| 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 //
27
28
29 // -------------------------------------------
30 // GEANT 4 class implementation file
31 //
32 // ---------------- G4FragmentingString -
33 // by Gunter Folger, September 200
34 // class for an excited string used in F
35 // -------------------------------------------
36
37
38 // G4FragmentingString
39 #include "G4FragmentingString.hh"
40 #include "G4ExcitedString.hh"
41
42 //--------------------------------------------
43
44 //--------------------------------------------
45
46 G4FragmentingString::G4FragmentingString(const
47 {
48 LeftParton=old.LeftParton;
49 RightParton=old.RightParton;
50 Ptleft=old.Ptleft;
51 Ptright=old.Ptright;
52 Pplus=old.Pplus;
53 Pminus=old.Pminus;
54 theStableParton=old.theStableParton;
55 theDecayParton=old.theDecayParton;
56 decaying=old.decaying;
57 Pstring=old.Pstring;
58 Pleft =old.Pleft;
59 Pright =old.Pright;
60 }
61
62 G4FragmentingString & G4FragmentingString::ope
63 {
64 if (this != &old)
65 {
66 LeftParton=old.LeftParton;
67 RightParton=old.RightParton;
68 Ptleft=old.Ptleft;
69 Ptright=old.Ptright;
70 Pplus=old.Pplus;
71 Pminus=old.Pminus;
72 theStableParton=old.theStableParton;
73 theDecayParton=old.theDecayParton;
74 decaying=old.decaying;
75 Pstring=old.Pstring;
76 Pleft =old.Pleft;
77 Pright =old.Pright;
78 }
79 return *this;
80 }
81
82 //--------------------------------------------
83
84 G4FragmentingString::G4FragmentingString(const
85 {
86 LeftParton=excited.GetLeftParton()->GetDefin
87 RightParton=excited.GetRightParton()->GetDef
88 Ptleft=excited.GetLeftParton()->Get4Momentum
89 Ptleft.setZ(0.);
90 Ptright=excited.GetRightParton()->Get4Moment
91 Ptright.setZ(0.);
92 theStableParton=0;
93 theDecayParton=0;
94
95 if (excited.GetDirection() > 0) {decay
96 else {decay
97
98 Pleft = excited.GetLeftParton()->Get4
99 Pright = excited.GetRightParton()->Get
100 Pstring= Pleft + Pright;
101
102 Pplus = Pstring.plus();
103 Pminus = Pstring.minus();
104 }
105
106 //--------------------------------------------
107
108 G4FragmentingString::G4FragmentingString(const
109 G4ParticleDefinition * newdecay,
110 const G4LorentzVector *momentum)
111 {
112 decaying=None;
113 // Momentum of produced hadron
114 G4LorentzVector Momentum = G4LorentzVe
115
116 if ( old.decaying == Left )
117 {
118 RightParton= old.RightParton;
119 Ptright = old.Ptright;
120 Pright = old.Pright;
121
122 LeftParton = newdecay;
123 Ptleft = old.Ptleft - mome
124 Ptleft.setZ(0.);
125 Pleft = old.Pleft - Momentum;
126
127 Pstring = Pleft + Pright;
128 Pplus = Pstring.plus();
129 Pminus = Pstring.minus();
130
131 theDecayParton=GetLeftParton()
132 theStableParton=GetRightParton
133 decaying = Left;
134 } else if ( old.decaying == Right )
135 {
136 RightParton = newdecay;
137 Ptright = old.Ptright - mo
138 Ptright.setZ(0.);
139 Pright = old.Pright - Momentum
140
141 LeftParton = old.LeftParton;
142 Ptleft = old.Ptleft;
143 Pleft = old.Pleft;
144
145 Pstring = Pleft + Pright;
146 Pplus = Pstring.plus();
147 Pminus = Pstring.minus();
148
149 theDecayParton=GetRightParton(
150 theStableParton=GetLeftParton(
151 decaying = Right;
152 } else
153 {
154 throw G4HadronicException(__FILE__, __LINE
155 "G4FragmentingString::G4Frag
156 }
157 }
158
159
160 //--------------------------------------------
161
162 G4FragmentingString::G4FragmentingString(const
163 G4ParticleDefinition * newdecay)
164 {
165 decaying=None;
166
167 Ptleft.setX(0.); Ptleft.setY(0.); Pt
168 Ptright.setX(0.); Ptright.setY(0.); Pt
169 Pplus=0.; Pminus=0.;
170 theStableParton=0; theDecayParton=0;
171
172 Pstring=G4LorentzVector(0.,0.,0.,0.);
173 Pleft =G4LorentzVector(0.,0.,0.,0.);
174 Pright =G4LorentzVector(0.,0.,0.,0.);
175
176 if ( old.decaying == Left )
177 {
178 RightParton= old.RightParton;
179 LeftParton = newdecay;
180 decaying=Left;
181 } else if ( old.decaying == Right )
182 {
183 RightParton = newdecay;
184 LeftParton = old.LeftParton;
185 decaying=Right;
186 } else
187 {
188 throw G4HadronicException(__FILE__, __LINE
189 "G4FragmentingString::G4Frag
190 }
191 }
192
193
194 //--------------------------------------------
195
196 G4FragmentingString::~G4FragmentingString()
197 {}
198
199 //--------------------------------------------
200
201 void G4FragmentingString::SetLeftPartonStable(
202 {
203 theStableParton=GetLeftParton();
204 theDecayParton=GetRightParton();
205 decaying=Right;
206 }
207
208 //--------------------------------------------
209
210 void G4FragmentingString::SetRightPartonStable
211 {
212 theStableParton=GetRightParton();
213 theDecayParton=GetLeftParton();
214 decaying=Left;
215 }
216
217 //--------------------------------------------
218
219 G4int G4FragmentingString::GetDecayDirection()
220 {
221 if (decaying == Left ) return +1;
222 else if (decaying == Right) return -1;
223 else throw G4HadronicException(__FILE__, __L
224 "G4FragmentingString::GetDecayD
225 return 0;
226 }
227
228 //--------------------------------------------
229
230 G4bool G4FragmentingString::IsAFourQuarkString
231 {
232 return LeftParton->GetParticleSubType()==
233 && RightParton->GetParticleSubType()==
234 }
235
236 //--------------------------------------------
237
238 G4bool G4FragmentingString::DecayIsQuark()
239 {
240 return theDecayParton->GetParticleSubType()=
241 }
242
243 G4bool G4FragmentingString::StableIsQuark()
244 {
245 return theStableParton->GetParticleSubType()
246 }
247
248 //--------------------------------------------
249
250 G4ThreeVector G4FragmentingString::StablePt()
251 {
252 if (decaying == Left ) return Ptright;
253 else if (decaying == Right ) return Ptleft;
254 else throw G4HadronicException(__FILE__, __L
255 return G4ThreeVector();
256 }
257
258 G4ThreeVector G4FragmentingString::DecayPt()
259 {
260 if (decaying == Left ) return Ptleft;
261 else if (decaying == Right ) return Ptright;
262 else throw G4HadronicException(__FILE__, __L
263 return G4ThreeVector();
264 }
265
266 //--------------------------------------------
267
268 G4double G4FragmentingString::LightConePlus()
269 {
270 return Pplus;
271 }
272
273 G4double G4FragmentingString::LightConeMinus()
274 {
275 return Pminus;
276 }
277
278 G4double G4FragmentingString::LightConeDecay()
279 {
280 if (decaying == Left ) return Pplus;
281 else if (decaying == Right ) return Pminus;
282 else throw G4HadronicException(__FILE__, __L
283 }
284
285 //--------------------------------------------
286
287 G4LorentzVector G4FragmentingString::Get4Momen
288 {
289 return Pstring;
290 }
291
292 G4double G4FragmentingString::Mass2() const
293 {
294 return Pstring.mag2();
295 }
296
297 G4double G4FragmentingString::Mass() const
298 {
299 return Pstring.mag();
300 }
301
302 G4double G4FragmentingString::MassT2() const
303 {
304 return Pplus*Pminus;
305 }
306
307 G4LorentzVector G4FragmentingString::GetPstrin
308 {
309 return Pstring;
310 }
311
312 G4LorentzVector G4FragmentingString::GetPleft(
313 {
314 return Pleft;
315 }
316
317 G4LorentzVector G4FragmentingString::GetPright
318 {
319 return Pright;
320 }
321
322 G4LorentzRotation G4FragmentingString::Transfo
323 {
324 G4LorentzVector momentum = Pstring;
325 G4LorentzRotation toAlignedCms(-1*momentu
326 momentum = toAlignedCms * Pleft;
327
328 toAlignedCms.rotateZ(-1*momentum.phi());
329 toAlignedCms.rotateY(-1*momentum.theta())
330
331 Pleft *= toAlignedCms;
332 Pright *= toAlignedCms;
333 Pstring *= toAlignedCms;
334
335 Ptleft = G4ThreeVector(Pleft.vect());
336 Ptleft.setZ(0.);
337 Ptright = G4ThreeVector(Pright.vect());
338 Pplus = Pstring.plus();
339 Pminus = Pstring.minus();
340
341 return toAlignedCms;
342 }
343