73 G4HadronPhysicsQGSP_BERT::tpdata = 0;
118 , QuasiElastic(quasiElastic)
122 void G4HadronPhysicsQGSP_BERT::CreateModels()
124 G4bool quasiElasticFTF=
false;
125 G4bool quasiElasticQGS=
true;
134 tpdata->theNeutrons->RegisterMe(tpdata->theQGSPNeutron=
new G4QGSPNeutronBuilder(quasiElasticQGS));
135 tpdata->theQGSPNeutron->SetMinEnergy(minQGSP);
136 tpdata->theNeutrons->RegisterMe(tpdata->theFTFPNeutron=
new G4FTFPNeutronBuilder(quasiElasticFTF));
137 tpdata->theFTFPNeutron->SetMinEnergy(minFTFP);
138 tpdata->theFTFPNeutron->SetMaxEnergy(maxFTFP);
141 tpdata->theBertiniNeutron->SetMinEnergy(minBERT);
142 tpdata->theBertiniNeutron->SetMaxEnergy(maxBERT);
146 tpdata->theQGSPPro->SetMinEnergy(minQGSP);
148 tpdata->theFTFPPro->SetMinEnergy(minFTFP);
149 tpdata->theFTFPPro->SetMaxEnergy(maxFTFP);
151 tpdata->theBertiniPro->SetMaxEnergy(maxBERT);
154 tpdata->thePiK->RegisterMe(tpdata->theQGSPPiK=
new G4QGSPPiKBuilder(quasiElasticQGS));
155 tpdata->theQGSPPiK->SetMinEnergy(minQGSP);
156 tpdata->thePiK->RegisterMe(tpdata->theFTFPPiK=
new G4FTFPPiKBuilder(quasiElasticFTF));
157 tpdata->theFTFPPiK->SetMinEnergy(minFTFP);
158 tpdata->theFTFPPiK->SetMaxEnergy(maxFTFP);
160 tpdata->theBertiniPiK->SetMaxEnergy(maxBERT);
170 delete tpdata->theBertiniNeutron;
171 delete tpdata->theQGSPNeutron;
172 delete tpdata->theFTFPNeutron;
173 delete tpdata->theNeutrons;
174 delete tpdata->theBertiniPiK;
175 delete tpdata->theQGSPPiK;
176 delete tpdata->theFTFPPiK;
177 delete tpdata->thePiK;
178 delete tpdata->theBertiniPro;
179 delete tpdata->theQGSPPro;
180 delete tpdata->theFTFPPro;
181 delete tpdata->thePro;
182 delete tpdata->theFTFPAntiBaryon;
183 delete tpdata->theAntiBaryon;
184 delete tpdata->theHyperon;
185 delete tpdata->xsNeutronInelasticXS;
186 delete tpdata->xsNeutronCaptureXS;
188 delete tpdata; tpdata = 0;
209 if ( tpdata == 0 ) tpdata =
new ThreadPrivate;
211 tpdata->theNeutrons->Build();
212 tpdata->thePro->Build();
213 tpdata->thePiK->Build();
214 tpdata->theHyperon->Build();
215 tpdata->theAntiBaryon->Build();
224 for (
size_t i=0; i < static_cast<size_t>(pv->
size()); ++i ) {
225 if (
fCapture == ((*pv)[i])->GetProcessSubType() ) {
234 capture->
AddDataSet(tpdata->xsNeutronCaptureXS);
G4int AddDiscreteProcess(G4VProcess *aProcess, G4int ord=ordDefault)
G4_DECLARE_PHYSCONSTR_FACTORY(G4HadronPhysicsQGSP_BERT)
static void ConstructParticle()
static void ConstructParticle()
static void ConstructParticle()
G4ProcessManager * GetProcessManager() const
static void ConstructParticle()
void RegisterMe(G4HadronicInteraction *a)
virtual void ConstructParticle()
static G4HadronicProcess * FindInelasticProcess(const G4ParticleDefinition *)
void AddDataSet(G4VCrossSectionDataSet *aDataSet)
G4HadronPhysicsQGSP_BERT(G4int verbose=1)
static G4Neutron * Neutron()
virtual void ConstructProcess()
virtual ~G4HadronPhysicsQGSP_BERT()
G4ProcessVector * GetProcessList() const