CMSSW Reference Manual

00001 #include "QBBCCMS.hh"
00002 #include "SimG4Core/PhysicsLists/interface/CMSEmStandardPhysics.h"
00003 #include "FWCore/MessageLogger/interface/MessageLogger.h"
00004 
00005 #include "G4DecayPhysics.hh"
00006 #include "G4EmExtraPhysics.hh"
00007 #include "G4QStoppingPhysics.hh"
00008 #include "G4LHEPStoppingPhysics.hh" 
00009 
00010 #include "G4DataQuestionaire.hh"
00011 #include "G4HadronInelasticQBBC.hh"
00012 #include "G4HadronInelasticQLHEP.hh"
00013 #include "G4HadronElasticPhysics.hh"
00014 #include "G4HadronDElasticPhysics.hh"
00015 #include "G4HadronHElasticPhysics.hh"
00016 #include "G4IonBinaryCascadePhysics.hh"
00017 #include "G4IonPhysics.hh"
00018 #include "G4NeutronTrackingCut.hh"
00019 
00020 QBBCCMS::QBBCCMS(G4LogicalVolumeToDDLogicalPartMap& map, 
00021                  const edm::ParameterSet & p) : PhysicsList(map, p) {
00022 
00023   G4DataQuestionaire it(photon);
00024   
00025   int  ver     = p.getUntrackedParameter<int>("Verbosity",0);
00026   bool emPhys  = p.getUntrackedParameter<bool>("EMPhysics",true);
00027   bool hadPhys = p.getUntrackedParameter<bool>("HadPhysics",true);
00028   bool ftf     = p.getUntrackedParameter<bool>("FlagFTF",false);
00029   bool bert    = p.getUntrackedParameter<bool>("FlagBERT",false);
00030   bool chips   = p.getUntrackedParameter<bool>("FlagCHIPS",false);
00031   bool hp      = p.getUntrackedParameter<bool>("FlagHP",false);
00032   bool glauber = p.getUntrackedParameter<bool>("FlagGlauber",false);
00033   edm::LogInfo("PhysicsList") << "You are using the simulation engine: "
00034                               << "QBBC 3.1 with Flags for EM Physics "
00035                               << emPhys << " and for Hadronic Physics "
00036                               << hadPhys << " Flags for FTF " << ftf
00037                               << " BERT " << bert << " CHIPS " << chips
00038                               << " HP " << hp << " Glauber " << glauber <<"\n";
00039 
00040   if (emPhys) {
00041     // EM Physics
00042     RegisterPhysics( new CMSEmStandardPhysics("standard EM",ver));
00043 
00044     // Synchroton Radiation & GN Physics
00045     RegisterPhysics(new G4EmExtraPhysics("extra EM"));
00046   }
00047 
00048   // Decays
00049   RegisterPhysics(new G4DecayPhysics("decay",ver));
00050 
00051   if (hadPhys) {
00052     // Hadron Elastic scattering
00053     RegisterPhysics(new G4HadronElasticPhysics("hElastic",ver,false,true));
00054 
00055     // Hadron Physics
00056     RegisterPhysics( new G4HadronInelasticQBBC("inelastic", ver, ftf,
00057                                                bert, chips, hp, glauber));
00058 
00059     // Stopping Physics
00060     RegisterPhysics(new G4QStoppingPhysics("stopping",ver));
00061 
00062     // Ion Physics
00063     RegisterPhysics(new G4IonBinaryCascadePhysics("ionBIC"));
00064 
00065     // Neutron tracking cut
00066     RegisterPhysics( new G4NeutronTrackingCut("Neutron tracking cut", ver));
00067   }
00068 }
00069
QBBCCMS.cc
