GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Tracking » Ideal Tracking Tasks
Ideal Tracking Tasks [message #14523] Wed, 20 March 2013 13:21 Go to next message
Ralf Kliemt is currently offline  Ralf Kliemt
Messages: 507
Registered: May 2007
Location: GSI, Darmstadt
first-grade participant

From: *gsi.de
Dear trackers.

There are now two ideal tracking tasks. The concept is to match hits from STT or FTS to MC tracks and adding MVD+GEM hits afterwards.

The packages & classes are:

sttmvdtracking/PndSttMvdGemTrackingIdeal
fts/PndFtsTrackerIdeal



Set your track output branch with:

SetTrackOutput("...")


Smearing can be introduced by:

SetMomentumSmearing(spx,spy,spz)
SetRelativeMomentumSmearing(dp/p)
SetVertexSmearing(svx,svy,svz)
SetTrackingEfficiency(epsilon)



One can deactivate the additional detectors by:

SetMvdActivity(kFALSE)
SetGemActivity(kFALSE)




Cheers.
Ralf

[Updated on: Wed, 20 March 2013 13:27] by Moderator

Report message to a moderator

Re: Ideal Tracking Tasks [message #14531 is a reply to message #14523] Thu, 21 March 2013 18:50 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
Hi Ralf,
I have added a ideal reconstruction macro under pandaroot/macro/run/recoideal_complete.C

I put there:

 PndSttMvdGemTrackingIdeal* trackStt = new PndSttMvdGemTrackingIdeal();
  trackStt->SetRelativeMomentumSmearing(0.05);
  trackStt->SetVertexSmearing(0.05, 0.05, 0.05);
  trackStt->SetTrackingEfficiency(1.);
  trackStt->SetTrackOutput("SttMvdGemIdealTrack");
  fRun->AddTask(trackStt);


What I find strange is that if I run muons (13) at 1 GeV, with a momentum smearing of 5%, once I plot:

cbmsim.Draw("-1./SttMvdGemIdealTrack.fTrackParamLast.fQp>>aa(200,-1.5,0.5)","SttMvdGemIdealTrack.fRefIndex==0")


I obtain a width of 37 MeV instead of 50. Vertices seem fine.

Could you please check also in your case?
The second thing I found is that half of the track are failing with the kalman, which sounds a bit strange. But this I could check more deeply later.
Re: Ideal Tracking Tasks [message #14538 is a reply to message #14523] Mon, 25 March 2013 11:13 Go to previous messageGo to next message
Ralf Kliemt is currently offline  Ralf Kliemt
Messages: 507
Registered: May 2007
Location: GSI, Darmstadt
first-grade participant

From: *gsi.de
Hi Stefano.

I found that to be a bug indeed. Fixed in SVN.

Ralf
Re: Ideal Tracking Tasks [message #14543 is a reply to message #14538] Tue, 26 March 2013 09:29 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
Hi Ralf,
thanks for the fix, I have tried and now the resution seems correct.
Still I have problems in the kalman, half of the tracks fail the fit. Is it possible that the candidate is not sorted in the proper way? If you run the recoideal_complete.C, after the fit in substance I do not see any momentum improvement.
Re: Ideal Tracking Tasks [message #14544 is a reply to message #14523] Tue, 26 March 2013 10:24 Go to previous messageGo to next message
Ralf Kliemt is currently offline  Ralf Kliemt
Messages: 507
Registered: May 2007
Location: GSI, Darmstadt
first-grade participant

From: *gsi.de
Hm.

Doesn't sorting occur in the interface of GenFit?

Ralf
Re: Ideal Tracking Tasks [message #14545 is a reply to message #14523] Tue, 26 March 2013 11:12 Go to previous messageGo to next message
Ralf Kliemt is currently offline  Ralf Kliemt
Messages: 507
Registered: May 2007
Location: GSI, Darmstadt
first-grade participant

From: *gsi.de
Hi.

I modified the sorting parameter to be rho=(x*x+y*y)
Should be better now.

Ralf
Re: Ideal Tracking Tasks [message #14546 is a reply to message #14545] Tue, 26 March 2013 17:38 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
Hi Ralf,
now it breaks:

===PndRecoKalmanFit::Init() finished ===================================================
-I- PndMCTrackAssociator::Init
[INFO   ]  Branch: STTHitMix  not found in Tree  
[INFO   ]  Branch: STTHitMix  not found in Tree  
-W- PndMCTrackAssociator::Init: Using STTHit array
[INFO   ]  Branch: MVDHitsStripMix  not found in Tree  
[INFO   ]  Branch: MVDHitsStripMix  not found in Tree  
-W- PndMCTrackAssociator::Init: Using MVDHitsStrip array
[INFO   ]  Branch: MVDHitsPixelMix  not found in Tree  
[INFO   ]  Branch: MVDHitsPixelMix  not found in Tree  
-W- PndMCTrackAssociator::Init: Using MVDHitsPixel
[INFO   ] The number of entries in chain is 1000 

 *** Break *** segmentation violation



===========================================================
There was a crash.
This is the entire stack trace of all threads:
===========================================================
#0  0x00007fb69d05b43e in __libc_waitpid (pid=<optimized out>, stat_loc=0x7fffb1524580, options=0) at ../sysdeps/unix/sysv/linux/waitpid.c:32
#1  0x00007fb69cfe129e in do_system (line=0x10d73550 "/home/spataro/sep12/fairbuild/etc/gdb-backtrace.sh 28302 1>&2") at ../sysdeps/posix/system.c:149
#2  0x00007fb69df18787 in TUnixSystem::StackTrace() () from /home/spataro/sep12/fairbuild/lib/root/libCore.so.5.34
#3  0x00007fb69df1b063 in TUnixSystem::DispatchSignals(ESignals) () from /home/spataro/sep12/fairbuild/lib/root/libCore.so.5.34
#4  <signal handler called>
#5  0x00007fb691d20c82 in FairHit::GetX (this=0x0) at /home/spataro/sep12/pandaroot/base/FairHit.h:37
#6  0x00007fb68d01fb91 in PndSttMvdGemTrackingIdeal::Exec (this=0x24ecda0, option=0x7fb691d33360 "") at /home/spataro/sep12/pandaroot/sttmvdtracking/PndSttMvdGemTrackingIdeal.cxx:190
#7  0x00007fb69deaed23 in TTask::ExecuteTasks(char const*) () from /home/spataro/sep12/fairbuild/lib/root/libCore.so.5.34
#8  0x00007fb69deaf356 in TTask::ExecuteTask(char const*) () from /home/spataro/sep12/fairbuild/lib/root/libCore.so.5.34
#9  0x00007fb691c5ad34 in FairRunAna::Run (this=0x1989cf0, Ev_start=0, Ev_end=1000) at /home/spataro/sep12/pandaroot/base/FairRunAna.cxx:549
#10 0x00007fb691ccd991 in G__G__BaseDict_998_0_6 (result7=0x7fffb1532e40, funcname=0x19849a0 "", libp=0x7fffb1528430, hash=0) at /home/spataro/sep12/build/base/G__BaseDict.cxx:14861
#11 0x00007fb69c44d9ab in Cint::G__ExceptionWrapper(int (*)(G__value*, char const*, G__param*, int), G__value*, char*, G__param*, int) () from /home/spataro/sep12/fairbuild/lib/root/libCint.so.5.34
#12 0x00007fb69c4f32e1 in G__execute_call () from /home/spataro/sep12/fairbuild/lib/root/libCint.so.5.34
#13 0x00007fb69c4f36ce in G__call_cppfunc () from /home/spataro/sep12/fairbuild/lib/root/libCint.so.5.34
#14 0x00007fb69c4d4077 in G__interpret_func () from /home/spataro/sep12/fairbuild/lib/root/libCint.so.5.34
#15 0x00007fb69c4c0546 in G__getfunction () from /home/spataro/sep12/fairbuild/lib/root/libCint.so.5.34
#16 0x00007fb69c5ac180 in G__getstructmem(int, G__FastAllocString&, char*, int, char*, int*, G__var_array*, int) () from /home/spataro/sep12/fairbuild/lib/root/libCint.so.5.34
#17 0x00007fb69c5a364e in G__getvariable () from /home/spataro/sep12/fairbuild/lib/root/libCint.so.5.34
#18 0x00007fb69c499b92 in G__getitem () from /home/spataro/sep12/fairbuild/lib/root/libCint.so.5.34
#19 0x00007fb69c4a02a4 in G__getexpr () from /home/spataro/sep12/fairbuild/lib/root/libCint.so.5.34
#20 0x00007fb69c523c76 in G__exec_statement () from /home/spataro/sep12/fairbuild/lib/root/libCint.so.5.34
#21 0x00007fb69c4d21c7 in G__interpret_func () from /home/spataro/sep12/fairbuild/lib/root/libCint.so.5.34
#22 0x00007fb69c4c05ae in G__getfunction () from /home/spataro/sep12/fairbuild/lib/root/libCint.so.5.34
#23 0x00007fb69c49a331 in G__getitem () from /home/spataro/sep12/fairbuild/lib/root/libCint.so.5.34
#24 0x00007fb69c4a02a4 in G__getexpr () from /home/spataro/sep12/fairbuild/lib/root/libCint.so.5.34
#25 0x00007fb69c4abfac in G__calc_internal () from /home/spataro/sep12/fairbuild/lib/root/libCint.so.5.34
#26 0x00007fb69c5336df in G__process_cmd () from /home/spataro/sep12/fairbuild/lib/root/libCint.so.5.34
#27 0x00007fb69dee2b2a in TCint::ProcessLine(char const*, TInterpreter::EErrorCode*) () from /home/spataro/sep12/fairbuild/lib/root/libCore.so.5.34
#28 0x00007fb69ded9523 in TCint::ProcessLineSynch(char const*, TInterpreter::EErrorCode*) () from /home/spataro/sep12/fairbuild/lib/root/libCore.so.5.34
#29 0x00007fb69de49b94 in TApplication::ExecuteFile(char const*, int*, bool) () from /home/spataro/sep12/fairbuild/lib/root/libCore.so.5.34
#30 0x00007fb69de488cf in TApplication::ProcessLine(char const*, bool, int*) () from /home/spataro/sep12/fairbuild/lib/root/libCore.so.5.34
#31 0x00007fb69daa080a in TRint::Run(bool) () from /home/spataro/sep12/fairbuild/lib/root/libRint.so.5.34
#32 0x0000000000400f6c in main ()
===========================================================


The lines below might hint at the cause of the crash.
If they do not help you then please submit a bug report at
http://root.cern.ch/bugs. Please post the ENTIRE stack trace
from above as an attachment in addition to anything else
that might help us fixing this issue.
===========================================================
#5  0x00007fb691d20c82 in FairHit::GetX (this=0x0) at /home/spataro/sep12/pandaroot/base/FairHit.h:37
#6  0x00007fb68d01fb91 in PndSttMvdGemTrackingIdeal::Exec (this=0x24ecda0, option=0x7fb691d33360 "") at /home/spataro/sep12/pandaroot/sttmvdtracking/PndSttMvdGemTrackingIdeal.cxx:190
#7  0x00007fb69deaed23 in TTask::ExecuteTasks(char const*) () from /home/spataro/sep12/fairbuild/lib/root/libCore.so.5.34
#8  0x00007fb69deaf356 in TTask::ExecuteTask(char const*) () from /home/spataro/sep12/fairbuild/lib/root/libCore.so.5.34
#9  0x00007fb691c5ad34 in FairRunAna::Run (this=0x1989cf0, Ev_start=0, Ev_end=1000) at /home/spataro/sep12/pandaroot/base/FairRunAna.cxx:549
===========================================================


Root > Function recoideal_complete() busy flag cleared
Re: Ideal Tracking Tasks [message #14547 is a reply to message #14523] Wed, 27 March 2013 09:49 Go to previous messageGo to next message
Ralf Kliemt is currently offline  Ralf Kliemt
Messages: 507
Registered: May 2007
Location: GSI, Darmstadt
first-grade participant

From: *gsi.de
Strange.

I fixed it. However I see often messages that something in GEANE aborts or GenFit cannot converge. Might be just due to one of the tracks going into a "bad" region...

Ralf
Re: Ideal Tracking Tasks [message #14548 is a reply to message #14547] Wed, 27 March 2013 10:33 Go to previous message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
Much better now, thanks, now most (98%) of the primary tracks tracks are fitted, and I obtain some reasonable momentum resolution.
Many thanks.
Previous Topic: Modifications in Kalman
Next Topic: Check of GENFIT effects on reconstructed barrel tracks
Goto Forum:
  


Current Time: Mon Nov 11 08:54:42 CET 2024

Total time taken to generate the page: 0.00872 seconds