00001 #include "AtNuFluxRW.h"
00002
00003 #include <iostream>
00004 #include <cmath>
00005 #include <cassert>
00006
00007 ClassImp(AtNuFluxRW)
00008
00009 AtNuFluxRW* AtNuFluxRW::fRWInstance = 0;
00010
00011 AtNuFluxRW* AtNuFluxRW::Instance()
00012 {
00013 if( !fRWInstance ){
00014 fRWInstance = new AtNuFluxRW();
00015 }
00016
00017 if( !fRWInstance ) assert(fRWInstance);
00018
00019 return fRWInstance;
00020 }
00021
00022 AtNuFluxRW::AtNuFluxRW()
00023 {
00024 this->LoadFlukaFluxRatios();
00025 this->LoadHondaFluxRatios();
00026 this->LoadBartolFluxRatios();
00027 this->LoadBartolFlavourRatios();
00028 }
00029
00030 AtNuFluxRW::~AtNuFluxRW()
00031 {
00032 if( fHondaFluxRatioNumu ) delete fHondaFluxRatioNumu;
00033 if( fHondaFluxRatioNumuBar ) delete fHondaFluxRatioNumuBar;
00034
00035 if( fFlukaFluxRatioNumu ) delete fFlukaFluxRatioNumu;
00036 if( fFlukaFluxRatioNumuBar ) delete fFlukaFluxRatioNumuBar;
00037
00038 if( fBartolFluxRatioNumu ) delete fBartolFluxRatioNumu;
00039 if( fBartolFluxRatioNumuBar ) delete fBartolFluxRatioNumuBar;
00040
00041 if( fBartolFlavourRatioNumu ) delete fBartolFlavourRatioNumu;
00042 if( fBartolFlavourRatioNumuBar ) delete fBartolFlavourRatioNumuBar;
00043 }
00044
00045 Double_t AtNuFluxRW::LookupRatio(TH2D* hFlux2D, Double_t energy, Double_t coszen)
00046 {
00047 if( hFlux2D==0 ) return 1.0;
00048
00049 Int_t Ebin = hFlux2D->GetXaxis()->FindBin(energy);
00050 if( Ebin<1 ) Ebin = 1;
00051 if( Ebin>hFlux2D->GetXaxis()->GetNbins() ) Ebin = hFlux2D->GetXaxis()->GetNbins();
00052
00053 Int_t Abin = hFlux2D->GetYaxis()->FindBin(coszen);
00054 if( Abin<1 ) Abin = 1;
00055 if( Abin>hFlux2D->GetYaxis()->GetNbins() ) Abin = hFlux2D->GetYaxis()->GetNbins();
00056
00057 return hFlux2D->GetBinContent(Ebin,Abin);
00058 }
00059
00060 Double_t AtNuFluxRW::GetHondaFluxRatio(Int_t idnu, Double_t energy, Double_t coszen)
00061 {
00062
00063 if( idnu==+14 || idnu==+16 ){
00064 return LookupRatio(fHondaFluxRatioNumu,energy,coszen);
00065 }
00066
00067
00068 if( idnu==-14 || idnu==-16 ){
00069 return LookupRatio(fHondaFluxRatioNumuBar,energy,coszen);
00070 }
00071
00072
00073 return 1.0;
00074 }
00075
00076 Double_t AtNuFluxRW::GetFlukaFluxRatio(Int_t idnu, Double_t energy, Double_t coszen)
00077 {
00078
00079 if( idnu==+14 || idnu==+16 ){
00080 return LookupRatio(fFlukaFluxRatioNumu,energy,coszen);
00081 }
00082
00083
00084 if( idnu==-14 || idnu==-16 ){
00085 return LookupRatio(fFlukaFluxRatioNumuBar,energy,coszen);
00086 }
00087
00088
00089 return 1.0;
00090 }
00091
00092 Double_t AtNuFluxRW::GetBartolFluxRatio(Int_t idnu, Double_t energy, Double_t coszen)
00093 {
00094
00095 if( idnu==+14 || idnu==+16 ){
00096 return LookupRatio(fBartolFluxRatioNumu,energy,coszen);
00097 }
00098
00099
00100 if( idnu==-14 || idnu==-16 ){
00101 return LookupRatio(fBartolFluxRatioNumuBar,energy,coszen);
00102 }
00103
00104
00105 return 1.0;
00106 }
00107
00108 Double_t AtNuFluxRW::GetBartolFlavourRatio(Int_t idnu, Double_t energy, Double_t coszen)
00109 {
00110
00111 if( idnu==+14 ){
00112 return LookupRatio(fBartolFlavourRatioNumu,energy,coszen);
00113 }
00114
00115
00116 if( idnu==-14 ){
00117 return LookupRatio(fBartolFlavourRatioNumuBar,energy,coszen);
00118 }
00119
00120
00121 return 1.0;
00122 }
00123
00124 TH2D* AtNuFluxRW::ApplyRatio(TH2D* hFlux2D, TH2D* hRatio2D)
00125 {
00126 TH2D* hCorrectedFlux2D = (TH2D*)(hFlux2D->Clone());
00127 hCorrectedFlux2D->Reset();
00128
00129 TString hName = hFlux2D->GetName();
00130 hName.Append("_Corrected");
00131 hCorrectedFlux2D->SetName(hName.Data());
00132
00133 for( Int_t nEnergy=0; nEnergy<hFlux2D->GetXaxis()->GetNbins(); nEnergy++ ){
00134 for( Int_t nAngle=0; nAngle<hFlux2D->GetYaxis()->GetNbins(); nAngle++ ){
00135 Double_t energy = hFlux2D->GetXaxis()->GetBinCenter(nEnergy+1);
00136 Double_t angle = hFlux2D->GetYaxis()->GetBinCenter(nAngle+1);
00137 Double_t content = hFlux2D->GetBinContent(nEnergy+1,nAngle+1);
00138 Double_t ratio = LookupRatio(hRatio2D,energy,angle);
00139 hCorrectedFlux2D->SetBinContent(nEnergy+1,nAngle+1,content*ratio);
00140 }
00141 }
00142
00143 return hCorrectedFlux2D;
00144 }
00145
00146 TH2D* AtNuFluxRW::ApplyBartolFluxRatios(Int_t idnu, TH2D* hFlux2D)
00147 {
00148 if( idnu==+14 || idnu==+16 ){
00149 return ApplyRatio(hFlux2D,fBartolFluxRatioNumu);
00150 }
00151
00152 if( idnu==-14 || idnu==-16 ){
00153 return ApplyRatio(hFlux2D,fBartolFluxRatioNumuBar);
00154 }
00155
00156 return ApplyRatio(hFlux2D);
00157 }
00158
00159 TH2D* AtNuFluxRW::ApplyBartolFlavourRatios(Int_t idnu, TH2D* hFlux2D)
00160 {
00161 if( idnu==+14 ){
00162 return ApplyRatio(hFlux2D,fBartolFlavourRatioNumu);
00163 }
00164
00165 if( idnu==-14 ){
00166 return ApplyRatio(hFlux2D,fBartolFlavourRatioNumuBar);
00167 }
00168
00169 return ApplyRatio(hFlux2D);
00170 }
00171
00172 void AtNuFluxRW::LoadHondaFluxRatios()
00173 {
00174
00175
00176
00177
00178
00179 Int_t EnergyBins = 50;
00180 Double_t EnergyEdges[51] = {
00181 0.1,
00182 0.125893,
00183 0.158489,
00184 0.199526,
00185 0.251189,
00186 0.316228,
00187 0.398107,
00188 0.501187,
00189 0.630957,
00190 0.794328,
00191 1,
00192 1.25893,
00193 1.58489,
00194 1.99526,
00195 2.51189,
00196 3.16228,
00197 3.98107,
00198 5.01187,
00199 6.30957,
00200 7.94328,
00201 10,
00202 12.5893,
00203 15.8489,
00204 19.9526,
00205 25.1189,
00206 31.6228,
00207 39.8107,
00208 50.1187,
00209 63.0957,
00210 79.4328,
00211 100,
00212 125.893,
00213 158.489,
00214 199.526,
00215 251.189,
00216 316.228,
00217 398.107,
00218 501.187,
00219 630.957,
00220 794.328,
00221 1000,
00222 1258.93,
00223 1584.89,
00224 1995.26,
00225 2511.89,
00226 3162.28,
00227 3981.07,
00228 5011.87,
00229 6309.57,
00230 7943.28,
00231 10000,
00232 };
00233
00234 Int_t AngleBins = 20;
00235 Double_t AngleEdges[21] = {
00236 -1,
00237 -0.9,
00238 -0.8,
00239 -0.7,
00240 -0.6,
00241 -0.5,
00242 -0.4,
00243 -0.3,
00244 -0.2,
00245 -0.1,
00246 0,
00247 0.1,
00248 0.2,
00249 0.3,
00250 0.4,
00251 0.5,
00252 0.6,
00253 0.7,
00254 0.8,
00255 0.9,
00256 1,
00257 };
00258
00259 Double_t tempFluxRatioNumu[50][20] = {
00260 { 0.720359, 0.793301, 0.899326, 0.923632, 0.88869, 0.890731,
00261 0.895898, 0.941403, 1.00064, 1.07675, 1.05964, 0.95933,
00262 0.889697, 0.824311, 0.772392, 0.747324, 0.708427, 0.700375,
00263 0.683478, 0.669825 },
00264 { 0.811819, 0.858748, 0.928462, 0.962145, 0.951695, 0.935997,
00265 0.972904, 1.0088, 1.0369, 1.11609, 1.0956, 1.04031,
00266 0.97184, 0.910671, 0.87014, 0.831496, 0.824536, 0.79331,
00267 0.784999, 0.772206 },
00268 { 0.901643, 0.938847, 0.986094, 1.00437, 0.980795, 0.992929,
00269 1.00864, 1.02331, 1.05646, 1.11394, 1.13605, 1.07346,
00270 1.00401, 0.957333, 0.934315, 0.909731, 0.905439, 0.893116,
00271 0.888796, 0.887069 },
00272 { 0.976218, 0.986363, 1.00512, 1.00628, 1.02265, 1.00834,
00273 1.01513, 1.03982, 1.04327, 1.06515, 1.06517, 1.04328,
00274 1.01365, 0.979517, 0.966544, 0.969086, 0.96095, 0.961108,
00275 0.969822, 0.949648 },
00276 { 1.01872, 1.01627, 1.02932, 1.03316, 1.03011, 1.02138,
00277 1.01611, 1.01047, 1.02486, 1.05918, 1.02993, 1.03229,
00278 1.01486, 0.998504, 1.01745, 1.00019, 1.0165, 1.00709,
00279 0.998735, 1.01248 },
00280 { 1.04757, 1.03383, 1.02004, 1.02277, 1.02824, 1.02772,
00281 1.01849, 1.0285, 1.00833, 1.03599, 1.05975, 1.01999,
00282 1.02301, 1.01764, 1.02007, 1.02611, 1.0395, 1.04548,
00283 1.05044, 1.04758 },
00284 { 1.04947, 1.05989, 1.0242, 1.03583, 1.03409, 1.03361,
00285 1.01578, 1.04418, 1.04013, 0.997514, 1.03242, 1.03338,
00286 1.04458, 1.03751, 1.05477, 1.04753, 1.05663, 1.05877,
00287 1.05929, 1.06716 },
00288 { 1.07568, 1.0654, 1.03961, 1.03068, 1.04431, 1.02462,
00289 1.02666, 1.05695, 1.03015, 1.09041, 1.04809, 1.03949,
00290 1.05707, 1.04793, 1.07593, 1.06837, 1.06618, 1.08101,
00291 1.08198, 1.09384 },
00292 { 1.06138, 1.06782, 1.04616, 1.04508, 1.04514, 1.03517,
00293 1.05694, 1.03784, 1.06191, 1.0501, 1.09622, 1.06553,
00294 1.07682, 1.06914, 1.07927, 1.082, 1.07267, 1.07801,
00295 1.08774, 1.11025 },
00296 { 1.09412, 1.06854, 1.04611, 1.06524, 1.05229, 1.07934,
00297 1.06993, 1.07127, 1.06371, 1.08644, 1.07538, 1.09415,
00298 1.08042, 1.093, 1.08035, 1.10018, 1.08149, 1.12184,
00299 1.12454, 1.08708 },
00300 { 1.11039, 1.10018, 1.09326, 1.08349, 1.06437, 1.09692,
00301 1.0888, 1.08211, 1.08283, 1.0978, 1.10365, 1.10362,
00302 1.10807, 1.09426, 1.11138, 1.09579, 1.11842, 1.10886,
00303 1.12435, 1.10071 },
00304 { 1.1361, 1.10531, 1.12103, 1.10344, 1.08005, 1.12129,
00305 1.1061, 1.10955, 1.10453, 1.09128, 1.06459, 1.08851,
00306 1.11315, 1.1012, 1.09238, 1.13819, 1.13852, 1.10766,
00307 1.12178, 1.12104 },
00308 { 1.13875, 1.14397, 1.13262, 1.11985, 1.13819, 1.1096,
00309 1.13455, 1.16484, 1.10533, 1.16723, 1.18577, 1.11948,
00310 1.13795, 1.14243, 1.12613, 1.15064, 1.14688, 1.12089,
00311 1.1336, 1.12485 },
00312 { 1.18039, 1.15135, 1.15731, 1.17255, 1.16296, 1.15735,
00313 1.16709, 1.14067, 1.14979, 1.12991, 1.19347, 1.10999,
00314 1.14749, 1.16314, 1.16297, 1.21201, 1.16578, 1.1863,
00315 1.18932, 1.20266 },
00316 { 1.19846, 1.18625, 1.17147, 1.13215, 1.1852, 1.17618,
00317 1.17061, 1.17287, 1.09525, 1.16026, 1.16056, 1.17877,
00318 1.12111, 1.18438, 1.19401, 1.18742, 1.16441, 1.18542,
00319 1.18331, 1.21773 },
00320 { 1.18228, 1.19926, 1.22963, 1.19432, 1.1596, 1.14937,
00321 1.15854, 1.18083, 1.17196, 1.19205, 1.19763, 1.19293,
00322 1.1861, 1.1911, 1.20507, 1.21879, 1.13275, 1.22823,
00323 1.16163, 1.23397 },
00324 { 1.23037, 1.20027, 1.18809, 1.17104, 1.17272, 1.19541,
00325 1.15537, 1.16861, 1.17049, 1.27104, 1.21755, 1.16005,
00326 1.20321, 1.19903, 1.16094, 1.13938, 1.21375, 1.20707,
00327 1.22883, 1.22645 },
00328 { 1.22187, 1.18382, 1.15814, 1.1994, 1.13506, 1.19049,
00329 1.13951, 1.17874, 1.16259, 1.21896, 1.13312, 1.13138,
00330 1.20175, 1.1914, 1.19718, 1.22886, 1.20109, 1.21387,
00331 1.17268, 1.18271 },
00332 { 1.2064, 1.21185, 1.17636, 1.13548, 1.18321, 1.18928,
00333 1.14983, 1.15034, 1.06072, 1.14785, 1.10597, 1.11695,
00334 1.14628, 1.16716, 1.14778, 1.18728, 1.17039, 1.1376,
00335 1.14722, 1.19284 },
00336 { 1.12981, 1.20517, 1.10963, 1.19575, 1.15659, 1.14774,
00337 1.07403, 1.14681, 1.15808, 1.03298, 1.12651, 1.167,
00338 1.15558, 1.08149, 1.16461, 1.19247, 1.18738, 1.24405,
00339 1.19572, 1.0725 },
00340 { 1.1511, 1.15231, 1.14933, 1.14589, 1.14508, 1.1447,
00341 1.14424, 1.13888, 1.12893, 1.11346, 1.11727, 1.13231,
00342 1.13415, 1.14208, 1.14571, 1.14495, 1.14798, 1.14751,
00343 1.15554, 1.14934 },
00344 { 1.14329, 1.1414, 1.14604, 1.14064, 1.13485, 1.13441,
00345 1.13029, 1.12785, 1.1186, 1.10503, 1.1105, 1.11887,
00346 1.12794, 1.13417, 1.13425, 1.13675, 1.13933, 1.14128,
00347 1.14472, 1.14024 },
00348 { 1.13441, 1.14003, 1.13291, 1.13382, 1.13016, 1.13221,
00349 1.12699, 1.12164, 1.11497, 1.10273, 1.10124, 1.11713,
00350 1.1273, 1.12773, 1.13036, 1.12587, 1.13491, 1.13656,
00351 1.13941, 1.13375 },
00352 { 1.12494, 1.12898, 1.13212, 1.12787, 1.12731, 1.12319,
00353 1.12358, 1.11564, 1.1124, 1.09107, 1.09365, 1.11407,
00354 1.12216, 1.1227, 1.1247, 1.12954, 1.12917, 1.13139,
00355 1.12865, 1.1263 },
00356 { 1.12473, 1.12036, 1.12515, 1.12121, 1.12552, 1.11977,
00357 1.12017, 1.11554, 1.11155, 1.10417, 1.10427, 1.11238,
00358 1.11665, 1.11846, 1.123, 1.12356, 1.12608, 1.12536,
00359 1.12016, 1.12335 },
00360 { 1.11929, 1.12043, 1.12485, 1.11824, 1.12171, 1.11996,
00361 1.12309, 1.12174, 1.11513, 1.11046, 1.1087, 1.11655,
00362 1.11417, 1.11852, 1.12062, 1.12059, 1.12305, 1.11882,
00363 1.1239, 1.11862 },
00364 { 1.11655, 1.11637, 1.11695, 1.11822, 1.12575, 1.12162,
00365 1.1204, 1.12111, 1.11985, 1.1116, 1.1141, 1.12162,
00366 1.11969, 1.1153, 1.11833, 1.12399, 1.1219, 1.11628,
00367 1.11423, 1.11661 },
00368 { 1.11135, 1.11244, 1.1154, 1.11615, 1.1088, 1.12134,
00369 1.1183, 1.12032, 1.12219, 1.11493, 1.11796, 1.12142,
00370 1.12113, 1.12485, 1.11878, 1.11905, 1.11432, 1.11331,
00371 1.11249, 1.11162 },
00372 { 1.09605, 1.09877, 1.10767, 1.10823, 1.11795, 1.11448,
00373 1.12246, 1.12601, 1.124, 1.11552, 1.11984, 1.12043,
00374 1.12625, 1.12131, 1.11771, 1.10814, 1.10896, 1.10441,
00375 1.10048, 1.09849 },
00376 { 1.08968, 1.08597, 1.1019, 1.10062, 1.11139, 1.11071,
00377 1.11827, 1.11996, 1.12631, 1.11418, 1.11608, 1.12855,
00378 1.11568, 1.11974, 1.10931, 1.11349, 1.10524, 1.09375,
00379 1.08888, 1.08871 },
00380 { 1.08272, 1.07979, 1.07296, 1.09641, 1.09896, 1.10429,
00381 1.11562, 1.12488, 1.12579, 1.12145, 1.11877, 1.12778,
00382 1.11316, 1.11345, 1.10274, 1.09863, 1.08844, 1.08766,
00383 1.0787, 1.07528 },
00384 { 1.05715, 1.07162, 1.08069, 1.07815, 1.09084, 1.1002,
00385 1.10575, 1.11338, 1.128, 1.11169, 1.11447, 1.12315,
00386 1.1134, 1.10299, 1.10264, 1.09237, 1.07953, 1.06554,
00387 1.07204, 1.05715 },
00388 { 1.05785, 1.06469, 1.06695, 1.06177, 1.07666, 1.08208,
00389 1.10408, 1.11009, 1.1213, 1.10846, 1.11855, 1.12361,
00390 1.11132, 1.09757, 1.07594, 1.08181, 1.07201, 1.07232,
00391 1.05339, 1.05775 },
00392 { 1.04628, 1.04894, 1.04511, 1.05601, 1.06488, 1.07203,
00393 1.09156, 1.10076, 1.12087, 1.11899, 1.10666, 1.11161,
00394 1.10041, 1.0856, 1.0616, 1.06136, 1.06424, 1.04722,
00395 1.03971, 1.04331 },
00396 { 1.02578, 1.03207, 1.03124, 1.03775, 1.04371, 1.06494,
00397 1.06497, 1.09948, 1.10943, 1.10405, 1.10746, 1.10762,
00398 1.09036, 1.07305, 1.06849, 1.04576, 1.04339, 1.02975,
00399 1.02125, 1.03481 },
00400 { 1.02402, 1.01825, 1.02115, 1.03183, 1.03533, 1.04203,
00401 1.06889, 1.07223, 1.09812, 1.10191, 1.10085, 1.10919,
00402 1.09066, 1.08796, 1.03291, 1.03575, 1.05052, 1.01271,
00403 1.00669, 1.01745 },
00404 { 0.97601, 1.01154, 1.02887, 1.00822, 1.01574, 1.04434,
00405 1.03381, 1.06895, 1.09126, 1.08, 1.08503, 1.08315,
00406 1.06469, 1.04636, 1.0384, 1.04315, 1.04777, 0.979799,
00407 1.01124, 0.987096 },
00408 { 0.980698, 1.01862, 1.0055, 1.02462, 0.985128, 0.995888,
00409 1.03434, 1.0326, 1.06082, 1.07983, 1.07919, 1.07433,
00410 1.04541, 1.04537, 1.02427, 1.02044, 0.979039, 0.998043,
00411 1.01255, 0.982553 },
00412 { 0.985548, 1.00081, 1.00725, 1.02392, 1.00049, 1.0117,
00413 1.00962, 1.02788, 1.0681, 1.06668, 1.07183, 1.05702,
00414 1.03429, 1.05098, 1.0235, 1.02786, 0.997753, 1.00043,
00415 1.00099, 0.980654 },
00416 { 0.95215, 0.989738, 0.963026, 0.947634, 0.987205, 1.00397,
00417 1.01354, 1.02783, 1.05208, 1.07734, 1.05458, 1.02989,
00418 1.03, 1.00584, 0.996533, 0.982034, 1.01694, 1.00216,
00419 0.99215, 0.962998 },
00420 { 1.01709, 0.972295, 1.00578, 1.00179, 0.99772, 1.01169,
00421 0.993987, 0.999361, 1.01396, 1.03763, 1.01768, 1.03304,
00422 1.03365, 1.01499, 0.984115, 0.997397, 1.00277, 1.02304,
00423 0.970194, 0.987445 },
00424 { 0.975049, 0.981989, 1.01819, 0.926574, 0.978372, 1.0284,
00425 0.984688, 1.00161, 0.99902, 1.0347, 1.02974, 0.985215,
00426 1.00539, 0.957156, 0.989713, 0.961819, 1.01739, 0.962689,
00427 0.992251, 0.926103 },
00428 { 1.00874, 1.0326, 0.964631, 1.02387, 0.953105, 0.980413,
00429 1.00585, 1.0298, 1.04538, 1.0401, 1.0487, 0.984558,
00430 0.995697, 0.991493, 0.959058, 1.03056, 0.983598, 0.984286,
00431 0.983021, 0.962649 },
00432 { 0.923048, 0.984893, 0.913356, 0.990271, 0.932158, 0.940045,
00433 0.939188, 0.994956, 0.988205, 1.01424, 1.01875, 0.985862,
00434 0.965603, 0.962185, 1.00052, 0.961401, 0.994707, 0.965186,
00435 0.998041, 1.06699 },
00436 { 0.952689, 1.00768, 0.916614, 0.941868, 0.950947, 0.973441,
00437 1.0287, 0.979704, 1.00533, 1.01974, 1.01249, 1.01507,
00438 0.997751, 0.955584, 0.964103, 1.00437, 0.956159, 0.988511,
00439 0.97202, 1.10934 },
00440 { 1.00621, 0.969905, 0.953438, 0.963676, 0.936194, 0.982752,
00441 0.980977, 0.960463, 0.962255, 0.987635, 0.91714, 1.00001,
00442 0.959234, 1.01537, 0.972583, 0.99775, 0.946224, 1.06324,
00443 0.984148, 0.970734 },
00444 { 0.951995, 1.0146, 0.966451, 1.08403, 0.959334, 0.986243,
00445 0.986639, 0.922884, 0.966221, 1.02831, 0.989717, 0.984058,
00446 0.986356, 0.956551, 0.949887, 0.953624, 0.976427, 1.01375,
00447 0.970966, 1.0026 },
00448 { 0.919405, 1.07264, 0.993215, 1.09949, 1.14462, 0.959278,
00449 1.05022, 0.917194, 1.01166, 1.01724, 1.01443, 0.978508,
00450 0.969163, 0.999111, 1.01185, 0.951327, 1.12472, 1.09681,
00451 1.06077, 1.08825 },
00452 { 0.972354, 0.948074, 1.10467, 0.821436, 0.963849, 0.906684,
00453 0.976617, 0.946898, 0.916566, 0.95836, 1.00307, 0.882556,
00454 1.24185, 0.92046, 1.05831, 0.938688, 0.981977, 1.05946,
00455 0.941017, 1.02242 },
00456 { 1.12965, 0.880083, 0.988796, 1.02966, 1.01482, 0.97416,
00457 1.0202, 1.0415, 1.05672, 0.934347, 1.02207, 1.03901,
00458 0.907022, 1.09341, 1.06995, 1.39554, 0.95927, 0.855533,
00459 0.817907, 0.869338 }
00460 };
00461
00462 fHondaFluxRatioNumu = new TH2D("hHondaFluxRatioNumu","",EnergyBins,EnergyEdges,AngleBins,AngleEdges);
00463 for( Int_t nx=0; nx<EnergyBins; nx++ ){
00464 for( Int_t ny=0; ny<AngleBins; ny++ ){
00465 fHondaFluxRatioNumu->SetBinContent(nx+1,ny+1,tempFluxRatioNumu[nx][ny]);
00466 }
00467 }
00468
00469 Double_t tempFluxRatioNumuBar[50][20] = {
00470 { 0.705133, 0.7943, 0.913706, 0.919652, 0.898913, 0.885134,
00471 0.894307, 0.932343, 0.985395, 1.03157, 1.02455, 0.958662,
00472 0.875577, 0.819564, 0.77819, 0.735101, 0.703953, 0.694607,
00473 0.676097, 0.666778 },
00474 { 0.794329, 0.866493, 0.949698, 0.970911, 0.948826, 0.949538,
00475 0.962661, 0.997567, 1.0531, 1.11652, 1.11494, 1.0234,
00476 0.954849, 0.899839, 0.848932, 0.818466, 0.798481, 0.785352,
00477 0.770943, 0.760943 },
00478 { 0.883217, 0.930264, 0.978082, 0.996879, 0.980149, 0.974922,
00479 0.999462, 1.02874, 1.08286, 1.12143, 1.10254, 1.047,
00480 1.00639, 0.94077, 0.930293, 0.899776, 0.881568, 0.884135,
00481 0.867732, 0.870368 },
00482 { 0.946934, 0.985803, 1.01274, 1.01208, 1.01284, 1.01805,
00483 1.00304, 1.03294, 1.05554, 1.0758, 1.12276, 1.04798,
00484 1.02107, 0.983328, 0.968824, 0.960838, 0.948438, 0.960309,
00485 0.954063, 0.953743 },
00486 { 1.00819, 1.01944, 1.01065, 1.01869, 1.01894, 1.01414,
00487 1.0123, 1.03468, 1.03833, 1.03174, 1.04626, 1.04002,
00488 1.01678, 1.014, 1.0033, 0.994082, 0.997439, 0.997324,
00489 1.01387, 1.02144 },
00490 { 1.03904, 1.03371, 1.01483, 1.03546, 1.02987, 1.02239,
00491 1.01477, 1.02417, 1.03699, 1.00359, 1.02797, 1.06239,
00492 1.02949, 1.03873, 1.04383, 1.0426, 1.03902, 1.04585,
00493 1.05546, 1.05332 },
00494 { 1.07246, 1.03317, 1.0204, 1.01604, 1.02214, 1.02477,
00495 1.02828, 1.03273, 1.02501, 1.02918, 1.05316, 1.05169,
00496 1.03752, 1.03116, 1.04735, 1.03585, 1.06831, 1.06725,
00497 1.07953, 1.07072 },
00498 { 1.06979, 1.04662, 1.02333, 1.04155, 1.02398, 1.03551,
00499 1.0422, 1.0317, 1.04971, 1.01329, 0.995565, 1.02665,
00500 1.06475, 1.04793, 1.05069, 1.05534, 1.07181, 1.08366,
00501 1.08138, 1.06107 },
00502 { 1.05827, 1.06336, 1.04431, 1.04371, 1.03866, 1.04363,
00503 1.03875, 1.04221, 1.04986, 1.05065, 1.02655, 1.01735,
00504 1.05924, 1.06284, 1.077, 1.08544, 1.08657, 1.11089,
00505 1.08088, 1.07522 },
00506 { 1.09838, 1.06833, 1.04944, 1.03933, 1.05021, 1.07651,
00507 1.06947, 1.05206, 1.07456, 1.04922, 1.03429, 1.09086,
00508 1.07231, 1.06563, 1.09761, 1.10441, 1.07589, 1.0984,
00509 1.10692, 1.10414 },
00510 { 1.09735, 1.05926, 1.08221, 1.06877, 1.06173, 1.04431,
00511 1.06163, 1.05545, 1.0499, 1.06126, 1.07627, 1.09785,
00512 1.08965, 1.0831, 1.08786, 1.09269, 1.07797, 1.09777,
00513 1.10847, 1.07936 },
00514 { 1.10494, 1.10999, 1.09035, 1.10165, 1.0757, 1.08716,
00515 1.11784, 1.10386, 1.07887, 1.05409, 1.06081, 1.07985,
00516 1.06992, 1.07563, 1.09316, 1.11355, 1.08677, 1.10125,
00517 1.10858, 1.10307 },
00518 { 1.06633, 1.10916, 1.08729, 1.0818, 1.10057, 1.11147,
00519 1.11002, 1.11307, 1.1051, 1.03308, 1.10983, 1.12574,
00520 1.11438, 1.12275, 1.12051, 1.11036, 1.13027, 1.09293,
00521 1.11491, 1.11677 },
00522 { 1.16918, 1.06589, 1.11561, 1.0767, 1.1046, 1.09654,
00523 1.09506, 1.07667, 1.11686, 1.06134, 1.07642, 1.0963,
00524 1.11021, 1.10293, 1.11079, 1.13937, 1.11544, 1.12877,
00525 1.10403, 1.15176 },
00526 { 1.13041, 1.12844, 1.11995, 1.11652, 1.103, 1.10452,
00527 1.08008, 1.11132, 1.11674, 1.04467, 1.0198, 1.06417,
00528 1.13236, 1.12685, 1.10334, 1.10725, 1.13853, 1.11495,
00529 1.12575, 1.11826 },
00530 { 1.13562, 1.13881, 1.14526, 1.12451, 1.09084, 1.13981,
00531 1.13376, 1.13784, 1.11226, 1.20303, 1.18051, 1.11957,
00532 1.14898, 1.12047, 1.10479, 1.16525, 1.14292, 1.11555,
00533 1.1329, 1.11313 },
00534 { 1.10223, 1.12731, 1.11671, 1.2029, 1.11916, 1.13736,
00535 1.12954, 1.18114, 1.14328, 1.18401, 1.14757, 1.13574,
00536 1.1303, 1.11546, 1.16661, 1.08987, 1.09791, 1.17782,
00537 1.11611, 1.1323 },
00538 { 1.12881, 1.17714, 1.11728, 1.1694, 1.17252, 1.14896,
00539 1.17745, 1.14006, 1.1291, 1.22367, 1.1102, 1.13685,
00540 1.09323, 1.13254, 1.12325, 1.12837, 1.14189, 1.14157,
00541 1.14519, 1.14733 },
00542 { 1.16973, 1.16469, 1.17574, 1.13778, 1.1252, 1.15677,
00543 1.189, 1.11608, 1.11831, 1.10882, 1.15818, 1.20677,
00544 1.13892, 1.17439, 1.16413, 1.13288, 1.10678, 1.11631,
00545 1.08629, 1.1376 },
00546 { 1.12826, 1.10288, 1.21341, 1.14021, 1.1797, 1.1559,
00547 1.16067, 1.13227, 1.14451, 1.15918, 1.21151, 1.11502,
00548 1.15246, 1.06628, 1.1391, 1.15752, 1.10684, 1.18234,
00549 1.15648, 1.13509 },
00550 { 1.14934, 1.14353, 1.14833, 1.13995, 1.13809, 1.13893,
00551 1.13442, 1.13488, 1.1259, 1.11419, 1.11284, 1.13253,
00552 1.13684, 1.1381, 1.14091, 1.14185, 1.14264, 1.14149,
00553 1.14814, 1.14853 },
00554 { 1.16095, 1.1593, 1.16022, 1.15405, 1.15332, 1.15007,
00555 1.14307, 1.14341, 1.13892, 1.1242, 1.12678, 1.13922,
00556 1.14437, 1.14569, 1.15234, 1.15359, 1.15419, 1.1595,
00557 1.16099, 1.1588 },
00558 { 1.17996, 1.17918, 1.17836, 1.16848, 1.17065, 1.16373,
00559 1.15307, 1.15678, 1.14122, 1.13179, 1.13469, 1.14693,
00560 1.15474, 1.15924, 1.1636, 1.16438, 1.17139, 1.17444,
00561 1.18003, 1.17848 },
00562 { 1.20275, 1.20148, 1.20205, 1.19573, 1.19472, 1.18742,
00563 1.17482, 1.16685, 1.1545, 1.12237, 1.13158, 1.16585,
00564 1.17009, 1.17332, 1.17742, 1.1848, 1.19301, 1.20006,
00565 1.20284, 1.20302 },
00566 { 1.23199, 1.23005, 1.22551, 1.216, 1.20596, 1.20472,
00567 1.19453, 1.18632, 1.17557, 1.1646, 1.16352, 1.17666,
00568 1.18679, 1.19638, 1.20046, 1.2107, 1.21719, 1.22444,
00569 1.22943, 1.23402 },
00570 { 1.26679, 1.25808, 1.25414, 1.24784, 1.23604, 1.22859,
00571 1.21677, 1.20914, 1.19875, 1.17636, 1.18044, 1.19569,
00572 1.20666, 1.2204, 1.22899, 1.23921, 1.24441, 1.25865,
00573 1.26412, 1.2668 },
00574 { 1.30729, 1.3022, 1.29234, 1.28425, 1.26724, 1.25899,
00575 1.25015, 1.22934, 1.22055, 1.19618, 1.19642, 1.21462,
00576 1.23154, 1.24583, 1.25785, 1.26872, 1.27616, 1.2898,
00577 1.3003, 1.31289 },
00578 { 1.3458, 1.33346, 1.32683, 1.31017, 1.30516, 1.29163,
00579 1.2781, 1.26183, 1.23963, 1.21289, 1.20949, 1.2383,
00580 1.25983, 1.27544, 1.28205, 1.30027, 1.32064, 1.32759,
00581 1.33107, 1.35272 },
00582 { 1.39339, 1.37818, 1.36722, 1.35793, 1.34069, 1.32341,
00583 1.30464, 1.28222, 1.26134, 1.23086, 1.23415, 1.26367,
00584 1.28575, 1.30516, 1.32728, 1.33689, 1.36279, 1.36977,
00585 1.38391, 1.4021 },
00586 { 1.43986, 1.42004, 1.41459, 1.39846, 1.38149, 1.37117,
00587 1.35151, 1.32025, 1.29475, 1.24535, 1.24385, 1.29182,
00588 1.31567, 1.34317, 1.3619, 1.37828, 1.38325, 1.40652,
00589 1.42174, 1.43721 },
00590 { 1.47846, 1.46473, 1.44732, 1.44225, 1.42023, 1.4048,
00591 1.38004, 1.35857, 1.32158, 1.265, 1.26025, 1.32548,
00592 1.35721, 1.37355, 1.39434, 1.42195, 1.44374, 1.45207,
00593 1.46786, 1.47445 },
00594 { 1.52119, 1.51495, 1.5031, 1.47081, 1.45149, 1.43951,
00595 1.40974, 1.37628, 1.34304, 1.29201, 1.28083, 1.33835,
00596 1.38144, 1.41119, 1.43139, 1.44877, 1.48598, 1.49132,
00597 1.51611, 1.53203 },
00598 { 1.56421, 1.54083, 1.54527, 1.51762, 1.50602, 1.47605,
00599 1.44575, 1.40985, 1.37555, 1.30113, 1.30346, 1.37167,
00600 1.41186, 1.44981, 1.4752, 1.50151, 1.53919, 1.53593,
00601 1.55018, 1.56263 },
00602 { 1.60936, 1.61102, 1.58039, 1.55951, 1.53931, 1.51283,
00603 1.47672, 1.43874, 1.4091, 1.31694, 1.31639, 1.41254,
00604 1.45811, 1.48385, 1.51517, 1.53574, 1.55859, 1.57769,
00605 1.58624, 1.62263 },
00606 { 1.6816, 1.60386, 1.63108, 1.63122, 1.56379, 1.55547,
00607 1.52758, 1.47433, 1.42696, 1.33888, 1.34102, 1.43307,
00608 1.49232, 1.52119, 1.53433, 1.5739, 1.59573, 1.60716,
00609 1.62192, 1.64429 },
00610 { 1.73812, 1.66965, 1.67677, 1.6745, 1.6444, 1.57884,
00611 1.53448, 1.52977, 1.47095, 1.36072, 1.35836, 1.48096,
00612 1.50554, 1.56879, 1.60711, 1.5996, 1.64558, 1.65923,
00613 1.68725, 1.69653 },
00614 { 1.73741, 1.74946, 1.71231, 1.68992, 1.63998, 1.647,
00615 1.58153, 1.58315, 1.51281, 1.38858, 1.38779, 1.50532,
00616 1.55179, 1.57645, 1.6258, 1.67598, 1.69954, 1.7439,
00617 1.72433, 1.721 },
00618 { 1.72621, 1.7482, 1.75391, 1.71555, 1.69618, 1.6894,
00619 1.62061, 1.59407, 1.55012, 1.43062, 1.42841, 1.55488,
00620 1.61292, 1.63951, 1.69458, 1.70176, 1.69348, 1.74217,
00621 1.72865, 1.76461 },
00622 { 1.81684, 1.82032, 1.80797, 1.74448, 1.7608, 1.68778,
00623 1.68566, 1.64367, 1.57199, 1.44767, 1.43786, 1.56677,
00624 1.61475, 1.66566, 1.68043, 1.76808, 1.74716, 1.80031,
00625 1.78198, 1.83992 },
00626 { 1.87268, 1.81514, 1.81341, 1.75813, 1.70021, 1.77617,
00627 1.71203, 1.68026, 1.60029, 1.47949, 1.47206, 1.58531,
00628 1.66951, 1.68274, 1.74688, 1.77939, 1.76911, 1.80966,
00629 1.7994, 1.87301 },
00630 { 1.82806, 1.82168, 1.87607, 1.82485, 1.76238, 1.72373,
00631 1.70933, 1.74603, 1.63884, 1.52924, 1.53744, 1.66926,
00632 1.69867, 1.77375, 1.80192, 1.79363, 1.86116, 1.80724,
00633 1.84731, 1.8975 },
00634 { 1.90363, 1.87933, 1.79027, 1.82264, 1.88811, 1.79203,
00635 1.78229, 1.70044, 1.69824, 1.52854, 1.58538, 1.67855,
00636 1.70352, 1.85734, 1.837, 1.91832, 1.84225, 1.93317,
00637 1.90808, 2.03093 },
00638 { 1.96605, 1.9209, 1.89505, 1.86872, 1.85034, 1.96523,
00639 1.84052, 1.76761, 1.72045, 1.58447, 1.61228, 1.7973,
00640 1.79289, 1.79712, 1.82767, 1.9419, 1.91372, 1.9548,
00641 1.99333, 1.92723 },
00642 { 2.00469, 2.07578, 2.02333, 1.92862, 1.93283, 1.9165,
00643 1.90044, 1.82469, 1.77742, 1.64828, 1.64222, 1.77863,
00644 1.80037, 1.83077, 1.88771, 2.03012, 1.87732, 1.74876,
00645 1.91666, 1.96412 },
00646 { 1.90766, 2.06549, 1.98681, 1.88428, 1.84486, 1.94434,
00647 1.81099, 1.85154, 1.69809, 1.69709, 1.66751, 1.73816,
00648 1.7965, 2.00486, 1.79924, 1.90534, 1.82175, 1.94305,
00649 2.1263, 1.90907 },
00650 { 1.87915, 1.99101, 1.88399, 2.09318, 1.87379, 2.06238,
00651 1.97225, 1.88009, 1.71774, 1.76548, 1.71688, 1.74072,
00652 1.9428, 1.97727, 2.04047, 1.94044, 1.91399, 2.05775,
00653 1.99772, 1.98832 },
00654 { 2.17599, 1.72374, 1.80553, 1.99612, 1.90734, 1.89137,
00655 1.87229, 1.98052, 1.84842, 1.667, 1.71548, 1.80627,
00656 1.96131, 1.79501, 2.00485, 1.87175, 2.08943, 2.15724,
00657 2.09258, 1.83488 },
00658 { 1.61225, 2.3388, 1.9408, 2.00414, 1.94543, 1.86327,
00659 1.85944, 2.12909, 1.93555, 1.59297, 1.84128, 1.85451,
00660 2.00904, 2.01206, 2.0765, 1.90297, 1.85861, 1.78392,
00661 1.9965, 1.95516 },
00662 { 1.69503, 2.17423, 2.08416, 2.03891, 1.70623, 2.03858,
00663 1.97139, 1.72152, 1.83812, 1.66838, 1.84196, 1.91136,
00664 1.93389, 2.00585, 2.03476, 1.96987, 1.80645, 1.93991,
00665 1.57272, 1.90624 },
00666 { 1.80686, 2.75457, 1.8771, 1.99271, 1.82314, 1.73241,
00667 1.79203, 1.7553, 1.96554, 1.69902, 1.7991, 1.80475,
00668 1.86286, 2.06388, 1.76461, 1.90267, 2.58264, 2.0764,
00669 1.767, 1.58713 }
00670 };
00671
00672 fHondaFluxRatioNumuBar = new TH2D("hHondaFluxRatioNumuBar","",EnergyBins,EnergyEdges,AngleBins,AngleEdges);
00673 for( Int_t nx=0; nx<EnergyBins; nx++ ){
00674 for( Int_t ny=0; ny<AngleBins; ny++ ){
00675 fHondaFluxRatioNumuBar->SetBinContent(nx+1,ny+1,tempFluxRatioNumuBar[nx][ny]);
00676 }
00677 }
00678
00679 return;
00680 }
00681
00682 void AtNuFluxRW::LoadFlukaFluxRatios()
00683 {
00684
00685
00686
00687
00688
00689 Int_t EnergyBins = 30;
00690 Double_t EnergyEdges[31] = {
00691 0.1,
00692 0.125893,
00693 0.158489,
00694 0.199526,
00695 0.251189,
00696 0.316228,
00697 0.398107,
00698 0.501187,
00699 0.630957,
00700 0.794328,
00701 1,
00702 1.25893,
00703 1.58489,
00704 1.99526,
00705 2.51189,
00706 3.16228,
00707 3.98107,
00708 5.01187,
00709 6.30957,
00710 7.94328,
00711 10,
00712 12.5893,
00713 15.8489,
00714 19.9526,
00715 25.1189,
00716 31.6228,
00717 39.8107,
00718 50.1187,
00719 63.0957,
00720 79.4328,
00721 100,
00722 };
00723
00724 Int_t AngleBins = 20;
00725 Double_t AngleEdges[21] = {
00726 -1,
00727 -0.9,
00728 -0.8,
00729 -0.7,
00730 -0.6,
00731 -0.5,
00732 -0.4,
00733 -0.3,
00734 -0.2,
00735 -0.1,
00736 0,
00737 0.1,
00738 0.2,
00739 0.3,
00740 0.4,
00741 0.5,
00742 0.6,
00743 0.7,
00744 0.8,
00745 0.9,
00746 1,
00747 };
00748
00749 Double_t tempFluxRatioNumu[30][20] = {
00750 { 0.840422, 0.87705, 0.934688, 0.976589, 0.974876, 0.997133,
00751 0.983382, 0.992384, 1.04525, 1.15244, 1.14142, 1.04213,
00752 0.984293, 0.927369, 0.881852, 0.867163, 0.830181, 0.827924,
00753 0.816066, 0.808012 },
00754 { 0.926969, 0.940956, 0.967104, 1.01872, 1.04036, 1.03584,
00755 1.05436, 1.04804, 1.06509, 1.17172, 1.15764, 1.10472,
00756 1.05144, 1.00346, 0.972511, 0.946645, 0.947392, 0.918824,
00757 0.916891, 0.908118 },
00758 { 0.97895, 0.992491, 1.00549, 1.04175, 1.04758, 1.0708,
00759 1.06923, 1.04677, 1.07574, 1.16438, 1.19801, 1.12455,
00760 1.06355, 1.02451, 1.00821, 0.993298, 0.995052, 0.984066,
00761 0.984916, 0.98428 },
00762 { 1.01067, 1.00498, 0.997878, 1.02038, 1.06626, 1.06147,
00763 1.05772, 1.05668, 1.06735, 1.13218, 1.14159, 1.09961,
00764 1.06503, 1.02501, 1.01173, 1.01959, 1.01161, 1.01055,
00765 1.02264, 0.999025 },
00766 { 1.02615, 1.00925, 1.00138, 1.02907, 1.05666, 1.05908,
00767 1.04621, 1.02356, 1.05518, 1.13481, 1.11541, 1.09418,
00768 1.0616, 1.03085, 1.04407, 1.02651, 1.04104, 1.02927,
00769 1.0223, 1.03323 },
00770 { 1.03955, 1.01353, 0.98098, 1.00772, 1.0415, 1.05275,
00771 1.03635, 1.03546, 1.03972, 1.10803, 1.14554, 1.08155,
00772 1.06196, 1.03651, 1.03095, 1.03401, 1.04684, 1.05032,
00773 1.05692, 1.05151 },
00774 { 1.02609, 1.02158, 0.969828, 1.00484, 1.03046, 1.03964,
00775 1.01265, 1.03366, 1.0545, 1.04331, 1.09412, 1.07823,
00776 1.06086, 1.03345, 1.04625, 1.03258, 1.04186, 1.04311,
00777 1.04528, 1.05227 },
00778 { 1.02814, 1.00388, 0.964053, 0.978136, 1.01453, 1.00621,
00779 0.993931, 1.01695, 1.01253, 1.09057, 1.06135, 1.04659,
00780 1.0381, 1.01049, 1.03628, 1.02236, 1.02211, 1.03681,
00781 1.03906, 1.05217 },
00782 { 0.985539, 0.978484, 0.948651, 0.966236, 0.984528, 0.983144,
00783 0.990259, 0.965815, 1.00339, 0.997712, 1.05529, 1.02844,
00784 1.01774, 0.994417, 1.00558, 1.00162, 0.993757, 1.00108,
00785 1.01096, 1.03329 },
00786 { 0.982514, 0.948947, 0.924249, 0.95685, 0.956261, 0.990715,
00787 0.963801, 0.959075, 0.96147, 0.984893, 0.984523, 1.00368,
00788 0.976896, 0.977057, 0.972528, 0.981443, 0.967363, 1.00533,
00789 1.0076, 0.974384 },
00790 { 0.95993, 0.944604, 0.938522, 0.942052, 0.93041, 0.968062,
00791 0.944421, 0.929909, 0.942288, 0.955893, 0.972053, 0.962979,
00792 0.958358, 0.938964, 0.960561, 0.940459, 0.960973, 0.95571,
00793 0.967179, 0.947774 },
00794 { 0.940157, 0.908715, 0.928857, 0.922184, 0.905714, 0.948219,
00795 0.919465, 0.916848, 0.916737, 0.91332, 0.898587, 0.905837,
00796 0.920511, 0.906117, 0.906796, 0.93594, 0.938358, 0.915237,
00797 0.923603, 0.924227 },
00798 { 0.906709, 0.906795, 0.907644, 0.902078, 0.916423, 0.899413,
00799 0.907459, 0.92709, 0.881327, 0.934491, 0.957184, 0.890575,
00800 0.902336, 0.904077, 0.898647, 0.909496, 0.90791, 0.891161,
00801 0.898456, 0.893693 },
00802 { 0.909814, 0.884431, 0.896884, 0.910102, 0.901812, 0.902066,
00803 0.902181, 0.881045, 0.886066, 0.877366, 0.929636, 0.854076,
00804 0.882906, 0.89292, 0.894411, 0.925695, 0.891105, 0.913367,
00805 0.912787, 0.923982 },
00806 { 0.900552, 0.890072, 0.878621, 0.84789, 0.89215, 0.886942,
00807 0.883239, 0.88863, 0.825906, 0.881006, 0.884023, 0.889832,
00808 0.847269, 0.889936, 0.893078, 0.884741, 0.867299, 0.890661,
00809 0.888197, 0.914311 },
00810 { 0.881396, 0.896532, 0.901616, 0.870362, 0.854736, 0.849712,
00811 0.865889, 0.889759, 0.864096, 0.8929, 0.902701, 0.884834,
00812 0.897773, 0.890881, 0.888527, 0.900204, 0.836041, 0.917734,
00813 0.869007, 0.922719 },
00814 { 0.90194, 0.881014, 0.917437, 0.889416, 0.880151, 0.894467,
00815 0.874528, 0.883313, 0.864998, 0.943504, 0.906853, 0.856194,
00816 0.910856, 0.898539, 0.859745, 0.843681, 0.90646, 0.917322,
00817 0.895364, 0.895972 },
00818 { 0.877034, 0.853173, 0.87813, 0.888662, 0.836146, 0.881792,
00819 0.838243, 0.880982, 0.84918, 0.889483, 0.830551, 0.829126,
00820 0.893938, 0.874343, 0.880488, 0.899203, 0.881769, 0.914314,
00821 0.840759, 0.842651 },
00822 { 0.861046, 0.868173, 0.882498, 0.828324, 0.867634, 0.875882,
00823 0.83913, 0.853538, 0.764237, 0.835009, 0.804898, 0.804856,
00824 0.847978, 0.849124, 0.839592, 0.865624, 0.85032, 0.85092,
00825 0.818272, 0.84908 },
00826 { 0.814619, 0.868039, 0.829126, 0.866669, 0.849354, 0.844058,
00827 0.782655, 0.857184, 0.8326, 0.746175, 0.810921, 0.838383,
00828 0.861179, 0.787891, 0.854677, 0.874375, 0.858519, 0.929754,
00829 0.859232, 0.7724 },
00830 { 0.839912, 0.838151, 0.854668, 0.826393, 0.84484, 0.841833,
00831 0.838396, 0.861452, 0.806416, 0.795237, 0.795214, 0.807704,
00832 0.856431, 0.83557, 0.841697, 0.844222, 0.828263, 0.853307,
00833 0.838909, 0.838539 },
00834 { 0.841239, 0.835541, 0.844869, 0.82072, 0.840709, 0.834053,
00835 0.83395, 0.860885, 0.796785, 0.779402, 0.781053, 0.795339,
00836 0.859116, 0.835604, 0.833107, 0.839841, 0.820427, 0.84344,
00837 0.836738, 0.841507 },
00838 { 0.838774, 0.835037, 0.829114, 0.814417, 0.836523, 0.826442,
00839 0.834585, 0.862857, 0.784558, 0.762781, 0.762801, 0.784694,
00840 0.862474, 0.834314, 0.823419, 0.832056, 0.814128, 0.831288,
00841 0.835001, 0.840306 },
00842 { 0.832603, 0.825802, 0.819527, 0.807237, 0.826646, 0.813892,
00843 0.835197, 0.859745, 0.777169, 0.743549, 0.745735, 0.776156,
00844 0.860269, 0.835033, 0.814564, 0.828506, 0.810048, 0.819239,
00845 0.824896, 0.833933 },
00846 { 0.823636, 0.813502, 0.806879, 0.800983, 0.818893, 0.803106,
00847 0.833114, 0.84876, 0.766483, 0.727132, 0.72711, 0.767074,
00848 0.850802, 0.831847, 0.805585, 0.817778, 0.804411, 0.807003,
00849 0.814046, 0.82229 },
00850 { 0.804692, 0.799141, 0.794776, 0.790953, 0.803637, 0.789965,
00851 0.826226, 0.834223, 0.755414, 0.709901, 0.70875, 0.756415,
00852 0.830037, 0.822851, 0.790829, 0.80359, 0.794334, 0.790506,
00853 0.801947, 0.803224 },
00854 { 0.782066, 0.776863, 0.777641, 0.779029, 0.789286, 0.775381,
00855 0.809164, 0.8014, 0.741867, 0.688316, 0.689865, 0.743061,
00856 0.801416, 0.805482, 0.773585, 0.789048, 0.781586, 0.777162,
00857 0.775198, 0.780561 },
00858 { 0.752632, 0.749598, 0.762307, 0.757796, 0.758955, 0.756511,
00859 0.782339, 0.75787, 0.721343, 0.666258, 0.668067, 0.720849,
00860 0.758386, 0.786926, 0.755133, 0.766845, 0.756555, 0.760879,
00861 0.748924, 0.750998 },
00862 { 0.711723, 0.711968, 0.738055, 0.724103, 0.741876, 0.72855,
00863 0.750054, 0.705974, 0.695288, 0.641971, 0.644456, 0.693081,
00864 0.704613, 0.749283, 0.730701, 0.735672, 0.724576, 0.735884,
00865 0.71195, 0.711646 },
00866 { 0.669569, 0.671043, 0.714502, 0.677259, 0.711067, 0.698983,
00867 0.701589, 0.636847, 0.662751, 0.6144, 0.615449, 0.664066,
00868 0.631455, 0.702509, 0.697745, 0.711669, 0.680096, 0.709219,
00869 0.671501, 0.667969 }
00870 };
00871
00872 fFlukaFluxRatioNumu = new TH2D("hFlukaFluxRatioNumu","",EnergyBins,EnergyEdges,AngleBins,AngleEdges);
00873 for( Int_t nx=0; nx<EnergyBins; nx++ ){
00874 for( Int_t ny=0; ny<AngleBins; ny++ ){
00875 fFlukaFluxRatioNumu->SetBinContent(nx+1,ny+1,tempFluxRatioNumu[nx][ny]);
00876 }
00877 }
00878
00879 Double_t tempFluxRatioNumuBar[30][20] = {
00880 { 0.813955, 0.880554, 0.953715, 0.963695, 0.979551, 0.977534,
00881 0.972933, 0.977053, 1.01003, 1.09478, 1.0922, 1.02583,
00882 0.958933, 0.91476, 0.877688, 0.842443, 0.81492, 0.812994,
00883 0.799772, 0.794956 },
00884 { 0.906648, 0.953322, 0.995764, 1.02309, 1.03268, 1.04295,
00885 1.03914, 1.03285, 1.06927, 1.16189, 1.16785, 1.0793,
00886 1.03057, 0.99315, 0.946973, 0.92694, 0.914156, 0.907307,
00887 0.898265, 0.892748 },
00888 { 0.967687, 0.989562, 1.00662, 1.03495, 1.04794, 1.05033,
00889 1.0605, 1.04846, 1.09523, 1.16135, 1.14993, 1.09518,
00890 1.06802, 1.01583, 1.00861, 0.985684, 0.973646, 0.981224,
00891 0.967247, 0.972892 },
00892 { 0.991858, 1.01151, 1.01583, 1.02905, 1.0576, 1.07458,
00893 1.04532, 1.04399, 1.07336, 1.12285, 1.18183, 1.09804,
00894 1.07076, 1.03873, 1.02046, 1.01448, 1.00555, 1.01874,
00895 1.01239, 1.01322 },
00896 { 1.02112, 1.01764, 0.992235, 1.01803, 1.04397, 1.04991,
00897 1.0389, 1.03678, 1.06082, 1.08782, 1.11203, 1.09425,
00898 1.05631, 1.04908, 1.02949, 1.01806, 1.0214, 1.0204,
00899 1.03596, 1.04356 },
00900 { 1.03206, 1.01305, 0.982932, 1.02285, 1.03961, 1.0437,
00901 1.02768, 1.01968, 1.06451, 1.0659, 1.10116, 1.12273,
00902 1.0602, 1.0577, 1.05059, 1.04353, 1.03921, 1.04453,
00903 1.05402, 1.0504 },
00904 { 1.04662, 0.997828, 0.974034, 0.990587, 1.01547, 1.03001,
00905 1.02316, 1.01201, 1.04113, 1.07662, 1.11572, 1.10016,
00906 1.05026, 1.02914, 1.03481, 1.01451, 1.04395, 1.04388,
00907 1.05508, 1.04726 },
00908 { 1.02127, 0.992092, 0.956757, 0.996096, 0.993806, 1.01754,
00909 1.00917, 0.987223, 1.03681, 1.03304, 1.0268, 1.0475,
00910 1.04834, 1.01599, 1.01068, 1.00583, 1.0196, 1.03302,
00911 1.03161, 1.01409 },
00912 { 0.984216, 0.982635, 0.954217, 0.975922, 0.982448, 0.993121,
00913 0.973113, 0.963795, 0.998541, 1.02415, 1.01221, 0.995614,
00914 1.00338, 0.993034, 1.00459, 1.00257, 1.00316, 1.02745,
00915 1.00104, 0.997168 },
00916 { 0.98761, 0.958462, 0.934194, 0.942014, 0.962714, 0.989743,
00917 0.96351, 0.937555, 0.975999, 0.969581, 0.9667, 1.01405,
00918 0.969301, 0.955468, 0.986127, 0.985209, 0.959008, 0.979765,
00919 0.988938, 0.98684 },
00920 { 0.950241, 0.917507, 0.934231, 0.937153, 0.936406, 0.922287,
00921 0.917466, 0.902756, 0.907597, 0.939338, 0.958714, 0.9671,
00922 0.9429, 0.929066, 0.93942, 0.937493, 0.922468, 0.941181,
00923 0.949811, 0.928142 },
00924 { 0.919697, 0.924763, 0.908541, 0.931705, 0.912169, 0.921021,
00925 0.931278, 0.909739, 0.894162, 0.891666, 0.903067, 0.906338,
00926 0.885935, 0.889556, 0.908473, 0.922189, 0.895008, 0.909152,
00927 0.914509, 0.913367 },
00928 { 0.861536, 0.896088, 0.880956, 0.886718, 0.904485, 0.906834,
00929 0.893743, 0.888632, 0.879259, 0.840086, 0.907926, 0.900853,
00930 0.887874, 0.895991, 0.897546, 0.891255, 0.900171, 0.873147,
00931 0.890441, 0.897142 },
00932 { 0.925017, 0.841357, 0.883133, 0.857694, 0.884814, 0.867021,
00933 0.86094, 0.840496, 0.862228, 0.840947, 0.855974, 0.847824,
00934 0.863373, 0.859504, 0.865201, 0.894354, 0.867702, 0.880532,
00935 0.862199, 0.906589 },
00936 { 0.883853, 0.876899, 0.869666, 0.867285, 0.868957, 0.851473,
00937 0.838532, 0.859709, 0.849183, 0.811147, 0.798071, 0.808903,
00938 0.8743, 0.867519, 0.84239, 0.857428, 0.874025, 0.857625,
00939 0.869227, 0.870517 },
00940 { 0.893944, 0.884452, 0.879849, 0.858588, 0.856081, 0.866517,
00941 0.881239, 0.887897, 0.846837, 0.923903, 0.911188, 0.856679,
00942 0.901609, 0.868033, 0.839657, 0.903914, 0.878702, 0.85994,
00943 0.878115, 0.874336 },
00944 { 0.876379, 0.903858, 0.875289, 0.936114, 0.892408, 0.885688,
00945 0.878715, 0.932544, 0.891969, 0.902956, 0.876699, 0.888734,
00946 0.891627, 0.863273, 0.897917, 0.854565, 0.842009, 0.916724,
00947 0.8826, 0.896507 },
00948 { 0.878963, 0.929851, 0.866767, 0.897388, 0.911348, 0.879581,
00949 0.897608, 0.882406, 0.866203, 0.938416, 0.85367, 0.875452,
00950 0.843557, 0.859121, 0.857549, 0.874313, 0.870554, 0.88027,
00951 0.896112, 0.890564 },
00952 { 0.897591, 0.908135, 0.904751, 0.869448, 0.865583, 0.880213,
00953 0.893247, 0.849962, 0.843268, 0.861255, 0.898265, 0.910786,
00954 0.866673, 0.879821, 0.88321, 0.868881, 0.844031, 0.85674,
00955 0.843284, 0.872949 },
00956 { 0.860163, 0.853109, 0.92984, 0.872262, 0.90232, 0.877299,
00957 0.866028, 0.855133, 0.853504, 0.901105, 0.94167, 0.829629,
00958 0.869636, 0.793252, 0.864147, 0.884708, 0.84564, 0.906451,
00959 0.892607, 0.86357 },
00960 { 0.867831, 0.877496, 0.875074, 0.870433, 0.866616, 0.861246,
00961 0.843616, 0.852774, 0.831459, 0.87798, 0.876609, 0.831759,
00962 0.853022, 0.842539, 0.864818, 0.868834, 0.871198, 0.871366,
00963 0.880759, 0.867405 },
00964 { 0.869339, 0.879671, 0.876674, 0.878389, 0.873719, 0.867182,
00965 0.847906, 0.852922, 0.829993, 0.881279, 0.882081, 0.827769,
00966 0.853044, 0.846821, 0.866685, 0.874014, 0.878079, 0.877752,
00967 0.880541, 0.868112 },
00968 { 0.877093, 0.8836, 0.880889, 0.883811, 0.881232, 0.869681,
00969 0.850304, 0.857053, 0.824268, 0.881685, 0.881795, 0.825359,
00970 0.856492, 0.852012, 0.86989, 0.877108, 0.886399, 0.880586,
00971 0.884923, 0.873408 },
00972 { 0.88204, 0.890896, 0.885366, 0.891302, 0.885061, 0.87298,
00973 0.855152, 0.862005, 0.821338, 0.873522, 0.877379, 0.822899,
00974 0.859078, 0.855205, 0.870471, 0.885673, 0.888663, 0.886682,
00975 0.889588, 0.88343 },
00976 { 0.887653, 0.897806, 0.890016, 0.887905, 0.888866, 0.871687,
00977 0.855953, 0.866537, 0.8189, 0.865546, 0.864905, 0.819642,
00978 0.866854, 0.858229, 0.868566, 0.892356, 0.888768, 0.890158,
00979 0.897351, 0.889105 },
00980 { 0.892832, 0.905664, 0.893565, 0.884684, 0.894954, 0.865424,
00981 0.854828, 0.872997, 0.815889, 0.843188, 0.846302, 0.81376,
00982 0.871165, 0.85856, 0.865709, 0.897249, 0.882239, 0.896822,
00983 0.910005, 0.892842 },
00984 { 0.894064, 0.923836, 0.900625, 0.874809, 0.895932, 0.858614,
00985 0.857579, 0.876727, 0.809252, 0.816875, 0.817175, 0.805324,
00986 0.878275, 0.855732, 0.857834, 0.896981, 0.869297, 0.89787,
00987 0.922487, 0.89789 },
00988 { 0.886917, 0.932699, 0.901197, 0.849552, 0.898207, 0.847619,
00989 0.853964, 0.886825, 0.796962, 0.779662, 0.77746, 0.796111,
00990 0.885423, 0.852868, 0.841331, 0.894837, 0.856341, 0.900216,
00991 0.931032, 0.89148 },
00992 { 0.870615, 0.948943, 0.900949, 0.828719, 0.893579, 0.827786,
00993 0.848829, 0.887204, 0.78355, 0.73544, 0.737145, 0.784997,
00994 0.889646, 0.848916, 0.830208, 0.891051, 0.831688, 0.901784,
00995 0.952889, 0.876058 },
00996 { 0.843181, 0.959424, 0.895298, 0.794099, 0.883033, 0.813929,
00997 0.851788, 0.896845, 0.77083, 0.683374, 0.681968, 0.769085,
00998 0.893734, 0.844724, 0.808426, 0.880979, 0.785459, 0.891535,
00999 0.960569, 0.841634 }
01000 };
01001
01002 fFlukaFluxRatioNumuBar = new TH2D("hFlukaFluxRatioNumuBar","",EnergyBins,EnergyEdges,AngleBins,AngleEdges);
01003 for( Int_t nx=0; nx<EnergyBins; nx++ ){
01004 for( Int_t ny=0; ny<AngleBins; ny++ ){
01005 fFlukaFluxRatioNumuBar->SetBinContent(nx+1,ny+1,tempFluxRatioNumuBar[nx][ny]);
01006 }
01007 }
01008
01009 return;
01010 }
01011
01012 void AtNuFluxRW::LoadBartolFluxRatios()
01013 {
01014
01015
01016
01017
01018
01019 Int_t EnergyBins = 50;
01020 Double_t EnergyEdges[51] = {
01021 0.1,
01022 0.125893,
01023 0.158489,
01024 0.199526,
01025 0.251189,
01026 0.316228,
01027 0.398107,
01028 0.501187,
01029 0.630957,
01030 0.794328,
01031 1,
01032 1.25893,
01033 1.58489,
01034 1.99526,
01035 2.51189,
01036 3.16228,
01037 3.98107,
01038 5.01187,
01039 6.30957,
01040 7.94328,
01041 10,
01042 12.5893,
01043 15.8489,
01044 19.9526,
01045 25.1189,
01046 31.6228,
01047 39.8107,
01048 50.1187,
01049 63.0957,
01050 79.4328,
01051 100,
01052 125.893,
01053 158.489,
01054 199.526,
01055 251.189,
01056 316.228,
01057 398.107,
01058 501.187,
01059 630.957,
01060 794.328,
01061 1000,
01062 1258.93,
01063 1584.89,
01064 1995.26,
01065 2511.89,
01066 3162.28,
01067 3981.07,
01068 5011.87,
01069 6309.57,
01070 7943.28,
01071 10000,
01072 };
01073
01074 Int_t AngleBins = 12;
01075 Double_t AngleEdges[13] = {
01076 -1,
01077 -0.875,
01078 -0.625,
01079 -0.375,
01080 -0.2,
01081 -0.1,
01082 0,
01083 0.1,
01084 0.2,
01085 0.375,
01086 0.625,
01087 0.875,
01088 1,
01089 };
01090
01091 Double_t tempFluxRatioNumu[50][12] = {
01092 { 0.407901, 0.87419, 0.616542, 0.688266, 0.915618, 1.58532,
01093 1.74734, 0.970917, 0.612829, 0.513663, 0.477047, 0.467546 },
01094 { 0.445901, 0.872082, 0.645515, 0.676216, 0.932318, 1.55963,
01095 1.73012, 0.955304, 0.616437, 0.532747, 0.501567, 0.504905 },
01096 { 0.492986, 0.848229, 0.681116, 0.694671, 0.92421, 1.52528,
01097 1.64609, 0.951425, 0.656118, 0.581295, 0.543958, 0.547079 },
01098 { 0.528652, 0.853431, 0.708676, 0.705465, 0.916687, 1.47216,
01099 1.63378, 0.972574, 0.68485, 0.621828, 0.593854, 0.601002 },
01100 { 0.561735, 0.831392, 0.722444, 0.722006, 0.893503, 1.34281,
01101 1.51554, 0.941985, 0.702086, 0.636388, 0.620098, 0.626743 },
01102 { 0.594873, 0.834399, 0.75474, 0.716433, 0.880281, 1.25631,
01103 1.33919, 0.919446, 0.714011, 0.669663, 0.648387, 0.657226 },
01104 { 0.620277, 0.829408, 0.780855, 0.719463, 0.832236, 1.21858,
01105 1.26482, 0.875715, 0.717299, 0.689813, 0.677408, 0.676059 },
01106 { 0.644834, 0.822422, 0.797919, 0.725557, 0.838948, 1.05533,
01107 1.16951, 0.861823, 0.727926, 0.700102, 0.684015, 0.693801 },
01108 { 0.691114, 0.840625, 0.816643, 0.737783, 0.799762, 1.03566,
01109 1.043, 0.816473, 0.729166, 0.717071, 0.730919, 0.717953 },
01110 { 0.719017, 0.850992, 0.836471, 0.752321, 0.816899, 0.973956,
01111 1.02467, 0.806696, 0.748375, 0.760005, 0.743555, 0.771294 },
01112 { 0.764811, 0.84433, 0.865975, 0.779598, 0.828076, 0.949407,
01113 0.974656, 0.829445, 0.778136, 0.799467, 0.78081, 0.817831 },
01114 { 0.801737, 0.871698, 0.890679, 0.817441, 0.874749, 0.969302,
01115 1.00796, 0.891016, 0.82688, 0.835656, 0.832187, 0.846983 },
01116 { 0.851308, 0.88483, 0.902341, 0.815974, 0.878693, 0.920148,
01117 0.918465, 0.87357, 0.834139, 0.857196, 0.86531, 0.895173 },
01118 { 0.904781, 0.933462, 0.895744, 0.864345, 0.902824, 0.960066,
01119 0.919705, 0.937197, 0.871655, 0.860518, 0.910428, 0.912769 },
01120 { 0.928794, 0.942555, 0.916024, 0.880167, 0.980136, 0.943384,
01121 0.950294, 0.912527, 0.908798, 0.893887, 0.927494, 0.937122 },
01122 { 0.975788, 0.906205, 0.954758, 0.929771, 0.944432, 0.923941,
01123 0.917028, 0.925558, 0.92427, 0.913486, 0.938226, 0.962792 },
01124 { 0.959312, 0.941718, 0.95912, 0.935943, 0.948323, 0.858838,
01125 0.893549, 0.958074, 0.909989, 0.976111, 0.928092, 0.965191 },
01126 { 0.950023, 0.982083, 0.954113, 0.936179, 0.947785, 0.899885,
01127 0.963736, 0.970708, 0.912453, 0.91452, 0.966331, 0.984249 },
01128 { 0.929443, 0.96308, 0.977871, 0.93561, 1.04203, 0.935262,
01129 0.971336, 0.989438, 0.934431, 0.991, 0.982813, 0.950919 },
01130 { 1.02293, 0.960648, 0.953228, 0.946973, 0.942216, 1.02765,
01131 0.945596, 0.935716, 0.941016, 0.928756, 0.922628, 1.06962 },
01132 { 0.975442, 0.94712, 0.963152, 0.909081, 0.961021, 0.944187,
01133 0.944215, 0.959489, 0.912573, 0.96302, 0.946721, 0.976709 },
01134 { 0.954521, 0.951504, 0.93789, 0.920042, 0.968064, 0.949156,
01135 0.94715, 0.969824, 0.919684, 0.937861, 0.951884, 0.954588 },
01136 { 0.980221, 0.980739, 0.933387, 0.895646, 0.947115, 0.944987,
01137 0.944963, 0.946951, 0.895734, 0.936308, 0.980007, 0.979505 },
01138 { 0.966654, 0.958888, 0.955399, 0.89338, 0.960503, 0.949842,
01139 0.947058, 0.961756, 0.893413, 0.953794, 0.958538, 0.966057 },
01140 { 0.967193, 0.935157, 0.943601, 0.925294, 0.945592, 0.954411,
01141 0.95444, 0.944864, 0.925321, 0.942915, 0.933895, 0.968163 },
01142 { 0.97509, 0.9379, 0.936289, 0.895395, 0.949033, 0.952325,
01143 0.953872, 0.947776, 0.900497, 0.936483, 0.937825, 0.974925 },
01144 { 0.974844, 0.949273, 0.951143, 0.888057, 0.950144, 0.930221,
01145 0.928132, 0.948617, 0.89038, 0.953059, 0.949053, 0.975188 },
01146 { 0.968934, 0.935001, 0.942496, 0.906561, 0.94847, 0.935474,
01147 0.932941, 0.949119, 0.904026, 0.939556, 0.936168, 0.968738 },
01148 { 0.973721, 0.9599, 0.945159, 0.890205, 0.966573, 0.940351,
01149 0.936725, 0.969651, 0.890462, 0.947307, 0.960404, 0.97169 },
01150 { 0.982156, 0.935489, 0.932213, 0.896294, 0.950343, 0.957662,
01151 0.95603, 0.948461, 0.897865, 0.931526, 0.936307, 0.982302 },
01152 { 0.977345, 0.970102, 0.963545, 0.896852, 0.964495, 0.942714,
01153 0.944971, 0.962793, 0.903185, 0.965066, 0.967328, 0.982911 },
01154 { 1.00439, 0.957662, 0.950131, 0.917232, 0.966896, 0.965724,
01155 0.96331, 0.971072, 0.918149, 0.948911, 0.962587, 1.00431 },
01156 { 0.967905, 1.00098, 0.967671, 0.89003, 0.971013, 0.961357,
01157 0.952684, 0.969011, 0.891553, 0.967983, 0.99899, 0.970126 },
01158 { 0.979338, 0.951973, 0.963948, 0.913014, 0.979198, 0.958417,
01159 0.969097, 0.987354, 0.915201, 0.969012, 0.951276, 0.983347 },
01160 { 1.04077, 0.942393, 0.958535, 0.900534, 0.985788, 0.968012,
01161 0.965025, 0.987405, 0.902251, 0.955214, 0.94414, 1.03586 },
01162 { 1.02074, 0.977487, 0.990401, 0.907411, 0.996908, 0.98743,
01163 0.988375, 0.986959, 0.89185, 0.99036, 0.978413, 1.02841 },
01164 { 1.02476, 1.01422, 0.95784, 0.898782, 0.994976, 1.01447,
01165 1.00977, 1.00243, 0.896561, 0.945832, 1.02195, 1.01565 },
01166 { 1.05813, 0.998507, 1.02375, 0.925399, 1.04316, 1.02827,
01167 1.02887, 1.03004, 0.914676, 1.00111, 1.01798, 1.05783 },
01168 { 1.05609, 1.00462, 0.992902, 0.916825, 1.01045, 1.01208,
01169 1.00721, 1.02104, 0.899057, 0.975705, 1.0157, 1.06022 },
01170 { 1.08674, 0.994028, 0.982287, 0.918874, 1.03912, 1.02558,
01171 1.04771, 1.06151, 0.920488, 0.982144, 0.955515, 1.07647 },
01172 { 1.04069, 1.0178, 1.02958, 0.981677, 1.06222, 1.05969,
01173 1.08046, 1.0426, 0.95391, 1.03907, 1.01133, 1.06558 },
01174 { 1.0953, 1.04574, 1.02943, 0.979644, 1.09715, 1.06773,
01175 1.07287, 1.11252, 0.988495, 1.04784, 1.03353, 1.13861 },
01176 { 1.05675, 1.02863, 1.03059, 0.95608, 1.04195, 1.04743,
01177 1.03884, 1.10631, 0.981417, 1.01581, 1.0475, 1.10792 },
01178 { 1.14721, 1.06448, 1.0757, 0.985877, 1.09084, 1.09804,
01179 1.09319, 1.09343, 0.994542, 1.03235, 1.03548, 1.02076 },
01180 { 1.12154, 1.08814, 1.05434, 0.958804, 1.06365, 1.05436,
01181 1.06191, 1.05345, 0.975105, 1.04579, 1.06025, 1.00413 },
01182 { 1.0512, 1.02273, 1.03443, 0.9821, 1.14604, 1.11062,
01183 1.19599, 1.10277, 0.970488, 1.01207, 0.982797, 1.07774 },
01184 { 1.09247, 0.983922, 0.995253, 0.997355, 1.11813, 1.06071,
01185 1.10208, 1.09786, 0.967704, 1.02564, 1.0107, 1.05864 },
01186 { 1.09168, 0.931565, 0.963358, 0.969301, 1.10096, 1.06495,
01187 1.06789, 1.13827, 0.952162, 1.01003, 0.890653, 0.957078 },
01188 { 1.09717, 1.03783, 1.11627, 0.987227, 1.26891, 1.15563,
01189 1.10412, 1.3178, 0.859991, 1.04903, 0.988425, 1.05713 },
01190 { 0.999478, 1.03798, 1.0156, 0.893401, 1.05337, 1.16683,
01191 1.06669, 1.07133, 0.954657, 0.874146, 1.14809, 1.23959 }
01192 };
01193
01194 fBartolFluxRatioNumu = new TH2D("hBartolFluxRatioNumu","",EnergyBins,EnergyEdges,AngleBins,AngleEdges);
01195 for( Int_t nx=0; nx<EnergyBins; nx++ ){
01196 for( Int_t ny=0; ny<AngleBins; ny++ ){
01197 fBartolFluxRatioNumu->SetBinContent(nx+1,ny+1,tempFluxRatioNumu[nx][ny]);
01198 }
01199 }
01200
01201 Double_t tempFluxRatioNumuBar[50][12] = {
01202 { 0.419909, 0.86696, 0.617805, 0.694347, 0.936305, 1.65452,
01203 1.80188, 0.977291, 0.621024, 0.521226, 0.486601, 0.478059 },
01204 { 0.457087, 0.860854, 0.640657, 0.685283, 0.904763, 1.5367,
01205 1.67528, 0.954441, 0.62423, 0.539869, 0.514669, 0.515865 },
01206 { 0.493547, 0.854435, 0.681477, 0.696709, 0.898733, 1.47652,
01207 1.6406, 0.964271, 0.647458, 0.574885, 0.550217, 0.549446 },
01208 { 0.540152, 0.841968, 0.706792, 0.702489, 0.906074, 1.45667,
01209 1.51981, 0.957375, 0.664831, 0.609647, 0.586539, 0.594554 },
01210 { 0.558002, 0.842949, 0.736069, 0.709655, 0.877461, 1.38931,
01211 1.50214, 0.925314, 0.686758, 0.639664, 0.618067, 0.612849 },
01212 { 0.599399, 0.834899, 0.746364, 0.71837, 0.844918, 1.30243,
01213 1.38521, 0.866158, 0.694154, 0.642933, 0.641561, 0.658374 },
01214 { 0.618166, 0.832418, 0.776161, 0.70901, 0.830014, 1.16974,
01215 1.23157, 0.845438, 0.702651, 0.686819, 0.670706, 0.682772 },
01216 { 0.655249, 0.822341, 0.786444, 0.717981, 0.805452, 1.1002,
01217 1.19497, 0.849365, 0.703929, 0.702722, 0.690835, 0.722096 },
01218 { 0.686463, 0.83479, 0.809993, 0.734854, 0.807546, 1.00753,
01219 1.09195, 0.852882, 0.727053, 0.718116, 0.719951, 0.738093 },
01220 { 0.698685, 0.841827, 0.832503, 0.755738, 0.793793, 1.00293,
01221 1.06289, 0.794962, 0.756311, 0.750372, 0.74752, 0.746921 },
01222 { 0.766557, 0.853241, 0.855347, 0.776714, 0.83853, 0.947186,
01223 0.968055, 0.809326, 0.762216, 0.77733, 0.786487, 0.819667 },
01224 { 0.79006, 0.842004, 0.857574, 0.776131, 0.836566, 0.946153,
01225 0.961994, 0.839372, 0.808635, 0.800766, 0.804363, 0.834603 },
01226 { 0.867501, 0.907733, 0.880754, 0.804795, 0.848948, 0.973121,
01227 0.921603, 0.834636, 0.808941, 0.842156, 0.872722, 0.86865 },
01228 { 0.835037, 0.907365, 0.937813, 0.859058, 0.87842, 0.966943,
01229 0.967174, 0.895318, 0.848295, 0.893034, 0.874765, 0.863815 },
01230 { 0.904532, 0.899746, 0.912705, 0.908259, 0.898149, 0.995856,
01231 1.02335, 0.94406, 0.886822, 0.888308, 0.887971, 0.928682 },
01232 { 0.902963, 0.890525, 0.917504, 0.88469, 0.924986, 0.876871,
01233 0.895394, 0.917661, 0.886617, 0.898717, 0.899884, 0.929002 },
01234 { 0.974293, 0.894307, 0.913058, 0.866757, 0.933625, 0.867731,
01235 0.893765, 0.937024, 0.897496, 0.9144, 0.907913, 0.961545 },
01236 { 0.933304, 0.913667, 0.875373, 0.875868, 0.943681, 0.849525,
01237 0.940958, 0.933712, 0.916019, 0.899732, 0.923929, 0.930549 },
01238 { 0.887209, 0.905562, 0.933587, 0.878906, 0.936487, 0.887606,
01239 0.851895, 0.867064, 0.87355, 0.931544, 0.949696, 0.921185 },
01240 { 0.958679, 0.879565, 0.930212, 0.904138, 0.923778, 0.903275,
01241 0.864285, 0.950362, 0.927507, 0.955342, 0.884841, 0.946768 },
01242 { 0.9197, 0.910948, 0.915674, 0.86876, 0.930058, 0.919897,
01243 0.92128, 0.929724, 0.868734, 0.91318, 0.911031, 0.919353 },
01244 { 0.90697, 0.918557, 0.911189, 0.871546, 0.919664, 0.927289,
01245 0.926447, 0.922136, 0.87174, 0.911396, 0.918127, 0.907802 },
01246 { 0.878294, 0.875183, 0.875293, 0.839296, 0.916957, 0.916563,
01247 0.916527, 0.915745, 0.838976, 0.87636, 0.874607, 0.880959 },
01248 { 0.895263, 0.853895, 0.858013, 0.835479, 0.927295, 0.910455,
01249 0.906606, 0.925536, 0.837389, 0.859093, 0.855112, 0.894417 },
01250 { 0.890479, 0.885531, 0.869181, 0.833928, 0.902262, 0.902183,
01251 0.903061, 0.901446, 0.83324, 0.86894, 0.885738, 0.88942 },
01252 { 0.86012, 0.846431, 0.860264, 0.820006, 0.89339, 0.895805,
01253 0.892723, 0.895727, 0.820076, 0.859227, 0.844828, 0.859259 },
01254 { 0.864054, 0.817477, 0.852692, 0.833401, 0.883778, 0.895342,
01255 0.895162, 0.888089, 0.833661, 0.85345, 0.820111, 0.8614 },
01256 { 0.830286, 0.822156, 0.819251, 0.799991, 0.881222, 0.901546,
01257 0.904078, 0.882164, 0.801415, 0.822587, 0.820294, 0.827235 },
01258 { 0.869974, 0.817073, 0.811808, 0.808597, 0.883987, 0.884445,
01259 0.882085, 0.882358, 0.807124, 0.811401, 0.81458, 0.864946 },
01260 { 0.807137, 0.80434, 0.831008, 0.797565, 0.878628, 0.868712,
01261 0.869756, 0.880621, 0.801195, 0.835496, 0.808747, 0.808106 },
01262 { 0.817639, 0.785971, 0.801099, 0.775249, 0.861652, 0.880509,
01263 0.883827, 0.859116, 0.777157, 0.803615, 0.784214, 0.819032 },
01264 { 0.793224, 0.791831, 0.778408, 0.770102, 0.864941, 0.860756,
01265 0.868271, 0.867974, 0.768037, 0.780019, 0.791509, 0.788658 },
01266 { 0.789749, 0.762709, 0.791717, 0.760806, 0.854404, 0.872839,
01267 0.871276, 0.856818, 0.759316, 0.791583, 0.759794, 0.789373 },
01268 { 0.777213, 0.755061, 0.768007, 0.776599, 0.87385, 0.887033,
01269 0.887403, 0.871721, 0.768832, 0.76777, 0.758978, 0.774709 },
01270 { 0.782047, 0.717187, 0.755737, 0.744113, 0.871312, 0.882519,
01271 0.881111, 0.867596, 0.739795, 0.760053, 0.72427, 0.794074 },
01272 { 0.750635, 0.731569, 0.758517, 0.693996, 0.856942, 0.898877,
01273 0.900438, 0.851152, 0.694862, 0.75975, 0.736655, 0.763372 },
01274 { 0.761118, 0.732715, 0.772226, 0.688345, 0.847226, 0.873483,
01275 0.873986, 0.851444, 0.697681, 0.770092, 0.728985, 0.769174 },
01276 { 0.774607, 0.73181, 0.736005, 0.680346, 0.838504, 0.877505,
01277 0.878862, 0.835941, 0.672439, 0.733923, 0.739178, 0.763102 },
01278 { 0.742569, 0.703755, 0.715088, 0.664591, 0.830326, 0.874947,
01279 0.880915, 0.833094, 0.674989, 0.716294, 0.708699, 0.738579 },
01280 { 0.737158, 0.689996, 0.730289, 0.701398, 0.836754, 0.881603,
01281 0.886054, 0.844662, 0.708866, 0.724037, 0.690787, 0.738451 },
01282 { 0.769987, 0.711012, 0.74531, 0.710789, 0.823083, 0.840787,
01283 0.836304, 0.808084, 0.712756, 0.722211, 0.714257, 0.745536 },
01284 { 0.751414, 0.71811, 0.718008, 0.707951, 0.811322, 0.874588,
01285 0.843233, 0.820841, 0.696452, 0.70208, 0.691169, 0.712002 },
01286 { 0.736834, 0.701795, 0.695325, 0.689812, 0.815717, 0.85955,
01287 0.844726, 0.780836, 0.689945, 0.705017, 0.680911, 0.74261 },
01288 { 0.701682, 0.675908, 0.682918, 0.657468, 0.773127, 0.830455,
01289 0.833518, 0.772601, 0.672326, 0.680118, 0.738788, 0.725297 },
01290 { 0.735046, 0.691412, 0.687795, 0.691496, 0.868924, 0.800936,
01291 0.815144, 0.848892, 0.678939, 0.69569, 0.700727, 0.730359 },
01292 { 0.749089, 0.68572, 0.673398, 0.656938, 0.867693, 0.745819,
01293 0.766935, 0.856235, 0.64287, 0.671985, 0.681261, 0.715962 },
01294 { 0.686553, 0.737141, 0.705175, 0.640149, 0.808821, 0.859627,
01295 0.835334, 0.827694, 0.654643, 0.694448, 0.641269, 0.75019 },
01296 { 0.833357, 0.668819, 0.704183, 0.617385, 0.78526, 0.911363,
01297 0.788458, 0.819574, 0.620585, 0.675242, 0.73846, 0.733908 },
01298 { 0.821739, 0.658415, 0.719287, 0.686069, 0.831316, 0.882383,
01299 0.799232, 0.799462, 0.632223, 0.686581, 0.771696, 0.80366 },
01300 { 0.786634, 0.66479, 0.750568, 0.718219, 0.813377, 0.879924,
01301 0.830972, 0.885847, 0.657697, 0.708878, 0.658219, 0.946438 }
01302 };
01303
01304 fBartolFluxRatioNumuBar = new TH2D("hBartolFluxRatioNumuBar","",EnergyBins,EnergyEdges,AngleBins,AngleEdges);
01305 for( Int_t nx=0; nx<EnergyBins; nx++ ){
01306 for( Int_t ny=0; ny<AngleBins; ny++ ){
01307 fBartolFluxRatioNumuBar->SetBinContent(nx+1,ny+1,tempFluxRatioNumuBar[nx][ny]);
01308 }
01309 }
01310
01311 return;
01312 }
01313
01314 void AtNuFluxRW::LoadBartolFlavourRatios()
01315 {
01316
01317
01318
01319
01320
01321 Int_t EnergyBins = 50;
01322 Double_t EnergyEdges[51] = {
01323 0.1,
01324 0.125893,
01325 0.158489,
01326 0.199526,
01327 0.251189,
01328 0.316228,
01329 0.398107,
01330 0.501187,
01331 0.630957,
01332 0.794328,
01333 1,
01334 1.25893,
01335 1.58489,
01336 1.99526,
01337 2.51189,
01338 3.16228,
01339 3.98107,
01340 5.01187,
01341 6.30957,
01342 7.94328,
01343 10,
01344 12.5893,
01345 15.8489,
01346 19.9526,
01347 25.1189,
01348 31.6228,
01349 39.8107,
01350 50.1187,
01351 63.0957,
01352 79.4328,
01353 100,
01354 125.893,
01355 158.489,
01356 199.526,
01357 251.189,
01358 316.228,
01359 398.107,
01360 501.187,
01361 630.957,
01362 794.328,
01363 1000,
01364 1258.93,
01365 1584.89,
01366 1995.26,
01367 2511.89,
01368 3162.28,
01369 3981.07,
01370 5011.87,
01371 6309.57,
01372 7943.28,
01373 10000,
01374 };
01375
01376 Int_t AngleBins = 20;
01377 Double_t AngleEdges[21] = {
01378 -1,
01379 -0.9,
01380 -0.8,
01381 -0.7,
01382 -0.6,
01383 -0.5,
01384 -0.4,
01385 -0.3,
01386 -0.2,
01387 -0.1,
01388 0,
01389 0.1,
01390 0.2,
01391 0.3,
01392 0.4,
01393 0.5,
01394 0.6,
01395 0.7,
01396 0.8,
01397 0.9,
01398 1,
01399 };
01400
01401 Double_t tempFlavourRatioNumu[50][20] = {
01402 { 0.483771, 0.482108, 0.468501, 0.467736, 0.487089,
01403 0.490678, 0.485871, 0.478179, 0.488448, 0.495532,
01404 0.487556, 0.493933, 0.497356, 0.490054, 0.49226,
01405 0.499126, 0.486597, 0.493977, 0.492738, 0.488128 },
01406 { 0.493145, 0.486648, 0.466921, 0.485041, 0.50345,
01407 0.49147, 0.517959, 0.515119, 0.507591, 0.536394,
01408 0.499488, 0.519409, 0.532653, 0.515699, 0.520559,
01409 0.51232, 0.528917, 0.511838, 0.505854, 0.503749 },
01410 { 0.50345, 0.488795, 0.47627, 0.489905, 0.49574,
01411 0.510479, 0.525588, 0.529513, 0.522196, 0.54615,
01412 0.558027, 0.533383, 0.530281, 0.526716, 0.531178,
01413 0.5302, 0.529989, 0.526617, 0.526447, 0.511257 },
01414 { 0.494412, 0.476955, 0.470938, 0.481523, 0.50026,
01415 0.500957, 0.50404, 0.522552, 0.51449, 0.514808,
01416 0.546352, 0.522192, 0.531277, 0.527599, 0.524156,
01417 0.533451, 0.522281, 0.520109, 0.510244, 0.497974 },
01418 { 0.485705, 0.480266, 0.477473, 0.484447, 0.50488,
01419 0.517759, 0.515405, 0.518466, 0.554715, 0.551937,
01420 0.543441, 0.534474, 0.548143, 0.534233, 0.538937,
01421 0.526618, 0.522721, 0.521193, 0.501988, 0.499907 },
01422 { 0.49551, 0.487358, 0.474522, 0.478369, 0.505141,
01423 0.524233, 0.540578, 0.547802, 0.549412, 0.572686,
01424 0.6097, 0.552332, 0.552638, 0.569476, 0.552663,
01425 0.533432, 0.539793, 0.518334, 0.524457, 0.502119 },
01426 { 0.485826, 0.481618, 0.467989, 0.48277, 0.496884,
01427 0.526026, 0.533962, 0.561812, 0.573883, 0.589635,
01428 0.599014, 0.593377, 0.574603, 0.562757, 0.565321,
01429 0.544403, 0.533308, 0.525284, 0.508585, 0.492155 },
01430 { 0.500213, 0.476787, 0.461245, 0.4819, 0.516985,
01431 0.520957, 0.537474, 0.576207, 0.581938, 0.591559,
01432 0.565733, 0.578725, 0.56797, 0.545511, 0.561957,
01433 0.551516, 0.523841, 0.51313, 0.502237, 0.495146 },
01434 { 0.45806, 0.464809, 0.452513, 0.479841, 0.491237,
01435 0.503463, 0.551025, 0.558151, 0.593469, 0.559609,
01436 0.590319, 0.588835, 0.582197, 0.565318, 0.54909,
01437 0.523059, 0.528094, 0.498051, 0.479855, 0.484975 },
01438 { 0.450572, 0.450117, 0.458697, 0.474582, 0.493752,
01439 0.534661, 0.560879, 0.567347, 0.593331, 0.644814,
01440 0.623036, 0.576097, 0.576903, 0.573431, 0.55282,
01441 0.52368, 0.51218, 0.493418, 0.49566, 0.454642 },
01442 { 0.428516, 0.433975, 0.453646, 0.457713, 0.490208,
01443 0.535391, 0.551684, 0.566821, 0.599828, 0.617741,
01444 0.62229, 0.5984, 0.592081, 0.552577, 0.546255,
01445 0.521942, 0.499915, 0.466678, 0.471527, 0.433905 },
01446 { 0.414213, 0.408831, 0.428097, 0.469648, 0.472421,
01447 0.512722, 0.53499, 0.594557, 0.591875, 0.605494,
01448 0.618075, 0.60744, 0.580938, 0.560856, 0.529507,
01449 0.511561, 0.488459, 0.447543, 0.431825, 0.417872 },
01450 { 0.387498, 0.395137, 0.422848, 0.439331, 0.469314,
01451 0.497287, 0.557325, 0.585158, 0.601038, 0.688688,
01452 0.638093, 0.587215, 0.568578, 0.55804, 0.525674,
01453 0.495038, 0.458604, 0.448664, 0.420186, 0.411906 },
01454 { 0.362212, 0.387219, 0.399575, 0.41824, 0.451681,
01455 0.468658, 0.538262, 0.550843, 0.602911, 0.606472,
01456 0.677418, 0.580113, 0.578858, 0.543293, 0.50707,
01457 0.465367, 0.435033, 0.426077, 0.403338, 0.381792 },
01458 { 0.335272, 0.335697, 0.380492, 0.366469, 0.433521,
01459 0.468824, 0.500071, 0.54058, 0.524332, 0.599026,
01460 0.598978, 0.573493, 0.553698, 0.49643, 0.470337,
01461 0.440934, 0.392765, 0.401862, 0.346422, 0.33904 },
01462 { 0.300426, 0.31122, 0.347986, 0.366456, 0.409579,
01463 0.475004, 0.473914, 0.522711, 0.605484, 0.595353,
01464 0.626054, 0.581413, 0.503464, 0.488057, 0.471954,
01465 0.397811, 0.363316, 0.371814, 0.315188, 0.303112 },
01466 { 0.281652, 0.307674, 0.293763, 0.320416, 0.359096,
01467 0.439461, 0.438454, 0.494095, 0.579221, 0.649487,
01468 0.632328, 0.559321, 0.521527, 0.456911, 0.375399,
01469 0.371075, 0.346508, 0.336367, 0.316047, 0.288439 },
01470 { 0.228568, 0.253805, 0.271087, 0.29186, 0.329319,
01471 0.407703, 0.411685, 0.467095, 0.552315, 0.692302,
01472 0.625927, 0.502497, 0.502373, 0.424052, 0.395593,
01473 0.364602, 0.322787, 0.297383, 0.276654, 0.240011 },
01474 { 0.246336, 0.249708, 0.235099, 0.295655, 0.281298,
01475 0.332034, 0.399846, 0.43911, 0.482904, 0.606731,
01476 0.538001, 0.472843, 0.442162, 0.390709, 0.340242,
01477 0.3173, 0.308261, 0.25273, 0.248806, 0.230625 },
01478 { 0.1631, 0.207982, 0.205912, 0.242623, 0.303877,
01479 0.288542, 0.338373, 0.437408, 0.521054, 0.597131,
01480 0.642662, 0.497307, 0.449432, 0.363256, 0.291221,
01481 0.321655, 0.234413, 0.240831, 0.229599, 0.207458 },
01482 { 0.169401, 0.184324, 0.200897, 0.222153, 0.249569,
01483 0.282551, 0.328894, 0.385506, 0.460316, 0.534085,
01484 0.534748, 0.462801, 0.385854, 0.327066, 0.282911,
01485 0.247291, 0.221632, 0.199857, 0.185385, 0.168774 },
01486 { 0.150167, 0.161729, 0.17881, 0.196742, 0.220957,
01487 0.252632, 0.293562, 0.352919, 0.431064, 0.512889,
01488 0.514071, 0.429936, 0.35205, 0.295427, 0.250656,
01489 0.219221, 0.195592, 0.177684, 0.161431, 0.150135 },
01490 { 0.132826, 0.144069, 0.155376, 0.174172, 0.194853,
01491 0.223914, 0.263341, 0.318823, 0.398018, 0.486413,
01492 0.490415, 0.396536, 0.318519, 0.261074, 0.223107,
01493 0.194106, 0.174984, 0.15696, 0.144289, 0.134147 },
01494 { 0.119721, 0.129023, 0.139478, 0.154179, 0.171923,
01495 0.19694, 0.234137, 0.285904, 0.363016, 0.458716,
01496 0.460095, 0.363648, 0.285564, 0.232981, 0.197246,
01497 0.173353, 0.154122, 0.14036, 0.127371, 0.118951 },
01498 { 0.108839, 0.113704, 0.123831, 0.137892, 0.152174,
01499 0.17414, 0.208053, 0.255379, 0.331325, 0.430478,
01500 0.431342, 0.329597, 0.253939, 0.206669, 0.17496,
01501 0.152544, 0.136938, 0.123935, 0.11439, 0.108862 },
01502 { 0.0973561, 0.103436, 0.111873, 0.121842, 0.134978,
01503 0.154286, 0.182153, 0.226978, 0.296657, 0.400109,
01504 0.400045, 0.297045, 0.224974, 0.180749, 0.153104,
01505 0.135703, 0.12208, 0.111167, 0.104326, 0.0995421 },
01506 { 0.0899228, 0.0953126, 0.101151, 0.110354, 0.121576,
01507 0.137762, 0.160898, 0.19966, 0.266258, 0.369011,
01508 0.369623, 0.266974, 0.199365, 0.161833, 0.136238,
01509 0.121946, 0.109746, 0.101529, 0.0953669, 0.091064 },
01510 { 0.0845741, 0.0881361, 0.0937866, 0.100352, 0.108946,
01511 0.123213, 0.143108, 0.176464, 0.237217, 0.337551,
01512 0.338259, 0.236706, 0.177374, 0.143615, 0.122609,
01513 0.109579, 0.0998899, 0.0943006, 0.0888237, 0.0846918 },
01514 { 0.0785875, 0.0810312, 0.0883395, 0.0917064, 0.100236,
01515 0.111409, 0.129469, 0.1586, 0.21045, 0.309442,
01516 0.310171, 0.208367, 0.155461, 0.129071, 0.109489,
01517 0.100451, 0.0940174, 0.0864788, 0.0830078, 0.0785639 },
01518 { 0.0744995, 0.0780492, 0.0813225, 0.0873548, 0.0934551,
01519 0.101645, 0.115324, 0.137047, 0.187218, 0.280671,
01520 0.279426, 0.186033, 0.139027, 0.115385, 0.0993897,
01521 0.0921002, 0.0859814, 0.0797118, 0.07803, 0.0766232 },
01522 { 0.0727529, 0.0733828, 0.0763106, 0.0804042, 0.0875126,
01523 0.0928214, 0.10437, 0.124952, 0.163905, 0.255326,
01524 0.254744, 0.165525, 0.12291, 0.104189, 0.0909762,
01525 0.085886, 0.0830063, 0.0789492, 0.0763853, 0.071196 },
01526 { 0.0669383, 0.0742784, 0.0745275, 0.076596, 0.0791858,
01527 0.0841418, 0.0954737, 0.111703, 0.148699, 0.229631,
01528 0.232904, 0.145686, 0.112031, 0.0978979, 0.0866845,
01529 0.0800304, 0.0759834, 0.0719356, 0.0731054, 0.0695679 },
01530 { 0.065199, 0.0675056, 0.0724757, 0.0699149, 0.0762228,
01531 0.079739, 0.0894119, 0.103028, 0.133076, 0.207396,
01532 0.212309, 0.129558, 0.10342, 0.0890477, 0.0818844,
01533 0.0775523, 0.0766142, 0.0696595, 0.0669001, 0.0666938 },
01534 { 0.064123, 0.0696213, 0.0655516, 0.0681806, 0.0738569,
01535 0.0781168, 0.0826077, 0.0945174, 0.114229, 0.194846,
01536 0.189752, 0.118109, 0.0922713, 0.083229, 0.0775152,
01537 0.0723676, 0.0686127, 0.0670341, 0.0640538, 0.0624709 },
01538 { 0.0591662, 0.0617666, 0.0615119, 0.0625595, 0.0693798,
01539 0.070784, 0.0763717, 0.0857307, 0.105883, 0.173673,
01540 0.174281, 0.10342, 0.0882207, 0.0774242, 0.0723012,
01541 0.0681079, 0.0638167, 0.0624689, 0.0604721, 0.0609709 },
01542 { 0.0558325, 0.057686, 0.0588582, 0.0665832, 0.0666299,
01543 0.0685581, 0.0737307, 0.0764614, 0.0964915, 0.16324,
01544 0.159981, 0.0995177, 0.0820128, 0.0743402, 0.0655594,
01545 0.0677487, 0.0670363, 0.0641988, 0.0596194, 0.0597144 },
01546 { 0.0587265, 0.0583574, 0.0607261, 0.0672015, 0.0623116,
01547 0.0654502, 0.0648915, 0.0760605, 0.0939325, 0.144998,
01548 0.143789, 0.088923, 0.0783105, 0.066246, 0.0710467,
01549 0.0619893, 0.0600036, 0.0534692, 0.0569893, 0.0561573 },
01550 { 0.0547225, 0.0596272, 0.0556315, 0.0568475, 0.0582189,
01551 0.0596421, 0.0703799, 0.0723212, 0.0803562, 0.129116,
01552 0.131531, 0.0858948, 0.0701371, 0.0620491, 0.0645727,
01553 0.0578932, 0.0585832, 0.0598886, 0.0559077, 0.0550951 },
01554 { 0.0552642, 0.0551161, 0.053484, 0.0591354, 0.0591139,
01555 0.059932, 0.0655763, 0.0678727, 0.0825035, 0.120444,
01556 0.118474, 0.0800536, 0.0659849, 0.0669378, 0.0606173,
01557 0.058444, 0.0576221, 0.0528608, 0.0556769, 0.0539962 },
01558 { 0.0489014, 0.0498292, 0.0475769, 0.0534821, 0.0549172,
01559 0.0559145, 0.0594824, 0.0625635, 0.0704311, 0.107317,
01560 0.109912, 0.0691772, 0.0630019, 0.0573904, 0.0587197,
01561 0.0556742, 0.0597354, 0.0520465, 0.0506518, 0.0496374 },
01562 { 0.0461725, 0.0528534, 0.0555651, 0.0557047, 0.0520061,
01563 0.060214, 0.0477867, 0.0589459, 0.0727363, 0.0973311,
01564 0.100299, 0.0693673, 0.0626349, 0.0622551, 0.0511697,
01565 0.0478783, 0.0493531, 0.0501689, 0.046332, 0.0582558 },
01566 { 0.0551047, 0.0498504, 0.0567048, 0.045668, 0.0483489,
01567 0.049971, 0.0519016, 0.0608866, 0.0646732, 0.0869788,
01568 0.0993125, 0.0616546, 0.0548603, 0.0621514, 0.0532224,
01569 0.0500748, 0.0494954, 0.0488039, 0.0493598, 0.0464522 },
01570 { 0.0578186, 0.0491193, 0.0521767, 0.0522645, 0.0500484,
01571 0.0458721, 0.062098, 0.0493422, 0.0583262, 0.0843645,
01572 0.0793672, 0.0646252, 0.0600312, 0.0530195, 0.051245,
01573 0.0598124, 0.0521268, 0.0508563, 0.0496217, 0.0590737 },
01574 { 0.0546895, 0.0470785, 0.0492463, 0.0401263, 0.0393656,
01575 0.0574015, 0.0425785, 0.0581796, 0.061947, 0.0794304,
01576 0.087147, 0.0586097, 0.0555622, 0.0532719, 0.050666,
01577 0.0613289, 0.0468641, 0.0475148, 0.0397562, 0.0490451 },
01578 { 0.0479479, 0.0525138, 0.0429685, 0.0499666, 0.0432402,
01579 0.0467873, 0.0499847, 0.050154, 0.0664756, 0.068267,
01580 0.0582761, 0.0602447, 0.051397, 0.0445861, 0.0556882,
01581 0.0453375, 0.0450909, 0.0404872, 0.057043, 0.0631275 },
01582 { 0.029464, 0.0388906, 0.0500278, 0.0442822, 0.0396983,
01583 0.0320882, 0.0536825, 0.0407348, 0.0453387, 0.0661033,
01584 0.0665525, 0.0531895, 0.0453006, 0.0705818, 0.0526684,
01585 0.0527037, 0.0451247, 0.0524748, 0.0545379, 0.070268 },
01586 { 0.0551849, 0.0458357, 0.0453332, 0.0490514, 0.0481249,
01587 0.052368, 0.0423143, 0.0575992, 0.0634309, 0.0660974,
01588 0.057854, 0.0495832, 0.0527541, 0.0430467, 0.0392717,
01589 0.0502788, 0.0372179, 0.0250528, 0.029136, 0.0571204 },
01590 { 0.0578826, 0.0250609, 0.0365661, 0.0539183, 0.0531626,
01591 0.085322, 0.057191, 0.0291048, 0.0593329, 0.0592411,
01592 0.0792843, 0.0364724, 0.0579378, 0.0600581, 0.0386804,
01593 0.0158845, 0.0504172, 0.0436397, 0.0410874, 0.0747608 },
01594 { 0.047278, 0.0108266, 0.0581433, 0.0435173, 0.0282802,
01595 0.0556926, 0.0608881, 0.0341592, 0.0468512, 0.0709938,
01596 0.0744972, 0.0344647, 0.0514944, 0.0616371, 0.0306792,
01597 0.0458638, 0.0402185, 0.0402181, 0.0544692, 0.0563191 },
01598 { 0.0488874, 0.016854, 0.0287852, 0.0486452, 0.0284928,
01599 0.0474454, 0.0162495, 0.0377561, 0.0634744, 0.0407885,
01600 0.060515, 0.0521133, 0.0392173, 0.0298258, 0.0295891,
01601 0.0556119, 0.0537151, 0.0231033, 0.00995481, 0.0147139 }
01602 };
01603
01604 fBartolFlavourRatioNumu = new TH2D("hBartolFlavourRatioNumu","",EnergyBins,EnergyEdges,AngleBins,AngleEdges);
01605 for( Int_t nx=0; nx<EnergyBins; nx++ ){
01606 for( Int_t ny=0; ny<AngleBins; ny++ ){
01607 fBartolFlavourRatioNumu->SetBinContent(nx+1,ny+1,tempFlavourRatioNumu[nx][ny]);
01608 }
01609 }
01610
01611 Double_t tempFlavourRatioNumuBar[50][20] = {
01612 { 0.412082, 0.446207, 0.465096, 0.460359, 0.453281,
01613 0.437887, 0.435633, 0.436573, 0.443728, 0.431594,
01614 0.42742, 0.424616, 0.426706, 0.422855, 0.426578,
01615 0.422978, 0.409966, 0.415512, 0.420574, 0.40902 },
01616 { 0.429984, 0.450796, 0.467453, 0.46128, 0.460112,
01617 0.471711, 0.456325, 0.461657, 0.443761, 0.46554,
01618 0.454107, 0.458411, 0.457236, 0.450232, 0.430269,
01619 0.443141, 0.428164, 0.435857, 0.428362, 0.424918 },
01620 { 0.433105, 0.446292, 0.466956, 0.472778, 0.459664,
01621 0.461006, 0.458654, 0.454493, 0.474067, 0.450731,
01622 0.432456, 0.443847, 0.464821, 0.443334, 0.451395,
01623 0.438165, 0.430345, 0.437035, 0.424072, 0.42411 },
01624 { 0.41547, 0.44927, 0.462523, 0.455846, 0.464306,
01625 0.463399, 0.437637, 0.452316, 0.449527, 0.443117,
01626 0.457429, 0.448325, 0.455673, 0.45393, 0.435496,
01627 0.441375, 0.428283, 0.427373, 0.419763, 0.422682 },
01628 { 0.423694, 0.436987, 0.448645, 0.458416, 0.45381,
01629 0.443066, 0.458565, 0.456408, 0.452544, 0.454817,
01630 0.453412, 0.452277, 0.451623, 0.449577, 0.444839,
01631 0.432767, 0.436776, 0.427372, 0.418028, 0.416475 },
01632 { 0.415695, 0.429557, 0.435888, 0.460423, 0.453169,
01633 0.453298, 0.452958, 0.461324, 0.474823, 0.481373,
01634 0.44374, 0.470241, 0.460486, 0.455262, 0.446434,
01635 0.432897, 0.431226, 0.425379, 0.423152, 0.414652 },
01636 { 0.414974, 0.41543, 0.433437, 0.451531, 0.442472,
01637 0.436401, 0.457848, 0.46187, 0.463416, 0.467512,
01638 0.468489, 0.472127, 0.468293, 0.448049, 0.440055,
01639 0.436197, 0.439518, 0.424458, 0.426239, 0.410102 },
01640 { 0.404801, 0.416695, 0.426806, 0.449237, 0.451902,
01641 0.449419, 0.458291, 0.453783, 0.475248, 0.461036,
01642 0.450015, 0.474451, 0.467286, 0.452803, 0.443891,
01643 0.439374, 0.429032, 0.416778, 0.415954, 0.3826 },
01644 { 0.386581, 0.417341, 0.419535, 0.447522, 0.451983,
01645 0.442891, 0.45073, 0.462137, 0.473213, 0.506575,
01646 0.472663, 0.470637, 0.456902, 0.452977, 0.448125,
01647 0.432245, 0.425759, 0.421255, 0.405408, 0.390853 },
01648 { 0.38914, 0.396454, 0.408412, 0.420317, 0.434366,
01649 0.448955, 0.469308, 0.463137, 0.506365, 0.473309,
01650 0.490495, 0.483602, 0.466677, 0.458701, 0.441111,
01651 0.441069, 0.416925, 0.401729, 0.392452, 0.386146 },
01652 { 0.364399, 0.376711, 0.395371, 0.415891, 0.414879,
01653 0.426647, 0.451585, 0.458355, 0.449111, 0.478984,
01654 0.480986, 0.492351, 0.469532, 0.44699, 0.445114,
01655 0.413849, 0.397488, 0.38728, 0.377112, 0.36034 },
01656 { 0.353007, 0.373926, 0.38936, 0.3944, 0.40981,
01657 0.43751, 0.456229, 0.461252, 0.491707, 0.451962,
01658 0.46186, 0.457081, 0.447181, 0.413908, 0.418922,
01659 0.420504, 0.390046, 0.380007, 0.365469, 0.348283 },
01660 { 0.303415, 0.363492, 0.358637, 0.376342, 0.4085,
01661 0.414435, 0.422241, 0.46363, 0.486646, 0.462166,
01662 0.448602, 0.494764, 0.487072, 0.447774, 0.415816,
01663 0.401511, 0.398458, 0.358477, 0.345809, 0.333924 },
01664 { 0.341528, 0.309434, 0.333661, 0.365176, 0.39094,
01665 0.407794, 0.420635, 0.419859, 0.470557, 0.470652,
01666 0.481248, 0.476548, 0.445981, 0.417633, 0.378496,
01667 0.388377, 0.362502, 0.355671, 0.310379, 0.305583 },
01668 { 0.29191, 0.314489, 0.318564, 0.350551, 0.336322,
01669 0.379333, 0.400711, 0.429679, 0.437953, 0.455759,
01670 0.401319, 0.454079, 0.448966, 0.405334, 0.374829,
01671 0.3625, 0.35585, 0.30944, 0.299335, 0.274751 },
01672 { 0.293613, 0.294652, 0.311178, 0.322521, 0.340589,
01673 0.367259, 0.379857, 0.423452, 0.466335, 0.470098,
01674 0.43711, 0.46497, 0.403181, 0.396796, 0.365467,
01675 0.355684, 0.321824, 0.288057, 0.275021, 0.281563 },
01676 { 0.239426, 0.273545, 0.313211, 0.308953, 0.315088,
01677 0.373761, 0.414211, 0.420854, 0.464901, 0.481912,
01678 0.487575, 0.458679, 0.418835, 0.389034, 0.357844,
01679 0.334657, 0.298711, 0.299819, 0.265162, 0.243532 },
01680 { 0.214475, 0.257205, 0.26995, 0.291824, 0.334321,
01681 0.328389, 0.375506, 0.402932, 0.445746, 0.500525,
01682 0.425476, 0.460156, 0.380677, 0.377773, 0.317102,
01683 0.312775, 0.279455, 0.280969, 0.264316, 0.221291 },
01684 { 0.202441, 0.225885, 0.255472, 0.267166, 0.287013,
01685 0.327647, 0.366529, 0.410451, 0.403483, 0.479243,
01686 0.419496, 0.482984, 0.402491, 0.382375, 0.348764,
01687 0.257041, 0.254401, 0.236084, 0.216587, 0.213084 },
01688 { 0.179954, 0.193844, 0.222862, 0.230639, 0.288656,
01689 0.286669, 0.336902, 0.362264, 0.45466, 0.49592,
01690 0.477066, 0.426637, 0.374788, 0.308645, 0.311289,
01691 0.276291, 0.228307, 0.20893, 0.238115, 0.23494 },
01692 { 0.174504, 0.185824, 0.202585, 0.218628, 0.239263,
01693 0.268425, 0.300748, 0.34331, 0.396976, 0.447518,
01694 0.446782, 0.39663, 0.345322, 0.302963, 0.267699,
01695 0.241377, 0.218861, 0.200956, 0.186883, 0.173887 },
01696 { 0.159006, 0.17012, 0.182721, 0.199817, 0.222308,
01697 0.24704, 0.279989, 0.322346, 0.379958, 0.434932,
01698 0.436291, 0.379098, 0.324053, 0.279305, 0.246665,
01699 0.220629, 0.200324, 0.182691, 0.169508, 0.15846 },
01700 { 0.145755, 0.154856, 0.166888, 0.181873, 0.200191,
01701 0.224519, 0.25634, 0.301256, 0.358876, 0.421369,
01702 0.421391, 0.361286, 0.299859, 0.257893, 0.225611,
01703 0.200406, 0.180594, 0.16616, 0.155017, 0.14533 },
01704 { 0.133866, 0.143334, 0.152407, 0.165488, 0.182148,
01705 0.205248, 0.23422, 0.279239, 0.337267, 0.404417,
01706 0.406198, 0.337136, 0.276908, 0.234852, 0.204838,
01707 0.182104, 0.165932, 0.152539, 0.143214, 0.133841 },
01708 { 0.123566, 0.131634, 0.139287, 0.151834, 0.16804,
01709 0.18802, 0.214236, 0.256979, 0.316113, 0.390613,
01710 0.386573, 0.317894, 0.25559, 0.215569, 0.187054,
01711 0.168076, 0.153555, 0.139994, 0.131555, 0.125902 },
01712 { 0.115897, 0.121767, 0.131581, 0.139697, 0.151212,
01713 0.171003, 0.19538, 0.233541, 0.293581, 0.370069,
01714 0.369491, 0.293331, 0.233459, 0.194641, 0.172387,
01715 0.153263, 0.139037, 0.129894, 0.123654, 0.116666 },
01716 { 0.111004, 0.116278, 0.121219, 0.129864, 0.140257,
01717 0.155804, 0.178266, 0.213658, 0.271505, 0.349669,
01718 0.351192, 0.269197, 0.213176, 0.179079, 0.157679,
01719 0.141758, 0.12908, 0.121258, 0.116909, 0.112008 },
01720 { 0.106056, 0.113147, 0.114718, 0.12155, 0.131086,
01721 0.144228, 0.162865, 0.19482, 0.246457, 0.330992,
01722 0.327864, 0.247094, 0.193329, 0.160031, 0.142282,
01723 0.130932, 0.122047, 0.115993, 0.110271, 0.107866 },
01724 { 0.105538, 0.106666, 0.110969, 0.116491, 0.121385,
01725 0.134591, 0.149833, 0.175646, 0.22705, 0.309264,
01726 0.310382, 0.226335, 0.177202, 0.14874, 0.134385,
01727 0.122798, 0.114333, 0.113435, 0.106166, 0.104194 },
01728 { 0.104183, 0.10512, 0.107658, 0.112914, 0.115858,
01729 0.123622, 0.141049, 0.165866, 0.20915, 0.290458,
01730 0.290746, 0.208201, 0.162713, 0.136834, 0.127399,
01731 0.118285, 0.109107, 0.104548, 0.102973, 0.104227 },
01732 { 0.0985275, 0.103386, 0.105375, 0.107694, 0.112493,
01733 0.119188, 0.132633, 0.150316, 0.189637, 0.271861,
01734 0.271673, 0.188251, 0.154303, 0.131769, 0.118001,
01735 0.114698, 0.106518, 0.106363, 0.102147, 0.0984905 },
01736 { 0.0986039, 0.0983164, 0.104228, 0.100105, 0.106353,
01737 0.113724, 0.122595, 0.138818, 0.177965, 0.252492,
01738 0.252904, 0.172283, 0.139934, 0.122001, 0.113189,
01739 0.10818, 0.105292, 0.0996621, 0.0995675, 0.0976611 },
01740 { 0.0952388, 0.0959157, 0.102642, 0.104566, 0.107753,
01741 0.111432, 0.117709, 0.131835, 0.15903, 0.239543,
01742 0.240279, 0.160642, 0.130822, 0.11649, 0.107555,
01743 0.105715, 0.103211, 0.0982097, 0.103424, 0.0970446 },
01744 { 0.0966725, 0.0984501, 0.101105, 0.100576, 0.0991159,
01745 0.109583, 0.11332, 0.12501, 0.147122, 0.224034,
01746 0.220114, 0.150081, 0.124032, 0.115834, 0.105434,
01747 0.103624, 0.105005, 0.101644, 0.0939662, 0.0933233 },
01748 { 0.099124, 0.0934924, 0.0917282, 0.101118, 0.104825,
01749 0.098954, 0.115263, 0.116816, 0.137998, 0.213482,
01750 0.207205, 0.138229, 0.119599, 0.110535, 0.104773,
01751 0.101785, 0.0996848, 0.0961159, 0.0906688, 0.0951277 },
01752 { 0.0922214, 0.0924557, 0.09404, 0.101346, 0.104847,
01753 0.0987602, 0.105733, 0.110381, 0.124784, 0.192419,
01754 0.195219, 0.127529, 0.1085, 0.107653, 0.108114,
01755 0.0984656, 0.101836, 0.0979752, 0.0880972, 0.0926429 },
01756 { 0.0869145, 0.0962932, 0.0990207, 0.0951642, 0.0909166,
01757 0.101027, 0.100303, 0.116116, 0.121938, 0.181358,
01758 0.187503, 0.12106, 0.112041, 0.101635, 0.106341,
01759 0.100281, 0.0836268, 0.0945236, 0.0882399, 0.084966 },
01760 { 0.0833919, 0.0833637, 0.0852624, 0.0929501, 0.0855587,
01761 0.0982595, 0.103193, 0.0989452, 0.120303, 0.171511,
01762 0.169597, 0.116081, 0.100552, 0.0975868, 0.102001,
01763 0.0992304, 0.0995348, 0.0894016, 0.0808943, 0.0871339 },
01764 { 0.0849254, 0.100271, 0.0852303, 0.0851574, 0.0922799,
01765 0.0914357, 0.103621, 0.101908, 0.110508, 0.161002,
01766 0.163773, 0.110799, 0.0911506, 0.101228, 0.0935324,
01767 0.0913888, 0.0829404, 0.0949507, 0.0862389, 0.0890564 },
01768 { 0.0866945, 0.0736328, 0.080485, 0.0895841, 0.0964735,
01769 0.0917, 0.0944815, 0.101367, 0.0972533, 0.152603,
01770 0.15325, 0.107244, 0.0944059, 0.095214, 0.0950256,
01771 0.0914785, 0.0831241, 0.0837255, 0.0873364, 0.0860993 },
01772 { 0.0779581, 0.0820391, 0.0916887, 0.0887483, 0.0827866,
01773 0.0827026, 0.0920711, 0.0967812, 0.103653, 0.143713,
01774 0.149082, 0.118462, 0.098462, 0.100916, 0.0892969,
01775 0.0852595, 0.0940911, 0.0742122, 0.0952267, 0.0936388 },
01776 { 0.0850528, 0.0791047, 0.0665815, 0.0837461, 0.0895629,
01777 0.0938629, 0.0976764, 0.102929, 0.102621, 0.137563,
01778 0.128113, 0.103627, 0.0927869, 0.102425, 0.0949018,
01779 0.0892192, 0.103132, 0.0932398, 0.0760799, 0.0901593 },
01780 { 0.0739873, 0.0657621, 0.0733051, 0.0810056, 0.0823362,
01781 0.080394, 0.082885, 0.0953372, 0.105809, 0.136333,
01782 0.142636, 0.101022, 0.0982706, 0.0851251, 0.0798855,
01783 0.0851645, 0.0799862, 0.087646, 0.0741014, 0.0812106 },
01784 { 0.0848307, 0.0738556, 0.0615958, 0.0647394, 0.0657774,
01785 0.0811213, 0.0891888, 0.0819757, 0.109486, 0.110607,
01786 0.117478, 0.104223, 0.0886376, 0.0773295, 0.0819245,
01787 0.0876723, 0.0725171, 0.0808546, 0.0813469, 0.095204 },
01788 { 0.0693136, 0.0975141, 0.0688994, 0.071388, 0.0883577,
01789 0.0654523, 0.0748456, 0.07202, 0.0828481, 0.115133,
01790 0.137037, 0.0930515, 0.0759633, 0.101175, 0.0873124,
01791 0.0629688, 0.0896685, 0.0692059, 0.0921279, 0.077337 },
01792 { 0.0686093, 0.0533264, 0.0434495, 0.0906161, 0.0764346,
01793 0.0809739, 0.0952387, 0.0845709, 0.0882732, 0.114535,
01794 0.0980442, 0.0891977, 0.0897545, 0.0827314, 0.0916396,
01795 0.127231, 0.0481388, 0.0897086, 0.0760457, 0.101561 },
01796 { 0.0740321, 0.0793183, 0.0488967, 0.0851262, 0.0764973,
01797 0.0579957, 0.121692, 0.0571608, 0.0703457, 0.113346,
01798 0.104348, 0.105721, 0.0784515, 0.0510833, 0.0689322,
01799 0.070756, 0.143634, 0.12002, 0.0406569, 0.0434415 },
01800 { 0.0325607, 0.170293, 0.0706702, 0.0635282, 0.0902266,
01801 0.118271, 0.0864529, 0.100342, 0.103515, 0.08514,
01802 0.113142, 0.0732865, 0.112173, 0.0644814, 0.0631694,
01803 0.0676042, 0.0594253, 0.102948, 0.0544197, 0.110599 },
01804 { 0.0396907, 0.0307508, 0.110389, 0.115979, 0.0810728,
01805 0.0655106, 0.042684, 0.0834926, 0.0970844, 0.0794593,
01806 0.114272, 0.0885057, 0.0870238, 0.121492, 0.0818473,
01807 0.108992, 0.0381046, 0.0274913, 0.0954279, 0.0635297 },
01808 { 0.015924, 0.0312376, 0.0517243, 0.0438412, 0.0503261,
01809 0.0501387, 0.105611, 0.0440046, 0.108186, 0.0893532,
01810 0.067842, 0.0789161, 0.0685462, 0.043805, 0.0707264,
01811 0.0298818, 0.10735, 0.0788188, 0.0314362, 0.0646983 }
01812 };
01813
01814 fBartolFlavourRatioNumuBar = new TH2D("hBartolFlavourRatioNumuBar","",EnergyBins,EnergyEdges,AngleBins,AngleEdges);
01815 for( Int_t nx=0; nx<EnergyBins; nx++ ){
01816 for( Int_t ny=0; ny<AngleBins; ny++ ){
01817 fBartolFlavourRatioNumuBar->SetBinContent(nx+1,ny+1,tempFlavourRatioNumuBar[nx][ny]);
01818 }
01819 }
01820
01821 return;
01822 }