Message: Problem when combining QGSP_BERT_HP and Penelope Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Problem when combining QGSP_BERT_HP and Penelope 

Forum: Physics List
Date: 09 Jul, 2014
From: Thomas Zerguerras <Thomas Zerguerras>

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:


 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:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: Problem when combining QGSP_BERT_HP and Penelope   (Luciano Pandola - 17 Jul, 2014)
2 None: Re: Problem when combining QGSP_BERT_HP and Penelope   (Luciano Pandola - 17 Jul, 2014)
1 Ok: Re: Problem when combining QGSP_BERT_HP and Penelope   (Thomas Zerguerras - 17 Jul, 2014)
 Add Message Add Message
to: "Problem when combining QGSP_BERT_HP and Penelope"

 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 ]