Linear Collider Forum



Home » General » General Questions » sid01 vs sid02
sid01 vs sid02 [message #1820] Tue, 28 July 2009 08:22
bweinert
Messages: 16
Registered: September 2008
Location: University of Rochester
Hi would there be any reason why using an sid02 event would change the results I'm getting for track matching. When I used the sid01file(mu_Theta4-176_1-50GeV-0-5000_SLIC-v2r3p10_geant4-v 9r0p1_LCPhys_sid01.slcio) I get real, virtual, and matched tracks that penetrate all the way to the edge of muon system. When I use the sid02 event
(mu_50.0GeV_Theta90_SLIC-v2r5p1_geant4-v9r1p2_LCPhys_sid02.s lcio) it only barely makes it to the muon detector for real, virtual, and matched hits. I've tried the same thing with higher energy sid02 events and I get the same result. Is there any reason why this might be the case? It looks like the problem is with the virtual track. I use TrackStepper from util/step to create the virtual tracks. Should I use something else now that I'm using sid02?

-Ben

I create the virtual tracks by using:
VirtPos[0] = rpVect[0];

VirtPos[1] = rpVect[1];

VirtPos[2] = rpVect[2];

VirtPosList.add(VirtPos);


I create the real hits by using (an example for the had calorimeter):
List<CalorimeterHit> hits = null;
try {
hits = event.get( CalorimeterHit.class, hcalHitmapName);
}
catch (Exception e) {}
if(hits==null) return;

for( int i = 0; i<hits.size(); ++i ) {
CalorimeterHit ihit = hits.get(i);
double[] Pos = ihit.getPosition();
matchHitsXYZ(ihit.getPosition());}


I match them by using:
protected void matchHitsXYZ(double [] realPos)
{

AIDA aida = AIDA.defaultInstance();
int nhitsTotal = 0;

// loop through virtual hits
for (double [] virtPos : VirtPosList) {

if(subdetName.equals(ecalSubdetName)) {
if (((realPos[0])-3)<=virtPos[0] && virtPos[0]<=((realPos[0])+3))
{

if(((realPos[1])-3)<=virtPos[1] && virtPos[1]<=((realPos[1])+3)){

if( ((realPos[2])-3)<=virtPos[2] && virtPos[2]<=((realPos[2])+3)){

aida.cloud2D("Y vs. X Matched").fill(realPos[0], realPos[1]);
aida.cloud1D("Rho
Matched" ).fill(Math.sqrt((realPos[0]*realPos[0])+(realPos[1]*realPos [1])+(realPos[2]*realPos[2])));
MatchedPosList.add(realPos);
foundHits.add(new BasicHep3Vector(realPos[0],
realPos[1], realPos[2]));

nhitsTotal++;
MatchedPosList.add(realPos);
}
} } }
else if(subdetName.equals(hcalSubdetName)) {
if (((realPos[0])-10)<=virtPos[0] && virtPos[0]<=((realPos[0])+10))
{

if(((realPos[1])-10)<=virtPos[1] && virtPos[1]<=((realPos[1])+10)){

if( ((realPos[2])-10)<=virtPos[2] && virtPos[2]<=((realPos[2])+10)){
aida.cloud2D("Y vs. X Matched").fill(realPos[0], realPos[1]);
aida.cloud1D("Rho
Matched" ).fill(Math.sqrt((realPos[0]*realPos[0])+(realPos[1]*realPos [1])+(realPos[2]*realPos[2])));
MatchedPosList.add(realPos);
foundHits.add(new BasicHep3Vector(realPos[0],
realPos[1], realPos[2]));

nhitsTotal++;
MatchedPosList.add(realPos);

}
} }}

}
}
 
Read Message
Previous Topic:location of SiD02 data files
Next Topic:2 things: thread links in mails and email address
Goto Forum:
  


Current Time: Tue Jan 23 23:41:36 Pacific Standard Time 2018
.:: Contact :: Home ::.

Powered by: FUDforum 3.0.1.
Copyright ©2001-2010 FUDforum Bulletin Board Software