Message: Re: Need help setting up physics list Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Idea Re: Need help setting up physics list 

Forum: Physics List
Re: None Need help setting up physics list (Tobias Kunz)
Date: 03 Feb, 2012
From: Andrea Dotti <Andrea Dotti>

Hello, if your class inherits from G4VModularPhysicsList you may find builders quite useful. For example for your energy range I suggest you to use the FTF (high energy model) and BERT (low energy) models to cover hadroncis physics (if you are using G4 version 9.4 or more recent).

Thanks to builders you can add all hadronics with the following line in the constructor of your physics list (agin this works only if your class inherits from G4VModularPhysicsList):

   // Hadron Elastic scattering
  this-> RegisterPhysics( new G4HadronElasticPhysics("elastic",ver,false));

   // Hadron Physics
   G4bool quasiElastic;
  this->RegisterPhysics(  new HadronPhysicsFTFP_BERT("hadron",quasiElastic=false));

  // Stopping Physics
  this->RegisterPhysics( new G4QStoppingPhysics("stopping"));

  // Ion Physics
  this->RegisterPhysics( new G4IonPhysics("ion"));

  // Neutron tracking cut
  this->RegisterPhysics( new G4NeutronTrackingCut("Neutron tracking cut", ver));

Please take a look at the source files FTFP_BERT.hh and FTFP_BERT.icc to see an example of a pre-packaged physics list provided with G4.

If you need extra accuracy for low energy (<20MeV) neutrons (at the price of higher CPU consumption), take a look at QGSP_BERT_HP or Shielding (available only in recent G4 versions) physics lists.

 Add Message Add Message
to: "Re: Need help setting up 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 ]