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 #include <iomanip>
00042
00043 #include "HadronPhysicsQGSP_BERT_TRV.hh"
00044
00045 #include "globals.hh"
00046 #include "G4ios.hh"
00047 #include "G4SystemOfUnits.hh"
00048 #include "G4ParticleDefinition.hh"
00049 #include "G4ParticleTable.hh"
00050
00051 #include "G4MesonConstructor.hh"
00052 #include "G4BaryonConstructor.hh"
00053 #include "G4ShortLivedConstructor.hh"
00054
00055
00056 #include "G4PhysicsConstructorFactory.hh"
00057
00058 G4_DECLARE_PHYSCONSTR_FACTORY(HadronPhysicsQGSP_BERT_TRV);
00059
00060 HadronPhysicsQGSP_BERT_TRV::HadronPhysicsQGSP_BERT_TRV(G4int)
00061 : G4VPhysicsConstructor("hInelastic QGSP_BERT_TRV")
00062 , theNeutrons(0)
00063 , theLEPNeutron(0)
00064 , theQGSPNeutron(0)
00065 , theBertiniNeutron(0)
00066 , thePiK(0)
00067 , theLEPPiK(0)
00068 , theQGSPPiK(0)
00069 , theBertiniPiK(0)
00070 , thePro(0)
00071 , theLEPPro(0)
00072 , theQGSPPro(0)
00073 , theBertiniPro(0)
00074 , theMisc(0)
00075 , QuasiElastic(true)
00076 {}
00077
00078 HadronPhysicsQGSP_BERT_TRV::HadronPhysicsQGSP_BERT_TRV(const G4String& name, G4bool quasiElastic)
00079 : G4VPhysicsConstructor(name)
00080 , theNeutrons(0)
00081 , theLEPNeutron(0)
00082 , theQGSPNeutron(0)
00083 , theBertiniNeutron(0)
00084 , thePiK(0)
00085 , theLEPPiK(0)
00086 , theQGSPPiK(0)
00087 , theBertiniPiK(0)
00088 , thePro(0)
00089 , theLEPPro(0)
00090 , theQGSPPro(0)
00091 , theBertiniPro(0)
00092 , theMisc(0)
00093 , QuasiElastic(quasiElastic)
00094 {}
00095
00096 void HadronPhysicsQGSP_BERT_TRV::CreateModels()
00097 {
00098 theNeutrons=new G4NeutronBuilder;
00099 theNeutrons->RegisterMe(theQGSPNeutron=new G4QGSPNeutronBuilder(QuasiElastic));
00100 theQGSPNeutron->SetMinEnergy(10.0*GeV);
00101 theNeutrons->RegisterMe(theLEPNeutron=new G4LEPNeutronBuilder);
00102 theLEPNeutron->SetMinInelasticEnergy(9.5*GeV);
00103 theLEPNeutron->SetMaxInelasticEnergy(15*GeV);
00104
00105 theNeutrons->RegisterMe(theBertiniNeutron=new G4BertiniNeutronBuilder);
00106 theBertiniNeutron->SetMinEnergy(0.0*GeV);
00107 theBertiniNeutron->SetMaxEnergy(9.9*GeV);
00108
00109 thePro=new G4ProtonBuilder;
00110 thePro->RegisterMe(theQGSPPro=new G4QGSPProtonBuilder(QuasiElastic));
00111 theQGSPPro->SetMinEnergy(10.0*GeV);
00112 thePro->RegisterMe(theLEPPro=new G4LEPProtonBuilder);
00113 theLEPPro->SetMinEnergy(9.5*GeV);
00114 theLEPPro->SetMaxEnergy(15*GeV);
00115
00116 thePro->RegisterMe(theBertiniPro=new G4BertiniProtonBuilder);
00117 theBertiniPro->SetMaxEnergy(9.9*GeV);
00118
00119 thePiK=new G4PiKBuilder;
00120 thePiK->RegisterMe(theQGSPPiK=new G4QGSPPiKBuilder(QuasiElastic));
00121 theQGSPPiK->SetMinEnergy(10.0*GeV);
00122 thePiK->RegisterMe(theLEPPiK=new G4LEPPiKBuilder);
00123 theLEPPiK->SetMinEnergy(9.5*GeV);
00124 theLEPPiK->SetMaxEnergy(15*GeV);
00125
00126 thePiK->RegisterMe(theBertiniPiK=new G4BertiniPiKBuilder);
00127 theBertiniPiK->SetMaxEnergy(9.9*GeV);
00128
00129 theMisc=new G4MiscBuilder;
00130 }
00131
00132 HadronPhysicsQGSP_BERT_TRV::~HadronPhysicsQGSP_BERT_TRV()
00133 {
00134 delete theMisc;
00135 delete theQGSPNeutron;
00136 delete theLEPNeutron;
00137 delete theBertiniNeutron;
00138 delete theQGSPPro;
00139 delete theLEPPro;
00140 delete thePro;
00141 delete theBertiniPro;
00142 delete theQGSPPiK;
00143 delete theLEPPiK;
00144 delete theBertiniPiK;
00145 delete thePiK;
00146 }
00147
00148 void HadronPhysicsQGSP_BERT_TRV::ConstructParticle()
00149 {
00150 G4MesonConstructor pMesonConstructor;
00151 pMesonConstructor.ConstructParticle();
00152
00153 G4BaryonConstructor pBaryonConstructor;
00154 pBaryonConstructor.ConstructParticle();
00155
00156 G4ShortLivedConstructor pShortLivedConstructor;
00157 pShortLivedConstructor.ConstructParticle();
00158 }
00159
00160 #include "G4ProcessManager.hh"
00161 void HadronPhysicsQGSP_BERT_TRV::ConstructProcess()
00162 {
00163 CreateModels();
00164 theNeutrons->Build();
00165 thePro->Build();
00166 thePiK->Build();
00167 theMisc->Build();
00168 }
00169