mdurk
commented
2 years ago Open mdurk opened 2 years ago
mdurk
commented
2 years ago 
lawrenceleejr
commented
2 years ago Looks like we need to set
export G4RADIOACTIVEDATA=/usr/local/share/Geant4-10.7.3/data/RadioactiveDecay5.6/
root@DESKTOP-9038P3Q:/home/mdurkee/work/ScatteringRadiation-master/build# ./run_beam
Creating output file output_10micron_1e4events.root
**************************************************************
Geant4 version Name: geant4-10-07-patch-03 (19-November-2021)
Copyright : Geant4 Collaboration
References : NIM A 506 (2003), 250-303
: IEEE-TNS 53 (2006), 270-278
: NIM A 835 (2016), 186-225
WWW : http://geant4.org/
**************************************************************
<<< Geant4 Physics List simulation engine: QGSP_BERT_HP
Checking overlaps for volume World (G4Box) ... OK!
Checking overlaps for volume physDetector (G4Box) ... OK!
hInelastic QGSP_BERT_HP Thresholds:
1) between BERT and FTF/P over the interval 3 to 6 GeV.
2) between FTF/P and QGS/P over the interval 12 to 25 GeV.
-- quasiElastic: 1 for QGS and 0 for FTF
terminate called after throwing an instance of 'G4HadronicException'
what(): In /usr/share/geant4.10.07.p03/source/processes/hadronic/models/particle_hp/src/G4ParticleHPInelastic.cc, line 88:
===> Please setenv G4PARTICLEHPDATA (recommended) or, at least setenv G4NEUTRONHPDATA to point to the neutron cross-section files.export G4NEUTRONHPDATA=/usr/local/share/Geant4-10.7.3/data/G4NDL4.6/ Hadronic Processes for triton
Process: hadElastic
Model: hElasticLHEP: 0 eV /n ---> 100 TeV/n
Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV
Process: tInelastic
Model: Binary Light Ion Cascade: 0 eV /n ---> 6 GeV/n
Model: FTFP: 3 GeV/n ---> 100 TeV/n
Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV
================================================================
=======================================================================
====== Pre-compound/De-excitation Physics Parameters ========
=======================================================================
Type of pre-compound inverse x-section 3
Pre-compound model active 1
Pre-compound excitation low energy (MeV) 0.1
Pre-compound excitation high energy (MeV) 30
Type of de-excitation inverse x-section 3
Type of de-excitation factory Evaporation+GEM
Number of de-excitation channels 68
Min excitation energy (keV) 0.01
Min energy per nucleon for multifragmentation (MeV) 2e+05
Limit excitation energy for Fermi BreakUp (MeV) 20
Level density (1/MeV) 0.075
Use simple level density model 1
Use discrete excitation energy of the residual 0
Time limit for long lived isomeres (ns) 1442.7
Isomer production flag 1
Internal e- conversion flag 1
Store e- internal conversion data 1
Electron internal conversion ID 0
Correlated gamma emission flag 0
Max 2J for sampling of angular correlations 10
=======================================================================Checking overlaps for volume World (G4Box) ... OK!
Checking overlaps for volume physDetector (G4Box) ... OK!
hInelastic QGSP_BERT_HP Thresholds:
1) between BERT and FTF/P over the interval 3 to 6 GeV.
2) between FTF/P and QGS/P over the interval 12 to 25 GeV.
-- quasiElastic: 1 for QGS and 0 for FTF
/usr/local/share/Geant4-10.7.3/data/G4NDL4.6/
@@@ G4ParticleHPInelastic instantiated for particle neutron data directory variable is G4NEUTRONHPDATA pointing to /usr/local/share/Geant4-10.7.3/data/G4NDL4.6//Inelastic
@@@ G4ParticleHPInelasticData instantiated for particle neutron data directory variable is G4NEUTRONHPDATA pointing to /usr/local/share/Geant4-10.7.3/data/G4NDL4.6/
NeutronHP: /Capture file for Z = 6, A = 12 is not found and NeutronHP will use /usr/local/share/Geant4-10.7.3/data/G4NDL4.6//Capture/CrossSection/6_nat_Carbon
NeutronHP: /Elastic file for Z = 6, A = 12 is not found and NeutronHP will use /usr/local/share/Geant4-10.7.3/data/G4NDL4.6//Elastic/CrossSection/6_nat_Carbon
NeutronHP: /Inelastic file for Z = 6, A = 12 is not found and NeutronHP will use /usr/local/share/Geant4-10.7.3/data/G4NDL4.6//Inelastic/CrossSection/6_nat_Carbon
NeutronHP: /Capture file for Z = 73, A = 180 is not found and NeutronHP will use /usr/local/share/Geant4-10.7.3/data/G4NDL4.6//Capture/CrossSection/73_181_Tantalum
NeutronHP: /Elastic file for Z = 73, A = 180 is not found and NeutronHP will use /usr/local/share/Geant4-10.7.3/data/G4NDL4.6//Elastic/CrossSection/73_181_Tantalum
NeutronHP: /Inelastic file for Z = 73, A = 180 is not found and NeutronHP will use /usr/local/share/Geant4-10.7.3/data/G4NDL4.6//Inelastic/CrossSection/73_181_Tantalum
Will save outputs to output_10micron_1e4events.root
Available UI session types: [ Qt, GAG, tcsh, csh ]
QStandardPaths: XDG_RUNTIME_DIR not set, defaulting to '/tmp/runtime-root'
parameter value (point) is not listed in the candidate List.@kpachal -- @mdurk has gotten your latest commits to run after some env fiddling. Using QBBC he seems to get a root tree that describes no interactions. But after this switch to QGSP_BERT, he gets no output file. He's going to dig into it but if you have any quick ideas for why it might not produce an output file (despite it not complaining at runtime) that'd be awesome.
mdurk
commented
2 years ago I noticed that it is currently complaining at runtime
================================================================
=======================================================================
====== Pre-compound/De-excitation Physics Parameters ========
=======================================================================
Type of pre-compound inverse x-section 3
Pre-compound model active 1
Pre-compound excitation low energy (MeV) 0.1
Pre-compound excitation high energy (MeV) 30
Type of de-excitation inverse x-section 3
Type of de-excitation factory Evaporation+GEM
Number of de-excitation channels 68
Min excitation energy (keV) 0.01
Min energy per nucleon for multifragmentation (MeV) 2e+05
Limit excitation energy for Fermi BreakUp (MeV) 20
Level density (1/MeV) 0.075
Use simple level density model 1
Use discrete excitation energy of the residual 0
Time limit for long lived isomeres (ns) 1442.7
Isomer production flag 1
Internal e- conversion flag 1
Store e- internal conversion data 1
Electron internal conversion ID 0
Correlated gamma emission flag 0
Max 2J for sampling of angular correlations 10
=======================================================================
munmap_chunk(): invalid pointer
Aborted (core dumped)Any ideas as to what this pointer is complaining about are welcome.
mdurk
commented
2 years ago I went back to the QBBC physics and am having the same issue. I have traced the problem to the beam setup, specifically commenting out line 170 allows the program to complete [ uiManager->ApplyCommand("/gps/direction 0 0 1"); ]
// Set up the beam. Can either do point-like,
// for studying scattering probability of a single electron,
// or realistic, for studying impact of non-constant targets.
uiManager->ApplyCommand("/gps/particle e-");
uiManager->ApplyCommand("/gps/ene/type Mono");
uiManager->ApplyCommand("/gps/ene/mono 31 MeV");
// Same center and direction regardless of what beam type it is
uiManager->ApplyCommand("/gps/pos/centre 0. 0. -1. m");
//uiManager->ApplyCommand("/gps/direction 0 0 1");There is an issue with the memory allocation of the gps/direction I suppose
kpachal
commented
2 years ago Hm, that's very weird. Referring to the problem Larry saw earlier (the tree with no interactions), try again after the git pull and see if it fixes that issue. As for the segfault, does it only appear with the one physics list? If leaving the gps/direction line in place works for one physics list but not the other, it feels like the problem is there somehow ...
Note: git pull will change to a 3rd physics list, on advice from some collaborators about what does best for low energy nuclear processes.