loon: AtNuFluxRW Class Reference

AtNuFluxRW::AtNuFluxRW ( ) [private]

Definition at line 22 of file AtNuFluxRW.cxx.

References LoadBartolFlavourRatios(), LoadBartolFluxRatios(), LoadFlukaFluxRatios(), and LoadHondaFluxRatios().

00023 {
00024   this->LoadFlukaFluxRatios();
00025   this->LoadHondaFluxRatios();
00026   this->LoadBartolFluxRatios();
00027   this->LoadBartolFlavourRatios();
00028 }

AtNuFluxRW::~AtNuFluxRW ( ) [private]

Definition at line 30 of file AtNuFluxRW.cxx.

References fBartolFlavourRatioNumu, fBartolFlavourRatioNumuBar, fBartolFluxRatioNumu, fBartolFluxRatioNumuBar, fFlukaFluxRatioNumu, fFlukaFluxRatioNumuBar, fHondaFluxRatioNumu, and fHondaFluxRatioNumuBar.

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 }

TH2D * AtNuFluxRW::ApplyBartolFlavourRatios	(	Int_t	idnu,
		TH2D *	hFlux2D
	)

Definition at line 159 of file AtNuFluxRW.cxx.

References ApplyRatio(), fBartolFlavourRatioNumu, and fBartolFlavourRatioNumuBar.

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 }

TH2D * AtNuFluxRW::ApplyBartolFluxRatios	(	Int_t	idnu,
		TH2D *	hFlux2D
	)

Definition at line 146 of file AtNuFluxRW.cxx.

References ApplyRatio(), fBartolFluxRatioNumu, and fBartolFluxRatioNumuBar.

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 }

TH2D * AtNuFluxRW::ApplyRatio	(	TH2D *	hFlux2D,
		TH2D *	hRatio2D = `0`
	)			`[private]`

Definition at line 124 of file AtNuFluxRW.cxx.

References Anp::angle(), and LookupRatio().

Referenced by ApplyBartolFlavourRatios(), and ApplyBartolFluxRatios().

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 }

Double_t AtNuFluxRW::GetBartolFlavourRatio	(	Int_t	idnu,
		Double_t	energy,
		Double_t	coszen
	)

Definition at line 108 of file AtNuFluxRW.cxx.

References fBartolFlavourRatioNumu, fBartolFlavourRatioNumuBar, and LookupRatio().

Referenced by AtNuEventCounter::AddData(), and AtNuEventNtuple::WriteVariables().

00109 {
00110   // muon neutrinos
00111   if( idnu==+14 ){
00112     return LookupRatio(fBartolFlavourRatioNumu,energy,coszen);
00113   }
00114 
00115   // muon anti-neutrinos
00116   if( idnu==-14 ){
00117     return LookupRatio(fBartolFlavourRatioNumuBar,energy,coszen);
00118   }
00119 
00120   // anything else
00121   return 1.0;
00122 }

TH2D* AtNuFluxRW::GetBartolFlavourRatioNumu ( ) [inline]

Definition at line 27 of file AtNuFluxRW.h.

References fBartolFlavourRatioNumu.

00027 { return fBartolFlavourRatioNumu;    }

TH2D* AtNuFluxRW::GetBartolFlavourRatioNumuBar ( ) [inline]

Definition at line 28 of file AtNuFluxRW.h.

References fBartolFlavourRatioNumuBar.

00028 { return fBartolFlavourRatioNumuBar; }

Double_t AtNuFluxRW::GetBartolFluxRatio	(	Int_t	idnu,
		Double_t	energy,
		Double_t	coszen
	)

Definition at line 92 of file AtNuFluxRW.cxx.

References fBartolFluxRatioNumu, fBartolFluxRatioNumuBar, and LookupRatio().

Referenced by AtNuEventCounter::AddData(), AtNuOscFitTemplateMaker::Run(), and AtNuEventNtuple::WriteVariables().

00093 {
00094   // muon(tau) neutrinos
00095   if( idnu==+14 || idnu==+16 ){
00096     return LookupRatio(fBartolFluxRatioNumu,energy,coszen);
00097   }
00098 
00099   // muon(tau) anti-neutrinos
00100   if( idnu==-14 || idnu==-16 ){
00101     return LookupRatio(fBartolFluxRatioNumuBar,energy,coszen);
00102   }
00103 
00104   // anything else
00105   return 1.0;
00106 }

TH2D* AtNuFluxRW::GetBartolFluxRatioNumu ( ) [inline]

Definition at line 24 of file AtNuFluxRW.h.

References fBartolFluxRatioNumu.

00024 { return fBartolFluxRatioNumu;    }

TH2D* AtNuFluxRW::GetBartolFluxRatioNumuBar ( ) [inline]

Definition at line 25 of file AtNuFluxRW.h.

References fBartolFluxRatioNumuBar.

00025 { return fBartolFluxRatioNumuBar; }

Double_t AtNuFluxRW::GetFlukaFluxRatio	(	Int_t	idnu,
		Double_t	energy,
		Double_t	coszen
	)

Definition at line 76 of file AtNuFluxRW.cxx.

References fFlukaFluxRatioNumu, fFlukaFluxRatioNumuBar, and LookupRatio().

Referenced by AtNuEventCounter::AddData(), and AtNuOscFitTemplateMaker::Run().

00077 {
00078   // muon(tau) neutrinos
00079   if( idnu==+14 || idnu==+16 ){
00080     return LookupRatio(fFlukaFluxRatioNumu,energy,coszen);
00081   }  
00082   
00083   // muon(tau) neutrinos
00084   if( idnu==-14 || idnu==-16 ){
00085     return LookupRatio(fFlukaFluxRatioNumuBar,energy,coszen);
00086   }
00087 
00088   // anything else
00089   return 1.0;
00090 }

TH2D* AtNuFluxRW::GetFlukaFluxRatioNumu ( ) [inline]

Definition at line 21 of file AtNuFluxRW.h.

References fFlukaFluxRatioNumu.

00021 { return fFlukaFluxRatioNumu;     }

TH2D* AtNuFluxRW::GetFlukaFluxRatioNumuBar ( ) [inline]

Definition at line 22 of file AtNuFluxRW.h.

References fFlukaFluxRatioNumuBar.

00022 { return fFlukaFluxRatioNumuBar;  }

Double_t AtNuFluxRW::GetHondaFluxRatio	(	Int_t	idnu,
		Double_t	energy,
		Double_t	coszen
	)

Definition at line 60 of file AtNuFluxRW.cxx.

References fHondaFluxRatioNumu, fHondaFluxRatioNumuBar, and LookupRatio().

Referenced by AtNuEventCounter::AddData(), and AtNuOscFitTemplateMaker::Run().

00061 {
00062   // muon(tau) neutrinos
00063   if( idnu==+14 || idnu==+16 ){
00064     return LookupRatio(fHondaFluxRatioNumu,energy,coszen);
00065   }  
00066   
00067   // muon(tau) neutrinos
00068   if( idnu==-14 || idnu==-16 ){
00069     return LookupRatio(fHondaFluxRatioNumuBar,energy,coszen);
00070   }
00071 
00072   // anything else
00073   return 1.0;
00074 }

TH2D* AtNuFluxRW::GetHondaFluxRatioNumu ( ) [inline]

Definition at line 18 of file AtNuFluxRW.h.

References fHondaFluxRatioNumu.

00018 { return fHondaFluxRatioNumu;     }

TH2D* AtNuFluxRW::GetHondaFluxRatioNumuBar ( ) [inline]

Definition at line 19 of file AtNuFluxRW.h.

References fHondaFluxRatioNumuBar.

00019 { return fHondaFluxRatioNumuBar;  }

AtNuFluxRW * AtNuFluxRW::Instance ( ) [static]

Definition at line 11 of file AtNuFluxRW.cxx.

Referenced by AtNuEventCounter::AddData(), AtNuOscFitTemplateMaker::Run(), and AtNuEventNtuple::WriteVariables().

00012 {
00013   if( !fRWInstance ){
00014     fRWInstance = new AtNuFluxRW();
00015   }
00016 
00017   if( !fRWInstance ) assert(fRWInstance);
00018 
00019   return fRWInstance;
00020 }

void AtNuFluxRW::LoadBartolFlavourRatios ( ) [private]

Definition at line 1314 of file AtNuFluxRW.cxx.

References fBartolFlavourRatioNumu, and fBartolFlavourRatioNumuBar.

Referenced by AtNuFluxRW().

01315 {
01316   /*
01317   std::cout << "  AtNuFluxRW: Loading Bartol Flavour Ratios " << std::endl 
01318             << "    Reweight: Bartol3D SolarMax Numu(Bar) to Nue(Bar)" << std::endl;
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 }

void AtNuFluxRW::LoadBartolFluxRatios ( ) [private]

Definition at line 1012 of file AtNuFluxRW.cxx.

References fBartolFluxRatioNumu, and fBartolFluxRatioNumuBar.

Referenced by AtNuFluxRW().

01013 {
01014   /*
01015   std::cout << "  AtNuFluxRW: Loading Bartol Flux Ratios " << std::endl 
01016             << "    Reweight: Bartol96 SolarMin to Bartol3D SolarMax " << std::endl;
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 }

void AtNuFluxRW::LoadFlukaFluxRatios ( ) [private]

Definition at line 682 of file AtNuFluxRW.cxx.

References fFlukaFluxRatioNumu, and fFlukaFluxRatioNumuBar.

Referenced by AtNuFluxRW().

00683 {
00684   /*
00685   std::cout << "  AtNuFluxRW: Loading Fluka Flux Ratios " << std::endl 
00686             << "    Reweight: Bartol3D SolarMax to Fluka3D SolarMax " << std::endl;
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 }

void AtNuFluxRW::LoadHondaFluxRatios ( ) [private]

Definition at line 172 of file AtNuFluxRW.cxx.

References fHondaFluxRatioNumu, and fHondaFluxRatioNumuBar.

Referenced by AtNuFluxRW().

00173 {
00174   /*
00175   std::cout << "  AtNuFluxRW: Loading Honda Flux Ratios " << std::endl 
00176             << "    Reweight: Bartol3D SolarMax to Honda3D SolarMax " << std::endl;
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 }

Double_t AtNuFluxRW::LookupRatio	(	TH2D *	hFlux2D,
		Double_t	energy,
		Double_t	coszen
	)			`[private]`

Definition at line 45 of file AtNuFluxRW.cxx.

Referenced by ApplyRatio(), GetBartolFlavourRatio(), GetBartolFluxRatio(), GetFlukaFluxRatio(), and GetHondaFluxRatio().

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 }

TH2D* AtNuFluxRW::fBartolFlavourRatioNumu [private]

Definition at line 54 of file AtNuFluxRW.h.

Referenced by ApplyBartolFlavourRatios(), GetBartolFlavourRatio(), GetBartolFlavourRatioNumu(), LoadBartolFlavourRatios(), and ~AtNuFluxRW().

TH2D* AtNuFluxRW::fBartolFlavourRatioNumuBar [private]

Definition at line 55 of file AtNuFluxRW.h.

Referenced by ApplyBartolFlavourRatios(), GetBartolFlavourRatio(), GetBartolFlavourRatioNumuBar(), LoadBartolFlavourRatios(), and ~AtNuFluxRW().

TH2D* AtNuFluxRW::fBartolFluxRatioNumu [private]

Definition at line 51 of file AtNuFluxRW.h.

Referenced by ApplyBartolFluxRatios(), GetBartolFluxRatio(), GetBartolFluxRatioNumu(), LoadBartolFluxRatios(), and ~AtNuFluxRW().

TH2D* AtNuFluxRW::fBartolFluxRatioNumuBar [private]

Definition at line 52 of file AtNuFluxRW.h.

Referenced by ApplyBartolFluxRatios(), GetBartolFluxRatio(), GetBartolFluxRatioNumuBar(), LoadBartolFluxRatios(), and ~AtNuFluxRW().

TH2D* AtNuFluxRW::fFlukaFluxRatioNumu [private]

Definition at line 48 of file AtNuFluxRW.h.

Referenced by GetFlukaFluxRatio(), GetFlukaFluxRatioNumu(), LoadFlukaFluxRatios(), and ~AtNuFluxRW().

TH2D* AtNuFluxRW::fFlukaFluxRatioNumuBar [private]

Definition at line 49 of file AtNuFluxRW.h.

Referenced by GetFlukaFluxRatio(), GetFlukaFluxRatioNumuBar(), LoadFlukaFluxRatios(), and ~AtNuFluxRW().

TH2D* AtNuFluxRW::fHondaFluxRatioNumu [private]

Definition at line 45 of file AtNuFluxRW.h.

Referenced by GetHondaFluxRatio(), GetHondaFluxRatioNumu(), LoadHondaFluxRatios(), and ~AtNuFluxRW().

TH2D* AtNuFluxRW::fHondaFluxRatioNumuBar [private]

Definition at line 46 of file AtNuFluxRW.h.

Referenced by GetHondaFluxRatio(), GetHondaFluxRatioNumuBar(), LoadHondaFluxRatios(), and ~AtNuFluxRW().

AtNuFluxRW * AtNuFluxRW::fRWInstance = 0 [static, private]

Definition at line 57 of file AtNuFluxRW.h.

AtNuFluxRW Class Reference

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation


Public Member Functions
Double_t	GetHondaFluxRatio (Int_t idnu, Double_t energy, Double_t coszen)
Double_t	GetFlukaFluxRatio (Int_t idnu, Double_t energy, Double_t coszen)
Double_t	GetBartolFluxRatio (Int_t idnu, Double_t energy, Double_t coszen)
Double_t	GetBartolFlavourRatio (Int_t idnu, Double_t energy, Double_t coszen)
TH2D *	GetHondaFluxRatioNumu ()
TH2D *	GetHondaFluxRatioNumuBar ()
TH2D *	GetFlukaFluxRatioNumu ()
TH2D *	GetFlukaFluxRatioNumuBar ()
TH2D *	GetBartolFluxRatioNumu ()
TH2D *	GetBartolFluxRatioNumuBar ()
TH2D *	GetBartolFlavourRatioNumu ()
TH2D *	GetBartolFlavourRatioNumuBar ()
TH2D *	ApplyBartolFluxRatios (Int_t idnu, TH2D *hFlux2D)
TH2D *	ApplyBartolFlavourRatios (Int_t idnu, TH2D *hFlux2D)
Static Public Member Functions
static AtNuFluxRW *	Instance ()
Private Member Functions
	AtNuFluxRW ()
	~AtNuFluxRW ()
Double_t	LookupRatio (TH2D *hFlux2D, Double_t energy, Double_t coszen)
TH2D *	ApplyRatio (TH2D hFlux2D, TH2D hRatio2D=0)
void	LoadHondaFluxRatios ()
void	LoadFlukaFluxRatios ()
void	LoadBartolFluxRatios ()
void	LoadBartolFlavourRatios ()
Private Attributes
TH2D *	fHondaFluxRatioNumu
TH2D *	fHondaFluxRatioNumuBar
TH2D *	fFlukaFluxRatioNumu
TH2D *	fFlukaFluxRatioNumuBar
TH2D *	fBartolFluxRatioNumu
TH2D *	fBartolFluxRatioNumuBar
TH2D *	fBartolFlavourRatioNumu
TH2D *	fBartolFlavourRatioNumuBar
Static Private Attributes
static AtNuFluxRW *	fRWInstance = 0