Home » PANDA » PandaRoot » Bugs, Fixes, Releases » STTPoints do not match with STTHelixHit
|
Re: STTPoints do not match with STTHelixHit [message #8283 is a reply to message #8281] |
Thu, 23 April 2009 12:07 |
Lia Lavezzi
Messages: 291 Registered: May 2007 Location: Torino
|
first-grade participant |
From: *pv.infn.it
|
|
Hi Tobias,
I tried to run your macros but I get some crash:
1) runMvdSim.C: I switched from FairEvtGenGenerator to FairBoxGenerator since I don' t have the data file from EvtGen but it crashes with this error
### Run 0 start.
FairMCApplication::GeneratePrimaries()
-I FairPrimaryGenerator: 10 primary tracks from vertex (0, 0, 0)Event Time = 0(ns)
>>> Event 0
*** Break *** segmentation violation
Using host libthread_db library "/lib/tls/libthread_db.so.1".
Attaching to program: /proc/7719/exe, process 7719
[Thread debugging using libthread_db enabled]
[New Thread -1208269120 (LWP 7719)]
0x00a8d7a2 in _dl_sysinfo_int80 () from /lib/ld-linux.so.2
#1 0x067270f3 in __waitpid_nocancel () from /lib/tls/libc.so.6
#2 0x066d07b9 in do_system () from /lib/tls/libc.so.6
#3 0x009f298d in system () from /lib/tls/libpthread.so.0
#4 0x0032c180 in TUnixSystem::Exec (this=0x8b14638,
#5 0x0032c634 in TUnixSystem::StackTrace (this=0x8b14638)
at core/unix/src/TUnixSystem.cxx:2121
#6 0x0032a479 in TUnixSystem::DispatchSignals (this=0x8b14638,
sig=kSigSegmentationViolation) at core/unix/src/TUnixSystem.cxx:1089
#7 0x00328315 in SigHandler (sig=kSigSegmentationViolation)
at core/unix/src/TUnixSystem.cxx:351
#8 0x0032f38c in sighandler (sig=11) at core/unix/src/TUnixSystem.cxx:3344
#9 <signal handler called>
#10 0x060c2a74 in TGeoVoxelFinder::GetValidExtra (this=0xac48e18,
array1=0xac4b395 "\224\226\222\022\020 ",
array2=0xac4b1d8 "\210�����\234\230\020", list=0xac4b6bc,
ncheck=@0xbff1b2c8) at geom/geom/src/TGeoVoxelFinder.cxx:605
#11 0x060c4791 in TGeoVoxelFinder::GetNextCandidates (this=0xac48e18,
point=0xbff1b320, ncheck=@0xbff1b2c8)
at geom/geom/src/TGeoVoxelFinder.cxx:1166
#12 0x060c57f7 in TGeoVoxelFinder::GetNextVoxel (this=0xac48e18,
point=0xbff1b320, ncheck=@0xbff1b2c8)
at geom/geom/src/TGeoVoxelFinder.cxx:1484
#13 0x0608d10a in TGeoShapeAssembly::DistFromOutside (this=0x95de2b0,
point=0xbff1b3e0, dir=0xbff1b3c0, iact=3, step=1411.954476830846, safe=0x0)
at geom/geom/src/TGeoShapeAssembly.cxx:230
#14 0x0608d030 in TGeoShapeAssembly::DistFromOutside (this=0x95df438,
point=0xbff1b4d0, dir=0xbff1b4b0, iact=3, step=1411.954476830846, safe=0x0)
at geom/geom/src/TGeoShapeAssembly.cxx:213
#15 0x0608d030 in TGeoShapeAssembly::DistFromOutside (this=0x961b5b8,
point=0xbff1b5c0, dir=0xbff1b5a0, iact=3, step=1411.954476830846, safe=0x0)
at geom/geom/src/TGeoShapeAssembly.cxx:213
#16 0x0608d030 in TGeoShapeAssembly::DistFromOutside (this=0x965ef80,
point=0xbff1b6c0, dir=0xbff1b6a0, iact=3, step=1411.954476830846, safe=0x0)
at geom/geom/src/TGeoShapeAssembly.cxx:213
#17 0x0605c574 in TGeoNavigator::FindNextDaughterBoundary (this=0x91820c8,
point=0xbff1b830, dir=0xbff1b810, idaughter=@0xbff1b808, compmatrix=true)
at geom/geom/src/TGeoNavigator.cxx:909
#18 0x0605b0b6 in TGeoNavigator::FindNextBoundary (this=0x91820c8,
stepmax=25697.528866946272, path=0x6f4c01a "", frombdr=true)
at geom/geom/src/TGeoNavigator.cxx:613
#19 0x0604042c in TGeoManager::FindNextBoundary (this=0x8f167c0,
stepmax=-25697.528866946272, path=0x6f4c01a "", frombdr=true)
at geom/geom/src/TGeoManager.cxx:2295
#20 0x06f43dc9 in TG4RootNavigator::ComputeStep (this=0xc65c830,
pGlobalPoint=@0xbff1bb60, pDirection=@0xbff1bb80,
pCurrentProposedStepLength=256975.28866946368, pNewSafety=@0xbff1bbb0)
at montecarlo/g4root/src/TG4RootNavigator.cxx:143
#21 0x03876012 in G4Transportation::AlongStepGetPhysicalInteractionLength (
this=0xe890d00, track=@0xf38e100, currentMinimumStep=256975.28866946368,
currentSafety=@0xbff1bc50, selection=0xc663824)
at src/G4Transportation.cc:227
#22 0x05ab362c in G4VProcess::AlongStepGPIL (this=0xe890d00, track=@0xf38e100,
previousStepSize=32.659569449450601, currentMinimumStep=256975.28866946368,
proposedSafety=@0xbff1bc50, selection=0xc663824)
at /home/lavezzi/fairsoft2/fairsoft/transport/geant4/source/processes/management/include/G4VProcess.hh:425
#23 0x05ab235f in G4SteppingManager::DefinePhysicalStepLength (this=0xc663718)
at src/G4SteppingManager2.cc:218
#24 0x05ab6c74 in G4SteppingManager::Stepping (this=0xc663718)
at src/G4SteppingManager.cc:181
#25 0x05ac05ea in G4TrackingManager::ProcessOneTrack (this=0xc6636b0,
apValueG4Track=0xf38e100) at src/G4TrackingManager.cc:126
#26 0x05e3ef3b in G4EventManager::DoProcessing (this=0xc663668,
anEvent=0xc66f6a8) at src/G4EventManager.cc:185
#27 0x05e3f72b in G4EventManager::ProcessOneEvent (this=0xc663668,
anEvent=0xc66f6a8) at src/G4EventManager.cc:335
#28 0x05d9f66f in G4RunManager::DoEventLoop (this=0xc6635b8, n_event=10,
macroFile=0x0, n_select=-1) at src/G4RunManager.cc:235
#29 0x05d9ef99 in G4RunManager::BeamOn (this=0xc6635b8, n_event=10,
macroFile=0x0, n_select=-1) at src/G4RunManager.cc:140
#30 0x07cb3925 in TG4RunManager::ProcessRun (this=0x93fd4c0, nofEvents=10)
at run/src/TG4RunManager.cxx:388
#31 0x07cba5f2 in TGeant4::ProcessRun (this=0x93f8e30, nofEvents=10)
at run/src/TGeant4.cxx:1140
#32 0x07dafec4 in FairMCApplication::RunMC (this=0x9189258, nofEvents=10)
at /home/lavezzi/test_dev/geanefor/trunk/base/FairMCApplication.cxx:228
#33 0x07dc91f6 in FairRunSim::Run (this=0x9103990, NStart=10, NStop=0)
at /home/lavezzi/test_dev/geanefor/trunk/base/FairRunSim.cxx:289
#34 0x07dfebf5 in G__FairDict_532_0_5 (result7=0xbff23300,
funcname=0x91019a0 "\001", libp=0xbff1d480, hash=0)
at /home/lavezzi/test_dev/geanefor/trunk/cbuild/base/FairDict.cxx:9292
#35 0x00b0bc0d in Cint::G__ExceptionWrapper (
funcp=0x7dfeade <G__FairDict_532_0_5>, result7=0xbff23300,
funcname=0x91019a0 "\001", libp=0xbff1d480, hash=0)
at cint/cint/src/Api.cxx:364
#36 0x00bb2ddd in G__execute_call (result7=0xbff23300, libp=0xbff1d480,
ifunc=0x91019a0, ifn=0) at cint/cint/src/newlink.cxx:2305
#37 0x00bb3422 in G__call_cppfunc (result7=0xbff23300, libp=0xbff1d480,
ifunc=0x91019a0, ifn=0) at cint/cint/src/newlink.cxx:2471
#38 0x00b951b0 in G__interpret_func (result7=0xbff23300,
funcname=0xbff22f00 "Run", libp=0xbff1d480, hash=309, p_ifunc=0x91019a0,
funcmatch=1, memfunc_flag=1) at cint/cint/src/ifunc.cxx:5245
#39 0x00b7bc28 in G__getfunction (item=0xbff26146 "Run(nEvents)",
known3=0xbff2555c, memfunc_flag=1) at cint/cint/src/func.cxx:2534
#40 0x00c691cf in G__getstructmem (store_var_type=112,
varname=0xbff23550 "p5��d~�", membername=0xbff26146 "Run(nEvents)",
tagname=0xbff23760 "fRun", known2=0xbff2555c, varglobal=0xd06ca0, objptr=2)
at cint/cint/src/var.cxx:6623
#41 0x00c5ccb6 in G__getvariable (item=0xbff26140 "fRun->Run(nEvents)",
known=0xbff2555c, varglobal=0xd06ca0, varlocal=0x0)
at cint/cint/src/var.cxx:5252
#42 0x00b6f873 in G__getitem (item=0xbff26140 "fRun->Run(nEvents)")
at cint/cint/src/expr.cxx:1884
#43 0x00b6dee2 in G__getexpr (expression=0xbff27a30 "fRun->Run(nEvents)")
at cint/cint/src/expr.cxx:1470
#44 0x00bda2ca in G__exec_function (statement=0xbff27a30 "fRun->Run(nEvents)",
pc=0xbff27e5c, piout=0xbff27e54, plargestep=0xbff27e44, presult=0xbff27a00)
at cint/cint/src/parse.cxx:601
#45 0x00be636a in G__exec_statement (mparen=0xbff27ec0)
at cint/cint/src/parse.cxx:6972
#46 0x00b4bf98 in G__exec_tempfile_core (
file=0xbff30cc0 "/home/lavezzi/test_dev/geanefor/trunk/macro/mvd/Tobias/./runMvdSim.C", fp=0x0) at cint/cint/src/debug.cxx:251
#47 0x00b4d4fc in G__exec_tempfile (
file=0xbff30cc0 "/home/lavezzi/test_dev/geanefor/trunk/macro/mvd/Tobias/./runMvdSim.C") at cint/cint/src/debug.cxx:798
#48 0x00bf1238 in G__process_cmd (line=0x717552 "cbmroot.Drc.DrcBarPoint_]",
prompt=0x8b177ec "", more=0x8b177e4, err=0xbff3152c, rslt=0xbff31530)
at cint/cint/src/pause.cxx:3074
#49 0x00317724 in TCint::ProcessLine (this=0x8b177c8,
line=0x717552 "cbmroot.Drc.DrcBarPoint_]", error=0xbff33d24)
at core/meta/src/TCint.cxx:339
#50 0x00317af9 in TCint::ProcessLineSynch (this=0x8b177c8,
line=0x717552 "cbmroot.Drc.DrcBarPoint_]", error=0xbff33d24)
at core/meta/src/TCint.cxx:406
#51 0x0022b879 in TApplication::ExecuteFile (file=0xbff31cc3 "runMvdSim.C",
error=0xbff33d24) at core/base/src/TApplication.cxx:935
#52 0x0022b129 in TApplication::ProcessFile (this=0x8b38570,
file=0xbff31cc3 "runMvdSim.C", error=0xbff33d24)
at core/base/src/TApplication.cxx:825
#53 0x0022b065 in TApplication::ProcessLine (this=0x8b38570,
line=0xbff31cc0 ".x runMvdSim.C", sync=false, err=0xbff33d24)
at core/base/src/TApplication.cxx:798
#54 0x008c3df6 in TRint::Run (this=0x8b38570, retrn=false)
at core/rint/src/TRint.cxx:355
#55 0x08048df7 in main (argc=1, argv=0xbff33e04) at main/src/rmain.cxx:29
If I change the solenoid geo file from PandaSolenoidV833.root to FullSolenoid.root (I copied this from tutorials/lhetrack/run_sim_sttcombi_pgun.C) it goes to the end without problems... I don' t exactly know why;
2) runMvdDigi.C: it crashes with this message
-I- PndMvdNoiseProducer: Intialisation successfull
Warning in <TGeoManager::Init>: Deleting previous geometry: FAIRGeom/FAIR geometry
Info in <TGeoManager::CloseGeometry>: Geometry loaded from file...
Info in <TGeoManager::SetTopVolume>: Top volume is cave. Master volume is cave
Info in <TGeoManager::Voxelize>: Voxelizing...
Info in <TGeoManager::CloseGeometry>: 473268 nodes/ 3406 volume UID's in FAIR geometry
Info in <TGeoManager::CloseGeometry>: ----------------modeler ready----------------
*** Break *** segmentation violation
Using host libthread_db library "/lib/tls/libthread_db.so.1".
Attaching to program: /proc/7365/exe, process 7365
[Thread debugging using libthread_db enabled]
[New Thread -1208207680 (LWP 7365)]
0x00a8d7a2 in _dl_sysinfo_int80 () from /lib/ld-linux.so.2
#1 0x023ed0f3 in __waitpid_nocancel () from /lib/tls/libc.so.6
#2 0x023967b9 in do_system () from /lib/tls/libc.so.6
#3 0x009f298d in system () from /lib/tls/libpthread.so.0
#4 0x0032c180 in TUnixSystem::Exec (this=0x9b05638,
at core/unix/src/TUnixSystem.cxx:1941
#5 0x0032c634 in TUnixSystem::StackTrace (this=0x9b05638) at core/unix/src/TUnixSystem.cxx:2121
#6 0x0032a479 in TUnixSystem::DispatchSignals (this=0x9b05638, sig=kSigSegmentationViolation)
at core/unix/src/TUnixSystem.cxx:1089
#7 0x00328315 in SigHandler (sig=kSigSegmentationViolation) at core/unix/src/TUnixSystem.cxx:351
#8 0x0032f38c in sighandler (sig=11) at core/unix/src/TUnixSystem.cxx:3344
#9 <signal handler called>
#10 0x532e7475 in ?? ()
#11 0x06d1ba8a in TGeoManager::GetVolume ()
from /home/lavezzi/fairsoft2/fairsoft/tools/root_v5.20.00/lib/libGeom.so
#12 0x03be6924 in PndMvdGeoHandling::GetPath (this=0xef6e728, id=
{_vptr.TString = 0x6f53c8, static fgInitialCapac = 15, static fgResizeInc = 16, static fgFreeboard = 15, fData = 0xccc4ae8 "1_1/84_0/165_1/164_2/163_1/146_4/143_1/142_1/", static fgIsA = 0xa199168})
at /home/lavezzi/test_dev/geanefor/trunk/mvd/MvdTools/PndMvdGeoHandling.cxx:90
#13 0x03bd337b in PndMvdStripHitProducer::SelectSensorParams (this=0xa372e08, detname=Cannot access memory at address 0x532e7475
)
at /home/lavezzi/test_dev/geanefor/trunk/mvd/MvdDigi/PndMvdStripHitProducer.cxx:415
#14 0x03bd257a in PndMvdStripHitProducer::Exec (this=0xa372e08, opt=0x21a6c18 "")
at /home/lavezzi/test_dev/geanefor/trunk/mvd/MvdDigi/PndMvdStripHitProducer.cxx:224
#15 0x002a99de in TTask::ExecuteTasks (this=0xa0f3208, option=0x21a6c18 "") at core/base/src/TTask.cxx:298#16 0x002a97d9 in TTask::ExecuteTask (this=0xa0f3208, option=0x21a6c18 "") at core/base/src/TTask.cxx:261
#17 0x021455f6 in FairRunAna::Run (this=0xa10ff88, Ev_start=0, Ev_end=5)
at /home/lavezzi/test_dev/geanefor/trunk/base/FairRunAna.cxx:248
#18 0x0217fb92 in G__FairDict_532_0_5 (result7=0xbfef9630, funcname=0xa10dff8 "\001", libp=0xbfef37b0,
hash=0) at /home/lavezzi/test_dev/geanefor/trunk/cbuild/base/FairDict.cxx:9288
#19 0x00b0bc0d in Cint::G__ExceptionWrapper (funcp=0x217fade <G__FairDict_532_0_5>, result7=0xbfef9630,
funcname=0xa10dff8 "\001", libp=0xbfef37b0, hash=0) at cint/cint/src/Api.cxx:364
#20 0x00bb2ddd in G__execute_call (result7=0xbfef9630, libp=0xbfef37b0, ifunc=0xa10dff8, ifn=0)
at cint/cint/src/newlink.cxx:2305
#21 0x00bb3422 in G__call_cppfunc (result7=0xbfef9630, libp=0xbfef37b0, ifunc=0xa10dff8, ifn=0)
at cint/cint/src/newlink.cxx:2471
#22 0x00b951b0 in G__interpret_func (result7=0xbfef9630, funcname=0xbfef9230 "Run", libp=0xbfef37b0,
hash=309, p_ifunc=0xa10dff8, funcmatch=1, memfunc_flag=1) at cint/cint/src/ifunc.cxx:5245
#23 0x00b7bc28 in G__getfunction (item=0xbfefc476 "Run(nStart,nEvents)", known3=0xbfefb88c, memfunc_flag=1) at cint/cint/src/func.cxx:2534
#24 0x00c691cf in G__getstructmem (store_var_type=112, varname=0xbfef9880 "�\230��d~�",
membername=0xbfefc476 "Run(nStart,nEvents)", tagname=0xbfef9a90 "fRun", known2=0xbfefb88c,
varglobal=0xd06ca0, objptr=2) at cint/cint/src/var.cxx:6623
#25 0x00c5ccb6 in G__getvariable (item=0xbfefc470 "fRun->Run(nStart,nEvents)", known=0xbfefb88c,
varglobal=0xd06ca0, varlocal=0x0) at cint/cint/src/var.cxx:5252
#26 0x00b6f873 in G__getitem (item=0xbfefc470 "fRun->Run(nStart,nEvents)") at cint/cint/src/expr.cxx:1884
#27 0x00b6dee2 in G__getexpr (expression=0xbfefdd60 "fRun->Run(nStart,nEvents)")
at cint/cint/src/expr.cxx:1470
#28 0x00bda2ca in G__exec_function (statement=0xbfefdd60 "fRun->Run(nStart,nEvents)", pc=0xbfefe18c,
piout=0xbfefe184, plargestep=0xbfefe174, presult=0xbfefdd30) at cint/cint/src/parse.cxx:601
#29 0x00be636a in G__exec_statement (mparen=0xbfefe1f0) at cint/cint/src/parse.cxx:6972
#30 0x00b4bf98 in G__exec_tempfile_core (
file=0xbff06ff0 "/home/lavezzi/test_dev/geanefor/trunk/macro/mvd/Tobias/./runMvdDigi.C", fp=0x0)
at cint/cint/src/debug.cxx:251
#31 0x00b4d4fc in G__exec_tempfile (
file=0xbff06ff0 "/home/lavezzi/test_dev/geanefor/trunk/macro/mvd/Tobias/./runMvdDigi.C")
at cint/cint/src/debug.cxx:798
#32 0x00bf1238 in G__process_cmd (line=0x717553 "Info_]", prompt=0x9b087ec "", more=0x9b087e4,
err=0xbff0785c, rslt=0xbff07860) at cint/cint/src/pause.cxx:3074
#33 0x00317724 in TCint::ProcessLine (this=0x9b087c8, line=0x717553 "Info_]", error=0xbff0a054)
at core/meta/src/TCint.cxx:339
#34 0x00317af9 in TCint::ProcessLineSynch (this=0x9b087c8, line=0x717553 "Info_]", error=0xbff0a054)
at core/meta/src/TCint.cxx:406
#35 0x0022b879 in TApplication::ExecuteFile (file=0xbff07ff3 "runMvdDigi.C", error=0xbff0a054)
at core/base/src/TApplication.cxx:935
#36 0x0022b129 in TApplication::ProcessFile (this=0x9b29570, file=0xbff07ff3 "runMvdDigi.C",
error=0xbff0a054) at core/base/src/TApplication.cxx:825
#37 0x0022b065 in TApplication::ProcessLine (this=0x9b29570, line=0xbff07ff0 ".x runMvdDigi.C", sync=false, err=0xbff0a054) at core/base/src/TApplication.cxx:798
#38 0x008c3df6 in TRint::Run (this=0x9b29570, retrn=false) at core/rint/src/TRint.cxx:355
#39 0x08048df7 in main (argc=1, argv=0xbff0a134) at main/src/rmain.cxx:29
So I could not check the SttHelixHit myself in this case , can you please run your reconstruction adding:
sttHHProducer->WriteHistograms()
just after the run has finished (like in macro/stt/runreco.C) to fill the residuals histograms and post these figures? It should be more evident there how big is the difference...
Thank you and ciao,
Lia.
|
|
|
Re: STTPoints do not match with STTHelixHit [message #8286 is a reply to message #8283] |
Thu, 23 April 2009 16:22 |
StefanoSpataro
Messages: 2736 Registered: June 2005 Location: Torino
|
first-grade participant |
From: *to.infn.it
|
|
Hi,
I have tried to run run_sim_sttcombi_pgun.C (modified to have 3 pions per event and with the momentum range suggested by Tobias), 10 events, and then run_rigi_sttcombi.C.
Just plotting from points.root:
cbmsim.Draw("STTPoint.GetXtot():STTPoint.GetZtot()","","",1,nevt)
from digi.root:
cbmsim.Draw("SttHelixHit.fY:SttHelixHit.fZ","","",1,nevt)
(where nevt is the number of the event that I want to analyse) I can see the same behaviour: it seems that there are still some problems with the z coordinate for low momentum tracks in helixhit.
|
|
|
Re: STTPoints do not match with STTHelixHit [message #8294 is a reply to message #8286] |
Fri, 24 April 2009 13:31 |
Susanna Costanza
Messages: 33 Registered: January 2008 Location: Pavia
|
continuous participant |
From: *pv.infn.it
|
|
Hi Tobias and Stefano,
we made some tests with 3 pions per event and with 3 muons per event, splitting the momentum range in 2 ranges: [0.1,1.] GeV/c and [1.,2.5] GeV/c.
Here is what we see for muons:
--> [0.1,1]:
- for 67% tracks, STTPoints match STTHelixHits
- for 17% tracks, failure in xy fit (no z reconstruction)
- for 16% tracks, all fits are successful but STTPoints do not match STTHelixHits
--> [1.,2.5]:
- for 83.3% tracks, STTPoints match STTHelixHits
- for 16.7% tracks, all fits are successful but STTPoints do not match STTHelixHits.
It's clear that the number of mismatches is higher for lower momenta, due to fit failures.
Concerning pions, the problem is that there are low momentum tracks (electrons) also in the range [1.,2.5] GeV/c, so the failures are more evident.
At present, this is the situation; z reconstruction is under investigation.
Ciao,
Susanna & Lia
|
|
|
|
|
|
|
Goto Forum:
Current Time: Sun Nov 10 03:09:06 CET 2024
Total time taken to generate the page: 0.00908 seconds
|