| About FairGeanePro [message #15561] |
Thu, 10 October 2013 18:20  |
Jifeng Hu
Messages: 31
Registered: October 2012
|
continuous participant |
From: *physik.uni-giessen.de
|
|
Hi, All,
who has experiences of FairGeanePro propagation, please have a look at my usage and results, I think I got some wrong results.
For 2 GeV momentum pion from the origin (0,0,0), isotropically generated, the flying path is expected to be the EMC shower center position. To do that, the codes are listed (see also PndPidCorrelator.cxx ):
0. Bool_t PndPidCorrelator::GetEmcInfo(FairTrackParH* helix, PndPidCandidate* pidCand) {
1. FairGeanePro *fProEmc = new FairGeanePro();
2. for (Int_t ee = 0; ee<emcEntries; ee++){
3. PndEmcCluster *emcHit = (PndEmcCluster*)fEmcCluster->At(ee);
4. TVector3 emcPos = emcHit->where();
5. fProEmc->SetPoint(emcPos);
6. fProEmc->PropagateToPCA(1, 1);
7. FairTrackParH *fRes= new FairTrackParH();
8. Bool_t rc = fProEmc->Propagate(helix, fRes, fPidHyp*pidCand->GetCharge()); // First propagation at module
9. if (!rc) continue;
10. emcGLength = fProEmc->GetLengthAtPCA();
11. }
22.}
emcGLength is extracted as the flying path. see its distribution in barrel.eps for Barrel Calorimeter, and so on.
see the Panda TDR:
http://www-panda.gsi.de/archive/public/panda_tpr.pdf
the geometry tell us the flying path is wrong.
Someone can tell me where I incorrectly make use of FairGeanePro class.
-
Attachment: backward.eps
(Size: 10.72KB, Downloaded 389 times)
-
Attachment: forward.eps
(Size: 35.28KB, Downloaded 406 times)
-
Attachment: barrel.eps
(Size: 328.60KB, Downloaded 405 times)
-
Attachment: shashylik.eps
(Size: 17.37KB, Downloaded 417 times)
|
|
|
|
GSI Forum
Search
Help
Members
Login
Home
