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_BIC.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_BIC);
00060
00061 HadronPhysicsQGSP_BIC::HadronPhysicsQGSP_BIC(G4int)
00062 : G4VPhysicsConstructor("hInelastic QGSP_BIC")
00063 , theNeutrons(0)
00064 , theLEPNeutron(0)
00065 , theQGSPNeutron(0)
00066 , theBinaryNeutron(0)
00067 , thePiK(0)
00068 , theLEPPiK(0)
00069 , theQGSPPiK(0)
00070 , thePro(0)
00071 , theLEPPro(0)
00072 , theQGSPPro(0)
00073 , theBinaryPro(0)
00074 , theMisc(0)
00075 , QuasiElastic(true)
00076 {}
00077
00078 HadronPhysicsQGSP_BIC::HadronPhysicsQGSP_BIC(const G4String& name, G4bool quasiElastic)
00079 : G4VPhysicsConstructor(name)
00080 , theNeutrons(0)
00081 , theLEPNeutron(0)
00082 , theQGSPNeutron(0)
00083 , theBinaryNeutron(0)
00084 , thePiK(0)
00085 , theLEPPiK(0)
00086 , theQGSPPiK(0)
00087 , thePro(0)
00088 , theLEPPro(0)
00089 , theQGSPPro(0)
00090 , theBinaryPro(0)
00091 , theMisc(0)
00092 , QuasiElastic(quasiElastic)
00093 {}
00094
00095 void HadronPhysicsQGSP_BIC::CreateModels()
00096 {
00097 theNeutrons=new G4NeutronBuilder;
00098
00099 theNeutrons->RegisterMe(theQGSPNeutron=new G4QGSPNeutronBuilder(QuasiElastic));
00100 theNeutrons->RegisterMe(theLEPNeutron=new G4LEPNeutronBuilder);
00101 theLEPNeutron->SetMinInelasticEnergy(9.5*GeV);
00102 theLEPNeutron->SetMaxInelasticEnergy(25*GeV);
00103
00104 theNeutrons->RegisterMe(theBinaryNeutron=new G4BinaryNeutronBuilder);
00105 theBinaryNeutron->SetMaxEnergy(9.9*GeV);
00106
00107 thePro=new G4ProtonBuilder;
00108 thePro->RegisterMe(theQGSPPro=new G4QGSPProtonBuilder(QuasiElastic));
00109 thePro->RegisterMe(theLEPPro=new G4LEPProtonBuilder);
00110 theLEPPro->SetMinEnergy(9.5*GeV);
00111 theLEPPro->SetMaxEnergy(25*GeV);
00112
00113 thePro->RegisterMe(theBinaryPro=new G4BinaryProtonBuilder);
00114 theBinaryPro->SetMaxEnergy(9.9*GeV);
00115
00116 thePiK=new G4PiKBuilder;
00117 thePiK->RegisterMe(theQGSPPiK=new G4QGSPPiKBuilder(QuasiElastic));
00118 thePiK->RegisterMe(theLEPPiK=new G4LEPPiKBuilder);
00119 theLEPPiK->SetMaxEnergy(25*GeV);
00120
00121 theMisc=new G4MiscBuilder;
00122 }
00123
00124 HadronPhysicsQGSP_BIC::~HadronPhysicsQGSP_BIC()
00125 {
00126 delete theMisc;
00127 delete theQGSPNeutron;
00128 delete theLEPNeutron;
00129 delete theBinaryNeutron;
00130 delete theQGSPPro;
00131 delete theLEPPro;
00132 delete thePro;
00133 delete theBinaryPro;
00134 delete theQGSPPiK;
00135 delete theLEPPiK;
00136 delete thePiK;
00137 }
00138
00139 void HadronPhysicsQGSP_BIC::ConstructParticle()
00140 {
00141 G4MesonConstructor pMesonConstructor;
00142 pMesonConstructor.ConstructParticle();
00143
00144 G4BaryonConstructor pBaryonConstructor;
00145 pBaryonConstructor.ConstructParticle();
00146
00147 G4ShortLivedConstructor pShortLivedConstructor;
00148 pShortLivedConstructor.ConstructParticle();
00149 }
00150
00151 #include "G4ProcessManager.hh"
00152 void HadronPhysicsQGSP_BIC::ConstructProcess()
00153 {
00154 CreateModels();
00155 theNeutrons->Build();
00156 thePro->Build();
00157 thePiK->Build();
00158 theMisc->Build();
00159 }
00160