CMS 3D CMS Logo

FTFPCMS_BERT_EMZ.cc
Go to the documentation of this file.
1 #include "FTFPCMS_BERT_EMZ.h"
4 
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"
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_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 "
26  << minFTFP / CLHEP::GeV << " to " << maxBERT / CLHEP::GeV << ":"
27  << maxBERTpi / CLHEP::GeV << " GeV";
28 
29  if (emPhys) {
30  // EM Physics
31  RegisterPhysics(new G4EmStandardPhysics_option4(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 }
FTFPCMS_BERT_EMZ::FTFPCMS_BERT_EMZ
FTFPCMS_BERT_EMZ(const edm::ParameterSet &p)
Definition: FTFPCMS_BERT_EMZ.cc:14
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
FTFPCMS_BERT_EMZ.h
GeV
const double GeV
Definition: MathUtil.h:16
CMSHadronPhysicsFTFP_BERT
Definition: CMSHadronPhysicsFTFP_BERT.h:19
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128