source: trunk/source/physics_lists/lists/src/G4HadronInelasticQBBC.cc@ 1164

Last change on this file since 1164 was 988, checked in by garnier, 17 years ago

fichiers manquants
File size: 11.2 KB
Line 
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26// $Id: G4HadronInelasticQBBC.cc,v 1.16 2009/02/19 12:45:42 vnivanch Exp $
27// GEANT4 tag $Name: geant4-09-02-ref-02 $
28//
29//---------------------------------------------------------------------------
30//
31// ClassName: G4HadronInelasticQBBC
32//
33// Author: 11 April 2006 V. Ivanchenko
34//
35// Modified:
36// 05.07.2006 V.Ivanchenko fix problem of initialisation of HP
37// 15.04.2007 V.Ivanchenko include quasi-elastic and change FTF low energy
38//
39//----------------------------------------------------------------------------
40//
41
42#include "G4HadronInelasticQBBC.hh"
43
44#include "G4HadronInelasticProcess.hh"
45
46#include "G4ParticleDefinition.hh"
47#include "G4ProcessManager.hh"
48
49#include "G4MesonConstructor.hh"
50#include "G4BaryonConstructor.hh"
51
52#include "G4PiNuclearCrossSection.hh"
53#include "G4BGGNucleonInelasticXS.hh"
54#include "G4BGGPionInelasticXS.hh"
55
56#include "G4TheoFSGenerator.hh"
57#include "G4QStringChipsParticleLevelInterface.hh"
58#include "G4StringChipsParticleLevelInterface.hh"
59#include "G4GeneratorPrecompoundInterface.hh"
60#include "G4QGSMFragmentation.hh"
61#include "G4LundStringFragmentation.hh"
62#include "G4ExcitedStringDecay.hh"
63
64#include "G4BinaryCascade.hh"
65#include "G4HadronFissionProcess.hh"
66#include "G4HadronCaptureProcess.hh"
67#include "G4LFission.hh"
68#include "G4LCapture.hh"
69
70#include "G4NeutronHPInelastic.hh"
71#include "G4NeutronHPFission.hh"
72#include "G4NeutronHPCapture.hh"
73
74#include "G4UInelasticCrossSection.hh"
75
76G4HadronInelasticQBBC::G4HadronInelasticQBBC(const G4String& name,
77 G4int ver, G4bool ftf, G4bool bert, G4bool chips, G4bool hp, G4bool glauber)
78 : G4VPhysicsConstructor(name), verbose(ver), ftfFlag(ftf), bertFlag(bert),
79 chipsFlag(chips), hpFlag(hp), glFlag(glauber), wasActivated(false)
80{
81 if(verbose > 1) G4cout << "### HadronInelasticQBBC bertFlag= "
82 << bertFlag <<G4endl;
83 theHPXSecI = 0;
84 theHPXSecC = 0;
85 theHPXSecF = 0;
86 theCascade = 0;
87 preCompound = 0;
88 theCHIPSCascade = 0;
89 theQuasiElastic = 0;
90 theQGStringDecay = 0;
91 theQGStringModel = 0;
92 theFTFBStringDecay = 0;
93 theFTFBStringModel = 0;
94 theFTFCStringDecay = 0;
95 theFTFCStringModel = 0;
96}
97
98G4HadronInelasticQBBC::~G4HadronInelasticQBBC()
99{
100 delete theCascade;
101 delete preCompound;
102 delete theCHIPSCascade;
103 delete theQuasiElastic;
104 delete theQGStringDecay;
105 delete theQGStringModel;
106 delete theFTFBStringDecay;
107 delete theFTFCStringDecay;
108 delete theFTFBStringModel;
109 delete theFTFCStringModel;
110 delete theHPXSecI;
111 delete theHPXSecC;
112 delete theHPXSecF;
113}
114
115void G4HadronInelasticQBBC::ConstructParticle()
116{
117 G4MesonConstructor pMesonConstructor;
118 pMesonConstructor.ConstructParticle();
119
120 G4BaryonConstructor pBaryonConstructor;
121 pBaryonConstructor.ConstructParticle();
122}
123
124void G4HadronInelasticQBBC::ConstructProcess()
125{
126 if(wasActivated) return;
127 wasActivated = true;
128
129 if(verbose > 1) {
130 G4cout << "### HadronInelasticQBBC Construct Process"
131 << " ftfFlag= " << ftfFlag << " bertFlag= " << bertFlag
132 << G4endl;
133 }
134 G4double minEstring = 9.5*GeV;
135 G4double maxEcascade = 7.5*GeV;
136 G4double minFTF = 4.5*GeV;
137 G4double maxFTF = 25.*GeV;
138
139 //Binary
140 G4HadronicInteraction* theBIC = new G4BinaryCascade();
141 theBIC->SetMinEnergy(0.0);
142 theBIC->SetMaxEnergy(maxEcascade);
143
144 //Bertini
145 G4HadronicInteraction* theBERT = new G4CascadeInterface();
146 theBERT->SetMinEnergy(0.0);
147 theBERT->SetMaxEnergy(maxEcascade);
148
149 //CHIPS
150 G4HadronicInteraction* theCHIPS = new G4StringChipsParticleLevelInterface();
151 theCHIPS->SetMinEnergy(0.0);
152 theCHIPS->SetMaxEnergy(maxEcascade);
153
154 //QGS
155 theCascade = new G4BinaryCascade();
156 preCompound = new G4GeneratorPrecompoundInterface();
157
158 theCHIPSCascade = new G4QStringChipsParticleLevelInterface;
159 G4TheoFSGenerator* theQGSModel = new G4TheoFSGenerator("QGSP");
160 theQGStringModel = new G4QGSModel< G4QGSParticipants >;
161 theQGStringDecay = new G4ExcitedStringDecay(new G4QGSMFragmentation());
162 theQGStringModel->SetFragmentationModel(theQGStringDecay);
163 theQGSModel->SetTransport(preCompound);
164
165 theQuasiElastic = new G4QuasiElasticChannel();
166 theQGSModel->SetQuasiElasticChannel(theQuasiElastic);
167 theQGSModel->SetHighEnergyGenerator(theQGStringModel);
168 theQGSModel->SetMinEnergy(minEstring);
169 theQGSModel->SetMaxEnergy(100*TeV);
170
171 //FTFB
172 G4TheoFSGenerator* theFTFBModel = new G4TheoFSGenerator("FTFP");
173 theFTFBStringModel = new G4FTFModel();
174 theFTFBStringDecay = new G4ExcitedStringDecay(new G4LundStringFragmentation());
175 theFTFBStringModel->SetFragmentationModel(theFTFBStringDecay);
176
177 // theFTFBModel->SetTransport(theCascade);
178 theFTFBModel->SetTransport(preCompound);
179 theFTFBModel->SetHighEnergyGenerator(theFTFBStringModel);
180 theFTFBModel->SetMinEnergy(minFTF);
181 theFTFBModel->SetMaxEnergy(100*TeV);
182
183 //FTFP
184 G4TheoFSGenerator* theFTFCModel = new G4TheoFSGenerator("FTFP");
185 theFTFCStringModel = new G4FTFModel();
186 theFTFCStringDecay = new G4ExcitedStringDecay(new G4LundStringFragmentation());
187 theFTFCStringModel->SetFragmentationModel(theFTFCStringDecay);
188
189 theFTFCModel->SetTransport(preCompound);
190 theFTFCModel->SetHighEnergyGenerator(theFTFCStringModel);
191 theFTFCModel->SetMinEnergy(minFTF);
192 theFTFCModel->SetMaxEnergy(maxFTF);
193
194 theParticleIterator->reset();
195 while( (*theParticleIterator)() ) {
196 G4ParticleDefinition* particle = theParticleIterator->value();
197 G4String pname = particle->GetParticleName();
198 if(verbose > 1) G4cout << "### HadronInelasticQBBC: " << pname << G4endl;
199 if(pname == "anti_lambda" ||
200 pname == "anti_neutron" ||
201 pname == "anti_omega-" ||
202 pname == "anti_proton" ||
203 pname == "anti_sigma-" ||
204 pname == "anti_sigma+" ||
205 pname == "anti_xi-" ||
206 pname == "anti_xi0" ||
207 pname == "kaon-" ||
208 pname == "kaon+" ||
209 pname == "kaon0S" ||
210 pname == "kaon0L" ||
211 pname == "lambda" ||
212 pname == "neutron" ||
213 pname == "omega-" ||
214 pname == "pi-" ||
215 pname == "pi+" ||
216 pname == "proton" ||
217 pname == "sigma-" ||
218 pname == "sigma+" ||
219 pname == "xi-" ||
220 pname == "xi0") {
221
222 G4ProcessManager* pmanager = particle->GetProcessManager();
223 G4HadronInelasticProcess* hp =
224 new G4HadronInelasticProcess("hInelastic", particle);
225 pmanager->AddDiscreteProcess(hp);
226
227 if(pname == "proton") {
228 hp->AddDataSet(&theXSecP);
229
230 hp->RegisterMe(theQGSModel);
231 hp->RegisterMe(theFTFCModel);
232 //if(ftfFlag) hp->RegisterMe(theFTFCModel);
233 //else hp->RegisterMe(theQGSModel);
234
235 if(bertFlag) hp->RegisterMe(theBERT);
236 else hp->RegisterMe(theBIC);
237
238 if(glFlag)
239 hp->AddDataSet(new G4BGGNucleonInelasticXS(particle));
240
241 } else if(pname == "neutron") {
242 hp->AddDataSet(&theXSecN);
243 hp->RegisterMe(theQGSModel);
244 hp->RegisterMe(theFTFCModel);
245 //if(ftfFlag) hp->RegisterMe(theFTFCModel);
246 //else hp->RegisterMe(theQGSModel);
247
248 G4HadronCaptureProcess* theNeutronCapture =
249 new G4HadronCaptureProcess("nCapture");
250 G4HadronFissionProcess* theNeutronFission =
251 new G4HadronFissionProcess("nFission");
252 pmanager->AddDiscreteProcess(theNeutronCapture);
253 pmanager->AddDiscreteProcess(theNeutronFission);
254
255 G4double emin = 0.0;
256 if(hpFlag) {
257 emin = 19.5*MeV;
258 theHPXSecI = new G4NeutronHPInelasticData;
259 theHPXSecC = new G4NeutronHPCaptureData;
260 theHPXSecF = new G4NeutronHPFissionData;
261 hp->AddDataSet(theHPXSecI);
262 theNeutronCapture->AddDataSet(theHPXSecC);
263 theNeutronFission->AddDataSet(theHPXSecF);
264 hp->RegisterMe(new G4NeutronHPInelastic());
265 theNeutronCapture->RegisterMe(new G4NeutronHPCapture());
266 theNeutronFission->RegisterMe(new G4NeutronHPFission());
267 }
268
269 G4HadronicInteraction* theB;
270 if(bertFlag) theB = new G4CascadeInterface();
271 else theB = new G4BinaryCascade();
272 theB->SetMinEnergy(emin);
273 theB->SetMaxEnergy(maxEcascade);
274 hp->RegisterMe(theB);
275
276 if(glFlag)
277 hp->AddDataSet(new G4BGGNucleonInelasticXS(particle));
278
279 G4HadronicInteraction* theC = new G4LCapture();
280 theC->SetMinEnergy(emin);
281 theNeutronCapture->RegisterMe(theC);
282
283 G4HadronicInteraction* theF = new G4LFission();
284 theF->SetMinEnergy(emin);
285 theNeutronFission->RegisterMe(theF);
286
287 } else if(pname == "pi-" || pname == "pi+") {
288 hp->AddDataSet(&thePiCross);
289 hp->RegisterMe(theQGSModel);
290 hp->RegisterMe(theFTFCModel);
291 //if(ftfFlag) hp->RegisterMe(theFTFCModel);
292 //else hp->RegisterMe(theQGSModel);
293
294 hp->RegisterMe(theBERT);
295 //if(bertFlag) hp->RegisterMe(theBERT);
296 //else hp->RegisterMe(theBIC);
297
298 if(glFlag)
299 hp->AddDataSet(new G4BGGPionInelasticXS(particle));
300
301 } else if(pname == "kaon-" ||
302 pname == "kaon+" ||
303 pname == "kaon0S" ||
304 pname == "kaon0L") {
305 hp->RegisterMe(theFTFBModel);
306 hp->RegisterMe(theBERT);
307 //hp->AddDataSet(new G4UInelasticCrossSection(particle));
308
309 } else if(pname == "lambda" ||
310 pname == "sigma-" ||
311 pname == "sigma+" ||
312 pname == "xi-" ||
313 pname == "xi0") {
314
315 hp->RegisterMe(theFTFBModel);
316 hp->RegisterMe(theBERT);
317 //hp->AddDataSet(new G4UInelasticCrossSection(particle));
318
319 } else if(pname == "anti_proton" || pname == "anti_neutron") {
320 hp->RegisterMe(theFTFBModel);
321 hp->RegisterMe(theCHIPS);
322 //hp->AddDataSet(new G4UInelasticCrossSection(particle));
323
324 } else {
325 hp->RegisterMe(theFTFBModel);
326 hp->RegisterMe(theCHIPS);
327 //hp->AddDataSet(new G4UInelasticCrossSection(particle));
328 }
329
330 if(verbose > 1)
331 G4cout << "### HadronInelasticQBBC: " << hp->GetProcessName()
332 << " added for " << pname << G4endl;
333 }
334 }
335}
Note: See TracBrowser for help on using the repository browser.