Geant4  v4-10.4-release
 모두 클래스 네임스페이스들 파일들 함수 변수 타입정의 열거형 타입 열거형 멤버 Friends 매크로 그룹들 페이지들
G4Reggeons.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 // $Id: G4Reggeons.cc 99348 2016-09-19 08:39:04Z vuzhinsk $
28 //
29 
30 #include "G4Reggeons.hh"
31 #include "G4PhysicalConstants.hh"
32 #include "G4SystemOfUnits.hh"
33 #include "G4Pow.hh"
34 #include "G4Exp.hh"
35 #include "G4Log.hh"
36 
37 
39 {
40 // KP Orig
41  Alpha_pomeron = 1.12; //0.9808;
42  Alphaprime_pomeron = 0.22/GeV/GeV; //0.25/GeV/GeV;
43  S0_pomeron = 1.0*GeV*GeV; //2.7*GeV*GeV;
44 
45  Alpha_pomeronHard = 1.47;
46  Gamma_pomeronHard = 0.0/GeV/GeV;
47 
48  G4int PDGcode = particle->GetPDGEncoding();
49  G4int absPDGcode = std::abs(PDGcode);
50 
51 //-------------------------------------------------------
52 // KP Orig
53  G4double C_pomeron_NN = 1.5;
54  G4double C_pomeron_N = std::sqrt(C_pomeron_NN);
55 
56  G4double Gamma_pomeron_NN = 2.14/GeV/GeV; //(2.6+3.96)
57  G4double Gamma_pomeron_N = std::sqrt(Gamma_pomeron_NN);
58  G4double Gamma_pomeron_Pr(0.), Gamma_pomeron_Tr(0.);
59 
60  G4double Rsquare_pomeron_NN = 3.30/GeV/GeV; //3.56
61  G4double Rsquare_pomeron_N = Rsquare_pomeron_NN/2.;
62  G4double Rsquare_pomeron_Pr(0.), Rsquare_pomeron_Tr(0.);
63 //-------------------------------------------------------
64 
65  if( absPDGcode > 1000 ) { // Projectile is baryon or anti_baryon --------
66  Cpr_pomeron = C_pomeron_N; // Shower enhancement coefficient for projectile
67  Ctr_pomeron = C_pomeron_N; // Shower enhancement coefficient for target
69 
70  Gamma_pomeron_Pr = Gamma_pomeron_N; // vertex constant for projectile
71  Gamma_pomeron_Tr = Gamma_pomeron_N; // vertex constant for target
72  Gamma_pomeron = Gamma_pomeron_Pr * Gamma_pomeron_Tr;
73 
74  Rsquare_pomeron_Pr = Rsquare_pomeron_N; // R^2 of pomeron-projectile interaction
75  Rsquare_pomeron_Tr = Rsquare_pomeron_N; // R^2 of pomeron-target interaction
76  Rsquare_pomeron = Rsquare_pomeron_Pr + Rsquare_pomeron_Tr;
77 
78  Freggeon_Alpha = 0.7; // Intersept of f-trajectory
79  Freggeon_Alphaprime = 0.8/GeV/GeV; // Slope of f-trajectory
80  Freggeon_Gamma = sqr(2.871)/GeV/GeV; // Vertex constant of f-meson - nucleon interactions
81  Freggeon_Rsquare = 2*0.916/GeV/GeV; // R^2 of f-meson - nucleon interactions
82  Freggeon_C = 1.0; // Shower enhancement coefficient
83  FParity = +1; // Parity of the trajectory
84 
85  Wreggeon_Alpha = 0.4; // Intersept of omega-trajectory (w)
86  Wreggeon_Alphaprime = 0.9/GeV/GeV; // Slope of w-trajectory
87  Wreggeon_Gamma = sqr(2.241)/GeV/GeV; // Vertex constant of w-meson - nucleon interactions
88  Wreggeon_Rsquare = 2*0.945/GeV/GeV *0.5; // R^2 of w-meson - nucleon interactions
89  Wreggeon_C = 1.0; // Shower enhancement coefficient
90  if(PDGcode > 0) WParity = -1; // Parity +1 for Pbar P, and -1 for PP interactions
91  if(PDGcode < 0) WParity = +1;
92  }
93  else if ( absPDGcode == 211 || PDGcode == 111 ) { // Projectile is Pion
94  Cpr_pomeron = 1.352;
95  Ctr_pomeron = C_pomeron_N;
97 // KP
98  Gamma_pomeron_Pr = 0.89/GeV; // 0.85 -> 0.89 // Uzhi
99  Gamma_pomeron_Tr = Gamma_pomeron_N;
100  Gamma_pomeron = Gamma_pomeron_Pr * Gamma_pomeron_Tr;
101 
102  Rsquare_pomeron_Pr = 0.5/GeV/GeV;
103  Rsquare_pomeron_Tr = Rsquare_pomeron_N;
104  Rsquare_pomeron = Rsquare_pomeron_Pr + Rsquare_pomeron_Tr;
105 
106  Freggeon_Alpha = 0.7;
107  Freggeon_Gamma = 3.524/GeV/GeV;
108  Freggeon_Rsquare = 1.0/GeV/GeV;
109  Freggeon_Alphaprime = 0.8/GeV/GeV;
110  Freggeon_C = 1.0;
111  FParity = +1;
112 
113  Wreggeon_Alpha = 0.5;
114  Wreggeon_Gamma = 0.56/GeV/GeV; // 1.12 -> 0.56 Uzhi
115  Wreggeon_Rsquare = 9.19/GeV/GeV;
116  Wreggeon_Alphaprime = 0.9/GeV/GeV;
117  Wreggeon_C = 1.0;
118  if(PDGcode > 0) WParity = -1;
119  if(PDGcode < 0) WParity = +1;
120  }
121  else if ( absPDGcode == 321 || absPDGcode == 311 ||
122  PDGcode == 130 || PDGcode == 310 ) { // Projectile is Kaon
123 
124  Cpr_pomeron = 1.522;
125  Ctr_pomeron = C_pomeron_N;
127 
128  Gamma_pomeron_Pr = 1.312/GeV;
129  Gamma_pomeron_Tr = Gamma_pomeron_N;
130  Gamma_pomeron = Gamma_pomeron_Pr * Gamma_pomeron_Tr;
131 
132  Rsquare_pomeron_Pr = 0.31/GeV/GeV;
133  Rsquare_pomeron_Tr = Rsquare_pomeron_N;
134  Rsquare_pomeron = Rsquare_pomeron_Pr + Rsquare_pomeron_Tr;
135 
136  Freggeon_Alpha = 0.0;
137  Freggeon_Gamma = 0.0/GeV/GeV;
138  Freggeon_Rsquare = 1.0/GeV/GeV;
139  Freggeon_Alphaprime = 0.0/GeV/GeV;
140  Freggeon_C = 1.0;
141  FParity = +1;
142 
143  Wreggeon_Alpha = 0.0;
144  Wreggeon_Gamma = 0.0/GeV/GeV;
145  Wreggeon_Rsquare = 1.0/GeV/GeV;
146  Wreggeon_Alphaprime = 0.0/GeV/GeV;
147  Wreggeon_C = 1.0;
148  WParity = -1;
149  }
150  else if ( absPDGcode == 22 ) { // Projectile is Gamma
151 
152  Cpr_pomeron = 1.437;
153  Ctr_pomeron = C_pomeron_N;
155 
156  Gamma_pomeron_Pr = 1.415/GeV/GeV;
157  Gamma_pomeron_Tr = Gamma_pomeron_N;
158  Gamma_pomeron = Gamma_pomeron_Pr * Gamma_pomeron_Tr;
159 
160  Rsquare_pomeron_Pr = 0.51/GeV/GeV;
161  Rsquare_pomeron_Tr = Rsquare_pomeron_N;
162  Rsquare_pomeron = Rsquare_pomeron_Pr + Rsquare_pomeron_Tr;
163 
164  Freggeon_Alpha = 0.0;
165  Freggeon_Gamma = 0.0/GeV/GeV;
166  Freggeon_Rsquare = 1.0/GeV/GeV;
167  Freggeon_Alphaprime = 0.0/GeV/GeV;
168  Freggeon_C = 1.0;
169  FParity = +1;
170 
171  Wreggeon_Alpha = 0.0;
172  Wreggeon_Gamma = 0.0/GeV/GeV;
173  Wreggeon_Rsquare = 1.0/GeV/GeV;
174  Wreggeon_Alphaprime = 0.0/GeV/GeV;
175  Wreggeon_C = 1.0;
176  WParity = -1;
177  }
178  else { // Projectile is undefined, Nucleon assumed
179  Cpr_pomeron = C_pomeron_N;
180  Ctr_pomeron = C_pomeron_N;
182 
183  Gamma_pomeron_Pr = Gamma_pomeron_N;
184  Gamma_pomeron_Tr = Gamma_pomeron_N;
185  Gamma_pomeron = Gamma_pomeron_Pr * Gamma_pomeron_Tr;
186 
187  Rsquare_pomeron_Pr = Rsquare_pomeron_N;
188  Rsquare_pomeron_Tr = Rsquare_pomeron_N;
189  Rsquare_pomeron = Rsquare_pomeron_Pr + Rsquare_pomeron_Tr;
190 
191  Freggeon_Alpha = 0.723;
192  Freggeon_Gamma = 8.801/GeV/GeV;
193  Freggeon_Rsquare = 0.396/GeV/GeV;
194  Freggeon_Alphaprime = 1.324/GeV/GeV;
195  Freggeon_C = 1.0;
196  FParity = +1;
197 
198  Wreggeon_Alpha = 0.353;
199  Wreggeon_Gamma = 8.516/GeV/GeV;
200  Wreggeon_Rsquare = 24.40/GeV/GeV;
201  Wreggeon_Alphaprime = 1.5/GeV/GeV;
202  Wreggeon_C = 1.0;
203  WParity = -1;
204  }
205 /*
206 G4cout<<G4endl<<"Reggeon's parameters for Particle "<<particle->GetParticleName()<<" "<<PDGcode<<G4endl<<G4endl;
207 G4cout<<"Alpha_pomeron "<<Alpha_pomeron;
208 G4cout<<" Alphaprime_pomeron "<<Alphaprime_pomeron*GeV*GeV;
209 G4cout<<" S0_pomeron "<<S0_pomeron/GeV/GeV<<G4endl;
210 G4cout<<"Gamma_pomeron "<<Gamma_pomeron*GeV*GeV;
211 G4cout<<" Rsquare_pomeron "<<Rsquare_pomeron*GeV*GeV;
212 G4cout<<" C_pomeron "<<C_pomeron<<G4endl<<G4endl;
213 */
214 }
215 
217 
219 
221 
223 
225 {
226  Xtotal =0.; XtotalP=0.; XtotalR=0.;
227  Xelastic=0.; Xpr_Diff=0.; Xtr_Diff=0.; XDDiff=0.; G4double XDiff=0.;
228  Xinel =0.; Xnd=0.; XndP=0.; XndR=0.;
229 
230  G4double AmplitudeP(0.), AmplitudeR(0.);
231 
232  G4double B_max = 10.*fermi;
233  G4double dB = B_max/10000.;
234 
235  G4double B =-dB/2.;
236 
237  G4double chiP(0.), chiR(0.), chiRin(0.); // chiPin Pomeron inelastic phase is a data member
238  chiPin=0.;
239  for(G4int i=0; i<10000;i++)
240  {
241  B += dB;
242 
243  chiP = Chi_pomeron(1.,B); chiR = Chi_reggeon(1.,B);
244  chiPin = Chi_pomeron(2.,B); chiRin = Chi_reggeon(2.,B);
245 
246  AmplitudeP = (1.0/C_pomeron)*(1.0 - G4Exp(-chiP))*G4Exp(-chiR);
247  AmplitudeR = (1.0 - G4Exp(-chiR));
248 
249  Xtotal += 2 * (AmplitudeP + AmplitudeR) * B * dB;
250  XtotalP += 2 * (AmplitudeP + 0. ) * B * dB;
251  XtotalR += 2 * (0. + AmplitudeR) * B * dB;
252 
253  Xelastic += sqr(AmplitudeP + AmplitudeR) * B * dB;
254  Xpr_Diff += (Cpr_pomeron - 1.0) * sqr(AmplitudeP) * B * dB;
255  Xtr_Diff += (Ctr_pomeron - 1.0) * sqr(AmplitudeP) * B * dB;
256  XDiff += (Cpr_pomeron - 1.0) * (Ctr_pomeron - 1.0) * sqr(AmplitudeP) * B * dB;
257 
258 // ----------------------------------
259  AmplitudeP = (1.0/C_pomeron)*(1.0 - G4Exp(-chiPin))*G4Exp(-chiRin);
260  AmplitudeR = (1.0 - G4Exp(-chiRin));
261 
262  Xnd += (AmplitudeP + AmplitudeR) * B * dB;
263  XndP += (AmplitudeP + 0. ) * B * dB;
264  XndR += (0. + AmplitudeR) * B * dB;
265  }
266 
268  Xelastic *=twopi; Xpr_Diff *=twopi; Xtr_Diff *=twopi; XDiff *=twopi;
269  Xinel = Xtotal - Xelastic;
270  (void)Xinel; // To avoid compiler warning "variable not used"
271 
272  Xnd *=twopi; XndP *=twopi; XndR *=twopi;
273  XDDiff = XDiff-Xpr_Diff-Xtr_Diff;
274 
275 /*
276 G4cout<<"Total totalP totalR "<<Xtotal/millibarn <<" "<<XtotalP/millibarn <<" "<<XtotalR/millibarn<<" mb"<<G4endl;
277 G4cout<<"Elastic "<<Xelastic/millibarn <<G4endl;
278 G4cout<<"PrDiff TrDiff W_Diff "<<Xpr_Diff/millibarn<<" "<<Xtr_Diff/millibarn<<" "<<XDiff/millibarn<<G4endl;
279 G4cout<<"Inelastic "<<Xinel/millibarn <<G4endl;
280 G4cout<<"NonDiff Pom & Reg "<<Xnd/millibarn <<" "<<XndP/millibarn <<" "<<XndR/millibarn <<G4endl;
281 */
282 }
283 
285 {
287  G4double Eikonal = Mult * C_pomeron * Gamma_pomeron/R2 *
289  G4Exp(-sqr(B)/4.0/R2/hbarc_squared);
290  return Eikonal;
291 }
292 
294 {
297 
298  G4double Eikonal = Mult * FParity * Freggeon_C * Freggeon_Gamma/R2F *
300  G4Exp(-sqr(B)/4.0/R2F/hbarc_squared);
301 
302  Eikonal+= Mult * WParity * Wreggeon_C * Wreggeon_Gamma/R2W *
304  G4Exp(-sqr(B)/4.0/R2W/hbarc_squared);
305  return Eikonal;
306 }
307 
311 
316 
318 
322 
323 //----------------------------------------------------------------------------------------------
325  G4double & Pint,
326  G4double & Pprd, G4double & Ptrd, G4double & Pdd,
327  G4double & Pnd, G4double & Pnvr)
328 {
329  // Puprose of the method is a calculation of inelastic interaction probability (Pint),
330  // probability of projectile diffraction (Pprd),
331  // probability of target diffraction (Ptrd),
332  // probability of double diffraction (Pdd ),
333  // probability of non-diffractive inter. (Pnd ),
334  // probability of quark-exc. inter. (Pnvr),
335  // number of cutted pomerons (NcutPomerons).
336  // The input parameters are B - impact parameter, and Mode = All/WITHOUT_R/NON_DIFF
337  //
338  if( B > 2.* fermi ) { Pint=0.; Pprd=0.; Ptrd=0.; Pdd=0.; Pnd=0.; Pnvr=0.; return;}
339  // At large B for hN collisions it is better to return zero inter. probability
340 
341  G4double chiP = Chi_pomeron(1.,B); G4double chiR = Chi_reggeon(1.,B);
342  chiPin = Chi_pomeron(2.,B); G4double chiRin = Chi_reggeon(2.,B);
343  //chiPin is data member of the class
344 
345  G4double Exp_ChiR = G4Exp(-chiR);
346 
347  G4double AmplitudeP = (1.0/C_pomeron)*(1.0 - G4Exp(-chiP))*Exp_ChiR;
348  G4double AmplitudeR = (1.0 - Exp_ChiR);
349 
350  G4double AmplitudeP2, Apr_Diff, Atr_Diff, ADiff;
351 
352  //Aelastic = sqr(AmplitudeP + AmplitudeR);
353  AmplitudeP2 = sqr(AmplitudeP);
354  Apr_Diff = (Cpr_pomeron - 1.0) * AmplitudeP2;
355  Atr_Diff = (Ctr_pomeron - 1.0) * AmplitudeP2;
356  ADiff = (Cpr_pomeron - 1.0) * (Ctr_pomeron - 1.0) * AmplitudeP2;
357 
358 // ----------------------------------
359  Exp_ChiR = G4Exp(-chiRin);
360  AmplitudeP = (1.0/C_pomeron)*(1.0 - G4Exp(-chiPin))*Exp_ChiR;
361  AmplitudeR = (1.0 - Exp_ChiR);
362 
363  G4double And, AndP, AndR;
364 
365  And = (AmplitudeP + AmplitudeR);
366  AndP = (AmplitudeP + 0. );
367  AndR = (0. + AmplitudeR);
368 
369 // ----------------------------------
370  if( Mode == ALL)
371  {
372  Pint = Apr_Diff + Atr_Diff + ADiff + And;
373  Pprd = Apr_Diff/Pint; // Probability of projectile diffraction
374  Ptrd = Atr_Diff/Pint; // Probability of target diffraction
375  Pdd = ADiff /Pint; // Probability of double diffraction
376  Pnd = AndP /Pint; // Probability of non-diffractive inelastic
377  // interaction
378  Pnvr = AndR /Pint; // Probability of non-vacuum reggeon (nvr)
379  // inelastic interaction
380  }
381  else if( Mode == WITHOUT_R)
382  {
383  Pint = Apr_Diff + Atr_Diff + ADiff + AndP;
384  Pprd = Apr_Diff/Pint;
385  Ptrd = Atr_Diff/Pint;
386  Pdd = ADiff /Pint;
387  Pnd = AndP /Pint;
388  Pnvr = 0.;
389  }
390  else
391  {// Mode == NON_DIFF (of projectile)
392  Pint = Atr_Diff + AndP;
393  Pprd = 0.;
394  Ptrd = Atr_Diff/Pint;
395  Pdd = 0.;
396  Pnd = AndP /Pint;
397  Pnvr = 0.;
398  }
399 
400  return;
401 }
402 
403 G4int G4Reggeons::ncPomerons() // Non-complite Poisson distribution
404 {
405  if( chiPin < 0.001 ) return 0; // At small average multiplicity of cutted pomerons
406  // it is better to return 0 to avoid problems with
407  // calculation exactness.
408  G4double ksi = G4UniformRand() * (1.0-G4Exp(-chiPin)) * G4Exp(chiPin);
409  G4double Term = chiPin;
410  G4double Sum = Term;
411  G4int nCuts = 1;
412 
413  while( Sum < ksi)
414  {
415  nCuts++;
416  Term *= chiPin/(G4double) nCuts;
417  Sum += Term;
418  }
419 
420  return nCuts;
421 }
422 
G4double C_pomeron
Definition: G4Reggeons.hh:97
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
Definition: G4Exp.hh:183
G4double GetPrDiffX()
Definition: G4Reggeons.cc:313
G4double WParity
Definition: G4Reggeons.hh:95
G4double GetNDr_X()
Definition: G4Reggeons.cc:321
G4double chiPin
Definition: G4Reggeons.hh:103
G4double GetNDp_X()
Definition: G4Reggeons.cc:320
G4double Alpha_pomeron
Definition: G4Reggeons.hh:91
G4double GetTotalX()
Definition: G4Reggeons.cc:308
G4double S0_pomeron
Definition: G4Reggeons.hh:91
G4double GetInelX()
Definition: G4Reggeons.cc:317
G4double GetTrDiffX()
Definition: G4Reggeons.cc:314
G4double Gamma_pomeron
Definition: G4Reggeons.hh:91
G4double XndP
Definition: G4Reggeons.hh:107
G4double G4Log(G4double x)
Definition: G4Log.hh:230
double S(double temp)
G4double Get_Cprojectile()
Definition: G4Reggeons.cc:216
G4double Gamma_pomeronHard
Definition: G4Reggeons.hh:92
G4double FParity
Definition: G4Reggeons.hh:94
G4double XtotalP
Definition: G4Reggeons.hh:105
G4double Xtotal
Definition: G4Reggeons.hh:105
G4double Freggeon_Rsquare
Definition: G4Reggeons.hh:94
double G4double
Definition: G4Types.hh:76
static constexpr double fermi
Definition: G4SIunits.hh:103
G4double Freggeon_Alpha
Definition: G4Reggeons.hh:94
static G4Pow * GetInstance()
Definition: G4Pow.cc:57
G4double powA(G4double A, G4double y) const
Definition: G4Pow.hh:242
G4double Xtr_Diff
Definition: G4Reggeons.hh:106
G4double GetElasticX()
Definition: G4Reggeons.cc:312
#define G4UniformRand()
Definition: Randomize.hh:53
static constexpr double twopi
Definition: G4SIunits.hh:76
void SetS(G4double S)
Definition: G4Reggeons.cc:222
G4double Xelastic
Definition: G4Reggeons.hh:106
G4double GetTotalXp()
Definition: G4Reggeons.cc:309
G4Reggeons(const G4ParticleDefinition *)
Definition: G4Reggeons.cc:38
G4double Xnd
Definition: G4Reggeons.hh:107
G4double Xpr_Diff
Definition: G4Reggeons.hh:106
G4double Rsquare_pomeron
Definition: G4Reggeons.hh:91
G4double Cpr_pomeron
Definition: G4Reggeons.hh:98
G4double Wreggeon_C
Definition: G4Reggeons.hh:95
G4double Freggeon_Alphaprime
Definition: G4Reggeons.hh:94
G4double Alphaprime_pomeron
Definition: G4Reggeons.hh:91
G4double GetDDiffX()
Definition: G4Reggeons.cc:315
G4double GetTotalXr()
Definition: G4Reggeons.cc:310
int G4int
Definition: G4Types.hh:78
G4double XtotalR
Definition: G4Reggeons.hh:105
G4double GetND_X()
Definition: G4Reggeons.cc:319
G4double Chi_pomeron(G4double Mult, G4double B)
Definition: G4Reggeons.cc:284
G4double Freggeon_Gamma
Definition: G4Reggeons.hh:94
T sqr(const T &x)
Definition: templates.hh:145
G4double Chi_reggeon(G4double Mult, G4double B)
Definition: G4Reggeons.cc:293
void CalculateXs()
Definition: G4Reggeons.cc:224
G4double Sint
Definition: G4Reggeons.hh:101
G4double XndR
Definition: G4Reggeons.hh:107
static constexpr double hbarc_squared
double B(double temperature)
G4double Wreggeon_Alphaprime
Definition: G4Reggeons.hh:95
typedef void(XMLCALL *XML_ElementDeclHandler)(void *userData
static constexpr double GeV
Definition: G4SIunits.hh:217
G4double Xinel
Definition: G4Reggeons.hh:107
G4double Wreggeon_Rsquare
Definition: G4Reggeons.hh:95
void GetProbabilities(G4double B, G4int Mode, G4double &Pint, G4double &Pprd, G4double &Ptrd, G4double &Pdd, G4double &Pnd, G4double &Pnvr)
Definition: G4Reggeons.cc:324
G4double Wreggeon_Gamma
Definition: G4Reggeons.hh:95
G4double Get_Ctarget()
Definition: G4Reggeons.cc:218
G4double Alpha_pomeronHard
Definition: G4Reggeons.hh:92
G4double Freggeon_C
Definition: G4Reggeons.hh:94
G4double Wreggeon_Alpha
Definition: G4Reggeons.hh:95
G4double Ctr_pomeron
Definition: G4Reggeons.hh:99
G4double XDDiff
Definition: G4Reggeons.hh:106
G4int ncPomerons()
Definition: G4Reggeons.cc:403