Message: Re: neutron production via (p,n) reaction: problem with *_BERT_HP physics list Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Feedback Re: neutron production via (p,n) reaction: problem with *_BERT_HP physics list 

Forum: Physics List
Re: Question neutron production via (p,n) reaction: problem with *_BERT_HP physics list (J)
Re: Idea Re: neutron production via (p,n) reaction: problem with *_BERT_HP physics list (Vladimir Ivanchenko)
Date: 28 Oct, 2008
From: J <J>

Hi, the problem is still not solved... unfortunately. I use the G4VModularPhysicsList. When I comment out the setcuts method by:

#include "ErsterVersuchPhysicsList.hh"
#include "G4ParticleTypes.hh"

#include "HadronPhysicsQGSP_BERT_HP.hh"

ErsterVersuchPhysicsList::ErsterVersuchPhysicsList() : G4VModularPhysicsList() {

  RegisterPhysics( new HadronPhysicsQGSP_BERT_HP("hadron") );
}

ErsterVersuchPhysicsList::~ErsterVersuchPhysicsList() {;}

//void ErsterVersuchPhysicsList::SetCuts()	// setting a range cut value for all particles
//{
  // uppress error messages            
//  G4int temp = GetVerboseLevel();
//  SetVerboseLevel(0);                                                           
  //  " G4VUserPhysicsList::SetCutsWithDefault" method sets 
  //   the default cut value for all particle types 
//  SetCutsWithDefault();   

  // Retrieve verbose level
//  SetVerboseLevel(temp);  
//}

I cannot compile. With the SetCutsWithDefault() I get the same output (compiling works naturally)...

->

========= Table of registered couples ==============================

Index : 0     used in the geometry : Yes     recalculation needed : No 
 Material : Vacuum
 Range cuts        :  gamma 0 fm     e- 0 fm     e+ 0 fm 
 Energy thresholds :  gamma -1 MeV    e- -1 MeV    e+ -1 MeV
 Region(s) which use this couple : 
    DefaultRegionForTheWorld

Index : 1     used in the geometry : Yes     recalculation needed : No 
 Material : Lithium
 Range cuts        :  gamma 0 fm     e- 0 fm     e+ 0 fm 
 Energy thresholds :  gamma -1 MeV    e- -1 MeV    e+ -1 MeV
 Region(s) which use this couple : 
    DefaultRegionForTheWorld

====================================================================

I did not manage to use the example from examples/extended/hadronic/Hadr01 for the setcuts yet. I think this could be because of the G4VModularPhysicsList(), but I'm not sure. When I use the part of the code about the setcuts I get the right range cuts in my output, but the hadronic physics seem not to work. Should I implement the hadronic physics in another way?

A colleague also experienced some trouble with these physics lists and the cut ranges...

Kind regards, Jörg

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

1 News: Re: neutron production via (p,n) reaction: problem with *_BERT_HP physics list   (Jörg Schütrumpf - 28 Oct, 2008)
2 More: Re: neutron production via (p,n) reaction: problem with *_BERT_HP physics list   (Jörg Schütrumpf - 28 Oct, 2008)
1 Idea: Re: neutron production via (p,n) reaction: problem with *_BERT_HP physics list   (Vladimir Ivanchenko - 05 Nov, 2008)
 Add Message Add Message
to: "Re: neutron production via (p,n) reaction: problem with *_BERT_HP physics list"

 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 ]