Home » PANDA » PandaRoot » Bugs, Fixes, Releases » analysis for TOF studies, fullsim crashes
analysis for TOF studies, fullsim crashes [message #9183] |
Tue, 11 August 2009 16:44 |
Bernhard Roth
Messages: 8 Registered: February 2009 Location: Bochum
|
occasional visitor |
From: *ep1.ruhr-uni-bochum.de
|
|
Hi,
I am doing TOF studies and also want to compare the results of my analysis between our BaBar-like software and PANDAroot. Therefore I'm analyzing the reaction
pbar p -> Phi Phi Pi0; Phi -> K+ K-;
all PHSP and at 2GeV beam momentum.
As background I'm considering Pi+ Pi- K+ K- Pi0
All with full simulation.
In PANDAroot I want to use the scripts from the PandaRootTutTorino09 which is tutorials/analysis/go_stt_evtgen.sh.
And I'm using the latest stable version from 10.8.2009, running on a GSI 32bit machine.
As first step I had to change run_reco_sttcombi.C:
Toggle Spoiler
- PndTpcLheTrackFinder* trackFinder = new PndTpcLheTrackFinder();
+ PndLheTrackFinder* trackFinder = new PndLheTrackFinder();
fRun->AddTask(trackFinder);
- PndTpcLheTrackFitter* trackFitter = new PndTpcLheTrackFitter("fitting");
+ PndLheTrackFitter* trackFitter = new PndLheTrackFitter("fitting");
fRun->AddTask(trackFitter);
(output of svn diff run_reco_sttcombi.C)
so removed the "Tpc" letters
Then it crashes whether I run over 10 events or more at different points.
On more events it crashes in run_kalman_stt.C with the following error messages in 4-kalman.log:
Toggle Spoiler
#1 0x412078c0 in __DTOR_END__ () from /lib/libc.so.6
#2 0x41118442 in do_system () from /lib/libc.so.6
#3 0x41093c5f in system () from /lib/libpthread.so.0
#4 0x40237363 in TUnixSystem::Exec (this=0x80e2a78,
at core/unix/src/TUnixSystem.cxx:1941
#5 0x40237836 in TUnixSystem::StackTrace (this=0x80e2a78) at core/unix/src/TUnixSystem.cxx:2121
#6 0x402356f5 in TUnixSystem::DispatchSignals (this=0x80e2a78, sig=kSigFloatingException) at core/unix/src/TUnixSystem.cxx:1089
#7 0x402336b4 in SigHandler (sig=kSigFloatingException) at core/unix/src/TUnixSystem.cxx:351
#8 0x4023a6d3 in sighandler (sig=8) at core/unix/src/TUnixSystem.cxx:3344
#9 0x41092825 in __pthread_sighandler () from /lib/libpthread.so.0
#10 <signal handler called>
#11 0x45e8ed38 in FairGeanePro::Propagate (this=0x9943150, TStart=0xbfd30760, TEnd=0xbfd30e10, PDG=-13)
at /d/panda02/broth/proot/stable/geane/FairGeanePro.cxx:157
#12 0x4439ca5a in GeaneTrackRep::extrapolate (this=0xe5430a8, pl=@0xbfd31ce0, statePred=@0xbfd31e60, covPred=@0xbfd31d60)
at /d/panda02/broth/proot/stable/trackrep/GeaneTrackRep.cxx:167
#13 0x4431b285 in Kalman::processHit (this=0xbfd333a0, tr=0xe542370, ihit=0, irep=0, rejectOutlier=false)
at /d/panda02/broth/proot/stable/genfit/Kalman.cxx:248
#14 0x4431a44f in Kalman::fittingPass (this=0xbfd333a0, trk=0xe542370, direction=1)
at /d/panda02/broth/proot/stable/genfit/Kalman.cxx:140
#15 0x44319be8 in Kalman::processTrack (this=0xbfd333a0, trk=0xe542370) at /d/panda02/broth/proot/stable/genfit/Kalman.cxx:38
#16 0x45e16ae5 in PndLheKalmanTask::Exec (this=0xb206128, opt=0x43d530a0 "")
at /d/panda02/broth/proot/stable/lhetrack/PndLheKalmanTask.cxx:238
#17 0x401b16af in TTask::ExecuteTasks (this=0x86f11d8, option=0x43d530a0 "") at core/base/src/TTask.cxx:298
#18 0x401b14b1 in TTask::ExecuteTask (this=0x86f11d8, option=0x43d530a0 "") at core/base/src/TTask.cxx:261
#19 0x43cdc46d in FairRunAna::Run (this=0x86f1150, Ev_start=0, Ev_end=100) at /d/panda02/broth/proot/stable/base/FairRunAna.cxx:248
#20 0x43d1a746 in G__FairDict_532_0_5 (result7=0xbfd3a900, funcname=0x86ef1a8 "\001", libp=0xbfd34a80, hash=0)
at /d/panda02/broth/proot/stable/build/base/FairDict.cxx:9067
#21 0x407b6126 in Cint::G__ExceptionWrapper (funcp=0x43d1a642 <G__FairDict_532_0_5>, result7=0xbfd3a900, funcname=0x86ef1a8 "\001",
libp=0xbfd34a80, hash=0) at cint/cint/src/Api.cxx:364
#22 0x408757f5 in G__execute_call (result7=0xbfd3a900, libp=0xbfd34a80, ifunc=0x86ef1a8, ifn=0) at cint/cint/src/newlink.cxx:2305
#23 0x40875ed8 in G__call_cppfunc (result7=0xbfd3a900, libp=0xbfd34a80, ifunc=0x86ef1a8, ifn=0) at cint/cint/src/newlink.cxx:2471
#24 0x40855818 in G__interpret_func (result7=0xbfd3a900, funcname=0xbfd3a500 "Run", libp=0xbfd34a80, hash=309, p_ifunc=0x86ef1a8,
funcmatch=1, memfunc_flag=1) at cint/cint/src/ifunc.cxx:5245
#25 0x40834ca1 in G__getfunction (item=0xbfd3d356 "Run(0,nEvents)", known3=0xbfd3cb9c, memfunc_flag=1) at cint/cint/src/func.cxx:2534
#26 0x40940b6b in G__getstructmem (store_var_type=112, varname=0xbfd3ab90 "", membername=0xbfd3d356 "Run(0,nEvents)",
tagname=0xbfd3ada0 "fRun", known2=0xbfd3cb9c, varglobal=0x409ebec0, objptr=2) at cint/cint/src/var.cxx:6623
#27 0x409329d7 in G__getvariable (item=0xbfd3d350 "fRun->Run(0,nEvents)", known=0xbfd3cb9c, varglobal=0x409ebec0, varlocal=0xbfd3f8a0)
at cint/cint/src/var.cxx:5252
#28 0x40825f01 in G__getitem (item=0xbfd3d350 "fRun->Run(0,nEvents)") at cint/cint/src/expr.cxx:1884
#29 0x4082387b in G__getexpr (expression=0xbfd3ec20 "fRun->Run(0,nEvents)") at cint/cint/src/expr.cxx:1470
#30 0x4089fde5 in G__exec_function (statement=0xbfd3ec20 "fRun->Run(0,nEvents)", pc=0xbfd3f04c, piout=0xbfd3f044,
plargestep=0xbfd3f034, presult=0xbfd3ebf0) at cint/cint/src/parse.cxx:601
#31 0x408af2a6 in G__exec_statement (mparen=0xbfd3f56c) at cint/cint/src/parse.cxx:6972
#32 0x408587a3 in G__interpret_func (result7=0xbfd45d40, funcname=0xbfd45940 "run_kalman_stt", libp=0xbfd3fec0, hash=1506,
p_ifunc=0x82b4fd0, funcmatch=1, memfunc_flag=0) at cint/cint/src/ifunc.cxx:6080
#33 0x4083584a in G__getfunction (item=0xbfd46640 "run_kalman_stt(100)", known3=0xbfd45e8c, memfunc_flag=0)
at cint/cint/src/func.cxx:2745
#34 0x40826044 in G__getitem (item=0xbfd46640 "run_kalman_stt(100)") at cint/cint/src/expr.cxx:1896
#35 0x4082387b in G__getexpr (expression=0x82f39e8 "run_kalman_stt(100)") at cint/cint/src/expr.cxx:1470
#36 0x408114ff in G__calc_internal (exprwithspace=0xbfd49d90 "run_kalman_stt(100)") at cint/cint/src/expr.cxx:1061
#37 0x408b7d27 in G__process_cmd (line=0x406e9b10 "1", prompt=0x80e62bc "", more=0x80e62b4, err=0xbfd4a60c, rslt=0xbfd4a610)
at cint/cint/src/pause.cxx:2234
...
whole file is in
/d/panda02/broth/proot/stable/tutorials/analysis/data_4crashed
In case of running only over 10 events it finishes run_kalman_stt.C successfully, but crashes in the next step, in makeTCands_stt.C. The error message in 5-microwriter.log is here:
Toggle Spoiler
#1 0x412078c0 in __DTOR_END__ () from /lib/libc.so.6
#2 0x41118442 in do_system () from /lib/libc.so.6
#3 0x41093c5f in system () from /lib/libpthread.so.0
#4 0x40237363 in TUnixSystem::Exec (this=0x80e2a78,
at core/unix/src/TUnixSystem.cxx:1941
#5 0x40237836 in TUnixSystem::StackTrace (this=0x80e2a78) at core/unix/src/TUnixSystem.cxx:2121
#6 0x402356f5 in TUnixSystem::DispatchSignals (this=0x80e2a78, sig=kSigSegmentationViolation) at core/unix/src/TUnixSystem.cxx:1089
#7 0x402336b4 in SigHandler (sig=kSigSegmentationViolation) at core/unix/src/TUnixSystem.cxx:351
#8 0x4023a6d3 in sighandler (sig=11) at core/unix/src/TUnixSystem.cxx:3344
#9 0x41092825 in __pthread_sighandler () from /lib/libpthread.so.0
#10 <signal handler called>
#11 0x458b1024 in typeinfo name for PndMvdConvertApvTask () from /d/panda02/broth/proot/stable/build/lib/libMvd.so
#12 0x443171e6 in Track::getCardinalRep (this=0x96905f0) at Track.h:182
#13 0x44333b1b in Track::getPos (this=0x96905f0) at Track.h:202
#14 0x462bc2bc in PndMicroWriter::Exec (this=0xb204300, opt=0x43d530a0 "")
at /d/panda02/broth/proot/stable/PndTools/AnalysisTools/PndMicroWriter.cxx:331
#15 0x401b16af in TTask::ExecuteTasks (this=0x86f26b0, option=0x43d530a0 "") at core/base/src/TTask.cxx:298
#16 0x401b14b1 in TTask::ExecuteTask (this=0x86f26b0, option=0x43d530a0 "") at core/base/src/TTask.cxx:261
#17 0x43cdc46d in FairRunAna::Run (this=0x86f2628, Ev_start=0, Ev_end=10) at /d/panda02/broth/proot/stable/base/FairRunAna.cxx:248
#18 0x43d1a746 in G__FairDict_532_0_5 (result7=0xbfb80740, funcname=0x86f0680 "\001", libp=0xbfb7a8c0, hash=0)
at /d/panda02/broth/proot/stable/build/base/FairDict.cxx:9067
#19 0x407b6126 in Cint::G__ExceptionWrapper (funcp=0x43d1a642 <G__FairDict_532_0_5>, result7=0xbfb80740, funcname=0x86f0680 "\001",
libp=0xbfb7a8c0, hash=0) at cint/cint/src/Api.cxx:364
#20 0x408757f5 in G__execute_call (result7=0xbfb80740, libp=0xbfb7a8c0, ifunc=0x86f0680, ifn=0) at cint/cint/src/newlink.cxx:2305
#21 0x40875ed8 in G__call_cppfunc (result7=0xbfb80740, libp=0xbfb7a8c0, ifunc=0x86f0680, ifn=0) at cint/cint/src/newlink.cxx:2471
#22 0x40855818 in G__interpret_func (result7=0xbfb80740, funcname=0xbfb80340 "Run", libp=0xbfb7a8c0, hash=309, p_ifunc=0x86f0680,
funcmatch=1, memfunc_flag=1) at cint/cint/src/ifunc.cxx:5245
#23 0x40834ca1 in G__getfunction (item=0xbfb83196 "Run(0,nevts)", known3=0xbfb829dc, memfunc_flag=1) at cint/cint/src/func.cxx:2534
#24 0x40940b6b in G__getstructmem (store_var_type=112, varname=0xbfb809d0 "", membername=0xbfb83196 "Run(0,nevts)",
tagname=0xbfb80be0 "fRun", known2=0xbfb829dc, varglobal=0x409ebec0, objptr=2) at cint/cint/src/var.cxx:6623
#25 0x409329d7 in G__getvariable (item=0xbfb83190 "fRun->Run(0,nevts)", known=0xbfb829dc, varglobal=0x409ebec0, varlocal=0xbfb856e0)
at cint/cint/src/var.cxx:5252
#26 0x40825f01 in G__getitem (item=0xbfb83190 "fRun->Run(0,nevts)") at cint/cint/src/expr.cxx:1884
#27 0x4082387b in G__getexpr (expression=0xbfb84a60 "fRun->Run(0,nevts)") at cint/cint/src/expr.cxx:1470
#28 0x4089fde5 in G__exec_function (statement=0xbfb84a60 "fRun->Run(0,nevts)", pc=0xbfb84e8c, piout=0xbfb84e84, plargestep=0xbfb84e74,
presult=0xbfb84a30) at cint/cint/src/parse.cxx:601
#29 0x408af2a6 in G__exec_statement (mparen=0xbfb853ac) at cint/cint/src/parse.cxx:6972
#30 0x408587a3 in G__interpret_func (result7=0xbfb8bb80, funcname=0xbfb8b780 "makeTCands_stt", libp=0xbfb85d00, hash=1429,
p_ifunc=0x82b4fd0, funcmatch=1, memfunc_flag=0) at cint/cint/src/ifunc.cxx:6080
#31 0x4083584a in G__getfunction (item=0xbfb8c480 "makeTCands_stt(10)", known3=0xbfb8bccc, memfunc_flag=0)
at cint/cint/src/func.cxx:2745
#32 0x40826044 in G__getitem (item=0xbfb8c480 "makeTCands_stt(10)") at cint/cint/src/expr.cxx:1896
#33 0x4082387b in G__getexpr (expression=0x82f31f0 "makeTCands_stt(10)") at cint/cint/src/expr.cxx:1470
#34 0x408114ff in G__calc_internal (exprwithspace=0xbfb8fbd0 "makeTCands_stt(10)") at cint/cint/src/expr.cxx:1061
#35 0x408b7d27 in G__process_cmd (line=0x406e9b0f "EventSummary_]", prompt=0x80e62bc "", more=0x80e62b4, err=0xbfb9044c,
rslt=0xbfb90450) at cint/cint/src/pause.cxx:2234
#36 0x4021e958 in TCint::ProcessLine (this=0x80e6298, line=0x406e9b0f "EventSummary_]", error=0xbfb92c54)
at core/meta/src/TCint.cxx:339
#37 0x4021ed70 in TCint::ProcessLineSynch (this=0x80e6298, line=0x406e9b0f "EventSummary_]", error=0xbfb92c54)
at core/meta/src/TCint.cxx:406
#38 0x4013581a in TApplication::ExecuteFile (file=0xbfb90bf3 "makeTCands_stt.C(10)", error=0xbfb92c54)
at core/base/src/TApplication.cxx:938
#39 0x40135094 in TApplication::ProcessFile (this=0x81078b0, file=0xbfb90bf3 "makeTCands_stt.C(10)", error=0xbfb92c54)
at core/base/src/TApplication.cxx:825
...
whole file is in
/d/panda02/broth/proot/stable/tutorials/analysis/data_5crashed
Did I use the correct macro with correct input?
Can some expert help me in this context?
Cheers,
Bernhard
P.S.
To compare the results between PANDAroot and the BaBar-like software it would be great to have the same data-sets available.
For signal (ppbarsystem -> Phi Phi Pi0, see above) I have now 100k events and for background Pi+ Pi- K+ K- Pi0 10M events for the BaBar-like software available. Is there an opportunity to order such an amount of data for doing analysis with PANDAroot? (all fullsim)
|
|
|
|
Re: analysis for TOF studies, fullsim crashes [message #9185 is a reply to message #9184] |
Tue, 11 August 2009 17:23 |
StefanoSpataro
Messages: 2736 Registered: June 2005 Location: Torino
|
first-grade participant |
From: *56-82-r.retail.telecomitalia.it
|
|
Hi,
about the first kind of errors, the structure of lhetrack was strongly changed with respect to the Torino meeting (therefore from the tutorial).
This means that the task name is a bit different (your change in reco), but I cannot see your crash in kalman macro without full access to the code. The only thing that I can suggest is to comment out the pidmaker in the reco macro (this could give some problems after with the new kalman task).
About the makeTCand, the code in rho was not updated according to the new structure, this is simply the reason why it is not working. I know Klaus was working to update it, but then he left for vacation.
If you want to reproduce your studies, I would suggest to use the svn release for the tutorial, exactly that revision numer (which was before the change).
|
|
|
|
|
|
|
Re: analysis for TOF studies, fullsim crashes [message #9199 is a reply to message #9183] |
Wed, 12 August 2009 16:22 |
Bernhard Roth
Messages: 8 Registered: February 2009 Location: Bochum
|
occasional visitor |
From: *ep1.ruhr-uni-bochum.de
|
|
Hi,
I'm now in revision 5813 and still running the macros from PandaRootTutTorino09 for fullsim of decay
p pbar -> phi phi pi0; phi -> K+ K- at 2 GeV
It finishes successfully on 10 and 100 events, but crashes on 1000 events in run_kalman_stt.C.
The error message in 4-kalman.log is:
Toggle Spoiler
#1 0x412078c0 in __DTOR_END__ () from /lib/libc.so.6
#2 0x41118442 in do_system () from /lib/libc.so.6
#3 0x41093c5f in system () from /lib/libpthread.so.0
#4 0x40237363 in TUnixSystem::Exec (this=0x80e2a78,
at core/unix/src/TUnixSystem.cxx:1941
#5 0x40237836 in TUnixSystem::StackTrace (this=0x80e2a78) at core/unix/src/TUnixSystem.cxx:2121
#6 0x402356f5 in TUnixSystem::DispatchSignals (this=0x80e2a78, sig=kSigFloatingException)
at core/unix/src/TUnixSystem.cxx:1089
#7 0x402336b4 in SigHandler (sig=kSigFloatingException) at core/unix/src/TUnixSystem.cxx:351
#8 0x4023a6d3 in sighandler (sig=8) at core/unix/src/TUnixSystem.cxx:3344
#9 0x41092825 in __pthread_sighandler () from /lib/libpthread.so.0
#10 <signal handler called>
#11 0x4601c2e3 in xmm55_ (a=0x46462e20, b=0x46462ee8, c=0x46462ee8) at matx55/xmm55.F:42
#12 0x460c310e in trprfn_ (x1=0x464341d8, p1=0x464341e4, h1=0x464341f0, x2=0x46434214, p2=0x46434220, h2=0x4643422c,
ch=0x46434250, xl=0x462b380c, r=0xbf9e3ef0, mvar=0xbf9e3ee8, iflag=0xbf9e3ee4, itran=0xbf9e3ee0, ierr=0xbf9e3edc)
at erpremc/trprfn.F:376
#13 0x460bcbda in erprop_ () at erdecks/erprop.F:62
#14 0x460bf33f in ertrch_ () at erdecks/ertrch.F:315
#15 0x460c0116 in ertrgo_ () at erdecks/ertrgo.F:236
#16 0x460bdda3 in ertrak_ (x1=0x9b145d8, p1=0x9b145e4, x2=0x9b1454c, p2=0x9b14558, ipa=0xbf9e4404, chopt=0xd00f2d8,
__g77_length_chopt=2) at erdecks/ertrak.F:211
#17 0x4617d44b in TGeant3::Ertrak (this=0xae7eec0, x1=0x9b145d8, p1=0x9b145e4, x2=0x9b1454c, p2=0x9b14558, ipa=6,
chopt=0xd00f2d8 "LE") at TGeant3/TGeant3.cxx:5392
#18 0x4511df1f in FairGeanePro::FindPCA (this=0x9b144c8, pca=2, PDGCode=13, point=
{<TObject> = {_vptr.TObject = 0x4324be08, fUniqueID = 0, fBits = 33554432, static fgDtorOnly = 0, static fgObjectStat = false, static fgIsA = 0x82b4838}, fX = 0, fY = 0, fZ = 0, static fgIsA = 0x8883f38},
wire1=
{<TObject> = {_vptr.TObject = 0x451352a8, fUniqueID = 0, fBits = 50331648, static fgDtorOnly = 0, static fgObjectStat = false, static fgIsA = 0x82b4838}, fX = 2.4562582328855371e-248, fY = 9.69071187023315
98e-262, fZ = 149.63770051287975, static fgIsA = 0x8883f38}, wire2=
{<TObject> = {_vptr.TObject = 0x4324be08, fUniqueID = 0, fBits = 33554432, static fgDtorOnly = 0, static fgObjectStat = false, static fgIsA = 0x82b4838}, fX = -24.684997098388671, fY = -19.100335878295898,
fZ = 79.9178, static fgIsA = 0x8883f38}, maxdistance=213.69470098481315, Rad=@0xbf9e4fe8, vpf=@0xbf9e53c0, vwi=@0xbf9e5350, Di=@0xbf9e4fe0,
trklength=@0xbf9e4fdc) at /d/panda02/broth/proot/5813/geane/FairGeanePro.cxx:553
#19 0x4438934d in GeaneTrackRep::extrapolateToLine (this=0xd901ef0, point1=@0xbf9e53f0, point2=@0xbf9e5390,
poca=@0xbf9e53c0, dirInPoca=@0xbf9e5310, poca_onwire=@0xbf9e5350)
at /d/panda02/broth/proot/5813/trackrep/GeaneTrackRep.cxx:343
#20 0x4432a738 in WirepointHitPolicy::detPlane (this=0xd9078e0, hit=0xd9075e0, rep=0xd901ef0)
at /d/panda02/broth/proot/5813/genfit/WirepointHitPolicy.cxx:91
#21 0x44a2996e in RecoHitIfc<WirepointHitPolicy>::getDetPlane (this=0xd9075e0, rep=0xd901ef0) at RecoHitIfc.h:65
#22 0x4430ca1c in Kalman::processHit (this=0xbf9e6500, hit=0xd9075e0, rep=0xd901ef0)
at /d/panda02/broth/proot/5813/genfit/Kalman.cxx:241
#23 0x4430bb18 in Kalman::fittingPass (this=0xbf9e6500, trk=0xd72f1f0, direction=1)
at /d/panda02/broth/proot/5813/genfit/Kalman.cxx:144
#24 0x4430b2b2 in Kalman::processTrack (this=0xbf9e6500, trk=0xd72f1f0)
at /d/panda02/broth/proot/5813/genfit/Kalman.cxx:36
#25 0x4509f3c5 in PndLheKalmanTask::Exec (this=0xb19c7d0, opt=0x43d53040 "")
at /d/panda02/broth/proot/5813/lhetrack/PndLheKalmanTask.cxx:247
#26 0x401b16af in TTask::ExecuteTasks (this=0x86c1a40, option=0x43d53040 "") at core/base/src/TTask.cxx:298
#27 0x401b14b1 in TTask::ExecuteTask (this=0x86c1a40, option=0x43d53040 "") at core/base/src/TTask.cxx:261
#28 0x43cdc46d in FairRunAna::Run (this=0x86c19b8, Ev_start=0, Ev_end=1000)
at /d/panda02/broth/proot/5813/base/FairRunAna.cxx:248
#29 0x43d1a746 in G__FairDict_532_0_5 (result7=0xbf9eda60, funcname=0x86bfa10 "\001", libp=0xbf9e7be0, hash=0)
at /d/panda02/broth/proot/5813/build/base/FairDict.cxx:9067
#30 0x407b6126 in Cint::G__ExceptionWrapper (funcp=0x43d1a642 <G__FairDict_532_0_5>, result7=0xbf9eda60,
funcname=0x86bfa10 "\001", libp=0xbf9e7be0, hash=0) at cint/cint/src/Api.cxx:364
#31 0x408757f5 in G__execute_call (result7=0xbf9eda60, libp=0xbf9e7be0, ifunc=0x86bfa10, ifn=0)
at cint/cint/src/newlink.cxx:2305
...
is in /d/panda02/broth/proot/5813/tutorials/analysis/data/4-kalman.log
But even if it runs fine, the tuples of PndChargedCandidates are all empty (in cases of running over 10 or 100 events).
Bernhard
|
|
|
Goto Forum:
Current Time: Tue Nov 26 06:15:28 CET 2024
Total time taken to generate the page: 0.00694 seconds
|