| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 /// \file hadronic/Hadr02/src/HadronPhysicsHIJ
27 /// \brief Implementation of the HadronPhysics
28 //
29 //
30 //--------------------------------------------
31 //
32 // ClassName:
33 //
34 // Author: 2012 A. Dotti
35 // created from HadronPhysicsHIJING
36 //
37 // Modified:
38 //
39 //--------------------------------------------
40 //
41 #ifdef G4_USE_HIJING
42 # include "HadronPhysicsHIJING.hh"
43
44 # include "G4AntiBarionBuilder.hh"
45 # include "G4BaryonConstructor.hh"
46 # include "G4BertiniNeutronBuilder.hh"
47 # include "G4BertiniPiKBuilder.hh"
48 # include "G4BertiniProtonBuilder.hh"
49 # include "G4ChipsKaonMinusInelasticXS.hh"
50 # include "G4ChipsKaonPlusInelasticXS.hh"
51 # include "G4ChipsKaonZeroInelasticXS.hh"
52 # include "G4CrossSectionDataSetRegistry.hh"
53 # include "G4FTFPAntiBarionBuilder.hh"
54 # include "G4FTFPNeutronBuilder.hh"
55 # include "G4FTFPPiKBuilder.hh"
56 # include "G4FTFPProtonBuilder.hh"
57 # include "G4HyperonFTFPBuilder.hh"
58 # include "G4MesonConstructor.hh"
59 # include "G4NeutronBuilder.hh"
60 # include "G4ParticleDefinition.hh"
61 # include "G4ParticleTable.hh"
62 # include "G4PhysListUtil.hh"
63 # include "G4PiKBuilder.hh"
64 # include "G4ProcessManager.hh"
65 # include "G4ProtonBuilder.hh"
66 # include "G4QGSPNeutronBuilder.hh"
67 # include "G4QGSPPiKBuilder.hh"
68 # include "G4QGSPProtonBuilder.hh"
69 # include "G4ShortLivedConstructor.hh"
70 # include "G4ios.hh"
71 # include "globals.hh"
72
73 # include <iomanip>
74
75 HadronPhysicsHIJING::HadronPhysicsHIJING(G4int
76 {
77 fNeutrons = 0;
78 fHIJINGNeutron = 0;
79 fPiK = 0;
80 fHIJINGPiK = 0;
81 fPro = 0;
82 fHIJINGPro = 0;
83 fHyperon = 0;
84 fFTFPHyperon = 0;
85 fAntiBaryon = 0;
86 fHIJINGAntiBaryon = 0;
87 fCHIPSInelastic = 0;
88 }
89
90 //....oooOO0OOooo........oooOO0OOooo........oo
91
92 void HadronPhysicsHIJING::CreateModels()
93 {
94 G4bool quasiElasFTF = false;
95 G4bool quasiElasQGS = true;
96 G4double minFTFP = 6 * GeV;
97 G4double maxBERT = 8 * GeV;
98 G4double minHIJING = 12 * GeV;
99 G4double minQGSP = 12 * GeV;
100 G4double maxFTFP = 25 * GeV;
101
102 // Proton
103 fPro = new G4ProtonBuilder;
104 fHIJINGPro = new HIJINGProtonBuilder();
105 fHIJINGPro->SetMinEnergy(minHIJING);
106 fPro->RegisterMe(fHIJINGPro);
107 G4FTFPProtonBuilder* FTFPPro = new G4FTFPPro
108 FTFPPro->SetMinEnergy(minFTFP);
109 FTFPPro->SetMaxEnergy(maxFTFP);
110 fPro->RegisterMe(FTFPPro);
111 G4BertiniProtonBuilder* BertPro = new G4Bert
112 BertPro->SetMaxEnergy(maxBERT);
113
114 fNeutrons = new G4NeutronBuilder;
115 fHIJINGNeutron = new HIJINGNeutronBuilder();
116 fHIJINGNeutron->SetMinEnergy(minHIJING);
117 fNeutrons->RegisterMe(fHIJINGNeutron);
118 // G4QGSPNeutronBuilder* QGSPNeu = new G4QGS
119 // QGSPNeu->SetMinEnergy(minQGSP);
120 // fNeutrons->RegisterMe(QGSPNeu);
121 G4FTFPNeutronBuilder* FTFPNeu = new G4FTFPNe
122 FTFPNeu->SetMinEnergy(minFTFP);
123 FTFPNeu->SetMaxEnergy(maxFTFP);
124 fNeutrons->RegisterMe(FTFPNeu);
125 G4BertiniNeutronBuilder* BertNeu = new G4Ber
126 BertNeu->SetMaxEnergy(maxBERT);
127
128 fPiK = new G4PiKBuilder;
129 G4QGSPPiKBuilder* QGSPPiK = new G4QGSPPiKBui
130 fPiK->RegisterMe(QGSPPiK);
131 QGSPPiK->SetMinEnergy(minQGSP);
132 G4FTFPPiKBuilder* FTFPPiK = new G4FTFPPiKBui
133 fPiK->RegisterMe(FTFPPiK);
134 FTFPPiK->SetMaxEnergy(maxFTFP);
135 FTFPPiK->SetMinEnergy(minFTFP);
136 G4BertiniPiKBuilder* BertiniPiK = new G4Bert
137 fPiK->RegisterMe(BertiniPiK);
138 BertiniPiK->SetMaxEnergy(maxBERT);
139
140 // For Hyperons use FTF model
141 fHyperon = new G4HyperonBuilder;
142 fFTFPHyperon = new G4HyperonFTFPBuilder;
143 fHyperon->RegisterMe(fFTFPHyperon);
144
145 fAntiBaryon = new G4AntiBarionBuilder;
146 // FTFPAntiBaryon=new G4FTFPAntiBarionBuilde
147 G4FTFPAntiBarionBuilder* FTFPAB = new G4FTFP
148 fAntiBaryon->RegisterMe(FTFPAB);
149 }
150
151 //....oooOO0OOooo........oooOO0OOooo........oo
152
153 HadronPhysicsHIJING::~HadronPhysicsHIJING()
154 {
155 delete fNeutrons;
156 delete fHIJINGNeutron;
157
158 delete fPiK;
159 delete fHIJINGPiK;
160
161 delete fPro;
162 delete fHIJINGPro;
163
164 delete fHyperon;
165 delete fFTFPHyperon;
166 delete fAntiBaryon;
167 delete fHIJINGAntiBaryon;
168
169 delete fCHIPSInelastic;
170 }
171
172 //....oooOO0OOooo........oooOO0OOooo........oo
173
174 void HadronPhysicsHIJING::ConstructParticle()
175 {
176 G4MesonConstructor pMesonConstructor;
177 pMesonConstructor.ConstructParticle();
178
179 G4BaryonConstructor pBaryonConstructor;
180 pBaryonConstructor.ConstructParticle();
181
182 G4ShortLivedConstructor pShortLivedConstruct
183 pShortLivedConstructor.ConstructParticle();
184 }
185
186 //....oooOO0OOooo........oooOO0OOooo........oo
187
188 void HadronPhysicsHIJING::ConstructProcess()
189 {
190 CreateModels();
191 fNeutrons->Build();
192 fPro->Build();
193 fPiK->Build();
194
195 // use CHIPS cross sections also for Kaons
196 G4VCrossSectionDataSet* ChipsKaonMinus =
197 G4CrossSectionDataSetRegistry::Instance()-
198 G4ChipsKaonMinusInelasticXS::Default_Nam
199 G4VCrossSectionDataSet* ChipsKaonPlus =
200 G4CrossSectionDataSetRegistry::Instance()-
201 G4ChipsKaonPlusInelasticXS::Default_Name
202 G4VCrossSectionDataSet* ChipsKaonZero =
203 G4CrossSectionDataSetRegistry::Instance()-
204 G4ChipsKaonZeroInelasticXS::Default_Name
205 //
206
207 G4PhysListUtil::FindInelasticProcess(G4KaonM
208 G4PhysListUtil::FindInelasticProcess(G4KaonP
209 G4PhysListUtil::FindInelasticProcess(G4KaonZ
210 G4PhysListUtil::FindInelasticProcess(G4KaonZ
211
212 fHyperon->Build();
213 fAntiBaryon->Build();
214 }
215
216 //....oooOO0OOooo........oooOO0OOooo........oo
217
218 G4HadronicProcess* HadronPhysicsHIJING::FindIn
219 {
220 G4HadronicProcess* had = 0;
221 if (p) {
222 G4ProcessVector* pvec = p->GetProcessManag
223 size_t n = pvec->size();
224 if (0 < n) {
225 for (size_t i = 0; i < n; ++i) {
226 if (fHadronInelastic == ((*pvec)[i])->
227 had = static_cast<G4HadronicProcess*
228
229 break;
230 }
231 }
232 }
233 }
234 return had;
235 }
236
237 #endif // G4_USE_HIJING
238