CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
FTFPCMS_BERT_EML.cc
Go to the documentation of this file.
1 #include "FTFPCMS_BERT_EML.hh"
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 
14 #include "G4DataQuestionaire.hh"
15 #include "G4HadronPhysicsFTFP_BERT.hh"
16 
17 FTFPCMS_BERT_EML::FTFPCMS_BERT_EML(G4LogicalVolumeToDDLogicalPartMap& map,
19  sim::ChordFinderSetter *chordFinderSetter_,
20  const edm::ParameterSet & p) : PhysicsList(map, table_, chordFinderSetter_, p) {
21 
22  G4DataQuestionaire it(photon);
23 
24  int ver = p.getUntrackedParameter<int>("Verbosity",0);
25  bool emPhys = p.getUntrackedParameter<bool>("EMPhysics",true);
26  bool hadPhys = p.getUntrackedParameter<bool>("HadPhysics",true);
27  bool tracking= p.getParameter<bool>("TrackingCut");
28  double timeLimit = p.getParameter<double>("MaxTrackTime")*ns;
29  edm::LogInfo("PhysicsList") << "You are using the simulation engine: "
30  << "FTFP_BERT_EML: \n Flags for EM Physics "
31  << emPhys << ", for Hadronic Physics "
32  << hadPhys << " and tracking cut " << tracking
33  << " t(ns)= " << timeLimit/ns;
34 
35  if (emPhys) {
36  // EM Physics
37  RegisterPhysics( new CMSEmStandardPhysics95msc93("EM standard msc93",ver,""));
38 
39  // Synchroton Radiation & GN Physics
40  G4EmExtraPhysics* gn = new G4EmExtraPhysics(ver);
41  RegisterPhysics(gn);
42  }
43 
44  // Decays
45  this->RegisterPhysics( new G4DecayPhysics(ver) );
46 
47  if (hadPhys) {
48  G4HadronicProcessStore::Instance()->SetVerbose(ver);
49 
50  // Hadron Elastic scattering
51  RegisterPhysics( new G4HadronElasticPhysics(ver));
52 
53  // Hadron Physics
54  RegisterPhysics( new G4HadronPhysicsFTFP_BERT(ver));
55 
56  // Stopping Physics
57  RegisterPhysics( new G4StoppingPhysics(ver));
58 
59  // Ion Physics
60  RegisterPhysics( new G4IonPhysics(ver));
61 
62  // Neutron tracking cut
63  if (tracking) {
64  G4NeutronTrackingCut* ncut= new G4NeutronTrackingCut(ver);
65  ncut->SetTimeLimit(timeLimit);
66  RegisterPhysics(ncut);
67  }
68  }
69 
70  // Monopoles
71  RegisterPhysics( new CMSMonopolePhysics(table_,chordFinderSetter_,p));
72 }
73 
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
HepPDT::ParticleDataTable ParticleDataTable
Table table_