Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
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
47G4UPiNuclearCrossSection::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
65G4UPiNuclearCrossSection::~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
89G4bool
90G4UPiNuclearCrossSection::IsElementApplicable(const G4DynamicParticle*,
91 G4int Z, const G4Material*)
92{
93 return (1 < Z);
94}
95
96G4double
97G4UPiNuclearCrossSection::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
111G4double
112G4UPiNuclearCrossSection::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
140G4double 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
170void 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
193void 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
208void 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
599void 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
G4double
double G4double
Definition: G4Types.hh:64
G4int
int G4int
Definition: G4Types.hh:66
G4bool
bool G4bool
Definition: G4Types.hh:67
G4endl
#define G4endl
Definition: G4ios.hh:52
G4cout
G4DLLIMPORT std::ostream G4cout
G4DynamicParticle::GetDefinition
G4ParticleDefinition * GetDefinition() const
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
G4LPhysicsFreeVector::PutValues
void PutValues(size_t binNumber, G4double binValue, G4double dataValue)
G4NistManager::Instance
static G4NistManager * Instance()
Definition: G4NistManager.cc:68
G4NistManager::GetAtomicMassAmu
G4double GetAtomicMassAmu(const G4String &symb) const
Definition: G4NistManager.hh:343
G4PhysicsTable::push_back
void push_back(G4PhysicsVector *)
G4PhysicsTable::clearAndDestroy
void clearAndDestroy()
G4PionMinus::PionMinus
static G4PionMinus * PionMinus()
Definition: G4PionMinus.cc:98
G4PionPlus::PionPlus
static G4PionPlus * PionPlus()
Definition: G4PionPlus.cc:98
G4UPiNuclearCrossSection::GetInelasticCrossSection
G4double GetInelasticCrossSection(const G4DynamicParticle *aParticle, G4int Z, G4int A)
Definition: G4UPiNuclearCrossSection.cc:112
G4UPiNuclearCrossSection::BuildPhysicsTable
void BuildPhysicsTable(const G4ParticleDefinition &)
Definition: G4UPiNuclearCrossSection.cc:208
G4UPiNuclearCrossSection::GetElasticCrossSection
G4double GetElasticCrossSection(const G4DynamicParticle *aParticle, G4int Z, G4int A)
Definition: G4UPiNuclearCrossSection.cc:97
G4UPiNuclearCrossSection::IsElementApplicable
virtual G4bool IsElementApplicable(const G4DynamicParticle *aParticle, G4int Z, const G4Material *)
Definition: G4UPiNuclearCrossSection.cc:90
G4UPiNuclearCrossSection::CrossSectionDescription
virtual void CrossSectionDescription(std::ostream &) const
Definition: G4UPiNuclearCrossSection.cc:599
G4UPiNuclearCrossSection::DumpPhysicsTable
void DumpPhysicsTable(const G4ParticleDefinition &)
Definition: G4UPiNuclearCrossSection.cc:193