Linear Collider Forum



Today's Messages (off)  | Unanswered Messages (on)

Forum: General Questions
 Topic: Detector Specifications
Detector Specifications [message #1885] Thu, 15 October 2009 08:05
bweinert
Messages: 16
Registered: September 2008
Location: University of Rochester
Hi, I'm running into a problem with the layer thickness of the muon system. The program updates the number of layers, when I run SID01 vs SID02 events, but it doesn't seem to update the layer thickness. It leaves me with 11 layers of 20mm thickness, so the detector stops early. This problem only occurs for the virtual tracks, not the real tracks. I've been trying to track down where the code gets the layer thickness from and I haven't found an actual number yet. Is there anywhere that I should look? Thank you.

-Ben
 Topic: 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);

}
} }}

}
}
 Topic: Registering and using forum.linearcollider.org
Registering and using forum.linearcollider.org [message #877] Sun, 27 May 2007 16:01
tonyj
Messages: 138
Registered: January 2004

I have written up some additional instructions on how to register with and use the forum.linearcollider.org web site. You can find these instructions at:

http://confluence.slac.stanford.edu/x/d4Q

These instructions are meant to complement the help built-in to the forum web site.

Feedback and suggestions are welcome.

[Updated on: Wed, 06 June 2007 00:15]


Current Time: Fri Oct 19 08:03:17 Pacific Daylight Time 2018
.:: Contact :: Home ::.

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