| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 // Author: Dennis Wright (SLAC)
27 // Date: 28 January 2013
28 //
29 // Description: implementation of numerically
30 // class for gamma p -> p pi0 rea
31 //
32 // 20130219 Inherit from templated base, move
33
34 #include "G4GamP2PPi0AngDst.hh"
35
36 namespace {
37 static const G4double eBins[15] =
38 { 0.145, 0.169, 0.240, 0.300, 0.360, 0.420
39 0.700, 1.050, 1.238, 1.400, 1.575, 1.825
40 2.900 };
41
42 static const G4double angleBins[19] =
43 { -1.000, -0.985, -0.940, -0.866, -0.766,
44 0.000, 0.174, 0.342, 0.500, 0.643,
45 1.000 };
46
47 static const G4double integralTable[15][19]
48 {0.0000, 0.0075, 0.0300, 0.0670, 0.1170, 0
49 0.5000, 0.5870, 0.6710, 0.7500, 0.8215, 0
50 1.0000},
51 {0.0000, 0.0123, 0.0482, 0.1041, 0.1739, 0
52 0.5543, 0.6214, 0.6862, 0.7500, 0.8124, 0
53 1.0000},
54 {0.0000, 0.0059, 0.0242, 0.0564, 0.1039, 0
55 0.5499, 0.6534, 0.7473, 0.8274, 0.8910, 0
56 1.0000},
57 {0.0000, 0.0039, 0.0164, 0.0399, 0.0769, 0
58 0.4913, 0.5980, 0.6981, 0.7866, 0.8598, 0
59 1.0000},
60 {0.0000, 0.0025, 0.0113, 0.0294, 0.0605, 0
61 0.4686, 0.5798, 0.6848, 0.7780, 0.8550, 0
62 1.0000},
63 {0.0000, 0.0015, 0.0074, 0.0212, 0.0476, 0
64 0.4485, 0.5627, 0.6714, 0.7684, 0.8486, 0
65 1.0000},
66 {0.0000, 0.0009, 0.0051, 0.0160, 0.0383, 0
67 0.4135, 0.5260, 0.6356, 0.7360, 0.8220, 0
68 1.0000},
69 {0.0000, 0.0022, 0.0105, 0.0290, 0.0630, 0
70 0.5261, 0.6454, 0.7529, 0.8421, 0.9092, 0
71 1.0000},
72 {0.0000, 0.0080, 0.0333, 0.0785, 0.1424, 0
73 0.5071, 0.5757, 0.6547, 0.7432, 0.8313, 0
74 1.0000},
75 {0.0000, 0.0064, 0.0288, 0.0717, 0.1333, 0
76 0.4566, 0.5307, 0.6114, 0.6937, 0.7726, 0
77 1.0000},
78 {0.0000, 0.0129, 0.0517, 0.1132, 0.1849, 0
79 0.4878, 0.5668, 0.6343, 0.6844, 0.7265, 0
80 1.0000},
81 {0.0000, 0.0027, 0.0197, 0.0638, 0.1277, 0
82 0.4892, 0.5796, 0.6246, 0.6388, 0.6663, 0
83 1.0000},
84 {0.0000, 0.0025, 0.0142, 0.0401, 0.0761, 0
85 0.4282, 0.5179, 0.5649, 0.5863, 0.6227, 0
86 1.0000},
87 {0.0000, 0.0010, 0.0036, 0.0062, 0.0083, 0
88 0.0844, 0.1307, 0.1957, 0.2864, 0.4073, 0
89 1.0000},
90 {0.0000, 0.0013, 0.0044, 0.0079, 0.0111, 0
91 0.0895, 0.1359, 0.2010, 0.2917, 0.4121, 0
92 1.0000}
93 };
94 }
95
96 // Constructor passes arrays to templated base
97
98 G4GamP2PPi0AngDst::G4GamP2PPi0AngDst(G4int ver
99 : G4NumIntTwoBodyAngDst<15,19>("G4GamP2PPi0A
100 integralTable, 1.5, verbose) {;}
101