5 #include "G4DecayPhysics.hh"
6 #include "G4EmStandardPhysics_option4.hh"
7 #include "G4EmExtraPhysics.hh"
8 #include "G4IonPhysics.hh"
9 #include "G4StoppingPhysics.hh"
10 #include "G4HadronElasticPhysics.hh"
11 #include "G4NeutronTrackingCut.hh"
12 #include "G4HadronicProcessStore.hh"
15 int ver =
p.getUntrackedParameter<
int>(
"Verbosity", 0);
16 bool emPhys =
p.getUntrackedParameter<
bool>(
"EMPhysics",
true);
17 bool hadPhys =
p.getUntrackedParameter<
bool>(
"HadPhysics",
true);
18 bool tracking =
p.getParameter<
bool>(
"TrackingCut");
19 double timeLimit =
p.getParameter<
double>(
"MaxTrackTime") * CLHEP::ns;
20 double minFTFP =
p.getParameter<
double>(
"EminFTFP") *
CLHEP::GeV;
21 double maxBERT =
p.getParameter<
double>(
"EmaxBERT") *
CLHEP::GeV;
22 double maxBERTpi =
p.getParameter<
double>(
"EmaxBERTpi") *
CLHEP::GeV;
23 edm::LogVerbatim(
"PhysicsList") <<
"You are using the simulation engine: FTFP_BERT_EMZ \n Flags for EM Physics "
24 << emPhys <<
", for Hadronic Physics " << hadPhys <<
" and tracking cut " <<
tracking
25 <<
" t(ns)= " << timeLimit / CLHEP::ns <<
"\n transition energy Bertini/FTFP from "
31 RegisterPhysics(
new G4EmStandardPhysics_option4(ver));
34 G4EmExtraPhysics* gn =
new G4EmExtraPhysics(ver);
39 this->RegisterPhysics(
new G4DecayPhysics(ver));
42 G4HadronicProcessStore::Instance()->SetVerbose(ver);
45 RegisterPhysics(
new G4HadronElasticPhysics(ver));
51 RegisterPhysics(
new G4StoppingPhysics(ver));
54 RegisterPhysics(
new G4IonPhysics(ver));
58 G4NeutronTrackingCut* ncut =
new G4NeutronTrackingCut(ver);
59 ncut->SetTimeLimit(timeLimit);
60 RegisterPhysics(ncut);