GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Tracking » GEM tracking
Re: GEM tracking [message #8793 is a reply to message #8364] Thu, 04 June 2009 18:10 Go to previous messageGo to previous message
Radoslaw Karabowicz is currently offline  Radoslaw Karabowicz
Messages: 108
Registered: June 2004
Location: GSI
continuous participant
From: *gsi.de
Dear Christian,

I have to apologize for my late answering, but the test and the results analysis took a bit longer I expected. The change from linux to mac didn't help either.
But here I am with something interesting, I hope.

I took the box generator, and used:
FairBoxGenerator* boxGen = new FairBoxGenerator(211,1);
boxGen->SetThetaRange(theta -0.01,theta +0.01);
boxGen->SetPhiRange (phi -0.01,phi +0.01);
boxGen->SetPRange (momentum-0.01,momentum+0.01);
primGen->AddGenerator(boxGen);

I have been varying the theta parameter from 2 to 40, phi from 0 to 360, and momentum was {0.1,0.2,0.5,1.0,2.0,4.0,8.0}. I have generated 1000 pions in each of the bins and tried to reconstruct them using start parameters:

TVector3 StartPos = TVector3 (0.,0.,0.);
TVector3 StartPosErr = TVector3(0.01,0.01,0.01);

TVector3 StartMom = TVector3 (1.,0.,1.);
StartMom.SetMagThetaPhi(fMomentum+0.05,
TMath::DegToRad()*(Double_t)fTheta,
TMath::DegToRad()*(Double_t)fPhi);
TVector3 StartMomErr = TVector3(0.1*StartMom);

where fMomentum, fTheta and fPhi were taken from the real parameters of the simulated pions.

As the results I have been looking at the distribution of the reconstructed momentum:
index.php?t=getfile&id=5405&private=0

This plot shows the distribution of reconstructed momentum magnitude for pions shot at phi = 0, theta = 10. The distribution look OK, maybe a bit to broad, but at least at correct momentum.

Taking the different distributions for pions shot at different theta and phi, I was able to construct a map of the mean values of the reconstructed momentum magnitude as functions of theta and phi. The resulting map for 2 GeV pions is:
index.php?t=getfile&id=5406&private=0

It certainly looks strange, but the behaviour is consistently repeated with different particle momenta (here 0.5 GeV):
index.php?t=getfile&id=5407&private=0

Have you ever seen anything so strange? For some phi angles the tracking performs good for all theta, but for some angles the fitting gives WAY too small or too large momenta, consistently for all theta. Even worse, the larger the theta angle, the greater the discrepancy.

The only thing I am afraid of is that I messed with the angles somehow, i mean rads and degs, I am checking it now... (well, tomorrow).

yours,
radek
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: change in genfit exception handling
Next Topic: GEM Tracking inside LHETRACK
Goto Forum:
  


Current Time: Sun Jun 16 23:42:06 CEST 2024

Total time taken to generate the page: 0.01079 seconds