00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042 #include <iomanip>
00043
00044 #include "HadronPhysicsQGSP.hh"
00045
00046 #include "globals.hh"
00047 #include "G4ios.hh"
00048 #include "G4SystemOfUnits.hh"
00049 #include "G4ParticleDefinition.hh"
00050 #include "G4ParticleTable.hh"
00051
00052 #include "G4MesonConstructor.hh"
00053 #include "G4BaryonConstructor.hh"
00054 #include "G4ShortLivedConstructor.hh"
00055
00056
00057 #include "G4PhysicsConstructorFactory.hh"
00058
00059 G4_DECLARE_PHYSCONSTR_FACTORY(HadronPhysicsQGSP);
00060
00061 HadronPhysicsQGSP::HadronPhysicsQGSP(G4int)
00062 : G4VPhysicsConstructor("hInelastic QGSP")
00063 , theNeutrons(0)
00064 , theLEPNeutron(0)
00065 , theQGSPNeutron(0)
00066 , thePiK(0)
00067 , theLEPPiK(0)
00068 , theQGSPPiK(0)
00069 , thePro(0)
00070 , theLEPPro(0)
00071 , theQGSPPro(0)
00072 , theMisc(0)
00073 , QuasiElastic(true)
00074 , ProjectileDiffraction(false)
00075 {
00076 }
00077
00078 HadronPhysicsQGSP::HadronPhysicsQGSP(const G4String& name, G4bool quasiElastic)
00079 : G4VPhysicsConstructor(name)
00080 , theNeutrons(0)
00081 , theLEPNeutron(0)
00082 , theQGSPNeutron(0)
00083 , thePiK(0)
00084 , theLEPPiK(0)
00085 , theQGSPPiK(0)
00086 , thePro(0)
00087 , theLEPPro(0)
00088 , theQGSPPro(0)
00089 , theMisc(0)
00090 , QuasiElastic(quasiElastic)
00091 , ProjectileDiffraction(false)
00092 {
00093 }
00094
00095 void HadronPhysicsQGSP::CreateModels()
00096 {
00097 theNeutrons=new G4NeutronBuilder;
00098 theQGSPNeutron=new G4QGSPNeutronBuilder(QuasiElastic, ProjectileDiffraction);
00099 theNeutrons->RegisterMe(theQGSPNeutron);
00100 theNeutrons->RegisterMe(theLEPNeutron=new G4LEPNeutronBuilder);
00101 theLEPNeutron->SetMaxInelasticEnergy(25*GeV);
00102
00103 thePro=new G4ProtonBuilder;
00104 theQGSPPro=new G4QGSPProtonBuilder(QuasiElastic, ProjectileDiffraction);
00105 thePro->RegisterMe(theQGSPPro);
00106 thePro->RegisterMe(theLEPPro=new G4LEPProtonBuilder);
00107 theLEPPro->SetMaxEnergy(25*GeV);
00108
00109 thePiK=new G4PiKBuilder;
00110 theQGSPPiK=new G4QGSPPiKBuilder(QuasiElastic, false);
00111 thePiK->RegisterMe(theQGSPPiK);
00112 thePiK->RegisterMe(theLEPPiK=new G4LEPPiKBuilder);
00113 theLEPPiK->SetMaxEnergy(25*GeV);
00114
00115 theMisc=new G4MiscBuilder;
00116 }
00117
00118 HadronPhysicsQGSP::~HadronPhysicsQGSP()
00119 {
00120 delete theMisc;
00121 delete theQGSPNeutron;
00122 delete theLEPNeutron;
00123 delete theQGSPPro;
00124 delete theLEPPro;
00125 delete thePro;
00126 delete theQGSPPiK;
00127 delete theLEPPiK;
00128 delete thePiK;
00129 }
00130
00131
00132 void HadronPhysicsQGSP::ConstructParticle()
00133 {
00134 G4MesonConstructor pMesonConstructor;
00135 pMesonConstructor.ConstructParticle();
00136
00137 G4BaryonConstructor pBaryonConstructor;
00138 pBaryonConstructor.ConstructParticle();
00139
00140 G4ShortLivedConstructor pShortLivedConstructor;
00141 pShortLivedConstructor.ConstructParticle();
00142 }
00143
00144 #include "G4ProcessManager.hh"
00145 void HadronPhysicsQGSP::ConstructProcess()
00146 {
00147 CreateModels();
00148 theNeutrons->Build();
00149 thePro->Build();
00150 thePiK->Build();
00151 theMisc->Build();
00152 }
00153