AtNuFluxRW Class Reference

#include <AtNuFluxRW.h>

List of all members.

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 AtNuFluxRWInstance ()

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 AtNuFluxRWfRWInstance = 0

Detailed Description

Definition at line 8 of file AtNuFluxRW.h.


Constructor & Destructor Documentation

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]

Member Function Documentation

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 }


Member Data Documentation

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

Definition at line 57 of file AtNuFluxRW.h.


The documentation for this class was generated from the following files:

Generated on 25 Sep 2017 for loon by  doxygen 1.6.1