GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Bugs, Fixes, Releases » [FIXED] Ideal Tracking Bugs(?) when using FairLinks
[FIXED] Ideal Tracking Bugs(?) when using FairLinks [message #18211] Tue, 12 May 2015 16:38 Go to next message
André Zambanini is currently offline  André Zambanini
Messages: 17
Registered: February 2012
Location: FZ Jülich
occasional visitor
From: *ext.kfa-juelich.de
Hello everyone,

for my analysis I had a closer look at some events in the eventDisplay and found some strange behavior of the ideal track finding. Basically, I found two problems and my guess is, they are independent of each other, but I'm not sure.

General Information
First, some information beforehand.
FairSoft: mar15
FairRoot: master (fb738d60 from 26.03.2015)
PandaRoot: r27581

I'm simulating events with EvtGen (momentum 4.07 GeV/c) using the decay chain: pbar p -> Xi+ Xi(1690)- -> Lambdabar pi+ Lambda K- and the lambdas decaying to pi p. The main message from the decay for you is that I have a lot of displaced vertices with distances of several centimeters to the IP.

For the reconstruction I'm using ideal track finding, both with the old track finder (PndSttMvdGemTrackingIdeal) and the new one (PndMCIdealTrackFinderNewLinks). The simulation chain uses FairLinks (fRun->SetUseFairLinks(kTRUE)Wink.


Wrongly Assigned GEM Hits
The first thing Tobias and me noticed were wrongly assigned GEM hits. From what I understood from Tobias and Stefano, this is a known issue. To illustrate this a bit, see these screenshots here:

index.php?t=getfile&id=8385&private=0
index.php?t=getfile&id=8384&private=0

Both, the old and the new track finder seem to assign GEM hits which don't belong to the track. The white dots indicate the MC points associated with the selected track - which works fine for MVD (blue squares) and STT (purple) hits, but seems to be off for GEM (red) hits.

Along with this comes wrong track information, as you can see with the red and blue circles, which indicate the first and last track parameters, respectively.

Tobias and me had a closer look at the new track finder and found out, that the spurious hit assignment happens when more than one MC point belong to a hit. The quick fix introduced by us is to simply ignore those hits. This came to the PndMCMatchNewLinks class with r27667 in the trunk.

My conclusion here is, that for now this is okay but maybe someone should have a detailed look at this.

Track Reconstruction with Kalman Task
After resolving the spurious GEM hits, the reconstructed track parameters still looked quite odd. I did some comparison and found out, that in a few cases the Kalman task messes things up, both for genfit 1 and 2. The following screenshots all show the result of the new ideal track finder (PndMCMatchNewLinks):

index.php?t=getfile&id=8386&private=0
index.php?t=getfile&id=8387&private=0
index.php?t=getfile&id=8388&private=0

For this event, the Kalman filter produces strange tracks, independent of the genfit version. Most other events I had a look at seemed ok with genfit 2, while genfit 1 sometimes causes charge flips for instance (as visible in the screenshot above).
Ideally, I would leave the Kalman filter out, especially because for ideal tracking it seems a bit odd to use it. But this doesn't work for the PID, which seems to require the covariant matrices or something else filled by the Kalman.

So concluding, I have two questions: Any ideas, what could cause the Kalman to produce these results? And secondly, why is it required to have it in the first place, shouldn't ideal PID be based on MC information?

[Updated on: Mon, 06 July 2015 14:12]

Report message to a moderator

Re: Ideal Tracking Bugs(?) when using FairLinks [message #18213 is a reply to message #18211] Tue, 12 May 2015 18:57 Go to previous messageGo to next message
André Zambanini is currently offline  André Zambanini
Messages: 17
Registered: February 2012
Location: FZ Jülich
occasional visitor
From: *netcologne.de
What the GEM hit assignment is concerned, I tried to run with the newest PandaRoot version again. As it turns out, the old ideal track finder (PndSttMvdGemTrackingIdeal) seems to assign all GEM hits correctly now. I have tried ~20 events and the assignment problem didn't occur, previously the effect appeared every 2nd or 3rd event. A crosscheck with the new track finder (PndMCIdealTrackFinderNewLinks) shows, that this still requires the mentioned modification by skipping GEM hits with more than one point.

Something I found earlier as well for the PndSttMvdGemTrackingIdeal algorithm but was not at my focus: the STT hits are not always assigned, even though they match quite nicely between MVD and GEM hits of a track. I attached a screenshot showing this on an example event (seems to happen every 2-3 events for my channel).

index.php?t=getfile&id=8389&private=0

The same event with PndMCIdealTrackFinderNewLinks looks fine.

About the Kalman filter: I still find some events where the genfit 1 propagation results in last track parameters with different charge sign (as the example in my first post), but the genfit 2 propagation is much better. All the events I looked at have a quite good match between ideal track and ideal track after Kalman.
Still, the question remains: why do we need a Kalman filter for ideal tracks?

[Updated on: Mon, 18 May 2015 15:43]

Report message to a moderator

Re: Ideal Tracking Bugs(?) when using FairLinks [message #18214 is a reply to message #18213] Tue, 12 May 2015 23:51 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *ip71.fastwebnet.it
Can you provide a set of macro (.dec sim digi reco pid and event display) so that one can reproduce what you are seeing?

Quote:
Still, the question remains: why do we need a Kalman filter for ideal tracks?


Kalman over ideal pattern recognition is to have ideal efficiency but realistic momentum reconstruction.

Why are you using the master of FairRoot?

[Updated on: Tue, 12 May 2015 23:58]

Report message to a moderator

Re: Ideal Tracking Bugs(?) when using FairLinks [message #18218 is a reply to message #18213] Wed, 13 May 2015 12:15 Go to previous messageGo to next message
Lia Lavezzi
Messages: 291
Registered: May 2007
Location: Torino
first-grade participant

From: *to.infn.it
Quote:
Something I found earlier as well for the PndSttMvdGemTrackingIdeal algorithm but was not at my focus: the STT hits are not always assigned, even though they match quite nicely between MVD and GEM hits of a track. I attached a screenshot showing this on an exemplary event (seems to happen every 2-3 events for my channel).


In the PndSttMvdGemTrackingIdeal there is a cut on (I think) 25 stt hits: if you have more than 25 hits, only the first 25 are assigned to the track since if you have a looping particle and you assign all the hits to the track cand, the Kalman messes up the track. Is this the case? Does it happen when you have more than 25 stt hits?

Ciao,
Lia.
Re: Ideal Tracking Bugs(?) when using FairLinks [message #18224 is a reply to message #18211] Wed, 13 May 2015 12:55 Go to previous messageGo to next message
André Zambanini is currently offline  André Zambanini
Messages: 17
Registered: February 2012
Location: FZ Jülich
occasional visitor
From: *ext.kfa-juelich.de
Quote:
Why are you using the master of FairRoot?


The initial reason was because the eventDisplay wouldn't work on my machine and we tried to include the most recent changes in the FairLink system. But I don't think it is necessary anymore, I can try to go back to v15-03 again.

I'll also provide you with a minimal macro example as soon as I can get on with it.


Quote:
Is this the case? Does it happen when you have more than 25 stt hits?


More than 25 hits in one track or overall? If the former, than no, it doesn't assign any hits (see also the screenshot in my last reply, the highlighted points - white points - are the assigned MC points, but none of the STT hits are assigned).
Re: Ideal Tracking Bugs(?) when using FairLinks [message #18226 is a reply to message #18224] Wed, 13 May 2015 16:35 Go to previous messageGo to next message
Lia Lavezzi
Messages: 291
Registered: May 2007
Location: Torino
first-grade participant

From: *to.infn.it
Quote:
More than 25 hits in one track or overall?


I mean more than 25 in one track.

It is a very strange behaviour since nothing has been changed in the stt concerning the hit -> refindex -> mc point -> track id -> mc track chain... could it be just a problem of display? Did you print out the hits and see the stt hits are missing?

Lia.
Re: Ideal Tracking Bugs(?) when using FairLinks [message #18240 is a reply to message #18226] Fri, 15 May 2015 14:47 Go to previous messageGo to next message
Lia Lavezzi
Messages: 291
Registered: May 2007
Location: Torino
first-grade participant

From: *to.infn.it
The STT hit wrong assignement should be fixed with the newest version of the SttMvdGemIdeal tracking. So, everyone who is using the ideal track finder, please update the sttmvdtracking directory!
The cut on 25 hits was wrongly assigned to the total number of hits in the event instead of the number of hits in the track and from this the missing hit problem had origin. Moreover, with Stefano, we decided to remove totally the cut @ 25 hits because of a logic problem, as said by Stefano in an earlier message to the forum.

Ciao,
Lia.
Re: Ideal Tracking Bugs(?) when using FairLinks [message #18241 is a reply to message #18240] Fri, 15 May 2015 22:56 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *ip71.fastwebnet.it
Moreover, a bug in the ideal PR with new links was found, in the calculation of first and last point. Hopefully it will be fixed in the begin of next week.
Meanwhile, it would be good to check with your events if something strange still happens with the old PndSttMvdGemTrackingIdeal, since it uses refindex and should be independent from the new links developments.
Re: Ideal Tracking Bugs(?) when using FairLinks [message #18242 is a reply to message #18240] Sat, 16 May 2015 09:57 Go to previous messageGo to next message
André Zambanini is currently offline  André Zambanini
Messages: 17
Registered: February 2012
Location: FZ Jülich
occasional visitor
From: *netcologne.de
Lia Lavezzi wrote on Fri, 15 May 2015 14:47
The STT hit wrong assignement should be fixed with the newest version of the SttMvdGemIdeal tracking. [...]
The cut on 25 hits was wrongly assigned to the total number of hits in the event instead of the number of hits in the track and from this the missing hit problem had origin.


Ah, this explains it. Thanks for fixing it!


Stefano Spataro wrote on Fri, 15 May 2015 22:56
Moreover, a bug in the ideal PR with new links was found, in the calculation of first and last point. Hopefully it will be fixed in the begin of next week.
Meanwhile, it would be good to check with your events if something strange still happens with the old PndSttMvdGemTrackingIdeal, since it uses refindex and should be independent from the new links developments.


Yes, I'll check this beginning of next week. Thank you!

PS: I'm a bit inactive since we had a holiday on Thursday and thus a long weekend in Germany. Sorry for the inconvenience.
Re: Ideal Tracking Bugs(?) when using FairLinks [message #18243 is a reply to message #18211] Mon, 18 May 2015 15:12 Go to previous messageGo to next message
André Zambanini is currently offline  André Zambanini
Messages: 17
Registered: February 2012
Location: FZ Jülich
occasional visitor
From: *ext.kfa-juelich.de
The STT hits are now assigned properly using the SttMvdGemIdeal tracking. And I don't see any strange GEM hit assignments anymore. So it seems, that these issues have been resolved.

Let me know when the first/last track parameter issue has been resolved, then I'll crosscheck this as well.
Re: Ideal Tracking Bugs(?) when using FairLinks [message #18244 is a reply to message #18243] Mon, 18 May 2015 15:18 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
Just a comment: the first/last problem was found for the new ideal PR, not for the old one. If you use PndSttMvdGemIdealTracking you should see no problems there.
Re: Ideal Tracking Bugs(?) when using FairLinks [message #18252 is a reply to message #18243] Tue, 19 May 2015 16:14 Go to previous messageGo to next message
Lia Lavezzi
Messages: 291
Registered: May 2007
Location: Torino
first-grade participant

From: *to.infn.it
Hallo Andre' and everyone,
the fix for the first/last assignement has been uploaded to svn by Stefano right now. Please cross your fingers and check if everything is fine now Rolling Eyes
Cheers,
Lia.
Re: Ideal Tracking Bugs(?) when using FairLinks [message #18255 is a reply to message #18252] Wed, 20 May 2015 14:47 Go to previous messageGo to next message
Jennifer Pütz is currently offline  Jennifer Pütz
Messages: 47
Registered: April 2015
Location: FZ Juelich
continuous participant
From: *ikp.kfa-juelich.de
Hi everyone,

I checked the fixed version with a simple simulation (100 events) of pbar p -> Pi+ Pi-. And I found the same problem with the first/last point again for some of my events.
Below there is a screenshot for one of these events.

index.php?t=getfile&id=8402&private=0

The macros I used for the simualtion and the reconstruction are attached.

Cheers,

Jenny
Re: Ideal Tracking Bugs(?) when using FairLinks [message #18256 is a reply to message #18255] Wed, 20 May 2015 14:51 Go to previous messageGo to next message
Lia Lavezzi
Messages: 291
Registered: May 2007
Location: Torino
first-grade participant

From: *to.infn.it
Hi Jenny,
I see from the macros that you have the Kalman filter on. Have you tried to switch it off and to look at the events directly after the pattern recognition alone?

Cheers,
Lia.
Re: Ideal Tracking Bugs(?) when using FairLinks [message #18257 is a reply to message #18256] Wed, 20 May 2015 14:55 Go to previous messageGo to next message
Jennifer Pütz is currently offline  Jennifer Pütz
Messages: 47
Registered: April 2015
Location: FZ Juelich
continuous participant
From: *ikp.kfa-juelich.de
Hi Lia,

I didn't try it. But I will and let you know if this is working.

Cheers.
Re: Ideal Tracking Bugs(?) when using FairLinks [message #18258 is a reply to message #18256] Wed, 20 May 2015 15:26 Go to previous messageGo to next message
Jennifer Pütz is currently offline  Jennifer Pütz
Messages: 47
Registered: April 2015
Location: FZ Juelich
continuous participant
From: *ikp.kfa-juelich.de
Hi,

I switched the Kalman Filter off and the events looking fine.
Below there is another event, for which the ideal tracking without the Kalman Filter is working fine. But if I switch on the Kalman Filter the first and last hit are changing.

I'm looking at the track with the highlighted points.
index.php?t=getfile&id=8403&private=0

Without Kalman:

index.php?t=getfile&id=8404&private=0

With Kalman:
index.php?t=getfile&id=8405&private=0

Cheers,

Jenny
Re: Ideal Tracking Bugs(?) when using FairLinks [message #18259 is a reply to message #18258] Wed, 20 May 2015 15:34 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 37.227.45*
The problem of the Kalman is well known and cannot be fixed.
You can try to use genfit2, which in theory should solve the problem (and if not we can address the developers).
Re: Ideal Tracking Bugs(?) when using FairLinks [message #18260 is a reply to message #18259] Wed, 20 May 2015 15:40 Go to previous message
Lia Lavezzi
Messages: 291
Registered: May 2007
Location: Torino
first-grade participant

From: *to.infn.it
I agree with Stefano. The fixes we made in these days were in the pattern recognition stage, not in the Kalman.
Let's hope genfit2 has that bug fixed (but I never tried it).
Anyway, glad to hear that at least before the Kalman everything is fine now.
Cheers,
Lia.
Previous Topic: Pion Measurement Asymmetry in GEANT3
Next Topic: [FIXED] Problem with PndVertexFitter for particles with neutral charge
Goto Forum:
  


Current Time: Sun Apr 28 02:02:03 CEST 2024

Total time taken to generate the page: 0.01015 seconds