| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 // GEANT4 Class header file
27 //
28 //
29 // File name: G4AtimaFluctuations
30 //
31 // Author: Jose Luis Rodriguez Sanchez
32 //
33 // Creation date: 16.01.2018
34 //
35 // Modifications:
36 // 09/10/2018: Deleted Titeica's correction at
37 //
38 //
39 // Class Description:
40 //
41 // This model calculates the energy-loss fluct
42 // developed at GSI, Darmstadt, Germany.
43 // For details: http://web-docs.gsi.de/~weick/
44 //
45 // Helmut Weick, GSI (responsible for fortran
46 // Andrej Prochazka, GSI (responsible for C ve
47 // Christoph Scheidenberger, GSI (project coor
48 //
49 // -------------------------------------------
50 //
51 //....oooOO0OOooo........oooOO0OOooo........oo
52 //....oooOO0OOooo........oooOO0OOooo........oo
53
54 #include "G4AtimaFluctuations.hh"
55 #include "G4PhysicalConstants.hh"
56 #include "G4SystemOfUnits.hh"
57 #include "Randomize.hh"
58 #include "G4Poisson.hh"
59 #include "G4MaterialCutsCouple.hh"
60 #include "G4Material.hh"
61 #include "G4NistManager.hh"
62 #include "G4DynamicParticle.hh"
63 #include "G4Log.hh"
64 #include "G4Exp.hh"
65 #include "G4Pow.hh"
66
67 //....oooOO0OOooo........oooOO0OOooo........oo
68
69 using namespace std;
70
71 G4double G4AtimaFluctuations::stepE = 0.0;
72 G4double G4AtimaFluctuations::tableE[] = {0.0}
73
74 G4AtimaFluctuations::G4AtimaFluctuations(const
75 : G4VEmFluctuationModel(nam),
76 particle(nullptr),
77 particleMass(CLHEP::proton_mass_c2),
78 charge(1.0),
79 chargeSquare(1.0),
80 effChargeSquare(1.0),
81 minLoss(0.001*CLHEP::eV)
82 {
83 g4calc = G4Pow::GetInstance();
84 kineticEnergy = 0.0;
85 beta2 = 0.0;
86 MLN10 = 2.30258509299;
87 atomic_mass_unit = 931.4940954; // MeV/c^2
88 dedx_constant = 0.3070749187; //4*pi*Na*me*
89 electron_mass = 0.510998928; // MeV/c^2
90 fine_structure = 1.0/137.035999139;
91 domega2dx_constant = dedx_constant*electron_
92 if(tableE[0] == 0.0) {
93 G4double logmin = 0.;
94 G4double logmax = 5.;
95 stepE = (logmax-logmin)/(G4double)(199);
96 for(G4int i=0; i<200; ++i){
97 tableE[i] = G4Exp(MLN10*(logmin + ((G4do
98 }
99 }
100 }
101
102 //....oooOO0OOooo........oooOO0OOooo........oo
103
104 G4AtimaFluctuations::~G4AtimaFluctuations() =
105
106 //....oooOO0OOooo........oooOO0OOooo........oo
107
108 void G4AtimaFluctuations::InitialiseMe(const G
109 {
110 particle = part;
111 particleMass = part->GetPDGMass();
112 charge = part->GetPDGCharge()/eplus;
113 chargeSquare = charge*charge;
114 effChargeSquare= chargeSquare;
115 }
116
117 //....oooOO0OOooo........oooOO0OOooo........oo
118
119 G4double
120 G4AtimaFluctuations::SampleFluctuations(const
121 const
122 const
123 const
124 const
125 const
126 {
127 // G4cout << "### meanLoss= " << meanLoss <
128 if(meanLoss <= minLoss) return meanLoss;
129
130 // G4cout << "G4AtimaFluctuations::SampleFluc
131 // << dp->GetKineticEnergy()
132 // << " Elim(MeV)= " << parameter*charge*
133 // << " " << parameter << " " << charge <<
134
135 const G4Material* material = couple->GetMate
136 G4double siga = Dispersion(material,dp,tcut,
137
138 CLHEP::HepRandomEngine* rndmEngine = G4Rando
139 //G4cout << "meanLoss= " << meanLoss << " lo
140 return G4RandGauss::shoot(rndmEngine,meanLos
141 }
142
143 //....oooOO0OOooo........oooOO0OOooo........oo
144
145 G4double G4AtimaFluctuations::Dispersion(const
146 const
147 const G4double, const G4doubl
148 const
149 {
150 kineticEnergy = dp->GetKineticEnergy();
151 const G4ParticleDefinition* p = dp->GetDefin
152 if(particle != p) { InitialiseMe(p); }
153
154 G4double ap = p->GetPDGMass()/atomic_mass_un
155 G4double zp = p->GetPDGCharge();
156 G4double ep = kineticEnergy/ap;// ep in MeV/
157 G4double zt = mat->GetIonisation()->GetZeffe
158 G4double at = G4NistManager::Instance()->Get
159
160 G4int z = G4lrint(zp);
161 if(z > 109) { z = 109; }
162
163 G4double gamma=1.0 + ep/atomic_mass_unit;
164 beta2=1.0-1.0/(gamma*gamma);
165 G4double beta = std::sqrt(beta2);
166 //z_eff_Pierce_Blann(Z, beta);
167 G4double zp_eff = zp*(1.0-G4Exp(-0.95*137.03
168 //
169 G4double f = domega2dx_constant*zp_eff*zp_ef
170 //
171 G4double cor =
172 24.89 * g4calc->powA(zt,1.2324)/(electron_
173 G4Log( 2.0*electron_mass*1e6*beta2/(33.05*
174 cor = std::max(cor, 0.0 );
175
176 //Lindhard corrections
177 if(ep<tableE[0])ep = tableE[0];
178 G4double da = (ap - element_atomic_weights[z
179
180 G4double v3 = EnergyTable_interpolate(tableE
181 G4double v4 = EnergyTable_interpolate(tableE
182 G4double dif = v4 - v3;
183 G4double X = v3+(dif*da/0.05);
184 X *= gamma*gamma;
185 //
186 G4double sse = 0.;
187 if(ep<30.0){
188 //Energy straggling Firsov
189 G4double factor = 4.8184e-3*g4calc->powA(z
190 sse = std::min(f*(X+cor), factor*beta2/fin
191 }else{
192 // sse = f*(X+cor);
193 sse = f*X;//09/10/2018
194 }
195 //
196 return sse*length/(cm)*mat->GetDensity()/(g/
197 }
198
199 //....oooOO0OOooo........oooOO0OOooo........oo
200
201 void G4AtimaFluctuations::SetParticleAndCharge
202
203 {
204 if(part != particle) {
205 particle = part;
206 particleMass = part->GetPDGMass();
207 charge = part->GetPDGCharge()/eplu
208 chargeSquare = charge*charge;
209 }
210 effChargeSquare = q2;
211 }
212
213 //....oooOO0OOooo........oooOO0OOooo........oo
214
215 G4double G4AtimaFluctuations::EnergyTable_inte
216 G4double r;
217 G4int num=200;
218 G4double lxval = G4Log(xval)/MLN10;
219 if(xval<table[0] || xval>table[num-1])return
220 if(xval==table[num-1])return y[num-1];
221 G4int i = (G4int)(lxval/stepE);
222 i = std::min(std::max(i, 0), num-2);
223 G4double linstep = table[i+1] - table[i];
224 G4double x = 1.0 - ((xval - table[i])/linste
225 r = (x*y[i]) + ((1-x)*y[i+1]);
226 return r;
227 }
228
229 //....oooOO0OOooo........oooOO0OOooo........oo
230
231 const G4double G4AtimaFluctuations::element_at
232 0.0,
233 1.00794, //H
234 4.0026, //He
235 6.941, //Li
236 9.01218, //Be
237 10.811, //B
238 12.0107, //C
239 14.0067, //N
240 15.9994, //O
241 18.9984, //F
242 20.1797, //Ne
243 22.9898, //Na
244 24.305, //Mg
245 26.9815, //Al
246 28.0855, //Si
247 30.9738, //P
248 32.065, //S
249 35.453, //Cl
250 39.948, //Ar
251 39.0983, //K
252 40.078, //Ca
253 44.9559, //Sc
254 47.867, //Ti
255 50.9415, //V
256 51.9961, //Cr
257 54.938, //Mn
258 55.845, //Fe
259 58.9332, //Co
260 58.6934, //Ni
261 63.546, //Cu
262 65.409, //Zn
263 69.723, //Ga
264 72.64, //Ge
265 74.9216, //As
266 78.96, //Se
267 79.904, //Br
268 83.798, //Kr
269 85.4678, //Rb
270 87.62, //Sr
271 88.9059, //Y
272 91.224, //Zr
273 92.9064, //Nb
274 95.94, //Mo
275 97.9072, //Tc
276 101.07, //Ru
277 102.906, //Rh
278 106.42, //Pd
279 107.868, //Ag
280 112.411, //Cd
281 114.818, //In
282 118.71, //Sn
283 121.76, //Sb
284 127.6, //Te
285 126.904, //I
286 131.293, //Xe
287 132.905, //Cs
288 137.327, //Ba
289 138.905, //La
290 140.116, //Ce
291 140.908, //Pr
292 144.24, //Nd
293 144.913, //Pm
294 150.36, //Sm
295 151.964, //Eu
296 157.25, //Gd
297 158.925, //Tb
298 162.5, //Dy
299 164.93, //Ho
300 167.259, //Er
301 168.934, //Tm
302 173.04, //Yb
303 174.967, //Lu
304 178.49, //Hf
305 180.948, //Ta
306 183.84, //W
307 186.207, //Re
308 190.23, //Os
309 192.217, //Ir
310 195.078, //Pt
311 196.967, //Au
312 200.59, //Hg
313 204.383, //Tl
314 207.2, //Pb
315 208.98, //Bi
316 208.982, //Po
317 209.987, //At
318 222.018, //Rn
319 223.02, //Fr
320 226.025, //Ra
321 227.028, //Ac
322 232.038, //Th
323 231.036, //Pa
324 238.029, //U
325 237.048, //Np
326 244.064, //Pu
327 243.061, //Am
328 247.07, //Cm
329 247.07, //Bk
330 251.08, //Cf
331 252.083, //Es
332 257.095, //Fm
333 258.098, //Md
334 259.101, //No
335 262.11, //Lr
336 261.109, //Rf
337 262.114, //Db
338 266.122, //Sg
339 264.125, //Bh
340 269.134, //Hs
341 268.139 //Mt
342 };
343
344 //....oooOO0OOooo........oooOO0OOooo........oo
345
346 //arrays dimensions are [z][energy], z=1 start
347
348 //LS X coefficient (dE straggling) for A=atomi
349 const G4double G4AtimaFluctuations::ls_X_coeff
350 {
351 {0.999182,0.999124,0.999062,0.998995,0.99892
352 {0.999599,0.999546,0.999489,0.999428,0.99936
353 {1.00017,1.00012,1.00007,1.00002,0.999959,0.
354 {1.00086,1.00082,1.00078,1.00073,1.00069,1.0
355 {1.00162,1.0016,1.00157,1.00154,1.00151,1.00
356 {1.00243,1.00243,1.00242,1.00241,1.00239,1.0
357 {1.00324,1.00327,1.00328,1.0033,1.0033,1.003
358 {1.00405,1.0041,1.00415,1.00419,1.00422,1.00
359 {1.00483,1.00491,1.00499,1.00506,1.00513,1.0
360 {1.00559,1.0057,1.0058,1.00591,1.00601,1.006
361 {1.0063,1.00644,1.00658,1.00672,1.00686,1.00
362 {1.00697,1.00715,1.00733,1.0075,1.00767,1.00
363 {1.00761,1.00782,1.00803,1.00824,1.00845,1.0
364 {1.0082,1.00845,1.00869,1.00894,1.00918,1.00
365 {1.00876,1.00904,1.00932,1.0096,1.00988,1.01
366 {1.00928,1.00959,1.0099,1.01021,1.01053,1.01
367 {1.00976,1.0101,1.01045,1.0108,1.01115,1.011
368 {1.01022,1.01059,1.01096,1.01134,1.01173,1.0
369 {1.01064,1.01104,1.01144,1.01185,1.01227,1.0
370 {1.01103,1.01146,1.01189,1.01233,1.01279,1.0
371 {1.01139,1.01185,1.01231,1.01278,1.01327,1.0
372 {1.01173,1.01221,1.0127,1.0132,1.01372,1.014
373 {1.01205,1.01255,1.01307,1.0136,1.01414,1.01
374 {1.01235,1.01287,1.01341,1.01397,1.01454,1.0
375 {1.01262,1.01317,1.01373,1.01431,1.01491,1.0
376 {1.01288,1.01345,1.01404,1.01464,1.01526,1.0
377 {1.01312,1.01371,1.01432,1.01494,1.01559,1.0
378 {1.01335,1.01395,1.01458,1.01523,1.0159,1.01
379 {1.01356,1.01418,1.01483,1.01549,1.01618,1.0
380 {1.01375,1.01439,1.01506,1.01574,1.01646,1.0
381 {1.01393,1.01459,1.01527,1.01598,1.01671,1.0
382 {1.0141,1.01478,1.01547,1.0162,1.01695,1.017
383 {1.01426,1.01495,1.01566,1.0164,1.01717,1.01
384 {1.01441,1.01511,1.01584,1.0166,1.01738,1.01
385 {1.01455,1.01526,1.016,1.01678,1.01758,1.018
386 {1.01467,1.0154,1.01616,1.01695,1.01776,1.01
387 {1.01479,1.01553,1.0163,1.0171,1.01794,1.018
388 {1.01491,1.01566,1.01644,1.01725,1.0181,1.01
389 {1.01501,1.01577,1.01656,1.01739,1.01825,1.0
390 {1.0151,1.01587,1.01668,1.01752,1.01839,1.01
391 {1.01519,1.01597,1.01679,1.01764,1.01852,1.0
392 {1.01527,1.01606,1.01689,1.01775,1.01865,1.0
393 {1.01535,1.01615,1.01698,1.01785,1.01876,1.0
394 {1.01542,1.01622,1.01706,1.01794,1.01886,1.0
395 {1.01548,1.01629,1.01714,1.01803,1.01896,1.0
396 {1.01554,1.01636,1.01721,1.01811,1.01905,1.0
397 {1.01559,1.01641,1.01728,1.01819,1.01913,1.0
398 {1.01563,1.01647,1.01734,1.01825,1.01921,1.0
399 {1.01567,1.01651,1.01739,1.01831,1.01928,1.0
400 {1.01571,1.01655,1.01744,1.01837,1.01934,1.0
401 {1.01574,1.01659,1.01748,1.01842,1.01939,1.0
402 {1.01577,1.01662,1.01752,1.01846,1.01944,1.0
403 {1.01579,1.01665,1.01755,1.01849,1.01949,1.0
404 {1.0158,1.01667,1.01758,1.01853,1.01952,1.02
405 {1.01582,1.01669,1.0176,1.01855,1.01955,1.02
406 {1.01583,1.0167,1.01761,1.01857,1.01958,1.02
407 {1.01583,1.01671,1.01763,1.01859,1.0196,1.02
408 {1.01583,1.01671,1.01763,1.0186,1.01962,1.02
409 {1.01583,1.01671,1.01763,1.01861,1.01962,1.0
410 {1.01582,1.01671,1.01763,1.01861,1.01963,1.0
411 {1.01581,1.0167,1.01763,1.0186,1.01963,1.020
412 {1.0158,1.01668,1.01762,1.0186,1.01962,1.020
413 {1.01578,1.01667,1.0176,1.01858,1.01961,1.02
414 {1.01575,1.01665,1.01758,1.01856,1.0196,1.02
415 {1.01573,1.01662,1.01756,1.01854,1.01958,1.0
416 {1.0157,1.01659,1.01753,1.01852,1.01955,1.02
417 {1.01567,1.01656,1.0175,1.01849,1.01952,1.02
418 {1.01563,1.01652,1.01746,1.01845,1.01949,1.0
419 {1.01559,1.01648,1.01742,1.01841,1.01945,1.0
420 {1.01555,1.01644,1.01738,1.01837,1.01941,1.0
421 {1.0155,1.01639,1.01733,1.01832,1.01936,1.02
422 {1.01545,1.01634,1.01728,1.01827,1.0193,1.02
423 {1.01539,1.01628,1.01722,1.01821,1.01925,1.0
424 {1.01533,1.01622,1.01716,1.01815,1.01918,1.0
425 {1.01527,1.01616,1.0171,1.01808,1.01912,1.02
426 {1.01521,1.01609,1.01703,1.01801,1.01904,1.0
427 {1.01514,1.01602,1.01695,1.01794,1.01897,1.0
428 {1.01506,1.01595,1.01688,1.01786,1.01889,1.0
429 {1.01499,1.01587,1.0168,1.01777,1.0188,1.019
430 {1.01491,1.01579,1.01671,1.01768,1.01871,1.0
431 {1.01482,1.0157,1.01662,1.01759,1.01861,1.01
432 {1.01474,1.01561,1.01653,1.0175,1.01851,1.01
433 {1.01465,1.01551,1.01643,1.01739,1.01841,1.0
434 {1.01455,1.01542,1.01633,1.01729,1.0183,1.01
435 {1.01445,1.01531,1.01622,1.01718,1.01818,1.0
436 {1.01435,1.01521,1.01611,1.01706,1.01806,1.0
437 {1.01424,1.0151,1.01599,1.01694,1.01794,1.01
438 {1.01413,1.01498,1.01587,1.01682,1.01781,1.0
439 {1.01402,1.01486,1.01575,1.01669,1.01767,1.0
440 {1.0139,1.01474,1.01562,1.01655,1.01753,1.01
441 {1.01378,1.01461,1.01549,1.01641,1.01739,1.0
442 {1.01365,1.01448,1.01535,1.01627,1.01723,1.0
443 {1.01352,1.01434,1.01521,1.01612,1.01708,1.0
444 {1.01339,1.0142,1.01506,1.01596,1.01692,1.01
445 {1.01325,1.01406,1.01491,1.0158,1.01675,1.01
446 {1.0131,1.01391,1.01475,1.01564,1.01658,1.01
447 {1.01296,1.01375,1.01459,1.01547,1.0164,1.01
448 {1.01281,1.01359,1.01442,1.0153,1.01622,1.01
449 {1.01265,1.01343,1.01425,1.01512,1.01603,1.0
450 {1.01249,1.01326,1.01408,1.01493,1.01583,1.0
451 {1.01233,1.01309,1.01389,1.01474,1.01563,1.0
452 {1.01216,1.01291,1.01371,1.01455,1.01543,1.0
453 {1.01199,1.01273,1.01352,1.01435,1.01522,1.0
454 {1.01181,1.01255,1.01332,1.01414,1.015,1.015
455 {1.01163,1.01236,1.01312,1.01393,1.01478,1.0
456 {1.01145,1.01216,1.01292,1.01372,1.01455,1.0
457 {1.01126,1.01197,1.01271,1.0135,1.01432,1.01
458 {1.01107,1.01176,1.0125,1.01327,1.01409,1.01
459 };
460
461 //....oooOO0OOooo........oooOO0OOooo........oo
462
463 //LS X coefficient for A=atomic weight * 1.05
464 const G4double G4AtimaFluctuations::ls_X_coeff
465 {
466 {0.999182,0.999124,0.999062,0.998995,0.99892
467 {0.999599,0.999546,0.999489,0.999428,0.99936
468 {1.00017,1.00012,1.00007,1.00002,0.999959,0.
469 {1.00086,1.00082,1.00078,1.00073,1.00069,1.0
470 {1.00162,1.0016,1.00157,1.00154,1.00151,1.00
471 {1.00243,1.00243,1.00242,1.00241,1.00239,1.0
472 {1.00324,1.00327,1.00328,1.0033,1.0033,1.003
473 {1.00405,1.0041,1.00415,1.00419,1.00422,1.00
474 {1.00483,1.00491,1.00499,1.00506,1.00513,1.0
475 {1.00559,1.0057,1.0058,1.00591,1.00601,1.006
476 {1.0063,1.00644,1.00658,1.00672,1.00686,1.00
477 {1.00697,1.00715,1.00733,1.0075,1.00767,1.00
478 {1.00761,1.00782,1.00803,1.00824,1.00845,1.0
479 {1.0082,1.00845,1.00869,1.00894,1.00918,1.00
480 {1.00876,1.00904,1.00932,1.0096,1.00988,1.01
481 {1.00928,1.00959,1.0099,1.01021,1.01053,1.01
482 {1.00976,1.0101,1.01045,1.0108,1.01115,1.011
483 {1.01022,1.01059,1.01096,1.01134,1.01173,1.0
484 {1.01064,1.01104,1.01144,1.01185,1.01227,1.0
485 {1.01103,1.01146,1.01189,1.01233,1.01279,1.0
486 {1.01139,1.01185,1.01231,1.01278,1.01327,1.0
487 {1.01173,1.01221,1.0127,1.0132,1.01372,1.014
488 {1.01205,1.01255,1.01307,1.0136,1.01414,1.01
489 {1.01235,1.01287,1.01341,1.01397,1.01454,1.0
490 {1.01262,1.01317,1.01373,1.01431,1.01491,1.0
491 {1.01288,1.01345,1.01404,1.01464,1.01526,1.0
492 {1.01312,1.01371,1.01432,1.01494,1.01559,1.0
493 {1.01335,1.01395,1.01458,1.01523,1.0159,1.01
494 {1.01356,1.01418,1.01483,1.01549,1.01618,1.0
495 {1.01375,1.01439,1.01506,1.01574,1.01646,1.0
496 {1.01393,1.01459,1.01527,1.01598,1.01671,1.0
497 {1.0141,1.01478,1.01547,1.0162,1.01695,1.017
498 {1.01426,1.01495,1.01566,1.0164,1.01717,1.01
499 {1.01441,1.01511,1.01584,1.0166,1.01738,1.01
500 {1.01455,1.01526,1.016,1.01678,1.01758,1.018
501 {1.01467,1.0154,1.01616,1.01695,1.01776,1.01
502 {1.01479,1.01553,1.0163,1.0171,1.01794,1.018
503 {1.01491,1.01566,1.01644,1.01725,1.0181,1.01
504 {1.01501,1.01577,1.01656,1.01739,1.01825,1.0
505 {1.0151,1.01587,1.01668,1.01752,1.01839,1.01
506 {1.01519,1.01597,1.01679,1.01764,1.01852,1.0
507 {1.01527,1.01606,1.01689,1.01775,1.01864,1.0
508 {1.01535,1.01615,1.01698,1.01785,1.01876,1.0
509 {1.01542,1.01622,1.01706,1.01794,1.01886,1.0
510 {1.01548,1.01629,1.01714,1.01803,1.01896,1.0
511 {1.01554,1.01636,1.01721,1.01811,1.01905,1.0
512 {1.01559,1.01641,1.01728,1.01819,1.01913,1.0
513 {1.01563,1.01647,1.01734,1.01825,1.01921,1.0
514 {1.01567,1.01651,1.01739,1.01831,1.01928,1.0
515 {1.01571,1.01655,1.01744,1.01837,1.01934,1.0
516 {1.01574,1.01659,1.01748,1.01842,1.01939,1.0
517 {1.01577,1.01662,1.01752,1.01846,1.01944,1.0
518 {1.01579,1.01665,1.01755,1.01849,1.01949,1.0
519 {1.0158,1.01667,1.01758,1.01853,1.01952,1.02
520 {1.01582,1.01669,1.0176,1.01855,1.01955,1.02
521 {1.01583,1.0167,1.01761,1.01857,1.01958,1.02
522 {1.01583,1.01671,1.01763,1.01859,1.0196,1.02
523 {1.01583,1.01671,1.01763,1.0186,1.01962,1.02
524 {1.01583,1.01671,1.01763,1.01861,1.01962,1.0
525 {1.01582,1.01671,1.01763,1.01861,1.01963,1.0
526 {1.01581,1.0167,1.01763,1.0186,1.01963,1.020
527 {1.0158,1.01668,1.01762,1.0186,1.01962,1.020
528 {1.01578,1.01667,1.0176,1.01858,1.01961,1.02
529 {1.01575,1.01665,1.01758,1.01856,1.0196,1.02
530 {1.01573,1.01662,1.01756,1.01854,1.01958,1.0
531 {1.0157,1.01659,1.01753,1.01852,1.01955,1.02
532 {1.01567,1.01656,1.0175,1.01849,1.01952,1.02
533 {1.01563,1.01652,1.01746,1.01845,1.01949,1.0
534 {1.01559,1.01648,1.01742,1.01841,1.01945,1.0
535 {1.01555,1.01644,1.01738,1.01837,1.01941,1.0
536 {1.0155,1.01639,1.01733,1.01832,1.01936,1.02
537 {1.01545,1.01634,1.01728,1.01827,1.0193,1.02
538 {1.01539,1.01628,1.01722,1.01821,1.01925,1.0
539 {1.01533,1.01622,1.01716,1.01815,1.01918,1.0
540 {1.01527,1.01616,1.0171,1.01808,1.01912,1.02
541 {1.01521,1.01609,1.01703,1.01801,1.01904,1.0
542 {1.01514,1.01602,1.01695,1.01794,1.01897,1.0
543 {1.01506,1.01595,1.01688,1.01786,1.01889,1.0
544 {1.01499,1.01587,1.0168,1.01777,1.0188,1.019
545 {1.01491,1.01579,1.01671,1.01768,1.01871,1.0
546 {1.01482,1.0157,1.01662,1.01759,1.01861,1.01
547 {1.01474,1.01561,1.01653,1.0175,1.01851,1.01
548 {1.01465,1.01552,1.01643,1.01739,1.01841,1.0
549 {1.01455,1.01542,1.01633,1.01729,1.0183,1.01
550 {1.01445,1.01531,1.01622,1.01718,1.01818,1.0
551 {1.01435,1.01521,1.01611,1.01706,1.01806,1.0
552 {1.01424,1.0151,1.01599,1.01694,1.01794,1.01
553 {1.01413,1.01498,1.01588,1.01682,1.01781,1.0
554 {1.01402,1.01486,1.01575,1.01669,1.01767,1.0
555 {1.0139,1.01474,1.01562,1.01655,1.01753,1.01
556 {1.01378,1.01461,1.01549,1.01641,1.01739,1.0
557 {1.01365,1.01448,1.01535,1.01627,1.01723,1.0
558 {1.01352,1.01434,1.01521,1.01612,1.01708,1.0
559 {1.01339,1.0142,1.01506,1.01596,1.01692,1.01
560 {1.01325,1.01406,1.01491,1.01581,1.01675,1.0
561 {1.0131,1.01391,1.01475,1.01564,1.01658,1.01
562 {1.01296,1.01375,1.01459,1.01547,1.0164,1.01
563 {1.01281,1.01359,1.01442,1.0153,1.01622,1.01
564 {1.01265,1.01343,1.01425,1.01512,1.01603,1.0
565 {1.01249,1.01326,1.01408,1.01493,1.01583,1.0
566 {1.01233,1.01309,1.0139,1.01474,1.01563,1.01
567 {1.01216,1.01292,1.01371,1.01455,1.01543,1.0
568 {1.01199,1.01274,1.01352,1.01435,1.01522,1.0
569 {1.01181,1.01255,1.01333,1.01414,1.015,1.015
570 {1.01163,1.01236,1.01313,1.01393,1.01478,1.0
571 {1.01145,1.01217,1.01292,1.01372,1.01456,1.0
572 {1.01126,1.01197,1.01271,1.0135,1.01432,1.01
573 {1.01107,1.01177,1.0125,1.01328,1.01409,1.01
574 };
575