/
radDamage.cc
219 lines (184 loc) · 83.6 KB
/
radDamage.cc
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
#include "radDamage.hh"
#include <iostream>
#include <cmath>
using namespace std;
/*E[MeV], theta[rad]*/
double radDamage::getNEIL(int partType,double energy, double theta){
int nDmg(-1);
if(abs(partType) == 11){
nDmg=3;
}else if(partType == 2112){
nDmg=0;
}else if(partType == 2212){
nDmg=1;
}else if(partType == 211 || partType == -211){
nDmg=2;
}else{
return -999;
}
double interpolatedValue = interpolate(xValNEIL[nDmg],yValNEIL[nDmg],energy);
double pi=acos(-1);
if( abs(theta/(pi/2) - 1) < 0.01 ) theta=pi/2*1.01;
interpolatedValue /= abs(cos(theta));
if(interpolatedValue < 0 || std::isnan(interpolatedValue) || std::isinf(interpolatedValue)){
cout<<__LINE__<<"\t"<<__PRETTY_FUNCTION__<<endl
<<"\tBad interpolated value: "<<interpolatedValue
<<"\t for nDmg "<<nDmg<<"\t and energy "<<energy<<"\t and theta"<<theta<<endl;
return 0;
}
return interpolatedValue;
}
/*E[MeV], theta[rad] => mrem*/
double radDamage::getMREM(int partType, double energy, double theta){
int nDmg(-1);
if(abs(partType) == 11){
nDmg=2;
}else if(partType == 2112){
nDmg=0;
}else if(partType == 22){
nDmg=1;
}else{
return -999;
}
//Tables are in MeV
double interpolatedValue = interpolate(xValMREM[nDmg],yValMREM[nDmg],energy);
double pi=acos(-1);
if( abs(theta/(pi/2) - 1) < 0.01 ) theta=pi/2*1.01;
if(nDmg==0)
interpolatedValue = 1*interpolatedValue/abs(cos(theta));
else if(nDmg==1 && interpolatedValue > 0)
interpolatedValue = 1/(3600*interpolatedValue*abs(cos(theta)));
else if(nDmg==2 && interpolatedValue > 0 )
interpolatedValue = 1/(3600*interpolatedValue*abs(cos(theta)));
if(interpolatedValue < 0 || std::isnan(interpolatedValue) || std::isinf(interpolatedValue)){
cout<<__LINE__<<"\t"<<__PRETTY_FUNCTION__<<endl
<<"\tBad interpolated value: "<<interpolatedValue
<<"\t for nDmg "<<nDmg<<"\t and energy "<<energy<<"\t and theta"<<theta<<endl;
return 0;
}
return interpolatedValue;
}
double radDamage::interpolate(vector<double> xV, vector<double> yV,double energy){
if( energy <= xV.front() ) return 0;
else if( energy >= xV.back() ) return yV.back();
int lowIndex = int(lower_bound(xV.begin(),xV.end(),energy) - xV.begin() ) - 1;
int highIndex = lowIndex + 1;
double lowVal = yV[lowIndex];
double highVal = yV[highIndex];
double lowX = xV[lowIndex];
double highX = xV[highIndex];
double returnVal = ( lowVal + (energy - lowX) * (highVal - lowVal) / (highX - lowX) );
return returnVal;
}
void radDamage::InitNEIL(){
double dmgNeutronX[] = {0.0000000001025 , 0.0000000001075 , 0.0000000001125 , 0.0000000001175 , 0.00000000012375 , 0.00000000013125 , 0.00000000013875 , 0.00000000014625 , 0.000000000155 , 0.000000000165 , 0.000000000175 , 0.000000000185 , 0.000000000195 , 0.000000000205 , 0.000000000215 , 0.000000000225 , 0.000000000235 , 0.0000000002475 , 0.0000000002625 , 0.000000000275 , 0.00000000029 , 0.00000000031 , 0.00000000033 , 0.00000000035 , 0.00000000037 , 0.00000000039 , 0.0000000004125 , 0.0000000004375 , 0.0000000004625 , 0.0000000004875 , 0.0000000005125 , 0.0000000005375 , 0.0000000005625 , 0.0000000005875 , 0.000000000615 , 0.000000000645 , 0.000000000675 , 0.000000000705 , 0.00000000074 , 0.00000000078 , 0.00000000082 , 0.00000000086 , 0.0000000009 , 0.00000000094 , 0.00000000098 , 0.000000001025 , 0.000000001075 , 0.000000001125 , 0.000000001175 , 0.0000000012375 , 0.0000000013125 , 0.0000000013875 , 0.0000000014625 , 0.00000000155 , 0.00000000165 , 0.00000000175 , 0.00000000185 , 0.00000000195 , 0.00000000205 , 0.00000000215 , 0.00000000225 , 0.00000000235 , 0.000000002475 , 0.000000002625 , 0.00000000275 , 0.0000000029 , 0.0000000031 , 0.0000000033 , 0.0000000035 , 0.0000000037 , 0.0000000039 , 0.000000004125 , 0.000000004375 , 0.000000004625 , 0.000000004875 , 0.000000005125 , 0.000000005375 , 0.000000005625 , 0.000000005875 , 0.00000000615 , 0.00000000645 , 0.00000000675 , 0.00000000705 , 0.000000007400001 , 0.0000000078 , 0.0000000082 , 0.0000000086 , 0.000000009 , 0.0000000094 , 0.0000000098 , 0.00000001025 , 0.00000001075 , 0.00000001125 , 0.00000001175 , 0.000000012375 , 0.000000013125 , 0.000000013875 , 0.000000014625 , 0.0000000155 , 0.0000000165 , 0.0000000175 , 0.0000000185 , 0.0000000195 , 0.0000000205 , 0.0000000215 , 0.0000000225 , 0.0000000235 , 0.00000002475 , 0.00000002625 , 0.0000000275 , 0.000000029 , 0.000000031 , 0.000000033 , 0.000000035 , 0.000000037 , 0.000000039 , 0.00000004125 , 0.00000004375 , 0.00000004625 , 0.00000004875 , 0.00000005125 , 0.00000005375 , 0.00000005625 , 0.00000005875 , 0.0000000615 , 0.0000000645 , 0.0000000675 , 0.0000000705 , 0.000000074 , 0.000000078 , 0.000000082 , 0.000000086 , 0.00000009 , 0.000000094 , 0.000000098 , 0.0000001025 , 0.0000001075 , 0.0000001125 , 0.0000001175 , 0.00000012375 , 0.00000013125 , 0.00000013875 , 0.00000014625 , 0.000000155 , 0.000000165 , 0.000000175 , 0.000000185 , 0.000000195 , 0.000000205 , 0.000000215 , 0.000000225 , 0.000000235 , 0.0000002475 , 0.0000002625 , 0.000000275 , 0.00000029 , 0.00000031 , 0.00000033 , 0.00000035 , 0.00000037 , 0.00000039 , 0.0000004125 , 0.0000004375 , 0.0000004625 , 0.0000004875 , 0.0000005125 , 0.0000005375 , 0.0000005625 , 0.0000005875 , 0.000000615 , 0.000000645 , 0.000000675 , 0.000000705 , 0.00000074 , 0.00000078 , 0.00000082 , 0.00000086 , 0.0000009 , 0.00000094 , 0.00000098 , 0.000001025 , 0.000001075 , 0.000001125 , 0.000001175 , 0.0000012375 , 0.0000013125 , 0.0000013875 , 0.0000014625 , 0.00000155 , 0.00000165 , 0.00000175 , 0.00000185 , 0.00000195 , 0.00000205 , 0.00000215 , 0.00000225 , 0.00000235 , 0.000002475 , 0.000002625 , 0.00000275 , 0.0000029 , 0.0000031 , 0.0000033 , 0.0000035 , 0.0000037 , 0.0000039 , 0.000004125 , 0.000004375 , 0.000004625 , 0.000004875 , 0.000005125 , 0.000005375 , 0.000005625 , 0.000005875 , 0.00000615 , 0.00000645 , 0.00000675 , 0.00000705 , 0.0000074 , 0.0000078 , 0.0000082 , 0.0000086 , 0.000009 , 0.0000094 , 0.0000098 , 0.00001025 , 0.00001075 , 0.00001125 , 0.00001175 , 0.000012375 , 0.000013125 , 0.000013875 , 0.000014625 , 0.0000155 , 0.0000165 , 0.0000175 , 0.0000185 , 0.0000195 , 0.0000205 , 0.0000215 , 0.0000225 , 0.0000235 , 0.00002475 , 0.00002625 , 0.0000275 , 0.000029 , 0.000031 , 0.000033 , 0.000035 , 0.000037 , 0.000039 , 0.00004125 , 0.00004375 , 0.00004625 , 0.00004875 , 0.00005125 , 0.00005375 , 0.00005625 , 0.00005875 , 0.0000615 , 0.0000645 , 0.0000675 , 0.0000705 , 0.000074 , 0.000078 , 0.000082 , 0.000086 , 0.00009 , 0.000094 , 0.000098 , 0.0001025 , 0.0001075 , 0.0001125 , 0.0001175 , 0.00012375 , 0.00013125 , 0.00013875 , 0.00014625 , 0.000155 , 0.000165 , 0.000175 , 0.000185 , 0.000195 , 0.000205 , 0.000215 , 0.000225 , 0.000235 , 0.0002475 , 0.0002625 , 0.000275 , 0.00029 , 0.00031 , 0.00033 , 0.00035 , 0.00037 , 0.00039 , 0.0004125 , 0.0004375 , 0.0004625 , 0.0004875 , 0.0005125 , 0.0005375 , 0.0005625 , 0.0005875 , 0.000615 , 0.000645 , 0.000675 , 0.000705 , 0.00074 , 0.00078 , 0.00082 , 0.00086 , 0.0009 , 0.00094 , 0.00098 , 0.001025 , 0.001075 , 0.001125 , 0.001175 , 0.0012375 , 0.0013125 , 0.0013875 , 0.0014625 , 0.00155 , 0.00165 , 0.00175 , 0.00185 , 0.00195 , 0.00205 , 0.00215 , 0.00225 , 0.00235 , 0.002475 , 0.002625 , 0.00275 , 0.0029 , 0.0031 , 0.0033 , 0.0035 , 0.0037 , 0.0039 , 0.004125 , 0.004375 , 0.004625 , 0.004875 , 0.005125 , 0.005375 , 0.005625 , 0.005875 , 0.00615 , 0.00645 , 0.00675 , 0.00705 , 0.0074 , 0.0078 , 0.0082 , 0.0086 , 0.009 , 0.0094 , 0.0098 , 0.01025 , 0.01075 , 0.01125 , 0.01175 , 0.012375 , 0.013125 , 0.013875 , 0.014625 , 0.0155 , 0.0165 , 0.0175 , 0.0185 , 0.0195 , 0.0205 , 0.0215 , 0.0225 , 0.0235 , 0.02475 , 0.02625 , 0.0275 , 0.029 , 0.031 , 0.033 , 0.035 , 0.037 , 0.039 , 0.04125 , 0.04375 , 0.04625 , 0.04875 , 0.05125 , 0.05375 , 0.05625 , 0.05875 , 0.06150001 , 0.0645 , 0.0675 , 0.0705 , 0.074 , 0.078 , 0.082 , 0.086 , 0.09 , 0.094 , 0.098 , 0.1025 , 0.1075 , 0.1125 , 0.1175 , 0.12375 , 0.13125 , 0.13875 , 0.14625 , 0.155 , 0.165 , 0.175 , 0.185 , 0.195 , 0.205 , 0.215 , 0.225 , 0.235 , 0.2475 , 0.2625 , 0.275 , 0.29 , 0.31 , 0.33 , 0.35 , 0.37 , 0.39 , 0.4125 , 0.4375 , 0.4625 , 0.4875 , 0.5125 , 0.5375 , 0.5625 , 0.5875 , 0.615 , 0.645 , 0.675 , 0.705 , 0.74 , 0.78 , 0.82 , 0.86 , 0.9 , 0.94 , 0.98 , 1.05 , 1.15 , 1.25 , 1.35 , 1.45 , 1.55 , 1.65 , 1.75 , 1.85 , 1.95 , 2.05 , 2.15 , 2.25 , 2.35 , 2.45 , 2.55 , 2.65 , 2.75 , 2.85 , 2.95 , 3.05 , 3.15 , 3.25 , 3.35 , 3.45 , 3.55 , 3.65 , 3.75 , 3.85 , 3.95 , 4.05 , 4.15 , 4.25 , 4.35 , 4.45 , 4.55 , 4.65 , 4.75 , 4.85 , 4.95 , 5.05 , 5.15 , 5.25 , 5.35 , 5.45 , 5.55 , 5.65 , 5.75 , 5.85 , 5.95 , 6.05 , 6.15 , 6.25 , 6.35 , 6.45 , 6.55 , 6.65 , 6.75 , 6.85 , 6.95 , 7.05 , 7.15 , 7.25 , 7.35 , 7.45 , 7.55 , 7.65 , 7.75 , 7.85 , 7.95 , 8.05 , 8.15 , 8.25 , 8.35 , 8.45 , 8.55 , 8.65 , 8.75 , 8.85 , 8.95 , 9.05 , 9.15 , 9.25 , 9.35 , 9.45 , 9.55 , 9.65 , 9.75 , 9.85 , 9.95 , 10.05 , 10.15 , 10.25 , 10.35 , 10.45 , 10.55 , 10.65 , 10.75 , 10.85 , 10.95 , 11.05 , 11.15 , 11.25 , 11.35 , 11.45 , 11.55 , 11.65 , 11.75 , 11.85 , 11.95 , 12.05 , 12.15 , 12.25 , 12.35 , 12.45 , 12.55 , 12.65 , 12.75 , 12.85 , 12.95 , 13.05 , 13.15 , 13.25 , 13.35 , 13.45 , 13.55 , 13.65 , 13.75 , 13.85 , 13.95 , 14.05 , 14.15 , 14.25 , 14.35 , 14.45 , 14.55 , 14.65 , 14.75 , 14.85 , 14.95 , 15.05 , 15.15 , 15.25 , 15.35 , 15.45 , 15.55 , 15.65 , 15.75 , 15.85 , 15.95 , 16.05 , 16.15 , 16.25 , 16.35 , 16.45 , 16.55 , 16.65 , 16.75 , 16.85 , 16.95 , 17.05 , 17.15 , 17.25 , 17.35 , 17.45 , 17.55 , 17.65 , 17.75 , 17.85 , 17.95 , 18.05 , 18.15 , 18.25 , 18.35 , 18.45 , 18.55 , 18.65 , 18.75 , 18.85 , 18.95 , 19.05 , 19.15 , 19.25 , 19.35 , 19.45 , 19.55 , 19.65 , 19.75 , 19.85 , 19.95 , 20 , 25 , 30 , 40 , 50 , 60 , 70 , 80 , 90 , 100 , 130 , 160 , 200 , 250 , 300 , 350 , 400 , 450 , 500 , 600 , 700 , 800 , 805 , 815 , 825 , 835 , 845 , 855 , 865 , 875 , 885 , 895 , 905 , 915 , 925 , 935 , 945 , 955 , 965 , 975 , 985 , 995 , 1005 , 1015 , 1025 , 1035 , 1045 , 1055 , 1065 , 1075 , 1085 , 1095 , 1105 , 1115 , 1125 , 1135 , 1145 , 1155 , 1165 , 1175 , 1185 , 1195 , 1205 , 1215 , 1225 , 1235 , 1245 , 1255 , 1265 , 1275 , 1285 , 1295 , 1305 , 1315 , 1325 , 1335 , 1345 , 1355 , 1365 , 1375 , 1385 , 1395 , 1405 , 1415 , 1425 , 1435 , 1445 , 1455 , 1465 , 1485 , 1495 , 1505 , 1515 , 1535 , 1545 , 1555 , 1565 , 1585 , 1595 , 1605 , 1615 , 1635 , 1645 , 1655 , 1665 , 1695 , 1705 , 1715 , 1735 , 1755 , 1765 , 1785 , 1795 , 1815 , 1835 , 1845 , 1855 , 1885 , 1895 , 1905 , 1915 , 1945 , 1955 , 1965 , 1985 , 2005 , 2015 , 2035 , 2045 , 2065 , 2085 , 2095 , 2105 , 2135 , 2145 , 2155 , 2165 , 2195 , 2205 , 2215 , 2235 , 2255 , 2265 , 2285 , 2295 , 2315 , 2335 , 2345 , 2355 , 2385 , 2395 , 2405 , 2415 , 2435 , 2445 , 2455 , 2465 , 2485 , 2495 , 2505 , 2515 , 2535 , 2545 , 2555 , 2565 , 2585 , 2595 , 2605 , 2615 , 2635 , 2645 , 2655 , 2665 , 2685 , 2695 , 2705 , 2715 , 2735 , 2745 , 2755 , 2765 , 2785 , 2795 , 2805 , 2815 , 2835 , 2845 , 2855 , 2865 , 2885 , 2895 , 2905 , 2915 , 2935 , 2945 , 2955 , 2965 , 2985 , 2995 , 3005 , 3015 , 3035 , 3045 , 3055 , 3065 , 3085 , 3095 , 3105 , 3115 , 3135 , 3145 , 3155 , 3165 , 3185 , 3195 , 3205 , 3215 , 3235 , 3245 , 3255 , 3265 , 3285 , 3295 , 3305 , 3315 , 3335 , 3345 , 3355 , 3365 , 3385 , 3395 , 3405 , 3415 , 3435 , 3445 , 3455 , 3465 , 3485 , 3495 , 3505 , 3515 , 3535 , 3545 , 3555 , 3565 , 3585 , 3595 , 3605 , 3615 , 3635 , 3645 , 3655 , 3665 , 3685 , 3695 , 3705 , 3715 , 3735 , 3745 , 3755 , 3765 , 3785 , 3795 , 3805 , 3815 , 3835 , 3845 , 3855 , 3865 , 3885 , 3895 , 3905 , 3915 , 3935 , 3945 , 3955 , 3965 , 3985 , 3995 , 4005 , 4015 , 4025 , 4035 , 4045 , 4055 , 4065 , 4075 , 4085 , 4095 , 4105 , 4115 , 4125 , 4135 , 4145 , 4155 , 4165 , 4175 , 4185 , 4195 , 4205 , 4215 , 4225 , 4235 , 4245 , 4255 , 4265 , 4275 , 4285 , 4295 , 4305 , 4315 , 4325 , 4335 , 4345 , 4355 , 4365 , 4375 , 4385 , 4395 , 4405 , 4415 , 4425 , 4435 , 4445 , 4455 , 4465 , 4475 , 4485 , 4495 , 4505 , 4515 , 4525 , 4535 , 4545 , 4555 , 4565 , 4575 , 4585 , 4595 , 4605 , 4615 , 4625 , 4635 , 4645 , 4655 , 4665 , 4675 , 4685 , 4695 , 4705 , 4715 , 4725 , 4735 , 4745 , 4755 , 4765 , 4775 , 4785 , 4795 , 4805 , 4815 , 4825 , 4835 , 4845 , 4855 , 4865 , 4875 , 4885 , 4895 , 4905 , 4915 , 4925 , 4935 , 4945 , 4955 , 4965 , 4975 , 4985 , 4995 , 5005 , 5015 , 5025 , 5035 , 5045 , 5055 , 5065 , 5075 , 5085 , 5095 , 5105 , 5115 , 5125 , 5135 , 5145 , 5155 , 5165 , 5175 , 5185 , 5195 , 5205 , 5215 , 5225 , 5235 , 5245 , 5255 , 5265 , 5275 , 5285 , 5295 , 5305 , 5315 , 5325 , 5335 , 5345 , 5355 , 5365 , 5375 , 5385 , 5395 , 5405 , 5415 , 5425 , 5435 , 5445 , 5455 , 5465 , 5475 , 5485 , 5495 , 5505 , 5515 , 5525 , 5535 , 5545 , 5555 , 5565 , 5575 , 5585 , 5595 , 5605 , 5615 , 5625 , 5635 , 5645 , 5655 , 5665 , 5675 , 5685 , 5695 , 5705 , 5715 , 5725 , 5735 , 5745 , 5755 , 5765 , 5775 , 5785 , 5795 , 5805 , 5815 , 5825 , 5835 , 5845 , 5855 , 5865 , 5875 , 5885 , 5895 , 5905 , 5915 , 5925 , 5935 , 5945 , 5955 , 5965 , 5975 , 5985 , 5995 , 6005 , 6015 , 6025 , 6035 , 6045 , 6055 , 6065 , 6075 , 6085 , 6095 , 6105 , 6115 , 6125 , 6135 , 6145 , 6155 , 6165 , 6175 , 6185 , 6195 , 6205 , 6215 , 6225 , 6235 , 6245 , 6255 , 6265 , 6275 , 6285 , 6295 , 6305 , 6315 , 6325 , 6335 , 6345 , 6355 , 6365 , 6375 , 6385 , 6395 , 6405 , 6415 , 6425 , 6435 , 6445 , 6455 , 6465 , 6475 , 6485 , 6495 , 6505 , 6515 , 6525 , 6535 , 6545 , 6555 , 6565 , 6575 , 6585 , 6595 , 6605 , 6615 , 6625 , 6635 , 6645 , 6655 , 6665 , 6675 , 6685 , 6695 , 6705 , 6715 , 6725 , 6735 , 6745 , 6755 , 6765 , 6775 , 6785 , 6795 , 6805 , 6815 , 6825 , 6835 , 6845 , 6855 , 6865 , 6875 , 6885 , 6895 , 6905 , 6915 , 6925 , 6935 , 6945 , 6955 , 6965 , 6975 , 6985 , 6995 , 7005 , 7015 , 7025 , 7035 , 7045 , 7055 , 7065 , 7075 , 7085 , 7095 , 7105 , 7115 , 7125 , 7135 , 7145 , 7155 , 7165 , 7185 , 7195 , 7205 , 7225 , 7235 , 7245 , 7255 , 7275 , 7285 , 7295 , 7315 , 7325 , 7335 , 7345 , 7365 , 7375 , 7385 , 7395 , 7415 , 7425 , 7435 , 7445 , 7465 , 7475 , 7485 , 7495 , 7515 , 7525 , 7535 , 7545 , 7565 , 7575 , 7585 , 7595 , 7615 , 7625 , 7635 , 7645 , 7665 , 7675 , 7685 , 7695 , 7715 , 7725 , 7735 , 7745 , 7765 , 7775 , 7785 , 7795 , 7815 , 7825 , 7835 , 7845 , 7865 , 7875 , 7885 , 7895 , 7915 , 7925 , 7935 , 7945 , 7965 , 7975 , 7985 , 7995 , 8015 , 8025 , 8035 , 8045 , 8065 , 8075 , 8085 , 8095 , 8115 , 8125 , 8135 , 8145 , 8165 , 8175 , 8185 , 8195 , 8215 , 8225 , 8235 , 8245 , 8265 , 8275 , 8285 , 8295 , 8315 , 8325 , 8335 , 8345 , 8365 , 8375 , 8385 , 8395 , 8415 , 8425 , 8435 , 8445 , 8465 , 8475 , 8485 , 8495 , 8515 , 8525 , 8535 , 8545 , 8565 , 8575 , 8585 , 8595 , 8615 , 8625 , 8635 , 8645 , 8665 , 8675 , 8685 , 8695 , 8715 , 8725 , 8735 , 8745 , 8765 , 8775 , 8785 , 8795 , 8815 , 8825 , 8835 , 8845 , 8865 , 8875 , 8885 , 8895 , 8915 , 8925 , 8935 , 8945 , 8965 , 8975 , 8995};
double dmgNeutronY[] = {0.0157512621578947 , 0.0153688431263158 , 0.0150285462842105 , 0.0147047370105263 , 0.0143243217263158 , 0.0139080129052632 , 0.0135286802105263 , 0.0131729914 , 0.0127967515473684 , 0.0124163857368421 , 0.0120576040842105 , 0.0117162806315789 , 0.0114202379473684 , 0.0111292886 , 0.0108712486526316 , 0.0106266001789474 , 0.0103932208315789 , 0.0101305860105263 , 0.00983410248421053 , 0.00960907055789474 , 0.00935597590526316 , 0.00904828568421053 , 0.00876959593684211 , 0.00851557343157895 , 0.00828972061052632 , 0.00807475530526316 , 0.00784674569473684 , 0.00762212181052632 , 0.00740983974736842 , 0.00721935081052632 , 0.00703882648421053 , 0.00687455121052632 , 0.00671839629473684 , 0.00657488923157895 , 0.00642530669473684 , 0.00627377464210526 , 0.00613365655789474 , 0.00600156346315789 , 0.00586411377894737 , 0.00570888346315789 , 0.00556681261052632 , 0.00543789627368421 , 0.00531176024210526 , 0.0052003202 , 0.00509114390526316 , 0.0049782726 , 0.00486066822105263 , 0.00475148697894737 , 0.00464878755789474 , 0.00452995098947368 , 0.00439795495789474 , 0.00427768217894737 , 0.00417096701052632 , 0.00404945709473684 , 0.00392496454736842 , 0.00381013552631579 , 0.00370610190526316 , 0.0036089532 , 0.00352023908421053 , 0.00343666113684211 , 0.00335997070736842 , 0.00328698381368421 , 0.00320340156947368 , 0.00311006530947368 , 0.00303850935368421 , 0.00296192223578947 , 0.00286374334210526 , 0.00277605648736842 , 0.00269463793157895 , 0.00262123876736842 , 0.00255236110631579 , 0.00248224760947368 , 0.00240946636315789 , 0.00234377203368421 , 0.00228261072842105 , 0.00222596005578947 , 0.00217384273263158 , 0.00212541774947368 , 0.00208162877684211 , 0.00203351835052632 , 0.00198478935578947 , 0.00194130439473684 , 0.00189812806736842 , 0.00185341156736842 , 0.0018045802368421 , 0.00176037899368421 , 0.00171843355473684 , 0.00168019100105263 , 0.00164359193157895 , 0.00160998044105263 , 0.00157389597789474 , 0.00153709093578947 , 0.00150286087684211 , 0.00147191681052632 , 0.00143326263684211 , 0.00139183245473684 , 0.00135348658947368 , 0.00131812368105263 , 0.00128019073684211 , 0.00124091799473684 , 0.00120442718421053 , 0.00117152572736842 , 0.00114099360631579 , 0.00111252186210526 , 0.00108650662210526 , 0.00106224988210526 , 0.00104035309157895 , 0.00101308903473684 , 0.00098412999368421 , 0.000961924738947368 , 0.000935915761052631 , 0.000905785214736842 , 0.000877309187368421 , 0.000852308511578947 , 0.000828499573684211 , 0.000807260841052632 , 0.000784695791578947 , 0.000761863082105263 , 0.000741004702105263 , 0.000721502903157895 , 0.000703625514736842 , 0.000687064397894737 , 0.000671623130526316 , 0.000657806598947369 , 0.000642776782105263 , 0.000627336334736842 , 0.000613643191578947 , 0.000600074087368421 , 0.00058607248 , 0.000570672667368421 , 0.000556732916842105 , 0.000543410094736842 , 0.000531299737894737 , 0.000519651992631579 , 0.000509002023157895 , 0.000497611581052631 , 0.000486396013684211 , 0.000475293138947368 , 0.000464580987368421 , 0.000453118812631579 , 0.000439661581052632 , 0.000427860183157895 , 0.00041659296 , 0.00040491510631579 , 0.00039236336 , 0.000380972423157895 , 0.000370517730526316 , 0.000360741497894737 , 0.000351797208421053 , 0.000343747877894737 , 0.000336171652736842 , 0.000328605628631579 , 0.000320370493052632 , 0.000311118761368421 , 0.000303933356947368 , 0.000296212853368421 , 0.000286487684210526 , 0.000277523227157895 , 0.000269453527052632 , 0.000261979663368421 , 0.000255112527052632 , 0.000248348442105263 , 0.000241028595684211 , 0.000234274809368421 , 0.000228383666947368 , 0.000222596597894737 , 0.000217600179894737 , 0.000212758768631579 , 0.000208173833684211 , 0.000203413164631579 , 0.000198468702 , 0.000194058941368421 , 0.000189812871789474 , 0.000185238235473684 , 0.000180509022105263 , 0.000176171351894737 , 0.000171842860736842 , 0.000167936743578947 , 0.000164461684736842 , 0.000160976952105263 , 0.00015752305031579 , 0.000154171080315789 , 0.000150809866210526 , 0.000147448189894737 , 0.000143356664842105 , 0.000139347496315789 , 0.000135574799157895 , 0.000131811873368421 , 0.000128070631894737 , 0.000124153688842105 , 0.000120401588947368 , 0.000117214429578947 , 0.000114089127052632 , 0.000111446935789474 , 0.000109021145894737 , 0.000106605376105263 , 0.000104189440947368 , 0.000101318180631579 , 0.0000985250649473684 , 0.0000962686094736842 , 0.0000936049987368421 , 0.0000905610156842105 , 0.0000877607762105263 , 0.0000851280118947368 , 0.0000828663413684211 , 0.000080647888631579 , 0.0000785959898947368 , 0.0000764402826315789 , 0.000074284493368421 , 0.0000721657362105263 , 0.0000704530647368421 , 0.0000688534803157895 , 0.0000672539284210526 , 0.0000656913266315789 , 0.0000642726908421053 , 0.0000627810109473684 , 0.0000612904826315789 , 0.0000599796693684211 , 0.0000585836541052632 , 0.0000570149161052632 , 0.0000557227418947368 , 0.0000544983147368421 , 0.0000532729664210526 , 0.0000520485717894737 , 0.0000508663917894737 , 0.0000498106408421053 , 0.0000486767127368421 , 0.0000475428666315789 , 0.0000464418082105263 , 0.0000453130951578947 , 0.0000439910621052632 , 0.0000427389308421053 , 0.0000416667281052631 , 0.0000404495392631579 , 0.0000392930401052632 , 0.000038208466 , 0.0000371229212631579 , 0.0000360650806315789 , 0.0000352180318947368 , 0.0000344161252631579 , 0.0000336132747052632 , 0.0000328361277578947 , 0.0000320394166631579 , 0.0000311039481157895 , 0.0000303462728736842 , 0.0000295896307473684 , 0.0000286007050105263 , 0.0000277834988526316 , 0.0000270155777052632 , 0.0000262476598105263 , 0.0000255002897052632 , 0.0000248331102631579 , 0.0000241237526526316 , 0.0000234185672210526 , 0.0000228203124315789 , 0.0000222692764210526 , 0.0000217182996421053 , 0.0000212320590736842 , 0.0000207889210842105 , 0.0000203027364947368 , 0.0000198503889684211 , 0.0000194433124526316 , 0.0000190362326842105 , 0.0000185613368736842 , 0.0000180329741157895 , 0.0000176114267894737 , 0.0000172095620736842 , 0.0000168086152631579 , 0.0000164210664315789 , 0.0000161003167894737 , 0.0000157498000631579 , 0.0000153612181157895 , 0.0000150168903263158 , 0.0000147033529789474 , 0.0000143188475157895 , 0.0000139333746736842 , 0.0000135726250526316 , 0.0000132128493157895 , 0.0000171374971368421 , 0.0000679163179789474 , 0.000131544770873684 , 0.000193809931526316 , 0.000225048217894737 , 0.000243256728421053 , 0.000261423653684211 , 0.000275397372631579 , 0.000286685087368421 , 0.000300740664210526 , 0.000317659829473684 , 0.000331809987368421 , 0.000348631135789474 , 0.000371201658947368 , 0.000393728606315789 , 0.000416253913684211 , 0.000438779715789474 , 0.000461261447368421 , 0.000486487528421053 , 0.000514566841052632 , 0.000542582175789474 , 0.000570596195789474 , 0.000598545912631579 , 0.000626431156842105 , 0.000654316895789474 , 0.000682212341052631 , 0.000712843249473684 , 0.000746188243157895 , 0.00077947864 , 0.000812760645263158 , 0.00085148858 , 0.000895711136842105 , 0.000939857224210526 , 0.00098400904631579 , 0.00102802172736842 , 0.00107195331684211 , 0.00111588014736842 , 0.00116509121052632 , 0.00121978902105263 , 0.00127438482105263 , 0.00132888850526316 , 0.001396637 , 0.00147795910526316 , 0.00155898351578947 , 0.00163999477894737 , 0.00173423262105263 , 0.00184157034736842 , 0.00194869896842105 , 0.00205575641052632 , 0.00216231407368421 , 0.00226922807368421 , 0.00238217517894737 , 0.00440461631578947 , 0.00258932454736842 , 0.00271699857894737 , 0.00287350269473684 , 0.00300416553684211 , 0.00315874455789474 , 0.00336605235789474 , 0.0035727464 , 0.00377771304210526 , 0.00398245275789474 , 0.00418621358947368 , 0.00441340610526316 , 0.00466609626315789 , 0.00491677512631579 , 0.00549160215789474 , 0.00541858412631579 , 0.00566222795789474 , 0.00590743705263158 , 0.00615081413684211 , 0.00641816272631579 , 0.00671005333684211 , 0.00699941355789474 , 0.00728532933684211 , 0.00761818767368421 , 0.00799787373684211 , 0.00837209887368421 , 0.00874529174736842 , 0.00911482617894737 , 0.0094766198631579 , 0.009837989 , 0.0102436503368421 , 0.0106949641263158 , 0.0111459523157895 , 0.0115852735789474 , 0.0121281263157895 , 0.0127679776842105 , 0.0133970077894737 , 0.0140284158947368 , 0.0155678135789474 , 0.015544694 , 0.0163355285263158 , 0.0171014743157895 , 0.0178215776842105 , 0.0185411056842105 , 0.0192593831578947 , 0.0199780206315789 , 0.0206895837894737 , 0.0215001503157895 , 0.0224461548421053 , 0.0231658078947368 , 0.0240007248421053 , 0.0250928786315789 , 0.0259890618947368 , 0.0267513493684211 , 0.0274325266315789 , 0.032468826 , 0.0279203647368421 , 0.0273873185263158 , 0.0259769833684211 , 0.0228832506315789 , 0.0162605327368421 , 0.0148522341052632 , 0.507765693263158 , 0.114514184210526 , 0.0779325315789474 , 0.067063391368421 , 0.0645463022105263 , 0.0598783783157895 , 0.0579702190526316 , 0.0563604932631579 , 0.0549398236842105 , 0.0537357332631579 , 0.0518848330526316 , 0.0499239576842105 , 0.0479017951578947 , 0.0452645898947369 , 0.0416229438947368 , 0.0376386267368421 , 0.0329657387368421 , 0.0264540102105263 , 0.0182649629473684 , 0.0112142661894737 , 0.0108952630105263 , 0.0484812250526316 , 0.217306914315789 , 0.69872448368421 , 1.19578253473684 , 1.16173878105263 , 0.971039121052632 , 0.828813665263158 , 0.733421109473684 , 0.670863415789474 , 0.617226047368421 , 0.576815617894737 , 0.554727975789474 , 0.539327132631579 , 0.52369232 , 0.519757351578947 , 0.509721517894737 , 0.503617456842105 , 0.531991170526316 , 0.540275653684211 , 0.545067665263158 , 0.558118343157895 , 0.574040401052632 , 0.598675961052632 , 0.755428550526316 , 1.28023130105263 , 0.599748601052632 , 0.542504049473684 , 0.553577689473684 , 0.578580974736842 , 0.601908811578948 , 0.672385488421053 , 0.920888551578947 , 1.45850912947368 , 0.831776736842105 , 0.943442582105263 , 1.17162367368421 , 1.17404088 , 0.802001543157895 , 0.657767798947369 , 0.967958267368421 , 0.940954031578947 , 1.07911352 , 1.12778311789474 , 1.76573729578947 , 0.836553629473684 , 1.41107035789474 , 1.39300664210526 , 1.02222715263158 , 1.15917350526316 , 1.12620987368421 , 1.10609297263158 , 1.26720092631579 , 1.38368153684211 , 1.23797910526316 , 1.15309143157895 , 1.43721563157895 , 1.06086187263158 , 1.27768070526316 , 1.42735447368421 , 1.28107795789474 , 1.21718761052632 , 1.25852553684211 , 1.21288823157895 , 0.750575436842105 , 1.23476252631579 , 1.18674070526316 , 1.45221824210526 , 1.43071410526316 , 1.14650597894737 , 1.78209818947368 , 1.44193953684211 , 1.50880461052632 , 1.49839917894737 , 1.672525 , 2.00336235789474 , 1.71740642105263 , 1.57875108421053 , 1.61158341052632 , 1.83127870526316 , 1.60486684210526 , 1.33566231578947 , 1.27684972631579 , 1.56710534736842 , 1.61574035789474 , 1.92938202105263 , 1.7942914 , 1.46796386315789 , 1.66561352631579 , 1.36398327368421 , 1.90003656842105 , 1.64147604210526 , 1.51442834736842 , 1.31441954736842 , 1.61103783157895 , 1.7782958 , 1.57655067368421 , 1.521453 , 1.77042192631579 , 1.46983914736842 , 1.80308335789474 , 1.81792343157895 , 1.79343238947368 , 1.7573862 , 1.79361070526316 , 1.8661588 , 1.86090195789474 , 1.99248669473684 , 1.8064966631579 , 1.66279608421053 , 1.74282183157895 , 1.76833778947368 , 1.78170988421053 , 1.77103870526316 , 1.76906190526316 , 1.58906486315789 , 1.68466082105263 , 1.87059383157895 , 1.87624228421053 , 1.68638874736842 , 1.61880791578947 , 1.75085283157895 , 1.82364149473684 , 1.76694875789474 , 1.66855772631579 , 1.70720475789474 , 1.75447877894737 , 1.7796052 , 1.76028574736842 , 1.78194009473684 , 1.7644398 , 1.77225547368421 , 1.70878303157895 , 1.69626821052632 , 1.67419409473684 , 1.7255118 , 1.76545288421053 , 1.76466037894737 , 1.73597844210526 , 1.68005977894737 , 1.69176962105263 , 1.73501585263158 , 1.76539863157895 , 1.77789982105263 , 1.76704234736842 , 1.75104012631579 , 1.75406268421053 , 1.76728756842105 , 1.78440068421053 , 1.81245296842105 , 1.81556608421053 , 1.82245722105263 , 1.79079383157895 , 1.75741704210526 , 1.83698125263158 , 1.85646310526316 , 1.84974588421053 , 1.85311922105263 , 1.86177897894737 , 1.82794374736842 , 1.79346717894737 , 1.80612885263158 , 1.83920676842105 , 1.80241610526316 , 1.77951367368421 , 1.80427749473684 , 1.87047569473684 , 1.80832013684211 , 1.78691104210526 , 1.81745627368421 , 1.83032054736842 , 1.84308125263158 , 1.84468901052632 , 1.79054329473684 , 1.76528918947368 , 1.77382794736842 , 1.80761515789474 , 1.81332495789474 , 1.77325673684211 , 1.80104069473684 , 1.84240529473684 , 1.84073691578947 , 1.82802410526316 , 1.80794287368421 , 1.80574167368421 , 1.8165948 , 1.85797098947368 , 1.89144098947368 , 1.89935962105263 , 1.89939144210526 , 1.89267822105263 , 1.87956034736842 , 1.86928132631579 , 1.91770804210526 , 1.97013645263158 , 1.95208435789474 , 1.92920005263158 , 1.93450167368421 , 1.94002696842105 , 1.93682362105263 , 1.93414658947368 , 1.93505345263158 , 1.93228105263158 , 1.91723203157895 , 1.90592123157895 , 1.91364171578947 , 1.92579314736842 , 1.95119133684211 , 1.97437574736842 , 1.96538073684211 , 1.94994502105263 , 1.93757338947368 , 1.9263564 , 1.94533028421053 , 1.98586841052632 , 2.00810714736842 , 1.99900036842105 , 1.98861403157895 , 1.97444967368421 , 1.96007648421053 , 1.95484478947368 , 1.95298971578947 , 1.95271033684211 , 1.95494907368421 , 1.95761381052632 , 1.96312390526316 , 1.9691570631579 , 1.97529223157895 , 2.071 , 2.049 , 2.041 , 2.012 , 1.805 , 1.644 , 1.499 , 1.378 , 1.264 , 1.168 , 0.974 , 0.865 , 0.791 , 0.733 , 0.696 , 0.693 , 0.685 , 0.685 , 0.681 , 0.679 , 0.677 , 0.672 , 0.600399052631579 , 0.598006210526316 , 0.595860210526316 , 0.594162842105263 , 0.593172315789474 , 0.592183894736842 , 0.591197473684211 , 0.590213263157895 , 0.589231157894737 , 0.588251157894737 , 0.587273263157895 , 0.586297473684211 , 0.585401157894737 , 0.584521894736842 , 0.583644210526316 , 0.582768 , 0.581893368421053 , 0.581020105263158 , 0.580148421052632 , 0.579278105263158 , 0.578381894736842 , 0.577317684210526 , 0.576255684210526 , 0.575195789473684 , 0.574138105263158 , 0.573421157894737 , 0.572815368421053 , 0.572210526315789 , 0.571606736842105 , 0.571090631578947 , 0.570754842105263 , 0.570420315789474 , 0.570086736842105 , 0.569754631578947 , 0.569423473684211 , 0.569093578947368 , 0.568764947368421 , 0.568437368421053 , 0.568107684210526 , 0.567774105263158 , 0.567441263157895 , 0.567109263157895 , 0.566777894736842 , 0.566447368421053 , 0.566117473684211 , 0.565788421052632 , 0.565460105263158 , 0.56497947368421 , 0.564266631578947 , 0.56355452631579 , 0.562843368421053 , 0.562132947368421 , 0.561423473684211 , 0.560714736842105 , 0.560006947368421 , 0.559299789473684 , 0.558698210526316 , 0.558098842105263 , 0.557500105263158 , 0.556902 , 0.55630452631579 , 0.555707578947368 , 0.555111263157895 , 0.554515578947368 , 0.553920526315789 , 0.553326 , 0.552732105263158 , 0.551546210526316 , 0.550954105263158 , 0.550362736842105 , 0.549771789473684 , 0.548591894736842 , 0.548002842105263 , 0.547414421052632 , 0.546826526315789 , 0.545652631578947 , 0.545066631578947 , 0.544481052631579 , 0.543896210526316 , 0.542728315789474 , 0.542145157894737 , 0.541562631578947 , 0.540980842105263 , 0.539238736842105 , 0.538659052631579 , 0.538080210526316 , 0.536924105263158 , 0.535794736842105 , 0.535359894736842 , 0.534491157894737 , 0.534057368421053 , 0.533190631578947 , 0.532325157894737 , 0.531893052631579 , 0.531461157894737 , 0.530167789473684 , 0.529737263157895 , 0.529307157894737 , 0.528877368421053 , 0.527590105263158 , 0.527161684210526 , 0.526733578947368 , 0.525878421052632 , 0.525024631578947 , 0.524598210526316 , 0.523746421052632 , 0.523320947368421 , 0.522471157894737 , 0.521622736842105 , 0.521198947368421 , 0.520775578947368 , 0.51950747368421 , 0.519085368421053 , 0.518663684210526 , 0.518242315789474 , 0.516980210526316 , 0.516560105263158 , 0.516140421052632 , 0.515301894736842 , 0.514464842105263 , 0.514046736842105 , 0.513725894736842 , 0.513616842105263 , 0.513398842105263 , 0.513180842105263 , 0.513071789473684 , 0.512962842105263 , 0.512636315789474 , 0.512527473684211 , 0.512418631578947 , 0.51231 , 0.512092526315789 , 0.511983894736842 , 0.511875157894737 , 0.511766631578947 , 0.51154947368421 , 0.511440842105263 , 0.511332421052632 , 0.511223894736842 , 0.511007052631579 , 0.510898631578947 , 0.510790210526316 , 0.510681894736842 , 0.510465368421053 , 0.510357052631579 , 0.510248736842105 , 0.510140421052632 , 0.509924210526316 , 0.509816105263158 , 0.509708 , 0.509599894736842 , 0.509383789473684 , 0.509275789473684 , 0.509167789473684 , 0.509059894736842 , 0.508844105263158 , 0.508736210526316 , 0.508628421052631 , 0.508520631578947 , 0.508305052631579 , 0.508197368421053 , 0.508089578947368 , 0.507982 , 0.507766736842105 , 0.507659052631579 , 0.507551473684211 , 0.507444 , 0.507228947368421 , 0.507121578947369 , 0.507014105263158 , 0.506906736842105 , 0.506692 , 0.506584736842105 , 0.506477473684211 , 0.506370105263158 , 0.506155684210526 , 0.506048526315789 , 0.505941263157895 , 0.505834210526316 , 0.50562 , 0.505512947368421 , 0.505405894736842 , 0.505298947368421 , 0.505084947368421 , 0.504978105263158 , 0.504871157894737 , 0.504764315789474 , 0.504550736842105 , 0.504443894736842 , 0.504337052631579 , 0.504230421052632 , 0.504016947368421 , 0.503910421052632 , 0.503803684210526 , 0.503697052631579 , 0.503484 , 0.503377473684211 , 0.503270947368421 , 0.503164631578947 , 0.502951684210526 , 0.502845263157895 , 0.502739052631579 , 0.502632631578947 , 0.502420105263158 , 0.502313789473684 , 0.502207578947368 , 0.502101473684211 , 0.501889052631579 , 0.501783052631579 , 0.501676947368421 , 0.501570842105263 , 0.501358631578947 , 0.501252736842105 , 0.501146842105263 , 0.501040842105263 , 0.500829157894737 , 0.500723157894737 , 0.500617368421053 , 0.500511578947368 , 0.500300105263158 , 0.500194315789474 , 0.500088736842105 , 0.499983052631579 , 0.499771789473684 , 0.499666210526316 , 0.499560631578947 , 0.499455052631579 , 0.499244105263158 , 0.499138631578947 , 0.499033263157895 , 0.498927789473684 , 0.498717157894737 , 0.498611789473684 , 0.498506526315789 , 0.498401157894737 , 0.498190736842105 , 0.498085578947368 , 0.497980421052632 , 0.497875263157895 , 0.497665052631579 , 0.49756 , 0.497455052631579 , 0.49735 , 0.49714 , 0.497035052631579 , 0.496930210526316 , 0.496825263157895 , 0.496615684210526 , 0.496510947368421 , 0.496406105263158 , 0.496301368421053 , 0.496091894736842 , 0.495987263157895 , 0.495882631578947 , 0.495778 , 0.495568947368421 , 0.495464421052632 , 0.495359789473684 , 0.495255473684211 , 0.495150947368421 , 0.49504652631579 , 0.494942105263158 , 0.494837684210526 , 0.494733368421053 , 0.494628947368421 , 0.494524631578947 , 0.494420526315789 , 0.494316210526316 , 0.494212 , 0.494107789473684 , 0.494003684210526 , 0.493899473684211 , 0.493795368421053 , 0.493691263157895 , 0.493587263157895 , 0.493483157894737 , 0.493379157894737 , 0.493275157894737 , 0.493171263157895 , 0.493067263157895 , 0.492963368421053 , 0.492859473684211 , 0.492755684210526 , 0.492651894736842 , 0.492548105263158 , 0.492444315789474 , 0.492340526315789 , 0.492236842105263 , 0.492133052631579 , 0.492029473684211 , 0.491925789473684 , 0.491822210526316 , 0.491718631578947 , 0.491615052631579 , 0.49151147368421 , 0.491408 , 0.491304526315789 , 0.491201052631579 , 0.491097578947368 , 0.490994210526316 , 0.490890842105263 , 0.490787368421053 , 0.490684105263158 , 0.490580842105263 , 0.490477473684211 , 0.490374315789474 , 0.490271052631579 , 0.490167789473684 , 0.490064736842105 , 0.48996147368421 , 0.489858421052632 , 0.489755263157895 , 0.489652210526316 , 0.489549157894737 , 0.489446105263158 , 0.489343052631579 , 0.489240210526316 , 0.489137263157895 , 0.489034315789474 , 0.488931368421053 , 0.488828526315789 , 0.488725684210526 , 0.488622947368421 , 0.488520105263158 , 0.488417263157895 , 0.488314631578947 , 0.488211894736842 , 0.488109263157895 , 0.48800652631579 , 0.487903894736842 , 0.487801368421053 , 0.487698736842105 , 0.487596105263158 , 0.487493578947368 , 0.487391157894737 , 0.487288631578947 , 0.487186210526316 , 0.487083684210526 , 0.486981263157895 , 0.486878947368421 , 0.486776631578947 , 0.486674315789474 , 0.486572 , 0.486469684210526 , 0.486367473684211 , 0.486265263157895 , 0.486162947368421 , 0.486060947368421 , 0.485958631578947 , 0.485856631578947 , 0.485754526315789 , 0.485652421052632 , 0.485550421052632 , 0.485448421052632 , 0.485346421052632 , 0.485244421052632 , 0.485142526315789 , 0.485040631578947 , 0.484938736842105 , 0.484836842105263 , 0.484734947368421 , 0.484633052631579 , 0.484531368421053 , 0.484429578947368 , 0.484327894736842 , 0.484226105263158 , 0.484124526315789 , 0.484022736842105 , 0.483921263157895 , 0.483819578947368 , 0.483718 , 0.483616421052632 , 0.483514842105263 , 0.483413368421053 , 0.483311894736842 , 0.483210421052632 , 0.483109052631579 , 0.483007578947368 , 0.482906210526316 , 0.482804842105263 , 0.482703473684211 , 0.482602210526316 , 0.482500947368421 , 0.482399684210526 , 0.482298526315789 , 0.482197368421053 , 0.482096105263158 , 0.481994947368421 , 0.481893789473684 , 0.481792631578947 , 0.481691578947368 , 0.481590631578947 , 0.481489578947368 , 0.481388631578947 , 0.481287578947368 , 0.481186631578947 , 0.481085684210526 , 0.480984842105263 , 0.480884 , 0.480783052631579 , 0.480682315789474 , 0.48058147368421 , 0.480480736842105 , 0.48038 , 0.480279263157895 , 0.480178526315789 , 0.480077894736842 , 0.479977157894737 , 0.479876631578947 , 0.479776 , 0.479675473684211 , 0.479574842105263 , 0.479474421052632 , 0.479373894736842 , 0.479273368421053 , 0.479172947368421 , 0.479072631578947 , 0.478972210526316 , 0.478871789473684 , 0.478771473684211 , 0.478671157894737 , 0.478570842105263 , 0.478470631578947 , 0.478370421052632 , 0.478270210526316 , 0.47817 , 0.478069789473684 , 0.477969684210526 , 0.47786947368421 , 0.477769473684211 , 0.477669368421053 , 0.477569368421053 , 0.477469263157895 , 0.477369368421053 , 0.477269368421053 , 0.477169368421053 , 0.47706947368421 , 0.476969578947368 , 0.476869789473684 , 0.476769894736842 , 0.476670105263158 , 0.476570315789474 , 0.476470631578947 , 0.476370842105263 , 0.476271052631579 , 0.476171473684211 , 0.476071789473684 , 0.475972105263158 , 0.47587252631579 , 0.475772947368421 , 0.475673368421053 , 0.475573789473684 , 0.475474315789474 , 0.475374842105263 , 0.475275368421053 , 0.475175894736842 , 0.475076526315789 , 0.474977052631579 , 0.474877789473684 , 0.474778421052632 , 0.474679157894737 , 0.474579789473684 , 0.474480631578947 , 0.474381263157895 , 0.474282105263158 , 0.474182947368421 , 0.474083789473684 , 0.473984631578947 , 0.473885578947368 , 0.473786421052632 , 0.473687368421053 , 0.473588315789474 , 0.473489263157895 , 0.473390315789474 , 0.473291368421053 , 0.473192421052632 , 0.473093473684211 , 0.472994631578947 , 0.472895684210526 , 0.472796842105263 , 0.472698105263158 , 0.472599263157895 , 0.47250052631579 , 0.472401789473684 , 0.472303052631579 , 0.472204421052632 , 0.472105684210526 , 0.472007157894737 , 0.471908526315789 , 0.471809894736842 , 0.471711263157895 , 0.471612736842105 , 0.471514210526316 , 0.471415684210526 , 0.471317263157895 , 0.471218842105263 , 0.471120421052632 , 0.471022 , 0.470923578947368 , 0.470825368421053 , 0.470727052631579 , 0.470628736842105 , 0.470530421052632 , 0.470432210526316 , 0.470333894736842 , 0.470235789473684 , 0.470137578947368 , 0.470039368421053 , 0.469941263157895 , 0.469843263157895 , 0.469745157894737 , 0.469647157894737 , 0.469549052631579 , 0.469451157894737 , 0.469353157894737 , 0.469255052631579 , 0.469157157894737 , 0.469059263157895 , 0.468961368421053 , 0.468863578947368 , 0.468765789473684 , 0.468668 , 0.468570210526316 , 0.468472421052632 , 0.468374736842105 , 0.468277052631579 , 0.468179263157895 , 0.468081684210526 , 0.467984 , 0.467886421052632 , 0.467788736842105 , 0.467691263157895 , 0.467593789473684 , 0.467496315789474 , 0.467398842105263 , 0.467301368421053 , 0.467203894736842 , 0.46710652631579 , 0.467009052631579 , 0.466911684210526 , 0.466814421052632 , 0.466717157894737 , 0.466619789473684 , 0.466522526315789 , 0.466425263157895 , 0.466328105263158 , 0.466230842105263 , 0.466133684210526 , 0.466036631578947 , 0.46593947368421 , 0.465842421052632 , 0.465745263157895 , 0.465648315789474 , 0.465551157894737 , 0.465454210526316 , 0.465357368421053 , 0.465260315789474 , 0.465163368421053 , 0.465066526315789 , 0.464969578947368 , 0.464872842105263 , 0.464776 , 0.464679157894737 , 0.464582421052632 , 0.464485578947368 , 0.464388947368421 , 0.464292210526316 , 0.464195473684211 , 0.464098842105263 , 0.464002315789474 , 0.463905578947368 , 0.463809052631579 , 0.46371252631579 , 0.463616 , 0.463423052631579 , 0.463326526315789 , 0.463230105263158 , 0.463037263157895 , 0.462940947368421 , 0.462844631578947 , 0.462748315789474 , 0.462555684210526 , 0.462459473684211 , 0.462363368421053 , 0.462170842105263 , 0.462074842105263 , 0.461978631578947 , 0.461882631578947 , 0.461690421052632 , 0.461594421052632 , 0.461498421052632 , 0.461402421052632 , 0.461210631578947 , 0.461114736842105 , 0.461018842105263 , 0.460922947368421 , 0.460731263157895 , 0.460635578947368 , 0.460539684210526 , 0.460444 , 0.460252631578947 , 0.460156947368421 , 0.460061263157895 , 0.459965684210526 , 0.45977452631579 , 0.459678947368421 , 0.45958347368421 , 0.459488 , 0.459297052631579 , 0.459201684210526 , 0.459106315789474 , 0.459010947368421 , 0.458820210526316 , 0.458724842105263 , 0.458629578947368 , 0.458534421052632 , 0.458343894736842 , 0.458248842105263 , 0.458153578947368 , 0.45805852631579 , 0.457868315789474 , 0.457773157894737 , 0.457678210526316 , 0.457583157894737 , 0.457393263157895 , 0.457298315789474 , 0.457203263157895 , 0.457108526315789 , 0.456918631578947 , 0.456823894736842 , 0.456729157894737 , 0.456634315789474 , 0.456444947368421 , 0.456350210526316 , 0.45625547368421 , 0.456160842105263 , 0.455971578947368 , 0.455877052631579 , 0.455782421052632 , 0.455687894736842 , 0.455498842105263 , 0.455404421052632 , 0.45531 , 0.455215578947368 , 0.455026842105263 , 0.454932421052632 , 0.454838105263158 , 0.454743789473684 , 0.454555263157895 , 0.454461052631579 , 0.454366947368421 , 0.454272736842105 , 0.454084421052632 , 0.453990315789474 , 0.453896210526316 , 0.453802105263158 , 0.453614105263158 , 0.453520105263158 , 0.453426105263158 , 0.453332210526316 , 0.453144421052632 , 0.453050526315789 , 0.452956631578947 , 0.452862736842105 , 0.452675157894737 , 0.45258147368421 , 0.452487684210526 , 0.452394 , 0.452206736842105 , 0.452113052631579 , 0.452019368421053 , 0.451925789473684 , 0.451738736842105 , 0.451645157894737 , 0.451551684210526 , 0.451458210526316 , 0.45127147368421 , 0.451178 , 0.451084631578947 , 0.450991157894737 , 0.450804631578947 , 0.450711263157895 , 0.450618105263158 , 0.450524842105263 , 0.450338421052632 , 0.450245157894737 , 0.450152105263158 , 0.450059052631579 , 0.449872736842105 , 0.449779789473684 , 0.449686736842105 , 0.449593684210526 , 0.449407789473684 , 0.449314842105263 , 0.449222 , 0.449129157894737 , 0.44894347368421 , 0.44885052631579 , 0.448757789473684 , 0.448665052631579 , 0.448479473684211 , 0.448386842105263 , 0.448294105263158 , 0.448201473684211 , 0.448016315789474 , 0.447923684210526 , 0.447831052631579 , 0.44773852631579 , 0.447553578947368 , 0.447461052631579 , 0.447368736842105 , 0.447276210526316 , 0.447091473684211 , 0.446999052631579 , 0.446906736842105 , 0.446814526315789 , 0.44663 , 0.446537684210526 , 0.446353263157895};
double dmgProtonX[]={1.000E-03, 2.000E-03, 3.000E-03, 5.000E-03, 7.000E-03, 1.000E-02, 2.000E-02, 3.000E-02, 5.000E-02, 7.000E-02, 1.000E-01, 2.000E-01, 3.000E-01, 5.000E-01, 7.000E-01, 1.000E+00, 2.000E+00, 3.000E+00, 5.000E+00, 7.000E+00, 1.000E+01, 1.500E+01, 2.000E+01, 2.500E+01, 3.000E+01, 3.500E+01, 4.500E+01, 5.500E+01, 5.000E+01, 6.500E+01, 7.000E+01, 7.500E+01, 8.500E+01, 9.500E+01, 1.000E+02, 1.050E+02, 1.150E+02, 1.250E+02, 1.350E+02, 1.450E+02, 1.550E+02, 1.650E+02, 1.750E+02, 1.850E+02, 1.950E+02, 2.000E+02, 2.050E+02, 2.150E+02, 2.250E+02, 2.350E+02, 2.450E+02, 2.550E+02, 2.650E+02, 2.750E+02, 2.850E+02, 2.950E+02, 3.050E+02, 3.150E+02, 3.250E+02, 3.350E+02, 3.450E+02, 3.550E+02, 3.650E+02, 3.750E+02, 3.850E+02, 3.950E+02, 4.050E+02, 4.150E+02, 4.250E+02, 4.350E+02, 4.450E+02, 4.550E+02, 4.650E+02, 4.750E+02, 4.850E+02, 4.950E+02, 5.050E+02, 5.150E+02, 5.250E+02, 5.350E+02, 5.450E+02, 5.550E+02, 5.650E+02, 5.750E+02, 5.850E+02, 5.950E+02, 6.050E+02, 6.150E+02, 6.250E+02, 6.350E+02, 6.450E+02, 6.550E+02, 6.650E+02, 6.750E+02, 6.850E+02, 6.950E+02, 7.050E+02, 7.150E+02, 7.250E+02, 7.350E+02, 7.450E+02, 7.550E+02, 7.650E+02, 7.750E+02, 7.850E+02, 7.950E+02, 8.050E+02, 8.150E+02, 8.250E+02, 8.350E+02, 8.450E+02, 8.550E+02, 8.650E+02, 8.750E+02, 8.850E+02, 8.950E+02, 9.050E+02, 9.150E+02, 9.250E+02, 9.350E+02, 9.450E+02, 9.550E+02, 9.650E+02, 9.750E+02, 9.850E+02, 9.950E+02, 1.005E+03, 1.015E+03, 1.025E+03, 1.035E+03, 1.045E+03, 1.055E+03, 1.065E+03, 1.075E+03, 1.085E+03, 1.095E+03, 1.105E+03, 1.115E+03, 1.125E+03, 1.135E+03, 1.145E+03, 1.155E+03, 1.165E+03, 1.175E+03, 1.185E+03, 1.195E+03, 1.205E+03, 1.215E+03, 1.225E+03, 1.235E+03, 1.245E+03, 1.255E+03, 1.265E+03, 1.275E+03, 1.285E+03, 1.295E+03, 1.305E+03, 1.315E+03, 1.325E+03, 1.335E+03, 1.345E+03, 1.355E+03, 1.365E+03, 1.375E+03, 1.385E+03, 1.395E+03, 1.405E+03, 1.415E+03, 1.425E+03, 1.435E+03, 1.445E+03, 1.455E+03, 1.465E+03, 1.475E+03, 1.485E+03, 1.495E+03, 1.505E+03, 1.515E+03, 1.525E+03, 1.535E+03, 1.545E+03, 1.555E+03, 1.565E+03, 1.575E+03, 1.585E+03, 1.595E+03, 1.605E+03, 1.615E+03, 1.625E+03, 1.635E+03, 1.645E+03, 1.655E+03, 1.665E+03, 1.675E+03, 1.685E+03, 1.695E+03, 1.705E+03, 1.715E+03, 1.725E+03, 1.735E+03, 1.745E+03, 1.755E+03, 1.765E+03, 1.775E+03, 1.785E+03, 1.795E+03, 1.805E+03, 1.815E+03, 1.825E+03, 1.835E+03, 1.845E+03, 1.855E+03, 1.865E+03, 1.875E+03, 1.885E+03, 1.895E+03, 1.905E+03, 1.915E+03, 1.925E+03, 1.935E+03, 1.945E+03, 1.955E+03, 1.965E+03, 1.975E+03, 1.985E+03, 1.995E+03, 2.005E+03, 2.015E+03, 2.025E+03, 2.035E+03, 2.045E+03, 2.055E+03, 2.065E+03, 2.075E+03, 2.085E+03, 2.095E+03, 2.105E+03, 2.115E+03, 2.125E+03, 2.135E+03, 2.145E+03, 2.155E+03, 2.165E+03, 2.175E+03, 2.185E+03, 2.195E+03, 2.205E+03, 2.215E+03, 2.225E+03, 2.235E+03, 2.245E+03, 2.255E+03, 2.265E+03, 2.275E+03, 2.285E+03, 2.295E+03, 2.305E+03, 2.315E+03, 2.325E+03, 2.335E+03, 2.345E+03, 2.355E+03, 2.365E+03, 2.375E+03, 2.385E+03, 2.395E+03, 2.405E+03, 2.415E+03, 2.425E+03, 2.435E+03, 2.445E+03, 2.455E+03, 2.465E+03, 2.475E+03, 2.485E+03, 2.495E+03, 2.505E+03, 2.515E+03, 2.525E+03, 2.535E+03, 2.545E+03, 2.555E+03, 2.565E+03, 2.575E+03, 2.585E+03, 2.595E+03, 2.605E+03, 2.615E+03, 2.625E+03, 2.635E+03, 2.645E+03, 2.655E+03, 2.665E+03, 2.675E+03, 2.685E+03, 2.695E+03, 2.705E+03, 2.715E+03, 2.725E+03, 2.735E+03, 2.745E+03, 2.755E+03, 2.765E+03, 2.775E+03, 2.785E+03, 2.795E+03, 2.805E+03, 2.815E+03, 2.825E+03, 2.835E+03, 2.845E+03, 2.855E+03, 2.865E+03, 2.875E+03, 2.885E+03, 2.895E+03, 2.905E+03, 2.915E+03, 2.925E+03, 2.935E+03, 2.945E+03, 2.955E+03, 2.965E+03, 2.975E+03, 2.985E+03, 2.995E+03, 3.005E+03, 3.015E+03, 3.025E+03, 3.035E+03, 3.045E+03, 3.055E+03, 3.065E+03, 3.075E+03, 3.085E+03, 3.095E+03, 3.105E+03, 3.115E+03, 3.125E+03, 3.135E+03, 3.145E+03, 3.155E+03, 3.165E+03, 3.175E+03, 3.185E+03, 3.195E+03, 3.205E+03, 3.215E+03, 3.225E+03, 3.235E+03, 3.245E+03, 3.255E+03, 3.265E+03, 3.275E+03, 3.285E+03, 3.295E+03, 3.305E+03, 3.315E+03, 3.325E+03, 3.335E+03, 3.345E+03, 3.355E+03, 3.365E+03, 3.375E+03, 3.385E+03, 3.395E+03, 3.405E+03, 3.415E+03, 3.425E+03, 3.435E+03, 3.445E+03, 3.455E+03, 3.465E+03, 3.475E+03, 3.485E+03, 3.495E+03, 3.505E+03, 3.515E+03, 3.525E+03, 3.535E+03, 3.545E+03, 3.555E+03, 3.565E+03, 3.575E+03, 3.585E+03, 3.595E+03, 3.605E+03, 3.615E+03, 3.625E+03, 3.635E+03, 3.645E+03, 3.655E+03, 3.665E+03, 3.675E+03, 3.685E+03, 3.695E+03, 3.705E+03, 3.715E+03, 3.725E+03, 3.735E+03, 3.745E+03, 3.755E+03, 3.765E+03, 3.775E+03, 3.785E+03, 3.795E+03, 3.805E+03, 3.815E+03, 3.825E+03, 3.835E+03, 3.845E+03, 3.855E+03, 3.865E+03, 3.875E+03, 3.885E+03, 3.895E+03, 3.905E+03, 3.915E+03, 3.925E+03, 3.935E+03, 3.945E+03, 3.955E+03, 3.965E+03, 3.975E+03, 3.985E+03, 3.995E+03, 4.005E+03, 4.015E+03, 4.025E+03, 4.035E+03, 4.045E+03, 4.055E+03, 4.065E+03, 4.075E+03, 4.085E+03, 4.095E+03, 4.105E+03, 4.115E+03, 4.125E+03, 4.135E+03, 4.145E+03, 4.155E+03, 4.165E+03, 4.175E+03, 4.185E+03, 4.195E+03, 4.205E+03, 4.215E+03, 4.225E+03, 4.235E+03, 4.245E+03, 4.255E+03, 4.265E+03, 4.275E+03, 4.285E+03, 4.295E+03, 4.305E+03, 4.315E+03, 4.325E+03, 4.335E+03, 4.345E+03, 4.355E+03, 4.365E+03, 4.375E+03, 4.385E+03, 4.395E+03, 4.405E+03, 4.415E+03, 4.425E+03, 4.435E+03, 4.445E+03, 4.455E+03, 4.465E+03, 4.475E+03, 4.485E+03, 4.495E+03, 4.505E+03, 4.515E+03, 4.525E+03, 4.535E+03, 4.545E+03, 4.555E+03, 4.565E+03, 4.575E+03, 4.585E+03, 4.595E+03, 4.605E+03, 4.615E+03, 4.625E+03, 4.635E+03, 4.645E+03, 4.655E+03, 4.665E+03, 4.675E+03, 4.685E+03, 4.695E+03, 4.705E+03, 4.715E+03, 4.725E+03, 4.735E+03, 4.745E+03, 4.755E+03, 4.765E+03, 4.775E+03, 4.785E+03, 4.795E+03, 4.805E+03, 4.815E+03, 4.825E+03, 4.835E+03, 4.845E+03, 4.855E+03, 4.865E+03, 4.875E+03, 4.885E+03, 4.895E+03, 4.905E+03, 4.915E+03, 4.925E+03, 4.935E+03, 4.945E+03, 4.955E+03, 4.965E+03, 4.975E+03, 4.985E+03, 4.995E+03, 5.005E+03, 5.015E+03, 5.025E+03, 5.035E+03, 5.045E+03, 5.055E+03, 5.065E+03, 5.075E+03, 5.085E+03, 5.095E+03, 5.105E+03, 5.115E+03, 5.125E+03, 5.135E+03, 5.145E+03, 5.155E+03, 5.165E+03, 5.175E+03, 5.185E+03, 5.195E+03, 5.205E+03, 5.215E+03, 5.225E+03, 5.235E+03, 5.245E+03, 5.255E+03, 5.265E+03, 5.275E+03, 5.285E+03, 5.295E+03, 5.305E+03, 5.315E+03, 5.325E+03, 5.335E+03, 5.345E+03, 5.355E+03, 5.365E+03, 5.375E+03, 5.385E+03, 5.395E+03, 5.405E+03, 5.415E+03, 5.425E+03, 5.435E+03, 5.445E+03, 5.455E+03, 5.465E+03, 5.475E+03, 5.485E+03, 5.495E+03, 5.505E+03, 5.515E+03, 5.525E+03, 5.535E+03, 5.545E+03, 5.555E+03, 5.565E+03, 5.575E+03, 5.585E+03, 5.595E+03, 5.605E+03, 5.615E+03, 5.625E+03, 5.635E+03, 5.645E+03, 5.655E+03, 5.665E+03, 5.675E+03, 5.685E+03, 5.695E+03, 5.705E+03, 5.715E+03, 5.725E+03, 5.735E+03, 5.745E+03, 5.755E+03, 5.765E+03, 5.775E+03, 5.785E+03, 5.795E+03, 5.805E+03, 5.815E+03, 5.825E+03, 5.835E+03, 5.845E+03, 5.855E+03, 5.865E+03, 5.875E+03, 5.885E+03, 5.895E+03, 5.905E+03, 5.915E+03, 5.925E+03, 5.935E+03, 5.945E+03, 5.955E+03, 5.965E+03, 5.975E+03, 5.985E+03, 5.995E+03, 6.005E+03, 6.015E+03, 6.025E+03, 6.035E+03, 6.045E+03, 6.055E+03, 6.065E+03, 6.075E+03, 6.085E+03, 6.095E+03, 6.105E+03, 6.115E+03, 6.125E+03, 6.135E+03, 6.145E+03, 6.155E+03, 6.165E+03, 6.175E+03, 6.185E+03, 6.195E+03, 6.205E+03, 6.215E+03, 6.225E+03, 6.235E+03, 6.245E+03, 6.255E+03, 6.265E+03, 6.275E+03, 6.285E+03, 6.295E+03, 6.305E+03, 6.315E+03, 6.325E+03, 6.335E+03, 6.345E+03, 6.355E+03, 6.365E+03, 6.375E+03, 6.385E+03, 6.395E+03, 6.405E+03, 6.415E+03, 6.425E+03, 6.435E+03, 6.445E+03, 6.455E+03, 6.465E+03, 6.475E+03, 6.485E+03, 6.495E+03, 6.505E+03, 6.515E+03, 6.525E+03, 6.535E+03, 6.545E+03, 6.555E+03, 6.565E+03, 6.575E+03, 6.585E+03, 6.595E+03, 6.605E+03, 6.615E+03, 6.625E+03, 6.635E+03, 6.645E+03, 6.655E+03, 6.665E+03, 6.675E+03, 6.685E+03, 6.695E+03, 6.705E+03, 6.715E+03, 6.725E+03, 6.735E+03, 6.745E+03, 6.755E+03, 6.765E+03, 6.775E+03, 6.785E+03, 6.795E+03, 6.805E+03, 6.815E+03, 6.825E+03, 6.835E+03, 6.845E+03, 6.855E+03, 6.865E+03, 6.875E+03, 6.885E+03, 6.895E+03, 6.905E+03, 6.915E+03, 6.925E+03, 6.935E+03, 6.945E+03, 6.955E+03, 6.965E+03, 6.975E+03, 6.985E+03, 6.995E+03, 7.005E+03, 7.015E+03, 7.025E+03, 7.035E+03, 7.045E+03, 7.055E+03, 7.065E+03, 7.075E+03, 7.085E+03, 7.095E+03, 7.105E+03, 7.115E+03, 7.125E+03, 7.135E+03, 7.145E+03, 7.155E+03, 7.165E+03, 7.175E+03, 7.185E+03, 7.195E+03, 7.205E+03, 7.215E+03, 7.225E+03, 7.413E+03, 7.245E+03, 7.255E+03, 7.265E+03, 7.275E+03, 7.285E+03, 7.295E+03, 7.305E+03, 7.315E+03, 7.325E+03, 7.335E+03, 7.345E+03, 7.355E+03, 7.365E+03, 7.375E+03, 7.385E+03, 7.395E+03, 7.405E+03, 7.415E+03, 7.425E+03, 7.435E+03, 7.445E+03, 7.455E+03, 7.465E+03, 7.475E+03, 7.485E+03, 7.495E+03, 7.505E+03, 7.515E+03, 7.525E+03, 7.535E+03, 7.545E+03, 7.555E+03, 7.565E+03, 7.575E+03, 7.585E+03, 7.595E+03, 7.605E+03, 7.615E+03, 7.625E+03, 7.635E+03, 7.645E+03, 7.655E+03, 7.665E+03, 7.675E+03, 7.685E+03, 7.695E+03, 7.705E+03, 7.715E+03, 7.725E+03, 7.735E+03, 7.745E+03, 7.755E+03, 7.765E+03, 7.775E+03, 7.785E+03, 7.795E+03, 7.805E+03, 7.815E+03, 7.825E+03, 7.835E+03, 7.845E+03, 7.855E+03, 7.865E+03, 7.875E+03, 7.885E+03, 7.895E+03, 7.905E+03, 7.915E+03, 7.925E+03, 7.935E+03, 7.945E+03, 7.955E+03, 7.965E+03, 7.975E+03, 7.985E+03, 7.995E+03, 8.005E+03, 8.015E+03, 8.025E+03, 8.035E+03, 8.045E+03, 8.055E+03, 8.065E+03, 8.075E+03, 8.085E+03, 8.095E+03, 8.105E+03, 8.115E+03, 8.125E+03, 8.135E+03, 8.145E+03, 8.155E+03, 8.165E+03, 8.175E+03, 8.185E+03, 8.195E+03, 8.205E+03, 8.215E+03, 8.225E+03, 8.235E+03, 8.245E+03, 8.255E+03, 8.265E+03, 8.275E+03, 8.285E+03, 8.295E+03, 8.305E+03, 8.315E+03, 8.325E+03, 8.335E+03, 8.345E+03, 8.355E+03, 8.365E+03, 8.375E+03, 8.385E+03, 8.395E+03, 8.405E+03, 8.415E+03, 8.425E+03, 8.435E+03, 8.445E+03, 8.455E+03, 8.465E+03, 8.475E+03, 8.485E+03, 8.495E+03, 8.505E+03, 8.515E+03, 8.525E+03, 8.535E+03, 8.545E+03, 8.555E+03, 8.565E+03, 8.575E+03, 8.585E+03, 8.595E+03, 8.605E+03, 8.615E+03, 8.625E+03, 8.635E+03, 8.645E+03, 8.655E+03, 8.665E+03, 8.675E+03, 8.685E+03, 8.695E+03, 8.705E+03, 8.715E+03, 8.725E+03, 8.735E+03, 8.745E+03, 8.755E+03, 8.765E+03, 8.775E+03, 8.785E+03, 8.795E+03, 8.805E+03, 8.815E+03, 8.825E+03, 8.835E+03, 8.845E+03, 8.855E+03, 8.865E+03, 8.875E+03, 8.885E+03, 8.895E+03, 8.905E+03, 8.915E+03, 8.925E+03, 8.935E+03, 8.945E+03, 8.955E+03, 8.965E+03, 8.975E+03, 8.985E+03, 8.995E+03, 9.005E+03};
double dmgProtonY[]={7.590E+03,5.170E+03,3.976E+03,2.778E+03,2.167E+03,1.650E+03,9.515E+02,6.824E+02,4.447E+02,3.338E+02,2.456E+02,1.336E+02,9.127E+01,5.862E+01,4.307E+01,3.133E+01,1.618E+01,1.101E+01,6.756E+00,5.140E+00,3.871E+00,3.376E+00,2.632E+00,2.558E+00,2.346E+00,2.125E+00,1.881E+00,1.715E+00,1.907E+00,1.580E+00,1.552E+00,1.468E+00,1.379E+00,1.307E+00,1.244E+00,1.276E+00,1.202E+00,1.161E+00,1.124E+00,1.099E+00,1.075E+00,1.051E+00,1.028E+00,1.013E+00,9.968E-01,9.525E-01,9.811E-01,9.653E-01,9.516E-01,9.421E-01,9.326E-01,9.242E-01,9.147E-01,9.063E-01,8.968E-01,8.905E-01,8.842E-01,8.779E-01,8.726E-01,8.663E-01,8.611E-01,8.558E-01,8.495E-01,8.442E-01,8.389E-01,8.326E-01,8.274E-01,8.221E-01,8.168E-01,8.116E-01,8.063E-01,8.011E-01,7.958E-01,7.905E-01,7.853E-01,7.821E-01,7.800E-01,7.768E-01,7.737E-01,7.716E-01,7.684E-01,7.653E-01,7.621E-01,7.600E-01,7.568E-01,7.547E-01,7.516E-01,7.484E-01,7.463E-01,7.442E-01,7.421E-01,7.400E-01,7.379E-01,7.358E-01,7.347E-01,7.326E-01,7.305E-01,7.284E-01,7.263E-01,7.242E-01,7.232E-01,7.211E-01,7.189E-01,7.168E-01,7.147E-01,7.137E-01,7.116E-01,7.095E-01,7.084E-01,7.074E-01,7.063E-01,7.042E-01,7.032E-01,7.021E-01,7.011E-01,6.989E-01,6.979E-01,6.968E-01,6.958E-01,6.947E-01,6.926E-01,6.916E-01,6.905E-01,6.895E-01,6.874E-01,6.863E-01,6.853E-01,6.842E-01,6.832E-01,6.811E-01,6.800E-01,6.789E-01,6.789E-01,6.779E-01,6.768E-01,6.758E-01,6.758E-01,6.747E-01,6.737E-01,6.726E-01,6.726E-01,6.716E-01,6.705E-01,6.695E-01,6.695E-01,6.684E-01,6.674E-01,6.663E-01,6.663E-01,6.653E-01,6.642E-01,6.632E-01,6.632E-01,6.621E-01,6.611E-01,6.600E-01,6.600E-01,6.589E-01,6.579E-01,6.568E-01,6.568E-01,6.558E-01,6.547E-01,6.547E-01,6.537E-01,6.526E-01,6.526E-01,6.516E-01,6.516E-01,6.505E-01,6.495E-01,6.495E-01,6.484E-01,6.474E-01,6.474E-01,6.463E-01,6.463E-01,6.453E-01,6.442E-01,6.442E-01,6.432E-01,6.432E-01,6.421E-01,6.411E-01,6.411E-01,6.400E-01,6.389E-01,6.389E-01,6.379E-01,6.379E-01,6.368E-01,6.358E-01,6.358E-01,6.347E-01,6.347E-01,6.337E-01,6.326E-01,6.326E-01,6.316E-01,6.316E-01,6.305E-01,6.295E-01,6.295E-01,6.284E-01,6.284E-01,6.284E-01,6.274E-01,6.274E-01,6.263E-01,6.263E-01,6.253E-01,6.253E-01,6.242E-01,6.242E-01,6.232E-01,6.232E-01,6.221E-01,6.221E-01,6.211E-01,6.211E-01,6.211E-01,6.200E-01,6.200E-01,6.189E-01,6.189E-01,6.179E-01,6.179E-01,6.168E-01,6.168E-01,6.158E-01,6.158E-01,6.158E-01,6.147E-01,6.147E-01,6.137E-01,6.137E-01,6.126E-01,6.126E-01,6.116E-01,6.116E-01,6.105E-01,6.105E-01,6.105E-01,6.095E-01,6.095E-01,6.084E-01,6.084E-01,6.074E-01,6.074E-01,6.063E-01,6.063E-01,6.053E-01,6.053E-01,6.053E-01,6.053E-01,6.042E-01,6.042E-01,6.042E-01,6.042E-01,6.042E-01,6.042E-01,6.042E-01,6.032E-01,6.032E-01,6.032E-01,6.032E-01,6.032E-01,6.032E-01,6.032E-01,6.021E-01,6.021E-01,6.021E-01,6.021E-01,6.021E-01,6.021E-01,6.021E-01,6.011E-01,6.011E-01,6.011E-01,6.011E-01,6.011E-01,6.011E-01,6.000E-01,6.000E-01,6.000E-01,6.000E-01,6.000E-01,6.000E-01,6.000E-01,5.989E-01,5.989E-01,5.989E-01,5.989E-01,5.989E-01,5.989E-01,5.989E-01,5.979E-01,5.979E-01,5.979E-01,5.979E-01,5.979E-01,5.979E-01,5.979E-01,5.968E-01,5.968E-01,5.968E-01,5.968E-01,5.968E-01,5.968E-01,5.958E-01,5.958E-01,5.958E-01,5.958E-01,5.958E-01,5.958E-01,5.958E-01,5.947E-01,5.947E-01,5.947E-01,5.947E-01,5.947E-01,5.947E-01,5.947E-01,5.937E-01,5.937E-01,5.937E-01,5.937E-01,5.937E-01,5.937E-01,5.937E-01,5.926E-01,5.926E-01,5.926E-01,5.926E-01,5.926E-01,5.926E-01,5.926E-01,5.916E-01,5.916E-01,5.916E-01,5.916E-01,5.916E-01,5.916E-01,5.905E-01,5.905E-01,5.905E-01,5.905E-01,5.905E-01,5.905E-01,5.905E-01,5.895E-01,5.895E-01,5.895E-01,5.895E-01,5.895E-01,5.895E-01,5.895E-01,5.884E-01,5.884E-01,5.884E-01,5.884E-01,5.884E-01,5.884E-01,5.884E-01,5.874E-01,5.874E-01,5.874E-01,5.874E-01,5.874E-01,5.874E-01,5.874E-01,5.863E-01,5.863E-01,5.863E-01,5.863E-01,5.863E-01,5.863E-01,5.863E-01,5.853E-01,5.853E-01,5.853E-01,5.853E-01,5.853E-01,5.853E-01,5.853E-01,5.842E-01,5.842E-01,5.842E-01,5.842E-01,5.842E-01,5.842E-01,5.842E-01,5.832E-01,5.832E-01,5.832E-01,5.832E-01,5.832E-01,5.832E-01,5.832E-01,5.821E-01,5.821E-01,5.821E-01,5.821E-01,5.821E-01,5.821E-01,5.821E-01,5.811E-01,5.811E-01,5.811E-01,5.811E-01,5.811E-01,5.811E-01,5.811E-01,5.800E-01,5.800E-01,5.800E-01,5.800E-01,5.800E-01,5.800E-01,5.800E-01,5.789E-01,5.789E-01,5.789E-01,5.789E-01,5.789E-01,5.789E-01,5.789E-01,5.779E-01,5.779E-01,5.779E-01,5.779E-01,5.779E-01,5.779E-01,5.779E-01,5.768E-01,5.768E-01,5.768E-01,5.768E-01,5.768E-01,5.768E-01,5.768E-01,5.758E-01,5.758E-01,5.758E-01,5.758E-01,5.758E-01,5.758E-01,5.758E-01,5.747E-01,5.747E-01,5.747E-01,5.747E-01,5.747E-01,5.747E-01,5.747E-01,5.737E-01,5.737E-01,5.737E-01,5.737E-01,5.737E-01,5.737E-01,5.737E-01,5.726E-01,5.726E-01,5.726E-01,5.726E-01,5.726E-01,5.726E-01,5.726E-01,5.716E-01,5.716E-01,5.716E-01,5.716E-01,5.716E-01,5.716E-01,5.716E-01,5.705E-01,5.705E-01,5.705E-01,5.705E-01,5.705E-01,5.705E-01,5.705E-01,5.695E-01,5.695E-01,5.695E-01,5.695E-01,5.695E-01,5.695E-01,5.695E-01,5.684E-01,5.684E-01,5.684E-01,5.684E-01,5.684E-01,5.684E-01,5.684E-01,5.674E-01,5.674E-01,5.674E-01,5.674E-01,5.674E-01,5.674E-01,5.674E-01,5.674E-01,5.663E-01,5.663E-01,5.663E-01,5.663E-01,5.663E-01,5.663E-01,5.663E-01,5.653E-01,5.653E-01,5.653E-01,5.653E-01,5.653E-01,5.653E-01,5.653E-01,5.642E-01,5.642E-01,5.642E-01,5.642E-01,5.642E-01,5.642E-01,5.642E-01,5.632E-01,5.632E-01,5.632E-01,5.632E-01,5.632E-01,5.632E-01,5.632E-01,5.621E-01,5.621E-01,5.621E-01,5.621E-01,5.621E-01,5.621E-01,5.621E-01,5.611E-01,5.611E-01,5.611E-01,5.611E-01,5.611E-01,5.611E-01,5.611E-01,5.611E-01,5.600E-01,5.600E-01,5.600E-01,5.600E-01,5.600E-01,5.600E-01,5.600E-01,5.589E-01,5.589E-01,5.589E-01,5.589E-01,5.589E-01,5.589E-01,5.589E-01,5.579E-01,5.579E-01,5.579E-01,5.579E-01,5.579E-01,5.579E-01,5.579E-01,5.568E-01,5.568E-01,5.568E-01,5.568E-01,5.568E-01,5.568E-01,5.568E-01,5.568E-01,5.558E-01,5.558E-01,5.558E-01,5.558E-01,5.558E-01,5.558E-01,5.558E-01,5.547E-01,5.547E-01,5.547E-01,5.547E-01,5.547E-01,5.547E-01,5.547E-01,5.537E-01,5.537E-01,5.537E-01,5.537E-01,5.537E-01,5.537E-01,5.537E-01,5.537E-01,5.526E-01,5.526E-01,5.526E-01,5.526E-01,5.526E-01,5.526E-01,5.526E-01,5.516E-01,5.516E-01,5.516E-01,5.516E-01,5.516E-01,5.516E-01,5.516E-01,5.505E-01,5.505E-01,5.505E-01,5.505E-01,5.505E-01,5.505E-01,5.505E-01,5.505E-01,5.495E-01,5.495E-01,5.495E-01,5.495E-01,5.495E-01,5.495E-01,5.495E-01,5.484E-01,5.484E-01,5.484E-01,5.484E-01,5.484E-01,5.484E-01,5.484E-01,5.474E-01,5.474E-01,5.474E-01,5.474E-01,5.474E-01,5.474E-01,5.474E-01,5.474E-01,5.463E-01,5.463E-01,5.463E-01,5.463E-01,5.463E-01,5.463E-01,5.463E-01,5.453E-01,5.453E-01,5.453E-01,5.453E-01,5.453E-01,5.453E-01,5.453E-01,5.453E-01,5.442E-01,5.442E-01,5.442E-01,5.442E-01,5.442E-01,5.442E-01,5.442E-01,5.432E-01,5.432E-01,5.432E-01,5.432E-01,5.432E-01,5.432E-01,5.432E-01,5.432E-01,5.421E-01,5.421E-01,5.421E-01,5.421E-01,5.421E-01,5.421E-01,5.421E-01,5.411E-01,5.411E-01,5.411E-01,5.411E-01,5.411E-01,5.411E-01,5.411E-01,5.411E-01,5.400E-01,5.400E-01,5.400E-01,5.400E-01,5.400E-01,5.400E-01,5.400E-01,5.389E-01,5.389E-01,5.389E-01,5.389E-01,5.389E-01,5.389E-01,5.389E-01,5.389E-01,5.379E-01,5.379E-01,5.379E-01,5.379E-01,5.379E-01,5.379E-01,5.379E-01,5.368E-01,5.368E-01,5.368E-01,5.368E-01,5.368E-01,5.368E-01,5.368E-01,5.368E-01,5.358E-01,5.358E-01,5.358E-01,5.358E-01,5.358E-01,5.358E-01,5.358E-01,5.347E-01,5.347E-01,5.347E-01,5.347E-01,5.347E-01,5.347E-01,5.347E-01,5.347E-01,5.337E-01,5.337E-01,5.337E-01,5.337E-01,5.337E-01,5.337E-01,5.337E-01,5.326E-01,5.326E-01,5.326E-01,5.326E-01,5.326E-01,5.326E-01,5.326E-01,5.326E-01,5.316E-01,5.316E-01,5.316E-01,5.316E-01,5.316E-01,5.316E-01,5.316E-01,5.316E-01,5.305E-01,5.305E-01,5.305E-01,5.305E-01,5.305E-01,5.305E-01,5.305E-01,5.295E-01,5.295E-01,5.295E-01,5.295E-01,5.295E-01,5.295E-01,5.295E-01,5.295E-01,5.284E-01,5.284E-01,5.284E-01,5.284E-01,5.284E-01,5.284E-01,5.284E-01,5.284E-01,5.274E-01,5.274E-01,5.274E-01,5.274E-01,5.274E-01,5.274E-01,5.274E-01,5.263E-01,5.263E-01,5.263E-01,5.263E-01,5.263E-01,5.263E-01,5.263E-01,5.263E-01,5.253E-01,5.253E-01,5.253E-01,5.253E-01,5.253E-01,5.253E-01,5.253E-01,5.253E-01,5.242E-01,5.242E-01,5.242E-01,5.242E-01,5.242E-01,5.242E-01,5.242E-01,5.242E-01,5.232E-01,5.232E-01,5.232E-01,5.232E-01,5.232E-01,5.232E-01,5.232E-01,5.221E-01,5.221E-01,5.221E-01,5.221E-01,5.221E-01,5.221E-01,5.221E-01,5.221E-01,5.211E-01,5.211E-01,5.211E-01,5.211E-01,5.211E-01,5.211E-01,5.211E-01,5.211E-01,5.200E-01,5.200E-01,5.200E-01,5.200E-01,5.200E-01,5.200E-01,5.200E-01,5.200E-01,5.189E-01,5.189E-01,5.189E-01,5.189E-01,5.189E-01,5.189E-01,5.189E-01,5.179E-01,5.179E-01,5.179E-01,5.179E-01,5.179E-01,5.179E-01,5.179E-01,5.179E-01,5.168E-01,5.168E-01,5.168E-01,5.168E-01,5.168E-01,5.168E-01,5.168E-01,5.168E-01,5.158E-01,5.158E-01,5.158E-01,5.158E-01,5.158E-01,5.158E-01,5.158E-01,5.158E-01,5.147E-01,5.147E-01,5.147E-01,5.147E-01,5.147E-01,5.147E-01,5.147E-01,5.147E-01,5.137E-01,5.137E-01,5.137E-01,5.137E-01,5.137E-01,5.137E-01,5.137E-01,5.137E-01,5.126E-01,5.126E-01,5.126E-01,5.126E-01,5.126E-01,5.126E-01,5.126E-01,5.126E-01,5.116E-01,5.116E-01,5.116E-01,5.116E-01,5.116E-01,5.116E-01,5.116E-01,5.116E-01,5.105E-01,5.105E-01,5.105E-01,5.105E-01,5.105E-01,5.105E-01,5.105E-01,5.095E-01,5.095E-01,5.095E-01,5.095E-01,5.095E-01,5.095E-01,5.095E-01,5.095E-01,5.084E-01,5.084E-01,5.084E-01,5.084E-01,5.084E-01,5.084E-01,5.084E-01,5.084E-01};
double dmgPionX[]={1.500E+01,2.500E+01,3.500E+01,4.500E+01,5.500E+01,6.500E+01,7.500E+01,8.500E+01,9.500E+01,1.050E+02,1.150E+02,1.250E+02,1.350E+02,1.450E+02,1.550E+02,1.650E+02,1.750E+02,1.850E+02,1.950E+02,2.050E+02,2.150E+02,2.250E+02,2.350E+02,2.450E+02,2.550E+02,2.650E+02,2.750E+02,2.850E+02,2.950E+02,3.050E+02,3.150E+02,3.250E+02,3.350E+02,3.450E+02,3.550E+02,3.650E+02,3.750E+02,3.850E+02,3.950E+02,4.050E+02,4.150E+02,4.250E+02,4.350E+02,4.450E+02,4.550E+02,4.650E+02,4.750E+02,4.850E+02,4.950E+02,5.050E+02,5.150E+02,5.250E+02,5.350E+02,5.450E+02,5.550E+02,5.650E+02,5.750E+02,5.850E+02,5.950E+02,6.050E+02,6.150E+02,6.250E+02,6.350E+02,6.450E+02,6.550E+02,6.650E+02,6.750E+02,6.850E+02,6.950E+02,7.050E+02,7.150E+02,7.250E+02,7.350E+02,7.450E+02,7.550E+02,7.650E+02,7.750E+02,7.850E+02,7.950E+02,8.050E+02,8.150E+02,8.250E+02,8.350E+02,8.450E+02,8.550E+02,8.650E+02,8.750E+02,8.850E+02,8.950E+02,9.050E+02,9.150E+02,9.250E+02,9.350E+02,9.450E+02,9.550E+02,9.650E+02,9.750E+02,9.850E+02,9.950E+02,1.005E+03,1.015E+03,1.025E+03,1.035E+03,1.045E+03,1.055E+03,1.065E+03,1.075E+03,1.085E+03,1.095E+03,1.105E+03,1.115E+03,1.125E+03,1.135E+03,1.145E+03,1.155E+03,1.165E+03,1.175E+03,1.185E+03,1.195E+03,1.205E+03,1.215E+03,1.225E+03,1.235E+03,1.245E+03,1.255E+03,1.265E+03,1.275E+03,1.285E+03,1.295E+03,1.305E+03,1.315E+03,1.325E+03,1.335E+03,1.345E+03,1.355E+03,1.365E+03,1.375E+03,1.385E+03,1.395E+03,1.405E+03,1.415E+03,1.425E+03,1.435E+03,1.445E+03,1.455E+03,1.465E+03,1.475E+03,1.485E+03,1.495E+03,1.505E+03,1.515E+03,1.525E+03,1.535E+03,1.545E+03,1.555E+03,1.565E+03,1.575E+03,1.585E+03,1.595E+03,1.605E+03,1.615E+03,1.625E+03,1.635E+03,1.645E+03,1.655E+03,1.665E+03,1.675E+03,1.685E+03,1.695E+03,1.705E+03,1.715E+03,1.725E+03,1.735E+03,1.745E+03,1.755E+03,1.765E+03,1.775E+03,1.785E+03,1.795E+03,1.805E+03,1.815E+03,1.825E+03,1.835E+03,1.845E+03,1.855E+03,1.865E+03,1.875E+03,1.885E+03,1.895E+03,1.905E+03,1.915E+03,1.925E+03,1.935E+03,1.945E+03,1.955E+03,1.965E+03,1.975E+03,1.985E+03,1.995E+03,2.005E+03,2.015E+03,2.025E+03,2.035E+03,2.045E+03,2.055E+03,2.065E+03,2.075E+03,2.085E+03,2.095E+03,2.105E+03,2.115E+03,2.125E+03,2.135E+03,2.145E+03,2.155E+03,2.165E+03,2.175E+03,2.185E+03,2.195E+03,2.205E+03,2.215E+03,2.225E+03,2.235E+03,2.245E+03,2.255E+03,2.265E+03,2.275E+03,2.285E+03,2.295E+03,2.305E+03,2.315E+03,2.325E+03,2.335E+03,2.345E+03,2.355E+03,2.365E+03,2.375E+03,2.385E+03,2.395E+03,2.405E+03,2.415E+03,2.425E+03,2.435E+03,2.445E+03,2.455E+03,2.465E+03,2.475E+03,2.485E+03,2.495E+03,2.505E+03,2.515E+03,2.525E+03,2.535E+03,2.545E+03,2.555E+03,2.565E+03,2.575E+03,2.585E+03,2.595E+03,2.605E+03,2.615E+03,2.625E+03,2.635E+03,2.645E+03,2.655E+03,2.665E+03,2.675E+03,2.685E+03,2.695E+03,2.705E+03,2.715E+03,2.725E+03,2.735E+03,2.745E+03,2.755E+03,2.765E+03,2.775E+03,2.785E+03,2.795E+03,2.805E+03,2.815E+03,2.825E+03,2.835E+03,2.845E+03,2.855E+03,2.865E+03,2.875E+03,2.885E+03,2.895E+03,2.905E+03,2.915E+03,2.925E+03,2.935E+03,2.945E+03,2.955E+03,2.965E+03,2.975E+03,2.985E+03,2.995E+03,3.005E+03,3.015E+03,3.025E+03,3.035E+03,3.045E+03,3.055E+03,3.065E+03,3.075E+03,3.085E+03,3.095E+03,3.105E+03,3.115E+03,3.125E+03,3.135E+03,3.145E+03,3.155E+03,3.165E+03,3.175E+03,3.185E+03,3.195E+03,3.205E+03,3.215E+03,3.225E+03,3.235E+03,3.245E+03,3.255E+03,3.265E+03,3.275E+03,3.285E+03,3.295E+03,3.305E+03,3.315E+03,3.325E+03,3.335E+03,3.345E+03,3.355E+03,3.365E+03,3.375E+03,3.385E+03,3.395E+03,3.405E+03,3.415E+03,3.425E+03,3.435E+03,3.445E+03,3.455E+03,3.465E+03,3.475E+03,3.485E+03,3.495E+03,3.505E+03,3.515E+03,3.525E+03,3.535E+03,3.545E+03,3.555E+03,3.565E+03,3.575E+03,3.585E+03,3.595E+03,3.605E+03,3.615E+03,3.625E+03,3.635E+03,3.645E+03,3.655E+03,3.665E+03,3.675E+03,3.685E+03,3.695E+03,3.705E+03,3.715E+03,3.725E+03,3.735E+03,3.745E+03,3.755E+03,3.765E+03,3.775E+03,3.785E+03,3.795E+03,3.805E+03,3.815E+03,3.825E+03,3.835E+03,3.845E+03,3.855E+03,3.865E+03,3.875E+03,3.885E+03,3.895E+03,3.905E+03,3.915E+03,3.925E+03,3.935E+03,3.945E+03,3.955E+03,3.965E+03,3.975E+03,3.985E+03,3.995E+03,4.005E+03,4.015E+03,4.025E+03,4.035E+03,4.045E+03,4.055E+03,4.065E+03,4.075E+03,4.085E+03,4.095E+03,4.105E+03,4.115E+03,4.125E+03,4.135E+03,4.145E+03,4.155E+03,4.165E+03,4.175E+03,4.185E+03,4.195E+03,4.205E+03,4.215E+03,4.225E+03,4.235E+03,4.245E+03,4.255E+03,4.265E+03,4.275E+03,4.285E+03,4.295E+03,4.305E+03,4.315E+03,4.325E+03,4.335E+03,4.345E+03,4.355E+03,4.365E+03,4.375E+03,4.385E+03,4.395E+03,4.405E+03,4.415E+03,4.425E+03,4.435E+03,4.445E+03,4.455E+03,4.465E+03,4.475E+03,4.485E+03,4.495E+03,4.505E+03,4.515E+03,4.525E+03,4.535E+03,4.545E+03,4.555E+03,4.565E+03,4.575E+03,4.585E+03,4.595E+03,4.605E+03,4.615E+03,4.625E+03,4.635E+03,4.645E+03,4.655E+03,4.665E+03,4.675E+03,4.685E+03,4.695E+03,4.705E+03,4.715E+03,4.725E+03,4.735E+03,4.745E+03,4.755E+03,4.765E+03,4.775E+03,4.785E+03,4.795E+03,4.805E+03,4.815E+03,4.825E+03,4.835E+03,4.845E+03,4.855E+03,4.865E+03,4.875E+03,4.885E+03,4.895E+03,4.905E+03,4.915E+03,4.925E+03,4.935E+03,4.945E+03,4.955E+03,4.965E+03,4.975E+03,4.985E+03,4.995E+03,5.005E+03,5.015E+03,5.025E+03,5.035E+03,5.045E+03,5.055E+03,5.065E+03,5.075E+03,5.085E+03,5.095E+03,5.105E+03,5.115E+03,5.125E+03,5.135E+03,5.145E+03,5.155E+03,5.165E+03,5.175E+03,5.185E+03,5.195E+03,5.205E+03,5.215E+03,5.225E+03,5.235E+03,5.245E+03,5.255E+03,5.265E+03,5.275E+03,5.285E+03,5.295E+03,5.305E+03,5.315E+03,5.325E+03,5.335E+03,5.345E+03,5.355E+03,5.365E+03,5.375E+03,5.385E+03,5.395E+03,5.405E+03,5.415E+03,5.425E+03,5.435E+03,5.445E+03,5.455E+03,5.465E+03,5.475E+03,5.485E+03,5.495E+03,5.505E+03,5.515E+03,5.525E+03,5.535E+03,5.545E+03,5.555E+03,5.565E+03,5.575E+03,5.585E+03,5.595E+03,5.605E+03,5.615E+03,5.625E+03,5.635E+03,5.645E+03,5.655E+03,5.665E+03,5.675E+03,5.685E+03,5.695E+03,5.705E+03,5.715E+03,5.725E+03,5.735E+03,5.745E+03,5.755E+03,5.765E+03,5.775E+03,5.785E+03,5.795E+03,5.805E+03,5.815E+03,5.825E+03,5.835E+03,5.845E+03,5.855E+03,5.865E+03,5.875E+03,5.885E+03,5.895E+03,5.905E+03,5.915E+03,5.925E+03,5.935E+03,5.945E+03,5.955E+03,5.965E+03,5.975E+03,5.985E+03,5.995E+03,6.005E+03,6.015E+03,6.025E+03,6.035E+03,6.045E+03,6.055E+03,6.065E+03,6.075E+03,6.085E+03,6.095E+03,6.105E+03,6.115E+03,6.125E+03,6.135E+03,6.145E+03,6.155E+03,6.165E+03,6.175E+03,6.185E+03,6.195E+03,6.205E+03,6.215E+03,6.225E+03,6.235E+03,6.245E+03,6.255E+03,6.265E+03,6.275E+03,6.285E+03,6.295E+03,6.305E+03,6.315E+03,6.325E+03,6.335E+03,6.345E+03,6.355E+03,6.365E+03,6.375E+03,6.385E+03,6.395E+03,6.405E+03,6.415E+03,6.425E+03,6.435E+03,6.445E+03,6.455E+03,6.465E+03,6.475E+03,6.485E+03,6.495E+03,6.505E+03,6.515E+03,6.525E+03,6.535E+03,6.545E+03,6.555E+03,6.565E+03,6.575E+03,6.585E+03,6.595E+03,6.605E+03,6.615E+03,6.625E+03,6.635E+03,6.645E+03,6.655E+03,6.665E+03,6.675E+03,6.685E+03,6.695E+03,6.705E+03,6.715E+03,6.725E+03,6.735E+03,6.745E+03,6.755E+03,6.765E+03,6.775E+03,6.785E+03,6.795E+03,6.805E+03,6.815E+03,6.825E+03,6.835E+03,6.845E+03,6.855E+03,6.865E+03,6.875E+03,6.885E+03,6.895E+03,6.905E+03,6.915E+03,6.925E+03,6.935E+03,6.945E+03,6.955E+03,6.965E+03,6.975E+03,6.985E+03,6.995E+03,7.005E+03,7.015E+03,7.025E+03,7.035E+03,7.045E+03,7.055E+03,7.065E+03,7.075E+03,7.085E+03,7.095E+03,7.105E+03,7.115E+03,7.125E+03,7.135E+03,7.145E+03,7.155E+03,7.165E+03,7.175E+03,7.185E+03,7.195E+03,7.205E+03,7.215E+03,7.225E+03,7.235E+03,7.245E+03,7.255E+03,7.265E+03,7.275E+03,7.285E+03,7.295E+03,7.305E+03,7.315E+03,7.325E+03,7.335E+03,7.345E+03,7.355E+03,7.365E+03,7.375E+03,7.385E+03,7.395E+03,7.405E+03,7.415E+03,7.425E+03,7.435E+03,7.445E+03,7.455E+03,7.465E+03,7.475E+03,7.485E+03,7.495E+03,7.505E+03,7.515E+03,7.525E+03,7.535E+03,7.545E+03,7.555E+03,7.565E+03,7.575E+03,7.585E+03,7.595E+03,7.605E+03,7.615E+03,7.625E+03,7.635E+03,7.645E+03,7.655E+03,7.665E+03,7.675E+03,7.685E+03,7.695E+03,7.705E+03,7.715E+03,7.725E+03,7.735E+03,7.745E+03,7.755E+03,7.765E+03,7.775E+03,7.785E+03,7.795E+03,7.805E+03,7.815E+03,7.825E+03,7.835E+03,7.845E+03,7.855E+03,7.865E+03,7.875E+03,7.885E+03,7.895E+03,7.905E+03,7.915E+03,7.925E+03,7.935E+03,7.945E+03,7.955E+03,7.965E+03,7.975E+03,7.985E+03,7.995E+03,8.005E+03,8.015E+03,8.025E+03,8.035E+03,8.045E+03,8.055E+03,8.065E+03,8.075E+03,8.085E+03,8.095E+03,8.105E+03,8.115E+03,8.125E+03,8.135E+03,8.145E+03,8.155E+03,8.165E+03,8.175E+03,8.185E+03,8.195E+03,8.205E+03,8.215E+03,8.225E+03,8.235E+03,8.245E+03,8.255E+03,8.265E+03,8.275E+03,8.285E+03,8.295E+03,8.305E+03,8.315E+03,8.325E+03,8.335E+03,8.345E+03,8.355E+03,8.365E+03,8.375E+03,8.385E+03,8.395E+03,8.405E+03,8.415E+03,8.425E+03,8.435E+03,8.445E+03,8.455E+03,8.465E+03,8.475E+03,8.485E+03,8.495E+03,8.505E+03,8.515E+03,8.525E+03,8.535E+03,8.545E+03,8.555E+03,8.565E+03,8.575E+03,8.585E+03,8.595E+03,8.605E+03,8.615E+03,8.625E+03,8.635E+03,8.645E+03,8.655E+03,8.665E+03,8.675E+03,8.685E+03,8.695E+03,8.705E+03,8.715E+03,8.725E+03,8.735E+03,8.745E+03,8.755E+03,8.765E+03,8.775E+03,8.785E+03,8.795E+03,8.805E+03,8.815E+03,8.825E+03,8.835E+03,8.845E+03,8.855E+03,8.865E+03,8.875E+03,8.885E+03,8.895E+03,8.905E+03,8.915E+03,8.925E+03,8.935E+03,8.945E+03,8.955E+03,8.965E+03,8.975E+03,8.985E+03,8.995E+03,9.005E+03};
double dmgPionY[]={5.589E-01,4.947E-01,4.663E-01,4.937E-01,5.274E-01,5.632E-01,6.042E-01,6.411E-01,7.011E-01,7.411E-01,7.663E-01,7.926E-01,8.242E-01,8.558E-01,8.863E-01,9.084E-01,9.221E-01,9.347E-01,9.432E-01,9.463E-01,9.453E-01,9.400E-01,9.316E-01,9.189E-01,9.042E-01,8.874E-01,8.684E-01,8.484E-01,8.274E-01,8.074E-01,7.895E-01,7.705E-01,7.526E-01,7.347E-01,7.168E-01,7.011E-01,6.853E-01,6.705E-01,6.558E-01,6.421E-01,6.295E-01,6.179E-01,6.063E-01,5.947E-01,5.853E-01,5.768E-01,5.684E-01,5.611E-01,5.558E-01,5.537E-01,5.526E-01,5.537E-01,5.547E-01,5.600E-01,5.653E-01,5.674E-01,5.684E-01,5.642E-01,5.600E-01,5.516E-01,5.442E-01,5.347E-01,5.253E-01,5.158E-01,5.063E-01,4.989E-01,4.926E-01,4.895E-01,4.863E-01,4.874E-01,4.895E-01,4.926E-01,4.968E-01,5.032E-01,5.126E-01,5.221E-01,5.326E-01,5.432E-01,5.537E-01,5.621E-01,5.695E-01,5.747E-01,5.800E-01,5.842E-01,5.863E-01,5.874E-01,5.853E-01,5.789E-01,5.716E-01,5.653E-01,5.600E-01,5.537E-01,5.474E-01,5.421E-01,5.368E-01,5.316E-01,5.274E-01,5.263E-01,5.242E-01,5.232E-01,5.211E-01,5.200E-01,5.179E-01,5.158E-01,5.147E-01,5.126E-01,5.137E-01,5.147E-01,5.168E-01,5.189E-01,5.211E-01,5.221E-01,5.242E-01,5.263E-01,5.284E-01,5.295E-01,5.316E-01,5.337E-01,5.347E-01,5.368E-01,5.379E-01,5.400E-01,5.411E-01,5.432E-01,5.442E-01,5.463E-01,5.453E-01,5.453E-01,5.442E-01,5.432E-01,5.432E-01,5.421E-01,5.411E-01,5.400E-01,5.400E-01,5.389E-01,5.368E-01,5.347E-01,5.326E-01,5.305E-01,5.284E-01,5.274E-01,5.253E-01,5.232E-01,5.211E-01,5.189E-01,5.168E-01,5.158E-01,5.137E-01,5.126E-01,5.105E-01,5.084E-01,5.074E-01,5.053E-01,5.032E-01,5.021E-01,5.011E-01,4.989E-01,4.979E-01,4.968E-01,4.958E-01,4.947E-01,4.937E-01,4.926E-01,4.916E-01,4.905E-01,4.895E-01,4.895E-01,4.884E-01,4.874E-01,4.874E-01,4.863E-01,4.863E-01,4.853E-01,4.842E-01,4.842E-01,4.842E-01,4.832E-01,4.832E-01,4.832E-01,4.832E-01,4.832E-01,4.832E-01,4.832E-01,4.832E-01,4.821E-01,4.821E-01,4.821E-01,4.821E-01,4.811E-01,4.811E-01,4.811E-01,4.800E-01,4.800E-01,4.800E-01,4.800E-01,4.789E-01,4.789E-01,4.789E-01,4.779E-01,4.779E-01,4.779E-01,4.779E-01,4.768E-01,4.768E-01,4.768E-01,4.758E-01,4.758E-01,4.758E-01,4.747E-01,4.747E-01,4.747E-01,4.747E-01,4.737E-01,4.737E-01,4.737E-01,4.726E-01,4.726E-01,4.726E-01,4.716E-01,4.716E-01,4.716E-01,4.716E-01,4.705E-01,4.705E-01,4.705E-01,4.695E-01,4.695E-01,4.695E-01,4.695E-01,4.695E-01,4.695E-01,4.695E-01,4.695E-01,4.695E-01,4.695E-01,4.695E-01,4.695E-01,4.695E-01,4.695E-01,4.695E-01,4.695E-01,4.684E-01,4.684E-01,4.684E-01,4.684E-01,4.684E-01,4.684E-01,4.684E-01,4.684E-01,4.674E-01,4.674E-01,4.674E-01,4.674E-01,4.674E-01,4.674E-01,4.674E-01,4.663E-01,4.663E-01,4.663E-01,4.663E-01,4.663E-01,4.663E-01,4.663E-01,4.663E-01,4.653E-01,4.653E-01,4.653E-01,4.653E-01,4.653E-01,4.653E-01,4.653E-01,4.642E-01,4.642E-01,4.642E-01,4.642E-01,4.642E-01,4.642E-01,4.642E-01,4.642E-01,4.632E-01,4.632E-01,4.632E-01,4.632E-01,4.632E-01,4.632E-01,4.632E-01,4.632E-01,4.632E-01,4.632E-01,4.632E-01,4.632E-01,4.632E-01,4.621E-01,4.621E-01,4.621E-01,4.621E-01,4.621E-01,4.621E-01,4.621E-01,4.621E-01,4.621E-01,4.621E-01,4.621E-01,4.621E-01,4.611E-01,4.611E-01,4.611E-01,4.611E-01,4.611E-01,4.611E-01,4.611E-01,4.611E-01,4.611E-01,4.611E-01,4.611E-01,4.611E-01,4.600E-01,4.600E-01,4.600E-01,4.600E-01,4.600E-01,4.600E-01,4.600E-01,4.600E-01,4.600E-01,4.600E-01,4.600E-01,4.600E-01,4.589E-01,4.589E-01,4.589E-01,4.589E-01,4.589E-01,4.589E-01,4.589E-01,4.589E-01,4.589E-01,4.589E-01,4.589E-01,4.579E-01,4.579E-01,4.579E-01,4.579E-01,4.579E-01,4.579E-01,4.579E-01,4.579E-01,4.579E-01,4.568E-01,4.568E-01,4.568E-01,4.568E-01,4.568E-01,4.568E-01,4.568E-01,4.558E-01,4.558E-01,4.558E-01,4.558E-01,4.558E-01,4.558E-01,4.558E-01,4.547E-01,4.547E-01,4.547E-01,4.547E-01,4.547E-01,4.547E-01,4.547E-01,4.537E-01,4.537E-01,4.537E-01,4.537E-01,4.537E-01,4.537E-01,4.537E-01,4.526E-01,4.526E-01,4.526E-01,4.526E-01,4.526E-01,4.526E-01,4.526E-01,4.516E-01,4.516E-01,4.516E-01,4.516E-01,4.516E-01,4.516E-01,4.505E-01,4.505E-01,4.505E-01,4.505E-01,4.505E-01,4.505E-01,4.505E-01,4.495E-01,4.495E-01,4.495E-01,4.495E-01,4.495E-01,4.495E-01,4.495E-01,4.484E-01,4.484E-01,4.484E-01,4.484E-01,4.484E-01,4.484E-01,4.484E-01,4.474E-01,4.474E-01,4.474E-01,4.474E-01,4.474E-01,4.474E-01,4.474E-01,4.463E-01,4.463E-01,4.463E-01,4.463E-01,4.463E-01,4.463E-01,4.463E-01,4.453E-01,4.453E-01,4.453E-01,4.453E-01,4.453E-01,4.453E-01,4.453E-01,4.442E-01,4.442E-01,4.442E-01,4.442E-01,4.442E-01,4.442E-01,4.442E-01,4.432E-01,4.432E-01,4.432E-01,4.432E-01,4.432E-01,4.432E-01,4.432E-01,4.421E-01,4.421E-01,4.421E-01,4.421E-01,4.421E-01,4.421E-01,4.421E-01,4.411E-01,4.411E-01,4.411E-01,4.411E-01,4.411E-01,4.411E-01,4.411E-01,4.400E-01,4.400E-01,4.400E-01,4.400E-01,4.400E-01,4.400E-01,4.400E-01,4.389E-01,4.389E-01,4.389E-01,4.389E-01,4.389E-01,4.389E-01,4.389E-01,4.379E-01,4.379E-01,4.379E-01,4.379E-01,4.379E-01,4.379E-01,4.379E-01,4.368E-01,4.368E-01,4.368E-01,4.368E-01,4.368E-01,4.368E-01,4.368E-01,4.358E-01,4.358E-01,4.358E-01,4.358E-01,4.358E-01,4.358E-01,4.358E-01,4.358E-01,4.347E-01,4.347E-01,4.347E-01,4.347E-01,4.347E-01,4.347E-01,4.347E-01,4.347E-01,4.337E-01,4.337E-01,4.337E-01,4.337E-01,4.337E-01,4.337E-01,4.337E-01,4.337E-01,4.326E-01,4.326E-01,4.326E-01,4.326E-01,4.326E-01,4.326E-01,4.326E-01,4.316E-01,4.316E-01,4.316E-01,4.316E-01,4.316E-01,4.316E-01,4.316E-01,4.316E-01,4.305E-01,4.305E-01,4.305E-01,4.305E-01,4.305E-01,4.305E-01,4.305E-01,4.305E-01,4.295E-01,4.295E-01,4.295E-01,4.295E-01,4.295E-01,4.295E-01,4.295E-01,4.295E-01,4.284E-01,4.284E-01,4.284E-01,4.284E-01,4.284E-01,4.284E-01,4.284E-01,4.284E-01,4.274E-01,4.274E-01,4.274E-01,4.274E-01,4.274E-01,4.274E-01,4.274E-01,4.274E-01,4.263E-01,4.263E-01,4.263E-01,4.263E-01,4.263E-01,4.263E-01,4.263E-01,4.253E-01,4.253E-01,4.253E-01,4.253E-01,4.253E-01,4.253E-01,4.253E-01,4.242E-01,4.242E-01,4.242E-01,4.242E-01,4.242E-01,4.242E-01,4.242E-01,4.232E-01,4.232E-01,4.232E-01,4.232E-01,4.232E-01,4.232E-01,4.232E-01,4.221E-01,4.221E-01,4.221E-01,4.221E-01,4.221E-01,4.221E-01,4.221E-01,4.211E-01,4.211E-01,4.211E-01,4.211E-01,4.211E-01,4.211E-01,4.211E-01,4.211E-01,4.200E-01,4.200E-01,4.200E-01,4.200E-01,4.200E-01,4.200E-01,4.200E-01,4.189E-01,4.189E-01,4.189E-01,4.189E-01,4.189E-01,4.189E-01,4.189E-01,4.179E-01,4.179E-01,4.179E-01,4.179E-01,4.179E-01,4.179E-01,4.179E-01,4.168E-01,4.168E-01,4.168E-01,4.168E-01,4.168E-01,4.168E-01,4.168E-01,4.168E-01,4.158E-01,4.158E-01,4.158E-01,4.158E-01,4.158E-01,4.158E-01,4.158E-01,4.147E-01,4.147E-01,4.147E-01,4.147E-01,4.147E-01,4.147E-01,4.147E-01,4.137E-01,4.137E-01,4.137E-01,4.137E-01,4.137E-01,4.137E-01,4.137E-01,4.137E-01,4.126E-01,4.126E-01,4.126E-01,4.126E-01,4.126E-01,4.126E-01,4.126E-01,4.116E-01,4.116E-01,4.116E-01,4.116E-01,4.116E-01,4.116E-01,4.116E-01,4.105E-01,4.105E-01,4.105E-01,4.105E-01,4.105E-01,4.105E-01,4.105E-01,4.105E-01,4.095E-01,4.095E-01,4.095E-01,4.095E-01,4.095E-01,4.095E-01,4.095E-01,4.084E-01,4.084E-01,4.084E-01,4.084E-01,4.084E-01,4.084E-01,4.084E-01,4.084E-01,4.074E-01,4.074E-01,4.074E-01,4.074E-01,4.074E-01,4.074E-01,4.074E-01,4.063E-01,4.063E-01,4.063E-01,4.063E-01,4.063E-01,4.063E-01,4.063E-01,4.063E-01,4.053E-01,4.053E-01,4.053E-01,4.053E-01,4.053E-01,4.053E-01,4.053E-01,4.042E-01,4.042E-01,4.042E-01,4.042E-01,4.042E-01,4.042E-01,4.042E-01,4.042E-01,4.032E-01,4.032E-01,4.032E-01,4.032E-01,4.032E-01,4.032E-01,4.032E-01,4.021E-01,4.021E-01,4.021E-01,4.021E-01,4.021E-01,4.021E-01,4.021E-01,4.021E-01,4.011E-01,4.011E-01,4.011E-01,4.011E-01,4.011E-01,4.011E-01,4.011E-01,4.000E-01,4.000E-01,4.000E-01,4.000E-01,4.000E-01,4.000E-01,4.000E-01,4.000E-01,3.989E-01,3.989E-01,3.989E-01,3.989E-01,3.989E-01,3.989E-01,3.989E-01,3.979E-01,3.979E-01,3.979E-01,3.979E-01,3.979E-01,3.979E-01,3.979E-01,3.979E-01,3.968E-01,3.968E-01,3.968E-01,3.968E-01,3.968E-01,3.968E-01,3.968E-01,3.968E-01,3.958E-01,3.958E-01,3.958E-01,3.958E-01,3.958E-01,3.958E-01,3.958E-01,3.947E-01,3.947E-01,3.947E-01,3.947E-01,3.947E-01,3.947E-01,3.947E-01,3.947E-01,3.937E-01,3.937E-01,3.937E-01,3.937E-01,3.937E-01,3.937E-01,3.937E-01,3.937E-01,3.926E-01,3.926E-01,3.926E-01,3.926E-01,3.926E-01,3.926E-01,3.926E-01,3.926E-01,3.916E-01,3.916E-01,3.916E-01,3.916E-01,3.916E-01,3.916E-01,3.916E-01,3.905E-01,3.905E-01,3.905E-01,3.905E-01,3.905E-01,3.905E-01,3.905E-01,3.905E-01,3.895E-01,3.895E-01,3.895E-01,3.895E-01,3.895E-01,3.895E-01,3.895E-01,3.895E-01,3.884E-01,3.884E-01,3.884E-01,3.884E-01,3.884E-01,3.884E-01,3.884E-01,3.884E-01,3.874E-01,3.874E-01,3.874E-01,3.874E-01,3.874E-01,3.874E-01,3.874E-01,3.874E-01,3.863E-01,3.863E-01,3.863E-01,3.863E-01,3.863E-01,3.863E-01,3.863E-01,3.863E-01,3.853E-01,3.853E-01,3.853E-01,3.853E-01,3.853E-01,3.853E-01,3.853E-01,3.853E-01,3.842E-01,3.842E-01,3.842E-01,3.842E-01,3.842E-01,3.842E-01,3.842E-01,3.832E-01,3.832E-01,3.832E-01,3.832E-01,3.832E-01,3.832E-01,3.832E-01,3.832E-01,3.821E-01,3.821E-01,3.821E-01,3.821E-01,3.821E-01,3.821E-01,3.821E-01,3.821E-01,3.811E-01,3.811E-01,3.811E-01,3.811E-01,3.811E-01,3.811E-01,3.811E-01,3.811E-01,3.800E-01,3.800E-01,3.800E-01,3.800E-01,3.800E-01,3.800E-01,3.800E-01,3.800E-01,3.789E-01,3.789E-01};
double dmgElectronX[]={0.1 ,0.11 ,0.13 ,0.14 ,0.15 ,0.17 ,0.18 ,0.2 ,0.22 ,0.24 ,0.26 ,0.29 ,0.32 ,0.35 ,0.38 ,0.42 ,0.46 ,0.5 ,0.55 ,0.6 ,0.66 ,0.72 ,0.79 ,0.87 ,0.95 ,1.05 ,1.15 ,1.26 ,1.38 ,1.51 ,1.66 ,1.82 ,2 ,2.19 ,2.4 ,2.63 ,2.88 ,3.16 ,3.47 ,3.8 ,4.17 ,4.57 ,5.01 ,5.5 ,6.03 ,6.61 ,7.24 ,7.94 ,8.71 ,9.55 ,10.47 ,11.48 ,12.59 ,13.8 ,15.14 ,16.6 ,18.2 ,19.95 ,21.88 ,23.99 ,26.3 ,28.84 ,31.62 ,34.67 ,38.02 ,41.69 ,45.71 ,50.12 ,54.95 ,60.26 ,66.07 ,72.44 ,79.43 ,87.1 ,95.5 ,104.71,114.82,125.89,138.04,151.36,165.96,181.97,199.53,218.78,239.88,263.03,288.4 ,316.23,346.74,380.19,416.87,457.09,501.19,549.54,602.56,660.69,724.44,794.33,870.96,954.99};
double dmgElectronY[]={0,0.00003,0.0001,0.00023,0.00041,0.00064,0.00092,0.00125,0.00163,0.00207,0.00255,0.00309,0.00368,0.00455,0.00542,0.00628,0.00715,0.00803,0.00895,0.00988,0.01081,0.01178,0.01282,0.01387,0.01491,0.01606,0.01732,0.01857,0.01983,0.02109,0.02234,0.0236,0.02488,0.0263,0.02775,0.02919,0.03064,0.03214,0.03368,0.03522,0.03675,0.03829,0.03984,0.04142,0.043,0.04459,0.04616,0.04771,0.04927,0.05082,0.05232,0.05375,0.05519,0.05663,0.05806,0.0595,0.06094,0.06235,0.0636,0.06483,0.06605,0.06728,0.06838,0.0694,0.07042,0.07143,0.07245,0.07344,0.07425,0.07505,0.07584,0.07657,0.07718,0.07779,0.0784,0.07886,0.07917,0.07948,0.0798,0.08011,0.08042,0.08074,0.08098,0.08098,0.08098,0.08098,0.08098,0.08098,0.08098,0.08098,0.08098,0.08098,0.08098,0.08098,0.08098,0.08098,0.08098,0.08098,0.08098,0.08098};
int dmgNeutronN = sizeof(dmgNeutronX)/sizeof(dmgNeutronX[0]);
int dmgProtonN = sizeof(dmgProtonX)/sizeof(dmgProtonX[0]);
int dmgPionN = sizeof(dmgPionX)/sizeof(dmgPionX[0]);
int dmgElectronN = sizeof(dmgElectronX)/sizeof(dmgElectronX[0]);
/*
0=neutron damage
1=proton damage
2=pion damage
3=electron damage
*/
vector<double> x,y;
for(int i=0;i<dmgNeutronN;i++){
x.push_back(dmgNeutronX[i]);
y.push_back(dmgNeutronY[i]);
}
xValNEIL.push_back(x);
yValNEIL.push_back(y);
x.clear();
y.clear();
for(int i=0;i<dmgProtonN;i++){
x.push_back(dmgProtonX[i]);
y.push_back(dmgProtonY[i]);
}
xValNEIL.push_back(x);
yValNEIL.push_back(y);
x.clear();
y.clear();
for(int i=0;i<dmgPionN;i++){
x.push_back(dmgPionX[i]);
y.push_back(dmgPionY[i]);
}
xValNEIL.push_back(x);
yValNEIL.push_back(y);
x.clear();
y.clear();
for(int i=0;i<dmgElectronN;i++){
x.push_back(dmgElectronX[i]);
y.push_back(dmgElectronY[i]);
}
xValNEIL.push_back(x);
yValNEIL.push_back(y);
x.clear();
y.clear();
}
void radDamage::InitMREM(){
double photonERef[]={0.010, 0.015, 0.020, 0.030, 0.040,
0.050, 0.060, 0.080, 0.100, 0.150,
0.200, 0.300, 0.400, 0.500, 0.600,
0.800, 1.000, 1.500, 2.000, 3.000,
4.000, 5.000, 6.000, 8.000, 10.000,
20.000, 30.000, 40.000, 50.000, 60.000,
80.000, 100.000, 200.000, 500.000, 1000.000,
2000.000, 5200.000, 10000.000, 20000.000};
double photonCRef[]={ 360.0, 900.0, 1700.0, 3900.0, 6400.0,
8300.0, 9000.0, 8300.0, 6800.0, 4200.0,
2900.0, 1800.0, 1300.0, 1100.0, 880.0,
680.0, 560.0, 410.0, 330.0, 250.0,
210.0, 180.0, 160.0, 130.0, 110.0,
64.0, 44.0, 34.0, 28.0, 23.0,
17.0, 14.0, 9.2, 5.8, 4.9,
4.3, 3.7, 3.4, 3.2};
double neutronERef[]={.250E-07, .100E-06, .100E-05, .100E-04, .100E-03,
.100E-02, .100E-01, .100E+00, .500E+00, .100E+01, .250E+01,
.500E+01, .700E+01, .100E+02, .140E+02, .200E+02, .400E+02,
.600E+02, .100E+03, .200E+03, .300E+03, .400E+03};
double neutronCRef[]={.1021E-05, .1021E-05, .1240E-05, .1240E-05, .1197E-05,
.1021E-05, .9921E-06, .6039E-05, .2572E-04, .3655E-04, .3472E-04,
.4340E-04, .4085E-04, .4085E-04, .5787E-04, .6313E-04, .6944E-04,
.6313E-04, .4960E-04, .5342E-04, .6313E-04, .6944E-04};
double electronERef[]={ 100.000, 200.000, 500.000, 1000.000,
5000.000, 10000.000, 20000.000};
double electronCRef[]={ 6.2, 5.3, 3.8, 3.0,
2.0, 1.8, 1.5};
int neutronNr = sizeof(neutronERef)/sizeof(neutronERef[0]);
int photonNr = sizeof(photonERef)/sizeof(photonERef[0]);
int electronNr = sizeof(electronERef)/sizeof(electronERef[0]);
//0=neutrons, 1=photons, 2=e/p
vector<double> x,y;
for(int i=0;i<neutronNr;i++){
x.push_back(neutronERef[i]);
y.push_back(neutronCRef[i]);
}
xValMREM.push_back(x);
yValMREM.push_back(y);
x.clear();
y.clear();
for(int i=0;i<photonNr;i++){
x.push_back(photonERef[i]);
y.push_back(photonCRef[i]);
}
xValMREM.push_back(x);
yValMREM.push_back(y);
x.clear();
y.clear();
for(int i=0;i<electronNr;i++){
x.push_back(electronERef[i]);
y.push_back(electronCRef[i]);
}
xValMREM.push_back(x);
yValMREM.push_back(y);
x.clear();
y.clear();
}