ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4HadronPhysicsFTFQGSP_BERT.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4HadronPhysicsFTFQGSP_BERT.cc
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 //
27 //---------------------------------------------------------------------------
28 // Author: Alberto Ribon
29 // Date: October 2017
30 //
31 // Hadron physics for the new, experimental physics list FTFQGSP_BERT,
32 // with QGS fragmentation of strings, instead of the Lund string
33 // fragmentation. Note that the string excitation is still done with FTF,
34 // exactly as for FTFP_BERT.
35 // Given that it is an experimental, and perhaps temporary, new type of
36 // hadron physics, corresponding builders are not created and everything
37 // is implemented directly in this class.
38 //----------------------------------------------------------------------------
39 //
40 #include <iomanip>
41 
42 #include "G4HadronPhysicsFTFQGSP_BERT.hh"
43 
44 #include "globals.hh"
45 #include "G4ios.hh"
46 #include "G4SystemOfUnits.hh"
47 #include "G4ParticleDefinition.hh"
48 #include "G4ParticleTable.hh"
49 
50 #include "G4MesonConstructor.hh"
51 #include "G4BaryonConstructor.hh"
52 #include "G4ShortLivedConstructor.hh"
53 
54 #include "G4ComponentGGHadronNucleusXsc.hh"
55 #include "G4CrossSectionDataSetRegistry.hh"
56 
57 #include "G4HadronCaptureProcess.hh"
58 #include "G4NeutronRadCapture.hh"
59 #include "G4NeutronInelasticXS.hh"
60 #include "G4NeutronCaptureXS.hh"
61 
62 #include "G4ProcessManager.hh"
63 #include "G4BGGNucleonInelasticXS.hh"
64 #include "G4BGGPionInelasticXS.hh"
65 #include "G4CrossSectionPairGG.hh"
66 #include "G4ComponentAntiNuclNuclearXS.hh"
67 #include "G4CrossSectionInelastic.hh"
68 
69 #include "G4PhysListUtil.hh"
70 
71 #include "G4HadronicParameters.hh"
72 
73 // factory
74 #include "G4PhysicsConstructorFactory.hh"
75 //
76 G4_DECLARE_PHYSCONSTR_FACTORY(G4HadronPhysicsFTFQGSP_BERT);
77 
78 G4HadronPhysicsFTFQGSP_BERT::G4HadronPhysicsFTFQGSP_BERT(G4int)
79  : G4VPhysicsConstructor("hInelastic FTFQGSP_BERT")
80  , theNeutronCaptureModel(0)
81  , thePreEquilib(0)
82  , theCascade(0)
83  , theStringModel(0)
84  , theStringDecay(0)
85  , theQGSMFragmentation(0)
86  , theHandler(0)
87  , theModel1(0)
88  , theModel2(0)
89  , theModel3(0)
90  , theBertini1(0)
91  , theBertini2(0)
92  , theNeutronCaptureProcess(0)
93  , theNeutronInelastic(0)
94  , theProtonInelastic(0)
95  , thePionMinusInelastic(0)
96  , thePionPlusInelastic(0)
97  , theKaonMinusInelastic(0)
98  , theKaonPlusInelastic(0)
99  , theKaonZeroLInelastic(0)
100  , theKaonZeroSInelastic(0)
101  , theLambdaInelastic(0)
102  , theAntiLambdaInelastic(0)
103  , theSigmaMinusInelastic(0)
104  , theAntiSigmaMinusInelastic(0)
105  , theSigmaPlusInelastic(0)
106  , theAntiSigmaPlusInelastic(0)
107  , theXiZeroInelastic(0)
108  , theAntiXiZeroInelastic(0)
109  , theXiMinusInelastic(0)
110  , theAntiXiMinusInelastic(0)
111  , theOmegaMinusInelastic(0)
112  , theAntiOmegaMinusInelastic(0)
113  , theAntiProtonInelastic(0)
114  , theAntiNeutronInelastic(0)
115  , theAntiDeuteronInelastic(0)
116  , theAntiTritonInelastic(0)
117  , theAntiHe3Inelastic(0)
118  , theAntiAlphaInelastic(0)
119  , thePiXS(0)
120  , theKaonAndHyperonXS(0)
121  , theAntiNucleonXS(0)
122  , theNeutronInelasticXS(0)
123  , theNeutronCaptureXS(0)
124 {}
125 
126 
127 G4HadronPhysicsFTFQGSP_BERT::G4HadronPhysicsFTFQGSP_BERT(const G4String& name, G4bool /* quasiElastic */)
128  : G4VPhysicsConstructor(name)
129  , theNeutronCaptureModel(0)
130  , thePreEquilib(0)
131  , theCascade(0)
132  , theStringModel(0)
133  , theStringDecay(0)
134  , theQGSMFragmentation(0)
135  , theHandler(0)
136  , theModel1(0)
137  , theModel2(0)
138  , theModel3(0)
139  , theBertini1(0)
140  , theBertini2(0)
141  , theNeutronCaptureProcess(0)
142  , theNeutronInelastic(0)
143  , theProtonInelastic(0)
144  , thePionMinusInelastic(0)
145  , thePionPlusInelastic(0)
146  , theKaonMinusInelastic(0)
147  , theKaonPlusInelastic(0)
148  , theKaonZeroLInelastic(0)
149  , theKaonZeroSInelastic(0)
150  , theLambdaInelastic(0)
151  , theAntiLambdaInelastic(0)
152  , theSigmaMinusInelastic(0)
153  , theAntiSigmaMinusInelastic(0)
154  , theSigmaPlusInelastic(0)
155  , theAntiSigmaPlusInelastic(0)
156  , theXiZeroInelastic(0)
157  , theAntiXiZeroInelastic(0)
158  , theXiMinusInelastic(0)
159  , theAntiXiMinusInelastic(0)
160  , theOmegaMinusInelastic(0)
161  , theAntiOmegaMinusInelastic(0)
162  , theAntiProtonInelastic(0)
163  , theAntiNeutronInelastic(0)
164  , theAntiDeuteronInelastic(0)
165  , theAntiTritonInelastic(0)
166  , theAntiHe3Inelastic(0)
167  , theAntiAlphaInelastic(0)
168  , thePiXS(0)
169  , theKaonAndHyperonXS(0)
170  , theAntiNucleonXS(0)
171  , theNeutronInelasticXS(0)
172  , theNeutronCaptureXS(0)
173 {}
174 
175 
176 void G4HadronPhysicsFTFQGSP_BERT::CreateModels()
177 {
178  G4double minFTFP = G4HadronicParameters::Instance()->GetMinEnergyTransitionFTF_Cascade();
179  G4double maxBERT = G4HadronicParameters::Instance()->GetMaxEnergyTransitionFTF_Cascade();
180  G4cout << " FTFQGSP_BERT : similar to FTFP_BERT but with" << G4endl
181  << " QGS string fragmentation (instead of Lund string fragmentation)." << G4endl;
182 
183  theStringModel = new G4FTFModel;
184 
185  theStringDecay = new G4ExcitedStringDecay( theQGSMFragmentation = new G4QGSMFragmentation );
186  theStringModel->SetFragmentationModel( theStringDecay );
187  thePreEquilib = new G4PreCompoundModel( theHandler = new G4ExcitationHandler );
188  theCascade = new G4GeneratorPrecompoundInterface( thePreEquilib );
189 
190  // FTF for neutrons, protons, pions, and kaons
191  theModel1 = new G4TheoFSGenerator( "FTFP" );
192  theModel1->SetMinEnergy( minFTFP );
193  theModel1->SetMaxEnergy( G4HadronicParameters::Instance()->GetMaxEnergy() );
194  theModel1->SetTransport( theCascade );
195  theModel1->SetHighEnergyGenerator( theStringModel );
196 
197  // BERT for neutrons, protons, pions, and kaons
198  theBertini1 = new G4CascadeInterface;
199  theBertini1->SetMinEnergy( 0.0 );
200  theBertini1->SetMaxEnergy( maxBERT );
201 
202  // FTF for hyperons
203  theModel2 = new G4TheoFSGenerator( "FTFP" );
204  theModel2->SetMinEnergy( G4HadronicParameters::Instance()->GetMinEnergyTransitionFTF_Cascade() );
205  theModel2->SetMaxEnergy( G4HadronicParameters::Instance()->GetMaxEnergy() );
206  theModel2->SetTransport( theCascade );
207  theModel2->SetHighEnergyGenerator( theStringModel );
208 
209  // BERT for hyperons
210  theBertini2 = new G4CascadeInterface;
211  theBertini2->SetMinEnergy( 0.0 );
212  theBertini2->SetMaxEnergy( G4HadronicParameters::Instance()->GetMaxEnergyTransitionFTF_Cascade() );
213 
214  // FTF for Antibaryons
215  theModel3 = new G4TheoFSGenerator( "FTFP" );
216  theModel3->SetMinEnergy( 0.0 );
217  theModel3->SetMaxEnergy( G4HadronicParameters::Instance()->GetMaxEnergy() );
218  theModel3->SetTransport( theCascade );
219  theModel3->SetHighEnergyGenerator( theStringModel );
220 
221  // Neutron Capture
222  theNeutronCaptureModel = new G4NeutronRadCapture;
223  theNeutronCaptureModel->SetMinEnergy( 0.0 );
224  theNeutronCaptureModel->SetMaxEnergy( G4HadronicParameters::Instance()->GetMaxEnergy() );
225 
226  // Cross sections
227  thePiXS = new G4BGGPionInelasticXS( G4PionMinus::Definition() );
228  theAntiNucleonXS = new G4CrossSectionInelastic( new G4ComponentAntiNuclNuclearXS );
229  theKaonAndHyperonXS = new G4CrossSectionInelastic( new G4ComponentGGHadronNucleusXsc );
230  theNeutronInelasticXS = G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet( G4NeutronInelasticXS::Default_Name() );
231  theNeutronCaptureXS = G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet( G4NeutronCaptureXS::Default_Name() );
232 }
233 
234 
235 G4HadronPhysicsFTFQGSP_BERT::~G4HadronPhysicsFTFQGSP_BERT()
236 {
237  delete theStringDecay;
238  delete theStringModel;
239  delete thePreEquilib;
240  delete theCascade;
241  delete theQGSMFragmentation;
242 }
243 
244 
245 void G4HadronPhysicsFTFQGSP_BERT::ConstructParticle()
246 {
247  G4MesonConstructor pMesonConstructor;
248  pMesonConstructor.ConstructParticle();
249 
250  G4BaryonConstructor pBaryonConstructor;
251  pBaryonConstructor.ConstructParticle();
252 
253  G4ShortLivedConstructor pShortLivedConstructor;
254  pShortLivedConstructor.ConstructParticle();
255 }
256 
257 
258 #include "G4ProcessManager.hh"
259 void G4HadronPhysicsFTFQGSP_BERT::ConstructProcess()
260 {
261  CreateModels();
262 
263  G4ProcessManager * aProcMan = 0;
264 
265  theNeutronInelastic = new G4NeutronInelasticProcess();
266  theNeutronInelastic->RegisterMe( theModel1 );
267  theNeutronInelastic->RegisterMe( theBertini1 );
268  theNeutronInelastic->AddDataSet( new G4BGGNucleonInelasticXS( G4Neutron::Neutron() ) );
269  theNeutronInelastic->AddDataSet( theNeutronInelasticXS );
270  aProcMan = G4Neutron::Neutron()->GetProcessManager();
271  aProcMan->AddDiscreteProcess( theNeutronInelastic );
272  theNeutronCaptureProcess = new G4HadronCaptureProcess();
273  theNeutronCaptureProcess->RegisterMe( theNeutronCaptureModel );
274  aProcMan->AddDiscreteProcess( theNeutronCaptureProcess );
275 
276  theProtonInelastic = new G4ProtonInelasticProcess();
277  theProtonInelastic->RegisterMe( theModel1 );
278  theProtonInelastic->RegisterMe( theBertini1 );
279  theProtonInelastic->AddDataSet( new G4BGGNucleonInelasticXS( G4Proton::Proton() ) );
280  aProcMan = G4Proton::Proton()->GetProcessManager();
281  aProcMan->AddDiscreteProcess( theProtonInelastic );
282 
283  thePionMinusInelastic = new G4PionMinusInelasticProcess();
284  thePionMinusInelastic->RegisterMe( theModel1 );
285  thePionMinusInelastic->RegisterMe( theBertini1 );
286  thePionMinusInelastic->AddDataSet( thePiXS );
287  aProcMan = G4PionMinus::PionMinus()->GetProcessManager();
288  aProcMan->AddDiscreteProcess( thePionMinusInelastic );
289 
290  thePionPlusInelastic = new G4PionPlusInelasticProcess();
291  thePionPlusInelastic->RegisterMe( theModel1 );
292  thePionPlusInelastic->RegisterMe( theBertini1 );
293  thePionPlusInelastic->AddDataSet( thePiXS );
294  aProcMan = G4PionPlus::PionPlus()->GetProcessManager();
295  aProcMan->AddDiscreteProcess( thePionPlusInelastic );
296 
297  theKaonMinusInelastic = new G4KaonMinusInelasticProcess();
298  theKaonMinusInelastic->RegisterMe( theModel1 );
299  theKaonMinusInelastic->RegisterMe( theBertini1 );
300  theKaonMinusInelastic->AddDataSet( theKaonAndHyperonXS );
301  aProcMan = G4KaonMinus::KaonMinus()->GetProcessManager();
302  aProcMan->AddDiscreteProcess( theKaonMinusInelastic );
303 
304  theKaonPlusInelastic = new G4KaonPlusInelasticProcess();
305  theKaonPlusInelastic->RegisterMe( theModel1 );
306  theKaonPlusInelastic->RegisterMe( theBertini1 );
307  theKaonPlusInelastic->AddDataSet( theKaonAndHyperonXS );
308  aProcMan = G4KaonPlus::KaonPlus()->GetProcessManager();
309  aProcMan->AddDiscreteProcess( theKaonPlusInelastic );
310 
311  theKaonZeroLInelastic = new G4KaonZeroLInelasticProcess();
312  theKaonZeroLInelastic->RegisterMe( theModel1 );
313  theKaonZeroLInelastic->RegisterMe( theBertini1 );
314  theKaonZeroLInelastic->AddDataSet( theKaonAndHyperonXS );
315  aProcMan = G4KaonZeroLong::KaonZeroLong()->GetProcessManager();
316  aProcMan->AddDiscreteProcess( theKaonZeroLInelastic );
317 
318  theKaonZeroSInelastic = new G4KaonZeroSInelasticProcess();
319  theKaonZeroSInelastic->RegisterMe( theModel1 );
320  theKaonZeroSInelastic->RegisterMe( theBertini1 );
321  theKaonZeroSInelastic->AddDataSet( theKaonAndHyperonXS );
322  aProcMan = G4KaonZeroShort::KaonZeroShort()->GetProcessManager();
323  aProcMan->AddDiscreteProcess( theKaonZeroSInelastic );
324 
325  theLambdaInelastic = new G4LambdaInelasticProcess();
326  theLambdaInelastic->RegisterMe( theModel2 );
327  theLambdaInelastic->RegisterMe( theBertini2 );
328  theLambdaInelastic->AddDataSet( theKaonAndHyperonXS );
329  aProcMan = G4Lambda::Lambda()->GetProcessManager();
330  aProcMan->AddDiscreteProcess( theLambdaInelastic );
331 
332  theAntiLambdaInelastic = new G4AntiLambdaInelasticProcess();
333  theAntiLambdaInelastic->RegisterMe( theModel3 );
334  theAntiLambdaInelastic->AddDataSet( theKaonAndHyperonXS );
335  aProcMan = G4AntiLambda::AntiLambda()->GetProcessManager();
336  aProcMan->AddDiscreteProcess( theAntiLambdaInelastic );
337 
338  theSigmaMinusInelastic = new G4SigmaMinusInelasticProcess();
339  theSigmaMinusInelastic->RegisterMe( theModel2 );
340  theSigmaMinusInelastic->RegisterMe( theBertini2 );
341  theSigmaMinusInelastic->AddDataSet( theKaonAndHyperonXS );
342  aProcMan = G4SigmaMinus::SigmaMinus()->GetProcessManager();
343  aProcMan->AddDiscreteProcess( theSigmaMinusInelastic );
344 
345  theAntiSigmaMinusInelastic = new G4AntiSigmaMinusInelasticProcess();
346  theAntiSigmaMinusInelastic->RegisterMe( theModel3 );
347  theAntiSigmaMinusInelastic->AddDataSet( theKaonAndHyperonXS );
348  aProcMan = G4AntiSigmaMinus::AntiSigmaMinus()->GetProcessManager();
349  aProcMan->AddDiscreteProcess( theAntiSigmaMinusInelastic );
350 
351  theSigmaPlusInelastic = new G4SigmaPlusInelasticProcess();
352  theSigmaPlusInelastic->RegisterMe( theModel2 );
353  theSigmaPlusInelastic->RegisterMe( theBertini2 );
354  theSigmaPlusInelastic->AddDataSet( theKaonAndHyperonXS );
355  aProcMan = G4SigmaPlus::SigmaPlus()->GetProcessManager();
356  aProcMan->AddDiscreteProcess( theSigmaPlusInelastic );
357 
358  theAntiSigmaPlusInelastic = new G4AntiSigmaPlusInelasticProcess();
359  theAntiSigmaPlusInelastic->RegisterMe( theModel3 );
360  theAntiSigmaPlusInelastic->AddDataSet( theKaonAndHyperonXS );
361  aProcMan = G4AntiSigmaPlus::AntiSigmaPlus()->GetProcessManager();
362  aProcMan->AddDiscreteProcess( theAntiSigmaPlusInelastic );
363 
364  theXiMinusInelastic = new G4XiMinusInelasticProcess();
365  theXiMinusInelastic->RegisterMe( theModel2 );
366  theXiMinusInelastic->RegisterMe( theBertini2 );
367  theXiMinusInelastic->AddDataSet( theKaonAndHyperonXS );
368  aProcMan = G4XiMinus::XiMinus()->GetProcessManager();
369  aProcMan->AddDiscreteProcess( theXiMinusInelastic );
370 
371  theAntiXiMinusInelastic = new G4AntiXiMinusInelasticProcess();
372  theAntiXiMinusInelastic->RegisterMe( theModel3 );
373  theAntiXiMinusInelastic->AddDataSet( theKaonAndHyperonXS );
374  aProcMan = G4AntiXiMinus::AntiXiMinus()->GetProcessManager();
375  aProcMan->AddDiscreteProcess( theAntiXiMinusInelastic );
376 
377  theXiZeroInelastic = new G4XiZeroInelasticProcess();
378  theXiZeroInelastic->RegisterMe( theModel2 );
379  theXiZeroInelastic->RegisterMe( theBertini2 );
380  theXiZeroInelastic->AddDataSet( theKaonAndHyperonXS );
381  aProcMan = G4XiZero::XiZero()->GetProcessManager();
382  aProcMan->AddDiscreteProcess( theXiZeroInelastic );
383 
384  theAntiXiZeroInelastic = new G4AntiXiZeroInelasticProcess();
385  theAntiXiZeroInelastic->RegisterMe( theModel3 );
386  theAntiXiZeroInelastic->AddDataSet( theKaonAndHyperonXS );
387  aProcMan = G4AntiXiZero::AntiXiZero()->GetProcessManager();
388  aProcMan->AddDiscreteProcess( theAntiXiZeroInelastic );
389 
390  theOmegaMinusInelastic = new G4OmegaMinusInelasticProcess();
391  theOmegaMinusInelastic->RegisterMe( theModel2 );
392  theOmegaMinusInelastic->RegisterMe( theBertini2 );
393  theOmegaMinusInelastic->AddDataSet( theKaonAndHyperonXS );
394  aProcMan = G4OmegaMinus::OmegaMinus()->GetProcessManager();
395  aProcMan->AddDiscreteProcess( theOmegaMinusInelastic );
396 
397  theAntiOmegaMinusInelastic = new G4AntiOmegaMinusInelasticProcess();
398  theAntiOmegaMinusInelastic->RegisterMe( theModel3 );
399  theAntiOmegaMinusInelastic->AddDataSet( theKaonAndHyperonXS );
400  aProcMan = G4AntiOmegaMinus::AntiOmegaMinus()->GetProcessManager();
401  aProcMan->AddDiscreteProcess( theAntiOmegaMinusInelastic );
402 
403  theAntiProtonInelastic = new G4AntiProtonInelasticProcess();
404  theAntiProtonInelastic->RegisterMe( theModel3 );
405  theAntiProtonInelastic->AddDataSet( theAntiNucleonXS );
406  aProcMan = G4AntiProton::AntiProton()->GetProcessManager();
407  aProcMan->AddDiscreteProcess( theAntiProtonInelastic );
408 
409  theAntiNeutronInelastic = new G4AntiNeutronInelasticProcess();
410  theAntiNeutronInelastic->RegisterMe( theModel3 );
411  theAntiNeutronInelastic->AddDataSet( theAntiNucleonXS );
412  aProcMan = G4AntiNeutron::AntiNeutron()->GetProcessManager();
413  aProcMan->AddDiscreteProcess( theAntiNeutronInelastic );
414 
415  theAntiDeuteronInelastic = new G4AntiDeuteronInelasticProcess();
416  theAntiDeuteronInelastic->RegisterMe( theModel3 );
417  theAntiDeuteronInelastic->AddDataSet( theAntiNucleonXS );
418  aProcMan = G4AntiDeuteron::AntiDeuteron()->GetProcessManager();
419  aProcMan->AddDiscreteProcess( theAntiDeuteronInelastic );
420 
421  theAntiTritonInelastic = new G4AntiTritonInelasticProcess();
422  theAntiTritonInelastic->RegisterMe( theModel3 );
423  theAntiTritonInelastic->AddDataSet( theAntiNucleonXS );
424  aProcMan = G4AntiTriton::AntiTriton()->GetProcessManager();
425  aProcMan->AddDiscreteProcess( theAntiTritonInelastic );
426 
427  theAntiHe3Inelastic = new G4AntiHe3InelasticProcess();
428  theAntiHe3Inelastic->RegisterMe( theModel3 );
429  theAntiHe3Inelastic->AddDataSet( theAntiNucleonXS );
430  aProcMan = G4AntiHe3::AntiHe3()->GetProcessManager();
431  aProcMan->AddDiscreteProcess( theAntiHe3Inelastic );
432 
433  theAntiAlphaInelastic = new G4AntiAlphaInelasticProcess();
434  theAntiAlphaInelastic->RegisterMe( theModel3 );
435  theAntiAlphaInelastic->AddDataSet( theAntiNucleonXS );
436  aProcMan = G4AntiAlpha::AntiAlpha()->GetProcessManager();
437  aProcMan->AddDiscreteProcess( theAntiAlphaInelastic );
438 }
439