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