Belle II Software development
TrgEclDataBase.cc
1/**************************************************************************
2 * basf2 (Belle II Analysis Software Framework) *
3 * Author: The Belle II Collaboration *
4 * *
5 * See git log for contributors and copyright holders. *
6 * This file is licensed under LGPL-3.0, see LICENSE.md. *
7 **************************************************************************/
8
9#include "trg/ecl/TrgEclDataBase.h"
10
11#include <TMatrixT.h>
12
13#include <iostream>
14#include <math.h>
15
16using namespace std;
17using namespace Belle2;
18
19TrgEclDataBase::TrgEclDataBase()
20{
21 _AmpCoefficient.clear();
22 _TimingCoefficient.clear();
23
24
25}
26void TrgEclDataBase::MakeFitterCoefficient(const std::vector<int>& SignalPDF, std::vector<int> NoiseMatrix)
27{
28 double DeltaTTT = 0.125;
29 const int NFitBin = 12; // The # of bin for fit
30 int Nsmalldt = 10; // # of bin between 2 Sampling point
31 int Nsmalldt2 = Nsmalldt * 2;
32 double Smalldt = DeltaTTT / Nsmalldt; // [us]
33 double SmallOffset = 1.0; // factor of table shift
34 // (NBitLutOffset)
35 // o For conversion from real to integer
36 // o Larger gives good E/T resolution.
37 // o Smaller is good, and ensitive for high E.
38 // o 12 is ok but 11 is bad
39 // int NBitLutOffset = 13;
40 int NBitLutOffset = 13;
41
42
43 //==========================
44 // noise average
45 //==========================
46 double noise_ave = 0.0; // [MeV]
47 for (int jjj = 0; jjj < 78; jjj++) {
48 noise_ave += NoiseMatrix[jjj];
49 }
50 noise_ave /= 78.;
51
52 //=========================
53 // noise matrix
54 //=========================
55 std::vector<std::vector<double>> noise_cm0 ;
56 noise_cm0.clear();
57 noise_cm0.resize(12, std::vector<double>(12, 0.0));
58 // sum noise matrix
59 for (int jjj = 0; jjj < NFitBin; jjj++) {
60 for (int kkk = 0; kkk < jjj + 1; kkk++) {
61 noise_cm0[jjj][kkk] += (NoiseMatrix[jjj] - noise_ave) * (NoiseMatrix[kkk] - noise_ave);
62 }
63 }
64
65
66
67 TMatrixT<double> noise_cm1(NFitBin, NFitBin); // NFitBin x NFitBin noise matrix
68 for (int kkk = 0; kkk < NFitBin; kkk++) {
69 for (int jjj = 0; jjj < kkk + 1; jjj++) {
70 for (int iii = 0; iii < kkk + 1; iii++) {
71 noise_cm1(jjj, iii) = noise_cm0[jjj][iii];
72 }
73 }
74 }
75 // noise_cm1.Print();
76 double determinant = noise_cm1.Determinant();
77 std::cout << "Determinant = " << determinant << std::endl;
78 if (determinant < 0.0) {
79 std::cout << "Error matrix not positive !!!!!! " << std::endl;
80 }
81 std::cout << "Noise Matrix inversion" << std::endl;
82 noise_cm1.Invert();
83 // noise_cm1.Print();
84
85
86 //===========================
87 // coefficient preparation
88 //===========================
89 // coefficient for A,B,P
90 std::vector<std::vector<double>> sg0;
91 std::vector<std::vector<double>> sg1;
92 std::vector<std::vector<double>> sg2;
93 sg0.clear();
94 sg1.clear();
95 sg2.clear();
96 sg0.resize(12, std::vector<double>(20, 0.0));
97 sg1.resize(12, std::vector<double>(20, 0.0));
98 sg2.resize(12, std::vector<double>(20, 0.0));
99 std::vector<std::vector<double>> f0;
100 std::vector<std::vector<double>> f1;
101 f0.clear();
102 f1.clear();
103 f0.resize(12, std::vector<double>(12, 0.0));
104 f1.resize(12, std::vector<double>(12, 0.0));
105
106 // coefficient for inverse covariance matrix
107 double g0g0[20] = {0}, g0g1[20] = {0}, g0g2[20] = {0};
108 double g1g1[20] = {0}, g1g2[20] = {0}, g2g2[20] = {0};
109 // determinant for inverse covariance matrix
110 // double dgg0[20] = {0};
111 // double dgg1[20] = {0};
112 double dgg2[20] = {0};
113
114 _AmpCoefficient.resize(12, std::vector<int>(10, 0));
115 _TimingCoefficient.resize(12, std::vector<int>(10, 0));
116
117 std::vector<std::vector<double>> fg31;
118 std::vector<std::vector<double>> fg32;
119 fg31.clear();
120 fg32.clear();
121 fg31.resize(20, std::vector<double>(12, 0.0));
122 fg32.resize(20, std::vector<double>(12, 0.0));
123
124 for (int kkk = 0; kkk < Nsmalldt2; kkk++) { // Nsmalldt2 loops for dt intervals
125 double ttt = -(SmallOffset + 4) * DeltaTTT + Smalldt * (kkk + 1);
126
127 for (int jjj = 0; jjj < NFitBin; jjj++) {
128 // (DeltaTTT=0.1 case)
129 // ttt ~= (-0.549 -> 0.751) [us] when kkk= 0
130 // ttt = (-0.500 -> 0.800) [us] when kkk= 47
131 // ttt = (-0.450 -> 0.850) [us] when kkk= 95
132 // ttt = (-0.400 -> 0.900) [us] when kkk=143
133 // ttt = (-0.350 -> 0.950) [us] when kkk=191
134 if (ttt > 0.0) {
135 // signal PDF calculation
136 f0[jjj][kkk] = interFADC(ttt, SignalPDF);
137 // signal PDF derivative calculation
138 double ddt = 0.005 * DeltaTTT;
139 // double ttt01 = ttt + ddt;
140 double ttt02 = ttt - ddt;
141 if (ttt > ddt) {
142 f1[jjj][kkk] = (f1[jjj][kkk] - interFADC(ttt02, SignalPDF)) / (2.0 * ddt);
143 } else {
144 f1[jjj][kkk] = f1[jjj][kkk] / (ddt + ttt);
145 }
146 } else {
147 f0[jjj][kkk] = 0.0;
148 f1[jjj][kkk] = 0.0;
149 }
150 ttt += DeltaTTT;
151
152 }
153 //
154 //
155 //
156 for (int jjj = 0; jjj < NFitBin; jjj++) {
157 double SigPDFj0 = f0[jjj][kkk];
158 double SigPDFj1 = f1[jjj][kkk];
159 for (int iii = 0; iii < NFitBin; iii++) {
160 double NoiseICM = noise_cm1(jjj, iii); // Noise Inverse Covariance Matrix
161 double SigPDFi0 = f0[iii][kkk];
162 double SigPDFi1 = f1[iii][kkk];
163 //
164 sg0[jjj][kkk] += NoiseICM * SigPDFi0;
165 sg1[jjj][kkk] += NoiseICM * SigPDFi1;
166 sg2[jjj][kkk] += NoiseICM;
167 //
168 g0g0[kkk] += NoiseICM * SigPDFj0 * SigPDFi0;
169 g0g1[kkk] += NoiseICM * SigPDFj0 * SigPDFi1;
170 g1g1[kkk] += NoiseICM * SigPDFi1 * SigPDFj1;
171 g0g2[kkk] += NoiseICM * SigPDFj0;
172 g1g2[kkk] += NoiseICM * SigPDFj1;
173 g2g2[kkk] += NoiseICM;
174 } // iii-loop end
175 } // jjj-loop end
176 //
177 // dgg0[kkk] = g0g0[kkk] * g1g1[kkk] - g0g1[kkk] * g0g1[kkk];
178 // dgg1[kkk] = g0g0[kkk] * g2g2[kkk] - g0g2[kkk] * g0g2[kkk];
179 dgg2[kkk] =
180 + g0g0[kkk] * g1g1[kkk] * g2g2[kkk]
181 - g0g1[kkk] * g0g1[kkk] * g2g2[kkk]
182 + g0g1[kkk] * g1g2[kkk] * g0g2[kkk] * 2
183 - g0g2[kkk] * g0g2[kkk] * g1g1[kkk]
184 - g1g2[kkk] * g1g2[kkk] * g0g0[kkk];
185 //
186 for (int iii = 0; iii < NFitBin; iii++) {
187 fg31[iii][kkk] = (1.0 / dgg2[kkk]) *
188 ((g1g1[kkk] * g2g2[kkk] - g1g2[kkk] * g1g2[kkk]) * sg0[iii][kkk] +
189 (g1g2[kkk] * g0g2[kkk] - g0g1[kkk] * g2g2[kkk]) * sg1[iii][kkk] +
190 (g0g1[kkk] * g1g2[kkk] - g1g1[kkk] * g0g2[kkk]) * sg2[iii][kkk]);
191 fg32[iii][kkk] = (1.0 / dgg2[kkk]) *
192 ((g1g2[kkk] * g0g2[kkk] - g0g1[kkk] * g2g2[kkk]) * sg0[iii][kkk] +
193 (g0g0[kkk] * g2g2[kkk] - g0g2[kkk] * g0g2[kkk]) * sg1[iii][kkk] +
194 (g0g1[kkk] * g0g2[kkk] - g0g0[kkk] * g1g2[kkk]) * sg2[iii][kkk]);
195 // fg33[iii][kkk] = (1.0/dgg2[kkk]) *
196 // ((g0g1[kkk]*g1g2[kkk]-g1g1[kkk]*g0g2[kkk]) * sg0[iii][kkk]+
197 // (g0g1[kkk]*g0g2[kkk]-g0g0[kkk]*g1g2[kkk]) * sg1[iii][kkk]+
198 // (g0g0[kkk]*g1g1[kkk]-g0g1[kkk]*g0g1[kkk]) * sg2[iii][kkk] );
199 }
200 } // kkk-loop end
201 //--------------------
202 // make integer
203 //--------------------
204 // int LutIntOffset = 10000; // this is correction from [us] to [0.1ns]
205 // int NBitLutOffset = 13; // this is minimum for both low and high energy.
206 int LutIntOffset = (int) pow(2, NBitLutOffset);
207 for (int kkk = 0; kkk < Nsmalldt2; kkk++) {
208 for (int iii = 0; iii < NFitBin; iii++) {
209 // int_f0[iii][kkk] = int (f0[iii][kkk]*LutIntOffset);
210 // int_f1[iii][kkk] = int (f1[iii][kkk]*LutIntOffset);
211 _AmpCoefficient[iii][kkk] = int (fg31[iii][kkk] * LutIntOffset);
212 _TimingCoefficient[iii][kkk] = int (fg32[iii][kkk] * LutIntOffset);
213 // int_fg33[iii][kkk] = int (fg33[iii][kkk]*LutIntOffset);
214 }
215 }
216
217
218
219}
220
221
222
223
224double TrgEclDataBase::GetTCFLatency(int TCId)
225{
226
227 //----------------------------------------------------------
228 // TC latency = Flight time of slowest TC - Other TCs
229 // Flight time was calculated supposing particle is gamma.
230 //----------------------------------------------------------
231
232
233 std::vector <double> TCLatency = { 0.715606, 0.417428, 4.32053e-07, 0.0101334, 0.41738, 0.715606, 0.417427, 2.63428e-06, 0.0101316, 0.417379, 0.715606, 0.417428, 0, 0.0101295, 0.41738, 0.715606, 0.417427, 5.84426e-06, 0.0101284, 0.41738, 0.715606, 0.417428, 4.32053e-07, 0.0101334, 0.41738, 0.715606, 0.417427, 2.63428e-06, 0.0101316, 0.417379, 0.715606, 0.417428, 0, 0.0101295, 0.41738, 0.715606, 0.417427, 5.84426e-06, 0.0101284, 0.41738, 0.715606, 0.417428, 4.32053e-07, 0.0101334, 0.41738, 0.715606, 0.417427, 2.63428e-06, 0.0101316, 0.417379, 0.715606, 0.417428, 0, 0.0101295, 0.41738, 0.715606, 0.417427, 5.84426e-06, 0.0101284, 0.41738, 0.715606, 0.417428, 4.32053e-07, 0.0101334, 0.41738, 0.715606, 0.417427, 2.63428e-06, 0.0101316, 0.417379, 0.715606, 0.417428, 0, 0.0101295, 0.41738, 0.715606, 0.417427, 5.84426e-06, 0.0101284, 0.41738, 0.0919657, 1.0329, 1.78857, 2.38025, 2.82433, 3.13196, 3.30927, 3.34687, 3.23671, 2.99469, 2.61992, 2.24861, 0.0919505, 1.03291, 1.78857, 2.38026, 2.82432, 3.13194, 3.30926, 3.34688, 3.23672, 2.99469, 2.61991, 2.24862, 0.0919648, 1.03289, 1.78857, 2.38026, 2.82433, 3.13195, 3.30928, 3.34687, 3.23672, 2.9947, 2.61993, 2.2486, 0.0919576, 1.0329, 1.78859, 2.38026, 2.82433, 3.13197, 3.30926, 3.34688, 3.23671, 2.99469, 2.61991, 2.2486, 0.0919594, 1.0329, 1.78858, 2.38025, 2.82433, 3.13195, 3.30927, 3.34688, 3.23671, 2.9947, 2.61992, 2.24861, 0.0919549, 1.0329, 1.78859, 2.38024, 2.82433, 3.13196, 3.30927, 3.34687, 3.23671, 2.9947, 2.61992, 2.24861, 0.0919564, 1.03289, 1.78858, 2.38024, 2.82433, 3.13197, 3.30927, 3.34688, 3.2367, 2.99469, 2.61991, 2.24861, 0.0919696, 1.0329, 1.78859, 2.38026, 2.82434, 3.13195, 3.30926, 3.34687, 3.23671, 2.99469, 2.61991, 2.24862, 0.0919613, 1.0329, 1.7886, 2.38024, 2.82432, 3.13197, 3.30925, 3.34686, 3.23672, 2.99471, 2.61992, 2.24861, 0.0919657, 1.0329, 1.78857, 2.38025, 2.82433, 3.13196, 3.30927, 3.34687, 3.23671, 2.99469, 2.61992, 2.24861, 0.0919505, 1.03291, 1.78857, 2.38026, 2.82432, 3.13194, 3.30926, 3.34688, 3.23672, 2.99469, 2.61991, 2.24862, 0.0919648, 1.03289, 1.78857, 2.38026, 2.82433, 3.13195, 3.30928, 3.34687, 3.23672, 2.9947, 2.61993, 2.2486, 0.0919576, 1.0329, 1.78859, 2.38026, 2.82433, 3.13197, 3.30926, 3.34688, 3.23671, 2.99469, 2.61991, 2.2486, 0.0919594, 1.0329, 1.78858, 2.38025, 2.82433, 3.13195, 3.30927, 3.34688, 3.23671, 2.9947, 2.61992, 2.24861, 0.0919549, 1.0329, 1.78859, 2.38024, 2.82433, 3.13196, 3.30927, 3.34687, 3.23671, 2.9947, 2.61992, 2.24861, 0.0919564, 1.03289, 1.78858, 2.38024, 2.82433, 3.13197, 3.30927, 3.34688, 3.2367, 2.99469, 2.61991, 2.24861, 0.0919696, 1.0329, 1.78859, 2.38026, 2.82434, 3.13195, 3.30926, 3.34687, 3.23671, 2.99469, 2.61991, 2.24862, 0.0919613, 1.0329, 1.7886, 2.38024, 2.82432, 3.13197, 3.30925, 3.34686, 3.23672, 2.99471, 2.61992, 2.24861, 0.0919657, 1.0329, 1.78857, 2.38025, 2.82433, 3.13196, 3.30927, 3.34687, 3.23671, 2.99469, 2.61992, 2.24861, 0.0919505, 1.03291, 1.78857, 2.38026, 2.82432, 3.13194, 3.30926, 3.34688, 3.23672, 2.99469, 2.61991, 2.24862, 0.0919648, 1.03289, 1.78857, 2.38026, 2.82433, 3.13195, 3.30928, 3.34687, 3.23672, 2.9947, 2.61993, 2.2486, 0.0919576, 1.0329, 1.78859, 2.38026, 2.82433, 3.13197, 3.30926, 3.34688, 3.23671, 2.99469, 2.61991, 2.2486, 0.0919594, 1.0329, 1.78858, 2.38025, 2.82433, 3.13195, 3.30927, 3.34688, 3.23671, 2.9947, 2.61992, 2.24861, 0.0919549, 1.0329, 1.78859, 2.38024, 2.82433, 3.13196, 3.30927, 3.34687, 3.23671, 2.9947, 2.61992, 2.24861, 0.0919564, 1.03289, 1.78858, 2.38024, 2.82433, 3.13197, 3.30927, 3.34688, 3.2367, 2.99469, 2.61991, 2.24861, 0.0919696, 1.0329, 1.78859, 2.38026, 2.82434, 3.13195, 3.30926, 3.34687, 3.23671, 2.99469, 2.61991, 2.24862, 0.0919613, 1.0329, 1.7886, 2.38024, 2.82432, 3.13197, 3.30925, 3.34686, 3.23672, 2.99471, 2.61992, 2.24861, 0.0919657, 1.0329, 1.78857, 2.38025, 2.82433, 3.13196, 3.30927, 3.34687, 3.23671, 2.99469, 2.61992, 2.24861, 0.0919505, 1.03291, 1.78857, 2.38026, 2.82432, 3.13194, 3.30926, 3.34688, 3.23672, 2.99469, 2.61991, 2.24862, 0.0919648, 1.03289, 1.78857, 2.38026, 2.82433, 3.13195, 3.30928, 3.34687, 3.23672, 2.9947, 2.61993, 2.2486, 0.0919576, 1.0329, 1.78859, 2.38026, 2.82433, 3.13197, 3.30926, 3.34688, 3.23671, 2.99469, 2.61991, 2.2486, 0.0919594, 1.0329, 1.78858, 2.38025, 2.82433, 3.13195, 3.30927, 3.34688, 3.23671, 2.9947, 2.61992, 2.24861, 0.0919549, 1.0329, 1.78859, 2.38024, 2.82433, 3.13196, 3.30927, 3.34687, 3.23671, 2.9947, 2.61992, 2.24861, 0.0919564, 1.03289, 1.78858, 2.38024, 2.82433, 3.13197, 3.30927, 3.34688, 3.2367, 2.99469, 2.61991, 2.24861, 0.0919696, 1.0329, 1.78859, 2.38026, 2.82434, 3.13195, 3.30926, 3.34687, 3.23671, 2.99469, 2.61991, 2.24862, 0.0919613, 1.0329, 1.7886, 2.38024, 2.82432, 3.13197, 3.30925, 3.34686, 3.23672, 2.99471, 2.61992, 2.24861, 3.37164, 2.58506, 2.56517, 3.37178, 3.37164, 2.58506, 2.56517, 3.37178, 3.37164, 2.58505, 2.56516, 3.37178, 3.37164, 2.58505, 2.56517, 3.37178, 3.37164, 2.58506, 2.56517, 3.37178, 3.37164, 2.58506, 2.56517, 3.37178, 3.37164, 2.58505, 2.56516, 3.37178, 3.37164, 2.58505, 2.56517, 3.37178, 3.37164, 2.58506, 2.56517, 3.37178, 3.37164, 2.58506, 2.56517, 3.37178, 3.37164, 2.58505, 2.56516, 3.37178, 3.37164, 2.58505, 2.56517, 3.37178, 3.37164, 2.58506, 2.56517, 3.37178, 3.37164, 2.58506, 2.56517, 3.37178, 3.37164, 2.58505, 2.56516, 3.37178, 3.37164, 2.58505, 2.56517, 3.37178};
234
235
236 return TCLatency[TCId - 1];
237}
238
239
240void
241TrgEclDataBase::getCoeffSigPDF(std::vector<std::vector<double>>& CoeffSigPDF0, std::vector<std::vector<double>>& CoeffSigPDF1)
242{
243 // double par_f0[20][12] ={
244
245 CoeffSigPDF0 = {
246 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.13653, 0.54824, 0.86612, 0.99103, 0.98471, 0.80303},
247 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00016, 0.17169, 0.58834, 0.88675, 0.99502, 0.97827, 0.77270},
248 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00113, 0.20964, 0.62692, 0.90539, 0.99778, 0.96963, 0.74131},
249 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00401, 0.24979, 0.66379, 0.92209, 0.99942, 0.95834, 0.70919},
250 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00994, 0.29156, 0.69881, 0.93692, 1.00000, 0.94415, 0.67666},
251 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.01980, 0.33439, 0.73187, 0.94996, 0.99960, 0.92700, 0.64398},
252 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.03415, 0.37779, 0.76290, 0.96127, 0.99829, 0.90698, 0.61140},
253 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.05317, 0.42126, 0.79184, 0.97096, 0.99614, 0.88428, 0.57915},
254 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.07679, 0.46440, 0.81870, 0.97909, 0.99323, 0.85917, 0.54742},
255 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.10472, 0.50684, 0.84345, 0.98575, 0.98951, 0.83197, 0.51637},
256 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.13653, 0.54824, 0.86612, 0.99103, 0.98471, 0.80303, 0.48613},
257 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00016, 0.17169, 0.58834, 0.88675, 0.99502, 0.97827, 0.77270, 0.45682},
258 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00113, 0.20964, 0.62692, 0.90539, 0.99778, 0.96963, 0.74131, 0.42853},
259 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00401, 0.24979, 0.66379, 0.92209, 0.99942, 0.95834, 0.70919, 0.40133},
260 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00994, 0.29156, 0.69881, 0.93692, 1.00000, 0.94415, 0.67666, 0.37526},
261 { 0.00000, 0.00000, 0.00000, 0.00000, 0.01980, 0.33439, 0.73187, 0.94996, 0.99960, 0.92700, 0.64398, 0.35035},
262 { 0.00000, 0.00000, 0.00000, 0.00000, 0.03415, 0.37779, 0.76290, 0.96127, 0.99829, 0.90698, 0.61140, 0.32664},
263 { 0.00000, 0.00000, 0.00000, 0.00000, 0.05317, 0.42126, 0.79184, 0.97096, 0.99614, 0.88428, 0.57915, 0.30411},
264 { 0.00000, 0.00000, 0.00000, 0.00000, 0.07679, 0.46440, 0.81870, 0.97909, 0.99323, 0.85917, 0.54742, 0.28278},
265 { 0.00000, 0.00000, 0.00000, 0.00000, 0.10472, 0.50684, 0.84345, 0.98575, 0.98951, 0.83197, 0.51637, 0.26261}
266 };
267
268
269 CoeffSigPDF1 = {
270 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00190, 2.68625, 3.26360, 1.73168, 0.36951, -0.44162, -2.37567},
271 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.03251, 2.93224, 3.14993, 1.56999, 0.26902, -0.59639, -2.47343},
272 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.13798, 3.13162, 3.02000, 1.41261, 0.17507, -0.79291, -2.54433},
273 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.33813, 3.28419, 2.87710, 1.26028, 0.08751, -1.01714, -2.58993},
274 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.62203, 3.39140, 2.72428, 1.11361, 0.00612, -1.25359, -2.61222},
275 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.96356, 3.45582, 2.56427, 0.97304, -0.06932, -1.48875, -2.61338},
276 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 1.33403, 3.48076, 2.39950, 0.83889, -0.13904, -1.71217, -2.59575},
277 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 1.70835, 3.46999, 2.23212, 0.71138, -0.20331, -1.91641, -2.56167},
278 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 2.06701, 3.42752, 2.06399, 0.59064, -0.26342, -2.09664, -2.51342},
279 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 2.39605, 3.35740, 1.89672, 0.47669, -0.33434, -2.25010, -2.45320},
280 { 0.00000, 0.00000, 0.00000, 0.00000, 0.00190, 2.68625, 3.26360, 1.73168, 0.36951, -0.44162, -2.37567, -2.38308},
281 { 0.00000, 0.00000, 0.00000, 0.00000, 0.03251, 2.93224, 3.14993, 1.56999, 0.26902, -0.59639, -2.47343, -2.30496},
282 { 0.00000, 0.00000, 0.00000, 0.00000, 0.13798, 3.13162, 3.02000, 1.41261, 0.17507, -0.79291, -2.54433, -2.22059},
283 { 0.00000, 0.00000, 0.00000, 0.00000, 0.33813, 3.28419, 2.87710, 1.26028, 0.08751, -1.01714, -2.58993, -2.13155},
284 { 0.00000, 0.00000, 0.00000, 0.00000, 0.62203, 3.39140, 2.72428, 1.11361, 0.00612, -1.25359, -2.61222, -2.03923},
285 { 0.00000, 0.00000, 0.00000, 0.00000, 0.96356, 3.45582, 2.56427, 0.97304, -0.06932, -1.48875, -2.61338, -1.94489},
286 { 0.00000, 0.00000, 0.00000, 0.00000, 1.33403, 3.48076, 2.39950, 0.83889, -0.13904, -1.71217, -2.59575, -1.84959},
287 { 0.00000, 0.00000, 0.00000, 0.00000, 1.70835, 3.46999, 2.23212, 0.71138, -0.20331, -1.91641, -2.56167, -1.75428},
288 { 0.00000, 0.00000, 0.00000, 0.00000, 2.06701, 3.42752, 2.06399, 0.59064, -0.26342, -2.09664, -2.51342, -1.65976},
289 { 0.00000, 0.00000, 0.00000, 0.00000, 2.39605, 3.35740, 1.89672, 0.47669, -0.33434, -2.25010, -2.45320, -1.56669}
290 };
291}
292void
293TrgEclDataBase::getCoeffNoise(int TCId, std::vector<std::vector<double>>& CoeffNoise31,
294 std::vector<std::vector<double>>& CoeffNoise32, std::vector<std::vector<double>>& CoeffNoise33)
295{
296
297 TCId = 0; //This value will be used later
298 // cppcheck-suppress knownConditionTrueFalse
299 if (TCId != 0) {return;}
300
301 CoeffNoise31 = {
302 { -0.34005, -0.16500, -0.12133, -0.14945, -0.02294, -0.27032, -0.04882, 0.07051, 0.23329, 0.31887, 0.28760, 0.20763},
303 { -0.35547, -0.16021, -0.11596, -0.15064, -0.01830, -0.29183, -0.02377, 0.08174, 0.24060, 0.31865, 0.28428, 0.19091},
304 { -0.37239, -0.15395, -0.11024, -0.15206, -0.01374, -0.31398, 0.00497, 0.09328, 0.24787, 0.31939, 0.27816, 0.17271},
305 { -0.39133, -0.14602, -0.10408, -0.15367, -0.01018, -0.33528, 0.03779, 0.10590, 0.25532, 0.32123, 0.26914, 0.15118},
306 { -0.41201, -0.13631, -0.09786, -0.15545, -0.00889, -0.35288, 0.07388, 0.11939, 0.26266, 0.32458, 0.25714, 0.12574},
307 { -0.43337, -0.12515, -0.09225, -0.15738, -0.01133, -0.36298, 0.11091, 0.13305, 0.26934, 0.32987, 0.24251, 0.09676},
308 { -0.45385, -0.11331, -0.08798, -0.15932, -0.01864, -0.36208, 0.14546, 0.14600, 0.27477, 0.33730, 0.22617, 0.06548},
309 { -0.47182, -0.10183, -0.08556, -0.16112, -0.03113, -0.34844, 0.17409, 0.15743, 0.27853, 0.34663, 0.20950, 0.03374},
310 { -0.48610, -0.09163, -0.08514, -0.16262, -0.04805, -0.32289, 0.19455, 0.16682, 0.28056, 0.35706, 0.19393, 0.00349},
311 { -0.49663, -0.08296, -0.08634, -0.16328, -0.06861, -0.28817, 0.20715, 0.17407, 0.28207, 0.36616, 0.18112, -0.02459},
312 { -0.50436, -0.07532, -0.08854, -0.16229, -0.09272, -0.24710, 0.21418, 0.17941, 0.28545, 0.36998, 0.17275, -0.05144},
313 { -0.51043, -0.06845, -0.09099, -0.16057, -0.11723, -0.20292, 0.21831, 0.18315, 0.29200, 0.36702, 0.16858, -0.07848},
314 { -0.51602, -0.06259, -0.09285, -0.16106, -0.13553, -0.15904, 0.22204, 0.18578, 0.30153, 0.35894, 0.16623, -0.10741},
315 { -0.52170, -0.05813, -0.09361, -0.16614, -0.14256, -0.11761, 0.22658, 0.18749, 0.31337, 0.34769, 0.16369, -0.13907},
316 { -0.52741, -0.05553, -0.09316, -0.17688, -0.13571, -0.08045, 0.23188, 0.18820, 0.32685, 0.33462, 0.16028, -0.17268},
317 { -0.53284, -0.05524, -0.09163, -0.19315, -0.11464, -0.04969, 0.23726, 0.18771, 0.34119, 0.32090, 0.15621, -0.20608},
318 { -0.53767, -0.05736, -0.08931, -0.21367, -0.08143, -0.02722, 0.24201, 0.18605, 0.35540, 0.30775, 0.15204, -0.23659},
319 { -0.54169, -0.06148, -0.08656, -0.23634, -0.04034, -0.01368, 0.24573, 0.18348, 0.36848, 0.29623, 0.14827, -0.26211},
320 { -0.54485, -0.06684, -0.08369, -0.25891, 0.00343, -0.00800, 0.24836, 0.18053, 0.37956, 0.28703, 0.14513, -0.28174},
321 { -0.54718, -0.07254, -0.08040, -0.28056, 0.04593, -0.00725, 0.25020, 0.17864, 0.38661, 0.28083, 0.14211, -0.29638}
322 };
323
324// double par_fg32[20][12]={
325
326 CoeffNoise32 = {
327 { -0.07680, 0.00074, 0.04672, 0.00598, 0.02151, -0.09937, 0.07807, 0.11668, 0.05597, -0.01337, -0.00295, -0.13318},
328 { -0.07820, 0.01135, 0.04339, 0.00338, 0.02373, -0.10586, 0.10140, 0.10029, 0.04966, -0.01285, -0.01141, -0.12488},
329 { -0.07873, 0.02017, 0.04065, 0.00144, 0.02418, -0.10718, 0.11882, 0.08745, 0.04387, -0.01080, -0.02133, -0.11854},
330 { -0.07876, 0.02748, 0.03832, 0.00014, 0.02204, -0.10229, 0.13109, 0.07806, 0.03874, -0.00789, -0.03124, -0.11567},
331 { -0.07820, 0.03359, 0.03604, -0.00066, 0.01712, -0.09097, 0.13884, 0.07119, 0.03404, -0.00441, -0.04056, -0.11602},
332 { -0.07662, 0.03857, 0.03348, -0.00103, 0.00954, -0.07328, 0.14202, 0.06579, 0.02942, -0.00050, -0.04891, -0.11848},
333 { -0.07351, 0.04231, 0.03042, -0.00103, -0.00023, -0.04991, 0.14031, 0.06102, 0.02462, 0.00361, -0.05593, -0.12169},
334 { -0.06856, 0.04467, 0.02687, -0.00067, -0.01136, -0.02238, 0.13371, 0.05633, 0.01954, 0.00755, -0.06126, -0.12443},
335 { -0.06185, 0.04566, 0.02302, 0.00004, -0.02284, 0.00706, 0.12294, 0.05148, 0.01431, 0.01082, -0.06472, -0.12592},
336 { -0.05384, 0.04559, 0.01917, 0.00140, -0.03417, 0.03616, 0.10946, 0.04638, 0.00984, 0.01183, -0.06573, -0.12607},
337 { -0.04527, 0.04493, 0.01560, 0.00378, -0.04534, 0.06326, 0.09507, 0.04098, 0.00750, 0.00846, -0.06357, -0.12541},
338 { -0.03692, 0.04387, 0.01270, 0.00626, -0.05416, 0.08670, 0.08150, 0.03544, 0.00754, 0.00091, -0.05908, -0.12476},
339 { -0.02945, 0.04232, 0.01081, 0.00692, -0.05696, 0.10514, 0.07005, 0.03007, 0.00922, -0.00869, -0.05435, -0.12510},
340 { -0.02302, 0.04019, 0.01001, 0.00474, -0.05209, 0.11850, 0.06100, 0.02498, 0.01176, -0.01864, -0.05067, -0.12677},
341 { -0.01742, 0.03737, 0.01013, -0.00029, -0.03978, 0.12702, 0.05386, 0.02007, 0.01461, -0.02812, -0.04828, -0.12918},
342 { -0.01229, 0.03380, 0.01092, -0.00757, -0.02131, 0.13074, 0.04792, 0.01522, 0.01731, -0.03666, -0.04685, -0.13123},
343 { -0.00739, 0.02957, 0.01214, -0.01612, 0.00124, 0.12981, 0.04260, 0.01040, 0.01944, -0.04391, -0.04597, -0.13180},
344 { -0.00254, 0.02496, 0.01355, -0.02481, 0.02538, 0.12480, 0.03757, 0.00570, 0.02068, -0.04970, -0.04532, -0.13026},
345 { 0.00234, 0.02033, 0.01500, -0.03273, 0.04877, 0.11685, 0.03273, 0.00129, 0.02077, -0.05400, -0.04478, -0.12657},
346 { 0.00727, 0.01598, 0.01669, -0.03977, 0.06995, 0.10736, 0.02806, -0.00210, 0.01862, -0.05649, -0.04449, -0.12108}
347 };
348
349 //double par_fg33[20][12] ={
350 CoeffNoise33 = {
351 { 0.23489, 0.15972, 0.15032, 0.10007, 0.11840, 0.19195, 0.06781, -0.03192, -0.02200, 0.00700, 0.00703, 0.01672},
352 { 0.24011, 0.15531, 0.15130, 0.10212, 0.11674, 0.20005, 0.05379, -0.02769, -0.02199, 0.00520, 0.00965, 0.01542},
353 { 0.24547, 0.15134, 0.15183, 0.10393, 0.11574, 0.20645, 0.04122, -0.02526, -0.02233, 0.00264, 0.01342, 0.01555},
354 { 0.25151, 0.14749, 0.15188, 0.10552, 0.11582, 0.21085, 0.02913, -0.02492, -0.02321, -0.00042, 0.01791, 0.01844},
355 { 0.25820, 0.14351, 0.15168, 0.10697, 0.11738, 0.21253, 0.01719, -0.02618, -0.02446, -0.00402, 0.02307, 0.02412},
356 { 0.26507, 0.13941, 0.15159, 0.10832, 0.12081, 0.21034, 0.00604, -0.02832, -0.02572, -0.00830, 0.02876, 0.03201},
357 { 0.27142, 0.13544, 0.15195, 0.10955, 0.12636, 0.20331, -0.00315, -0.03069, -0.02666, -0.01327, 0.03460, 0.04116},
358 { 0.27658, 0.13197, 0.15298, 0.11063, 0.13388, 0.19126, -0.00922, -0.03277, -0.02709, -0.01877, 0.04005, 0.05051},
359 { 0.28015, 0.12931, 0.15468, 0.11151, 0.14283, 0.17510, -0.01160, -0.03425, -0.02700, -0.02439, 0.04457, 0.05909},
360 { 0.28215, 0.12752, 0.15688, 0.11195, 0.15277, 0.15643, -0.01063, -0.03498, -0.02702, -0.02887, 0.04736, 0.06644},
361 { 0.28303, 0.12639, 0.15932, 0.11157, 0.16349, 0.13693, -0.00746, -0.03493, -0.02832, -0.03028, 0.04754, 0.07271},
362 { 0.28355, 0.12567, 0.16158, 0.11099, 0.17333, 0.11819, -0.00377, -0.03439, -0.03135, -0.02837, 0.04567, 0.07889},
363 { 0.28448, 0.12519, 0.16325, 0.11164, 0.17954, 0.10105, -0.00121, -0.03380, -0.03581, -0.02438, 0.04353, 0.08653},
364 { 0.28608, 0.12496, 0.16415, 0.11453, 0.18039, 0.08548, -0.00029, -0.03331, -0.04133, -0.01932, 0.04223, 0.09645},
365 { 0.28815, 0.12521, 0.16432, 0.12005, 0.17506, 0.07165, -0.00066, -0.03276, -0.04761, -0.01365, 0.04196, 0.10828},
366 { 0.29040, 0.12624, 0.16389, 0.12817, 0.16336, 0.06037, -0.00169, -0.03193, -0.05434, -0.00778, 0.04245, 0.12086},
367 { 0.29259, 0.12824, 0.16303, 0.13834, 0.14605, 0.05262, -0.00282, -0.03071, -0.06109, -0.00221, 0.04325, 0.13271},
368 { 0.29455, 0.13110, 0.16193, 0.14962, 0.12497, 0.04890, -0.00372, -0.02916, -0.06738, 0.00257, 0.04403, 0.14258},
369 { 0.29624, 0.13453, 0.16076, 0.16096, 0.10251, 0.04890, -0.00424, -0.02750, -0.07277, 0.00618, 0.04458, 0.14984},
370 { 0.29765, 0.13817, 0.15938, 0.17191, 0.08070, 0.05153, -0.00443, -0.02647, -0.07611, 0.00809, 0.04503, 0.15455}
371 };
372
373
374 return;
375}
376
377
378
379
380void TrgEclDataBase:: readNoiseLMatrix(std::vector<std::vector<double>>& MatrixParallel,
381 std::vector<std::vector<double>>& MatrixSerial)
382{
383
384 //double L
385 MatrixParallel = { { 0.00222966, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
386 { 0.00175567, 0.0013744, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
387 { 0.000978697, 0.00159799, 0.00120828, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
388 { 0.000411268, 0.00106236, 0.00150164, 0.00119109, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
389 { 0.000104886, 0.000533208, 0.00101586, 0.00149401, 0.00118828, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
390 { 1.6782e-05, 0.000148717, 0.000548641, 0.00100194, 0.00149732, 0.00118411, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
391 { 6.3743e-06, 1.90824e-05, 0.000163147, 0.000544966, 0.00100262, 0.00149888, 0.0011813, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
392 { 2.35494e-06, 7.33264e-06, 1.93629e-05, 0.000164577, 0.000544718, 0.00100488, 0.00149767, 0.00118081, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
393 { 8.45505e-07, 2.74031e-06, 7.45357e-06, 1.9282e-05, 0.000164886, 0.000545786, 0.00100482, 0.00149732, 0.00118077, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
394 { 2.96716e-07, 9.92614e-07, 2.7925e-06, 7.42397e-06, 1.92962e-05, 0.000165394, 0.000546186, 0.00100452, 0.00149734, 0.00118075, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
395 { 1.02205e-07, 3.50798e-07, 1.01349e-06, 2.78241e-06, 7.42938e-06, 1.93363e-05, 0.000165698, 0.000546122, 0.00100451, 0.00149735, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
396 { 3.46628e-08, 1.21525e-07, 3.58736e-07, 1.01011e-06, 2.78446e-06, 7.44503e-06, 1.93478e-05, 0.000165743, 0.000546108, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
397 { 1.16028e-08, 4.1411e-08, 1.24434e-07, 3.57618e-07, 1.01087e-06, 2.79044e-06, 7.44973e-06, 1.93465e-05, 0.000165746, 0.000546112, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
398 { 3.84053e-09, 1.39171e-08, 4.24471e-08, 1.24068e-07, 3.57887e-07, 1.01307e-06, 2.79233e-06, 7.44926e-06, 1.93464e-05, 0.00016575, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
399 { 1.25893e-09, 4.62224e-09, 1.42781e-08, 4.23287e-08, 1.24163e-07, 3.58676e-07, 1.01379e-06, 2.79218e-06, 7.44922e-06, 1.93466e-05, 0.000165752, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
400 { 4.09187e-10, 1.51963e-09, 4.74579e-09, 1.424e-08, 4.23611e-08, 1.24439e-07, 3.58941e-07, 1.01374e-06, 2.79217e-06, 7.44932e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
401 { 1.32005e-10, 4.95189e-10, 1.56129e-09, 4.73367e-09, 1.4251e-08, 4.2456e-08, 1.24533e-07, 3.58927e-07, 1.01374e-06, 2.79221e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
402 { 4.23033e-11, 1.60111e-10, 5.09061e-10, 1.55744e-09, 4.73733e-09, 1.42831e-08, 4.24891e-08, 1.24529e-07, 3.58925e-07, 1.01375e-06, 2.79223e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
403 { 1.34764e-11, 5.14127e-11, 1.6468e-10, 5.07849e-10, 1.55865e-09, 4.74806e-09, 1.42945e-08, 4.24878e-08, 1.24529e-07, 3.5893e-07, 1.01376e-06, 2.79223e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
404 { 4.27025e-12, 1.64076e-11, 5.29043e-11, 1.643e-10, 5.08245e-10, 1.5622e-09, 4.7519e-09, 1.42941e-08, 4.24877e-08, 1.2453e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
405 { 1.34658e-12, 5.20739e-12, 1.68905e-11, 5.27855e-11, 1.64429e-10, 5.09406e-10, 1.56348e-09, 4.75179e-09, 1.42941e-08, 4.24882e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
406 { 4.22769e-13, 1.64447e-12, 5.36263e-12, 1.68536e-11, 5.2827e-11, 1.64806e-10, 5.0983e-10, 1.56345e-09, 4.75177e-09, 1.42943e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
407 { 1.322e-13, 5.16973e-13, 1.69406e-12, 5.35117e-12, 1.68669e-11, 5.29484e-11, 1.64944e-10, 5.0982e-10, 1.56344e-09, 4.75184e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
408 { 4.11877e-14, 1.61852e-13, 5.32724e-13, 1.69052e-12, 5.3554e-12, 1.69057e-11, 5.29933e-11, 1.64941e-10, 5.09818e-10, 1.56347e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
409 { 1.2789e-14, 5.04811e-14, 1.66829e-13, 5.31633e-13, 1.69186e-12, 5.36777e-12, 1.69202e-11, 5.29925e-11, 1.64941e-10, 5.09825e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
410 { 3.95867e-15, 1.56904e-14, 5.20466e-14, 1.66494e-13, 5.32054e-13, 1.69578e-12, 5.37239e-12, 1.692e-11, 5.29923e-11, 1.64943e-10, 5.0983e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
411 { 1.22183e-15, 4.86128e-15, 1.61807e-14, 5.19439e-14, 1.66626e-13, 5.33288e-13, 1.69725e-12, 5.37232e-12, 1.69199e-11, 5.2993e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
412 { 3.76104e-16, 1.5017e-15, 5.01428e-15, 1.61493e-14, 5.19852e-14, 1.67013e-13, 5.33754e-13, 1.69723e-12, 5.37231e-12, 1.69201e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
413 { 1.15484e-16, 4.62624e-16, 1.54928e-15, 5.0047e-15, 1.61622e-14, 5.21062e-14, 1.6716e-13, 5.33748e-13, 1.69722e-12, 5.37238e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
414 { 3.53776e-17, 1.42156e-16, 4.77368e-16, 1.54636e-15, 5.00869e-15, 1.61999e-14, 5.21521e-14, 1.67158e-13, 5.33747e-13, 1.69725e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
415 { 1.08141e-17, 4.35782e-17, 1.46712e-16, 4.7648e-16, 1.54759e-15, 5.02038e-15, 1.62142e-14, 5.21517e-14, 1.67158e-13, 5.33754e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
416 { 3.29886e-18, 1.33294e-17, 4.49821e-17, 1.46442e-16, 4.76862e-16, 1.55121e-15, 5.02484e-15, 1.62141e-14, 5.21515e-14, 1.6716e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
417 { 1.0044e-18, 4.06863e-18, 1.37609e-17, 4.49005e-17, 1.4656e-16, 4.77977e-16, 1.55259e-15, 5.02481e-15, 1.6214e-14, 5.21523e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
418 { 3.05258e-19, 1.23947e-18, 4.20092e-18, 1.37362e-17, 4.49366e-17, 1.46903e-16, 4.78405e-16, 1.55259e-15, 5.0248e-15, 1.62143e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
419 { 9.26167e-20, 3.76903e-19, 1.27994e-18, 4.19347e-18, 1.37472e-17, 4.50419e-17, 1.47035e-16, 4.78403e-16, 1.55258e-15, 5.02487e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
420 { 2.80553e-20, 1.14412e-19, 3.89258e-19, 1.2777e-18, 4.19684e-18, 1.37795e-17, 4.50825e-17, 1.47034e-16, 4.78402e-16, 1.5526e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
421 { 8.48559e-21, 3.46739e-20, 1.18176e-19, 3.88582e-19, 1.27873e-18, 4.2067e-18, 1.37919e-17, 4.50823e-17, 1.47034e-16, 4.78409e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
422 { 2.56286e-21, 1.04922e-20, 3.58187e-20, 1.17973e-19, 3.88895e-19, 1.28173e-18, 4.21051e-18, 1.37919e-17, 4.50822e-17, 1.47036e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
423 { 7.72993e-22, 3.17024e-21, 1.08397e-20, 3.57576e-20, 1.18068e-19, 3.89809e-19, 1.2829e-18, 4.2105e-18, 1.37919e-17, 4.50829e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
424 { 2.32843e-22, 9.56572e-22, 3.27558e-21, 1.08214e-20, 3.57865e-20, 1.18346e-19, 3.90165e-19, 1.28289e-18, 4.21049e-18, 1.37921e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
425 { 7.00514e-23, 2.88252e-22, 9.8845e-22, 3.27008e-21, 1.08301e-20, 3.58707e-20, 1.18454e-19, 3.90164e-19, 1.28289e-18, 4.21055e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
426 { 2.10504e-23, 8.67527e-23, 2.97884e-22, 9.86805e-22, 3.27273e-21, 1.08556e-20, 3.59036e-20, 1.18454e-19, 3.90163e-19, 1.28291e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
427 { 6.31856e-24, 2.60782e-23, 8.96594e-23, 2.97392e-22, 9.87602e-22, 3.28044e-21, 1.08656e-20, 3.59035e-20, 1.18453e-19, 3.90169e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
428 { 1.89458e-24, 7.8303e-24, 2.69541e-23, 8.95124e-23, 2.97633e-22, 9.89932e-22, 3.28346e-21, 1.08656e-20, 3.59035e-20, 1.18455e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
429 { 5.67497e-25, 2.3486e-24, 8.09394e-24, 2.69102e-23, 8.95848e-23, 2.98335e-22, 9.90843e-22, 3.28346e-21, 1.08656e-20, 3.5904e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
430 { 1.69821e-25, 7.03706e-25, 2.42786e-24, 8.08084e-24, 2.6932e-23, 8.97964e-23, 2.9861e-22, 9.90844e-22, 3.28345e-21, 1.08657e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
431 { 5.07706e-26, 2.10641e-25, 7.27505e-25, 2.42395e-24, 8.08739e-24, 2.69957e-23, 8.98793e-23, 2.98611e-22, 9.90841e-22, 3.28349e-21, 1.08658e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
432 { 1.51651e-26, 6.29919e-26, 2.1778e-25, 7.26342e-25, 2.42592e-24, 8.10651e-24, 2.70206e-23, 8.98794e-23, 2.9861e-22, 9.90855e-22, 3.28353e-21, 1.08658e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
433 { 4.5259e-27, 1.88205e-26, 6.51309e-26, 2.17434e-25, 7.26931e-25, 2.43166e-24, 8.11402e-24, 2.70206e-23, 8.98792e-23, 2.98614e-22, 9.90865e-22, 3.28352e-21, 1.08658e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
434 { 1.3496e-27, 5.61826e-27, 1.94608e-26, 6.50279e-26, 2.1761e-25, 7.28651e-25, 2.43391e-24, 8.11403e-24, 2.70206e-23, 8.98805e-23, 2.98617e-22, 9.90864e-22, 3.28352e-21, 1.08658e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
435 { 4.02128e-28, 1.67575e-27, 5.80974e-27, 1.94302e-26, 6.50807e-26, 2.18125e-25, 7.29327e-25, 2.43392e-24, 8.11401e-24, 2.7021e-23, 8.98813e-23, 2.98617e-22, 9.90864e-22, 3.28352e-21, 1.08658e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
436 { 1.19727e-28, 4.99423e-28, 1.73296e-27, 5.80066e-27, 1.9446e-26, 6.52349e-26, 2.18328e-25, 7.29329e-25, 2.43391e-24, 8.11413e-24, 2.70212e-23, 8.98813e-23, 2.98617e-22, 9.90864e-22, 3.28352e-21, 1.08658e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
437 { 3.56205e-29, 1.48728e-28, 5.16502e-28, 1.73027e-27, 5.80537e-27, 1.94921e-26, 6.52956e-26, 2.18329e-25, 7.29327e-25, 2.43394e-24, 8.1142e-24, 2.70212e-23, 8.98813e-23, 2.98617e-22, 9.90864e-22, 3.28352e-21, 1.08658e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
438 { 1.05902e-29, 4.42583e-29, 1.53822e-28, 5.15702e-28, 1.73167e-27, 5.81913e-27, 1.95102e-26, 6.52957e-26, 2.18328e-25, 7.29338e-25, 2.43397e-24, 8.1142e-24, 2.70212e-23, 8.98813e-23, 2.98617e-22, 9.90864e-22, 3.28352e-21, 1.08658e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
439 { 3.14637e-30, 1.3161e-29, 4.57766e-29, 1.53585e-28, 5.16122e-28, 1.73578e-27, 5.82456e-27, 1.95103e-26, 6.52956e-26, 2.18331e-25, 7.29345e-25, 2.43397e-24, 8.1142e-24, 2.70212e-23, 8.98813e-23, 2.98617e-22, 9.90864e-22, 3.28352e-21, 1.08658e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
440 { 9.34181e-31, 3.91094e-30, 1.36131e-29, 4.57064e-29, 1.5371e-28, 5.17346e-28, 1.7374e-27, 5.82457e-27, 1.95103e-26, 6.52965e-26, 2.18333e-25, 7.29344e-25, 2.43397e-24, 8.1142e-24, 2.70212e-23, 8.98813e-23, 2.98617e-22, 9.90864e-22, 3.28352e-21, 1.08658e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, 0, },
441 { 2.7719e-31, 1.16141e-30, 4.0455e-30, 1.35923e-29, 4.57436e-29, 1.54075e-28, 5.17829e-28, 1.7374e-27, 5.82456e-27, 1.95105e-26, 6.52971e-26, 2.18333e-25, 7.29344e-25, 2.43397e-24, 8.1142e-24, 2.70212e-23, 8.98813e-23, 2.98617e-22, 9.90864e-22, 3.28352e-21, 1.08658e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, 0, },
442 { 8.21977e-32, 3.44679e-31, 1.20143e-30, 4.03935e-30, 1.36034e-29, 4.58522e-29, 1.54219e-28, 5.17831e-28, 1.7374e-27, 5.82465e-27, 1.95107e-26, 6.52971e-26, 2.18333e-25, 7.29344e-25, 2.43397e-24, 8.1142e-24, 2.70212e-23, 8.98813e-23, 2.98617e-22, 9.90864e-22, 3.28352e-21, 1.08658e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, 0, },
443 { 2.43604e-32, 1.02229e-31, 3.5657e-31, 1.19961e-30, 4.04264e-30, 1.36357e-29, 4.58951e-29, 1.5422e-28, 5.1783e-28, 1.73742e-27, 5.8247e-27, 1.95107e-26, 6.52971e-26, 2.18333e-25, 7.29344e-25, 2.43397e-24, 8.1142e-24, 2.70212e-23, 8.98813e-23, 2.98617e-22, 9.90864e-22, 3.28352e-21, 1.08658e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, 0, },
444 { 7.21544e-33, 3.03023e-32, 1.0576e-31, 3.56032e-31, 1.20059e-30, 4.05223e-30, 1.36485e-29, 4.58953e-29, 1.54219e-28, 5.17837e-28, 1.73744e-27, 5.8247e-27, 1.95107e-26, 6.52971e-26, 2.18333e-25, 7.29344e-25, 2.43397e-24, 8.1142e-24, 2.70212e-23, 8.98813e-23, 2.98617e-22, 9.90864e-22, 3.28352e-21, 1.08658e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, 0, },
445 { 2.13599e-33, 8.97689e-33, 3.13503e-32, 1.05601e-31, 3.56322e-31, 1.20344e-30, 4.05603e-30, 1.36485e-29, 4.58952e-29, 1.54221e-28, 5.17842e-28, 1.73744e-27, 5.8247e-27, 1.95107e-26, 6.52971e-26, 2.18333e-25, 7.29344e-25, 2.43397e-24, 8.1142e-24, 2.70212e-23, 8.98813e-23, 2.98617e-22, 9.90864e-22, 3.28352e-21, 1.08658e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, 0, },
446 { 6.31983e-34, 2.65787e-33, 9.28771e-33, 3.13033e-32, 1.05687e-31, 3.57168e-31, 1.20457e-30, 4.05605e-30, 1.36485e-29, 4.58958e-29, 1.54223e-28, 5.17842e-28, 1.73744e-27, 5.8247e-27, 1.95107e-26, 6.52971e-26, 2.18333e-25, 7.29344e-25, 2.43397e-24, 8.1142e-24, 2.70212e-23, 8.98813e-23, 2.98617e-22, 9.90864e-22, 3.28352e-21, 1.08658e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, 0, },
447 { 1.8689e-34, 7.86516e-34, 2.75e-33, 9.27386e-33, 3.13289e-32, 1.05939e-31, 3.57504e-31, 1.20457e-30, 4.05604e-30, 1.36487e-29, 4.58963e-29, 1.54223e-28, 5.17842e-28, 1.73744e-27, 5.8247e-27, 1.95107e-26, 6.52971e-26, 2.18333e-25, 7.29344e-25, 2.43397e-24, 8.1142e-24, 2.70212e-23, 8.98813e-23, 2.98617e-22, 9.90864e-22, 3.28352e-21, 1.08658e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072, 0, },
448 { 5.52395e-35, 2.32624e-34, 8.1381e-34, 2.74592e-33, 9.28142e-33, 3.14033e-32, 1.06038e-31, 3.57506e-31, 1.20457e-30, 4.0561e-30, 1.36488e-29, 4.58963e-29, 1.54223e-28, 5.17842e-28, 1.73744e-27, 5.8247e-27, 1.95107e-26, 6.52971e-26, 2.18333e-25, 7.29344e-25, 2.43397e-24, 8.1142e-24, 2.70212e-23, 8.98813e-23, 2.98617e-22, 9.90864e-22, 3.28352e-21, 1.08658e-20, 3.59043e-20, 1.18456e-19, 3.90172e-19, 1.28292e-18, 4.21059e-18, 1.37922e-17, 4.50833e-17, 1.47037e-16, 4.78413e-16, 1.55262e-15, 5.02491e-15, 1.62144e-14, 5.21527e-14, 1.67162e-13, 5.33759e-13, 1.69726e-12, 5.37243e-12, 1.69203e-11, 5.29935e-11, 1.64944e-10, 5.09829e-10, 1.56348e-09, 4.75188e-09, 1.42944e-08, 4.24886e-08, 1.24531e-07, 3.58933e-07, 1.01376e-06, 2.79222e-06, 7.44937e-06, 1.93467e-05, 0.000165753, 0.000546116, 0.00100453, 0.00149734, 0.00118072 }
449 }
450 ;
451 //double L
452 MatrixSerial = { { 0.013249, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
453 { 0.00475532, 0.0123662, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
454 { 0.00274508, 0.0040392, 0.012316, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
455 { 0.00210822, 0.00213035, 0.00394697, 0.0122871, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
456 { 0.000376401, 0.00211398, 0.00217582, 0.00399753, 0.012256, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
457 { 1.98224e-05, 0.000395649, 0.00213374, 0.00220255, 0.00397452, 0.0122547, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
458 { 0.00059039, -0.000205793, 0.000340816, 0.00209816, 0.00223999, 0.0039839, 0.0122367, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
459 { 0.00165779, -4.95393e-06, -0.000346551, 0.000233604, 0.00221429, 0.00226548, 0.00389532, 0.012141, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
460 { 0.00126864, 0.00128829, -7.01616e-05, -0.000397125, 0.000287709, 0.00222591, 0.00222532, 0.00384044, 0.0121394, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
461 { 0.000687338, 0.00109489, 0.00127108, -7.94679e-05, -0.00038827, 0.000289372, 0.00222296, 0.00224361, 0.00384555, 0.0121344, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
462 { 0.000317359, 0.000614368, 0.00109251, 0.00127564, -8.75235e-05, -0.000390029, 0.000296407, 0.00225782, 0.00225064, 0.00383409, 0.0121295, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
463 { 0.00013335, 0.000288736, 0.000614988, 0.00109745, 0.00127107, -8.78134e-05, -0.000389504, 0.000298635, 0.0022587, 0.00224948, 0.0038354, 0.0121291, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
464 { 5.2629e-05, 0.000122631, 0.00028945, 0.000617872, 0.00109573, 0.00127143, -8.97687e-05, -0.000397417, 0.000297259, 0.00226201, 0.00225458, 0.00383378, 0.0121278, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
465 { 1.98609e-05, 4.87488e-05, 0.000123036, 0.00029082, 0.000617307, 0.00109606, 0.00127165, -8.56127e-05, -0.000397023, 0.000296772, 0.00226256, 0.00225423, 0.00383371, 0.0121277, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
466 { 7.24665e-06, 1.84921e-05, 4.89358e-05, 0.000123619, 0.00029065, 0.000617492, 0.00109675, 0.001287, -8.23599e-05, -0.000402957, 0.000290846, 0.00226433, 0.00225779, 0.00383321, 0.0121252, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
467 { 2.57546e-06, 6.7736e-06, 1.85699e-05, 4.91674e-05, 0.000123572, 0.000290737, 0.000617964, 0.00110859, 0.00128981, -8.76795e-05, -0.000408429, 0.000292502, 0.00226756, 0.00225733, 0.0038293, 0.0121242, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
468 { 8.96199e-07, 2.41469e-06, 6.8039e-06, 1.86576e-05, 4.91554e-05, 0.000123609, 0.00029098, 0.000624391, 0.00111019, 0.00128724, -8.97364e-05, -0.000407831, 0.000293576, 0.00226741, 0.00225637, 0.0038293, 0.0121242, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
469 { 3.06502e-07, 8.42315e-07, 2.42598e-06, 6.83599e-06, 1.86549e-05, 4.917e-05, 0.000123717, 0.000293947, 0.000625148, 0.00110915, 0.00128718, -8.98748e-05, -0.000408129, 0.000293623, 0.00226838, 0.00225742, 0.00382911, 0.0121239, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
470 { 1.03318e-07, 2.88652e-07, 8.46391e-07, 2.4374e-06, 6.83552e-06, 1.86605e-05, 4.92145e-05, 0.000124964, 0.000294271, 0.000624761, 0.00110947, 0.001287, -9.01261e-05, -0.000408095, 0.00029403, 0.00226884, 0.00225732, 0.00382887, 0.0121239, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
471 { 3.4402e-08, 9.74643e-08, 2.90087e-07, 8.50369e-07, 2.43738e-06, 6.83755e-06, 1.86777e-05, 4.97063e-05, 0.000125093, 0.000294133, 0.000624996, 0.00110937, 0.00128703, -9.01484e-05, -0.000408414, 0.000293807, 0.00226888, 0.00225754, 0.00382893, 0.0121238, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
472 { 1.13347e-08, 3.24992e-08, 9.79591e-08, 2.91448e-07, 8.50401e-07, 2.4381e-06, 6.84398e-06, 1.88632e-05, 4.97555e-05, 0.000125045, 0.000294255, 0.000624947, 0.00110944, 0.001287, -9.03853e-05, -0.000408601, 0.000293842, 0.00226908, 0.0022576, 0.00382881, 0.0121238, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
473 { 3.70048e-09, 1.07209e-08, 3.26671e-08, 9.84181e-08, 2.9147e-07, 8.50652e-07, 2.44042e-06, 6.91165e-06, 1.88813e-05, 4.97394e-05, 0.0001251, 0.000294234, 0.000624997, 0.00110943, 0.00128716, -9.01351e-05, -0.000408641, 0.000293698, 0.00226905, 0.00225766, 0.00382887, 0.0121237, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
474 { 1.19844e-09, 3.5038e-09, 1.07771e-08, 3.282e-08, 9.84289e-08, 2.91556e-07, 8.5147e-07, 2.46447e-06, 6.91811e-06, 1.8876e-05, 4.97619e-05, 0.000125092, 0.00029426, 0.00062499, 0.00110961, 0.00128751, -9.02068e-05, -0.000408924, 0.00029364, 0.00226921, 0.00225776, 0.00382875, 0.0121237, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
475 { 3.85382e-10, 1.1358e-09, 3.52239e-09, 1.08275e-08, 3.28245e-08, 9.84579e-08, 2.91839e-07, 8.59835e-07, 2.46672e-06, 6.91636e-06, 1.88847e-05, 4.97588e-05, 0.000125103, 0.000294257, 0.0006251, 0.00110984, 0.00128746, -9.03819e-05, -0.00040896, 0.000293728, 0.00226926, 0.0022577, 0.00382869, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
476 { 1.23146e-10, 3.65539e-10, 1.14189e-09, 3.53884e-09, 1.08292e-08, 3.28341e-08, 9.8554e-08, 2.94699e-07, 8.6061e-07, 2.46616e-06, 6.9196e-06, 1.88836e-05, 4.97635e-05, 0.000125102, 0.00029431, 0.000625218, 0.00110982, 0.00128741, -9.03843e-05, -0.000408957, 0.00029373, 0.00226927, 0.00225773, 0.0038287, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
477 { 3.9128e-11, 1.16891e-10, 3.67517e-10, 1.14722e-09, 3.53948e-09, 1.08324e-08, 3.28664e-08, 9.9518e-08, 2.94961e-07, 8.60429e-07, 2.46733e-06, 6.9192e-06, 1.88854e-05, 4.9763e-05, 0.000125125, 0.000294363, 0.000625207, 0.00110981, 0.00128742, -9.0404e-05, -0.000408969, 0.000293747, 0.00226931, 0.00225774, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
478 { 1.23691e-11, 3.71649e-11, 1.17529e-10, 3.6923e-10, 1.14745e-09, 3.54052e-09, 1.08431e-08, 3.31873e-08, 9.96054e-08, 2.94904e-07, 8.60839e-07, 2.46719e-06, 6.91987e-06, 1.88852e-05, 4.97723e-05, 0.000125147, 0.000294359, 0.00062521, 0.00110982, 0.00128741, -9.04038e-05, -0.00040897, 0.000293745, 0.00226931, 0.00225774, 0.00382869, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
479 { 3.89207e-12, 1.17555e-11, 3.73691e-11, 1.18076e-10, 3.6931e-10, 1.14778e-09, 3.54403e-09, 1.09488e-08, 3.32162e-08, 9.95873e-08, 2.95045e-07, 8.60793e-07, 2.46743e-06, 6.91982e-06, 1.88888e-05, 4.97808e-05, 0.000125145, 0.000294361, 0.000625213, 0.00110981, 0.00128742, -9.04147e-05, -0.000408988, 0.00029374, 0.00226932, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
480 { 1.21953e-12, 3.70095e-12, 1.18205e-11, 3.7543e-11, 1.18103e-10, 3.69418e-10, 1.14893e-09, 3.57854e-09, 1.09582e-08, 3.32105e-08, 9.96353e-08, 2.9503e-07, 8.60878e-07, 2.46741e-06, 6.92113e-06, 1.88919e-05, 4.97801e-05, 0.000125147, 0.000294363, 0.000625213, 0.00110982, 0.00128741, -9.04157e-05, -0.000408989, 0.000293741, 0.00226932, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
481 { 3.80658e-13, 1.16021e-12, 3.72154e-12, 1.18754e-11, 3.75522e-11, 1.18138e-10, 3.69787e-10, 1.1601e-09, 3.58161e-09, 1.09565e-08, 3.32266e-08, 9.96301e-08, 2.95059e-07, 8.60871e-07, 2.46788e-06, 6.92227e-06, 1.88917e-05, 4.97808e-05, 0.000125147, 0.000294363, 0.000625216, 0.00110982, 0.00128742, -9.04109e-05, -0.000408996, 0.000293736, 0.00226932, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
482 { 1.18398e-13, 3.62304e-13, 1.1667e-12, 3.73883e-12, 1.18785e-11, 3.75631e-11, 1.18256e-10, 3.7338e-10, 1.16109e-09, 3.58105e-09, 1.09618e-08, 3.32249e-08, 9.96401e-08, 2.95057e-07, 8.61036e-07, 2.46828e-06, 6.92219e-06, 1.8892e-05, 4.9781e-05, 0.000125148, 0.000294364, 0.000625216, 0.00110983, 0.00128743, -9.04186e-05, -0.000409002, 0.000293738, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
483 { 3.67061e-14, 1.12735e-13, 3.64339e-13, 1.17212e-12, 3.73983e-12, 1.1882e-11, 3.76008e-11, 1.19404e-10, 3.73696e-10, 1.16092e-09, 3.5828e-09, 1.09612e-08, 3.32282e-08, 9.96394e-08, 2.95113e-07, 8.61175e-07, 2.46826e-06, 6.9223e-06, 1.88921e-05, 4.97811e-05, 0.000125148, 0.000294364, 0.000625221, 0.00110983, 0.00128742, -9.04205e-05, -0.000409002, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
484 { 1.13456e-14, 3.49636e-14, 1.13371e-13, 3.66029e-13, 1.17244e-12, 3.74092e-12, 1.18939e-11, 3.79656e-11, 1.19504e-10, 3.73643e-10, 1.16149e-09, 3.58262e-09, 1.09623e-08, 3.3228e-08, 9.96585e-08, 2.95161e-07, 8.61166e-07, 2.4683e-06, 6.92232e-06, 1.88921e-05, 4.97813e-05, 0.000125148, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04209e-05, -0.000409002, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
485 { 3.49711e-15, 1.08108e-14, 3.51616e-14, 1.13896e-13, 3.66134e-13, 1.17278e-12, 3.74469e-12, 1.20092e-11, 3.79974e-11, 1.19488e-10, 3.73827e-10, 1.16143e-09, 3.58299e-09, 1.09623e-08, 3.32344e-08, 9.96744e-08, 2.95157e-07, 8.6118e-07, 2.46831e-06, 6.92234e-06, 1.88922e-05, 4.97813e-05, 0.000125149, 0.000294367, 0.000625221, 0.00110983, 0.00128742, -9.04214e-05, -0.000409002, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
486 { 1.07514e-15, 3.33332e-15, 1.08722e-14, 3.53246e-14, 1.1393e-13, 3.66241e-13, 1.17397e-12, 3.78097e-12, 1.20192e-11, 3.79924e-11, 1.19547e-10, 3.73808e-10, 1.16155e-09, 3.58296e-09, 1.09644e-08, 3.32397e-08, 9.96734e-08, 2.95163e-07, 8.61184e-07, 2.46831e-06, 6.92236e-06, 1.88922e-05, 4.97817e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04214e-05, -0.000409003, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
487 { 3.29743e-16, 1.0251e-15, 3.35233e-15, 1.09226e-14, 3.53352e-14, 1.13963e-13, 3.66611e-13, 1.18534e-12, 3.78411e-12, 1.20177e-11, 3.80111e-11, 1.19541e-10, 3.73846e-10, 1.16154e-09, 3.58365e-09, 1.09661e-08, 3.32393e-08, 9.96751e-08, 2.95164e-07, 8.61185e-07, 2.46832e-06, 6.92237e-06, 1.88923e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04219e-05, -0.000409003, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
488 { 1.00904e-16, 3.14481e-16, 1.03096e-15, 3.36785e-15, 1.0926e-14, 3.53455e-14, 1.14079e-13, 3.70159e-13, 1.18632e-12, 3.78364e-12, 1.20236e-11, 3.80093e-11, 1.19553e-10, 3.73844e-10, 1.16176e-09, 3.58422e-09, 1.0966e-08, 3.32399e-08, 9.96755e-08, 2.95164e-07, 8.61189e-07, 2.46832e-06, 6.92242e-06, 1.88923e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04217e-05, -0.000409003, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
489 { 3.08122e-17, 9.62584e-17, 3.16284e-16, 1.03573e-15, 3.36891e-15, 1.09291e-14, 3.53814e-14, 1.15182e-13, 3.70464e-13, 1.18617e-12, 3.78551e-12, 1.2023e-11, 3.80132e-11, 1.19552e-10, 3.73916e-10, 1.16194e-09, 3.58418e-09, 1.09662e-08, 3.324e-08, 9.96758e-08, 2.95165e-07, 8.61189e-07, 2.46834e-06, 6.92243e-06, 1.88923e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04219e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
490 { 9.39029e-18, 2.94009e-17, 9.68118e-17, 3.17748e-16, 1.03606e-15, 3.36989e-15, 1.09403e-14, 3.57234e-14, 1.15277e-13, 3.70421e-13, 1.18676e-12, 3.78533e-12, 1.20243e-11, 3.80129e-11, 1.19576e-10, 3.73975e-10, 1.16193e-09, 3.58425e-09, 1.09663e-08, 3.32401e-08, 9.96761e-08, 2.95166e-07, 8.61196e-07, 2.46835e-06, 6.92243e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
491 { 2.85646e-18, 8.96222e-18, 2.95703e-17, 9.72596e-17, 3.17852e-16, 1.03636e-15, 3.37332e-15, 1.1046e-14, 3.57528e-14, 1.15264e-13, 3.70605e-13, 1.18671e-12, 3.78572e-12, 1.20242e-11, 3.80203e-11, 1.19594e-10, 3.73971e-10, 1.16195e-09, 3.58426e-09, 1.09663e-08, 3.32402e-08, 9.96762e-08, 2.95168e-07, 8.61197e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
492 { 8.67394e-19, 2.72683e-18, 9.01398e-18, 2.9707e-17, 9.72919e-17, 3.17944e-16, 1.03742e-15, 3.40591e-15, 1.1055e-14, 3.57488e-14, 1.15321e-13, 3.70587e-13, 1.18683e-12, 3.7857e-12, 1.20265e-11, 3.80262e-11, 1.19593e-10, 3.73978e-10, 1.16196e-09, 3.58427e-09, 1.09663e-08, 3.32403e-08, 9.9677e-08, 2.95169e-07, 8.61197e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
493 { 2.62958e-19, 8.28197e-19, 2.74261e-18, 9.05565e-18, 2.9717e-17, 9.73202e-17, 3.18269e-16, 1.04744e-15, 3.40869e-15, 1.10538e-14, 3.57666e-14, 1.15315e-13, 3.70625e-13, 1.18682e-12, 3.78643e-12, 1.20284e-11, 3.80259e-11, 1.19595e-10, 3.73979e-10, 1.16196e-09, 3.58428e-09, 1.09663e-08, 3.32405e-08, 9.96771e-08, 2.95168e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
494 { 7.95939e-20, 2.51123e-19, 8.32999e-19, 2.75528e-18, 9.05873e-18, 2.97257e-17, 9.74196e-17, 3.21341e-16, 1.04829e-15, 3.40833e-15, 1.10593e-14, 3.57649e-14, 1.15327e-13, 3.70623e-13, 1.18705e-12, 3.78702e-12, 1.20283e-11, 3.80266e-11, 1.19596e-10, 3.7398e-10, 1.16197e-09, 3.58429e-09, 1.09664e-08, 3.32406e-08, 9.96771e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
495 { 2.40564e-20, 7.60253e-20, 2.52582e-19, 8.36848e-19, 2.75623e-18, 9.06137e-18, 2.9756e-17, 9.83596e-17, 3.21602e-16, 1.04818e-15, 3.41002e-15, 1.10588e-14, 3.57686e-14, 1.15327e-13, 3.70695e-13, 1.18723e-12, 3.78698e-12, 1.20285e-11, 3.80267e-11, 1.19596e-10, 3.73982e-10, 1.16197e-09, 3.58431e-09, 1.09664e-08, 3.32406e-08, 9.96772e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
496 { 7.26062e-21, 2.29817e-20, 7.64677e-20, 2.53749e-19, 8.3714e-19, 2.75704e-18, 9.07064e-18, 3.00431e-17, 9.84395e-17, 3.21569e-16, 1.0487e-15, 3.40987e-15, 1.10599e-14, 3.57683e-14, 1.15349e-13, 3.70753e-13, 1.18722e-12, 3.78706e-12, 1.20286e-11, 3.80268e-11, 1.19597e-10, 3.73982e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96772e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
497 { 2.18846e-21, 6.93739e-21, 2.31157e-20, 7.68208e-20, 2.53838e-19, 8.37384e-19, 2.75986e-18, 9.15812e-18, 3.00674e-17, 9.84296e-17, 3.21729e-16, 1.04865e-15, 3.41022e-15, 1.10599e-14, 3.57753e-14, 1.15367e-13, 3.70749e-13, 1.18725e-12, 3.78707e-12, 1.20286e-11, 3.8027e-11, 1.19597e-10, 3.73985e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96772e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
498 { 6.58803e-22, 2.09136e-21, 6.97789e-21, 2.32224e-20, 7.68482e-20, 2.53912e-19, 8.38242e-19, 2.78647e-18, 9.16553e-18, 3.00645e-17, 9.84787e-17, 3.21715e-16, 1.04876e-15, 3.4102e-15, 1.1062e-14, 3.57808e-14, 1.15366e-13, 3.70756e-13, 1.18725e-12, 3.78708e-12, 1.20286e-11, 3.8027e-11, 1.19598e-10, 3.73986e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96772e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
499 { 1.98084e-22, 6.29662e-22, 2.10358e-21, 7.01009e-21, 2.32308e-20, 7.68706e-20, 2.54173e-19, 8.46322e-19, 2.78872e-18, 9.16464e-18, 3.00795e-17, 9.84742e-17, 3.21748e-16, 1.04876e-15, 3.41086e-15, 1.10637e-14, 3.57805e-14, 1.15368e-13, 3.70758e-13, 1.18725e-12, 3.78709e-12, 1.20286e-11, 3.80273e-11, 1.19598e-10, 3.73985e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96772e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
500 { 5.94904e-23, 1.89349e-22, 6.33349e-22, 2.11329e-21, 7.01264e-21, 2.32375e-20, 7.69494e-20, 2.56622e-19, 8.47004e-19, 2.78845e-18, 9.16922e-18, 3.00781e-17, 9.84844e-17, 3.21746e-16, 1.04896e-15, 3.41139e-15, 1.10636e-14, 3.57812e-14, 1.15368e-13, 3.70759e-13, 1.18726e-12, 3.7871e-12, 1.20287e-11, 3.80274e-11, 1.19598e-10, 3.73986e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96773e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
501 { 1.78471e-23, 5.68743e-23, 1.90459e-22, 6.3627e-22, 2.11406e-21, 7.01468e-21, 2.32614e-20, 7.76908e-20, 2.56829e-19, 8.46925e-19, 2.78985e-18, 9.1688e-18, 3.00812e-17, 9.84837e-17, 3.21809e-16, 1.04912e-15, 3.41136e-15, 1.10638e-14, 3.57813e-14, 1.15369e-13, 3.7076e-13, 1.18726e-12, 3.78713e-12, 1.20288e-11, 3.80274e-11, 1.19598e-10, 3.73986e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96773e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
502 { 5.34856e-24, 1.70645e-23, 5.72081e-23, 1.91337e-22, 6.36505e-22, 2.11468e-21, 7.02188e-21, 2.34855e-20, 7.77533e-20, 2.56805e-19, 8.47348e-19, 2.78972e-18, 9.16975e-18, 3.0081e-17, 9.8503e-17, 3.21858e-16, 1.04911e-15, 3.41142e-15, 1.10639e-14, 3.57814e-14, 1.15369e-13, 3.7076e-13, 1.18727e-12, 3.78713e-12, 1.20288e-11, 3.80274e-11, 1.19598e-10, 3.73986e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96773e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
503 { 1.6013e-24, 5.11462e-24, 1.71647e-23, 5.74719e-23, 1.91408e-22, 6.3669e-22, 2.11685e-21, 7.08951e-21, 2.35043e-20, 7.77462e-20, 2.56933e-19, 8.47309e-19, 2.79001e-18, 9.16969e-18, 3.00869e-17, 9.85181e-17, 3.21855e-16, 1.04913e-15, 3.41144e-15, 1.10639e-14, 3.57815e-14, 1.15369e-13, 3.70763e-13, 1.18727e-12, 3.78713e-12, 1.20288e-11, 3.80274e-11, 1.19598e-10, 3.73986e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96773e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
504 { 4.78953e-25, 1.53143e-24, 5.14471e-24, 1.72439e-23, 5.74935e-23, 1.91464e-22, 6.37345e-22, 2.13723e-21, 7.0952e-21, 2.35022e-20, 7.77851e-20, 2.56922e-19, 8.47397e-19, 2.78999e-18, 9.17148e-18, 3.00915e-17, 9.85172e-17, 3.21861e-16, 1.04914e-15, 3.41144e-15, 1.10639e-14, 3.57816e-14, 1.1537e-13, 3.70764e-13, 1.18727e-12, 3.78714e-12, 1.20288e-11, 3.80274e-11, 1.19598e-10, 3.73986e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96773e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
505 { 1.43125e-25, 4.58107e-25, 1.54045e-24, 5.16843e-24, 1.72504e-23, 5.75102e-23, 1.91661e-22, 6.43481e-22, 2.13895e-21, 7.09457e-21, 2.3514e-20, 7.77816e-20, 2.56948e-19, 8.47392e-19, 2.79054e-18, 9.17289e-18, 3.00912e-17, 9.85191e-17, 3.21863e-16, 1.04914e-15, 3.41146e-15, 1.1064e-14, 3.57819e-14, 1.1537e-13, 3.70764e-13, 1.18727e-12, 3.78714e-12, 1.20288e-11, 3.80274e-11, 1.19598e-10, 3.73986e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96773e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
506 { 4.27325e-26, 1.3691e-25, 4.60808e-25, 1.54755e-24, 5.17039e-24, 1.72554e-23, 5.75694e-23, 1.93506e-22, 6.43996e-22, 2.13876e-21, 7.09812e-21, 2.35129e-20, 7.77897e-20, 2.56947e-19, 8.47558e-19, 2.79096e-18, 9.1728e-18, 3.00918e-17, 9.85195e-17, 3.21864e-16, 1.04914e-15, 3.41146e-15, 1.1064e-14, 3.57819e-14, 1.1537e-13, 3.70764e-13, 1.18727e-12, 3.78714e-12, 1.20288e-11, 3.80274e-11, 1.19598e-10, 3.73986e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96773e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, 0, },
507 { 1.27478e-26, 4.0881e-26, 1.37718e-25, 4.62931e-25, 1.54814e-24, 5.17189e-24, 1.72732e-23, 5.81234e-23, 1.9366e-22, 6.43941e-22, 2.13983e-21, 7.0978e-21, 2.35154e-20, 7.77892e-20, 2.56997e-19, 8.47687e-19, 2.79094e-18, 9.17298e-18, 3.0092e-17, 9.85198e-17, 3.21865e-16, 1.04914e-15, 3.41149e-15, 1.10641e-14, 3.57819e-14, 1.1537e-13, 3.70764e-13, 1.18727e-12, 3.78714e-12, 1.20288e-11, 3.80274e-11, 1.19598e-10, 3.73986e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96773e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, 0, },
508 { 3.7998e-27, 1.21967e-26, 4.11225e-26, 1.38353e-25, 4.63109e-25, 1.54859e-24, 5.17722e-24, 1.74394e-23, 5.81699e-23, 1.93644e-22, 6.44263e-22, 2.13973e-21, 7.09854e-21, 2.35152e-20, 7.78044e-20, 2.57036e-19, 8.4768e-19, 2.79099e-18, 9.17302e-18, 3.0092e-17, 9.85201e-17, 3.21865e-16, 1.04915e-15, 3.41149e-15, 1.10641e-14, 3.57819e-14, 1.1537e-13, 3.70764e-13, 1.18727e-12, 3.78714e-12, 1.20288e-11, 3.80274e-11, 1.19598e-10, 3.73986e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96773e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, 0, },
509 { 1.13175e-27, 3.63586e-27, 1.22688e-26, 4.1312e-26, 1.38406e-25, 4.63243e-25, 1.55019e-24, 5.22703e-24, 1.74533e-23, 5.8165e-23, 1.93741e-22, 6.44234e-22, 2.13996e-21, 7.09849e-21, 2.35198e-20, 7.78163e-20, 2.57034e-19, 8.47696e-19, 2.791e-18, 9.17304e-18, 3.00921e-17, 9.85202e-17, 3.21868e-16, 1.04915e-15, 3.41149e-15, 1.10641e-14, 3.57819e-14, 1.1537e-13, 3.70764e-13, 1.18727e-12, 3.78714e-12, 1.20288e-11, 3.80274e-11, 1.19598e-10, 3.73986e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96773e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, 0, },
510 { 3.36833e-28, 1.08302e-27, 3.65738e-27, 1.23253e-26, 4.1328e-26, 1.38446e-25, 4.63721e-25, 1.5651e-24, 5.2312e-24, 1.74518e-23, 5.81941e-23, 1.93732e-22, 6.44301e-22, 2.13994e-21, 7.09989e-21, 2.35234e-20, 7.78156e-20, 2.57039e-19, 8.477e-19, 2.79101e-18, 9.17308e-18, 3.00922e-17, 9.8521e-17, 3.21868e-16, 1.04915e-15, 3.41149e-15, 1.10641e-14, 3.5782e-14, 1.1537e-13, 3.70764e-13, 1.18727e-12, 3.78714e-12, 1.20288e-11, 3.80274e-11, 1.19598e-10, 3.73986e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96773e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, 0, },
511 { 1.00177e-28, 3.22356e-28, 1.08943e-27, 3.67422e-27, 1.23301e-26, 4.134e-26, 1.38589e-25, 4.68182e-25, 1.56635e-24, 5.23077e-24, 1.74606e-23, 5.81915e-23, 1.93753e-22, 6.44297e-22, 2.14036e-21, 7.10097e-21, 2.35232e-20, 7.78171e-20, 2.5704e-19, 8.47702e-19, 2.79102e-18, 9.17309e-18, 3.00924e-17, 9.85211e-17, 3.21868e-16, 1.04916e-15, 3.41149e-15, 1.10641e-14, 3.57819e-14, 1.1537e-13, 3.70764e-13, 1.18727e-12, 3.78714e-12, 1.20288e-11, 3.80274e-11, 1.19598e-10, 3.73986e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96773e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, 0, },
512 { 2.97726e-29, 9.58791e-29, 3.24267e-28, 1.09445e-27, 3.67566e-27, 1.23337e-26, 4.13827e-26, 1.39922e-25, 4.68554e-25, 1.56622e-24, 5.23339e-24, 1.74598e-23, 5.81976e-23, 1.93751e-22, 6.44424e-22, 2.14069e-21, 7.1009e-21, 2.35237e-20, 7.78174e-20, 2.57041e-19, 8.47705e-19, 2.79102e-18, 9.17316e-18, 3.00925e-17, 9.85212e-17, 3.21868e-16, 1.04916e-15, 3.4115e-15, 1.10641e-14, 3.5782e-14, 1.1537e-13, 3.70764e-13, 1.18727e-12, 3.78714e-12, 1.20288e-11, 3.80274e-11, 1.19598e-10, 3.73986e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96773e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, 0, },
513 { 8.84255e-30, 2.84977e-29, 9.64478e-29, 3.2576e-28, 1.09488e-27, 3.67673e-27, 1.23464e-26, 4.17806e-26, 1.40033e-25, 4.68516e-25, 1.56701e-24, 5.23316e-24, 1.74616e-23, 5.81972e-23, 1.93789e-22, 6.44522e-22, 2.14067e-21, 7.10105e-21, 2.35238e-20, 7.78176e-20, 2.57042e-19, 8.47706e-19, 2.79105e-18, 9.17317e-18, 3.00925e-17, 9.85212e-17, 3.21868e-16, 1.04916e-15, 3.4115e-15, 1.10641e-14, 3.5782e-14, 1.1537e-13, 3.70764e-13, 1.18727e-12, 3.78714e-12, 1.20288e-11, 3.80274e-11, 1.19598e-10, 3.73986e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96773e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, 0, },
514 { 2.62457e-30, 8.46455e-30, 2.86668e-29, 9.68918e-29, 3.25889e-28, 1.09519e-27, 3.68053e-27, 1.24651e-26, 4.18138e-26, 1.40022e-25, 4.68751e-25, 1.56694e-24, 5.2337e-24, 1.74615e-23, 5.82086e-23, 1.93819e-22, 6.44516e-22, 2.14071e-21, 7.10107e-21, 2.35238e-20, 7.78179e-20, 2.57042e-19, 8.47713e-19, 2.79105e-18, 9.17317e-18, 3.00925e-17, 9.85212e-17, 3.21868e-16, 1.04916e-15, 3.4115e-15, 1.10641e-14, 3.5782e-14, 1.15371e-13, 3.70764e-13, 1.18727e-12, 3.78714e-12, 1.20288e-11, 3.80274e-11, 1.19598e-10, 3.73986e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96773e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, 0, },
515 { 7.7852e-31, 2.51256e-30, 8.51483e-30, 2.87988e-29, 9.69303e-29, 3.25983e-28, 1.09633e-27, 3.71591e-27, 1.2475e-26, 4.18105e-26, 1.40092e-25, 4.6873e-25, 1.5671e-24, 5.23367e-24, 1.7465e-23, 5.82175e-23, 1.93817e-22, 6.44529e-22, 2.14072e-21, 7.10109e-21, 2.35239e-20, 7.7818e-20, 2.57044e-19, 8.47714e-19, 2.79105e-18, 9.17318e-18, 3.00925e-17, 9.85212e-17, 3.21868e-16, 1.04916e-15, 3.4115e-15, 1.10641e-14, 3.5782e-14, 1.15371e-13, 3.70764e-13, 1.18727e-12, 3.78714e-12, 1.20288e-11, 3.80274e-11, 1.19598e-10, 3.73986e-10, 1.16198e-09, 3.58432e-09, 1.09664e-08, 3.32406e-08, 9.96773e-08, 2.95169e-07, 8.61198e-07, 2.46835e-06, 6.92244e-06, 1.88924e-05, 4.97818e-05, 0.000125149, 0.000294367, 0.000625222, 0.00110983, 0.00128742, -9.04221e-05, -0.000409004, 0.00029374, 0.00226933, 0.00225775, 0.00382868, 0.0121236, },
516 };
517
518
519}
520
521
522//
523//
524
525
526
527double TrgEclDataBase::GetCMPhi(int tcid)
528{
529
530
531 int index = tcid - 1;
532 if (index == -1 || index > 576) {return 0;}
533
534 const std::vector<double> CMPhi = {20.6916, 10.7535, 9.79258, 22.0817, 23.2945, 46.4763, 35.7108, 34.1959, 46.2791, 47.9495, 70.9428, 59.9706, 58.1014, 69.8504, 71.7511, 93.9289, 83.2808, 81.2993, 92.6644, 94.5637, 115.592, 105.612, 103.734, 114.735, 116.446, 136.23, 127.091, 125.477, 136.179, 137.575, 156.174, 147.928, 146.679, 157.176, 158.181, 175.745, 168.366, 167.535, 177.931, 178.515, 195.249, 188.658, 188.259, 198.664, 198.829, 214.992, 209.059, 209.07, 219.595, 219.379, 235.286, 229.821, 230.186, 240.934, 240.417, 256.461, 251.192, 251.805, 262.868, 262.176, 278.838, 273.388, 274.09, 285.528, 284.845, 302.654, 296.552, 297.133, 308.953, 308.505, 327.911, 320.683, 320.907, 333.044, 333.052, 354.422, 345.563, 345.234, 357.574, 358.198, 8.10196, 8.02207, 7.9584, 7.90825, 7.86997, 7.84221, 7.82458, 7.82007, 7.83209, 7.85461, 7.88748, 7.9192, 18.8526, 18.6718, 18.5285, 18.4173, 18.3338, 18.2755, 18.241, 18.2332, 18.2552, 18.3014, 18.372, 18.4416, 29.5551, 29.2811, 29.0642, 28.8961, 28.7707, 28.6839, 28.6335, 28.6227, 28.6542, 28.7225, 28.8281, 28.9329, 40.1849, 39.8292, 39.5473, 39.329, 39.1663, 39.0542, 38.99, 38.9762, 39.0161, 39.104, 39.2409, 39.3769, 50.7216, 50.2986, 49.9629, 49.7024, 49.5084, 49.3749, 49.2989, 49.2827, 49.3297, 49.4342, 49.5973, 49.7597, 61.1489, 60.675, 60.2979, 60.0049, 59.7865, 59.6363, 59.5509, 59.5328, 59.5856, 59.7028, 59.8865, 60.0694, 71.4555, 70.9479, 70.5429, 70.2276, 69.9924, 69.8306, 69.7386, 69.7194, 69.7759, 69.9024, 70.1003, 70.2971, 81.6352, 81.1113, 80.6921, 80.3651, 80.1206, 79.9524, 79.857, 79.837, 79.8956, 80.0271, 80.2327, 80.4372, 91.6863, 91.1628, 90.7428, 90.4143, 90.1684, 89.999, 89.9029, 89.8827, 89.9418, 90.0743, 90.2813, 90.4868, 101.611, 101.104, 100.695, 100.375, 100.135, 99.9698, 99.8758, 99.856, 99.9139, 100.043, 100.246, 100.446, 111.416, 110.939, 110.554, 110.251, 110.024, 109.867, 109.777, 109.759, 109.814, 109.937, 110.128, 110.318, 121.11, 120.675, 120.323, 120.046, 119.838, 119.694, 119.611, 119.594, 119.645, 119.758, 119.934, 120.108, 130.704, 130.322, 130.012, 129.768, 129.584, 129.456, 129.383, 129.367, 129.413, 129.513, 129.669, 129.822, 140.211, 139.891, 139.631, 139.426, 139.27, 139.162, 139.1, 139.086, 139.125, 139.21, 139.342, 139.471, 149.648, 149.397, 149.192, 149.03, 148.907, 148.821, 148.77, 148.759, 148.791, 148.859, 148.963, 149.066, 159.029, 158.851, 158.706, 158.591, 158.503, 158.441, 158.405, 158.396, 158.42, 158.469, 158.544, 158.617, 168.371, 168.271, 168.189, 168.123, 168.072, 168.035, 168.013, 168.008, 168.022, 168.052, 168.095, 168.137, 177.692, 177.671, 177.653, 177.638, 177.625, 177.615, 177.607, 177.604, 177.61, 177.619, 177.631, 177.641, 187.008, 187.066, 187.113, 187.149, 187.175, 187.191, 187.198, 187.198, 187.195, 187.184, 187.163, 187.141, 196.337, 196.474, 196.585, 196.671, 196.734, 196.776, 196.797, 196.801, 196.789, 196.758, 196.705, 196.652, 205.697, 205.91, 206.081, 206.216, 206.316, 206.382, 206.417, 206.423, 206.404, 206.353, 206.27, 206.186, 215.104, 215.388, 215.617, 215.797, 215.931, 216.021, 216.068, 216.077, 216.05, 215.981, 215.869, 215.757, 224.574, 224.924, 225.206, 225.427, 225.591, 225.702, 225.761, 225.772, 225.738, 225.653, 225.515, 225.378, 234.124, 234.531, 234.859, 235.116, 235.306, 235.435, 235.504, 235.517, 235.477, 235.378, 235.218, 235.058, 243.767, 244.222, 244.588, 244.873, 245.086, 245.229, 245.306, 245.321, 245.276, 245.166, 244.988, 244.81, 253.515, 254.006, 254.401, 254.708, 254.936, 255.09, 255.174, 255.189, 255.141, 255.022, 254.831, 254.64, 263.379, 263.894, 264.306, 264.626, 264.864, 265.024, 265.111, 265.127, 265.076, 264.953, 264.755, 264.555, 273.367, 273.89, 274.308, 274.632, 274.871, 275.033, 275.12, 275.137, 275.086, 274.961, 274.761, 274.56, 283.483, 283.998, 284.408, 284.725, 284.96, 285.117, 285.202, 285.219, 285.169, 285.048, 284.852, 284.655, 293.728, 294.217, 294.606, 294.906, 295.128, 295.276, 295.356, 295.371, 295.325, 295.211, 295.026, 294.84, 304.096, 304.543, 304.897, 305.17, 305.37, 305.505, 305.577, 305.59, 305.548, 305.446, 305.279, 305.11, 314.58, 314.968, 315.274, 315.51, 315.682, 315.797, 315.858, 315.87, 315.834, 315.746, 315.603, 315.458, 325.166, 325.48, 325.726, 325.915, 326.052, 326.144, 326.192, 326.201, 326.174, 326.104, 325.989, 325.873, 335.836, 336.061, 336.237, 336.372, 336.47, 336.535, 336.568, 336.574, 336.556, 336.507, 336.426, 336.342, 346.566, 346.693, 346.793, 346.868, 346.923, 346.958, 346.975, 346.977, 346.968, 346.942, 346.898, 346.852, 357.332, 357.356, 357.374, 357.386, 357.394, 357.397, 357.397, 357.398, 357.399, 357.397, 357.391, 357.344, 22.8924, 21.4447, 9.50308, 10.7985, 46.8294, 45.0749, 33.2578, 34.918, 70.2369, 68.3176, 56.7095, 58.6126, 92.9993, 91.0699, 79.7134, 81.704, 115.12, 113.316, 102.211, 104.139, 136.694, 135.116, 124.225, 125.97, 157.875, 156.585, 145.846, 147.326, 178.846, 177.871, 167.208, 168.377, 199.806, 199.141, 188.469, 189.318, 220.954, 220.56, 209.797, 210.348, 242.477, 242.283, 231.354, 231.661, 264.531, 264.433, 253.28, 253.433, 287.221, 287.083, 275.674, 275.797, 310.561, 310.233, 298.577, 298.816, 334.449, 333.798, 321.947, 322.451, 358.254, 357.566, 345.656, 346.536};
535 // cppcheck-suppress containerOutOfBounds
536 return CMPhi[index];
537
538
539}
540
541
542double TrgEclDataBase::GetCMTheta(int tcid)
543{
544
545 int index = tcid - 1;
546 if (index == -1 || index >= 576) {return 0;}
547 const double CMTheta[576] = {17.7315, 26.6566, 35.2717, 35.2375, 26.8222, 18.3848, 27.1, 35.6647, 35.8264, 27.4757, 19.3899, 27.9295, 36.4172, 36.7099, 28.4402, 20.5497, 28.9843, 37.3876, 37.7295, 29.5395, 21.6747, 30.0836, 38.4111, 38.7217, 30.5977, 22.6118, 31.0639, 39.3335, 39.5456, 31.4681, 23.2526, 31.7978, 40.0311, 40.0973, 32.0444, 23.5314, 32.2004, 40.4201, 40.3129, 32.2625, 23.4214, 32.2285, 40.4566, 40.1687, 32.0997, 22.933, 31.879, 40.137, 39.6803, 31.5729, 22.1146, 31.1893, 39.4964, 38.9036, 30.74, 21.0555, 30.2382, 38.6109, 37.9349, 29.701, 19.8894, 29.1468, 37.5961, 36.9083, 28.5972, 18.7918, 28.0743, 36.6016, 35.9841, 27.6015, 17.9622, 27.2007, 35.7921, 35.3251, 26.8926, 18.1528, 26.6928, 35.3174, 35.6425, 27.3093, 42.6511, 50.1664, 58.5442, 67.7097, 77.5183, 87.754, 98.1393, 108.338, 117.941, 126.69, 134.503, 139.718, 42.7522, 50.257, 58.6227, 67.7736, 77.566, 87.7843, 98.1521, 108.334, 117.922, 126.658, 134.461, 139.669, 42.93, 50.4166, 58.7605, 67.8859, 77.6499, 87.8377, 98.1745, 108.327, 117.889, 126.602, 134.387, 139.585, 43.1774, 50.6388, 58.9526, 68.0427, 77.7669, 87.9122, 98.2057, 108.317, 117.842, 126.525, 134.284, 139.467, 43.4846, 50.915, 59.1913, 68.2376, 77.9126, 88.0048, 98.2447, 108.305, 117.784, 126.428, 134.156, 139.32, 43.8404, 51.2351, 59.4682, 68.4639, 78.0816, 88.1125, 98.2899, 108.29, 117.716, 126.315, 134.007, 139.149, 44.2317, 51.5875, 59.7733, 68.7132, 78.2681, 88.2313, 98.3398, 108.274, 117.641, 126.191, 133.843, 138.96, 44.6454, 51.9604, 60.0963, 68.9775, 78.4655, 88.3571, 98.3928, 108.257, 117.562, 126.059, 133.668, 138.76, 45.0678, 52.3415, 60.4268, 69.2477, 78.6679, 88.4861, 98.447, 108.24, 117.481, 125.924, 133.489, 138.554, 45.486, 52.7194, 60.7543, 69.516, 78.8686, 88.6141, 98.5008, 108.223, 117.4, 125.79, 133.312, 138.35, 45.8873, 53.0822, 61.0694, 69.7739, 79.0618, 88.7373, 98.5527, 108.206, 117.322, 125.66, 133.14, 138.153, 46.2609, 53.4201, 61.363, 70.0144, 79.2419, 88.8523, 98.6011, 108.191, 117.25, 125.539, 132.98, 137.97, 46.5968, 53.7244, 61.6274, 70.2311, 79.4043, 88.9559, 98.6448, 108.177, 117.184, 125.43, 132.835, 137.804, 46.8865, 53.9867, 61.8555, 70.4182, 79.5445, 89.0454, 98.6825, 108.165, 117.127, 125.336, 132.711, 137.661, 47.1231, 54.2012, 62.042, 70.5712, 79.6591, 89.1186, 98.7133, 108.155, 117.081, 125.258, 132.608, 137.543, 47.3007, 54.3622, 62.1821, 70.6861, 79.7453, 89.1737, 98.7365, 108.147, 117.046, 125.2, 132.531, 137.455, 47.4156, 54.4665, 62.2727, 70.7606, 79.801, 89.2093, 98.7516, 108.142, 117.023, 125.163, 132.482, 137.398, 47.4653, 54.5114, 62.312, 70.7926, 79.8253, 89.2247, 98.7581, 108.14, 117.014, 125.146, 132.46, 137.373, 47.4486, 54.4964, 62.2987, 70.7819, 79.8171, 89.2196, 98.7559, 108.141, 117.017, 125.152, 132.467, 137.382, 47.3656, 54.4212, 62.2334, 70.7282, 79.7769, 89.1939, 98.7451, 108.145, 117.033, 125.179, 132.503, 137.423, 47.2186, 54.2878, 62.1174, 70.633, 79.7056, 89.1483, 98.7259, 108.151, 117.062, 125.227, 132.567, 137.496, 47.0107, 54.0994, 61.9536, 70.4987, 79.6049, 89.084, 98.6988, 108.159, 117.103, 125.295, 132.657, 137.599, 46.7467, 53.8601, 61.7455, 70.3281, 79.477, 89.0024, 98.6644, 108.17, 117.154, 125.381, 132.771, 137.73, 46.4328, 53.5759, 61.4984, 70.1256, 79.3253, 88.9055, 98.6236, 108.183, 117.216, 125.483, 132.906, 137.885, 46.0766, 53.2535, 61.2183, 69.8961, 79.1534, 88.7959, 98.5774, 108.198, 117.285, 125.599, 133.059, 138.06, 45.6877, 52.9017, 60.9128, 69.646, 78.9659, 88.6763, 98.5271, 108.214, 117.361, 125.724, 133.225, 138.251, 45.2763, 52.5299, 60.5904, 69.3817, 78.7683, 88.5502, 98.474, 108.231, 117.44, 125.857, 133.4, 138.452, 44.8546, 52.1493, 60.26, 69.1115, 78.566, 88.4213, 98.4198, 108.249, 117.522, 125.992, 133.58, 138.658, 44.4348, 51.7706, 59.932, 68.8431, 78.3653, 88.2933, 98.366, 108.266, 117.602, 126.126, 133.757, 138.861, 44.0306, 51.4064, 59.6166, 68.5852, 78.1724, 88.1705, 98.3143, 108.282, 117.68, 126.255, 133.927, 139.057, 43.6556, 51.0689, 59.3246, 68.3466, 77.9941, 88.0568, 98.2665, 108.298, 117.751, 126.374, 134.084, 139.237, 43.3228, 50.7696, 59.0656, 68.1352, 77.8361, 87.9562, 98.2243, 108.311, 117.814, 126.479, 134.223, 139.397, 43.0446, 50.5196, 58.8495, 67.9587, 77.7043, 87.8723, 98.189, 108.322, 117.867, 126.566, 134.339, 139.53, 42.8312, 50.3279, 58.684, 67.8236, 77.6035, 87.8082, 98.1621, 108.331, 117.907, 126.633, 134.428, 139.632, 42.6913, 50.2024, 58.5755, 67.7352, 77.5373, 87.7661, 98.1444, 108.337, 117.934, 126.677, 134.486, 139.699, 43.068, 50.6411, 59.0737, 68.2876, 78.1345, 88.3931, 98.7829, 108.961, 118.52, 127.215, 134.972, 140.314, 155.334, 145.937, 145.824, 155.423, 154.984, 145.644, 145.63, 155.185, 154.464, 145.2, 145.254, 154.739, 153.865, 144.682, 144.764, 154.168, 153.283, 144.172, 144.242, 153.566, 152.799, 143.744, 143.767, 153.023, 152.478, 143.456, 143.404, 152.615, 152.359, 143.343, 143.2, 152.391, 152.455, 143.419, 143.181, 152.38, 152.755, 143.674, 143.348, 152.581, 153.223, 144.077, 143.682, 152.971, 153.799, 144.576, 144.139, 153.502, 154.4, 145.099, 144.658, 154.102, 154.933, 145.565, 145.161, 154.681, 155.306, 145.893, 145.566, 155.145, 156.461, 146.543, 145.802, 155.409, };
548
549 return CMTheta[index];
550
551
552}
553
554
555//
556//
557//
558double TrgEclDataBase::GetCMEnergy(int tcid)
559{
560
561 int index = tcid - 1;
562 if (index == -1 || index > 576) {return 0;}
563
564 const std::vector<double> CM1GeV = {1.32241, 1.31221, 1.2956, 1.29427, 1.31074, 1.3186, 1.30827, 1.29087, 1.28717, 1.30489, 1.31247, 1.30075, 1.28168, 1.27634, 1.29601, 1.30501, 1.29088, 1.26958, 1.26355, 1.28555, 1.29738, 1.28024, 1.25653, 1.25082, 1.27514, 1.29072, 1.27045, 1.24451, 1.24004, 1.26633, 1.28602, 1.26293, 1.23526, 1.23272, 1.26038, 1.28393, 1.25875, 1.23004, 1.22983, 1.2581, 1.28476, 1.25845, 1.22955, 1.23176, 1.2598, 1.28838, 1.26209, 1.23385, 1.23826, 1.26525, 1.29429, 1.26918, 1.24236, 1.24846, 1.27372, 1.30163, 1.27871, 1.25394, 1.26094, 1.28398, 1.30931, 1.28933, 1.26694, 1.27387, 1.29454, 1.31615, 1.29941, 1.2794, 1.28526, 1.30375, 1.32108, 1.30737, 1.28932, 1.29322, 1.31012, 1.32282, 1.31189, 1.29505, 1.29508, 1.31148, 1.27686, 1.25364, 1.22346, 1.18599, 1.1418, 1.09253, 1.04086, 0.99008, 0.943691, 0.903755, 0.870748, 0.850432, 1.27534, 1.25197, 1.22166, 1.18411, 1.1399, 1.09068, 1.03911, 0.988488, 0.942277, 0.902522, 0.869687, 0.849489, 1.27267, 1.24903, 1.2185, 1.18081, 1.13655, 1.08742, 1.03604, 0.985685, 0.939787, 0.90035, 0.867818, 0.847828, 1.26894, 1.24492, 1.21409, 1.1762, 1.1319, 1.08287, 1.03175, 0.981777, 0.936313, 0.89732, 0.865208, 0.845506, 1.26428, 1.2398, 1.20858, 1.17045, 1.12609, 1.07721, 1.02642, 0.976906, 0.931981, 0.893539, 0.861949, 0.842606, 1.25885, 1.23383, 1.20218, 1.16378, 1.11935, 1.07063, 1.02022, 0.971248, 0.926946, 0.889139, 0.858154, 0.839225, 1.25284, 1.22723, 1.19511, 1.15641, 1.11192, 1.06338, 1.01338, 0.965, 0.92138, 0.884272, 0.85395, 0.835478, 1.24644, 1.22021, 1.18759, 1.14859, 1.10403, 1.05568, 1.00613, 0.95837, 0.91547, 0.879097, 0.849476, 0.831486, 1.23985, 1.213, 1.17988, 1.14057, 1.09595, 1.04781, 0.998696, 0.951575, 0.909406, 0.873782, 0.844875, 0.827378, 1.23328, 1.20582, 1.17221, 1.1326, 1.08792, 1.03998, 0.991319, 0.944824, 0.903377, 0.868491, 0.840288, 0.823278, 1.22694, 1.19889, 1.16482, 1.12492, 1.08019, 1.03245, 0.984214, 0.938318, 0.897562, 0.863382, 0.835854, 0.819311, 1.22099, 1.1924, 1.1579, 1.11775, 1.07298, 1.02543, 0.977585, 0.932244, 0.892127, 0.858603, 0.831702, 0.815592, 1.21561, 1.18654, 1.15166, 1.11128, 1.06648, 1.0191, 0.97161, 0.926768, 0.887223, 0.854286, 0.827947, 0.812227, 1.21095, 1.18146, 1.14626, 1.10569, 1.06086, 1.01362, 0.966447, 0.922033, 0.882981, 0.850548, 0.824692, 0.809309, 1.20712, 1.1773, 1.14184, 1.10112, 1.05626, 1.00915, 0.962225, 0.918159, 0.879508, 0.847486, 0.822025, 0.806915, 1.20424, 1.17417, 1.13851, 1.09768, 1.0528, 1.00579, 0.959047, 0.915243, 0.876893, 0.845179, 0.820013, 0.805108, 1.20237, 1.17213, 1.13635, 1.09544, 1.05057, 1.00361, 0.956989, 0.913354, 0.875198, 0.843683, 0.818709, 0.803937, 1.20156, 1.17126, 1.13542, 1.09448, 1.0496, 1.00266, 0.956099, 0.912537, 0.874465, 0.843036, 0.818144, 0.80343, 1.20184, 1.17155, 1.13573, 1.09481, 1.04992, 1.00298, 0.956397, 0.91281, 0.874711, 0.843254, 0.818334, 0.803601, 1.20318, 1.17302, 1.13729, 1.09641, 1.05154, 1.00455, 0.957876, 0.914168, 0.87593, 0.84433, 0.819274, 0.804445, 1.20557, 1.17561, 1.14004, 1.09926, 1.0544, 1.00733, 0.960503, 0.916578, 0.878092, 0.846239, 0.820939, 0.805941, 1.20894, 1.17927, 1.14393, 1.10329, 1.05844, 1.01126, 0.964215, 0.919984, 0.881147, 0.848934, 0.823288, 0.80805, 1.2132, 1.18391, 1.14886, 1.10839, 1.06356, 1.01625, 0.968923, 0.924302, 0.885018, 0.852346, 0.826261, 0.810717, 1.21824, 1.1894, 1.15471, 1.11444, 1.06964, 1.02217, 0.97451, 0.929424, 0.889607, 0.856388, 0.829779, 0.813871, 1.22393, 1.1956, 1.16131, 1.12128, 1.07653, 1.02888, 0.980834, 0.93522, 0.894795, 0.860953, 0.833747, 0.817426, 1.2301, 1.20234, 1.16849, 1.12874, 1.08402, 1.03618, 0.987726, 0.941532, 0.90044, 0.865915, 0.838057, 0.821283, 1.23658, 1.20942, 1.17606, 1.13659, 1.09193, 1.04389, 0.994993, 0.948184, 0.906384, 0.871135, 0.842584, 0.825332, 1.24318, 1.21665, 1.18378, 1.14461, 1.10002, 1.05177, 1.00243, 0.954984, 0.912454, 0.876459, 0.847195, 0.829452, 1.2497, 1.22379, 1.19142, 1.15257, 1.10803, 1.05958, 1.0098, 0.961723, 0.918464, 0.881724, 0.851751, 0.833518, 1.25593, 1.23063, 1.19874, 1.16019, 1.11573, 1.06709, 1.01687, 0.968188, 0.924225, 0.886765, 0.856106, 0.837403, 1.26167, 1.23693, 1.2055, 1.16724, 1.12284, 1.07403, 1.02342, 0.974167, 0.929548, 0.891417, 0.860122, 0.84098, 1.26674, 1.2425, 1.21148, 1.17347, 1.12914, 1.08018, 1.02921, 0.979454, 0.934251, 0.895524, 0.863663, 0.844132, 1.27095, 1.24713, 1.21646, 1.17867, 1.13439, 1.0853, 1.03404, 0.983863, 0.93817, 0.898943, 0.866608, 0.846753, 1.27416, 1.25067, 1.22026, 1.18264, 1.13841, 1.08922, 1.03774, 0.987234, 0.941165, 0.901554, 0.868855, 0.84875, 1.27626, 1.25298, 1.22275, 1.18524, 1.14104, 1.09179, 1.04016, 0.989441, 0.943125, 0.903262, 0.870324, 0.850056, 1.27593, 1.25233, 1.22172, 1.18382, 1.13925, 1.08973, 1.03798, 0.987338, 0.94127, 0.901723, 0.869107, 0.84843, 0.798408, 0.826887, 0.828557, 0.799225, 0.795186, 0.822604, 0.825698, 0.797044, 0.790324, 0.816059, 0.820141, 0.792916, 0.784607, 0.808329, 0.812819, 0.787518, 0.778923, 0.800626, 0.804911, 0.781705, 0.774121, 0.794093, 0.797623, 0.776364, 0.77089, 0.789651, 0.792008, 0.772275, 0.769678, 0.7879, 0.788839, 0.770015, 0.770651, 0.789071, 0.788539, 0.769894, 0.773678, 0.793011, 0.791147, 0.771929, 0.778336, 0.79919, 0.796317, 0.775841, 0.783962, 0.806745, 0.803343, 0.78108, 0.78972, 0.814567, 0.81122, 0.786882, 0.794718, 0.821443, 0.81876, 0.79237, 0.798153, 0.826252, 0.824762, 0.796675, 0.796741, 0.826444, 0.828226, 0.799095, };
565 // cppcheck-suppress containerOutOfBounds
566 return CM1GeV[index];
567
568
569}
570
571int TrgEclDataBase::Get3DBhabhaLUT(int tcid)
572{
573
574 const std::vector<int> BhabhaLUT = {0, 147805, 221373, 286893, 287085, 221565, 148205, 221773,
575 295469, 295661, 221965, 156797, 230349, 295853, 304237, 230541,
576 173549, 238909, 304413, 312797, 247293, 182093, 247469, 312973,
577 321341, 255821, 190605, 255997, 321500, 329868, 256173, 190925,
578 264525, 330044, 330204, 264685, 199437, 264845, 330380, 330540,
579 265021, 191549, 265181, 330700, 330876, 265341, 191869, 265501,
580 331036, 331212, 265661, 183997, 257645, 323180, 323356, 257805,
581 176141, 249789, 323533, 315517, 249965, 168317, 241949, 315693,
582 307693, 242141, 160509, 234141, 307869, 299869, 234333, 152717,
583 226333, 300061, 292061, 226525, 153133, 226733, 292253, 300653,
584 226925, 352397, 409741, 483468, 557196, 639115, 721035, 802954,
585 884874, 966793, 1040521, 1106057, 1147017, 352573, 409917, 483644,
586 557356, 639275, 721195, 803114, 885034, 966953, 1040681, 1098025,
587 1147176, 352749, 410076, 483804, 557532, 639451, 721371, 803290,
588 885210, 967129, 1040857, 1098201, 1147352, 352909, 418444, 483980,
589 557692, 639611, 721531, 803450, 885370, 967289, 1041017, 1098361,
590 1139320, 353085, 418604, 484140, 557868, 639787, 721691, 803610,
591 885530, 967449, 1032985, 1098537, 1139496, 361437, 418780, 484300,
592 558028, 639947, 721867, 803786, 885706, 967625, 1033161, 1098697,
593 1139656, 361597, 427132, 492668, 566380, 640107, 722027, 803946,
594 885866, 967785, 1033321, 1098857, 1139816, 369964, 427292, 492828,
595 566539, 640267, 722187, 804106, 886026, 967945, 1033481, 1099016,
596 1139976, 370124, 427452, 492988, 566699, 648619, 722346, 804266,
597 886186, 959913, 1033641, 1090984, 1140136, 370284, 435804, 501340,
598 575051, 648779, 730698, 812618, 886345, 960073, 1033801, 1091144,
599 1132104, 378620, 435964, 501500, 575211, 648939, 730858, 812778,
600 886505, 960233, 1033961, 1091304, 1132264, 378780, 436124, 501644,
601 575371, 649099, 731018, 812938, 886665, 960393, 1034121, 1091464,
602 1132424, 387132, 444460, 509996, 575531, 649259, 731162, 813082,
603 886809, 960537, 1026089, 1091624, 1132584, 387276, 444620, 510155,
604 575675, 657595, 731322, 813242, 886969, 960697, 1026233, 1091768,
605 1132728, 387436, 444764, 510299, 584027, 657755, 731482, 813402,
606 887129, 960857, 1026392, 1091928, 1132888, 387580, 444924, 510459,
607 584187, 657915, 731626, 813546, 887273, 961001, 1026536, 1092088,
608 1124856, 387724, 445068, 510603, 584331, 658059, 731786, 813706,
609 887433, 961161, 1026696, 1084040, 1125000, 387884, 453420, 510763,
610 584491, 658218, 731946, 813866, 887593, 961321, 1026856, 1084200,
611 1125160, 388028, 445372, 510907, 584635, 658362, 732090, 814010,
612 887737, 961465, 1027000, 1084344, 1125304, 388172, 445516, 511067,
613 584795, 658523, 732250, 814170, 887897, 961625, 1027160, 1092696,
614 1125464, 388332, 445676, 511211, 584939, 658667, 732394, 814314,
615 888041, 961769, 1027304, 1092840, 1125608, 388476, 445820, 511371,
616 576907, 658827, 732554, 814474, 888201, 961929, 1027464, 1093000,
617 1133960, 388636, 445980, 511515, 577051, 650795, 732714, 814634,
618 888361, 962089, 1027625, 1093160, 1134104, 380588, 446140, 503484,
619 577211, 650939, 732858, 814794, 888521, 962233, 1027769, 1093304,
620 1134264, 380748, 438092, 503644, 577371, 651099, 733018, 814938,
621 888665, 962393, 1036121, 1093464, 1134424, 380908, 438252, 503788,
622 577531, 651259, 733178, 815098, 888825, 962553, 1036281, 1093624,
623 1134584, 372860, 438412, 503948, 569499, 651419, 733338, 807066,
624 888985, 962713, 1036441, 1093784, 1134744, 373020, 430380, 495916,
625 569659, 651579, 725307, 807226, 889146, 971065, 1036601, 1102136,
626 1143096, 364988, 430540, 496076, 569820, 643547, 725467, 807386,
627 889306, 971225, 1036761, 1102297, 1143256, 365165, 422508, 496252,
628 569980, 643707, 725627, 807546, 889466, 971385, 1036921, 1102457,
629 1143416, 365325, 422684, 488220, 561948, 643867, 725803, 807722,
630 889642, 971561, 1037081, 1102617, 1143576, 357309, 422844, 488380,
631 562124, 644043, 725963, 807882, 889802, 971721, 1037257, 1102793,
632 1143736, 357469, 423004, 488556, 562284, 644203, 726123, 808042,
633 889962, 971881, 1045609, 1102953, 1152104, 357645, 414989, 488716,
634 562444, 644363, 726299, 808218, 890138, 972057, 1045785, 1103113,
635 1152264, 357821, 415165, 488892, 562620, 644539, 726459, 808378,
636 890298, 972217, 1045945, 1103289, 1152441, 357981, 423517, 489052,
637 562780, 644699, 726619, 816730, 898650, 980569, 1046105, 1111641,
638 1152600, 1270136, 1196376, 1196200, 1269944, 1270520, 1196760, 1196568,
639 1270328, 1262696, 1188936, 1188760, 1270712, 1263064, 1189304, 1189128,
640 1262888, 1255224, 1181464, 1181288, 1263240, 1255576, 1181816, 1181640,
641 1255400, 1247720, 1173976, 1173800, 1255736, 1248056, 1174312, 1174136,
642 1247880, 1248392, 1174648, 1174472, 1248216, 1256920, 1183192, 1174824,
643 1256744, 1257256, 1183528, 1183352, 1257096, 1265816, 1192072, 1183704,
644 1265624, 1266168, 1192440, 1192264, 1265992, 1274744, 1201000, 1192632,
645 1274552, 1275112, 1201384, 1201192, 1274920, 1283688, 1209960, 1201576,
646 1275320, 0, 0, 0, 0, 0, 0, 0,
647 0, 0, 0, 0, 0, 0, 0, 0,
648 0, 0, 0, 0, 0, 0, 0, 0,
649 0, 0, 0, 0, 0, 0, 0, 0,
650 0, 0, 0, 0, 0, 0, 0, 0,
651 0, 0, 0, 0, 0, 0, 0, 0,
652 0, 0, 0, 0, 0, 0, 0, 0,
653 0, 0, 0, 0, 0, 0, 0, 0,
654 0, 0, 0, 0, 0, 0, 0, 0,
655 0, 0, 0, 0, 0, 0, 0, 0,
656 0, 0, 0, 0, 0, 0, 0, 0,
657 0, 0, 0, 0, 0, 0, 0, 0,
658 0, 0, 0, 0, 0, 0, 0, 0,
659 0, 0, 0, 0, 0, 0, 0, 0,
660 0, 0, 0, 0, 0, 0, 0, 0,
661 0, 0, 0, 0, 0, 0, 0, 0,
662 0, 0, 0, 0, 0, 0, 0, 0,
663 0, 0, 0, 0, 0, 0, 0, 0,
664 0, 0, 0, 0, 0, 0, 0, 0,
665 0, 0, 0, 0, 0, 0, 0, 0,
666 0, 0, 0, 0, 0, 0, 0, 0,
667 0, 0, 0, 0, 0, 0, 0, 0,
668 0, 0, 0, 0, 0, 0, 0, 0,
669 0, 0, 0, 0, 0, 0, 0, 0,
670 0, 0, 0, 0, 0, 0, 0, 0,
671 0, 0, 0, 0, 0, 0, 0, 0,
672 0, 0, 0, 0, 0, 0, 0, 0,
673 0, 0, 0, 0, 0, 0, 0, 0,
674 0, 0, 0, 0, 0, 0, 0, 0,
675 0, 0, 0, 0, 0, 0, 0, 0,
676 0, 0, 0, 0, 0, 0, 0, 0,
677 0, 0, 0, 0, 0, 0, 0, 0,
678 0, 0, 0, 0, 0, 0, 0, 0,
679 0, 0, 0, 0, 0, 0, 0, 0,
680 0, 0, 0, 0, 0, 0, 0, 0,
681 0, 0, 0, 0, 0, 0, 0, 0,
682 0, 0, 0, 0, 0, 0, 0, 0,
683 0, 0, 0, 0, 0, 0, 0, 0,
684 0, 0, 0, 0, 0, 0, 0, 0,
685 0, 0, 0, 0, 0, 0, 0, 0,
686 0, 0, 0, 0, 0, 0, 0, 0,
687 0, 0, 0, 0, 0, 0, 0, 0,
688 0, 0, 0, 0, 0, 0, 0, 0,
689 0, 0, 0, 0, 0, 0, 0, 0,
690 0, 0, 0, 0, 0, 0, 0, 0,
691 0, 0, 0, 0, 0, 0, 0, 0,
692 0, 0, 0, 0, 0, 0, 0, 0,
693 0, 0, 0, 0, 0, 0, 0, 0,
694 0, 0, 0, 0, 0, 0, 0, 0,
695 0, 0, 0, 0, 0, 0, 0, 0,
696 0, 0, 0, 0, 0, 0, 0, 0,
697 0, 0, 0, 0, 0, 0, 0, 0,
698 0, 0, 0, 0, 0, 0, 0, 0,
699 0, 0, 0, 0, 0, 0, 0, 0,
700 0, 0, 0, 0, 0, 0, 0, 0,
701 0, 0, 0, 0, 0, 0, 0, 0
702 };
703 return BhabhaLUT[tcid];
704
705}
706
707
708
709//
710//
711
712
713double
714TrgEclDataBase::interFADC(double timing, std::vector<int> SignalAmp)
715{
716
717 // std::vector<double> SignalAmp;
718 std::vector<double> SignalTime;
719
720 // SignalAmp.clear();
721 SignalTime.clear();
722
723 timing = timing * 1000;
724
725 int startbin = (int)(timing) / 125;
726 int endbin = startbin + 1;
727
728 SignalTime = {0, 125, 250, 375, 500, 625, 750, 875};
729
730 double E_out = 0;
731 if (timing < 0 || timing > 875) {return E_out;}
732
733 if (timing > SignalTime[startbin] && timing < SignalTime[endbin]) {
734 E_out = ((((double)SignalAmp[endbin] - (double)SignalAmp[startbin])) / (SignalTime[endbin] - SignalTime[startbin])) *
735 (timing - SignalTime[startbin]) + (double)SignalAmp[startbin];
736 return E_out;
737 } else {
738 E_out = 0;
739 }
740
741
742 return E_out;
743
744}
Belle2::TOP::SignalPDF
Parametrization of signal PDF in a single pixel.
Definition: SignalPDF.h:25
Belle2::TrgEclDataBase::GetCMEnergy
double GetCMEnergy(int)
TC CM Phi
Definition: TrgEclDataBase.cc:558
Belle2::TrgEclDataBase::Get3DBhabhaLUT
int Get3DBhabhaLUT(int)
TC CM Phi
Definition: TrgEclDataBase.cc:571
Belle2::TrgEclDataBase::GetCMPhi
double GetCMPhi(int)
TC CM Phi
Definition: TrgEclDataBase.cc:527
Belle2::TrgEclDataBase::getCoeffNoise
void getCoeffNoise(int, std::vector< std::vector< double > > &, std::vector< std::vector< double > > &, std::vector< std::vector< double > > &)
read coefficient for fit
Definition: TrgEclDataBase.cc:293
Belle2::TrgEclDataBase::readNoiseLMatrix
void readNoiseLMatrix(std::vector< std::vector< double > > &, std::vector< std::vector< double > > &)
Noise Matrix.
Definition: TrgEclDataBase.cc:380
Belle2::TrgEclDataBase::GetTCFLatency
double GetTCFLatency(int)
TC flight time latency
Definition: TrgEclDataBase.cc:224
Belle2::TrgEclDataBase::TrgEclDataBase
TrgEclDataBase()
TrgEclDataBase Constructor.
Definition: TrgEclDataBase.cc:19
Belle2::TrgEclDataBase::getCoeffSigPDF
void getCoeffSigPDF(std::vector< std::vector< double > > &, std::vector< std::vector< double > > &)
read coefficient for fit
Definition: TrgEclDataBase.cc:241
Belle2::TrgEclDataBase::_AmpCoefficient
std::vector< std::vector< int > > _AmpCoefficient
Amplitude Coefficient.
Definition: TrgEclDataBase.h:53
Belle2::TrgEclDataBase::interFADC
double interFADC(double, std::vector< int >)
Util to interpalate Signal Shape.
Definition: TrgEclDataBase.cc:714
Belle2::TrgEclDataBase::_TimingCoefficient
std::vector< std::vector< int > > _TimingCoefficient
Timing Coefficient.
Definition: TrgEclDataBase.h:55
Belle2::TrgEclDataBase::GetCMTheta
double GetCMTheta(int)
TC CM Phi
Definition: TrgEclDataBase.cc:542
Belle2::TrgEclDataBase::MakeFitterCoefficient
void MakeFitterCoefficient(const std::vector< int > &, std::vector< int >)
Make Fitter Coefficients from Signal PDF and Noise covariance metrix.
Definition: TrgEclDataBase.cc:26
Belle2
Abstract base class for different kinds of events.
Definition: MillepedeAlgorithm.h:17
std
STL namespace.