|Message: Problem when combining QGSP_BERT_HP and Penelope||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
I am currently writing an application to simulate low energy neutron interaction in 4 volumes which consist of four Boxes made of Polyethylene, CR39, PMMA and Water respectively. I use Geant4 9.6p03 and its associated data files. I want to use the QGSP_BERT_HP physics list with activation of thermal neutron physics and the Penelope Physics List for the electromagnetic processes at Low Energy. When running, my application crashes : -------- EEEE ------- G4Exception-START -------- EEEE ------- *** G4Exception : em2037 issued by : G4PenelopeOscillatorManager::BuildOscillatorTable() Error in Sterheimer factor calculation 5.82907e+11 56.3994 *** Fatal Exception *** core dump *** -------- EEEE -------- G4Exception-END --------- EEEE ------- *** G4Exception: Aborting execution *** Since I want to use the neutron thermal scattering in my application, I added TS_H_of_Polyethylene and TS_H_of_Water in my list of materials following: a=1.01*g/mole; elH=new G4Element(name="Hydrogen",symbol="H2",z=1.,a); //Definition of Hydrogen for activation of thermal neutron physics elHPolyEthy = new G4Element(name="TS_H_of_Polyethylene","H",z=1., 1.01) ; elHWater = new G4Element(name="TS_H_of_Water","H",z=1., 1.01) ; When I define Polyethylene with elHPolyEthy following: density = 0.94*g/cm3; Polyethylene = new G4Material(name="Polyethylene",density,nel=2, kStateSolid, 293.15*kelvin); Polyethylene->AddElement(elC, fractionmass=0.856289); Polyethylene->AddElement(elHPolyEthy, fractionmass=0.143711); and use the so-defined material in my Detector Construction, the program crashes and gives the error message mentionned above When I replace elHPolyEthy by elH, the program runs. The same crash occurs when defining Water following: density=1.0*g/cm3 ; Water= new G4Material(name="Water", density, nel=2, kStateLiquid, 293.15*kelvin) ; Water->AddElement(elHWater, fractionmass=0.111894) ; Water->AddElement(elO, fractionmass=0.888106) ; If I replace elHWater by elH, the program runs. The application crashes the same way when I use the NIST manager tool to define Polyethylene and Water, i.e respectively by: G4NistManager* man = G4NistManager::Instance() ; man->SetVerbose(1) ; .... Water = man->FindOrBuildMaterial("G4_WATER") ; Polyethylene = man->FindOrBuildMaterial("G4_POLYETHYLENE") ; ... As a last test, I run the program without Penelope in the Physics List but with QGSP_BERT_HP with thermal neutron and by defining Polyethylene and Water with elHPolyEthy and elHWater respectively. The program runs successfully. So right now, I am deep in trouble since I need to use both the Neutron Physics List( QGSP_BERT_HP ) including thermal scattering and the Penelope Physics to simulate Low Energy Electromagnetic processes with gamma, e-, e+, proton etc . But it seems there is an incompatibility hindering the combination of these both lists. As a remark, I use a Geantino as a primary particle for these tests, i.e the problem should occur when loading materials and physics. Has such a problem already been reported ? How can I combine both Low energy electromagnetic physics with QGSP_BERT_HP including thermal neutron scattering ? Many thanks for your help and suggestion ...
|Inline Depth:||Outline Depth:||Add message:|