34 #define INCLXX_IN_GEANT4_MODE 1
109 return 77./(plab+1.5);
115 return 6.3555*std::exp(-3.2481*alp-0.377*alp*alp);
117 else if (plab < 0.800) {
118 return (33.+196.*std::sqrt(std::pow(
std::abs(plab-0.95),5)));
121 return 31./std::sqrt(plab);
126 return 34.*std::pow(plab/0.4, (-2.104));
128 else if (plab < 0.800) {
129 return (23.5+1000.*std::pow(plab-0.7, 4));
132 return (1250./(50.+plab)-4.*std::pow(plab-1.3, 2));
138 return 34.*std::pow(plab/0.4, (-2.104));
140 else if (plab < 0.800) {
141 return (23.5+1000.*std::pow(plab-0.7, 4));
144 return (1250./(50.+plab)-4.*std::pow(plab-1.3, 2));
153 const G4double momentumGeV = 0.001 * pLab;
158 if(isospin==2 || isospin==-2) {
160 xs = (41.0 + (60.0*momentumGeV - 54.0)*std::exp(-1.2*momentumGeV) - 77.0/(momentumGeV + 1.5));
161 }
else if(pLab >= 1500.0 && pLab < 2000.0) {
162 xs = (41.0 + 60.0*(momentumGeV - 0.9)*std::exp(-1.2*momentumGeV) - 1250.0/(momentumGeV+50.0)+ 4.0*std::pow(momentumGeV - 1.3, 2));
163 }
else if(pLab < 1500.0) {
164 xs = (23.5 + 24.6/(1.0 + std::exp(-10.0*momentumGeV + 12.0))
165 -1250.0/(momentumGeV +50.0)+4.0*std::pow(momentumGeV - 1.3,2));
167 }
else if(isospin==0) {
169 xs = (42.0 - 77.0/(momentumGeV + 1.5));
170 }
else if(pLab >= 1000.0 && pLab < 2000.0) {
171 xs = (24.2 + 8.9*momentumGeV - 31.1/std::sqrt(momentumGeV));
172 }
else if(pLab < 1000.0) {
173 xs = (33.0 + 196.0*std::sqrt(std::pow(
std::abs(momentumGeV - 0.95),5))
174 -31.1/std::sqrt(momentumGeV));
178 if(xs < 0.0)
return 0.0;
185 return -2.33730e-06*std::pow(x, 3)+1.13819e-02*std::pow(x,2)
186 -1.83993e+01*x+9893.4;
187 }
else if(x > 1750.0 && x <= 2175.0) {
188 return 1.13531e-06*std::pow(x, 3)-6.91694e-03*std::pow(x, 2)
189 +1.39907e+01*x-9360.76;
191 return -3.18087*std::log(x)+52.9784;
198 return 0.00120683*(x-1372.52)*(x-1372.52)+26.2058;
199 }
else if(x > 1475.0 && x <= 1565.0) {
200 return 1.15873e-05*x*x+49965.6/((x-1519.59)*(x-1519.59)+2372.55);
201 }
else if(x > 1565.0 && x <= 2400.0) {
202 return 34.0248+43262.2/((x-1681.65)*(x-1681.65)+1689.35);
203 }
else if(x > 2400.0 && x <= 7500.0) {
204 return 3.3e-7*(x-7500.0)*(x-7500.0)+24.5;
224 return inelastic +
elastic(p1, p2);
237 if(x>10000.)
return 0.0;
245 }
else if(particle2->
isPion()) {
256 G4double q2=(y-1076.0*1076.0)*(y-800.0*800.0)/y/4.0;
260 G4double q3 = std::pow(std::sqrt(q2),3);
262 G4double spnResult = 326.5/(std::pow((x-1215.0-ramass)*2.0/(110.0-ramass), 2)+1.0);
263 spnResult = spnResult*(1.0-5.0*ramass/1215.0);
265 spnResult = spnResult*f3*cg/6.0;
267 if(x < 1200.0 && spnResult < 5.0) {
273 if((ind2t3 == 1 && ipit3 == 2) || (ind2t3 == -1 && ipit3 == -2))
275 else if((ind2t3 == 1 && ipit3 == -2) || (ind2t3 == -1 && ipit3 == 2))
288 if(isospin==4 || isospin==-4)
return 0.0;
302 if(Ecm <= 938.3 + deltaMass) {
306 if(Ecm < 938.3 + deltaMass + 2.0) {
307 Ecm = 938.3 + deltaMass + 2.0;
319 result *= 3.*(32.0 + isospin * isospin * (deltaIsospin * deltaIsospin - 5))/64.0;
320 result /= 1.0 + 0.25 * isospin * isospin;
350 return 5.5e-6 * x/(7.7 +
x);
352 return (5.34 + 0.67*(x - 2.0)) * 1.0
e-6;
357 return b/(1.0 + std::exp(-(x - 0.45)/0.05));
358 }
else if(pl < 1100.0) {
359 return (9.87 - 4.88 * x) * 1.0e-6;
361 return (3.68 + 0.76*x) * 1.0e-6;