| Ideal Tracking Tasks [message #14523] |
Wed, 20 March 2013 13:21  |
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  |
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 #14546 is a reply to message #14545] |
Tue, 26 March 2013 17:38 |
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
|
|
|
|
|
|
|
|
GSI Forum
Search
Help
Members
Login
Home
