Home » PANDA » PandaRoot » Tracking » geane propagation: floating point exception
Re: geane propagation: floating point exception [message #5160 is a reply to message #5130] |
Thu, 04 October 2007 17:06 |
Andrea Fontana
Messages: 32 Registered: May 2007
|
continuous participant |
From: *pv.infn.it
|
|
Hi Sebastian,
we are writing with an update on our progress with genfit.
First of all, the error of 0 momentum was due to a series of
reasons:
1) the PDG code of the track was hardwired to 211 (pion)
in GeaneTrackRep and in the macros we have 13 (muon).
2) the simulated momentum of 0.5 is too low. By increasing it
to 1.5 it works. We know that we need to track also slow
particles, but this will come...
3) finally we have protected the CbmTrackParP class by rejecting
the tracks with 0 momentum that are now skipped. This has been
commited to the svn repository.
Moreover, to generalize the use of the detector planes and
to make them independent from the particle direction (in
the demo they were orthogonal to z) we have modified the
DemoRecoHit class. We have created a new class DemoRecoHit
which we invite you to look at (it is in the attachment as I
cannot commit to recotasks/demo): here we defined a plane
at each hit with the correct orientation orthogonal to the
track momentum.
This is now general and can be used also with an isotropic
generator: for brevity we have done it only in the third
constructor DemoRecoHit(CbmMCPoint *point).
Having fixed this, we are now extrapolating muons of 1.5 GeV
in the TPC isotropically generated. The Kalman filter seems
to be working, but it fails on some tracks. If we shoot along
one axis (take X, to cross the TPC) the result is much better
and we are now trying to understand this.
We attach the plots of the reconstructed total momentum:
in the file X.ps we shoot along X, in Isotropic.ps isotropically.
Other thoughts in our minds:
- is there a way to reject one hit in the extrapolation/filtering
if something goes wrong? We have seen the statusFlag in
AbsTrackRep. Maybe we could also modify all the extrapolate
methods to return a bool that can be used to decide.
- We had to comment the call to FillGeotrack in DemoKalmanTask.
What is it used for?
For now this is all. Anyway we are very happy!!!
Ciao,
Andrea and Lia
-
Attachment: X.ps
(Size: 10.85KB, Downloaded 416 times)
-
Attachment: Isotropic.ps
(Size: 11.78KB, Downloaded 448 times)
-
Attachment: DemoRecoHit.cxx
(Size: 3.65KB, Downloaded 437 times)
|
|
|
Goto Forum:
Current Time: Sat Oct 12 11:09:33 CEST 2024
Total time taken to generate the page: 0.00930 seconds
|