G4PhysListFactory Class Reference

#include <G4PhysListFactory.hh>

Public Member Functions
	G4PhysListFactory (G4int ver=1)
	~G4PhysListFactory ()
G4VModularPhysicsList *	GetReferencePhysList (const G4String &)
G4VModularPhysicsList *	ReferencePhysList ()
G4bool	IsReferencePhysList (const G4String &) const
const std::vector< G4String > &	AvailablePhysLists () const
const std::vector< G4String > &	AvailablePhysListsEM () const
void	SetVerbose (G4int val)

Detailed Description

Definition at line 45 of file G4PhysListFactory.hh.

Constructor & Destructor Documentation

G4PhysListFactory()

G4PhysListFactory::G4PhysListFactory

(

G4int

ver = 1

)

Definition at line 82 of file G4PhysListFactory.cc.

  : defName("FTFP_BERT"),verbose(ver),theMessenger(nullptr)
{
  nlists_hadr = 36;
  G4String ss[36] = {
    "FTFP_BERT","FTFP_BERT_TRV","FTFP_BERT_ATL","FTFP_BERT_HP","FTFQGSP_BERT",
    "FTFP_INCLXX","FTFP_INCLXX_HP","FTF_BIC","LBE","QBBC",
    "QGSP_BERT","QGSP_BERT_HP","QGSP_BIC","QGSP_BIC_HP","QGSP_BIC_AllHP",
    "QGSP_FTFP_BERT","QGSP_INCLXX","QGSP_INCLXX_HP","QGS_BIC",
    "Shielding","ShieldingLEND","ShieldingLIQMD","ShieldingM","NuBeam",
    "Shielding_HP","ShieldingLIQMD_HP","ShieldingM_HP",
    "FTFP_BERT_HPT","FTFP_INCLXX_HPT","QGSP_BERT_HPT","QGSP_BIC_HPT",
    "QGSP_BIC_AllHPT","QGSP_INCLXX_HPT","Shielding_HPT","ShieldingLIQMD_HPT",
    "ShieldingM_HPT"};
  for(size_t i=0; i<nlists_hadr; ++i) {
    listnames_hadr.push_back(ss[i]);
  }
 
  nlists_em = 12;
  G4String s1[12] = {"","_EMV","_EMX","_EMY","_EMZ","_LIV","_PEN",
             "__GS","__SS","_EM0","_WVI","__LE"};
  for(size_t i=0; i<nlists_em; ++i) {
    listnames_em.push_back(s1[i]);
  }
}

~G4PhysListFactory()

G4PhysListFactory::~G4PhysListFactory

(

)

Definition at line 108 of file G4PhysListFactory.cc.

{
  delete theMessenger;
}

Member Function Documentation

AvailablePhysLists()

const std::vector< G4String > & G4PhysListFactory::AvailablePhysLists

(

)

const

Definition at line 271 of file G4PhysListFactory.cc.

{
  return listnames_hadr;
}

AvailablePhysListsEM()

const std::vector< G4String > & G4PhysListFactory::AvailablePhysListsEM

(

)

const

Definition at line 277 of file G4PhysListFactory.cc.

{
  return listnames_em;
}

GetReferencePhysList()

G4VModularPhysicsList * G4PhysListFactory::GetReferencePhysList

(

const G4String &

name

)

Definition at line 134 of file G4PhysListFactory.cc.

{
  // analysis on the string 
  size_t n = name.size();
 
  // last characters in the string
  size_t em_opt = 0;
  G4String em_name = "";
 
  // check EM options
  if(n > 4) {
    em_name = name.substr(n - 4, 4);
    for(size_t i=1; i<nlists_em; ++i) { 
      if(listnames_em[i] == em_name) { 
    em_opt = i;
        n -= 4;
        break; 
      }
    }
    if(0 == em_opt) { em_name = ""; }
  }
 
  // hadronic pHysics List
  G4String had_name = name.substr(0, n);
 
  if(0 < verbose) {
    G4cout << "G4PhysListFactory::GetReferencePhysList <" << had_name
       << em_name << ">  EMoption= " << em_opt << G4endl;
  }
  G4VModularPhysicsList* p = nullptr;
  if(had_name == "FTFP_BERT")               {p = new FTFP_BERT(verbose);}
  else if(had_name == "FTFP_BERT_HP")       {p = new FTFP_BERT_HP(verbose);}
  else if(had_name == "FTFP_BERT_TRV")      {p = new FTFP_BERT_TRV(verbose);}
  else if(had_name == "FTFP_BERT_ATL")      {p = new FTFP_BERT_ATL(verbose);}
  else if(had_name == "FTFQGSP_BERT")       {p = new FTFQGSP_BERT(verbose);}
  else if(had_name == "FTFP_INCLXX")        {p = new FTFP_INCLXX(verbose);}
  else if(had_name == "FTFP_INCLXX_HP")     {p = new FTFP_INCLXX_HP(verbose);}
  else if(had_name == "FTF_BIC")            {p = new FTF_BIC(verbose);}
  else if(had_name == "LBE")                {p = new LBE();}
  else if(had_name == "QBBC")               {p = new QBBC(verbose);}
  else if(had_name == "QGSP_BERT")          {p = new QGSP_BERT(verbose);}
  else if(had_name == "QGSP_BERT_HP")       {p = new QGSP_BERT_HP(verbose);}
  else if(had_name == "QGSP_BIC")           {p = new QGSP_BIC(verbose);}
  else if(had_name == "QGSP_BIC_HP")        {p = new QGSP_BIC_HP(verbose);}
  else if(had_name == "QGSP_BIC_AllHP")     {p = new QGSP_BIC_AllHP(verbose);}
  else if(had_name == "QGSP_FTFP_BERT")     {p = new QGSP_FTFP_BERT(verbose);}
  else if(had_name == "QGSP_INCLXX")        {p = new QGSP_INCLXX(verbose);}
  else if(had_name == "QGSP_INCLXX_HP")     {p = new QGSP_INCLXX_HP(verbose);}
  else if(had_name == "QGS_BIC")            {p = new QGS_BIC(verbose);}
  else if(had_name == "Shielding")          {p = new Shielding(verbose);}
  else if(had_name == "ShieldingLEND")      {p = new ShieldingLEND(verbose);}
  else if(had_name == "ShieldingLIQMD")     {p = new Shielding(verbose,"HP","",true);}  
  else if(had_name == "ShieldingM")         {p = new Shielding(verbose,"HP","M");}
  else if(had_name == "NuBeam")             {p = new NuBeam(verbose);}
  else if(had_name == "Shielding_HP")       {p = new Shielding(verbose);}
  else if(had_name == "ShieldingLIQMD_HP")  {p = new Shielding(verbose,"HP","",true);}
  else if(had_name == "ShieldingM_HP")      {p = new Shielding(verbose,"HP","M");}
  else if(had_name == "FTFP_BERT_HPT")      {p = new FTFP_BERT_HP(verbose);
                                             p->RegisterPhysics(new G4ThermalNeutrons);}
  else if(had_name == "FTFP_INCLXX_HPT")    {p = new FTFP_INCLXX_HP(verbose);
                                             p->RegisterPhysics(new G4ThermalNeutrons);}
  else if(had_name == "QGSP_BERT_HPT")      {p = new QGSP_BERT_HP(verbose);
                                             p->RegisterPhysics(new G4ThermalNeutrons);}
  else if(had_name == "QGSP_BIC_HPT")       {p = new QGSP_BIC_HPT(verbose);}
  else if(had_name == "QGSP_BIC_AllHPT")    {p = new QGSP_BIC_AllHP(verbose);
                                             p->RegisterPhysics(new G4ThermalNeutrons);}
  else if(had_name == "QGSP_INCLXX_HPT")    {p = new QGSP_INCLXX_HP(verbose);
                                             p->RegisterPhysics(new G4ThermalNeutrons);}
  else if(had_name == "Shielding_HPT")      {p = new Shielding(verbose);
                                             p->RegisterPhysics(new G4ThermalNeutrons);}
  else if(had_name == "ShieldingLIQMD_HPT") {p = new Shielding(verbose,"HP","",true);
                                             p->RegisterPhysics(new G4ThermalNeutrons);}
  else if(had_name == "ShieldingM_HPT")     {p = new Shielding(verbose,"HP","M");
                                             p->RegisterPhysics(new G4ThermalNeutrons);}
  else {
    p = new FTFP_BERT(verbose);
    G4ExceptionDescription ed;
    ed << "PhysicsList " << had_name << " is not known;"
       << " the default FTFP_BERT is created";
    G4Exception("G4PhysListFactory: ","pl0003",JustWarning,ed,"");
  }
  if(nullptr != p) {
    if(0 < em_opt && had_name != "LBE") {
      if(1 == em_opt) { 
    p->ReplacePhysics(new G4EmStandardPhysics_option1(verbose)); 
      } else if(2 == em_opt) {
    p->ReplacePhysics(new G4EmStandardPhysics_option2(verbose)); 
      } else if(3 == em_opt) {
    p->ReplacePhysics(new G4EmStandardPhysics_option3(verbose)); 
      } else if(4 == em_opt) {
    p->ReplacePhysics(new G4EmStandardPhysics_option4(verbose)); 
      } else if(5 == em_opt) {
    p->ReplacePhysics(new G4EmLivermorePhysics(verbose)); 
      } else if(6 == em_opt) {
    p->ReplacePhysics(new G4EmPenelopePhysics(verbose)); 
      } else if(7 == em_opt) {
    p->ReplacePhysics(new G4EmStandardPhysicsGS(verbose)); 
      } else if(8 == em_opt) {
    p->ReplacePhysics(new G4EmStandardPhysicsSS(verbose)); 
      } else if(9 == em_opt) {
    p->ReplacePhysics(new G4EmStandardPhysics(verbose)); 
      } else if(10 == em_opt) {
    p->ReplacePhysics(new G4EmStandardPhysicsWVI(verbose)); 
      } else if(11 == em_opt) {
    p->ReplacePhysics(new G4EmLowEPPhysics(verbose)); 
      }
    }
    theMessenger = new G4PhysListFactoryMessenger(p);
  }
  if(0 < verbose) G4cout << G4endl;
  return p;
}

Referenced by ReferencePhysList().

IsReferencePhysList()

G4bool G4PhysListFactory::IsReferencePhysList

(

const G4String &

name

)

const

Definition at line 247 of file G4PhysListFactory.cc.

{
  G4bool res = false;
  size_t n = name.size();
  if(n > 4) {
    G4String em_name = name.substr(n - 4, 4);
    for(size_t i=1; i<nlists_em; ++i) { 
      if(listnames_em[i] == em_name) { 
        n -= 4;
        break; 
      }
    }
  }
  G4String had_name = name.substr(0, n);
  for(size_t i=0; i<nlists_hadr; ++i) {
    if(had_name == listnames_hadr[i]) {
      res = true;
      break;
    }
  }
  return res;
}

ReferencePhysList()

G4VModularPhysicsList * G4PhysListFactory::ReferencePhysList

(

)

Definition at line 114 of file G4PhysListFactory.cc.

{
  // instantiate PhysList by environment variable "PHYSLIST"
  G4String name = "";
  char* path = std::getenv("PHYSLIST");
  if (path) {
    name = G4String(path);
  } else {
    name = defName;
    G4cout << "### G4PhysListFactory WARNING: "
       << " environment variable PHYSLIST is not defined"
       << G4endl
       << "    Default Physics Lists " << name 
       << " is instantiated" 
       << G4endl;
  }
  return GetReferencePhysList(name);
}

SetVerbose()

void G4PhysListFactory::SetVerbose ( G4int val )

inline

Definition at line 68 of file G4PhysListFactory.hh.

{ verbose = val; }

---

The documentation for this class was generated from the following files:

• G4PhysListFactory.hh
• G4PhysListFactory.cc