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