| GeaneTrackRep with GEANE. [message #6205] |
Thu, 03 April 2008 11:03  |
Lia Lavezzi
Messages: 291
Registered: May 2007
Location: Torino
|
first-grade participant |
From: *PV.INFN.IT
|
|
Hallo Sebastian,
Susanna and I are having a look to the recotask/demo, running the runDemo.C with GEANE, just to test the GeaneTrackRep (we are using the new external packages).
We have some remarks:
1 - the Propagate(...) function of CbmGeanePro returns a value kTRUE or kFALSE to state whether the propagation was successful or not. This flag should be added in the GeaneTrackRep::extrapolate(...) function line where the actual propagation is performed by geane, i.e.:
_geane->Propagate(&par,&result,211);
just to stop the procedure when the GEANE propagation failed.
2 - in the function GeaneTrackRep::getMom(const DetPlane& pl) the TVector3 result is filled with the
variables statePred[3][0] and statePred[4][0] (which are v and w), instead of
statePred[1][0] and statePred[2][0] (which correspond to the derivatives v' and w').
The magnitude of the vector is correct because it is normalized to the momentum value, but the
components are wrong.
3 - pay attention also to the fact that the PDG code in the reconstruction is 211 (pion) while in the
simulation is 13 (muon).
In attachment you can find an example of GeaneTrackRep with these corrections (concerning the getMom
function all the calculations are explicit, but actually the simple change from statePred[3][0] and
statePred[4][0] to statePred[1][0] and statePred[2][0] should be enough).
We hope this could be useful, please have a look to this and tell us what do you think 
Ciao,
Lia and Susanna.
|
|
|
|
| Re: GeaneTrackRep with GEANE. [message #6206 is a reply to message #6205] |
Thu, 03 April 2008 11:13  |
Sebastian Neubert
Messages: 282
Registered: March 2006
Location: Munich
|
first-grade participant |
From: 213.215.163*
|
|
Hi Lia and Susanna!
Thank you very much for looking at this! I really appreciate that very much!
I will look at your code and put it to the repository. Concerning the PDG code, this of course has to become a variable of course, which you can intialize when you build a GeaneTrackRep. I will take care of this.
Kind Regards, Sebastian.
Sebastian Neubert
Technische Universität München
Department Physik E18
sneubert@e18.physik.tu-muenchen.de
tel: +49-8928912592
|
|
|
|
GSI Forum
Search
Help
Members
Login
Home
