Message: problem of using the para hydrogen scattering kernel in geant4 Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None problem of using the para hydrogen scattering kernel in geant4 

Forum: Hadronic Processes
Date: 23 Apr, 2014
From: Lu, Wei <Lu, Wei>

Hi, I took the comments out in the file "G4NeutronHPThermalScatteringNames.cc" and enabled the scattering kernels for para hydrogen and others. From the output it seems to pick the up the scattering kernel for the para hydrogen as shown below:

"Neutron HP Thermal Scattering Data: Following material-element pairs and/or elements are registered. Element TS_H_of_Para_Hydrogen, internal thermal scattering id 1 Element TS_Beryllium_Metal, internal thermal scattering id 2 Material G4_WATER - Element H, internal thermal scattering id 0"

However, the program crashed later when I ran it on a Linux machine. The error information is attached below. Interestingly, it does not crash when I ran it on my Mac though the results are weird (a big jump where the scattering kernel starts at 4 eV, please find in the attached picture). I think the difference between the versions of Geant4 on Mac and Linux is one is 4.10.00.p01 and the other is 4.10.00. The Mac one is installed under "lib" instead of "lib64" on Linux.

My questions is 1) the para hydrogen scattering kernel reliable for use? 2) if so, how to fix the crashing problem in Linux? 3) is the result from the Mac reasonable? It is certainly not reasonable so far, I will have to long run to verify that. Thanks a lot.

Wei

=========================================================== There was a crash. This is the entire stack trace of all threads: ===========================================================

#0  0x00000038a5aac8be in waitpid () from /lib64/libc.so.6
#1  0x00000038a5a3e909 in do_system () from /lib64/libc.so.6
#2  0x00007fdb587c3188 in TUnixSystem::StackTrace() () from /SNS/users/7wi/root/lib/libCore.so
#3  0x00007fdb587c2003 in TUnixSystem::DispatchSignals(ESignals) () from /SNS/users/7wi/root/lib/libCore.so
#4  <signal handler called>
#5  0x00007fdb5b45e497 in G4NeutronHPThermalScattering::create_sE_and_EPM_from_pE_and_vE_P_E_isoAng(double, double, std::vector<E_P_E_isoAng*, std::allocator<E_P_E_isoAng*> >*) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4processes.so
#6  0x00007fdb5b462cdd in G4NeutronHPThermalScattering::ApplyYourself(G4HadProjectile const&, G4Nucleus&) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4processes.so
#7  0x00007fdb5b547a9e in G4HadronElasticProcess::PostStepDoIt(G4Track const&, G4Step const&) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4processes.so
#8  0x00007fdb5c02ff9b in G4SteppingManager::InvokePSDIP(unsigned long) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4tracking.so
#9  0x00007fdb5c0303ad in G4SteppingManager::InvokePostStepDoItProcs() () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4tracking.so
#10 0x00007fdb5c02e225 in G4SteppingManager::Stepping() () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4tracking.so
#11 0x00007fdb5c036d0c in G4TrackingManager::ProcessOneTrack(G4Track*) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4tracking.so
#12 0x00007fdb5c263a5e in G4EventManager::DoProcessing(G4Event*) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4event.so
#13 0x00007fdb5c4d2b47 in G4RunManager::ProcessOneEvent(int) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4run.so
#14 0x00007fdb5c4d2285 in G4RunManager::DoEventLoop(int, char const*, int) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4run.so
#15 0x00007fdb5c4d2222 in G4RunManager::BeamOn(int, char const*, int) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4run.so
#16 0x00007fdb5c4e582a in G4RunMessenger::SetNewValue(G4UIcommand*, G4String) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4run.so
#17 0x00007fdb5964988b in G4UIcommand::DoIt(G4String) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4intercoms.so
#18 0x00007fdb59655fc7 in G4UImanager::ApplyCommand(char const*) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4intercoms.so
#19 0x00007fdb5963b0d7 in G4UIbatch::ExecCommand(G4String const&) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4intercoms.so
#20 0x00007fdb5963c09b in G4UIbatch::SessionStart() () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4intercoms.so
#21 0x00007fdb59653b33 in G4UImanager::ExecuteMacroFile(char const*) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4intercoms.so
#22 0x00007fdb5965002f in G4UIcontrolMessenger::SetNewValue(G4UIcommand*, G4String) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4intercoms.so
#23 0x00007fdb5964988b in G4UIcommand::DoIt(G4String) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4intercoms.so
#24 0x00007fdb59655fc7 in G4UImanager::ApplyCommand(char const*) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4intercoms.so
#25 0x0000000000414319 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  0x00007fdb5b45e497 in G4NeutronHPThermalScattering::create_sE_and_EPM_from_pE_and_vE_P_E_isoAng(double, double, std::vector<E_P_E_isoAng*, std::allocator<E_P_E_isoAng*> >*) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4processes.so
#6  0x00007fdb5b462cdd in G4NeutronHPThermalScattering::ApplyYourself(G4HadProjectile const&, G4Nucleus&) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4processes.so
#7  0x00007fdb5b547a9e in G4HadronElasticProcess::PostStepDoIt(G4Track const&, G4Step const&) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4processes.so
#8  0x00007fdb5c02ff9b in G4SteppingManager::InvokePSDIP(unsigned long) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4tracking.so
#9  0x00007fdb5c0303ad in G4SteppingManager::InvokePostStepDoItProcs() () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4tracking.so
#10 0x00007fdb5c02e225 in G4SteppingManager::Stepping() () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4tracking.so
#11 0x00007fdb5c036d0c in G4TrackingManager::ProcessOneTrack(G4Track*) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4tracking.so
#12 0x00007fdb5c263a5e in G4EventManager::DoProcessing(G4Event*) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4event.so
#13 0x00007fdb5c4d2b47 in G4RunManager::ProcessOneEvent(int) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4run.so
#14 0x00007fdb5c4d2285 in G4RunManager::DoEventLoop(int, char const*, int) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4run.so
#15 0x00007fdb5c4d2222 in G4RunManager::BeamOn(int, char const*, int) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4run.so
#16 0x00007fdb5c4e582a in G4RunMessenger::SetNewValue(G4UIcommand*, G4String) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4run.so
#17 0x00007fdb5964988b in G4UIcommand::DoIt(G4String) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4intercoms.so
#18 0x00007fdb59655fc7 in G4UImanager::ApplyCommand(char const*) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4intercoms.so
#19 0x00007fdb5963b0d7 in G4UIbatch::ExecCommand(G4String const&) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4intercoms.so
#20 0x00007fdb5963c09b in G4UIbatch::SessionStart() () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4intercoms.so
#21 0x00007fdb59653b33 in G4UImanager::ExecuteMacroFile(char const*) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4intercoms.so
#22 0x00007fdb5965002f in G4UIcontrolMessenger::SetNewValue(G4UIcommand*, G4String) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4intercoms.so
#23 0x00007fdb5964988b in G4UIcommand::DoIt(G4String) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4intercoms.so
#24 0x00007fdb59655fc7 in G4UImanager::ApplyCommand(char const*) () from /SNS/users/7wi/geant4.10.00.p01/lib64/libG4intercoms.so
#25 0x0000000000414319 in main ()
===========================================================

Segmentation fault (core dumped)

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2014/04/23/11.13-30674-spectrum.tiff.pdf

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: problem of using the para hydrogen scattering kernel in geant4   (Koi, Tatsumi - 23 Apr, 2014)
 Add Message Add Message
to: "problem of using the para hydrogen scattering kernel in geant4"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]