| 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 // INCL++ intra-nuclear cascade model
27 // Alain Boudard, CEA-Saclay, France
28 // Joseph Cugnon, University of Liege, Belgium
29 // Jean-Christophe David, CEA-Saclay, France
30 // Pekka Kaitaniemi, CEA-Saclay, France, and H
31 // Sylvie Leray, CEA-Saclay, France
32 // Davide Mancusi, CEA-Saclay, France
33 //
34 #define INCLXX_IN_GEANT4_MODE 1
35
36 #include "globals.hh"
37
38 #include "G4INCLNNToMissingStrangenessChannel.
39 #include "G4INCLKinematicsUtils.hh"
40 #include "G4INCLBinaryCollisionAvatar.hh"
41 #include "G4INCLRandom.hh"
42 #include "G4INCLGlobals.hh"
43 #include "G4INCLLogger.hh"
44 #include <algorithm>
45 #include "G4INCLPhaseSpaceGenerator.hh"
46
47 namespace G4INCL {
48
49 const G4double NNToMissingStrangenessChannel
50
51 NNToMissingStrangenessChannel::NNToMissingSt
52 : particle1(p1), particle2(p2)
53 {}
54
55 NNToMissingStrangenessChannel::~NNToMissingS
56
57 void NNToMissingStrangenessChannel::fillFina
58
59 const G4double sqrtS = KinematicsUtils::to
60 const G4double pLab = 0.001*KinematicsUtil
61 // assert(sqrtS > 3100); // ! > 3047.34. Not s
62
63 const G4int iso = ParticleTable::getIsospi
64 // assert(iso == -2 || iso == 0 || iso == 2);
65 G4int iso_system = 0;
66 G4int available_iso = 0;
67 G4int nbr_pions = 0;
68 G4int min_pions = 0;
69 G4int max_pions = 0;
70
71 G4double rdm = Random::shoot();
72
73 //G4int nbr_particle = 3;
74
75 if(rdm < 0.35){
76 // N-K-Lambda chosen
77 particle2->setType(Lambda);
78 available_iso = 2;
79 min_pions = 3;
80 max_pions = G4int((sqrtS-ParticleTable::
81 }
82 else if((iso == 0 && rdm < 0.55) || rdm <
83 // N-N-K-Kb chosen
84 //nbr_particle++;
85 available_iso = 4;
86 min_pions = 1;
87 max_pions = G4int((sqrtS-2.*ParticleTabl
88 }
89 else{
90 // N-K-Sigma chosen
91 available_iso = 4;
92 min_pions = 3;
93 max_pions = G4int((sqrtS-ParticleTable::
94 }
95 /* Gaussian noise + mean
96 G4double intermediaire = min_pions + Rando
97 nbr_pions = std::min(max_pions,std::max(mi
98
99 available_iso += nbr_pions*2;
100 //nbr_particle += nbr_pions;
101
102 ParticleList list;
103 ParticleType PionType = PiZero;
104 const ThreeVector &rcol1 = particle1->getP
105 const ThreeVector zero;
106 Particle *kaon = new Particle(KPlus,zero,r
107
108 for(Int_t i=0; i<nbr_pions; i++){
109 Particle *pion = new Particle(PionType,z
110 if(available_iso-std::abs(iso-iso_system
111 rdm = Random::shoot();
112 if(((iso == 0) && rdm < 0.3) || ((iso
113 pion->setType(PiMinus);
114 iso_system -= 2;
115 available_iso -= 2;
116 }
117 else if(((iso == 0) && rdm < 0.7) || (
118 pion->setType(PiZero);
119 available_iso -= 2;
120 }
121 else{
122 pion->setType(PiPlus);
123 iso_system += 2;
124 available_iso -= 2;
125 }
126 }
127 else if(available_iso-std::abs(iso-iso_s
128 rdm = Random::shoot();
129 // pn pp too hig
130 if(((iso == 0) && (rdm*0.7 < 0.3)) ||
131 pion->setType(ParticleTable::getPion
132 iso_system += Math::sign(iso-iso_sys
133 available_iso -= 2;
134 }
135 else{
136 PionType = PiZero;
137 available_iso -= 2;
138 }
139 }
140 else if(available_iso-std::abs(iso-iso_s
141 pion->setType(ParticleTable::getPionTy
142 iso_system += Math::sign(iso-iso_syste
143 available_iso -= 2;
144 }
145 else INCL_ERROR("Pion Generation Problem
146 list.push_back(pion);
147 }
148
149 if(particle2->isLambda()){ // N-K-Lambda
150 // assert(available_iso == 2);
151 if(std::abs(iso-iso_system) == 2){
152 particle1->setType(ParticleTable::getN
153 kaon->setType(ParticleTable::getKaonTy
154 }
155 else if(std::abs(iso-iso_system) == 0){
156 rdm = G4int(Random::shoot()*2.)*2-1;
157 particle1->setType(ParticleTable::getN
158 kaon->setType(ParticleTable::getKaonTy
159 }
160 else INCL_ERROR("Isospin non-conservatio
161 }
162 else if(min_pions == 1){ // N-N-K-Kb chose
163 // assert(available_iso == 4);
164 Particle *antikaon = new Particle(KMinus
165 if(std::abs(iso-iso_system) == 4){
166 particle1->setType(ParticleTable::getN
167 particle2->setType(ParticleTable::getN
168 kaon->setType(ParticleTable::getKaonTy
169 antikaon->setType(ParticleTable::getAn
170 }
171 else if(std::abs(iso-iso_system) == 2){
172 rdm = G4int(Random::shoot()*2.)*2-1;
173 particle1->setType(ParticleTable::getN
174 particle2->setType(ParticleTable::getN
175 kaon->setType(ParticleTable::getKaonTy
176 antikaon->setType(ParticleTable::getAn
177 }
178 else if(std::abs(iso-iso_system) == 0){
179 rdm = G4int(Random::shoot()*2.)*2-1;
180 G4double rdm2 = G4int(Random::shoot()*
181 particle1->setType(ParticleTable::getN
182 particle2->setType(ParticleTable::getN
183 kaon->setType(ParticleTable::getKaonTy
184 antikaon->setType(ParticleTable::getAn
185 }
186 else INCL_ERROR("Isospin non-conservatio
187 list.push_back(antikaon);
188 nbr_pions += 1; // need for addCreatedPa
189 }
190 else{// N-K-Sigma
191 // assert(available_iso == 4);
192 if(std::abs(iso-iso_system) == 4){
193 particle1->setType(ParticleTable::getN
194 particle2->setType(ParticleTable::getS
195 kaon->setType(ParticleTable::getKaonTy
196 }
197 else if(std::abs(iso-iso_system) == 2){
198 rdm = Random::shoot();
199 // pp(pn) sp:0.42(0.23) sz:0.51(0.54)
200 if(((iso == 0) && (rdm*0.77 < 0.23)) |
201 particle2->setType(ParticleTable::ge
202 rdm = G4int(Random::shoot()*2.)*2-1;
203 particle1->setType(ParticleTable::ge
204 kaon->setType(ParticleTable::getKaon
205 }
206 else{
207 particle2->setType(SigmaZero);
208 particle1->setType(ParticleTable::ge
209 kaon->setType(ParticleTable::getKaon
210 }
211 }
212 else if(std::abs(iso-iso_system) == 0){
213 if(((iso == 0) && rdm < 0.23) || ((iso
214 particle2->setType(SigmaPlus);
215 particle1->setType(Neutron);
216 kaon->setType(KZero);
217 }
218 else if(((iso == 0) && rdm < 0.77) ||
219 particle2->setType(SigmaZero);
220 rdm = G4int(Random::shoot()*2.)*2-1;
221 particle1->setType(ParticleTable::ge
222 kaon->setType(ParticleTable::getKaon
223 }
224 else{
225 particle2->setType(SigmaMinus);
226 particle1->setType(Proton);
227 kaon->setType(KPlus);
228 }
229 }
230 else INCL_ERROR("Isospin non-conservatio
231 }
232
233 list.push_back(particle1);
234 list.push_back(particle2);
235 list.push_back(kaon);
236
237 PhaseSpaceGenerator::generateBiased(sqrtS,
238
239 fs->addModifiedParticle(particle1);
240 fs->addModifiedParticle(particle2);
241 fs->addCreatedParticle(kaon);
242 for(Int_t i=0; i<nbr_pions; i++) fs->addCr
243
244 }
245 }
246