Geant4  v4-10.4-release
 모두 클래스 네임스페이스들 파일들 함수 변수 타입정의 열거형 타입 열거형 멤버 Friends 매크로 그룹들 페이지들
G4INCLNNToNSKpiChannel.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 // 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 Helsinki Institute of Physics, Finland
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 
39 #include "G4INCLKinematicsUtils.hh"
41 #include "G4INCLRandom.hh"
42 #include "G4INCLGlobals.hh"
43 #include "G4INCLLogger.hh"
44 #include <algorithm>
46 
47 namespace G4INCL {
48 
49  const G4double NNToNSKpiChannel::angularSlope = 2.; // What is the exact effect? Sould be check
50 
52  : particle1(p1), particle2(p2)
53  {}
54 
56 
58 
59  // pp (36) pn (36)
60  //
61  // pp -> p pi+ S- K+ (9)
62  // pp -> p pi+ S0 K0 (9)
63  // pp -> p pi0 S+ K0 (4)
64  // pp -> n pi+ S+ K0 (2)
65  // pp -> p pi0 S0 K+ (4)
66  // pp -> n pi+ S0 K+ (2)
67  // pp -> p pi- S+ K+ (2)
68  // pp -> n pi0 S+ K+ (4)
69 
70  // pn -> p pi0 S- K+ (4)
71  // pn -> n pi+ S- K+ (2)
72  // pn -> p pi0 S0 K0 (2)
73  // pn -> n pi+ S0 K0 (1)
74  // pn -> p pi+ S- K0 (9)
75 
77 
79 
80  ParticleType KaonType;
81  ParticleType PionType;
82 
83  G4double rdm = Random::shoot();
84 
85  if(iso == 2){
86  if(rdm * 36. < 9.){
87  KaonType = KPlus;
88  PionType = PiPlus;
90  }
91  else if(rdm * 36. < 18.){
92  KaonType = KZero;
93  PionType = PiPlus;
95  }
96  else if(rdm * 36. < 22.){
97  KaonType = KZero;
98  PionType = PiZero;
100  }
101  else if(rdm * 36. < 24.){
102  KaonType = KZero;
103  PionType = PiPlus;
106  }
107  else if(rdm * 36. < 28.){
108  KaonType = KPlus;
109  PionType = PiZero;
111  }
112  else if(rdm * 36. < 30.){
113  KaonType = KPlus;
114  PionType = PiPlus;
117  }
118  else if(rdm * 36. < 32.){
119  KaonType = KPlus;
120  PionType = PiMinus;
122  }
123  else{
124  KaonType = KPlus;
125  PionType = PiZero;
128  }
129 
130  }
131  else if(iso == -2){
132  if(rdm * 36. < 9.){
133  KaonType = KZero;
134  PionType = PiMinus;
136  }
137  else if(rdm * 36. < 18.){
138  KaonType = KPlus;
139  PionType = PiMinus;
141  }
142  else if(rdm * 36. < 22.){
143  KaonType = KPlus;
144  PionType = PiZero;
146  }
147  else if(rdm * 36. < 24.){
148  KaonType = KPlus;
149  PionType = PiMinus;
152  }
153  else if(rdm * 36. < 28.){
154  KaonType = KZero;
155  PionType = PiZero;
157  }
158  else if(rdm * 36. < 30.){
159  KaonType = KZero;
160  PionType = PiMinus;
163  }
164  else if(rdm * 36. < 32.){
165  KaonType = KZero;
166  PionType = PiPlus;
168  }
169  else{
170  KaonType = KZero;
171  PionType = PiZero;
174  }
175 
176  }
177  else if(rdm*36. < 4.){
178  KaonType = KPlus;
179  PionType = PiZero;
182  }
183  else if(rdm*36. < 6.){
184  KaonType = KZero;
185  PionType = PiZero;
188  }
189  else if(rdm*36. < 8.){
190  KaonType = KPlus;
191  PionType = PiPlus;
194  }
195  else if(rdm*36. < 9.){
196  KaonType = KZero;
197  PionType = PiMinus;
200  }
201  else if(rdm*36. < 18.){
202  KaonType = KZero;
203  PionType = PiZero;
206  }
207  else if(rdm*36. < 27.){
208  KaonType = KPlus;
209  PionType = PiZero;
212  }
213  else if(rdm*36. < 28.){
214  KaonType = KZero;
215  PionType = PiPlus;
218  }
219  else if(rdm*36. < 30.){
220  KaonType = KPlus;
221  PionType = PiMinus;
224  }
225  else if(rdm*36. < 32.){
226  KaonType = KZero;
227  PionType = PiPlus;
230  }
231  else{
232  KaonType = KPlus;
233  PionType = PiMinus;
236  }
237 
238  ParticleList list;
239  list.push_back(particle1);
240  list.push_back(particle2);
241  const ThreeVector &rcol1 = particle1->getPosition();
242  const ThreeVector &rcol2 = particle2->getPosition();
243  const ThreeVector zero;
244  Particle *pion = new Particle(PionType,zero,rcol1);
245  Particle *kaon = new Particle(KaonType,zero,rcol2);
246  list.push_back(kaon);
247  list.push_back(pion);
248 
250 
251  INCL_DEBUG("NNToNSKpi " << (kaon->getMomentum().theta()) * 180. / G4INCL::Math::pi << '\n');
252 
255  fs->addCreatedParticle(kaon);
256  fs->addCreatedParticle(pion);
257 
258  }
259 }
void generateBiased(const G4double sqrtS, ParticleList &particles, const size_t index, const G4double slope)
Generate a biased event in the CM system.
G4bool pion(G4int ityp)
void addModifiedParticle(Particle *p)
const G4INCL::ThreeVector & getPosition() const
G4double totalEnergyInCM(Particle const *const p1, Particle const *const p2)
const G4INCL::ThreeVector & getMomentum() const
#define INCL_DEBUG(x)
static const G4double angularSlope
double G4double
Definition: G4Types.hh:76
G4double shoot()
Definition: G4INCLRandom.cc:93
void fillFinalState(FinalState *fs)
const G4double pi
void setType(ParticleType t)
NNToNSKpiChannel(Particle *, Particle *)
int G4int
Definition: G4Types.hh:78
G4int getIsospin(const ParticleType t)
Get the isospin of a particle.
G4INCL::ParticleType getType() const
G4double theta() const
void addCreatedParticle(Particle *p)