6 #include "G4DecayPhysics.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_EMM: \n Flags for EM Physics: "
24 << emPhys <<
"; Hadronic Physics: " << hadPhys <<
"; tracking cut: " <<
tracking
25 <<
"; time limit(ns)= " << timeLimit / CLHEP::ns
26 <<
"\n transition energy Bertini/FTFP from " << minFTFP /
CLHEP::GeV <<
" to "
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);