CMS 3D CMS Logo

FTFPCMS_BERT_EMM.cc
Go to the documentation of this file.
1 #include "FTFPCMS_BERT_EMM.h"
5 
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"
13 
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 "
27  << maxBERT / CLHEP::GeV << ":" << maxBERTpi / CLHEP::GeV << " GeV";
28 
29  if (emPhys) {
30  // EM Physics
31  RegisterPhysics(new CMSEmStandardPhysicsLPM(ver));
32 
33  // Synchroton Radiation & GN Physics
34  G4EmExtraPhysics* gn = new G4EmExtraPhysics(ver);
35  RegisterPhysics(gn);
36  }
37 
38  // Decays
39  this->RegisterPhysics(new G4DecayPhysics(ver));
40 
41  if (hadPhys) {
42  G4HadronicProcessStore::Instance()->SetVerbose(ver);
43 
44  // Hadron Elastic scattering
45  RegisterPhysics(new G4HadronElasticPhysics(ver));
46 
47  // Hadron Physics
48  RegisterPhysics(new CMSHadronPhysicsFTFP_BERT(minFTFP, maxBERT, maxBERTpi));
49 
50  // Stopping Physics
51  RegisterPhysics(new G4StoppingPhysics(ver));
52 
53  // Ion Physics
54  RegisterPhysics(new G4IonPhysics(ver));
55 
56  // Neutron tracking cut
57  if (tracking) {
58  G4NeutronTrackingCut* ncut = new G4NeutronTrackingCut(ver);
59  ncut->SetTimeLimit(timeLimit);
60  RegisterPhysics(ncut);
61  }
62  }
63 }
MessageLogger.h
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
CMSHadronPhysicsFTFP_BERT.h
tracking
Definition: TempMeasurements.h:8
PhysicsList
Definition: PhysicsList.h:7
edm::ParameterSet
Definition: ParameterSet.h:47
GeV
const double GeV
Definition: MathUtil.h:16
FTFPCMS_BERT_EMM::FTFPCMS_BERT_EMM
FTFPCMS_BERT_EMM(const edm::ParameterSet &p)
Definition: FTFPCMS_BERT_EMM.cc:14
CMSEmStandardPhysicsLPM
Definition: CMSEmStandardPhysicsLPM.h:7
CMSHadronPhysicsFTFP_BERT
Definition: CMSHadronPhysicsFTFP_BERT.h:19
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
CMSEmStandardPhysicsLPM.h
FTFPCMS_BERT_EMM.h