Belle II Software development
Modules | Classes | Functions | Variables
simulation
Packages

Modules

 simulation data objects
 
 simulation modules
 

Classes

class  BkgNeutronWeight
 The class to get the weighting factor for a 1-MeV-equivalent neutron flux on silicon. More...
 
class  BkgSensitiveDetector
 The Class for BeamBackground Sensitive Detector. More...
 
class  ROISimulationParameters
 The payload containing all PXD ROI parameters. More...
 
class  DetectorConstruction
 Class responsible to connect to geometry to simulation. More...
 
class  RandomEngine
 Interface class to make Geant4 use the Belle2 RandomGenerator. More...
 
class  G4LongLivedNeutral
 A class to hold long-lived neutral particle description. More...
 
class  G4LongLivedNeutralDecay
 This class is a decay process. More...
 
class  G4LongLivedNeutralPhysics
 LongLivedNeutral physics Class – to be registered in the physics list. More...
 
class  G4LongLivedNeutralTransportation
 Concrete class that does the geometrical transport. More...
 

Functions

static BkgNeutronWeight & getInstance ()
 Return a reference to the singleton BkgNeutronWeight instance.
 
double getWeight (double ke)
 Get weighting factor to convert a neutron to its 1-MeV equivalent.
 
 BkgNeutronWeight (void)
 Constructor is private for this singleton.
 
void addPoint (double, double)
 Function to add one (x,y) data point to the spline table.
 
void doSpline ()
 Function to calculate the cubic spline derivatives.
 
 BkgSensitiveDetector (const char *subDet, int iden=0)
 Constructor.
 
bool step (G4Step *aStep, G4TouchableHistory *) override
 Process each step and calculate variables defined in PXDSimHit.
 
virtual G4VParticleChange * AtRestDoIt (const G4Track &aTrack, const G4Step &aStep) override
 G4VProcess::AtRestDoIt() implementation for decay at rest.
 
G4double GetRemainderLifeTime () const
 Get Remainder of life time at rest decay.
 

Variables

static BkgNeutronWeight * m_Instance = 0
 static pointer to the singleton instance of this class
 

Detailed Description

Function Documentation

◆ addPoint()

void addPoint ( double  x,
double  y 
)
private

Function to add one (x,y) data point to the spline table.

Definition at line 1429 of file BkgNeutronWeight.cc.

1430 {
1431 m_x.push_back(x);
1432 m_y.push_back(y);
1433 m_b.push_back(0.0);
1434 m_c.push_back(0.0);
1435 m_d.push_back(0.0);
1436 }
Belle2::BkgNeutronWeight::m_b
std::vector< double > m_b
Vector of first derivatives in spline fit (calculated)
Definition: BkgNeutronWeight.h:54
Belle2::BkgNeutronWeight::m_c
std::vector< double > m_c
Vector of second derivatives in spline fit (calculated)
Definition: BkgNeutronWeight.h:57
Belle2::BkgNeutronWeight::m_y
std::vector< double > m_y
Vector of tabulated y coordinates (weighting factors)
Definition: BkgNeutronWeight.h:51
Belle2::BkgNeutronWeight::m_d
std::vector< double > m_d
Vector of third derivatives in spline fit (calculated)
Definition: BkgNeutronWeight.h:60
Belle2::BkgNeutronWeight::m_x
std::vector< double > m_x
Vector of tabulated x coordinates (logarithm of kinetic energy)
Definition: BkgNeutronWeight.h:48

◆ AtRestDoIt()

G4VParticleChange * AtRestDoIt ( const G4Track &  aTrack,
const G4Step &  aStep 
)
inlineoverridevirtual

G4VProcess::AtRestDoIt() implementation for decay at rest.

Proposes particle change by DecayIt method at rest

Parameters
aTrackPropagating particle track reference
aStepCurrent step reference

Definition at line 185 of file G4LongLivedNeutralDecay.h.

189 {
190 return DecayIt(aTrack, aStep);
191 }
Belle2::G4LongLivedNeutralDecay::DecayIt
virtual G4VParticleChange * DecayIt(const G4Track &aTrack, const G4Step &aStep)
G4Decay::DecayIt() implemention.
Definition: G4LongLivedNeutralDecay.cc:88

◆ BkgNeutronWeight()

BkgNeutronWeight ( void  )
private

Constructor is private for this singleton.

Definition at line 43 of file BkgNeutronWeight.cc.

44 {
45 addPoint(log(1.025E-10), 1.575E-02);
46 addPoint(log(1.075E-10), 1.537E-02);
47 addPoint(log(1.125E-10), 1.503E-02);
48 addPoint(log(1.175E-10), 1.470E-02);
49 addPoint(log(1.238E-10), 1.432E-02);
50 addPoint(log(1.313E-10), 1.391E-02);
51 addPoint(log(1.388E-10), 1.353E-02);
52 addPoint(log(1.463E-10), 1.317E-02);
53 addPoint(log(1.550E-10), 1.280E-02);
54 addPoint(log(1.650E-10), 1.242E-02);
55 addPoint(log(1.750E-10), 1.206E-02);
56 addPoint(log(1.850E-10), 1.172E-02);
57 addPoint(log(1.950E-10), 1.142E-02);
58 addPoint(log(2.050E-10), 1.113E-02);
59 addPoint(log(2.150E-10), 1.087E-02);
60 addPoint(log(2.250E-10), 1.063E-02);
61 addPoint(log(2.350E-10), 1.039E-02);
62 addPoint(log(2.475E-10), 1.013E-02);
63 addPoint(log(2.625E-10), 9.834E-03);
64 addPoint(log(2.750E-10), 9.609E-03);
65 addPoint(log(2.900E-10), 9.356E-03);
66 addPoint(log(3.100E-10), 9.048E-03);
67 addPoint(log(3.300E-10), 8.770E-03);
68 addPoint(log(3.500E-10), 8.516E-03);
69 addPoint(log(3.700E-10), 8.290E-03);
70 addPoint(log(3.900E-10), 8.075E-03);
71 addPoint(log(4.125E-10), 7.847E-03);
72 addPoint(log(4.375E-10), 7.622E-03);
73 addPoint(log(4.625E-10), 7.410E-03);
74 addPoint(log(4.875E-10), 7.219E-03);
75 addPoint(log(5.125E-10), 7.039E-03);
76 addPoint(log(5.375E-10), 6.875E-03);
77 addPoint(log(5.625E-10), 6.718E-03);
78 addPoint(log(5.875E-10), 6.575E-03);
79 addPoint(log(6.150E-10), 6.425E-03);
80 addPoint(log(6.450E-10), 6.274E-03);
81 addPoint(log(6.750E-10), 6.134E-03);
82 addPoint(log(7.050E-10), 6.002E-03);
83 addPoint(log(7.400E-10), 5.864E-03);
84 addPoint(log(7.800E-10), 5.709E-03);
85 addPoint(log(8.200E-10), 5.567E-03);
86 addPoint(log(8.600E-10), 5.438E-03);
87 addPoint(log(9.000E-10), 5.312E-03);
88 addPoint(log(9.400E-10), 5.200E-03);
89 addPoint(log(9.800E-10), 5.091E-03);
90 addPoint(log(1.025E-09), 4.978E-03);
91 addPoint(log(1.075E-09), 4.861E-03);
92 addPoint(log(1.125E-09), 4.751E-03);
93 addPoint(log(1.175E-09), 4.649E-03);
94 addPoint(log(1.238E-09), 4.530E-03);
95 addPoint(log(1.313E-09), 4.398E-03);
96 addPoint(log(1.388E-09), 4.278E-03);
97 addPoint(log(1.463E-09), 4.171E-03);
98 addPoint(log(1.550E-09), 4.049E-03);
99 addPoint(log(1.650E-09), 3.925E-03);
100 addPoint(log(1.750E-09), 3.810E-03);
101 addPoint(log(1.850E-09), 3.706E-03);
102 addPoint(log(1.950E-09), 3.609E-03);
103 addPoint(log(2.050E-09), 3.520E-03);
104 addPoint(log(2.150E-09), 3.437E-03);
105 addPoint(log(2.250E-09), 3.360E-03);
106 addPoint(log(2.350E-09), 3.287E-03);
107 addPoint(log(2.475E-09), 3.203E-03);
108 addPoint(log(2.625E-09), 3.110E-03);
109 addPoint(log(2.750E-09), 3.039E-03);
110 addPoint(log(2.900E-09), 2.962E-03);
111 addPoint(log(3.100E-09), 2.864E-03);
112 addPoint(log(3.300E-09), 2.776E-03);
113 addPoint(log(3.500E-09), 2.695E-03);
114 addPoint(log(3.700E-09), 2.621E-03);
115 addPoint(log(3.900E-09), 2.552E-03);
116 addPoint(log(4.125E-09), 2.482E-03);
117 addPoint(log(4.375E-09), 2.409E-03);
118 addPoint(log(4.625E-09), 2.344E-03);
119 addPoint(log(4.875E-09), 2.283E-03);
120 addPoint(log(5.125E-09), 2.226E-03);
121 addPoint(log(5.375E-09), 2.174E-03);
122 addPoint(log(5.625E-09), 2.125E-03);
123 addPoint(log(5.875E-09), 2.082E-03);
124 addPoint(log(6.150E-09), 2.034E-03);
125 addPoint(log(6.450E-09), 1.985E-03);
126 addPoint(log(6.750E-09), 1.941E-03);
127 addPoint(log(7.050E-09), 1.898E-03);
128 addPoint(log(7.400E-09), 1.853E-03);
129 addPoint(log(7.800E-09), 1.805E-03);
130 addPoint(log(8.200E-09), 1.760E-03);
131 addPoint(log(8.600E-09), 1.718E-03);
132 addPoint(log(9.000E-09), 1.680E-03);
133 addPoint(log(9.400E-09), 1.644E-03);
134 addPoint(log(9.800E-09), 1.610E-03);
135 addPoint(log(1.025E-08), 1.574E-03);
136 addPoint(log(1.075E-08), 1.537E-03);
137 addPoint(log(1.125E-08), 1.503E-03);
138 addPoint(log(1.175E-08), 1.472E-03);
139 addPoint(log(1.238E-08), 1.433E-03);
140 addPoint(log(1.313E-08), 1.392E-03);
141 addPoint(log(1.388E-08), 1.353E-03);
142 addPoint(log(1.463E-08), 1.318E-03);
143 addPoint(log(1.550E-08), 1.280E-03);
144 addPoint(log(1.650E-08), 1.241E-03);
145 addPoint(log(1.750E-08), 1.204E-03);
146 addPoint(log(1.850E-08), 1.172E-03);
147 addPoint(log(1.950E-08), 1.141E-03);
148 addPoint(log(2.050E-08), 1.113E-03);
149 addPoint(log(2.150E-08), 1.087E-03);
150 addPoint(log(2.250E-08), 1.062E-03);
151 addPoint(log(2.350E-08), 1.040E-03);
152 addPoint(log(2.475E-08), 1.013E-03);
153 addPoint(log(2.625E-08), 9.841E-04);
154 addPoint(log(2.750E-08), 9.619E-04);
155 addPoint(log(2.900E-08), 9.359E-04);
156 addPoint(log(3.100E-08), 9.058E-04);
157 addPoint(log(3.300E-08), 8.773E-04);
158 addPoint(log(3.500E-08), 8.523E-04);
159 addPoint(log(3.700E-08), 8.285E-04);
160 addPoint(log(3.900E-08), 8.073E-04);
161 addPoint(log(4.125E-08), 7.847E-04);
162 addPoint(log(4.375E-08), 7.619E-04);
163 addPoint(log(4.625E-08), 7.410E-04);
164 addPoint(log(4.875E-08), 7.215E-04);
165 addPoint(log(5.125E-08), 7.036E-04);
166 addPoint(log(5.375E-08), 6.871E-04);
167 addPoint(log(5.625E-08), 6.716E-04);
168 addPoint(log(5.875E-08), 6.578E-04);
169 addPoint(log(6.150E-08), 6.428E-04);
170 addPoint(log(6.450E-08), 6.273E-04);
171 addPoint(log(6.750E-08), 6.136E-04);
172 addPoint(log(7.050E-08), 6.001E-04);
173 addPoint(log(7.400E-08), 5.861E-04);
174 addPoint(log(7.800E-08), 5.707E-04);
175 addPoint(log(8.200E-08), 5.567E-04);
176 addPoint(log(8.600E-08), 5.434E-04);
177 addPoint(log(9.000E-08), 5.313E-04);
178 addPoint(log(9.400E-08), 5.197E-04);
179 addPoint(log(9.800E-08), 5.090E-04);
180 addPoint(log(1.025E-07), 4.976E-04);
181 addPoint(log(1.075E-07), 4.864E-04);
182 addPoint(log(1.125E-07), 4.753E-04);
183 addPoint(log(1.175E-07), 4.646E-04);
184 addPoint(log(1.238E-07), 4.531E-04);
185 addPoint(log(1.313E-07), 4.397E-04);
186 addPoint(log(1.388E-07), 4.279E-04);
187 addPoint(log(1.463E-07), 4.166E-04);
188 addPoint(log(1.550E-07), 4.049E-04);
189 addPoint(log(1.650E-07), 3.924E-04);
190 addPoint(log(1.750E-07), 3.810E-04);
191 addPoint(log(1.850E-07), 3.705E-04);
192 addPoint(log(1.950E-07), 3.607E-04);
193 addPoint(log(2.050E-07), 3.518E-04);
194 addPoint(log(2.150E-07), 3.437E-04);
195 addPoint(log(2.250E-07), 3.362E-04);
196 addPoint(log(2.350E-07), 3.286E-04);
197 addPoint(log(2.475E-07), 3.204E-04);
198 addPoint(log(2.625E-07), 3.111E-04);
199 addPoint(log(2.750E-07), 3.039E-04);
200 addPoint(log(2.900E-07), 2.962E-04);
201 addPoint(log(3.100E-07), 2.865E-04);
202 addPoint(log(3.300E-07), 2.775E-04);
203 addPoint(log(3.500E-07), 2.695E-04);
204 addPoint(log(3.700E-07), 2.620E-04);
205 addPoint(log(3.900E-07), 2.551E-04);
206 addPoint(log(4.125E-07), 2.483E-04);
207 addPoint(log(4.375E-07), 2.410E-04);
208 addPoint(log(4.625E-07), 2.343E-04);
209 addPoint(log(4.875E-07), 2.284E-04);
210 addPoint(log(5.125E-07), 2.226E-04);
211 addPoint(log(5.375E-07), 2.176E-04);
212 addPoint(log(5.625E-07), 2.128E-04);
213 addPoint(log(5.875E-07), 2.082E-04);
214 addPoint(log(6.150E-07), 2.034E-04);
215 addPoint(log(6.450E-07), 1.985E-04);
216 addPoint(log(6.750E-07), 1.941E-04);
217 addPoint(log(7.050E-07), 1.898E-04);
218 addPoint(log(7.400E-07), 1.852E-04);
219 addPoint(log(7.800E-07), 1.805E-04);
220 addPoint(log(8.200E-07), 1.762E-04);
221 addPoint(log(8.600E-07), 1.718E-04);
222 addPoint(log(9.000E-07), 1.679E-04);
223 addPoint(log(9.400E-07), 1.645E-04);
224 addPoint(log(9.800E-07), 1.610E-04);
225 addPoint(log(1.025E-06), 1.575E-04);
226 addPoint(log(1.075E-06), 1.542E-04);
227 addPoint(log(1.125E-06), 1.508E-04);
228 addPoint(log(1.175E-06), 1.474E-04);
229 addPoint(log(1.238E-06), 1.434E-04);
230 addPoint(log(1.313E-06), 1.393E-04);
231 addPoint(log(1.388E-06), 1.356E-04);
232 addPoint(log(1.463E-06), 1.318E-04);
233 addPoint(log(1.550E-06), 1.281E-04);
234 addPoint(log(1.650E-06), 1.242E-04);
235 addPoint(log(1.750E-06), 1.204E-04);
236 addPoint(log(1.850E-06), 1.172E-04);
237 addPoint(log(1.950E-06), 1.141E-04);
238 addPoint(log(2.050E-06), 1.114E-04);
239 addPoint(log(2.150E-06), 1.090E-04);
240 addPoint(log(2.250E-06), 1.066E-04);
241 addPoint(log(2.350E-06), 1.042E-04);
242 addPoint(log(2.475E-06), 1.013E-04);
243 addPoint(log(2.625E-06), 9.853E-05);
244 addPoint(log(2.750E-06), 9.627E-05);
245 addPoint(log(2.900E-06), 9.360E-05);
246 addPoint(log(3.100E-06), 9.056E-05);
247 addPoint(log(3.300E-06), 8.776E-05);
248 addPoint(log(3.500E-06), 8.513E-05);
249 addPoint(log(3.700E-06), 8.287E-05);
250 addPoint(log(3.900E-06), 8.065E-05);
251 addPoint(log(4.125E-06), 7.860E-05);
252 addPoint(log(4.375E-06), 7.644E-05);
253 addPoint(log(4.625E-06), 7.428E-05);
254 addPoint(log(4.875E-06), 7.217E-05);
255 addPoint(log(5.125E-06), 7.045E-05);
256 addPoint(log(5.375E-06), 6.885E-05);
257 addPoint(log(5.625E-06), 6.725E-05);
258 addPoint(log(5.875E-06), 6.569E-05);
259 addPoint(log(6.150E-06), 6.427E-05);
260 addPoint(log(6.450E-06), 6.278E-05);
261 addPoint(log(6.750E-06), 6.129E-05);
262 addPoint(log(7.050E-06), 5.998E-05);
263 addPoint(log(7.400E-06), 5.858E-05);
264 addPoint(log(7.800E-06), 5.701E-05);
265 addPoint(log(8.200E-06), 5.572E-05);
266 addPoint(log(8.600E-06), 5.450E-05);
267 addPoint(log(9.000E-06), 5.327E-05);
268 addPoint(log(9.400E-06), 5.205E-05);
269 addPoint(log(9.800E-06), 5.087E-05);
270 addPoint(log(1.025E-05), 4.981E-05);
271 addPoint(log(1.075E-05), 4.868E-05);
272 addPoint(log(1.125E-05), 4.754E-05);
273 addPoint(log(1.175E-05), 4.644E-05);
274 addPoint(log(1.238E-05), 4.531E-05);
275 addPoint(log(1.313E-05), 4.399E-05);
276 addPoint(log(1.388E-05), 4.274E-05);
277 addPoint(log(1.463E-05), 4.167E-05);
278 addPoint(log(1.550E-05), 4.045E-05);
279 addPoint(log(1.650E-05), 3.929E-05);
280 addPoint(log(1.750E-05), 3.821E-05);
281 addPoint(log(1.850E-05), 3.712E-05);
282 addPoint(log(1.950E-05), 3.607E-05);
283 addPoint(log(2.050E-05), 3.522E-05);
284 addPoint(log(2.150E-05), 3.442E-05);
285 addPoint(log(2.250E-05), 3.361E-05);
286 addPoint(log(2.350E-05), 3.284E-05);
287 addPoint(log(2.475E-05), 3.204E-05);
288 addPoint(log(2.625E-05), 3.110E-05);
289 addPoint(log(2.750E-05), 3.035E-05);
290 addPoint(log(2.900E-05), 2.959E-05);
291 addPoint(log(3.100E-05), 2.860E-05);
292 addPoint(log(3.300E-05), 2.778E-05);
293 addPoint(log(3.500E-05), 2.702E-05);
294 addPoint(log(3.700E-05), 2.625E-05);
295 addPoint(log(3.900E-05), 2.550E-05);
296 addPoint(log(4.125E-05), 2.483E-05);
297 addPoint(log(4.375E-05), 2.412E-05);
298 addPoint(log(4.625E-05), 2.342E-05);
299 addPoint(log(4.875E-05), 2.282E-05);
300 addPoint(log(5.125E-05), 2.227E-05);
301 addPoint(log(5.375E-05), 2.172E-05);
302 addPoint(log(5.625E-05), 2.123E-05);
303 addPoint(log(5.875E-05), 2.079E-05);
304 addPoint(log(6.150E-05), 2.030E-05);
305 addPoint(log(6.450E-05), 1.985E-05);
306 addPoint(log(6.750E-05), 1.944E-05);
307 addPoint(log(7.050E-05), 1.904E-05);
308 addPoint(log(7.400E-05), 1.856E-05);
309 addPoint(log(7.800E-05), 1.803E-05);
310 addPoint(log(8.200E-05), 1.761E-05);
311 addPoint(log(8.600E-05), 1.721E-05);
312 addPoint(log(9.000E-05), 1.681E-05);
313 addPoint(log(9.400E-05), 1.642E-05);
314 addPoint(log(9.800E-05), 1.610E-05);
315 addPoint(log(1.025E-04), 1.575E-05);
316 addPoint(log(1.075E-04), 1.536E-05);
317 addPoint(log(1.125E-04), 1.502E-05);
318 addPoint(log(1.175E-04), 1.470E-05);
319 addPoint(log(1.238E-04), 1.432E-05);
320 addPoint(log(1.313E-04), 1.393E-05);
321 addPoint(log(1.388E-04), 1.357E-05);
322 addPoint(log(1.463E-04), 1.321E-05);
323 addPoint(log(1.550E-04), 1.714E-05);
324 addPoint(log(1.650E-04), 6.792E-05);
325 addPoint(log(1.750E-04), 1.315E-04);
326 addPoint(log(1.850E-04), 1.938E-04);
327 addPoint(log(1.950E-04), 2.250E-04);
328 addPoint(log(2.050E-04), 2.433E-04);
329 addPoint(log(2.150E-04), 2.614E-04);
330 addPoint(log(2.250E-04), 2.754E-04);
331 addPoint(log(2.350E-04), 2.867E-04);
332 addPoint(log(2.475E-04), 3.007E-04);
333 addPoint(log(2.625E-04), 3.177E-04);
334 addPoint(log(2.750E-04), 3.318E-04);
335 addPoint(log(2.900E-04), 3.486E-04);
336 addPoint(log(3.100E-04), 3.712E-04);
337 addPoint(log(3.300E-04), 3.937E-04);
338 addPoint(log(3.500E-04), 4.163E-04);
339 addPoint(log(3.700E-04), 4.388E-04);
340 addPoint(log(3.900E-04), 4.613E-04);
341 addPoint(log(4.125E-04), 4.865E-04);
342 addPoint(log(4.375E-04), 5.146E-04);
343 addPoint(log(4.625E-04), 5.426E-04);
344 addPoint(log(4.875E-04), 5.706E-04);
345 addPoint(log(5.125E-04), 5.985E-04);
346 addPoint(log(5.375E-04), 6.264E-04);
347 addPoint(log(5.625E-04), 6.543E-04);
348 addPoint(log(5.875E-04), 6.822E-04);
349 addPoint(log(6.150E-04), 7.128E-04);
350 addPoint(log(6.450E-04), 7.462E-04);
351 addPoint(log(6.750E-04), 7.795E-04);
352 addPoint(log(7.050E-04), 8.128E-04);
353 addPoint(log(7.400E-04), 8.515E-04);
354 addPoint(log(7.800E-04), 8.957E-04);
355 addPoint(log(8.200E-04), 9.399E-04);
356 addPoint(log(8.600E-04), 9.840E-04);
357 addPoint(log(9.000E-04), 1.028E-03);
358 addPoint(log(9.400E-04), 1.072E-03);
359 addPoint(log(9.800E-04), 1.116E-03);
360 addPoint(log(1.025E-03), 1.165E-03);
361 addPoint(log(1.075E-03), 1.220E-03);
362 addPoint(log(1.125E-03), 1.274E-03);
363 addPoint(log(1.175E-03), 1.329E-03);
364 addPoint(log(1.238E-03), 1.397E-03);
365 addPoint(log(1.313E-03), 1.478E-03);
366 addPoint(log(1.388E-03), 1.559E-03);
367 addPoint(log(1.463E-03), 1.640E-03);
368 addPoint(log(1.550E-03), 1.734E-03);
369 addPoint(log(1.650E-03), 1.842E-03);
370 addPoint(log(1.750E-03), 1.949E-03);
371 addPoint(log(1.850E-03), 2.056E-03);
372 addPoint(log(1.950E-03), 2.162E-03);
373 addPoint(log(2.050E-03), 2.269E-03);
374 addPoint(log(2.150E-03), 2.382E-03);
375 addPoint(log(2.250E-03), 4.405E-03);
376 addPoint(log(2.350E-03), 2.589E-03);
377 addPoint(log(2.475E-03), 2.717E-03);
378 addPoint(log(2.625E-03), 2.874E-03);
379 addPoint(log(2.750E-03), 3.004E-03);
380 addPoint(log(2.900E-03), 3.159E-03);
381 addPoint(log(3.100E-03), 3.366E-03);
382 addPoint(log(3.300E-03), 3.573E-03);
383 addPoint(log(3.500E-03), 3.778E-03);
384 addPoint(log(3.700E-03), 3.982E-03);
385 addPoint(log(3.900E-03), 4.186E-03);
386 addPoint(log(4.125E-03), 4.413E-03);
387 addPoint(log(4.375E-03), 4.666E-03);
388 addPoint(log(4.625E-03), 4.917E-03);
389 addPoint(log(4.875E-03), 5.492E-03);
390 addPoint(log(5.125E-03), 5.419E-03);
391 addPoint(log(5.375E-03), 5.662E-03);
392 addPoint(log(5.625E-03), 5.907E-03);
393 addPoint(log(5.875E-03), 6.151E-03);
394 addPoint(log(6.150E-03), 6.418E-03);
395 addPoint(log(6.450E-03), 6.710E-03);
396 addPoint(log(6.750E-03), 6.999E-03);
397 addPoint(log(7.050E-03), 7.285E-03);
398 addPoint(log(7.400E-03), 7.618E-03);
399 addPoint(log(7.800E-03), 7.998E-03);
400 addPoint(log(8.200E-03), 8.372E-03);
401 addPoint(log(8.600E-03), 8.745E-03);
402 addPoint(log(9.000E-03), 9.115E-03);
403 addPoint(log(9.400E-03), 9.477E-03);
404 addPoint(log(9.800E-03), 9.838E-03);
405 addPoint(log(1.025E-02), 1.024E-02);
406 addPoint(log(1.075E-02), 1.069E-02);
407 addPoint(log(1.125E-02), 1.115E-02);
408 addPoint(log(1.175E-02), 1.159E-02);
409 addPoint(log(1.238E-02), 1.213E-02);
410 addPoint(log(1.313E-02), 1.277E-02);
411 addPoint(log(1.388E-02), 1.340E-02);
412 addPoint(log(1.463E-02), 1.403E-02);
413 addPoint(log(1.550E-02), 1.557E-02);
414 addPoint(log(1.650E-02), 1.554E-02);
415 addPoint(log(1.750E-02), 1.634E-02);
416 addPoint(log(1.850E-02), 1.710E-02);
417 addPoint(log(1.950E-02), 1.782E-02);
418 addPoint(log(2.050E-02), 1.854E-02);
419 addPoint(log(2.150E-02), 1.926E-02);
420 addPoint(log(2.250E-02), 1.998E-02);
421 addPoint(log(2.350E-02), 2.069E-02);
422 addPoint(log(2.475E-02), 2.150E-02);
423 addPoint(log(2.625E-02), 2.245E-02);
424 addPoint(log(2.750E-02), 2.317E-02);
425 addPoint(log(2.900E-02), 2.400E-02);
426 addPoint(log(3.100E-02), 2.509E-02);
427 addPoint(log(3.300E-02), 2.599E-02);
428 addPoint(log(3.500E-02), 2.675E-02);
429 addPoint(log(3.700E-02), 2.743E-02);
430 addPoint(log(3.900E-02), 3.247E-02);
431 addPoint(log(4.125E-02), 2.792E-02);
432 addPoint(log(4.375E-02), 2.739E-02);
433 addPoint(log(4.625E-02), 2.598E-02);
434 addPoint(log(4.875E-02), 2.288E-02);
435 addPoint(log(5.125E-02), 1.626E-02);
436 addPoint(log(5.375E-02), 1.485E-02);
437 addPoint(log(5.625E-02), 5.078E-01);
438 addPoint(log(5.875E-02), 1.145E-01);
439 addPoint(log(6.150E-02), 7.793E-02);
440 addPoint(log(6.450E-02), 6.706E-02);
441 addPoint(log(6.750E-02), 6.455E-02);
442 addPoint(log(7.050E-02), 5.988E-02);
443 addPoint(log(7.400E-02), 5.797E-02);
444 addPoint(log(7.800E-02), 5.636E-02);
445 addPoint(log(8.200E-02), 5.494E-02);
446 addPoint(log(8.600E-02), 5.374E-02);
447 addPoint(log(9.000E-02), 5.188E-02);
448 addPoint(log(9.400E-02), 4.992E-02);
449 addPoint(log(9.800E-02), 4.790E-02);
450 addPoint(log(1.025E-01), 4.526E-02);
451 addPoint(log(1.075E-01), 4.162E-02);
452 addPoint(log(1.125E-01), 3.764E-02);
453 addPoint(log(1.175E-01), 3.297E-02);
454 addPoint(log(1.238E-01), 2.645E-02);
455 addPoint(log(1.313E-01), 1.826E-02);
456 addPoint(log(1.388E-01), 1.121E-02);
457 addPoint(log(1.463E-01), 1.090E-02);
458 addPoint(log(1.550E-01), 4.848E-02);
459 addPoint(log(1.650E-01), 2.173E-01);
460 addPoint(log(1.750E-01), 6.987E-01);
461 addPoint(log(1.850E-01), 1.196E+00);
462 addPoint(log(1.950E-01), 1.162E+00);
463 addPoint(log(2.050E-01), 9.710E-01);
464 addPoint(log(2.150E-01), 8.288E-01);
465 addPoint(log(2.250E-01), 7.334E-01);
466 addPoint(log(2.350E-01), 6.709E-01);
467 addPoint(log(2.475E-01), 6.172E-01);
468 addPoint(log(2.625E-01), 5.768E-01);
469 addPoint(log(2.750E-01), 5.547E-01);
470 addPoint(log(2.900E-01), 5.393E-01);
471 addPoint(log(3.100E-01), 5.237E-01);
472 addPoint(log(3.300E-01), 5.198E-01);
473 addPoint(log(3.500E-01), 5.097E-01);
474 addPoint(log(3.700E-01), 5.036E-01);
475 addPoint(log(3.900E-01), 5.320E-01);
476 addPoint(log(4.125E-01), 5.403E-01);
477 addPoint(log(4.375E-01), 5.451E-01);
478 addPoint(log(4.625E-01), 5.581E-01);
479 addPoint(log(4.875E-01), 5.740E-01);
480 addPoint(log(5.125E-01), 5.987E-01);
481 addPoint(log(5.375E-01), 7.554E-01);
482 addPoint(log(5.625E-01), 1.280E+00);
483 addPoint(log(5.875E-01), 5.997E-01);
484 addPoint(log(6.150E-01), 5.425E-01);
485 addPoint(log(6.450E-01), 5.536E-01);
486 addPoint(log(6.750E-01), 5.786E-01);
487 addPoint(log(7.050E-01), 6.019E-01);
488 addPoint(log(7.400E-01), 6.724E-01);
489 addPoint(log(7.800E-01), 9.209E-01);
490 addPoint(log(8.200E-01), 1.459E+00);
491 addPoint(log(8.600E-01), 8.318E-01);
492 addPoint(log(9.000E-01), 9.434E-01);
493 addPoint(log(9.400E-01), 1.172E+00);
494 addPoint(log(9.800E-01), 1.174E+00);
495 addPoint(log(1.050E+00), 8.020E-01);
496 addPoint(log(1.150E+00), 6.578E-01);
497 addPoint(log(1.250E+00), 9.680E-01);
498 addPoint(log(1.350E+00), 9.410E-01);
499 addPoint(log(1.450E+00), 1.079E+00);
500 addPoint(log(1.550E+00), 1.128E+00);
501 addPoint(log(1.650E+00), 1.766E+00);
502 addPoint(log(1.750E+00), 8.366E-01);
503 addPoint(log(1.850E+00), 1.411E+00);
504 addPoint(log(1.950E+00), 1.393E+00);
505 addPoint(log(2.050E+00), 1.022E+00);
506 addPoint(log(2.150E+00), 1.159E+00);
507 addPoint(log(2.250E+00), 1.126E+00);
508 addPoint(log(2.350E+00), 1.106E+00);
509 addPoint(log(2.450E+00), 1.267E+00);
510 addPoint(log(2.550E+00), 1.384E+00);
511 addPoint(log(2.650E+00), 1.238E+00);
512 addPoint(log(2.750E+00), 1.153E+00);
513 addPoint(log(2.850E+00), 1.437E+00);
514 addPoint(log(2.950E+00), 1.061E+00);
515 addPoint(log(3.050E+00), 1.278E+00);
516 addPoint(log(3.150E+00), 1.427E+00);
517 addPoint(log(3.250E+00), 1.281E+00);
518 addPoint(log(3.350E+00), 1.217E+00);
519 addPoint(log(3.450E+00), 1.259E+00);
520 addPoint(log(3.550E+00), 1.213E+00);
521 addPoint(log(3.650E+00), 7.506E-01);
522 addPoint(log(3.750E+00), 1.235E+00);
523 addPoint(log(3.850E+00), 1.187E+00);
524 addPoint(log(3.950E+00), 1.452E+00);
525 addPoint(log(4.050E+00), 1.431E+00);
526 addPoint(log(4.150E+00), 1.147E+00);
527 addPoint(log(4.250E+00), 1.782E+00);
528 addPoint(log(4.350E+00), 1.442E+00);
529 addPoint(log(4.450E+00), 1.509E+00);
530 addPoint(log(4.550E+00), 1.498E+00);
531 addPoint(log(4.650E+00), 1.673E+00);
532 addPoint(log(4.750E+00), 2.003E+00);
533 addPoint(log(4.850E+00), 1.717E+00);
534 addPoint(log(4.950E+00), 1.579E+00);
535 addPoint(log(5.050E+00), 1.612E+00);
536 addPoint(log(5.150E+00), 1.831E+00);
537 addPoint(log(5.250E+00), 1.605E+00);
538 addPoint(log(5.350E+00), 1.336E+00);
539 addPoint(log(5.450E+00), 1.277E+00);
540 addPoint(log(5.550E+00), 1.567E+00);
541 addPoint(log(5.650E+00), 1.616E+00);
542 addPoint(log(5.750E+00), 1.929E+00);
543 addPoint(log(5.850E+00), 1.794E+00);
544 addPoint(log(5.950E+00), 1.468E+00);
545 addPoint(log(6.050E+00), 1.666E+00);
546 addPoint(log(6.150E+00), 1.364E+00);
547 addPoint(log(6.250E+00), 1.900E+00);
548 addPoint(log(6.350E+00), 1.641E+00);
549 addPoint(log(6.450E+00), 1.514E+00);
550 addPoint(log(6.550E+00), 1.314E+00);
551 addPoint(log(6.650E+00), 1.611E+00);
552 addPoint(log(6.750E+00), 1.778E+00);
553 addPoint(log(6.850E+00), 1.577E+00);
554 addPoint(log(6.950E+00), 1.521E+00);
555 addPoint(log(7.050E+00), 1.770E+00);
556 addPoint(log(7.150E+00), 1.470E+00);
557 addPoint(log(7.250E+00), 1.803E+00);
558 addPoint(log(7.350E+00), 1.818E+00);
559 addPoint(log(7.450E+00), 1.793E+00);
560 addPoint(log(7.550E+00), 1.757E+00);
561 addPoint(log(7.650E+00), 1.794E+00);
562 addPoint(log(7.750E+00), 1.866E+00);
563 addPoint(log(7.850E+00), 1.861E+00);
564 addPoint(log(7.950E+00), 1.992E+00);
565 addPoint(log(8.050E+00), 1.806E+00);
566 addPoint(log(8.150E+00), 1.663E+00);
567 addPoint(log(8.250E+00), 1.743E+00);
568 addPoint(log(8.350E+00), 1.768E+00);
569 addPoint(log(8.450E+00), 1.782E+00);
570 addPoint(log(8.550E+00), 1.771E+00);
571 addPoint(log(8.650E+00), 1.769E+00);
572 addPoint(log(8.750E+00), 1.589E+00);
573 addPoint(log(8.850E+00), 1.685E+00);
574 addPoint(log(8.950E+00), 1.871E+00);
575 addPoint(log(9.050E+00), 1.876E+00);
576 addPoint(log(9.150E+00), 1.686E+00);
577 addPoint(log(9.250E+00), 1.619E+00);
578 addPoint(log(9.350E+00), 1.751E+00);
579 addPoint(log(9.450E+00), 1.824E+00);
580 addPoint(log(9.550E+00), 1.767E+00);
581 addPoint(log(9.650E+00), 1.669E+00);
582 addPoint(log(9.750E+00), 1.707E+00);
583 addPoint(log(9.850E+00), 1.754E+00);
584 addPoint(log(9.950E+00), 1.780E+00);
585 addPoint(log(1.005E+01), 1.760E+00);
586 addPoint(log(1.015E+01), 1.782E+00);
587 addPoint(log(1.025E+01), 1.764E+00);
588 addPoint(log(1.035E+01), 1.772E+00);
589 addPoint(log(1.045E+01), 1.709E+00);
590 addPoint(log(1.055E+01), 1.696E+00);
591 addPoint(log(1.065E+01), 1.674E+00);
592 addPoint(log(1.075E+01), 1.726E+00);
593 addPoint(log(1.085E+01), 1.765E+00);
594 addPoint(log(1.095E+01), 1.765E+00);
595 addPoint(log(1.105E+01), 1.736E+00);
596 addPoint(log(1.115E+01), 1.680E+00);
597 addPoint(log(1.125E+01), 1.692E+00);
598 addPoint(log(1.135E+01), 1.735E+00);
599 addPoint(log(1.145E+01), 1.765E+00);
600 addPoint(log(1.155E+01), 1.778E+00);
601 addPoint(log(1.165E+01), 1.767E+00);
602 addPoint(log(1.175E+01), 1.751E+00);
603 addPoint(log(1.185E+01), 1.754E+00);
604 addPoint(log(1.195E+01), 1.767E+00);
605 addPoint(log(1.205E+01), 1.784E+00);
606 addPoint(log(1.215E+01), 1.812E+00);
607 addPoint(log(1.225E+01), 1.816E+00);
608 addPoint(log(1.235E+01), 1.822E+00);
609 addPoint(log(1.245E+01), 1.791E+00);
610 addPoint(log(1.255E+01), 1.757E+00);
611 addPoint(log(1.265E+01), 1.837E+00);
612 addPoint(log(1.275E+01), 1.856E+00);
613 addPoint(log(1.285E+01), 1.850E+00);
614 addPoint(log(1.295E+01), 1.853E+00);
615 addPoint(log(1.305E+01), 1.862E+00);
616 addPoint(log(1.315E+01), 1.828E+00);
617 addPoint(log(1.325E+01), 1.793E+00);
618 addPoint(log(1.335E+01), 1.806E+00);
619 addPoint(log(1.345E+01), 1.839E+00);
620 addPoint(log(1.355E+01), 1.802E+00);
621 addPoint(log(1.365E+01), 1.780E+00);
622 addPoint(log(1.375E+01), 1.804E+00);
623 addPoint(log(1.385E+01), 1.870E+00);
624 addPoint(log(1.395E+01), 1.808E+00);
625 addPoint(log(1.405E+01), 1.787E+00);
626 addPoint(log(1.415E+01), 1.817E+00);
627 addPoint(log(1.425E+01), 1.830E+00);
628 addPoint(log(1.435E+01), 1.843E+00);
629 addPoint(log(1.445E+01), 1.845E+00);
630 addPoint(log(1.455E+01), 1.791E+00);
631 addPoint(log(1.465E+01), 1.765E+00);
632 addPoint(log(1.475E+01), 1.774E+00);
633 addPoint(log(1.485E+01), 1.808E+00);
634 addPoint(log(1.495E+01), 1.813E+00);
635 addPoint(log(1.505E+01), 1.773E+00);
636 addPoint(log(1.515E+01), 1.801E+00);
637 addPoint(log(1.525E+01), 1.842E+00);
638 addPoint(log(1.535E+01), 1.841E+00);
639 addPoint(log(1.545E+01), 1.828E+00);
640 addPoint(log(1.555E+01), 1.808E+00);
641 addPoint(log(1.565E+01), 1.806E+00);
642 addPoint(log(1.575E+01), 1.817E+00);
643 addPoint(log(1.585E+01), 1.858E+00);
644 addPoint(log(1.595E+01), 1.891E+00);
645 addPoint(log(1.605E+01), 1.899E+00);
646 addPoint(log(1.615E+01), 1.899E+00);
647 addPoint(log(1.625E+01), 1.893E+00);
648 addPoint(log(1.635E+01), 1.880E+00);
649 addPoint(log(1.645E+01), 1.869E+00);
650 addPoint(log(1.655E+01), 1.918E+00);
651 addPoint(log(1.665E+01), 1.970E+00);
652 addPoint(log(1.675E+01), 1.952E+00);
653 addPoint(log(1.685E+01), 1.929E+00);
654 addPoint(log(1.695E+01), 1.935E+00);
655 addPoint(log(1.705E+01), 1.940E+00);
656 addPoint(log(1.715E+01), 1.937E+00);
657 addPoint(log(1.725E+01), 1.934E+00);
658 addPoint(log(1.735E+01), 1.935E+00);
659 addPoint(log(1.745E+01), 1.932E+00);
660 addPoint(log(1.755E+01), 1.917E+00);
661 addPoint(log(1.765E+01), 1.906E+00);
662 addPoint(log(1.775E+01), 1.914E+00);
663 addPoint(log(1.785E+01), 1.926E+00);
664 addPoint(log(1.795E+01), 1.951E+00);
665 addPoint(log(1.805E+01), 1.974E+00);
666 addPoint(log(1.815E+01), 1.965E+00);
667 addPoint(log(1.825E+01), 1.950E+00);
668 addPoint(log(1.835E+01), 1.938E+00);
669 addPoint(log(1.845E+01), 1.926E+00);
670 addPoint(log(1.855E+01), 1.945E+00);
671 addPoint(log(1.865E+01), 1.986E+00);
672 addPoint(log(1.875E+01), 2.008E+00);
673 addPoint(log(1.885E+01), 1.999E+00);
674 addPoint(log(1.895E+01), 1.989E+00);
675 addPoint(log(1.905E+01), 1.974E+00);
676 addPoint(log(1.915E+01), 1.960E+00);
677 addPoint(log(1.925E+01), 1.955E+00);
678 addPoint(log(1.935E+01), 1.953E+00);
679 addPoint(log(1.945E+01), 1.953E+00);
680 addPoint(log(1.955E+01), 1.955E+00);
681 addPoint(log(1.965E+01), 1.958E+00);
682 addPoint(log(1.975E+01), 1.963E+00);
683 addPoint(log(1.985E+01), 1.969E+00);
684 addPoint(log(1.995E+01), 1.975E+00);
685 addPoint(log(2.000E+01), 2.071E+00);
686 addPoint(log(2.500E+01), 2.049E+00);
687 addPoint(log(3.000E+01), 2.041E+00);
688 addPoint(log(4.000E+01), 2.012E+00);
689 addPoint(log(5.000E+01), 1.805E+00);
690 addPoint(log(6.000E+01), 1.644E+00);
691 addPoint(log(7.000E+01), 1.499E+00);
692 addPoint(log(8.000E+01), 1.378E+00);
693 addPoint(log(9.000E+01), 1.264E+00);
694 addPoint(log(1.000E+02), 1.168E+00);
695 addPoint(log(1.300E+02), 9.740E-01);
696 addPoint(log(1.600E+02), 8.650E-01);
697 addPoint(log(2.000E+02), 7.910E-01);
698 addPoint(log(2.500E+02), 7.330E-01);
699 addPoint(log(3.000E+02), 6.960E-01);
700 addPoint(log(3.500E+02), 6.930E-01);
701 addPoint(log(4.000E+02), 6.850E-01);
702 addPoint(log(4.500E+02), 6.850E-01);
703 addPoint(log(5.000E+02), 6.810E-01);
704 addPoint(log(6.000E+02), 6.790E-01);
705 addPoint(log(7.000E+02), 6.770E-01);
706 addPoint(log(8.000E+02), 6.720E-01);
707 addPoint(log(8.050E+02), 6.004E-01);
708 addPoint(log(8.150E+02), 5.980E-01);
709 addPoint(log(8.250E+02), 5.959E-01);
710 addPoint(log(8.350E+02), 5.942E-01);
711 addPoint(log(8.450E+02), 5.932E-01);
712 addPoint(log(8.550E+02), 5.922E-01);
713 addPoint(log(8.650E+02), 5.912E-01);
714 addPoint(log(8.750E+02), 5.902E-01);
715 addPoint(log(8.850E+02), 5.892E-01);
716 addPoint(log(8.950E+02), 5.883E-01);
717 addPoint(log(9.050E+02), 5.873E-01);
718 addPoint(log(9.150E+02), 5.863E-01);
719 addPoint(log(9.250E+02), 5.854E-01);
720 addPoint(log(9.350E+02), 5.845E-01);
721 addPoint(log(9.450E+02), 5.836E-01);
722 addPoint(log(9.550E+02), 5.828E-01);
723 addPoint(log(9.650E+02), 5.819E-01);
724 addPoint(log(9.750E+02), 5.810E-01);
725 addPoint(log(9.850E+02), 5.801E-01);
726 addPoint(log(9.950E+02), 5.793E-01);
727 addPoint(log(1.005E+03), 5.784E-01);
728 addPoint(log(1.015E+03), 5.773E-01);
729 addPoint(log(1.025E+03), 5.763E-01);
730 addPoint(log(1.035E+03), 5.752E-01);
731 addPoint(log(1.045E+03), 5.741E-01);
732 addPoint(log(1.055E+03), 5.734E-01);
733 addPoint(log(1.065E+03), 5.728E-01);
734 addPoint(log(1.075E+03), 5.722E-01);
735 addPoint(log(1.085E+03), 5.716E-01);
736 addPoint(log(1.095E+03), 5.711E-01);
737 addPoint(log(1.105E+03), 5.708E-01);
738 addPoint(log(1.115E+03), 5.704E-01);
739 addPoint(log(1.125E+03), 5.701E-01);
740 addPoint(log(1.135E+03), 5.698E-01);
741 addPoint(log(1.145E+03), 5.694E-01);
742 addPoint(log(1.155E+03), 5.691E-01);
743 addPoint(log(1.165E+03), 5.688E-01);
744 addPoint(log(1.175E+03), 5.684E-01);
745 addPoint(log(1.185E+03), 5.681E-01);
746 addPoint(log(1.195E+03), 5.678E-01);
747 addPoint(log(1.205E+03), 5.674E-01);
748 addPoint(log(1.215E+03), 5.671E-01);
749 addPoint(log(1.225E+03), 5.668E-01);
750 addPoint(log(1.235E+03), 5.664E-01);
751 addPoint(log(1.245E+03), 5.661E-01);
752 addPoint(log(1.255E+03), 5.658E-01);
753 addPoint(log(1.265E+03), 5.655E-01);
754 addPoint(log(1.275E+03), 5.650E-01);
755 addPoint(log(1.285E+03), 5.643E-01);
756 addPoint(log(1.295E+03), 5.636E-01);
757 addPoint(log(1.305E+03), 5.628E-01);
758 addPoint(log(1.315E+03), 5.621E-01);
759 addPoint(log(1.325E+03), 5.614E-01);
760 addPoint(log(1.335E+03), 5.607E-01);
761 addPoint(log(1.345E+03), 5.600E-01);
762 addPoint(log(1.355E+03), 5.593E-01);
763 addPoint(log(1.365E+03), 5.587E-01);
764 addPoint(log(1.375E+03), 5.581E-01);
765 addPoint(log(1.385E+03), 5.575E-01);
766 addPoint(log(1.395E+03), 5.569E-01);
767 addPoint(log(1.405E+03), 5.563E-01);
768 addPoint(log(1.415E+03), 5.557E-01);
769 addPoint(log(1.425E+03), 5.551E-01);
770 addPoint(log(1.435E+03), 5.545E-01);
771 addPoint(log(1.445E+03), 5.539E-01);
772 addPoint(log(1.455E+03), 5.533E-01);
773 addPoint(log(1.465E+03), 5.527E-01);
774 addPoint(log(1.485E+03), 5.515E-01);
775 addPoint(log(1.495E+03), 5.510E-01);
776 addPoint(log(1.505E+03), 5.504E-01);
777 addPoint(log(1.515E+03), 5.498E-01);
778 addPoint(log(1.535E+03), 5.486E-01);
779 addPoint(log(1.545E+03), 5.480E-01);
780 addPoint(log(1.555E+03), 5.474E-01);
781 addPoint(log(1.565E+03), 5.468E-01);
782 addPoint(log(1.585E+03), 5.457E-01);
783 addPoint(log(1.595E+03), 5.451E-01);
784 addPoint(log(1.605E+03), 5.445E-01);
785 addPoint(log(1.615E+03), 5.439E-01);
786 addPoint(log(1.635E+03), 5.427E-01);
787 addPoint(log(1.645E+03), 5.421E-01);
788 addPoint(log(1.655E+03), 5.416E-01);
789 addPoint(log(1.665E+03), 5.410E-01);
790 addPoint(log(1.695E+03), 5.392E-01);
791 addPoint(log(1.705E+03), 5.387E-01);
792 addPoint(log(1.715E+03), 5.381E-01);
793 addPoint(log(1.735E+03), 5.369E-01);
794 addPoint(log(1.755E+03), 5.358E-01);
795 addPoint(log(1.765E+03), 5.354E-01);
796 addPoint(log(1.785E+03), 5.345E-01);
797 addPoint(log(1.795E+03), 5.341E-01);
798 addPoint(log(1.815E+03), 5.332E-01);
799 addPoint(log(1.835E+03), 5.323E-01);
800 addPoint(log(1.845E+03), 5.319E-01);
801 addPoint(log(1.855E+03), 5.315E-01);
802 addPoint(log(1.885E+03), 5.302E-01);
803 addPoint(log(1.895E+03), 5.297E-01);
804 addPoint(log(1.905E+03), 5.293E-01);
805 addPoint(log(1.915E+03), 5.289E-01);
806 addPoint(log(1.945E+03), 5.276E-01);
807 addPoint(log(1.955E+03), 5.272E-01);
808 addPoint(log(1.965E+03), 5.267E-01);
809 addPoint(log(1.985E+03), 5.259E-01);
810 addPoint(log(2.005E+03), 5.250E-01);
811 addPoint(log(2.015E+03), 5.246E-01);
812 addPoint(log(2.035E+03), 5.237E-01);
813 addPoint(log(2.045E+03), 5.233E-01);
814 addPoint(log(2.065E+03), 5.225E-01);
815 addPoint(log(2.085E+03), 5.216E-01);
816 addPoint(log(2.095E+03), 5.212E-01);
817 addPoint(log(2.105E+03), 5.208E-01);
818 addPoint(log(2.135E+03), 5.195E-01);
819 addPoint(log(2.145E+03), 5.191E-01);
820 addPoint(log(2.155E+03), 5.187E-01);
821 addPoint(log(2.165E+03), 5.182E-01);
822 addPoint(log(2.195E+03), 5.170E-01);
823 addPoint(log(2.205E+03), 5.166E-01);
824 addPoint(log(2.215E+03), 5.161E-01);
825 addPoint(log(2.235E+03), 5.153E-01);
826 addPoint(log(2.255E+03), 5.145E-01);
827 addPoint(log(2.265E+03), 5.140E-01);
828 addPoint(log(2.285E+03), 5.137E-01);
829 addPoint(log(2.295E+03), 5.136E-01);
830 addPoint(log(2.315E+03), 5.134E-01);
831 addPoint(log(2.335E+03), 5.132E-01);
832 addPoint(log(2.345E+03), 5.131E-01);
833 addPoint(log(2.355E+03), 5.130E-01);
834 addPoint(log(2.385E+03), 5.126E-01);
835 addPoint(log(2.395E+03), 5.125E-01);
836 addPoint(log(2.405E+03), 5.124E-01);
837 addPoint(log(2.415E+03), 5.123E-01);
838 addPoint(log(2.435E+03), 5.121E-01);
839 addPoint(log(2.445E+03), 5.120E-01);
840 addPoint(log(2.455E+03), 5.119E-01);
841 addPoint(log(2.465E+03), 5.118E-01);
842 addPoint(log(2.485E+03), 5.115E-01);
843 addPoint(log(2.495E+03), 5.114E-01);
844 addPoint(log(2.505E+03), 5.113E-01);
845 addPoint(log(2.515E+03), 5.112E-01);
846 addPoint(log(2.535E+03), 5.110E-01);
847 addPoint(log(2.545E+03), 5.109E-01);
848 addPoint(log(2.555E+03), 5.108E-01);
849 addPoint(log(2.565E+03), 5.107E-01);
850 addPoint(log(2.585E+03), 5.105E-01);
851 addPoint(log(2.595E+03), 5.104E-01);
852 addPoint(log(2.605E+03), 5.102E-01);
853 addPoint(log(2.615E+03), 5.101E-01);
854 addPoint(log(2.635E+03), 5.099E-01);
855 addPoint(log(2.645E+03), 5.098E-01);
856 addPoint(log(2.655E+03), 5.097E-01);
857 addPoint(log(2.665E+03), 5.096E-01);
858 addPoint(log(2.685E+03), 5.094E-01);
859 addPoint(log(2.695E+03), 5.093E-01);
860 addPoint(log(2.705E+03), 5.092E-01);
861 addPoint(log(2.715E+03), 5.091E-01);
862 addPoint(log(2.735E+03), 5.088E-01);
863 addPoint(log(2.745E+03), 5.087E-01);
864 addPoint(log(2.755E+03), 5.086E-01);
865 addPoint(log(2.765E+03), 5.085E-01);
866 addPoint(log(2.785E+03), 5.083E-01);
867 addPoint(log(2.795E+03), 5.082E-01);
868 addPoint(log(2.805E+03), 5.081E-01);
869 addPoint(log(2.815E+03), 5.080E-01);
870 addPoint(log(2.835E+03), 5.078E-01);
871 addPoint(log(2.845E+03), 5.077E-01);
872 addPoint(log(2.855E+03), 5.076E-01);
873 addPoint(log(2.865E+03), 5.074E-01);
874 addPoint(log(2.885E+03), 5.072E-01);
875 addPoint(log(2.895E+03), 5.071E-01);
876 addPoint(log(2.905E+03), 5.070E-01);
877 addPoint(log(2.915E+03), 5.069E-01);
878 addPoint(log(2.935E+03), 5.067E-01);
879 addPoint(log(2.945E+03), 5.066E-01);
880 addPoint(log(2.955E+03), 5.065E-01);
881 addPoint(log(2.965E+03), 5.064E-01);
882 addPoint(log(2.985E+03), 5.062E-01);
883 addPoint(log(2.995E+03), 5.060E-01);
884 addPoint(log(3.005E+03), 5.059E-01);
885 addPoint(log(3.015E+03), 5.058E-01);
886 addPoint(log(3.035E+03), 5.056E-01);
887 addPoint(log(3.045E+03), 5.055E-01);
888 addPoint(log(3.055E+03), 5.054E-01);
889 addPoint(log(3.065E+03), 5.053E-01);
890 addPoint(log(3.085E+03), 5.051E-01);
891 addPoint(log(3.095E+03), 5.050E-01);
892 addPoint(log(3.105E+03), 5.049E-01);
893 addPoint(log(3.115E+03), 5.048E-01);
894 addPoint(log(3.135E+03), 5.046E-01);
895 addPoint(log(3.145E+03), 5.044E-01);
896 addPoint(log(3.155E+03), 5.043E-01);
897 addPoint(log(3.165E+03), 5.042E-01);
898 addPoint(log(3.185E+03), 5.040E-01);
899 addPoint(log(3.195E+03), 5.039E-01);
900 addPoint(log(3.205E+03), 5.038E-01);
901 addPoint(log(3.215E+03), 5.037E-01);
902 addPoint(log(3.235E+03), 5.035E-01);
903 addPoint(log(3.245E+03), 5.034E-01);
904 addPoint(log(3.255E+03), 5.033E-01);
905 addPoint(log(3.265E+03), 5.032E-01);
906 addPoint(log(3.285E+03), 5.030E-01);
907 addPoint(log(3.295E+03), 5.028E-01);
908 addPoint(log(3.305E+03), 5.027E-01);
909 addPoint(log(3.315E+03), 5.026E-01);
910 addPoint(log(3.335E+03), 5.024E-01);
911 addPoint(log(3.345E+03), 5.023E-01);
912 addPoint(log(3.355E+03), 5.022E-01);
913 addPoint(log(3.365E+03), 5.021E-01);
914 addPoint(log(3.385E+03), 5.019E-01);
915 addPoint(log(3.395E+03), 5.018E-01);
916 addPoint(log(3.405E+03), 5.017E-01);
917 addPoint(log(3.415E+03), 5.016E-01);
918 addPoint(log(3.435E+03), 5.014E-01);
919 addPoint(log(3.445E+03), 5.013E-01);
920 addPoint(log(3.455E+03), 5.011E-01);
921 addPoint(log(3.465E+03), 5.010E-01);
922 addPoint(log(3.485E+03), 5.008E-01);
923 addPoint(log(3.495E+03), 5.007E-01);
924 addPoint(log(3.505E+03), 5.006E-01);
925 addPoint(log(3.515E+03), 5.005E-01);
926 addPoint(log(3.535E+03), 5.003E-01);
927 addPoint(log(3.545E+03), 5.002E-01);
928 addPoint(log(3.555E+03), 5.001E-01);
929 addPoint(log(3.565E+03), 5.000E-01);
930 addPoint(log(3.585E+03), 4.998E-01);
931 addPoint(log(3.595E+03), 4.997E-01);
932 addPoint(log(3.605E+03), 4.996E-01);
933 addPoint(log(3.615E+03), 4.995E-01);
934 addPoint(log(3.635E+03), 4.992E-01);
935 addPoint(log(3.645E+03), 4.991E-01);
936 addPoint(log(3.655E+03), 4.990E-01);
937 addPoint(log(3.665E+03), 4.989E-01);
938 addPoint(log(3.685E+03), 4.987E-01);
939 addPoint(log(3.695E+03), 4.986E-01);
940 addPoint(log(3.705E+03), 4.985E-01);
941 addPoint(log(3.715E+03), 4.984E-01);
942 addPoint(log(3.735E+03), 4.982E-01);
943 addPoint(log(3.745E+03), 4.981E-01);
944 addPoint(log(3.755E+03), 4.980E-01);
945 addPoint(log(3.765E+03), 4.979E-01);
946 addPoint(log(3.785E+03), 4.977E-01);
947 addPoint(log(3.795E+03), 4.976E-01);
948 addPoint(log(3.805E+03), 4.975E-01);
949 addPoint(log(3.815E+03), 4.974E-01);
950 addPoint(log(3.835E+03), 4.971E-01);
951 addPoint(log(3.845E+03), 4.970E-01);
952 addPoint(log(3.855E+03), 4.969E-01);
953 addPoint(log(3.865E+03), 4.968E-01);
954 addPoint(log(3.885E+03), 4.966E-01);
955 addPoint(log(3.895E+03), 4.965E-01);
956 addPoint(log(3.905E+03), 4.964E-01);
957 addPoint(log(3.915E+03), 4.963E-01);
958 addPoint(log(3.935E+03), 4.961E-01);
959 addPoint(log(3.945E+03), 4.960E-01);
960 addPoint(log(3.955E+03), 4.959E-01);
961 addPoint(log(3.965E+03), 4.958E-01);
962 addPoint(log(3.985E+03), 4.956E-01);
963 addPoint(log(3.995E+03), 4.955E-01);
964 addPoint(log(4.005E+03), 4.954E-01);
965 addPoint(log(4.015E+03), 4.953E-01);
966 addPoint(log(4.025E+03), 4.952E-01);
967 addPoint(log(4.035E+03), 4.950E-01);
968 addPoint(log(4.045E+03), 4.949E-01);
969 addPoint(log(4.055E+03), 4.948E-01);
970 addPoint(log(4.065E+03), 4.947E-01);
971 addPoint(log(4.075E+03), 4.946E-01);
972 addPoint(log(4.085E+03), 4.945E-01);
973 addPoint(log(4.095E+03), 4.944E-01);
974 addPoint(log(4.105E+03), 4.943E-01);
975 addPoint(log(4.115E+03), 4.942E-01);
976 addPoint(log(4.125E+03), 4.941E-01);
977 addPoint(log(4.135E+03), 4.940E-01);
978 addPoint(log(4.145E+03), 4.939E-01);
979 addPoint(log(4.155E+03), 4.938E-01);
980 addPoint(log(4.165E+03), 4.937E-01);
981 addPoint(log(4.175E+03), 4.936E-01);
982 addPoint(log(4.185E+03), 4.935E-01);
983 addPoint(log(4.195E+03), 4.934E-01);
984 addPoint(log(4.205E+03), 4.933E-01);
985 addPoint(log(4.215E+03), 4.932E-01);
986 addPoint(log(4.225E+03), 4.931E-01);
987 addPoint(log(4.235E+03), 4.930E-01);
988 addPoint(log(4.245E+03), 4.929E-01);
989 addPoint(log(4.255E+03), 4.928E-01);
990 addPoint(log(4.265E+03), 4.927E-01);
991 addPoint(log(4.275E+03), 4.925E-01);
992 addPoint(log(4.285E+03), 4.924E-01);
993 addPoint(log(4.295E+03), 4.923E-01);
994 addPoint(log(4.305E+03), 4.922E-01);
995 addPoint(log(4.315E+03), 4.921E-01);
996 addPoint(log(4.325E+03), 4.920E-01);
997 addPoint(log(4.335E+03), 4.919E-01);
998 addPoint(log(4.345E+03), 4.918E-01);
999 addPoint(log(4.355E+03), 4.917E-01);
1000 addPoint(log(4.365E+03), 4.916E-01);
1001 addPoint(log(4.375E+03), 4.915E-01);
1002 addPoint(log(4.385E+03), 4.914E-01);
1003 addPoint(log(4.395E+03), 4.913E-01);
1004 addPoint(log(4.405E+03), 4.912E-01);
1005 addPoint(log(4.415E+03), 4.911E-01);
1006 addPoint(log(4.425E+03), 4.910E-01);
1007 addPoint(log(4.435E+03), 4.909E-01);
1008 addPoint(log(4.445E+03), 4.908E-01);
1009 addPoint(log(4.455E+03), 4.907E-01);
1010 addPoint(log(4.465E+03), 4.906E-01);
1011 addPoint(log(4.475E+03), 4.905E-01);
1012 addPoint(log(4.485E+03), 4.904E-01);
1013 addPoint(log(4.495E+03), 4.903E-01);
1014 addPoint(log(4.505E+03), 4.902E-01);
1015 addPoint(log(4.515E+03), 4.901E-01);
1016 addPoint(log(4.525E+03), 4.900E-01);
1017 addPoint(log(4.535E+03), 4.899E-01);
1018 addPoint(log(4.545E+03), 4.898E-01);
1019 addPoint(log(4.555E+03), 4.897E-01);
1020 addPoint(log(4.565E+03), 4.895E-01);
1021 addPoint(log(4.575E+03), 4.894E-01);
1022 addPoint(log(4.585E+03), 4.893E-01);
1023 addPoint(log(4.595E+03), 4.892E-01);
1024 addPoint(log(4.605E+03), 4.891E-01);
1025 addPoint(log(4.615E+03), 4.890E-01);
1026 addPoint(log(4.625E+03), 4.889E-01);
1027 addPoint(log(4.635E+03), 4.888E-01);
1028 addPoint(log(4.645E+03), 4.887E-01);
1029 addPoint(log(4.655E+03), 4.886E-01);
1030 addPoint(log(4.665E+03), 4.885E-01);
1031 addPoint(log(4.675E+03), 4.884E-01);
1032 addPoint(log(4.685E+03), 4.883E-01);
1033 addPoint(log(4.695E+03), 4.882E-01);
1034 addPoint(log(4.705E+03), 4.881E-01);
1035 addPoint(log(4.715E+03), 4.880E-01);
1036 addPoint(log(4.725E+03), 4.879E-01);
1037 addPoint(log(4.735E+03), 4.878E-01);
1038 addPoint(log(4.745E+03), 4.877E-01);
1039 addPoint(log(4.755E+03), 4.876E-01);
1040 addPoint(log(4.765E+03), 4.875E-01);
1041 addPoint(log(4.775E+03), 4.874E-01);
1042 addPoint(log(4.785E+03), 4.873E-01);
1043 addPoint(log(4.795E+03), 4.872E-01);
1044 addPoint(log(4.805E+03), 4.871E-01);
1045 addPoint(log(4.815E+03), 4.870E-01);
1046 addPoint(log(4.825E+03), 4.869E-01);
1047 addPoint(log(4.835E+03), 4.868E-01);
1048 addPoint(log(4.845E+03), 4.867E-01);
1049 addPoint(log(4.855E+03), 4.866E-01);
1050 addPoint(log(4.865E+03), 4.865E-01);
1051 addPoint(log(4.875E+03), 4.864E-01);
1052 addPoint(log(4.885E+03), 4.863E-01);
1053 addPoint(log(4.895E+03), 4.862E-01);
1054 addPoint(log(4.905E+03), 4.861E-01);
1055 addPoint(log(4.915E+03), 4.860E-01);
1056 addPoint(log(4.925E+03), 4.859E-01);
1057 addPoint(log(4.935E+03), 4.858E-01);
1058 addPoint(log(4.945E+03), 4.857E-01);
1059 addPoint(log(4.955E+03), 4.856E-01);
1060 addPoint(log(4.965E+03), 4.854E-01);
1061 addPoint(log(4.975E+03), 4.853E-01);
1062 addPoint(log(4.985E+03), 4.852E-01);
1063 addPoint(log(4.995E+03), 4.851E-01);
1064 addPoint(log(5.005E+03), 4.850E-01);
1065 addPoint(log(5.015E+03), 4.849E-01);
1066 addPoint(log(5.025E+03), 4.848E-01);
1067 addPoint(log(5.035E+03), 4.847E-01);
1068 addPoint(log(5.045E+03), 4.846E-01);
1069 addPoint(log(5.055E+03), 4.845E-01);
1070 addPoint(log(5.065E+03), 4.844E-01);
1071 addPoint(log(5.075E+03), 4.843E-01);
1072 addPoint(log(5.085E+03), 4.842E-01);
1073 addPoint(log(5.095E+03), 4.841E-01);
1074 addPoint(log(5.105E+03), 4.840E-01);
1075 addPoint(log(5.115E+03), 4.839E-01);
1076 addPoint(log(5.125E+03), 4.838E-01);
1077 addPoint(log(5.135E+03), 4.837E-01);
1078 addPoint(log(5.145E+03), 4.836E-01);
1079 addPoint(log(5.155E+03), 4.835E-01);
1080 addPoint(log(5.165E+03), 4.834E-01);
1081 addPoint(log(5.175E+03), 4.833E-01);
1082 addPoint(log(5.185E+03), 4.832E-01);
1083 addPoint(log(5.195E+03), 4.831E-01);
1084 addPoint(log(5.205E+03), 4.830E-01);
1085 addPoint(log(5.215E+03), 4.829E-01);
1086 addPoint(log(5.225E+03), 4.828E-01);
1087 addPoint(log(5.235E+03), 4.827E-01);
1088 addPoint(log(5.245E+03), 4.826E-01);
1089 addPoint(log(5.255E+03), 4.825E-01);
1090 addPoint(log(5.265E+03), 4.824E-01);
1091 addPoint(log(5.275E+03), 4.823E-01);
1092 addPoint(log(5.285E+03), 4.822E-01);
1093 addPoint(log(5.295E+03), 4.821E-01);
1094 addPoint(log(5.305E+03), 4.820E-01);
1095 addPoint(log(5.315E+03), 4.819E-01);
1096 addPoint(log(5.325E+03), 4.818E-01);
1097 addPoint(log(5.335E+03), 4.817E-01);
1098 addPoint(log(5.345E+03), 4.816E-01);
1099 addPoint(log(5.355E+03), 4.815E-01);
1100 addPoint(log(5.365E+03), 4.814E-01);
1101 addPoint(log(5.375E+03), 4.813E-01);
1102 addPoint(log(5.385E+03), 4.812E-01);
1103 addPoint(log(5.395E+03), 4.811E-01);
1104 addPoint(log(5.405E+03), 4.810E-01);
1105 addPoint(log(5.415E+03), 4.809E-01);
1106 addPoint(log(5.425E+03), 4.808E-01);
1107 addPoint(log(5.435E+03), 4.807E-01);
1108 addPoint(log(5.445E+03), 4.806E-01);
1109 addPoint(log(5.455E+03), 4.805E-01);
1110 addPoint(log(5.465E+03), 4.804E-01);
1111 addPoint(log(5.475E+03), 4.803E-01);
1112 addPoint(log(5.485E+03), 4.802E-01);
1113 addPoint(log(5.495E+03), 4.801E-01);
1114 addPoint(log(5.505E+03), 4.800E-01);
1115 addPoint(log(5.515E+03), 4.799E-01);
1116 addPoint(log(5.525E+03), 4.798E-01);
1117 addPoint(log(5.535E+03), 4.797E-01);
1118 addPoint(log(5.545E+03), 4.796E-01);
1119 addPoint(log(5.555E+03), 4.795E-01);
1120 addPoint(log(5.565E+03), 4.794E-01);
1121 addPoint(log(5.575E+03), 4.793E-01);
1122 addPoint(log(5.585E+03), 4.792E-01);
1123 addPoint(log(5.595E+03), 4.791E-01);
1124 addPoint(log(5.605E+03), 4.790E-01);
1125 addPoint(log(5.615E+03), 4.789E-01);
1126 addPoint(log(5.625E+03), 4.788E-01);
1127 addPoint(log(5.635E+03), 4.787E-01);
1128 addPoint(log(5.645E+03), 4.786E-01);
1129 addPoint(log(5.655E+03), 4.785E-01);
1130 addPoint(log(5.665E+03), 4.784E-01);
1131 addPoint(log(5.675E+03), 4.783E-01);
1132 addPoint(log(5.685E+03), 4.782E-01);
1133 addPoint(log(5.695E+03), 4.781E-01);
1134 addPoint(log(5.705E+03), 4.780E-01);
1135 addPoint(log(5.715E+03), 4.779E-01);
1136 addPoint(log(5.725E+03), 4.778E-01);
1137 addPoint(log(5.735E+03), 4.777E-01);
1138 addPoint(log(5.745E+03), 4.776E-01);
1139 addPoint(log(5.755E+03), 4.775E-01);
1140 addPoint(log(5.765E+03), 4.774E-01);
1141 addPoint(log(5.775E+03), 4.773E-01);
1142 addPoint(log(5.785E+03), 4.772E-01);
1143 addPoint(log(5.795E+03), 4.771E-01);
1144 addPoint(log(5.805E+03), 4.770E-01);
1145 addPoint(log(5.815E+03), 4.769E-01);
1146 addPoint(log(5.825E+03), 4.768E-01);
1147 addPoint(log(5.835E+03), 4.767E-01);
1148 addPoint(log(5.845E+03), 4.766E-01);
1149 addPoint(log(5.855E+03), 4.765E-01);
1150 addPoint(log(5.865E+03), 4.764E-01);
1151 addPoint(log(5.875E+03), 4.763E-01);
1152 addPoint(log(5.885E+03), 4.762E-01);
1153 addPoint(log(5.895E+03), 4.761E-01);
1154 addPoint(log(5.905E+03), 4.760E-01);
1155 addPoint(log(5.915E+03), 4.759E-01);
1156 addPoint(log(5.925E+03), 4.758E-01);
1157 addPoint(log(5.935E+03), 4.757E-01);
1158 addPoint(log(5.945E+03), 4.756E-01);
1159 addPoint(log(5.955E+03), 4.755E-01);
1160 addPoint(log(5.965E+03), 4.754E-01);
1161 addPoint(log(5.975E+03), 4.753E-01);
1162 addPoint(log(5.985E+03), 4.752E-01);
1163 addPoint(log(5.995E+03), 4.751E-01);
1164 addPoint(log(6.005E+03), 4.750E-01);
1165 addPoint(log(6.015E+03), 4.749E-01);
1166 addPoint(log(6.025E+03), 4.748E-01);
1167 addPoint(log(6.035E+03), 4.747E-01);
1168 addPoint(log(6.045E+03), 4.746E-01);
1169 addPoint(log(6.055E+03), 4.745E-01);
1170 addPoint(log(6.065E+03), 4.744E-01);
1171 addPoint(log(6.075E+03), 4.743E-01);
1172 addPoint(log(6.085E+03), 4.742E-01);
1173 addPoint(log(6.095E+03), 4.741E-01);
1174 addPoint(log(6.105E+03), 4.740E-01);
1175 addPoint(log(6.115E+03), 4.739E-01);
1176 addPoint(log(6.125E+03), 4.738E-01);
1177 addPoint(log(6.135E+03), 4.737E-01);
1178 addPoint(log(6.145E+03), 4.736E-01);
1179 addPoint(log(6.155E+03), 4.735E-01);
1180 addPoint(log(6.165E+03), 4.734E-01);
1181 addPoint(log(6.175E+03), 4.733E-01);
1182 addPoint(log(6.185E+03), 4.732E-01);
1183 addPoint(log(6.195E+03), 4.731E-01);
1184 addPoint(log(6.205E+03), 4.730E-01);
1185 addPoint(log(6.215E+03), 4.729E-01);
1186 addPoint(log(6.225E+03), 4.728E-01);
1187 addPoint(log(6.235E+03), 4.727E-01);
1188 addPoint(log(6.245E+03), 4.726E-01);
1189 addPoint(log(6.255E+03), 4.725E-01);
1190 addPoint(log(6.265E+03), 4.724E-01);
1191 addPoint(log(6.275E+03), 4.723E-01);
1192 addPoint(log(6.285E+03), 4.722E-01);
1193 addPoint(log(6.295E+03), 4.721E-01);
1194 addPoint(log(6.305E+03), 4.720E-01);
1195 addPoint(log(6.315E+03), 4.719E-01);
1196 addPoint(log(6.325E+03), 4.718E-01);
1197 addPoint(log(6.335E+03), 4.717E-01);
1198 addPoint(log(6.345E+03), 4.716E-01);
1199 addPoint(log(6.355E+03), 4.715E-01);
1200 addPoint(log(6.365E+03), 4.714E-01);
1201 addPoint(log(6.375E+03), 4.713E-01);
1202 addPoint(log(6.385E+03), 4.712E-01);
1203 addPoint(log(6.395E+03), 4.711E-01);
1204 addPoint(log(6.405E+03), 4.710E-01);
1205 addPoint(log(6.415E+03), 4.709E-01);
1206 addPoint(log(6.425E+03), 4.708E-01);
1207 addPoint(log(6.435E+03), 4.707E-01);
1208 addPoint(log(6.445E+03), 4.706E-01);
1209 addPoint(log(6.455E+03), 4.705E-01);
1210 addPoint(log(6.465E+03), 4.704E-01);
1211 addPoint(log(6.475E+03), 4.703E-01);
1212 addPoint(log(6.485E+03), 4.702E-01);
1213 addPoint(log(6.495E+03), 4.701E-01);
1214 addPoint(log(6.505E+03), 4.700E-01);
1215 addPoint(log(6.515E+03), 4.699E-01);
1216 addPoint(log(6.525E+03), 4.698E-01);
1217 addPoint(log(6.535E+03), 4.697E-01);
1218 addPoint(log(6.545E+03), 4.696E-01);
1219 addPoint(log(6.555E+03), 4.695E-01);
1220 addPoint(log(6.565E+03), 4.695E-01);
1221 addPoint(log(6.575E+03), 4.694E-01);
1222 addPoint(log(6.585E+03), 4.693E-01);
1223 addPoint(log(6.595E+03), 4.692E-01);
1224 addPoint(log(6.605E+03), 4.691E-01);
1225 addPoint(log(6.615E+03), 4.690E-01);
1226 addPoint(log(6.625E+03), 4.689E-01);
1227 addPoint(log(6.635E+03), 4.688E-01);
1228 addPoint(log(6.645E+03), 4.687E-01);
1229 addPoint(log(6.655E+03), 4.686E-01);
1230 addPoint(log(6.665E+03), 4.685E-01);
1231 addPoint(log(6.675E+03), 4.684E-01);
1232 addPoint(log(6.685E+03), 4.683E-01);
1233 addPoint(log(6.695E+03), 4.682E-01);
1234 addPoint(log(6.705E+03), 4.681E-01);
1235 addPoint(log(6.715E+03), 4.680E-01);
1236 addPoint(log(6.725E+03), 4.679E-01);
1237 addPoint(log(6.735E+03), 4.678E-01);
1238 addPoint(log(6.745E+03), 4.677E-01);
1239 addPoint(log(6.755E+03), 4.676E-01);
1240 addPoint(log(6.765E+03), 4.675E-01);
1241 addPoint(log(6.775E+03), 4.674E-01);
1242 addPoint(log(6.785E+03), 4.673E-01);
1243 addPoint(log(6.795E+03), 4.672E-01);
1244 addPoint(log(6.805E+03), 4.671E-01);
1245 addPoint(log(6.815E+03), 4.670E-01);
1246 addPoint(log(6.825E+03), 4.669E-01);
1247 addPoint(log(6.835E+03), 4.668E-01);
1248 addPoint(log(6.845E+03), 4.667E-01);
1249 addPoint(log(6.855E+03), 4.666E-01);
1250 addPoint(log(6.865E+03), 4.665E-01);
1251 addPoint(log(6.875E+03), 4.664E-01);
1252 addPoint(log(6.885E+03), 4.663E-01);
1253 addPoint(log(6.895E+03), 4.662E-01);
1254 addPoint(log(6.905E+03), 4.661E-01);
1255 addPoint(log(6.915E+03), 4.660E-01);
1256 addPoint(log(6.925E+03), 4.659E-01);
1257 addPoint(log(6.935E+03), 4.658E-01);
1258 addPoint(log(6.945E+03), 4.657E-01);
1259 addPoint(log(6.955E+03), 4.656E-01);
1260 addPoint(log(6.965E+03), 4.656E-01);
1261 addPoint(log(6.975E+03), 4.655E-01);
1262 addPoint(log(6.985E+03), 4.654E-01);
1263 addPoint(log(6.995E+03), 4.653E-01);
1264 addPoint(log(7.005E+03), 4.652E-01);
1265 addPoint(log(7.015E+03), 4.651E-01);
1266 addPoint(log(7.025E+03), 4.650E-01);
1267 addPoint(log(7.035E+03), 4.649E-01);
1268 addPoint(log(7.045E+03), 4.648E-01);
1269 addPoint(log(7.055E+03), 4.647E-01);
1270 addPoint(log(7.065E+03), 4.646E-01);
1271 addPoint(log(7.075E+03), 4.645E-01);
1272 addPoint(log(7.085E+03), 4.644E-01);
1273 addPoint(log(7.095E+03), 4.643E-01);
1274 addPoint(log(7.105E+03), 4.642E-01);
1275 addPoint(log(7.115E+03), 4.641E-01);
1276 addPoint(log(7.125E+03), 4.640E-01);
1277 addPoint(log(7.135E+03), 4.639E-01);
1278 addPoint(log(7.145E+03), 4.638E-01);
1279 addPoint(log(7.155E+03), 4.637E-01);
1280 addPoint(log(7.165E+03), 4.636E-01);
1281 addPoint(log(7.185E+03), 4.634E-01);
1282 addPoint(log(7.195E+03), 4.633E-01);
1283 addPoint(log(7.205E+03), 4.632E-01);
1284 addPoint(log(7.225E+03), 4.630E-01);
1285 addPoint(log(7.235E+03), 4.629E-01);
1286 addPoint(log(7.245E+03), 4.628E-01);
1287 addPoint(log(7.255E+03), 4.627E-01);
1288 addPoint(log(7.275E+03), 4.626E-01);
1289 addPoint(log(7.285E+03), 4.625E-01);
1290 addPoint(log(7.295E+03), 4.624E-01);
1291 addPoint(log(7.315E+03), 4.622E-01);
1292 addPoint(log(7.325E+03), 4.621E-01);
1293 addPoint(log(7.335E+03), 4.620E-01);
1294 addPoint(log(7.345E+03), 4.619E-01);
1295 addPoint(log(7.365E+03), 4.617E-01);
1296 addPoint(log(7.375E+03), 4.616E-01);
1297 addPoint(log(7.385E+03), 4.615E-01);
1298 addPoint(log(7.395E+03), 4.614E-01);
1299 addPoint(log(7.415E+03), 4.612E-01);
1300 addPoint(log(7.425E+03), 4.611E-01);
1301 addPoint(log(7.435E+03), 4.610E-01);
1302 addPoint(log(7.445E+03), 4.609E-01);
1303 addPoint(log(7.465E+03), 4.607E-01);
1304 addPoint(log(7.475E+03), 4.606E-01);
1305 addPoint(log(7.485E+03), 4.605E-01);
1306 addPoint(log(7.495E+03), 4.604E-01);
1307 addPoint(log(7.515E+03), 4.603E-01);
1308 addPoint(log(7.525E+03), 4.602E-01);
1309 addPoint(log(7.535E+03), 4.601E-01);
1310 addPoint(log(7.545E+03), 4.600E-01);
1311 addPoint(log(7.565E+03), 4.598E-01);
1312 addPoint(log(7.575E+03), 4.597E-01);
1313 addPoint(log(7.585E+03), 4.596E-01);
1314 addPoint(log(7.595E+03), 4.595E-01);
1315 addPoint(log(7.615E+03), 4.593E-01);
1316 addPoint(log(7.625E+03), 4.592E-01);
1317 addPoint(log(7.635E+03), 4.591E-01);
1318 addPoint(log(7.645E+03), 4.590E-01);
1319 addPoint(log(7.665E+03), 4.588E-01);
1320 addPoint(log(7.675E+03), 4.587E-01);
1321 addPoint(log(7.685E+03), 4.586E-01);
1322 addPoint(log(7.695E+03), 4.585E-01);
1323 addPoint(log(7.715E+03), 4.583E-01);
1324 addPoint(log(7.725E+03), 4.582E-01);
1325 addPoint(log(7.735E+03), 4.582E-01);
1326 addPoint(log(7.745E+03), 4.581E-01);
1327 addPoint(log(7.765E+03), 4.579E-01);
1328 addPoint(log(7.775E+03), 4.578E-01);
1329 addPoint(log(7.785E+03), 4.577E-01);
1330 addPoint(log(7.795E+03), 4.576E-01);
1331 addPoint(log(7.815E+03), 4.574E-01);
1332 addPoint(log(7.825E+03), 4.573E-01);
1333 addPoint(log(7.835E+03), 4.572E-01);
1334 addPoint(log(7.845E+03), 4.571E-01);
1335 addPoint(log(7.865E+03), 4.569E-01);
1336 addPoint(log(7.875E+03), 4.568E-01);
1337 addPoint(log(7.885E+03), 4.567E-01);
1338 addPoint(log(7.895E+03), 4.566E-01);
1339 addPoint(log(7.915E+03), 4.564E-01);
1340 addPoint(log(7.925E+03), 4.564E-01);
1341 addPoint(log(7.935E+03), 4.563E-01);
1342 addPoint(log(7.945E+03), 4.562E-01);
1343 addPoint(log(7.965E+03), 4.560E-01);
1344 addPoint(log(7.975E+03), 4.559E-01);
1345 addPoint(log(7.985E+03), 4.558E-01);
1346 addPoint(log(7.995E+03), 4.557E-01);
1347 addPoint(log(8.015E+03), 4.555E-01);
1348 addPoint(log(8.025E+03), 4.554E-01);
1349 addPoint(log(8.035E+03), 4.553E-01);
1350 addPoint(log(8.045E+03), 4.552E-01);
1351 addPoint(log(8.065E+03), 4.550E-01);
1352 addPoint(log(8.075E+03), 4.549E-01);
1353 addPoint(log(8.085E+03), 4.548E-01);
1354 addPoint(log(8.095E+03), 4.547E-01);
1355 addPoint(log(8.115E+03), 4.546E-01);
1356 addPoint(log(8.125E+03), 4.545E-01);
1357 addPoint(log(8.135E+03), 4.544E-01);
1358 addPoint(log(8.145E+03), 4.543E-01);
1359 addPoint(log(8.165E+03), 4.541E-01);
1360 addPoint(log(8.175E+03), 4.540E-01);
1361 addPoint(log(8.185E+03), 4.539E-01);
1362 addPoint(log(8.195E+03), 4.538E-01);
1363 addPoint(log(8.215E+03), 4.536E-01);
1364 addPoint(log(8.225E+03), 4.535E-01);
1365 addPoint(log(8.235E+03), 4.534E-01);
1366 addPoint(log(8.245E+03), 4.533E-01);
1367 addPoint(log(8.265E+03), 4.531E-01);
1368 addPoint(log(8.275E+03), 4.531E-01);
1369 addPoint(log(8.285E+03), 4.530E-01);
1370 addPoint(log(8.295E+03), 4.529E-01);
1371 addPoint(log(8.315E+03), 4.527E-01);
1372 addPoint(log(8.325E+03), 4.526E-01);
1373 addPoint(log(8.335E+03), 4.525E-01);
1374 addPoint(log(8.345E+03), 4.524E-01);
1375 addPoint(log(8.365E+03), 4.522E-01);
1376 addPoint(log(8.375E+03), 4.521E-01);
1377 addPoint(log(8.385E+03), 4.520E-01);
1378 addPoint(log(8.395E+03), 4.519E-01);
1379 addPoint(log(8.415E+03), 4.517E-01);
1380 addPoint(log(8.425E+03), 4.516E-01);
1381 addPoint(log(8.435E+03), 4.516E-01);
1382 addPoint(log(8.445E+03), 4.515E-01);
1383 addPoint(log(8.465E+03), 4.513E-01);
1384 addPoint(log(8.475E+03), 4.512E-01);
1385 addPoint(log(8.485E+03), 4.511E-01);
1386 addPoint(log(8.495E+03), 4.510E-01);
1387 addPoint(log(8.515E+03), 4.508E-01);
1388 addPoint(log(8.525E+03), 4.507E-01);
1389 addPoint(log(8.535E+03), 4.506E-01);
1390 addPoint(log(8.545E+03), 4.505E-01);
1391 addPoint(log(8.565E+03), 4.503E-01);
1392 addPoint(log(8.575E+03), 4.502E-01);
1393 addPoint(log(8.585E+03), 4.502E-01);
1394 addPoint(log(8.595E+03), 4.501E-01);
1395 addPoint(log(8.615E+03), 4.499E-01);
1396 addPoint(log(8.625E+03), 4.498E-01);
1397 addPoint(log(8.635E+03), 4.497E-01);
1398 addPoint(log(8.645E+03), 4.496E-01);
1399 addPoint(log(8.665E+03), 4.494E-01);
1400 addPoint(log(8.675E+03), 4.493E-01);
1401 addPoint(log(8.685E+03), 4.492E-01);
1402 addPoint(log(8.695E+03), 4.491E-01);
1403 addPoint(log(8.715E+03), 4.489E-01);
1404 addPoint(log(8.725E+03), 4.489E-01);
1405 addPoint(log(8.735E+03), 4.488E-01);
1406 addPoint(log(8.745E+03), 4.487E-01);
1407 addPoint(log(8.765E+03), 4.485E-01);
1408 addPoint(log(8.775E+03), 4.484E-01);
1409 addPoint(log(8.785E+03), 4.483E-01);
1410 addPoint(log(8.795E+03), 4.482E-01);
1411 addPoint(log(8.815E+03), 4.480E-01);
1412 addPoint(log(8.825E+03), 4.479E-01);
1413 addPoint(log(8.835E+03), 4.478E-01);
1414 addPoint(log(8.845E+03), 4.477E-01);
1415 addPoint(log(8.865E+03), 4.476E-01);
1416 addPoint(log(8.875E+03), 4.475E-01);
1417 addPoint(log(8.885E+03), 4.474E-01);
1418 addPoint(log(8.895E+03), 4.473E-01);
1419 addPoint(log(8.915E+03), 4.471E-01);
1420 addPoint(log(8.925E+03), 4.470E-01);
1421 addPoint(log(8.935E+03), 4.469E-01);
1422 addPoint(log(8.945E+03), 4.468E-01);
1423 addPoint(log(8.965E+03), 4.466E-01);
1424 addPoint(log(8.975E+03), 4.465E-01);
1425 addPoint(log(8.995E+03), 4.464E-01);
1426 doSpline();
1427 };
Belle2::BkgNeutronWeight::addPoint
void addPoint(double, double)
Function to add one (x,y) data point to the spline table.
Definition: BkgNeutronWeight.cc:1429
Belle2::BkgNeutronWeight::doSpline
void doSpline()
Function to calculate the cubic spline derivatives.
Definition: BkgNeutronWeight.cc:1438

◆ BkgSensitiveDetector()

BkgSensitiveDetector ( const char *  subDet,
int  iden = 0 
)
explicit

Constructor.

Parameters
subDetname of the subdetector (i.e. "IR","PXD","SVD",...)
idenidentifier of subdetector component (optional, 0 by default)

Definition at line 31 of file BkgSensitiveDetector.cc.

31 :
32 Simulation::SensitiveDetectorBase("BKG", Const::invalidDetector), m_trackID(0), m_startPos(0., 0., 0.), m_startMom(0., 0., 0.),
33 m_startTime(0), m_startEnergy(0), m_energyDeposit(0), m_trackLength(0.)
34 {
35 // registration of store arrays and relations
36
37 StoreArray<MCParticle> mcParticles;
38 StoreArray<BeamBackHit> beamBackHits;
39
40 beamBackHits.registerInDataStore();
41 mcParticles.registerRelationTo(beamBackHits);
42
43 // additional registration of MCParticle relation
44
45 RelationArray relation(mcParticles, beamBackHits);
46 registerMCParticleRelation(relation);
47
48 m_eclrepscale = 0;
49 std::string subDet = subDett;
50 if (subDet == "IR") m_subDet = 0;
51 else if (subDet == "PXD") m_subDet = 1;
52 else if (subDet == "SVD") m_subDet = 2;
53 else if (subDet == "CDC") m_subDet = 3;
54 else if (subDet == "ARICH") m_subDet = 4;
55 else if (subDet == "TOP") m_subDet = 5;
56 else if (subDet == "ECL") m_subDet = 6;
57 else if (subDet == "EKLM") m_subDet = 7;
58 else if (subDet == "BKLM") m_subDet = 8;
59 else m_subDet = 99;
60
61 if (m_subDet == 6) {
62 const int nc[] = {3, 4, 6, 9, 2, 9, 6, 4}; // the number of crystals in a sector per ring
63 const int indx[] = {96, 288, 864, 1151, 7776, 8064, 8544, 8736}; // regions with the same number of crystals in a ring
64 std::vector<int> indxv(indx, indx + sizeof(indx) / sizeof(int));
65 m_eclrepscale = nc[upper_bound(indxv.begin(), indxv.end(), iden) - indxv.begin()];
66 }
67
68 m_identifier = iden;
69
70 }
Belle2::BkgSensitiveDetector::m_energyDeposit
double m_energyDeposit
Energy deposited in volume.
Definition: BkgSensitiveDetector.h:69
Belle2::BkgSensitiveDetector::m_startPos
ROOT::Math::XYZVector m_startPos
Particle position at the entrance in volume.
Definition: BkgSensitiveDetector.h:57
Belle2::BkgSensitiveDetector::m_startEnergy
double m_startEnergy
Particle energy at the entrance in volume.
Definition: BkgSensitiveDetector.h:66
Belle2::BkgSensitiveDetector::m_trackLength
double m_trackLength
Length of the track in the volume.
Definition: BkgSensitiveDetector.h:72
Belle2::BkgSensitiveDetector::m_eclrepscale
int m_eclrepscale
Replica (=sector) scale in ECL.
Definition: BkgSensitiveDetector.h:75
Belle2::BkgSensitiveDetector::m_subDet
int m_subDet
Subdetector id number.
Definition: BkgSensitiveDetector.h:48
Belle2::BkgSensitiveDetector::m_startMom
ROOT::Math::XYZVector m_startMom
Particle momentum at the entrance in volume.
Definition: BkgSensitiveDetector.h:60
Belle2::BkgSensitiveDetector::m_identifier
int m_identifier
Identifier of subdetector component.
Definition: BkgSensitiveDetector.h:51
Belle2::Simulation::SensitiveDetectorBase::registerMCParticleRelation
static void registerMCParticleRelation(const std::string &name, RelationArray::EConsolidationAction ignoreAction=RelationArray::c_negativeWeight)
Register an relation involving MCParticles.
Definition: SensitiveDetectorBase.cc:22

◆ doSpline()

void doSpline ( )
private

Function to calculate the cubic spline derivatives.

Definition at line 1438 of file BkgNeutronWeight.cc.

1439 {
1440 // converted from fortran in Bevington's Data Analysis ... ugh!
1441 m_n = m_x.size();
1442 int n = m_n;
1443 int i = 0;
1444
1445 m_d[0] = m_x[1] - m_x[0];
1446 m_c[1] = (m_y[1] - m_y[0]) / m_d[0];
1447
1448 for (i = 1; i < n - 1; i++) {
1449 m_d[i] = m_x[i + 1] - m_x[i];
1450 m_b[i] = 2.0 * (m_d[i - 1] + m_d[i]);
1451 m_c[i + 1] = (m_y[i + 1] - m_y[i]) / m_d[i];
1452 m_c[i] = m_c[i + 1] - m_c[i];
1453 }
1454
1455 m_b[0] = -m_d[0];
1456 m_b[n - 1] = -m_d[n - 2];
1457
1458 m_c[0] = m_c[2] / (m_x[3] - m_x[1]) - m_c[1] / (m_x[2] - m_x[0]);
1459 m_c[0] = m_c[0] * m_d[0] * m_d[0] / (m_x[3] - m_x[0]);
1460 m_c[n - 1] = m_c[n - 2] / (m_x[n - 1] - m_x[n - 3]) - m_c[n - 3] / (m_x[n - 2] - m_x[n - 4]);
1461 m_c[n - 1] = -m_c[n - 1] * m_d[n - 2] * m_d[n - 2] / (m_x[n - 1] - m_x[n - 4]);
1462
1463 for (i = 1; i < n; i++) {
1464 double t = m_d[i - 1] / m_b[i - 1];
1465 m_b[i] = m_b[i] - t * m_d[i - 1];
1466 m_c[i] = m_c[i] - t * m_c[i - 1];
1467 }
1468
1469 m_c[n - 1] = m_c[n - 1] / m_b[n - 1];
1470
1471 for (i = n - 2; i >= 0; i--) {
1472 m_c[i] = (m_c[i] - m_d[i] * m_c[i + 1]) / m_b[i];
1473 }
1474
1475 m_b[n - 1] = (m_y[n - 1] - m_y[n - 2]) / m_d[n - 2] + m_d[n - 2] * (m_c[n - 2] + 2.0 * m_c[n - 1]);
1476
1477 for (i = 0; i < n - 1; i++) {
1478 m_b[i] = (m_y[i + 1] - m_y[i]) / m_d[i] - m_d[i] * (m_c[i + 1] + 2.0 * m_c[i]);
1479 m_d[i] = (m_c[i + 1] - m_c[i]) / m_d[i];
1480 m_c[i] = 3.0 * m_c[i];
1481 }
1482
1483 m_c[n - 1] = 3.0 * m_c[n - 1];
1484 m_d[n - 1] = m_d[n - 2];
1485
1486 }
Belle2::BkgNeutronWeight::m_n
int m_n
Number of tabulated points.
Definition: BkgNeutronWeight.h:45

◆ getInstance()

BkgNeutronWeight & getInstance ( )
static

Return a reference to the singleton BkgNeutronWeight instance.

Definition at line 20 of file BkgNeutronWeight.cc.

21 {
22 if (m_Instance == 0) m_Instance = new BkgNeutronWeight();
23 return *m_Instance;
24 }
Belle2::BkgNeutronWeight::m_Instance
static BkgNeutronWeight * m_Instance
static pointer to the singleton instance of this class
Definition: BkgNeutronWeight.h:63
Belle2::BkgNeutronWeight::BkgNeutronWeight
BkgNeutronWeight(void)
Constructor is private for this singleton.
Definition: BkgNeutronWeight.cc:43

◆ GetRemainderLifeTime()

G4double GetRemainderLifeTime ( ) const
inline

Get Remainder of life time at rest decay.

Returns remaining lifetime.

Definition at line 195 of file G4LongLivedNeutralDecay.h.

◆ getWeight()

double getWeight ( double  ke)

Get weighting factor to convert a neutron to its 1-MeV equivalent.

Definition at line 26 of file BkgNeutronWeight.cc.

27 {
28 if (m_Instance == 0) m_Instance = new BkgNeutronWeight();
29 if (ke <= 0) return 0;
30
31 double x = log(ke);
32 int i = 0;
33 int j = m_n;
34 while (j > i + 1) {
35 int k = (i + j) / 2;
36 if (x < m_x[k]) j = k;
37 if (x >= m_x[k]) i = k;
38 }
39 double dx = x - m_x[i];
40 return (m_y[i] + dx * (m_b[i] + dx * (m_c[i] + dx * m_d[i])));
41 }

◆ step()

G4bool step ( G4Step *  aStep,
G4TouchableHistory *   
)
overridevirtual

Process each step and calculate variables defined in PXDSimHit.

Parameters
aStepCurrent Geant4 step in the sensitive medium.
Returns
true if a hit was stored, o.w. false.

Implements SensitiveDetectorBase.

Definition at line 73 of file BkgSensitiveDetector.cc.

74 {
75
76 const G4StepPoint& preStep = *aStep->GetPreStepPoint();
77 G4Track& track = *aStep->GetTrack();
78
79 if (m_trackID != track.GetTrackID()) {
80 //TrackID changed, store track informations
81 m_trackID = track.GetTrackID();
82 //Get world position
83 const G4ThreeVector& worldPosition = preStep.GetPosition();
84 m_startPos.SetXYZ(worldPosition.x() * Unit::mm / Unit::cm, worldPosition.y() * Unit::mm / Unit::cm,
85 worldPosition.z() * Unit::mm / Unit::cm);
86 //Get momentum
87 const G4ThreeVector& momentum = preStep.GetMomentum() ;
88 m_startMom.SetXYZ(momentum.x() * Unit::MeV, momentum.y() * Unit::MeV,
89 momentum.z() * Unit::MeV);
90 //Get time
91 m_startTime = preStep.GetGlobalTime();
92 //Get energy
93 m_startEnergy = preStep.GetKineticEnergy() * Unit::MeV;
94 //Reset energy deposit;
95 m_energyDeposit = 0;
96 //Reset track lenght;
97 m_trackLength = 0;
98 }
99 //Update energy deposit
100 m_energyDeposit += aStep->GetTotalEnergyDeposit() * Unit::MeV;
101
102 m_trackLength += aStep->GetStepLength() * Unit::mm;
103
104 //Save Hit if track leaves volume or is killed
105 if (track.GetNextVolume() != track.GetVolume() || track.GetTrackStatus() >= fStopAndKill) {
106 int pdgCode = track.GetDefinition()->GetPDGEncoding();
107 double endEnergy = track.GetKineticEnergy() * Unit::MeV;
108 double neutWeight = 0;
109 if (pdgCode == Const::neutron.getPDGCode()) {
110 BkgNeutronWeight& wt = BkgNeutronWeight::getInstance();
111 neutWeight = wt.getWeight(m_startEnergy / Unit::MeV);
112 }
113 StoreArray<BeamBackHit> beamBackHits;
114 int nentr = beamBackHits.getEntries();
115 int id = m_identifier;
116 if (m_subDet == 6) { // ECL
117 int sector = preStep.GetTouchable()->GetCopyNumber(1);
118 if ((sector == 0) && (m_eclrepscale == 2) && ((id & 1) != 0)) id += 144; // barrel odd-numbered diodes in sector 0 -> sector 144
119 id += sector * m_eclrepscale;
120 }
121 beamBackHits.appendNew(m_subDet, id, pdgCode, m_trackID, m_startPos, m_startMom, m_startTime, m_startEnergy, endEnergy,
122 m_energyDeposit, m_trackLength, neutWeight);
123
124
125 // create relation to MCParticle
126 StoreArray<MCParticle> mcParticles;
127 RelationArray relBeamBackHitToMCParticle(mcParticles, beamBackHits);
128 relBeamBackHitToMCParticle.add(m_trackID, nentr);
129
130 //Reset TrackID
131 m_trackID = 0;
132 }
133
134 return true;
135 }
Belle2::Const::neutron
static const ParticleType neutron
neutron particle
Definition: Const.h:675
Belle2::Unit::mm
static const double mm
[millimeters]
Definition: Unit.h:70
Belle2::Unit::MeV
static const double MeV
[megaelectronvolt]
Definition: Unit.h:114
Belle2::Unit::cm
static const double cm
Standard units with the value = 1.
Definition: Unit.h:47
Belle2::BkgNeutronWeight::getInstance
static BkgNeutronWeight & getInstance()
Return a reference to the singleton BkgNeutronWeight instance.
Definition: BkgNeutronWeight.cc:20

Variable Documentation

◆ m_Instance

BkgNeutronWeight * m_Instance = 0
staticprivate

static pointer to the singleton instance of this class

Definition at line 63 of file BkgNeutronWeight.h.