donghee@donghee-laptop:~/GSI/pandaroot/macro/dvcs$ gdb root.exe GNU gdb 6.8-debian Copyright (C) 2008 Free Software Foundation, Inc. License GPLv3+: GNU GPL version 3 or later This is free software: you are free to change and redistribute it. There is NO WARRANTY, to the extent permitted by law. Type "show copying" and "show warranty" for details. This GDB was configured as "i486-linux-gnu"... (gdb) run -b -q My_dvcs_full_run.C Starting program: /home/donghee/GSI/fairsoft/tools/root_v5.20.00/bin/root.exe -b -q My_dvcs_full_run.C [Thread debugging using libthread_db enabled] [New Thread 0xb6c316c0 (LWP 14510)] ******************************************* * * * W E L C O M E to R O O T * * * * Version 5.20/00 24 June 2008 * * * * You are welcome to visit our Web site * * http://root.cern.ch * * * ******************************************* ROOT 5.20/00 (trunk@24524, Mar 27 2009, 22:10:00 on linux) CINT/ROOT C/C++ Interpreter version 5.16.29, Jan 08, 2008 Type ? for help. Commands must be C++ statements. Enclose multiple statements between { }. root [0] Processing My_dvcs_full_run.C... - RTDB container factory FairBaseContFact - RTDB container factory PndFieldContFact - RTDB container factory PndPassiveContFact PSaid instance created... access via gSaid->f() - RTDB container factory PndEmcContFact - RTDB container factory PndTpcContFact Warning in : class PndSttTrack already in TClassTable - RTDB container factory PndSttContFact Warning in : class PndMvdPixel already in TClassTable Warning in : class PndMvdStrip already in TClassTable Warning in : class PndMvdMCPoint already in TClassTable Warning in : class PndMvdDigi already in TClassTable Warning in : class PndMvdDigiPixel already in TClassTable Warning in : class PndMvdApvHit already in TClassTable - RTDB container factory PndMvdContFact Warning in : class PndMvdPidCand already in TClassTable - RTDB container factory PndTofContFact - RTDB container factory PndDrcContFact - RTDB container factory PndMdtContFact - RTDB container factory PndDchContFact - RTDB container factory PndLheContFact -I- FairRun::SetMaterials() Media file used: /home/donghee/GSI/pandaroot/geometry/media_pnd.geo MDT version used: torino -I PndPythiaGenerator: Opening input file /home/donghee/GSI/pandaroot/macro/dvcs/hepevt/Panda_event_1.txt ============== FairRunSim: Initialising simulation run ============== Info in : Geometry FAIRGeom, FAIR geometry created -I- FairGeoMedia Read media -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/TransMap.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/DipoleMap.root -I- PndFieldMap: Reading field map from ROOT file /home/donghee/GSI/pandaroot/input/SolenoidMap.root Loading Geant4 granular libraries ... Loading VGM libraries ... Loading libraries ... finished Info in : Top volume is cave. Master volume is cave Material aluminium is not defined Create Medium aluminium Material iron is not defined Create Medium iron ====== EMC:: ConstructASCIIGeometry() ====== ============================================= PndEmcReader: EMC geometry file == /home/donghee/GSI/pandaroot/geometry/emc_module1245.dat Emc module = 1 ******** Emc module = 2 ******** Emc module = 3 ******** Emc module = 4 ******** Emc module = 5 ******** Material silicon is not defined Create Medium silicon Material carbon is not defined Create Medium carbon -I container name PndGeoSttPar Material DCHmixture2bar is not defined Create Medium DCHmixture2bar Material tungsten is not defined Create Medium tungsten Material Aerogel is not defined Create Medium Aerogel Material C4F10 is not defined Create Medium C4F10 FairMCApplication::ConstructGeometry() : Now closing the geometry Info in : Fixing runtime shapes... Info in : ...Nothing to fix Info in : Counting nodes... Info in : Voxelizing... Info in : Building cache... Info in : --- Maximum geometry depth set to 100 Info in : 476454 nodes/ 3450 volume UID's in FAIR geometry Info in : ----------------modeler ready---------------- Info in : TG4RootNavigator created and registered to G4TransportationManager Running TVirtualMCApplication::ConstructGeometry ************************************************************* Geant4 version Name: geant4-09-01-patch-02 (9-May-2008) Copyright : Geant4 Collaboration Reference : NIM A 506 (2003), 250-303 WWW : http://cern.ch/geant4 ************************************************************* Info in : Creating G4 hierarchy ... Info in : Converting reflections in: FAIRGeom - FAIR geometry ... Info in : Done ===> GEANT4 materials created and mapped to TGeo ones... ===> GEANT4 physical volumes created and mapped to TGeo hierarchy... ### INFO: TG4RootDetectorConstruction::Construct() finished TG4PostDetConstruction::Initialize G4 Stat: instantiated 4353 logical volumes 81688 physical volumes Info in : ROOT detector construction class connected to G4RunManager Adding HadronPhysicsList QGSP_BERT_EMV <<< Geant4 Physics List engine packaging library: PACK 5.4 <<< Geant4 Physics List simulation engine: QGSP_BERT_EMV 1.0 Adding SpecialPhysicsList stepLimiter+specialCuts+specialControls Debug mode is switched on. Visualization Manager instantiating... Visualization Manager initialising... Registering graphics systems... You have successfully registered the following graphics systems. Current available graphics systems are: ASCIITree (ATree) DAWNFILE (DAWNFILE) G4HepRepFile (HepRepFile) G4HepRep (HepRepXML) RayTracer (RayTracer) VRML1FILE (VRML1FILE) VRML2FILE (VRML2FILE) Registering model factories... You have successfully registered the following model factories. Registered model factories: drawByCharge drawByParticleID Registered filter factories: None Geant4 has been created. -I g4Config() using g4conf macro: /home/donghee/GSI/pandaroot/gconfig/g4config.in SetCuts Macro: Setting Processes.. SetCuts Macro: Setting cuts.. ### Adding Neutron tracking cut for neutron ### cut value is 10 microseconds ### Hadron physics constructed. ### Processes mapped to VMC controls ok. ### Step limiter physics constructed. ### Special Cuts constructed. ### Special Controls constructed. ### User particles physics constructed. ### Processes mapped to VMC codes ok. -I- Initializing PndMvdDetector() ************************************************************* initialisation for run id 524362640 ************************************************************* -I FairParSet::init() 2 0 Error in : FairBaseParSet not initialized -I FairParSet::init() 2 0 Error in : PndGeoPassivePar not initialized -I FairParSet::init() 2 0 Error in : PndGeoSttPar not initialized -I FairParSet::init() 2 0 Error in : PndGeoTofPar not initialized -I FairParSet::init() 2 0 Error in : PndGeoMdtPar not initialized Error in : Error occured during initialization -I- FairMCApplication -> simulation RunID: 524362640 ********************************************** GEANT4 Geometry statistics: 4353 logical volumes 81688 physical volumes 21 materials 20 user limits 3448 sensitive detectors ********************************************** -I- FairMCApplication:: Monte carlo Engine Initialisation with TGeant4 create PndFieldPar container PndMultiFieldPar create PndFieldPar container PndTransPar create PndFieldPar container PndDipolePar create PndFieldPar container PndSolenoidPar RuntimeDb: write container FairBaseParSet *** FairBaseParSet written to ROOT file version: 36 RuntimeDb: write container PndGeoPassivePar *** PndGeoPassivePar written to ROOT file version: 36 RuntimeDb: write container PndGeoSttPar *** PndGeoSttPar written to ROOT file version: 32 RuntimeDb: write container PndGeoTofPar *** PndGeoTofPar written to ROOT file version: 36 RuntimeDb: write container PndGeoMdtPar RuntimeDb: write container PndMultiFieldPar *** PndMultiFieldPar written to ROOT file version: 37 RuntimeDb: write container PndTransPar RuntimeDb: write container PndDipolePar RuntimeDb: write container PndSolenoidPar -------------------------------------------------------------------------------- -------------- actual containers in runtime database ------------------------- FairBaseParSet Test class for parameter io PndGeoPassivePar Passive Geometry Parameters PndGeoSttPar Stt Geometry Parameters PndGeoTofPar tof Geometry Parameters PndGeoMdtPar MDT Geometry Parameters PndMultiFieldPar Multiple Field parameter container PndTransPar Trans. Field parameter container PndDipolePar Dipole Field parameter container PndSolenoidPar Solenoid Field parameter container -------------- runs, versions ------------------------------------------------ run id container 1st-inp 2nd-inp output run: 524362640 FairBaseParSet 524362640 -1 36 PndGeoPassivePar 524362640 -1 36 PndGeoSttPar 524362640 -1 32 PndGeoTofPar 524362640 -1 36 PndGeoMdtPar -1 -1 0 PndMultiFieldPar 524362640 -1 37 PndTransPar 524362640 -1 0 PndDipolePar 524362640 -1 0 PndSolenoidPar 524362640 -1 0 -------------- input/output -------------------------------------------------- first input: none second input: none output: OBJ: FairParRootFile /home/donghee/GSI/pandaroot/macro/dvcs/hepevt/Panda_event_1.full.param.root : 0 at: 0xe286bc8 Root file I/O /home/donghee/GSI/pandaroot/macro/dvcs/hepevt/Panda_event_1.full.param.root is open detector I/Os: FairGenericParIo phot: Total cross sections from Sandia parametrisation. Sampling according PhotoElectric model compt: Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV Sampling according Klein-Nishina model tables are built for gamma Lambda tables from 100 eV to 100 GeV in 90 bins. conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z; sampling secondary e+e- according Bethe-Heitler model tables are built for gamma Lambda tables from 1.022 MeV to 100 GeV in 100 bins. msc: Model variant of multiple scattering for e- Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 0 eIoni: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.8, integral: 1, fluct: 1 eBrem: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1 eIoni: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.8, integral: 1, fluct: 1 eBrem: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1 annihil: Sampling according eplus2gg model tables are built for e+ Lambda tables from 100 eV to 100 TeV in 120 bins. msc: Model variant of multiple scattering for proton Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 0 hIoni: tables are built for proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. NuclearStopping= 1 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 msc: Model variant of multiple scattering for GenericIon LateralDisplacementFlag= 0 Skin= 0 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 1 ionIoni: tables are built for GenericIon dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from BraggIon for protons below. NuclearStopping= 1 Stopping Power data for 8 ion/material pairs are used. Step function: finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1 hIoni: tables are built for anti_proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. NuclearStopping= 1 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 msc: Model variant of multiple scattering for mu+ Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 0 muIoni: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 muBrems: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model muPairProd: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model muIoni: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 muBrems: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model muPairProd: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model G4UHadronElasticProcess for neutron PDGcode= 2112 Elow(MeV)= 19 Elowest(eV)= 0 hIoni: tables are built for pi+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV Parametrisation from Bragg for protons below. NuclearStopping= 1 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 msc: Model variant of multiple scattering for pi- Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 0 hIoni: tables are built for pi- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV Parametrisation from Bragg for protons below. NuclearStopping= 1 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : air Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes recalculation needed : No Material : steel Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 20.8323 keV e- 1.31192 MeV e+ 1.23361 MeV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : vacuum Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV Region(s) which use this couple : DefaultRegionForTheWorld Index : 3 used in the geometry : Yes recalculation needed : No Material : Al+Be Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 4.45676 keV e- 496.074 keV e+ 478.087 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 4 used in the geometry : Yes recalculation needed : No Material : aluminium Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 6.88731 keV e- 596.68 keV e+ 568.011 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 5 used in the geometry : Yes recalculation needed : No Material : iron Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 20.8323 keV e- 1.28002 MeV e+ 1.21851 MeV Region(s) which use this couple : DefaultRegionForTheWorld Index : 6 used in the geometry : Yes recalculation needed : No Material : PWO Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 84.7768 keV e- 1.13176 MeV e+ 1.06419 MeV Region(s) which use this couple : DefaultRegionForTheWorld Index : 7 used in the geometry : Yes recalculation needed : No Material : lead Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 100.511 keV e- 1.37814 MeV e+ 1.28002 MeV Region(s) which use this couple : DefaultRegionForTheWorld Index : 8 used in the geometry : Yes recalculation needed : No Material : FscScint Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 2.36895 keV e- 355.791 keV e+ 347.138 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 9 used in the geometry : Yes recalculation needed : No Material : silicon Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 6.88731 keV e- 540.718 keV e+ 521.113 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 10 used in the geometry : Yes recalculation needed : No Material : carbon Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 3.29462 keV e- 568.011 keV e+ 554.196 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 11 used in the geometry : Yes recalculation needed : No Material : mylar Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 2.97898 keV e- 417.539 keV e+ 402.4 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 12 used in the geometry : Yes recalculation needed : No Material : argon Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV Region(s) which use this couple : DefaultRegionForTheWorld Index : 13 used in the geometry : Yes recalculation needed : No Material : copper Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 24.6072 keV e- 1.39521 MeV e+ 1.31192 MeV Region(s) which use this couple : DefaultRegionForTheWorld Index : 14 used in the geometry : Yes recalculation needed : No Material : polyvinyltoluene Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 2.36895 keV e- 355.791 keV e+ 347.138 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 15 used in the geometry : Yes recalculation needed : No Material : DCHmixture2bar Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 990 eV e- 1.79809 keV e+ 1.76756 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 16 used in the geometry : Yes recalculation needed : No Material : tungsten Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 105.401 keV e- 2.31122 MeV e+ 2.14666 MeV Region(s) which use this couple : DefaultRegionForTheWorld Index : 17 used in the geometry : Yes recalculation needed : No Material : C4F10 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 990 eV e- 11.6998 keV e+ 11.6998 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 18 used in the geometry : Yes recalculation needed : No Material : Aerogel Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 1.21271 keV e- 80.2205 keV e+ 79.239 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 19 used in the geometry : Yes recalculation needed : No Material : ArCO2 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 5.18033 keV e- 417.539 keV e+ 402.4 keV Region(s) which use this couple : DefaultRegionForTheWorld ==================================================================== ### Run 0 start. FairMCApplication::GeneratePrimaries() -I FairPrimaryGenerator: 3 primary tracks from vertex (0, 0, 0)Event Time = 0(ns) >>> Event 0 -I- PndStack: Filling MCTrack array... -I- PndStack: Number of primaries = 3 Total number of particles = 6078 Number of tracks in output = 171 -I- PndStack: Updating track indizes......stack and 7 collections updated. -I- PndStt: 77 points registered in this event. FairMCApplication::GeneratePrimaries() -I FairPrimaryGenerator: 3 primary tracks from vertex (0, 0, 0)Event Time = 0(ns) >>> Event 1 -I- PndStack: Filling MCTrack array... -I- PndStack: Number of primaries = 3 Total number of particles = 7306 Number of tracks in output = 66 -I- PndStack: Updating track indizes......stack and 7 collections updated. -I- PndStt: 26 points registered in this event. FairMCApplication::GeneratePrimaries() -I FairPrimaryGenerator: 3 primary tracks from vertex (0, 0, 0)Event Time = 0(ns) >>> Event 2 -I- PndStack: Filling MCTrack array... -I- PndStack: Number of primaries = 3 Total number of particles = 11623 Number of tracks in output = 63 -I- PndStack: Updating track indizes......stack and 7 collections updated. -I- PndStt: 145 points registered in this event. FairMCApplication::GeneratePrimaries() -I FairPrimaryGenerator: 3 primary tracks from vertex (0, 0, 0)Event Time = 0(ns) >>> Event 3 Info in : Cannot compute distance from inside the assembly (but from a component) -I- PndStack: Filling MCTrack array... -I- PndStack: Number of primaries = 3 Total number of particles = 6347 Number of tracks in output = 42 -I- PndStack: Updating track indizes......stack and 7 collections updated. -I- PndStt: 12 points registered in this event. FairMCApplication::GeneratePrimaries() -I FairPrimaryGenerator: 3 primary tracks from vertex (0, 0, 0)Event Time = 0(ns) >>> Event 4 -I- PndStack: Filling MCTrack array... -I- PndStack: Number of primaries = 3 Total number of particles = 8189 Number of tracks in output = 14 -I- PndStack: Updating track indizes......stack and 7 collections updated. -I- PndStt: 26 points registered in this event. FairMCApplication::GeneratePrimaries() -I FairPrimaryGenerator: 3 primary tracks from vertex (0, 0, 0)Event Time = 0(ns) >>> Event 5 -I- PndStack: Filling MCTrack array... -I- PndStack: Number of primaries = 3 Total number of particles = 5245 Number of tracks in output = 105 -I- PndStack: Updating track indizes......stack and 7 collections updated. -I- PndStt: 25 points registered in this event. FairMCApplication::GeneratePrimaries() -I FairPrimaryGenerator: 3 primary tracks from vertex (0, 0, 0)Event Time = 0(ns) >>> Event 6 -I- PndStack: Filling MCTrack array... -I- PndStack: Number of primaries = 3 Total number of particles = 3722 Number of tracks in output = 6 -I- PndStack: Updating track indizes......stack and 7 collections updated. -I- PndStt: 21 points registered in this event. FairMCApplication::GeneratePrimaries() -I FairPrimaryGenerator: 3 primary tracks from vertex (0, 0, 0)Event Time = 0(ns) >>> Event 7 -I- PndStack: Filling MCTrack array... -I- PndStack: Number of primaries = 3 Total number of particles = 11272 Number of tracks in output = 23 -I- PndStack: Updating track indizes......stack and 7 collections updated. -I- PndStt: 25 points registered in this event. FairMCApplication::GeneratePrimaries() -I FairPrimaryGenerator: 3 primary tracks from vertex (0, 0, 0)Event Time = 0(ns) >>> Event 8 -I- PndStack: Filling MCTrack array... -I- PndStack: Number of primaries = 3 Total number of particles = 5712 Number of tracks in output = 12 -I- PndStack: Updating track indizes......stack and 7 collections updated. -I- PndStt: 25 points registered in this event. FairMCApplication::GeneratePrimaries() -I FairPrimaryGenerator: 3 primary tracks from vertex (0, 0, 0)Event Time = 0(ns) >>> Event 9 Program received signal SIGFPE, Arithmetic exception. [Switching to Thread 0xb6c316c0 (LWP 14510)] G4QElasticCrossSection::GetTabValues (this=0xaf1f34a0, lp=1.826771653543307, PDG=2212, tgZ=82, tgN=125) at src/G4QElasticCrossSection.cc:1221 1221 if(a<6.5) Current language: auto; currently c++ (gdb) where #0 G4QElasticCrossSection::GetTabValues (this=0xaf1f34a0, lp=1.826771653543307, PDG=2212, tgZ=82, tgN=125) at src/G4QElasticCrossSection.cc:1221 #1 0xaf1266b4 in G4QElasticCrossSection::GetPTables (this=0xaf1f34a0, LP=2.57306636095633, ILP=1.4488188976377951, PDG=2212, tgZ=82, tgN=125) at src/G4QElasticCrossSection.cc:801 #2 0xaf126f5a in G4QElasticCrossSection::CalculateCrossSection (this=0xaf1f34a0, CS=false, F=-1, I=11, PDG=2212, tgZ=82, tgN=125, pIU=13105.950463856872) at src/G4QElasticCrossSection.cc:378 #3 0xaf127e11 in G4QElasticCrossSection::GetCrossSection (this=0xaf1f34a0, fCS=false, pMom=13105.950463856872, tgZ=82, tgN=125, pPDG=2212) at src/G4QElasticCrossSection.cc:227 #4 0xaf45e6db in G4HadronElastic::ApplyYourself (this=0xdf73ce0, aTrack=@0xbfacc560, targetNucleus=@0xdf93020) at src/G4HadronElastic.cc:205 #5 0xaf465051 in G4UHadronElasticProcess::PostStepDoIt (this=0xdf92bc8, track=@0xebb1810, step=@0xbd7cba0) at src/G4UHadronElasticProcess.cc:310 #6 0xaf036b38 in G4SteppingManager::InvokePSDIP (this=0xbd7ca88, np=4) at src/G4SteppingManager2.cc:503 #7 0xaf036e1f in G4SteppingManager::InvokePostStepDoItProcs (this=0xbd7ca88) at src/G4SteppingManager2.cc:479 #8 0xaf032d4c in G4SteppingManager::Stepping (this=0xbd7ca88) at src/G4SteppingManager.cc:210 #9 0xaf0403a0 in G4TrackingManager::ProcessOneTrack (this=0xbd7ca60, apValueG4Track=0xebb1810) at src/G4TrackingManager.cc:126 #10 0xaef14b91 in G4EventManager::DoProcessing (this=0xbd7ca18, anEvent=0xbd88a80) at src/G4EventManager.cc:185 #11 0xaef15328 in G4EventManager::ProcessOneEvent (this=0xbd7ca18, anEvent=0xbd88a80) at src/G4EventManager.cc:335 #12 0xaeeb9db1 in G4RunManager::DoEventLoop (this=0xbd7c968, n_event=10, macroFile=0x0, n_select=-1) at src/G4RunManager.cc:235 #13 0xaeeb8724 in G4RunManager::BeamOn (this=0xbd7c968, n_event=10, macroFile=0x0, n_select=-1) at src/G4RunManager.cc:140 #14 0xae66fc29 in TG4RunManager::ProcessRun () from /home/donghee/GSI/fairsoft/transport/geant4_vmc/lib/tgt_linux/libgeant4vmc.so #15 0xae6743f9 in TGeant4::ProcessRun () from /home/donghee/GSI/fairsoft/transport/geant4_vmc/lib/tgt_linux/libgeant4vmc.so #16 0xb2b12ff5 in FairMCApplication::RunMC (this=0x87284c8, nofEvents=10) at /home/donghee/GSI/pandaroot/base/FairMCApplication.cxx:228 #17 0xb2b27887 in FairRunSim::Run (this=0x8620c40, NStart=10, NStop=0) at /home/donghee/GSI/pandaroot/base/FairRunSim.cxx:289 #18 0xb2b64c9c in G__FairDict_532_0_5 (result7=0xbfad347c, funcname=0x861ef28 "\001", libp=0xbfacde10, hash=0) at /home/donghee/GSI/buildpanda/base/FairDict.cxx:9292 #19 0xb70c7767 in Cint::G__ExceptionWrapper () from /home/donghee/GSI/fairsoft/tools/root/lib/libCint.so.5.20 #20 0xb7182c98 in G__execute_call () from /home/donghee/GSI/fairsoft/tools/root/lib/libCint.so.5.20 #21 0xb7183d90 in G__call_cppfunc () from /home/donghee/GSI/fairsoft/tools/root/lib/libCint.so.5.20 #22 0xb7157123 in G__interpret_func () from /home/donghee/GSI/fairsoft/tools/root/lib/libCint.so.5.20 #23 0xb7145369 in G__getfunction () from /home/donghee/GSI/fairsoft/tools/root/lib/libCint.so.5.20 #24 0xb723eaf4 in G__getstructmem () from /home/donghee/GSI/fairsoft/tools/root/lib/libCint.so.5.20 #25 0xb7233f7a in G__getvariable () from /home/donghee/GSI/fairsoft/tools/root/lib/libCint.so.5.20 ---Type to continue, or q to quit--- #26 0xb711802f in G__getitem () from /home/donghee/GSI/fairsoft/tools/root/lib/libCint.so.5.20 #27 0xb711f938 in G__getexpr () from /home/donghee/GSI/fairsoft/tools/root/lib/libCint.so.5.20 #28 0xb71b4a8c in G__exec_statement () from /home/donghee/GSI/fairsoft/tools/root/lib/libCint.so.5.20 #29 0xb7159385 in G__interpret_func () from /home/donghee/GSI/fairsoft/tools/root/lib/libCint.so.5.20 #30 0xb71454bd in G__getfunction () from /home/donghee/GSI/fairsoft/tools/root/lib/libCint.so.5.20 #31 0xb711865d in G__getitem () from /home/donghee/GSI/fairsoft/tools/root/lib/libCint.so.5.20 #32 0xb711f938 in G__getexpr () from /home/donghee/GSI/fairsoft/tools/root/lib/libCint.so.5.20 #33 0xb712b512 in G__calc_internal () from /home/donghee/GSI/fairsoft/tools/root/lib/libCint.so.5.20 #34 0xb71c155f in G__process_cmd () from /home/donghee/GSI/fairsoft/tools/root/lib/libCint.so.5.20 #35 0xb79a0250 in TCint::ProcessLine () from /home/donghee/GSI/fairsoft/tools/root/lib/libCore.so.5.20 #36 0xb799ff0f in TCint::ProcessLineSynch () from /home/donghee/GSI/fairsoft/tools/root/lib/libCore.so.5.20 #37 0xb78cff80 in TApplication::ExecuteFile () from /home/donghee/GSI/fairsoft/tools/root/lib/libCore.so.5.20 #38 0xb78d0514 in TApplication::ProcessFile () from /home/donghee/GSI/fairsoft/tools/root/lib/libCore.so.5.20 #39 0xb78cd584 in TApplication::ProcessLine () from /home/donghee/GSI/fairsoft/tools/root/lib/libCore.so.5.20 #40 0xb6f15bb4 in TRint::Run () from /home/donghee/GSI/fairsoft/tools/root/lib/libRint.so.5.20 #41 0x08048eb5 in main () (gdb) (gdb)