GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung
Home » PANDA » PandaRoot » Bugs, Fixes, Releases » Crash in PndTpcRiemannTrackingTask
icon4.gif  Crash in PndTpcRiemannTrackingTask [message #12192] Mon, 04 July 2011 20:11 Go to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant
From: *40-79-r.retail.telecomitalia.it
Submitting TPC macros on the grid, only 20 events for each job, we have seen that almost 50% of the jobs are crashing because of PndTpcRiemannTrackingTask.
In particular, the last message is the following:

PndTpcRiemannTrackingTask::Exec; Event Number: 16
Looping over 58 riemann tracks to write out
RKTrackRep::setPDG particle -22 not known to TDatabasePDG -> abort


This comes from the MonteCarlo particle hypothesis (why is it used inside a pattern recognition task?). Somehow some "photon" tracks are reconstructed, multiplied by some charge and finally -22 is not taken by the particle DB -> crash.

If we run 500 or 1000 events we expect to have a very low success rate.
I would ask to tpc tracking experts to take a look and to fix it as soon as possible, becase this prevents us to run production on the grid.

Report message to a moderator

Re: Crash in PndTpcRiemannTrackingTask [message #12194 is a reply to message #12192] Tue, 05 July 2011 11:47 Go to previous messageGo to next message
Johannes Rauch is currently offline  Johannes Rauch
Messages: 41
Registered: September 2010
Location: TUM
continuous participant
 From: *natpool.mwn.de
Dear Stefano,

I commited a fix: If the pdg is -22 it is reset to 22, and an additional check is performed if the pdg is valid. If not, the track is skipped.

Please update the PndTpcRiemannTrackingTask to version 12579.

regards,

Johannes

Report message to a moderator

Re: Crash in PndTpcRiemannTrackingTask [message #12197 is a reply to message #12194] Tue, 05 July 2011 12:56 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
Are we sure that genfit+geane is able to propagate photons correctly? I would say that this was never tested before.

And I think it makes no sense for tracking of charged particles. How can be that a photon is creating a track, which is also reconstructed?

Report message to a moderator

Re: Crash in PndTpcRiemannTrackingTask [message #12234 is a reply to message #12192] Fri, 08 July 2011 13:14 Go to previous messageGo to next message
Johannes Rauch is currently offline  Johannes Rauch
Messages: 41
Registered: September 2010
Location: TUM
continuous participant
 From: *natpool.mwn.de
Hi,

I have been working on the bug that we adressed in the meeting yesterday (if the Riemann PR gets the track in a wrong direction, it flips the pdg id).

What I'm doing now is sorting the track in a way that the winding matches the charge gathered from the mc pdg (instead of flipping the charge resp. pdg id).

The problem here is that secondary tracks will get the pdg of the primary track they stem from, because for secondaries we have no more pdg ids available at this point.
We still have to think about a solution here.

Dear Analysts, could you please update the RiemannTrackingTask and check if efficiencies increase?

regards,

Johannes

Report message to a moderator

Re: Crash in PndTpcRiemannTrackingTask [message #12246 is a reply to message #12192] Tue, 12 July 2011 16:56 Go to previous messageGo to next message
Johannes Rauch is currently offline  Johannes Rauch
Messages: 41
Registered: September 2010
Location: TUM
continuous participant
 From: *natpool.mwn.de
Hi,

I just checked in a new version of the PndTpcRiemannTrackingTask.
You can now specify a pdg in the reco macro, which can then be flipped by the RiemannTracker according to the winding.
The default pdg is 211.
Please remember to switch off MCPid which is on by default.

For example:

tpcSPR->SetPDG(15);
tpcSPR->SetMCPid(false);


Please have a try.

regards,

Johannes

Report message to a moderator

Re: Crash in PndTpcRiemannTrackingTask [message #12259 is a reply to message #12192] Thu, 14 July 2011 10:25 Go to previous messageGo to next message
Dima Melnychuk is currently offline  Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...
first-grade participant
 From: *fuw.edu.pl
Hi,

I checked TPC performance for eta_c reconstruction after latest update of PndTpcRiemannTrackingTask.

I used kaon hypothesis in tpc reconstruction macro

tpcSPR->SetPDG(321);
tpcSPR->SetMCPid(false);

and eta_c peak appeared back but efficiency is still low in comparison with STT (5.4-5.9% for TPC vs 20-29% for STT).

So with 2000 events multiplicity of charged tracks:
index.php?t=getfile&id=6505&private=0
And maximum efficiency we can expect is 54%.

Mass distribution for phi and eta_c (using 4C-fit) efficiency - 5.9%
index.php?t=getfile&id=6506&private=0

Mass distribution for phi and eta_c (using vertex fit) efficiency - 5.4%
index.php?t=getfile&id=6507&private=0

Efficiencies in principle depend on cuts on chi2 of the fit and can be increased, but the 4-times difference between TPC and STT is with the same cuts.

Dima

Report message to a moderator

Re: Crash in PndTpcRiemannTrackingTask [message #12260 is a reply to message #12192] Thu, 14 July 2011 10:53 Go to previous messageGo to next message
Johannes Rauch is currently offline  Johannes Rauch
Messages: 41
Registered: September 2010
Location: TUM
continuous participant
 From: *natpool.mwn.de
Hi Dima,

its very strange that the TPC efficiency is so low. Yesterday I started to run the tpc macros in the /macro/run/tdrct/eta_c folder, to see whats happening, but I ran into a problem:
When runnin run_ana_eta_c_tpc.C, he complains:

Warning in <TTree::AddFriend>: FriendElement cbmsim in file evt_points_tpc.root has less entries 10 than its parent Tree: 11

How can I get it running?

regards,

Johannes

Report message to a moderator

Re: Crash in PndTpcRiemannTrackingTask [message #12261 is a reply to message #12192] Thu, 14 July 2011 10:58 Go to previous messageGo to next message
Dima Melnychuk is currently offline  Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...
first-grade participant
 From: *fuw.edu.pl
Hi Johannes,

This warning according to Tobias is harmless.
( http://forum.gsi.de/index.php?t=tree&th=3267&start=0&rid=78& amp;S=226ae0a4cb62b7d1cbdcce6c470541d1#page_top)
It happens when you use trunk version of base and additional entries in digitization, reconstruction should be simply ignored.

Dima

Report message to a moderator

Re: Crash in PndTpcRiemannTrackingTask [message #12263 is a reply to message #12261] Thu, 14 July 2011 11:06 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
However wil will run the old base packages on the grid, to avoid such a problem.

Report message to a moderator

Re: Crash in PndTpcRiemannTrackingTask [message #12280 is a reply to message #12246] Sat, 16 July 2011 23:27 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant
From: *16-79-r.retail.telecomitalia.it
Hi,
I have tried with the new july11 release (12718) to run eta_c events with tpc, ideal riemann and using the pion hypothesis.

The following plot shows the single kaon efficiency as a function of momentum, as showed in the last tracking meeting. In blue with ideal pid, in red with pion pid:

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

You can see that the low efficiency is connected really to the ideal id code.
I suppose tomorrow I will start new files on the grid setting the pion hypothesis in the tpcriemann, to have a "decent" efficiency.

Report message to a moderator

Re: Crash in PndTpcRiemannTrackingTask [message #12281 is a reply to message #12192] Sun, 17 July 2011 01:18 Go to previous messageGo to next message
Dima Melnychuk is currently offline  Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...
first-grade participant
 From: *play-internet.pl
Hi Stefano,

I have a clarification here.
Pion hypothesis gives better efficiencies than ideal pid, but not better then kaon hypothesis.

Here is the plot done with your macro when in reconstruction kaon hypothesis was used:
index.php?t=getfile&id=6517&private=0

I would say that efficiency is not worse that with pion hypotheis.

And this efficiency is with with reconstructed data, for which I obtained eta_c reconstruction efficiency 6%, i.e. 3-4 times lower than for STT.

So in the grid reconstruction macro specific hypotheis is prefered over ideal pid at the moment, but for eta_c channel it could be both pion or kaon.

Dima
  • Attachment: eff.png
    (Size: 9.47KB, Downloaded 549 times)

Report message to a moderator

Re: Crash in PndTpcRiemannTrackingTask [message #12282 is a reply to message #12281] Sun, 17 July 2011 10:05 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant
From: *16-79-r.retail.telecomitalia.it
I think the problem is independent from the particle hypothesis, simply the "ideal" option has some problems in the code and does not work properly.
However, it is too early to use specific particle hyp at the pattern recognition level, i think.

Report message to a moderator

Re: Crash in PndTpcRiemannTrackingTask [message #12283 is a reply to message #12282] Sun, 17 July 2011 16:58 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant
From: *16-79-r.retail.telecomitalia.it
Meanwhile I have tried to check also, tpc with pion riemann, phi and eta_c peaks with standard pion hypothesis (red) and with ideal hypothesis (blue) in PndPidCorrelator.

K+K- invariant mass

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

eta_c invariant mass

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

It seems that the ideal PidCorrelator (back extrapolation to IP using ideal particle hypothesis) produces a much sharper peak than with the default pion hypothesis. It seems a bit strange to me that the phi-phi invariant mass shows the same behaviour... however I will run the ideal pid correlator for tpc.

I will check the same also with stt.

Report message to a moderator

Re: Crash in PndTpcRiemannTrackingTask [message #12286 is a reply to message #12283] Mon, 18 July 2011 09:59 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
Sorry,
my fault,
in reality blue peak is with pion hypothesis, while in read the ideal hypothesis.
Something is going wrong in the PidCorrelator with ideal hyp, it seems in the backward extrapolation the phi angle is wrong. I will investigate. The same for stt...

Report message to a moderator

Re: Crash in PndTpcRiemannTrackingTask [message #12287 is a reply to message #12286] Mon, 18 July 2011 14:46 Go to previous message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant
From: *to.infn.it
Hi,
there was a problem in the ideal settings of the PndPidCorrelator, which is now fixed in the trunk (12725). Now phi and psi peaks ae almost the same w/ and w/o ideal pid.

Report message to a moderator

Previous Topic: source of low efficiency of eta_c reconstruction with TPC
Next Topic: Number of entries in TTree
Goto Forum:
  

[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Tue Dec 10 19:48:45 CET 2024

Total time taken to generate the page: 0.00867 seconds
.:: Contact :: Home :: Legal Notice :: Data Protection Policy ::.

Powered by: FUDforum 3.1.3.
Copyright ©2001-2023 FUDforum Bulletin Board Software