CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
FTFPCMS_BERT.cc
Go to the documentation of this file.
1 #include "FTFPCMS_BERT.h"
4 
5 #include "G4EmStandardPhysics.hh"
6 #include "G4DecayPhysics.hh"
7 #include "G4EmExtraPhysics.hh"
8 #include "G4IonPhysics.hh"
9 #include "G4StoppingPhysics.hh"
10 #include "G4HadronElasticPhysics.hh"
11 
13  int ver = p.getUntrackedParameter<int>("Verbosity", 0);
14  bool emPhys = p.getUntrackedParameter<bool>("EMPhysics", true);
15  bool hadPhys = p.getUntrackedParameter<bool>("HadPhysics", true);
16  double minFTFP = p.getParameter<double>("EminFTFP") * CLHEP::GeV;
17  double maxBERT = p.getParameter<double>("EmaxBERT") * CLHEP::GeV;
18  double maxBERTpi = p.getParameter<double>("EmaxBERTpi") * CLHEP::GeV;
19  edm::LogVerbatim("PhysicsList") << "CMS Physics List FTFP_BERT: "
20  << "\n Flags for EM Physics: " << emPhys << "; Hadronic Physics: " << hadPhys
21  << "\n transition energy Bertini/FTFP from " << minFTFP / CLHEP::GeV << " to "
22  << maxBERT / CLHEP::GeV << ":" << maxBERTpi / CLHEP::GeV << " GeV";
23 
24  if (emPhys) {
25  // EM Physics
26  RegisterPhysics(new G4EmStandardPhysics(ver));
27 
28  // Synchroton Radiation & GN Physics
29  G4EmExtraPhysics* gn = new G4EmExtraPhysics(ver);
30  RegisterPhysics(gn);
31  }
32 
33  // Decays
34  RegisterPhysics(new G4DecayPhysics(ver));
35 
36  if (hadPhys) {
37  // Hadron Elastic scattering
38  RegisterPhysics(new G4HadronElasticPhysics(ver));
39 
40  // Hadron Physics
41  RegisterPhysics(new CMSHadronPhysicsFTFP_BERT(minFTFP, maxBERT, maxBERTpi, minFTFP, maxBERT));
42 
43  // Stopping Physics
44  RegisterPhysics(new G4StoppingPhysics(ver));
45 
46  // Ion Physics
47  RegisterPhysics(new G4IonPhysics(ver));
48  }
49 }
Log< level::Info, true > LogVerbatim
FTFPCMS_BERT(const edm::ParameterSet &p)
Definition: FTFPCMS_BERT.cc:12