ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4UPiNuclearCrossSection.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4UPiNuclearCrossSection.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 // Calculation of the total, elastic and inelastic cross-sections
28 // based on Barashenkov parametrisations of pion data
29 //
30 // 16.08.06 V.Ivanchenko - first implementation on base of
31 // J.P Wellisch class G4PiNuclearCrossSection
32 // 22.01.07 V.Ivanchenko - add cross section interfaces with Z and A
33 // 05.03.07 V.Ivanchenko - fix weight for interpolation
34 // 13.03.07 V.Ivanchenko - cleanup at low energies
35 // 11.09.09 V.Ivanchenko - fixed bug in interpolation
36 //
37 
38 #include "G4UPiNuclearCrossSection.hh"
39 #include "G4SystemOfUnits.hh"
40 #include "G4LPhysicsFreeVector.hh"
41 #include "G4PionMinus.hh"
42 #include "G4PionPlus.hh"
43 #include "G4PhysicsTable.hh"
44 #include "G4NistManager.hh"
45 
46 G4int G4UPiNuclearCrossSection::theZ[NZ] =
47 {2,4,6,7,8,11,13,20,26,29,42,48,50,74,82,92};
48 G4int G4UPiNuclearCrossSection::idxZ[93] = {0};
49 G4double G4UPiNuclearCrossSection::theA[NZ] = {0.0};
50 G4double G4UPiNuclearCrossSection::APower[93] = {0.0};
51 
52 G4PhysicsTable* G4UPiNuclearCrossSection::piPlusElastic = nullptr;
53 G4PhysicsTable* G4UPiNuclearCrossSection::piPlusInelastic = nullptr;
54 G4PhysicsTable* G4UPiNuclearCrossSection::piMinusElastic = nullptr;
55 G4PhysicsTable* G4UPiNuclearCrossSection::piMinusInelastic = nullptr;
56 
57 #ifdef G4MULTITHREADED
58 G4Mutex G4UPiNuclearCrossSection::pionUXSMutex = G4MUTEX_INITIALIZER;
59 #endif
60 
61 G4UPiNuclearCrossSection::G4UPiNuclearCrossSection()
62  : G4VCrossSectionDataSet("G4UPiNuclearCrossSection")
63 {
64  isMaster = false;
65  piPlus = G4PionPlus::PionPlus();
66  piMinus = G4PionMinus::PionMinus();
67 
68  aPower = 0.75;
69  elow = 20.0*CLHEP::MeV;
70 }
71 
72 G4UPiNuclearCrossSection::~G4UPiNuclearCrossSection()
73 {
74  if(isMaster) {
75  if(piPlusElastic) {
76  piPlusElastic->clearAndDestroy();
77  delete piPlusElastic;
78  piPlusElastic = nullptr;
79  }
80  if(piPlusInelastic) {
81  piPlusInelastic->clearAndDestroy();
82  delete piPlusInelastic;
83  piPlusInelastic = nullptr;
84  }
85  if(piMinusElastic) {
86  piMinusElastic->clearAndDestroy();
87  delete piMinusElastic;
88  piMinusElastic = nullptr;
89  }
90  if(piMinusInelastic) {
91  piMinusInelastic->clearAndDestroy();
92  delete piMinusInelastic;
93  piMinusInelastic = nullptr;
94  }
95  }
96 }
97 
98 G4bool
99 G4UPiNuclearCrossSection::IsElementApplicable(const G4DynamicParticle*,
100  G4int Z, const G4Material*)
101 {
102  return (1 < Z);
103 }
104 
105 G4double G4UPiNuclearCrossSection::Interpolate(
106  G4int Z, G4int A, G4double e, const G4PhysicsTable* table) const
107 {
108  G4double res = 0.0;
109  G4double ekin = std::max(e, elow);
110  G4int iz = std::min(Z, 92);
111  G4int idx = idxZ[iz];
112  size_t jdx = (size_t)(std::max(ekin - elow, 0.0)*0.06);
113  //G4cout << "Interpolate: Z= " << iz << " A= " << A << " idx= " << idx
114  // << " jdx= " << jdx << " Ekin= " << ekin << G4endl;
115  if(idx < 0 || 2 == iz) {
116  res = ((*table)[std::abs(idx)])->Value(ekin, jdx);
117  //G4cout << "1: jdx= " << jdx << G4endl;
118  } else {
119  G4int iz2 = theZ[idx];
120  G4double x2 = (((*table)[idx])->Value(ekin, jdx))*APower[iz]/APower[iz2];
121  //G4cout << "2: jdx= " << jdx << G4endl;
122  G4int iz1 = theZ[idx-1];
123  G4double x1 = (((*table)[idx-1])->Value(ekin, jdx))*APower[iz]/APower[iz1];
124  G4double w1 = ((G4double)A - theA[idx-1])/(theA[idx] - theA[idx-1]);
125  res = w1*x2 + (1.0 - w1)*x1;
126  }
127  //G4cout << " res(nb)= " << res/CLHEP::barn << G4endl;
128  return res;
129 }
130 
131 void G4UPiNuclearCrossSection::AddDataSet(const G4String& p,
132  const G4double* tot,
133  const G4double* in,
134  const G4double* e,
135  G4int n)
136 {
137  G4LPhysicsFreeVector* pvin = new G4LPhysicsFreeVector(n,e[0]*GeV,e[n-1]*GeV);
138  pvin->SetSpline(true);
139  G4LPhysicsFreeVector* pvel = new G4LPhysicsFreeVector(n,e[0]*GeV,e[n-1]*GeV);
140  pvel->SetSpline(true);
141  for(G4int i=0; i<n; ++i) {
142  pvin->PutValues(i,e[i]*GeV,in[i]*millibarn);
143  pvel->PutValues(i,e[i]*GeV,std::max(0.0,(tot[i]-in[i])*millibarn));
144  }
145  if(p == "pi+") {
146  piPlusInelastic->push_back(pvin);
147  piPlusElastic->push_back(pvel);
148  } else {
149  piMinusInelastic->push_back(pvin);
150  piMinusElastic->push_back(pvel);
151  }
152 }
153 
154 void G4UPiNuclearCrossSection::DumpPhysicsTable(const G4ParticleDefinition& p)
155 {
156  if(&p == piPlus) {
157  G4cout << "### G4UPiNuclearCrossSection Elastic data for pi+" << G4endl;
158  G4cout << *piPlusElastic << G4endl;
159  G4cout << "### G4UPiNuclearCrossSection Inelastic data for pi+" << G4endl;
160  G4cout << *piPlusInelastic << G4endl;
161  } else if(&p == piMinus) {
162  G4cout << "### G4UPiNuclearCrossSection Elastic data for pi-" << G4endl;
163  G4cout << *piMinusElastic << G4endl;
164  G4cout << "### G4UPiNuclearCrossSection Inelastic data for pi-" << G4endl;
165  G4cout << *piMinusInelastic << G4endl;
166  }
167 }
168 
169 void G4UPiNuclearCrossSection::BuildPhysicsTable(const G4ParticleDefinition& p)
170 {
171  //G4cout << "G4UPiNuclearCrossSection::BuildPhysicsTable for "
172  // << p.GetParticleName() << " " << piPlusElastic << G4endl;
173  if(piPlusElastic) { return; }
174 
175  if(&p != piPlus && &p != piMinus) {
176  G4ExceptionDescription ed;
177  ed << "This cross section is applicable only to pions and not to "
178  << p.GetParticleName() << G4endl;
179  G4Exception("G4UPiNuclearCrossSection::BuildPhysicsTable", "had001",
180  FatalException, ed);
181  return;
182  }
183 
184  if(!piPlusElastic) {
185 #ifdef G4MULTITHREADED
186  G4MUTEXLOCK(&pionUXSMutex);
187  if(!piPlusElastic) {
188 #endif
189  isMaster = true;
190 #ifdef G4MULTITHREADED
191  }
192  G4MUTEXUNLOCK(&pionUXSMutex);
193 #endif
194  }
195  if(!isMaster) { return; }
196 
197  G4NistManager* nist = G4NistManager::Instance();
198  G4Pow* g4pow = G4Pow::GetInstance();
199  for(G4int i=0; i<NZ; ++i) {
200  theA[i] = nist->GetAtomicMassAmu(theZ[i]);
201  }
202  for(G4int i=1; i<93; ++i) {
203  APower[i] = g4pow->powA(nist->GetAtomicMassAmu(i),aPower);
204  }
205  G4int idx = 1;
206  for(G4int i=3; i<93; ++i) {
207  if(theZ[idx] == i) {
208  idxZ[i] = -idx;
209  ++idx;
210  } else {
211  idxZ[i] = idx;
212  }
213  }
214 
215  piPlusElastic = new G4PhysicsTable();
216  piPlusInelastic = new G4PhysicsTable();
217  piMinusElastic = new G4PhysicsTable();
218  piMinusInelastic = new G4PhysicsTable();
219 
220  LoadData();
221  //G4cout << "### Initialisation done " << G4endl;
222 }
223 
224 void G4UPiNuclearCrossSection::LoadData()
225 {
226  static const G4double e1[38] = {
227  0.02, 0.04, 0.06, 0.08, 0.1, 0.12, 0.13, 0.14, 0.15, 0.16,
228  0.17, 0.18, 0.19, 0.2, 0.22,0.24, 0.26, 0.28, 0.3, 0.35,
229  0.4, 0.45, 0.5, 0.55, 0.6, 0.7, 0.8, 0.9, 1, 2,
230  3, 5, 10, 20, 50, 100, 500, 1000};
231  static const G4double e2[39] = {
232  0.02, 0.04, 0.06, 0.08, 0.1, 0.11, 0.12, 0.13, 0.14, 0.15,
233  0.16, 0.17, 0.18, 0.2, 0.22, 0.24, 0.26, 0.28, 0.3, 0.35,
234  0.4, 0.45, 0.5, 0.55, 0.575,0.6, 0.7, 0.8, 0.9, 1,
235  2, 3, 5, 10, 20, 50, 100, 500, 1000};
236  static const G4double e3[31] = {
237  0.02, 0.04, 0.06, 0.08, 0.1, 0.12, 0.14, 0.16, 0.18, 0.2,
238  0.22, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.6, 0.7, 0.8,
239  0.9, 1, 2, 3, 5, 10, 20, 50, 100, 500, 1000};
240  static const G4double e4[32] = {
241  0.02, 0.04, 0.06, 0.08, 0.1, 0.12, 0.14, 0.16, 0.18, 0.2,
242  0.22, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.7,
243  0.8, 0.9, 1, 2, 3, 5, 10, 20, 50, 100, 500, 1000};
244  static const G4double e5[34] = {
245  0.02, 0.04, 0.05, 0.06, 0.07, 0.08, 0.09, 0.1, 0.12, 0.14,
246  0.16, 0.18, 0.2, 0.22, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5,
247  0.6, 0.7, 0.8, 0.9, 1, 2, 3, 5, 10, 20, 50, 100, 500, 1000};
248  static const G4double e6[35] = {
249  0.02, 0.04, 0.05, 0.06, 0.07, 0.08, 0.09, 0.1, 0.12, 0.14,
250  0.16, 0.18, 0.2, 0.22, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5,
251  0.55, 0.6, 0.7, 0.8, 0.9, 1, 2, 3, 5, 10, 20, 50, 100, 500, 1000};
252 
253  static const G4double he_t[38] = {
254  40, 70, 108, 152, 208, 276, 300, 320, 329, 333,
255  332, 328, 322, 310, 288, 260, 240, 216, 196, 144,
256  125, 112, 108.5, 109, 110.5, 117, 123, 128.5, 135, 110,
257  96, 87, 85, 83.5, 83.5, 83.5, 83.5, 83.5};
258  static const G4double he_in[38] = {
259  18, 38, 62, 98, 136, 176, 190, 200, 209, 212,
260  212, 208, 204, 196, 176, 164, 150, 134, 124, 97.5,
261  90, 85,82.5, 83.5, 86.5, 93, 97.5,100, 102, 83,
262  77, 75, 74, 72.5, 72.5, 72.5, 72.5, 72.5};
263  static const G4double be_m_t[38] = {
264  150, 210, 294, 396, 520, 600, 623, 635, 642, 640,
265  630, 615, 600, 576, 540, 504, 470, 435, 400, 340,
266  294, 258, 236, 230, 233, 244, 257, 270, 276, 250,
267  230, 215, 205, 194, 188, 186, 186, 186};
268  static const G4double be_m_in[38] = {
269  90, 126, 177, 240, 320, 380, 400, 410, 414, 410,
270  400, 387, 371, 360, 333, 312, 285, 260, 237, 216,
271  198, 187, 182, 180, 182, 187, 193, 203, 207, 179,
272  172, 165, 159, 155, 144, 144, 144, 144};
273  static const G4double be_p_t[38] = {
274  96, 150, 222, 320, 430, 514, 545, 565, 574, 574,
275  564, 552, 535, 522, 490, 462, 432, 398, 367, 314,
276  276, 248, 232, 230, 233, 244, 257, 270, 276, 250,
277  230, 215, 205, 194, 188, 186, 186, 186};
278  static const G4double be_p_in[38] = {
279  60, 95, 142, 194, 262, 319, 345, 361, 364, 364,
280  354, 350, 330, 319, 298, 280, 258, 237, 216, 200,
281  189, 183, 182, 180, 182, 187, 193, 203, 207, 179,
282  172, 165, 159, 155, 144, 144, 144, 144};
283 
284  static const G4double c_m_t[39] = {
285  204, 260, 366, 517, 630, 673, 694, 704, 710, 711,
286  706, 694, 676, 648, 616, 584, 548, 518, 489, 426,
287  376, 342, 323, 310, 312, 313, 319, 333, 342, 348,
288  310, 290, 268, 250, 245, 237, 234, 234, 234};
289  static const G4double c_m_in[39] = {
290  128, 160, 224, 315, 388, 416, 430, 438, 444, 445,
291  440, 432, 416, 400, 380, 354, 320, 304, 288, 264,
292  246, 240, 233, 232, 233, 234, 238, 246, 252, 256,
293  220, 210, 198, 187, 183, 176, 174, 174, 174};
294  static const G4double c_p_t[39] = {
295  140, 192, 294, 428, 594, 642, 662, 687, 685, 688,
296  684, 672, 656, 630, 598, 567, 533, 504, 474, 416,
297  369, 336, 319, 310, 312, 313, 319, 333, 342, 348,
298  310, 290, 268, 250, 245, 237, 234, 234, 234};
299  static const G4double c_p_in[39] = {
300  94, 132, 184, 260, 370, 398, 408, 420, 426, 428,
301  424, 416, 400, 386, 366, 340, 308, 294, 280, 257,
302  241, 236, 231, 232, 233, 234, 238, 246, 252, 256,
303  220, 210, 198, 187, 183, 176, 174, 174, 174};
304  static const G4double n_m_t[39] = {
305  246, 308, 424, 590, 729, 776, 800, 821, 822, 817,
306  800, 778, 768, 728, 690, 654, 615, 584, 556, 480,
307  430, 393, 373, 367, 368, 370, 375, 388, 390, 397,
308  364, 337, 310, 291, 275, 268, 268, 268, 268};
309  static const G4double n_m_in[39] = {
310  155, 188, 256, 360, 456, 492, 512, 526, 526, 520,
311  504, 491, 475, 450, 425, 396, 376, 360, 340, 300,
312  282, 270, 265, 265, 266, 268, 273, 280, 288, 288,
313  256, 237, 226, 218, 208, 202, 202, 202, 202};
314  static const G4double n_p_t[39] = {
315  150, 212, 328, 500, 680, 735, 762, 781, 782, 779,
316  770, 748, 740, 706, 672, 633, 600, 569, 541, 467,
317  419, 385, 368, 364, 366, 368, 375, 388, 390, 397,
318  364, 337, 310, 291, 275, 268, 268, 268, 268};
319  static const G4double n_p_in[39] = {
320  90, 140, 208, 300, 426, 467, 490, 504, 504, 500,
321  484, 474, 460, 437, 413, 381, 365, 350, 330, 292,
322  276, 267, 263, 264, 265, 267, 273, 280, 288, 288,
323  256, 237, 226, 218, 208, 202, 202, 202, 202};
324 
325  static const G4double o_m_t[31] = {
326  280, 360, 500, 685, 812, 861, 870, 865, 835, 800,
327  755, 700, 600, 537, 493, 468, 441, 436, 443, 449,
328  460, 463, 432, 385, 350, 325, 312, 307, 303, 303, 303};
329  static const G4double o_m_in[31] = {
330  190, 207, 300, 420, 500, 540, 550, 542, 520, 490,
331  460, 423, 360, 339, 321, 314, 312, 314, 319, 324,
332  328, 330, 300, 275, 250, 240, 229, 225, 222, 222, 222};
333  static const G4double o_p_t[31] = {
334  170, 240, 390, 570, 740, 818, 830, 822, 800, 765,
335  725, 675, 585, 525, 483, 458, 444, 447, 453, 449,
336  460, 463, 432, 385, 350, 325, 312, 307, 303, 303, 303};
337  static const G4double o_p_in[31] = {
338  100, 145, 240, 340, 470, 518, 530, 522, 505, 477,
339  448, 412, 350, 330, 316, 310, 308, 311, 317, 324,
340  328, 330, 300, 275, 250, 240, 229, 225, 222, 222, 222};
341  static const G4double na_m_t[31] = {
342  450, 545, 705, 910, 1020, 1075, 1087, 1080, 1042, 987,
343  943, 885, 790, 700, 650, 610, 585, 575, 585, 595,
344  600, 610, 556, 524, 494, 458, 445, 429, 427, 427, 427};
345  static const G4double na_m_in[31] = {
346  275, 315, 413, 545, 620, 660, 670, 662, 630, 593,
347  570, 520, 465, 420, 410, 395, 390, 400, 410, 418,
348  420, 422, 372, 348, 330, 320, 310, 294, 292, 292, 292};
349  static const G4double na_p_t[31] = {
350  210, 320, 530, 795, 960, 1035, 1050, 1040, 1007, 957,
351  918, 865, 773, 685, 636, 598, 575, 565, 578, 590,
352  598, 610, 556, 524, 494, 458, 445, 429, 427, 427, 427};
353  static const G4double na_p_in[31] = {
354  115, 210, 340, 495, 585, 630, 645, 637, 605, 572,
355  550, 505, 455, 410, 401, 388, 383, 393, 405, 414,
356  418, 422, 372, 348, 330, 320, 310, 294, 292, 292, 292};
357  static const G4double al_m_t[31] = {
358  532, 637, 832, 1057, 1207, 1230, 1210, 1174, 1133, 1095,
359  1038, 970, 890, 807, 750, 710, 675, 665, 670, 673,
360  678, 682, 618, 574, 546, 520, 507, 495, 488, 488, 488};
361  static const G4double al_m_in[31] = {
362  300, 360, 495, 665, 750, 765, 750, 730, 700, 660,
363  615, 570, 520, 490, 470, 450, 448, 450, 450, 452,
364  456, 460, 408, 392, 376, 356, 347, 338, 332, 332, 332};
365  static const G4double al_p_t[31] = {
366  225, 350, 616, 945, 1122, 1175, 1157, 1128, 1088, 1045,
367  988, 935, 870, 787, 730, 690, 660, 652, 660, 668,
368  678, 682, 618, 574, 546, 520, 507, 495, 488, 488, 488};
369  static const G4double al_p_in[31] = {
370  120, 238, 390, 610, 712, 735, 720, 703, 655, 635,
371  590, 550, 505, 475, 455, 438, 440, 445, 445, 450,
372  456, 460, 408, 392, 376, 356, 347, 338, 332, 332, 332};
373 
374  static const G4double ca_m_t[31] = {
375  800, 980, 1240, 1460, 1570, 1600, 1580, 1535, 1475, 1425,
376  1375,1295, 1200, 1083, 1000, 948, 915, 895, 900, 908,
377  915, 922, 856, 795, 740, 705, 682, 660, 660, 660, 660};
378  static const G4double ca_m_in[31] = {
379  470, 550, 620, 860, 955, 960, 920, 860, 820, 780,
380  740, 665, 637, 615, 600, 590, 580, 580, 600, 608,
381  610, 615, 550, 525, 510, 488, 470, 450, 450, 450, 450};
382  static const G4double ca_p_t[31] = {
383  275, 445, 790, 1195, 1440, 1485, 1475, 1435, 1385, 1335,
384  1295,1245,1160, 1050, 970, 923, 895, 877, 887, 897,
385  904, 913, 855, 795, 740, 705, 682, 660, 660, 660, 660};
386  static const G4double ca_p_in[31] = {
387  160, 315, 500, 745, 870, 905, 900, 860, 810, 770,
388  740, 710, 640, 617, 595, 585, 575, 575, 590, 600,
389  602, 608, 550, 525, 510, 488, 470, 450, 450, 450, 450};
390  // ca data may have typo
391 
392  static const G4double fe_m_t[32] = {
393  1175, 1363, 1670, 1950, 2050, 2040, 1975, 1886, 1834, 1773,
394  1720, 1635, 1474, 1380, 1269, 1225, 1182, 1162, 1159, 1162,
395  1178, 1190, 1197, 1102, 1135, 975, 945, 925, 905, 905,
396  905, 905};
397  static const G4double fe_m_in[32] = {
398  625, 725, 910, 1180, 1275, 1250, 1200, 1150, 1100, 1040,
399  995, 925, 825, 810, 780, 760, 745, 740, 740, 740,
400  750, 760, 765, 690, 660, 635, 615, 600, 585, 585,
401  585, 585};
402  static const G4double fe_p_t[32] = {
403  330, 575, 1010, 1500, 1837, 1875, 1820, 1751, 1691, 1636,
404  1690, 1450, 1396, 1305, 1219, 1190, 1148, 1138, 1134, 1144,
405  1163, 1175, 1183, 1198, 1135, 975, 945, 925, 905, 905,
406  905, 905};
407  static const G4double fe_p_in[32] = {
408  210, 410, 707, 1010, 1125, 1150, 1100, 1070, 1010, 960,
409  920, 776, 780, 760, 750, 740, 720, 725, 725, 730,
410  740, 750, 755, 690, 660, 635, 615, 600, 585, 585,
411  585, 585};
412  static const G4double cu_m_t[32] = {
413  1400, 1600, 1875, 2088, 2200, 2220, 2175, 2125, 2075, 2012,
414  1950, 1855, 1670, 1530, 1430, 1370, 1315, 1315, 1315, 1330,
415  1345, 1360, 1365, 1250, 1185, 1128, 1070, 1035, 1010, 1010,
416  1010, 1010};
417  static const G4double cu_m_in[32] = {
418  725, 840, 1020, 1200, 1295, 1300, 1267, 1240, 1213, 1175,
419  1125, 1042, 950, 900, 860, 840, 830, 832, 835, 840,
420  850, 860, 865, 785, 735, 705, 680, 650, 630, 630,
421  630, 630};
422  static const G4double cu_p_t[32] = {
423  355, 605, 1120, 1630, 1940, 2010, 2010, 1980, 1925, 1895,
424  1830, 1730, 1585, 1490, 1400, 1340, 1290, 1290, 1290, 1310,
425  1330, 1345, 1350, 1240, 1185, 1128, 1070, 1035, 1010, 1010,
426  1010, 1010};
427  static const G4double cu_p_in[32] = {
428  230, 425, 780, 1025, 1155, 1190, 1190, 1180, 1125, 1100,
429  1050, 1000, 900, 870, 835, 815, 810, 812, 815, 825,
430  840, 850, 855, 780, 735, 705, 680, 650, 630, 630,
431  630, 630};
432 
433  static const G4double mo_m_t[34] = {
434  2430, 2610, 2710, 2790, 2880, 2940, 2965, 2970, 2970, 2920,
435  2840, 2720, 2570, 2500, 2365, 2200, 2050, 1926, 1825, 1768,
436  1749, 1750, 1778, 1789, 1808, 1690, 1645, 1530, 1492, 1450,
437  1425, 1425, 1425, 1425};
438  static const G4double mo_m_in[34] = {
439  925, 1125, 1250, 1375, 1500, 1600, 1680, 1750, 1770, 1730,
440  1660, 1580, 1500, 1450, 1330, 1250, 1190, 1140, 1100, 1075,
441  1075, 1070, 1088, 1095, 1110, 1035, 1005, 940, 917, 880,
442  860, 860, 860, 860};
443  static const G4double mo_p_t[34] = {
444  410, 730, 1110, 1530, 1920, 2200, 2385, 2520, 2600, 2630,
445  2575, 2470, 2320, 2285, 2185, 2053, 1945, 1852, 1776, 1719,
446  1710, 1716, 1746, 1759, 1778, 1675, 1645, 1530, 1492, 1450,
447  1425, 1425, 1425, 1425};
448  static const G4double mo_p_in[34] = {
449  270, 540, 825, 975, 1140, 1285, 1400, 1480, 1555, 1580,
450  1525, 1470, 1360, 1340, 1255, 1160, 1120, 1085, 1060, 1045,
451  1045, 1045, 1065, 1075, 1090, 1025, 1005, 940, 917, 880,
452  860, 860, 860, 860};
453  static const G4double cd_m_t[34] = {
454  3060, 3125, 3170, 3220, 3255, 3280, 3290, 3260, 3270, 3200,
455  3120, 3080, 3090, 2920, 2810, 2640, 2362, 2230, 2115, 2050,
456  2020, 2025, 2040, 2070, 2100, 1900, 1795, 1740, 1675, 1645,
457  1625, 1620, 1620, 1620};
458  static const G4double cd_m_in[34]= {
459  1025, 1275, 1440, 1625, 1740, 1800, 1880, 1920, 1980, 1920,
460  1850, 1810, 1720, 1650, 1560, 1450, 1330, 1290, 1245, 1210,
461  1200, 1200, 1205, 1205, 1230, 1130, 1085, 1060, 1000, 985,
462  975, 970, 970, 970};
463  static const G4double cd_p_t[34] = {
464  455, 780, 1170, 1700, 2120, 2400, 2600, 2720, 2820, 2840,
465  2800, 2760, 2720, 2640, 2560, 2450, 2252, 2130, 2035, 1985,
466  1970, 1975, 2005, 2035, 2070, 1880, 1795, 1740, 1675, 1645,
467  1625, 1620, 1620, 1620};
468  static const G4double cd_p_in[34] = {
469  310, 580, 880, 1060, 1270, 1400, 1530, 1610, 1660, 1680,
470  1640, 1600, 1560, 1500, 1430, 1330, 1280, 1230, 1200, 1180,
471  1170, 1175, 1180, 1180, 1210, 1120, 1085, 1060, 1000, 985,
472  975, 970, 970, 970};
473 
474  static const G4double sn_m_t[35] = {
475  3000, 3180, 3250, 3300, 3300, 3410, 3470, 3450, 3410, 3350,
476  3280, 3200, 3120, 3050, 2900, 2630, 2500, 2325, 2190, 2100,
477  2060, 2055, 2055, 2055, 2067, 2085, 2000, 1900, 1835, 1770,
478  1720, 1700, 1695, 1695, 1695};
479  static const G4double sn_m_in[35] = {
480  1050, 1350, 1520, 1650, 1800, 1980, 2070, 2120, 2090, 2050,
481  1980, 1920, 1830, 1770, 1670, 1500, 1435, 1350, 1300, 1230,
482  1220, 1235, 1235, 1235, 1237, 1240, 1160, 1120, 1090, 1065,
483  1040, 1020, 1015, 1015, 1015};
484  static const G4double sn_p_t[35] = {
485  465, 800, 1200, 1760, 2170, 2480, 2730, 2885, 2970, 2980,
486  2970, 2890, 2840, 2790, 2620, 2450, 2335, 2205, 2080, 2020,
487  2010, 1990, 1990, 2015, 2030, 2045, 1980, 1890, 1835, 1770,
488  1720, 1700, 1695, 1695, 1695};
489  static const G4double sn_p_in[35] = {
490  315, 590, 880, 1220, 1460, 1580, 1700, 1770, 1810, 1810,
491  1800, 1730, 1680, 1630, 1530, 1400, 1335, 1270, 1210, 1180,
492  1190, 1190, 1190, 1205, 1210, 1210, 1150, 1115, 1090, 1065,
493  1040, 1020, 1015, 1015, 1015};
494  static const G4double w_m_t[35] = {
495  5200, 5115, 5025, 4975, 4900, 4850, 4780, 4725, 4600, 4490,
496  4355, 4255, 4125, 4040, 3830, 3580, 3330, 3110, 2955, 2860,
497  2852, 2845, 2885, 2900, 2915, 2940, 2800, 2660, 2570, 2490,
498  2460, 2425, 2420, 2420, 2420};
499  static const G4double w_m_in[35] = {
500  1450, 1850, 2100, 2350, 2550, 2700, 2825, 2900, 2850, 2750,
501  2630, 2525, 2400, 2300, 2200, 2070, 1880, 1770, 1715, 1680,
502  1680, 1680, 1685, 1690, 1700, 1720, 1635, 1560, 1530, 1460,
503  1440, 1410, 1410, 1410, 1410};
504  static const G4double w_p_t[35] = {
505  480, 900, 1500, 2350, 3020, 3420, 3650, 3775, 3875, 3830,
506  3750, 3700, 3630, 3550, 3550, 3290, 3070, 2890, 2840, 2730,
507  2725, 2720, 2770, 2805, 2828, 2865, 2770, 2640, 2570, 2490,
508  2460, 2425, 2420, 2420, 2420};
509  static const G4double w_p_in[35] = {
510  325, 680, 990, 1500, 1850, 2150, 2250, 2300, 2350, 2330,
511  2280, 2230, 2200, 2120, 2130, 1900, 1780, 1670, 1635, 1600,
512  1602, 1605, 1610, 1615, 1630, 1660, 1620, 1550, 1530, 1460,
513  1440, 1410, 1410, 1410, 1410};
514 
515  static const G4double pb_m_t[35] = {
516  5890, 5700, 5610, 5580, 5550, 5480, 5400, 5300, 5100, 4930,
517  4750, 4600, 4400, 4280, 4170, 3915, 3650, 3470, 3260, 3150,
518  3120, 3070, 3085, 3100, 3120, 3160, 3070, 2930, 2820, 2750,
519  2710, 2655, 2640, 2640, 2640};
520  static const G4double pb_m_in[35] = {
521  1575, 2025, 2300, 2575, 2850, 3000, 3115, 3180, 3080, 2940,
522  2800, 2670, 2550, 2450, 2370, 2220, 2110, 2000, 1920, 1880,
523  1850, 1800, 1805, 1810, 1820, 1840, 1800, 1720, 1640, 1620,
524  1570, 1530, 1530, 1530, 1530};
525  static const G4double pb_p_t[35] = {
526  515, 940, 1500, 2400, 3270, 3750, 4050, 4140, 4260, 4200,
527  4080, 3990, 3990, 3810, 3730, 3520, 3370, 3186, 3110, 3010,
528  2990, 2985, 3005, 3020, 3040, 3080, 3020, 2905, 2790, 2750,
529  2710, 2655, 2640, 2640, 2640};
530  static const G4double pb_p_in[35] = {
531  348, 707, 1040, 1650, 2100, 2400, 2580, 2640, 2650, 2520,
532  2410, 2300, 2250, 2190, 2130, 2000, 1930, 1870, 1830, 1790,
533  1770, 1765, 1775, 1780, 1790, 1800, 1775, 1710, 1620, 1620,
534  1570, 1530, 1530, 1530, 1530};
535  static const G4double u_m_t[35] = {
536  7080, 6830, 6650, 6530, 6400, 6280, 6100, 5840, 5660, 5520,
537  5330, 5160, 4990, 4810, 4630, 4323, 4130, 3870, 3700, 3550,
538  3490, 3465, 3467, 3475, 3495, 3515, 3440, 3360, 3150, 3040,
539  2985, 2955, 2940, 2940, 2940};
540  static const G4double u_m_in[35] = {
541  1740, 2220, 2500, 2820, 3080, 3300, 3420, 3500, 3420, 3330,
542  3200, 3060, 2940, 2850, 2710, 2470, 2380, 2250, 2160, 2080,
543  2040, 2045, 2047, 2050, 2055, 2060, 2010, 1980, 1830, 1780,
544  1735, 1710, 1700, 1700, 1700};
545  static const G4double u_p_t[35] = {
546  485, 960, 1580, 2700, 3550, 4050, 4320, 4420, 4620, 4660,
547  4580, 4470, 4350, 4295, 4187, 3938, 3755, 3573, 3450, 3342,
548  3310, 3295, 3310, 3330, 3375, 3405, 3350, 3338, 3135, 3040,
549  2985, 2955, 2940, 2940, 2940};
550  static const G4double u_p_in[35] = {
551  334, 720, 1020, 1560, 2100, 2300, 2550, 2700, 2880, 2880,
552  2760, 2660, 2550, 2510, 2430, 2270, 2130, 2060, 2000, 1970,
553  1950, 1950, 1960, 1960, 1970, 1980, 1950, 1978, 1830, 1780,
554  1735, 1710, 1700, 1700, 1700};
555 
556  AddDataSet("pi-",he_t, he_in, e1, 38);
557  AddDataSet("pi+",he_t, he_in, e1, 38);
558  AddDataSet("pi-",be_m_t, be_m_in, e1, 38);
559  AddDataSet("pi+",be_p_t, be_p_in, e1, 38);
560  AddDataSet("pi-",c_m_t, c_m_in, e2, 39);
561  AddDataSet("pi+",c_p_t, c_p_in, e2, 39);
562  AddDataSet("pi-",n_m_t, n_m_in, e2, 39);
563  AddDataSet("pi+",n_p_t, n_p_in, e2, 39);
564  AddDataSet("pi-",o_m_t, o_m_in, e3, 31);
565  AddDataSet("pi+",o_p_t, o_p_in, e3, 31);
566  AddDataSet("pi-",na_m_t, na_m_in, e3, 31);
567  AddDataSet("pi+",na_p_t, na_p_in, e3, 31);
568  AddDataSet("pi-",al_m_t, al_m_in, e3, 31);
569  AddDataSet("pi+",al_p_t, al_p_in, e3, 31);
570  AddDataSet("pi-",ca_m_t, ca_m_in, e3, 31);
571  AddDataSet("pi+",ca_p_t, ca_p_in, e3, 31);
572  AddDataSet("pi-",fe_m_t, fe_m_in, e4, 32);
573  AddDataSet("pi+",fe_p_t, fe_p_in, e4, 32);
574  AddDataSet("pi-",cu_m_t, cu_m_in, e4, 32);
575  AddDataSet("pi+",cu_p_t, cu_p_in, e4, 32);
576  AddDataSet("pi-",mo_m_t, mo_m_in, e5, 34);
577  AddDataSet("pi+",mo_p_t, mo_p_in, e5, 34);
578  AddDataSet("pi-",cd_m_t, cd_m_in, e5, 34);
579  AddDataSet("pi+",cd_p_t, cd_p_in, e5, 34);
580  AddDataSet("pi-",sn_m_t, sn_m_in, e6, 35);
581  AddDataSet("pi+",sn_p_t, sn_p_in, e6, 35);
582  AddDataSet("pi-",w_m_t, w_m_in, e6, 35);
583  AddDataSet("pi+",w_p_t, w_p_in, e6, 35);
584  AddDataSet("pi-",pb_m_t, pb_m_in, e6, 35);
585  AddDataSet("pi+",pb_p_t, pb_p_in, e6, 35);
586  AddDataSet("pi-",u_m_t, u_m_in, e6, 35);
587  AddDataSet("pi+",u_p_t, u_p_in, e6, 35);
588 }
589 
590 void G4UPiNuclearCrossSection::CrossSectionDescription(std::ostream& outFile) const
591 {
592  outFile << "G4UPiNuclearCrossSection calculates the total, elastic and\n"
593  << "inelastic cross sections for pion scattering from nuclei\n"
594  << "heavier than hydrogen. It is based on the Barashenkov\n"
595  << "parameterization and is valid for all incident energies.\n";
596 }
597