AtNuFluxRW.cxx

Go to the documentation of this file.

#include "AtNuFluxRW.h"

#include <iostream>
#include <cmath>
#include <cassert>

ClassImp(AtNuFluxRW)

AtNuFluxRW* AtNuFluxRW::fRWInstance = 0;

AtNuFluxRW* AtNuFluxRW::Instance()
{
  if( !fRWInstance ){
    fRWInstance = new AtNuFluxRW();
  }

  if( !fRWInstance ) assert(fRWInstance);

  return fRWInstance;
}

AtNuFluxRW::AtNuFluxRW()
{
  this->LoadFlukaFluxRatios();
  this->LoadHondaFluxRatios();
  this->LoadBartolFluxRatios();
  this->LoadBartolFlavourRatios();
}
  
AtNuFluxRW::~AtNuFluxRW()
{
  if( fHondaFluxRatioNumu )        delete fHondaFluxRatioNumu;
  if( fHondaFluxRatioNumuBar )     delete fHondaFluxRatioNumuBar;

  if( fFlukaFluxRatioNumu )        delete fFlukaFluxRatioNumu;
  if( fFlukaFluxRatioNumuBar )     delete fFlukaFluxRatioNumuBar;

  if( fBartolFluxRatioNumu )       delete fBartolFluxRatioNumu;
  if( fBartolFluxRatioNumuBar )    delete fBartolFluxRatioNumuBar;

  if( fBartolFlavourRatioNumu )    delete fBartolFlavourRatioNumu;
  if( fBartolFlavourRatioNumuBar ) delete fBartolFlavourRatioNumuBar;
}

Double_t AtNuFluxRW::LookupRatio(TH2D* hFlux2D, Double_t energy, Double_t coszen)
{
  if( hFlux2D==0 ) return 1.0;

  Int_t Ebin = hFlux2D->GetXaxis()->FindBin(energy);
  if( Ebin<1 ) Ebin = 1;
  if( Ebin>hFlux2D->GetXaxis()->GetNbins() ) Ebin = hFlux2D->GetXaxis()->GetNbins();

  Int_t Abin = hFlux2D->GetYaxis()->FindBin(coszen);
  if( Abin<1 ) Abin = 1;
  if( Abin>hFlux2D->GetYaxis()->GetNbins() ) Abin = hFlux2D->GetYaxis()->GetNbins();
    
  return hFlux2D->GetBinContent(Ebin,Abin);
}

Double_t AtNuFluxRW::GetHondaFluxRatio(Int_t idnu, Double_t energy, Double_t coszen)
{
  // muon(tau) neutrinos
  if( idnu==+14 || idnu==+16 ){
    return LookupRatio(fHondaFluxRatioNumu,energy,coszen);
  }  
  
  // muon(tau) neutrinos
  if( idnu==-14 || idnu==-16 ){
    return LookupRatio(fHondaFluxRatioNumuBar,energy,coszen);
  }

  // anything else
  return 1.0;
}
  
Double_t AtNuFluxRW::GetFlukaFluxRatio(Int_t idnu, Double_t energy, Double_t coszen)
{
  // muon(tau) neutrinos
  if( idnu==+14 || idnu==+16 ){
    return LookupRatio(fFlukaFluxRatioNumu,energy,coszen);
  }  
  
  // muon(tau) neutrinos
  if( idnu==-14 || idnu==-16 ){
    return LookupRatio(fFlukaFluxRatioNumuBar,energy,coszen);
  }

  // anything else
  return 1.0;
}

Double_t AtNuFluxRW::GetBartolFluxRatio(Int_t idnu, Double_t energy, Double_t coszen)
{
  // muon(tau) neutrinos
  if( idnu==+14 || idnu==+16 ){
    return LookupRatio(fBartolFluxRatioNumu,energy,coszen);
  }

  // muon(tau) anti-neutrinos
  if( idnu==-14 || idnu==-16 ){
    return LookupRatio(fBartolFluxRatioNumuBar,energy,coszen);
  }

  // anything else
  return 1.0;
}

Double_t AtNuFluxRW::GetBartolFlavourRatio(Int_t idnu, Double_t energy, Double_t coszen)
{
  // muon neutrinos
  if( idnu==+14 ){
    return LookupRatio(fBartolFlavourRatioNumu,energy,coszen);
  }

  // muon anti-neutrinos
  if( idnu==-14 ){
    return LookupRatio(fBartolFlavourRatioNumuBar,energy,coszen);
  }

  // anything else
  return 1.0;
}

TH2D* AtNuFluxRW::ApplyRatio(TH2D* hFlux2D, TH2D* hRatio2D)
{
  TH2D* hCorrectedFlux2D = (TH2D*)(hFlux2D->Clone());
  hCorrectedFlux2D->Reset();

  TString hName = hFlux2D->GetName();
  hName.Append("_Corrected");
  hCorrectedFlux2D->SetName(hName.Data());

  for( Int_t nEnergy=0; nEnergy<hFlux2D->GetXaxis()->GetNbins(); nEnergy++ ){
    for( Int_t nAngle=0; nAngle<hFlux2D->GetYaxis()->GetNbins(); nAngle++ ){
      Double_t energy = hFlux2D->GetXaxis()->GetBinCenter(nEnergy+1);
      Double_t angle = hFlux2D->GetYaxis()->GetBinCenter(nAngle+1);
      Double_t content = hFlux2D->GetBinContent(nEnergy+1,nAngle+1);
      Double_t ratio = LookupRatio(hRatio2D,energy,angle);
      hCorrectedFlux2D->SetBinContent(nEnergy+1,nAngle+1,content*ratio);
    }
  }

  return hCorrectedFlux2D;
}

TH2D* AtNuFluxRW::ApplyBartolFluxRatios(Int_t idnu, TH2D* hFlux2D)
{
  if( idnu==+14 || idnu==+16 ){
    return ApplyRatio(hFlux2D,fBartolFluxRatioNumu);
  }

  if( idnu==-14 || idnu==-16 ){
    return ApplyRatio(hFlux2D,fBartolFluxRatioNumuBar);
  }

  return ApplyRatio(hFlux2D);
}

TH2D* AtNuFluxRW::ApplyBartolFlavourRatios(Int_t idnu, TH2D* hFlux2D)
{
  if( idnu==+14 ){
    return ApplyRatio(hFlux2D,fBartolFlavourRatioNumu);
  }

  if( idnu==-14 ){
    return ApplyRatio(hFlux2D,fBartolFlavourRatioNumuBar);
  }

  return ApplyRatio(hFlux2D);
}

void AtNuFluxRW::LoadHondaFluxRatios()
{
  /*
  std::cout << "  AtNuFluxRW: Loading Honda Flux Ratios " << std::endl 
            << "    Reweight: Bartol3D SolarMax to Honda3D SolarMax " << std::endl;
  */

  Int_t EnergyBins = 50;
  Double_t EnergyEdges[51] = {
    0.1,
    0.125893,
    0.158489,
    0.199526,
    0.251189,
    0.316228,
    0.398107,
    0.501187,
    0.630957,
    0.794328,
    1,
    1.25893,
    1.58489,
    1.99526,
    2.51189,
    3.16228,
    3.98107,
    5.01187,
    6.30957,
    7.94328,
    10,
    12.5893,
    15.8489,
    19.9526,
    25.1189,
    31.6228,
    39.8107,
    50.1187,
    63.0957,
    79.4328,
    100,
    125.893,
    158.489,
    199.526,
    251.189,
    316.228,
    398.107,
    501.187,
    630.957,
    794.328,
    1000,
    1258.93,
    1584.89,
    1995.26,
    2511.89,
    3162.28,
    3981.07,
    5011.87,
    6309.57,
    7943.28,
    10000,
  };

  Int_t AngleBins = 20;
  Double_t AngleEdges[21] = {
    -1,
    -0.9,
    -0.8,
    -0.7,
    -0.6,
    -0.5,
    -0.4,
    -0.3,
    -0.2,
    -0.1,
    0,
    0.1,
    0.2,
    0.3,
    0.4,
    0.5,
    0.6,
    0.7,
    0.8,
    0.9,
    1,
  };

  Double_t tempFluxRatioNumu[50][20] = {
    { 0.720359, 0.793301, 0.899326, 0.923632, 0.88869, 0.890731,
      0.895898, 0.941403, 1.00064, 1.07675, 1.05964, 0.95933,
      0.889697, 0.824311, 0.772392, 0.747324, 0.708427, 0.700375,
      0.683478, 0.669825 },
    { 0.811819, 0.858748, 0.928462, 0.962145, 0.951695, 0.935997,
      0.972904, 1.0088, 1.0369, 1.11609, 1.0956, 1.04031,
      0.97184, 0.910671, 0.87014, 0.831496, 0.824536, 0.79331,
      0.784999, 0.772206 },
    { 0.901643, 0.938847, 0.986094, 1.00437, 0.980795, 0.992929,
      1.00864, 1.02331, 1.05646, 1.11394, 1.13605, 1.07346,
      1.00401, 0.957333, 0.934315, 0.909731, 0.905439, 0.893116,
      0.888796, 0.887069 },
    { 0.976218, 0.986363, 1.00512, 1.00628, 1.02265, 1.00834,
      1.01513, 1.03982, 1.04327, 1.06515, 1.06517, 1.04328,
      1.01365, 0.979517, 0.966544, 0.969086, 0.96095, 0.961108,
      0.969822, 0.949648 },
    { 1.01872, 1.01627, 1.02932, 1.03316, 1.03011, 1.02138,
      1.01611, 1.01047, 1.02486, 1.05918, 1.02993, 1.03229,
      1.01486, 0.998504, 1.01745, 1.00019, 1.0165, 1.00709,
      0.998735, 1.01248 },
    { 1.04757, 1.03383, 1.02004, 1.02277, 1.02824, 1.02772,
      1.01849, 1.0285, 1.00833, 1.03599, 1.05975, 1.01999,
      1.02301, 1.01764, 1.02007, 1.02611, 1.0395, 1.04548,
      1.05044, 1.04758 },
    { 1.04947, 1.05989, 1.0242, 1.03583, 1.03409, 1.03361,
      1.01578, 1.04418, 1.04013, 0.997514, 1.03242, 1.03338,
      1.04458, 1.03751, 1.05477, 1.04753, 1.05663, 1.05877,
      1.05929, 1.06716 },
    { 1.07568, 1.0654, 1.03961, 1.03068, 1.04431, 1.02462,
      1.02666, 1.05695, 1.03015, 1.09041, 1.04809, 1.03949,
      1.05707, 1.04793, 1.07593, 1.06837, 1.06618, 1.08101,
      1.08198, 1.09384 },
    { 1.06138, 1.06782, 1.04616, 1.04508, 1.04514, 1.03517,
      1.05694, 1.03784, 1.06191, 1.0501, 1.09622, 1.06553,
      1.07682, 1.06914, 1.07927, 1.082, 1.07267, 1.07801,
      1.08774, 1.11025 },
    { 1.09412, 1.06854, 1.04611, 1.06524, 1.05229, 1.07934,
      1.06993, 1.07127, 1.06371, 1.08644, 1.07538, 1.09415,
      1.08042, 1.093, 1.08035, 1.10018, 1.08149, 1.12184,
      1.12454, 1.08708 },
    { 1.11039, 1.10018, 1.09326, 1.08349, 1.06437, 1.09692,
      1.0888, 1.08211, 1.08283, 1.0978, 1.10365, 1.10362,
      1.10807, 1.09426, 1.11138, 1.09579, 1.11842, 1.10886,
      1.12435, 1.10071 },
    { 1.1361, 1.10531, 1.12103, 1.10344, 1.08005, 1.12129,
      1.1061, 1.10955, 1.10453, 1.09128, 1.06459, 1.08851,
      1.11315, 1.1012, 1.09238, 1.13819, 1.13852, 1.10766,
      1.12178, 1.12104 },
    { 1.13875, 1.14397, 1.13262, 1.11985, 1.13819, 1.1096,
      1.13455, 1.16484, 1.10533, 1.16723, 1.18577, 1.11948,
      1.13795, 1.14243, 1.12613, 1.15064, 1.14688, 1.12089,
      1.1336, 1.12485 },
    { 1.18039, 1.15135, 1.15731, 1.17255, 1.16296, 1.15735,
      1.16709, 1.14067, 1.14979, 1.12991, 1.19347, 1.10999,
      1.14749, 1.16314, 1.16297, 1.21201, 1.16578, 1.1863,
      1.18932, 1.20266 },
    { 1.19846, 1.18625, 1.17147, 1.13215, 1.1852, 1.17618,
      1.17061, 1.17287, 1.09525, 1.16026, 1.16056, 1.17877,
      1.12111, 1.18438, 1.19401, 1.18742, 1.16441, 1.18542,
      1.18331, 1.21773 },
    { 1.18228, 1.19926, 1.22963, 1.19432, 1.1596, 1.14937,
      1.15854, 1.18083, 1.17196, 1.19205, 1.19763, 1.19293,
      1.1861, 1.1911, 1.20507, 1.21879, 1.13275, 1.22823,
      1.16163, 1.23397 },
    { 1.23037, 1.20027, 1.18809, 1.17104, 1.17272, 1.19541,
      1.15537, 1.16861, 1.17049, 1.27104, 1.21755, 1.16005,
      1.20321, 1.19903, 1.16094, 1.13938, 1.21375, 1.20707,
      1.22883, 1.22645 },
    { 1.22187, 1.18382, 1.15814, 1.1994, 1.13506, 1.19049,
      1.13951, 1.17874, 1.16259, 1.21896, 1.13312, 1.13138,
      1.20175, 1.1914, 1.19718, 1.22886, 1.20109, 1.21387,
      1.17268, 1.18271 },
    { 1.2064, 1.21185, 1.17636, 1.13548, 1.18321, 1.18928,
      1.14983, 1.15034, 1.06072, 1.14785, 1.10597, 1.11695,
      1.14628, 1.16716, 1.14778, 1.18728, 1.17039, 1.1376,
      1.14722, 1.19284 },
    { 1.12981, 1.20517, 1.10963, 1.19575, 1.15659, 1.14774,
      1.07403, 1.14681, 1.15808, 1.03298, 1.12651, 1.167,
      1.15558, 1.08149, 1.16461, 1.19247, 1.18738, 1.24405,
      1.19572, 1.0725 },
    { 1.1511, 1.15231, 1.14933, 1.14589, 1.14508, 1.1447,
      1.14424, 1.13888, 1.12893, 1.11346, 1.11727, 1.13231,
      1.13415, 1.14208, 1.14571, 1.14495, 1.14798, 1.14751,
      1.15554, 1.14934 },
    { 1.14329, 1.1414, 1.14604, 1.14064, 1.13485, 1.13441,
      1.13029, 1.12785, 1.1186, 1.10503, 1.1105, 1.11887,
      1.12794, 1.13417, 1.13425, 1.13675, 1.13933, 1.14128,
      1.14472, 1.14024 },
    { 1.13441, 1.14003, 1.13291, 1.13382, 1.13016, 1.13221,
      1.12699, 1.12164, 1.11497, 1.10273, 1.10124, 1.11713,
      1.1273, 1.12773, 1.13036, 1.12587, 1.13491, 1.13656,
      1.13941, 1.13375 },
    { 1.12494, 1.12898, 1.13212, 1.12787, 1.12731, 1.12319,
      1.12358, 1.11564, 1.1124, 1.09107, 1.09365, 1.11407,
      1.12216, 1.1227, 1.1247, 1.12954, 1.12917, 1.13139,
      1.12865, 1.1263 },
    { 1.12473, 1.12036, 1.12515, 1.12121, 1.12552, 1.11977,
      1.12017, 1.11554, 1.11155, 1.10417, 1.10427, 1.11238,
      1.11665, 1.11846, 1.123, 1.12356, 1.12608, 1.12536,
      1.12016, 1.12335 },
    { 1.11929, 1.12043, 1.12485, 1.11824, 1.12171, 1.11996,
      1.12309, 1.12174, 1.11513, 1.11046, 1.1087, 1.11655,
      1.11417, 1.11852, 1.12062, 1.12059, 1.12305, 1.11882,
      1.1239, 1.11862 },
    { 1.11655, 1.11637, 1.11695, 1.11822, 1.12575, 1.12162,
      1.1204, 1.12111, 1.11985, 1.1116, 1.1141, 1.12162,
      1.11969, 1.1153, 1.11833, 1.12399, 1.1219, 1.11628,
      1.11423, 1.11661 },
    { 1.11135, 1.11244, 1.1154, 1.11615, 1.1088, 1.12134,
      1.1183, 1.12032, 1.12219, 1.11493, 1.11796, 1.12142,
      1.12113, 1.12485, 1.11878, 1.11905, 1.11432, 1.11331,
      1.11249, 1.11162 },
    { 1.09605, 1.09877, 1.10767, 1.10823, 1.11795, 1.11448,
      1.12246, 1.12601, 1.124, 1.11552, 1.11984, 1.12043,
      1.12625, 1.12131, 1.11771, 1.10814, 1.10896, 1.10441,
      1.10048, 1.09849 },
    { 1.08968, 1.08597, 1.1019, 1.10062, 1.11139, 1.11071,
      1.11827, 1.11996, 1.12631, 1.11418, 1.11608, 1.12855,
      1.11568, 1.11974, 1.10931, 1.11349, 1.10524, 1.09375,
      1.08888, 1.08871 },
    { 1.08272, 1.07979, 1.07296, 1.09641, 1.09896, 1.10429,
      1.11562, 1.12488, 1.12579, 1.12145, 1.11877, 1.12778,
      1.11316, 1.11345, 1.10274, 1.09863, 1.08844, 1.08766,
      1.0787, 1.07528 },
    { 1.05715, 1.07162, 1.08069, 1.07815, 1.09084, 1.1002,
      1.10575, 1.11338, 1.128, 1.11169, 1.11447, 1.12315,
      1.1134, 1.10299, 1.10264, 1.09237, 1.07953, 1.06554,
      1.07204, 1.05715 },
    { 1.05785, 1.06469, 1.06695, 1.06177, 1.07666, 1.08208,
      1.10408, 1.11009, 1.1213, 1.10846, 1.11855, 1.12361,
      1.11132, 1.09757, 1.07594, 1.08181, 1.07201, 1.07232,
      1.05339, 1.05775 },
    { 1.04628, 1.04894, 1.04511, 1.05601, 1.06488, 1.07203,
      1.09156, 1.10076, 1.12087, 1.11899, 1.10666, 1.11161,
      1.10041, 1.0856, 1.0616, 1.06136, 1.06424, 1.04722,
      1.03971, 1.04331 },
    { 1.02578, 1.03207, 1.03124, 1.03775, 1.04371, 1.06494,
      1.06497, 1.09948, 1.10943, 1.10405, 1.10746, 1.10762,
      1.09036, 1.07305, 1.06849, 1.04576, 1.04339, 1.02975,
      1.02125, 1.03481 },
    { 1.02402, 1.01825, 1.02115, 1.03183, 1.03533, 1.04203,
      1.06889, 1.07223, 1.09812, 1.10191, 1.10085, 1.10919,
      1.09066, 1.08796, 1.03291, 1.03575, 1.05052, 1.01271,
      1.00669, 1.01745 },
    { 0.97601, 1.01154, 1.02887, 1.00822, 1.01574, 1.04434,
      1.03381, 1.06895, 1.09126, 1.08, 1.08503, 1.08315,
      1.06469, 1.04636, 1.0384, 1.04315, 1.04777, 0.979799,
      1.01124, 0.987096 },
    { 0.980698, 1.01862, 1.0055, 1.02462, 0.985128, 0.995888,
      1.03434, 1.0326, 1.06082, 1.07983, 1.07919, 1.07433,
      1.04541, 1.04537, 1.02427, 1.02044, 0.979039, 0.998043,
      1.01255, 0.982553 },
    { 0.985548, 1.00081, 1.00725, 1.02392, 1.00049, 1.0117,
      1.00962, 1.02788, 1.0681, 1.06668, 1.07183, 1.05702,
      1.03429, 1.05098, 1.0235, 1.02786, 0.997753, 1.00043,
      1.00099, 0.980654 },
    { 0.95215, 0.989738, 0.963026, 0.947634, 0.987205, 1.00397,
      1.01354, 1.02783, 1.05208, 1.07734, 1.05458, 1.02989,
      1.03, 1.00584, 0.996533, 0.982034, 1.01694, 1.00216,
      0.99215, 0.962998 },
    { 1.01709, 0.972295, 1.00578, 1.00179, 0.99772, 1.01169,
      0.993987, 0.999361, 1.01396, 1.03763, 1.01768, 1.03304,
      1.03365, 1.01499, 0.984115, 0.997397, 1.00277, 1.02304,
      0.970194, 0.987445 },
    { 0.975049, 0.981989, 1.01819, 0.926574, 0.978372, 1.0284,
      0.984688, 1.00161, 0.99902, 1.0347, 1.02974, 0.985215,
      1.00539, 0.957156, 0.989713, 0.961819, 1.01739, 0.962689,
      0.992251, 0.926103 },
    { 1.00874, 1.0326, 0.964631, 1.02387, 0.953105, 0.980413,
      1.00585, 1.0298, 1.04538, 1.0401, 1.0487, 0.984558,
      0.995697, 0.991493, 0.959058, 1.03056, 0.983598, 0.984286,
      0.983021, 0.962649 },
    { 0.923048, 0.984893, 0.913356, 0.990271, 0.932158, 0.940045,
      0.939188, 0.994956, 0.988205, 1.01424, 1.01875, 0.985862,
      0.965603, 0.962185, 1.00052, 0.961401, 0.994707, 0.965186,
      0.998041, 1.06699 },
    { 0.952689, 1.00768, 0.916614, 0.941868, 0.950947, 0.973441,
      1.0287, 0.979704, 1.00533, 1.01974, 1.01249, 1.01507,
      0.997751, 0.955584, 0.964103, 1.00437, 0.956159, 0.988511,
      0.97202, 1.10934 },
    { 1.00621, 0.969905, 0.953438, 0.963676, 0.936194, 0.982752,
      0.980977, 0.960463, 0.962255, 0.987635, 0.91714, 1.00001,
      0.959234, 1.01537, 0.972583, 0.99775, 0.946224, 1.06324,
      0.984148, 0.970734 },
    { 0.951995, 1.0146, 0.966451, 1.08403, 0.959334, 0.986243,
      0.986639, 0.922884, 0.966221, 1.02831, 0.989717, 0.984058,
      0.986356, 0.956551, 0.949887, 0.953624, 0.976427, 1.01375,
      0.970966, 1.0026 },
    { 0.919405, 1.07264, 0.993215, 1.09949, 1.14462, 0.959278,
      1.05022, 0.917194, 1.01166, 1.01724, 1.01443, 0.978508,
      0.969163, 0.999111, 1.01185, 0.951327, 1.12472, 1.09681,
      1.06077, 1.08825 },
    { 0.972354, 0.948074, 1.10467, 0.821436, 0.963849, 0.906684,
      0.976617, 0.946898, 0.916566, 0.95836, 1.00307, 0.882556,
      1.24185, 0.92046, 1.05831, 0.938688, 0.981977, 1.05946,
      0.941017, 1.02242 },
    { 1.12965, 0.880083, 0.988796, 1.02966, 1.01482, 0.97416,
      1.0202, 1.0415, 1.05672, 0.934347, 1.02207, 1.03901,
      0.907022, 1.09341, 1.06995, 1.39554, 0.95927, 0.855533,
      0.817907, 0.869338 }
  };

  fHondaFluxRatioNumu = new TH2D("hHondaFluxRatioNumu","",EnergyBins,EnergyEdges,AngleBins,AngleEdges);
  for( Int_t nx=0; nx<EnergyBins; nx++ ){
    for( Int_t ny=0; ny<AngleBins; ny++ ){
      fHondaFluxRatioNumu->SetBinContent(nx+1,ny+1,tempFluxRatioNumu[nx][ny]);
    }
  }

  Double_t tempFluxRatioNumuBar[50][20] = {
    { 0.705133, 0.7943, 0.913706, 0.919652, 0.898913, 0.885134,
      0.894307, 0.932343, 0.985395, 1.03157, 1.02455, 0.958662,
      0.875577, 0.819564, 0.77819, 0.735101, 0.703953, 0.694607,
      0.676097, 0.666778 },
    { 0.794329, 0.866493, 0.949698, 0.970911, 0.948826, 0.949538,
      0.962661, 0.997567, 1.0531, 1.11652, 1.11494, 1.0234,
      0.954849, 0.899839, 0.848932, 0.818466, 0.798481, 0.785352,
      0.770943, 0.760943 },
    { 0.883217, 0.930264, 0.978082, 0.996879, 0.980149, 0.974922,
      0.999462, 1.02874, 1.08286, 1.12143, 1.10254, 1.047,
      1.00639, 0.94077, 0.930293, 0.899776, 0.881568, 0.884135,
      0.867732, 0.870368 },
    { 0.946934, 0.985803, 1.01274, 1.01208, 1.01284, 1.01805,
      1.00304, 1.03294, 1.05554, 1.0758, 1.12276, 1.04798,
      1.02107, 0.983328, 0.968824, 0.960838, 0.948438, 0.960309,
      0.954063, 0.953743 },
    { 1.00819, 1.01944, 1.01065, 1.01869, 1.01894, 1.01414,
      1.0123, 1.03468, 1.03833, 1.03174, 1.04626, 1.04002,
      1.01678, 1.014, 1.0033, 0.994082, 0.997439, 0.997324,
      1.01387, 1.02144 },
    { 1.03904, 1.03371, 1.01483, 1.03546, 1.02987, 1.02239,
      1.01477, 1.02417, 1.03699, 1.00359, 1.02797, 1.06239,
      1.02949, 1.03873, 1.04383, 1.0426, 1.03902, 1.04585,
      1.05546, 1.05332 },
    { 1.07246, 1.03317, 1.0204, 1.01604, 1.02214, 1.02477,
      1.02828, 1.03273, 1.02501, 1.02918, 1.05316, 1.05169,
      1.03752, 1.03116, 1.04735, 1.03585, 1.06831, 1.06725,
      1.07953, 1.07072 },
    { 1.06979, 1.04662, 1.02333, 1.04155, 1.02398, 1.03551,
      1.0422, 1.0317, 1.04971, 1.01329, 0.995565, 1.02665,
      1.06475, 1.04793, 1.05069, 1.05534, 1.07181, 1.08366,
      1.08138, 1.06107 },
    { 1.05827, 1.06336, 1.04431, 1.04371, 1.03866, 1.04363,
      1.03875, 1.04221, 1.04986, 1.05065, 1.02655, 1.01735,
      1.05924, 1.06284, 1.077, 1.08544, 1.08657, 1.11089,
      1.08088, 1.07522 },
    { 1.09838, 1.06833, 1.04944, 1.03933, 1.05021, 1.07651,
      1.06947, 1.05206, 1.07456, 1.04922, 1.03429, 1.09086,
      1.07231, 1.06563, 1.09761, 1.10441, 1.07589, 1.0984,
      1.10692, 1.10414 },
    { 1.09735, 1.05926, 1.08221, 1.06877, 1.06173, 1.04431,
      1.06163, 1.05545, 1.0499, 1.06126, 1.07627, 1.09785,
      1.08965, 1.0831, 1.08786, 1.09269, 1.07797, 1.09777,
      1.10847, 1.07936 },
    { 1.10494, 1.10999, 1.09035, 1.10165, 1.0757, 1.08716,
      1.11784, 1.10386, 1.07887, 1.05409, 1.06081, 1.07985,
      1.06992, 1.07563, 1.09316, 1.11355, 1.08677, 1.10125,
      1.10858, 1.10307 },
    { 1.06633, 1.10916, 1.08729, 1.0818, 1.10057, 1.11147,
      1.11002, 1.11307, 1.1051, 1.03308, 1.10983, 1.12574,
      1.11438, 1.12275, 1.12051, 1.11036, 1.13027, 1.09293,
      1.11491, 1.11677 },
    { 1.16918, 1.06589, 1.11561, 1.0767, 1.1046, 1.09654,
      1.09506, 1.07667, 1.11686, 1.06134, 1.07642, 1.0963,
      1.11021, 1.10293, 1.11079, 1.13937, 1.11544, 1.12877,
      1.10403, 1.15176 },
    { 1.13041, 1.12844, 1.11995, 1.11652, 1.103, 1.10452,
      1.08008, 1.11132, 1.11674, 1.04467, 1.0198, 1.06417,
      1.13236, 1.12685, 1.10334, 1.10725, 1.13853, 1.11495,
      1.12575, 1.11826 },
    { 1.13562, 1.13881, 1.14526, 1.12451, 1.09084, 1.13981,
      1.13376, 1.13784, 1.11226, 1.20303, 1.18051, 1.11957,
      1.14898, 1.12047, 1.10479, 1.16525, 1.14292, 1.11555,
      1.1329, 1.11313 },
    { 1.10223, 1.12731, 1.11671, 1.2029, 1.11916, 1.13736,
      1.12954, 1.18114, 1.14328, 1.18401, 1.14757, 1.13574,
      1.1303, 1.11546, 1.16661, 1.08987, 1.09791, 1.17782,
      1.11611, 1.1323 },
    { 1.12881, 1.17714, 1.11728, 1.1694, 1.17252, 1.14896,
      1.17745, 1.14006, 1.1291, 1.22367, 1.1102, 1.13685,
      1.09323, 1.13254, 1.12325, 1.12837, 1.14189, 1.14157,
      1.14519, 1.14733 },
    { 1.16973, 1.16469, 1.17574, 1.13778, 1.1252, 1.15677,
      1.189, 1.11608, 1.11831, 1.10882, 1.15818, 1.20677,
      1.13892, 1.17439, 1.16413, 1.13288, 1.10678, 1.11631,
      1.08629, 1.1376 },
    { 1.12826, 1.10288, 1.21341, 1.14021, 1.1797, 1.1559,
      1.16067, 1.13227, 1.14451, 1.15918, 1.21151, 1.11502,
      1.15246, 1.06628, 1.1391, 1.15752, 1.10684, 1.18234,
      1.15648, 1.13509 },
    { 1.14934, 1.14353, 1.14833, 1.13995, 1.13809, 1.13893,
      1.13442, 1.13488, 1.1259, 1.11419, 1.11284, 1.13253,
      1.13684, 1.1381, 1.14091, 1.14185, 1.14264, 1.14149,
      1.14814, 1.14853 },
    { 1.16095, 1.1593, 1.16022, 1.15405, 1.15332, 1.15007,
      1.14307, 1.14341, 1.13892, 1.1242, 1.12678, 1.13922,
      1.14437, 1.14569, 1.15234, 1.15359, 1.15419, 1.1595,
      1.16099, 1.1588 },
    { 1.17996, 1.17918, 1.17836, 1.16848, 1.17065, 1.16373,
      1.15307, 1.15678, 1.14122, 1.13179, 1.13469, 1.14693,
      1.15474, 1.15924, 1.1636, 1.16438, 1.17139, 1.17444,
      1.18003, 1.17848 },
    { 1.20275, 1.20148, 1.20205, 1.19573, 1.19472, 1.18742,
      1.17482, 1.16685, 1.1545, 1.12237, 1.13158, 1.16585,
      1.17009, 1.17332, 1.17742, 1.1848, 1.19301, 1.20006,
      1.20284, 1.20302 },
    { 1.23199, 1.23005, 1.22551, 1.216, 1.20596, 1.20472,
      1.19453, 1.18632, 1.17557, 1.1646, 1.16352, 1.17666,
      1.18679, 1.19638, 1.20046, 1.2107, 1.21719, 1.22444,
      1.22943, 1.23402 },
    { 1.26679, 1.25808, 1.25414, 1.24784, 1.23604, 1.22859,
      1.21677, 1.20914, 1.19875, 1.17636, 1.18044, 1.19569,
      1.20666, 1.2204, 1.22899, 1.23921, 1.24441, 1.25865,
      1.26412, 1.2668 },
    { 1.30729, 1.3022, 1.29234, 1.28425, 1.26724, 1.25899,
      1.25015, 1.22934, 1.22055, 1.19618, 1.19642, 1.21462,
      1.23154, 1.24583, 1.25785, 1.26872, 1.27616, 1.2898,
      1.3003, 1.31289 },
    { 1.3458, 1.33346, 1.32683, 1.31017, 1.30516, 1.29163,
      1.2781, 1.26183, 1.23963, 1.21289, 1.20949, 1.2383,
      1.25983, 1.27544, 1.28205, 1.30027, 1.32064, 1.32759,
      1.33107, 1.35272 },
    { 1.39339, 1.37818, 1.36722, 1.35793, 1.34069, 1.32341,
      1.30464, 1.28222, 1.26134, 1.23086, 1.23415, 1.26367,
      1.28575, 1.30516, 1.32728, 1.33689, 1.36279, 1.36977,
      1.38391, 1.4021 },
    { 1.43986, 1.42004, 1.41459, 1.39846, 1.38149, 1.37117,
      1.35151, 1.32025, 1.29475, 1.24535, 1.24385, 1.29182,
      1.31567, 1.34317, 1.3619, 1.37828, 1.38325, 1.40652,
      1.42174, 1.43721 },
    { 1.47846, 1.46473, 1.44732, 1.44225, 1.42023, 1.4048,
      1.38004, 1.35857, 1.32158, 1.265, 1.26025, 1.32548,
      1.35721, 1.37355, 1.39434, 1.42195, 1.44374, 1.45207,
      1.46786, 1.47445 },
    { 1.52119, 1.51495, 1.5031, 1.47081, 1.45149, 1.43951,
      1.40974, 1.37628, 1.34304, 1.29201, 1.28083, 1.33835,
      1.38144, 1.41119, 1.43139, 1.44877, 1.48598, 1.49132,
      1.51611, 1.53203 },
    { 1.56421, 1.54083, 1.54527, 1.51762, 1.50602, 1.47605,
      1.44575, 1.40985, 1.37555, 1.30113, 1.30346, 1.37167,
      1.41186, 1.44981, 1.4752, 1.50151, 1.53919, 1.53593,
      1.55018, 1.56263 },
    { 1.60936, 1.61102, 1.58039, 1.55951, 1.53931, 1.51283,
      1.47672, 1.43874, 1.4091, 1.31694, 1.31639, 1.41254,
      1.45811, 1.48385, 1.51517, 1.53574, 1.55859, 1.57769,
      1.58624, 1.62263 },
    { 1.6816, 1.60386, 1.63108, 1.63122, 1.56379, 1.55547,
      1.52758, 1.47433, 1.42696, 1.33888, 1.34102, 1.43307,
      1.49232, 1.52119, 1.53433, 1.5739, 1.59573, 1.60716,
      1.62192, 1.64429 },
    { 1.73812, 1.66965, 1.67677, 1.6745, 1.6444, 1.57884,
      1.53448, 1.52977, 1.47095, 1.36072, 1.35836, 1.48096,
      1.50554, 1.56879, 1.60711, 1.5996, 1.64558, 1.65923,
      1.68725, 1.69653 },
    { 1.73741, 1.74946, 1.71231, 1.68992, 1.63998, 1.647,
      1.58153, 1.58315, 1.51281, 1.38858, 1.38779, 1.50532,
      1.55179, 1.57645, 1.6258, 1.67598, 1.69954, 1.7439,
      1.72433, 1.721 },
    { 1.72621, 1.7482, 1.75391, 1.71555, 1.69618, 1.6894,
      1.62061, 1.59407, 1.55012, 1.43062, 1.42841, 1.55488,
      1.61292, 1.63951, 1.69458, 1.70176, 1.69348, 1.74217,
      1.72865, 1.76461 },
    { 1.81684, 1.82032, 1.80797, 1.74448, 1.7608, 1.68778,
      1.68566, 1.64367, 1.57199, 1.44767, 1.43786, 1.56677,
      1.61475, 1.66566, 1.68043, 1.76808, 1.74716, 1.80031,
      1.78198, 1.83992 },
    { 1.87268, 1.81514, 1.81341, 1.75813, 1.70021, 1.77617,
      1.71203, 1.68026, 1.60029, 1.47949, 1.47206, 1.58531,
      1.66951, 1.68274, 1.74688, 1.77939, 1.76911, 1.80966,
      1.7994, 1.87301 },
    { 1.82806, 1.82168, 1.87607, 1.82485, 1.76238, 1.72373,
      1.70933, 1.74603, 1.63884, 1.52924, 1.53744, 1.66926,
      1.69867, 1.77375, 1.80192, 1.79363, 1.86116, 1.80724,
      1.84731, 1.8975 },
    { 1.90363, 1.87933, 1.79027, 1.82264, 1.88811, 1.79203,
      1.78229, 1.70044, 1.69824, 1.52854, 1.58538, 1.67855,
      1.70352, 1.85734, 1.837, 1.91832, 1.84225, 1.93317,
      1.90808, 2.03093 },
    { 1.96605, 1.9209, 1.89505, 1.86872, 1.85034, 1.96523,
      1.84052, 1.76761, 1.72045, 1.58447, 1.61228, 1.7973,
      1.79289, 1.79712, 1.82767, 1.9419, 1.91372, 1.9548,
      1.99333, 1.92723 },
    { 2.00469, 2.07578, 2.02333, 1.92862, 1.93283, 1.9165,
      1.90044, 1.82469, 1.77742, 1.64828, 1.64222, 1.77863,
      1.80037, 1.83077, 1.88771, 2.03012, 1.87732, 1.74876,
      1.91666, 1.96412 },
    { 1.90766, 2.06549, 1.98681, 1.88428, 1.84486, 1.94434,
      1.81099, 1.85154, 1.69809, 1.69709, 1.66751, 1.73816,
      1.7965, 2.00486, 1.79924, 1.90534, 1.82175, 1.94305,
      2.1263, 1.90907 },
    { 1.87915, 1.99101, 1.88399, 2.09318, 1.87379, 2.06238,
      1.97225, 1.88009, 1.71774, 1.76548, 1.71688, 1.74072,
      1.9428, 1.97727, 2.04047, 1.94044, 1.91399, 2.05775,
      1.99772, 1.98832 },
    { 2.17599, 1.72374, 1.80553, 1.99612, 1.90734, 1.89137,
      1.87229, 1.98052, 1.84842, 1.667, 1.71548, 1.80627,
      1.96131, 1.79501, 2.00485, 1.87175, 2.08943, 2.15724,
      2.09258, 1.83488 },
    { 1.61225, 2.3388, 1.9408, 2.00414, 1.94543, 1.86327,
      1.85944, 2.12909, 1.93555, 1.59297, 1.84128, 1.85451,
      2.00904, 2.01206, 2.0765, 1.90297, 1.85861, 1.78392,
      1.9965, 1.95516 },
    { 1.69503, 2.17423, 2.08416, 2.03891, 1.70623, 2.03858,
      1.97139, 1.72152, 1.83812, 1.66838, 1.84196, 1.91136,
      1.93389, 2.00585, 2.03476, 1.96987, 1.80645, 1.93991,
      1.57272, 1.90624 },
    { 1.80686, 2.75457, 1.8771, 1.99271, 1.82314, 1.73241,
      1.79203, 1.7553, 1.96554, 1.69902, 1.7991, 1.80475,
      1.86286, 2.06388, 1.76461, 1.90267, 2.58264, 2.0764,
      1.767, 1.58713 }
  };

  fHondaFluxRatioNumuBar = new TH2D("hHondaFluxRatioNumuBar","",EnergyBins,EnergyEdges,AngleBins,AngleEdges);
  for( Int_t nx=0; nx<EnergyBins; nx++ ){
    for( Int_t ny=0; ny<AngleBins; ny++ ){
      fHondaFluxRatioNumuBar->SetBinContent(nx+1,ny+1,tempFluxRatioNumuBar[nx][ny]);
    }
  }

  return;
}

void AtNuFluxRW::LoadFlukaFluxRatios()
{
  /*
  std::cout << "  AtNuFluxRW: Loading Fluka Flux Ratios " << std::endl 
            << "    Reweight: Bartol3D SolarMax to Fluka3D SolarMax " << std::endl;
  */

  Int_t EnergyBins = 30;
  Double_t EnergyEdges[31] = {
    0.1,
    0.125893,
    0.158489,
    0.199526,
    0.251189,
    0.316228,
    0.398107,
    0.501187,
    0.630957,
    0.794328,
    1,
    1.25893,
    1.58489,
    1.99526,
    2.51189,
    3.16228,
    3.98107,
    5.01187,
    6.30957,
    7.94328,
    10,
    12.5893,
    15.8489,
    19.9526,
    25.1189,
    31.6228,
    39.8107,
    50.1187,
    63.0957,
    79.4328,
    100,
  };

  Int_t AngleBins = 20;
  Double_t AngleEdges[21] = {
    -1,
    -0.9,
    -0.8,
    -0.7,
    -0.6,
    -0.5,
    -0.4,
    -0.3,
    -0.2,
    -0.1,
    0,
    0.1,
    0.2,
    0.3,
    0.4,
    0.5,
    0.6,
    0.7,
    0.8,
    0.9,
    1,
  };

  Double_t tempFluxRatioNumu[30][20] = {
    { 0.840422, 0.87705, 0.934688, 0.976589, 0.974876, 0.997133,
      0.983382, 0.992384, 1.04525, 1.15244, 1.14142, 1.04213,
      0.984293, 0.927369, 0.881852, 0.867163, 0.830181, 0.827924,
      0.816066, 0.808012 },
    { 0.926969, 0.940956, 0.967104, 1.01872, 1.04036, 1.03584,
      1.05436, 1.04804, 1.06509, 1.17172, 1.15764, 1.10472,
      1.05144, 1.00346, 0.972511, 0.946645, 0.947392, 0.918824,
      0.916891, 0.908118 },
    { 0.97895, 0.992491, 1.00549, 1.04175, 1.04758, 1.0708,
      1.06923, 1.04677, 1.07574, 1.16438, 1.19801, 1.12455,
      1.06355, 1.02451, 1.00821, 0.993298, 0.995052, 0.984066,
      0.984916, 0.98428 },
    { 1.01067, 1.00498, 0.997878, 1.02038, 1.06626, 1.06147,
      1.05772, 1.05668, 1.06735, 1.13218, 1.14159, 1.09961,
      1.06503, 1.02501, 1.01173, 1.01959, 1.01161, 1.01055,
      1.02264, 0.999025 },
    { 1.02615, 1.00925, 1.00138, 1.02907, 1.05666, 1.05908,
      1.04621, 1.02356, 1.05518, 1.13481, 1.11541, 1.09418,
      1.0616, 1.03085, 1.04407, 1.02651, 1.04104, 1.02927,
      1.0223, 1.03323 },
    { 1.03955, 1.01353, 0.98098, 1.00772, 1.0415, 1.05275,
      1.03635, 1.03546, 1.03972, 1.10803, 1.14554, 1.08155,
      1.06196, 1.03651, 1.03095, 1.03401, 1.04684, 1.05032,
      1.05692, 1.05151 },
    { 1.02609, 1.02158, 0.969828, 1.00484, 1.03046, 1.03964,
      1.01265, 1.03366, 1.0545, 1.04331, 1.09412, 1.07823,
      1.06086, 1.03345, 1.04625, 1.03258, 1.04186, 1.04311,
      1.04528, 1.05227 },
    { 1.02814, 1.00388, 0.964053, 0.978136, 1.01453, 1.00621,
      0.993931, 1.01695, 1.01253, 1.09057, 1.06135, 1.04659,
      1.0381, 1.01049, 1.03628, 1.02236, 1.02211, 1.03681,
      1.03906, 1.05217 },
    { 0.985539, 0.978484, 0.948651, 0.966236, 0.984528, 0.983144,
      0.990259, 0.965815, 1.00339, 0.997712, 1.05529, 1.02844,
      1.01774, 0.994417, 1.00558, 1.00162, 0.993757, 1.00108,
      1.01096, 1.03329 },
    { 0.982514, 0.948947, 0.924249, 0.95685, 0.956261, 0.990715,
      0.963801, 0.959075, 0.96147, 0.984893, 0.984523, 1.00368,
      0.976896, 0.977057, 0.972528, 0.981443, 0.967363, 1.00533,
      1.0076, 0.974384 },
    { 0.95993, 0.944604, 0.938522, 0.942052, 0.93041, 0.968062,
      0.944421, 0.929909, 0.942288, 0.955893, 0.972053, 0.962979,
      0.958358, 0.938964, 0.960561, 0.940459, 0.960973, 0.95571,
      0.967179, 0.947774 },
    { 0.940157, 0.908715, 0.928857, 0.922184, 0.905714, 0.948219,
      0.919465, 0.916848, 0.916737, 0.91332, 0.898587, 0.905837,
      0.920511, 0.906117, 0.906796, 0.93594, 0.938358, 0.915237,
      0.923603, 0.924227 },
    { 0.906709, 0.906795, 0.907644, 0.902078, 0.916423, 0.899413,
      0.907459, 0.92709, 0.881327, 0.934491, 0.957184, 0.890575,
      0.902336, 0.904077, 0.898647, 0.909496, 0.90791, 0.891161,
      0.898456, 0.893693 },
    { 0.909814, 0.884431, 0.896884, 0.910102, 0.901812, 0.902066,
      0.902181, 0.881045, 0.886066, 0.877366, 0.929636, 0.854076,
      0.882906, 0.89292, 0.894411, 0.925695, 0.891105, 0.913367,
      0.912787, 0.923982 },
    { 0.900552, 0.890072, 0.878621, 0.84789, 0.89215, 0.886942,
      0.883239, 0.88863, 0.825906, 0.881006, 0.884023, 0.889832,
      0.847269, 0.889936, 0.893078, 0.884741, 0.867299, 0.890661,
      0.888197, 0.914311 },
    { 0.881396, 0.896532, 0.901616, 0.870362, 0.854736, 0.849712,
      0.865889, 0.889759, 0.864096, 0.8929, 0.902701, 0.884834,
      0.897773, 0.890881, 0.888527, 0.900204, 0.836041, 0.917734,
      0.869007, 0.922719 },
    { 0.90194, 0.881014, 0.917437, 0.889416, 0.880151, 0.894467,
      0.874528, 0.883313, 0.864998, 0.943504, 0.906853, 0.856194,
      0.910856, 0.898539, 0.859745, 0.843681, 0.90646, 0.917322,
      0.895364, 0.895972 },
    { 0.877034, 0.853173, 0.87813, 0.888662, 0.836146, 0.881792,
      0.838243, 0.880982, 0.84918, 0.889483, 0.830551, 0.829126,
      0.893938, 0.874343, 0.880488, 0.899203, 0.881769, 0.914314,
      0.840759, 0.842651 },
    { 0.861046, 0.868173, 0.882498, 0.828324, 0.867634, 0.875882,
      0.83913, 0.853538, 0.764237, 0.835009, 0.804898, 0.804856,
      0.847978, 0.849124, 0.839592, 0.865624, 0.85032, 0.85092,
      0.818272, 0.84908 },
    { 0.814619, 0.868039, 0.829126, 0.866669, 0.849354, 0.844058,
      0.782655, 0.857184, 0.8326, 0.746175, 0.810921, 0.838383,
      0.861179, 0.787891, 0.854677, 0.874375, 0.858519, 0.929754,
      0.859232, 0.7724 },
    { 0.839912, 0.838151, 0.854668, 0.826393, 0.84484, 0.841833,
      0.838396, 0.861452, 0.806416, 0.795237, 0.795214, 0.807704,
      0.856431, 0.83557, 0.841697, 0.844222, 0.828263, 0.853307,
      0.838909, 0.838539 },
    { 0.841239, 0.835541, 0.844869, 0.82072, 0.840709, 0.834053,
      0.83395, 0.860885, 0.796785, 0.779402, 0.781053, 0.795339,
      0.859116, 0.835604, 0.833107, 0.839841, 0.820427, 0.84344,
      0.836738, 0.841507 },
    { 0.838774, 0.835037, 0.829114, 0.814417, 0.836523, 0.826442,
      0.834585, 0.862857, 0.784558, 0.762781, 0.762801, 0.784694,
      0.862474, 0.834314, 0.823419, 0.832056, 0.814128, 0.831288,
      0.835001, 0.840306 },
    { 0.832603, 0.825802, 0.819527, 0.807237, 0.826646, 0.813892,
      0.835197, 0.859745, 0.777169, 0.743549, 0.745735, 0.776156,
      0.860269, 0.835033, 0.814564, 0.828506, 0.810048, 0.819239,
      0.824896, 0.833933 },
    { 0.823636, 0.813502, 0.806879, 0.800983, 0.818893, 0.803106,
      0.833114, 0.84876, 0.766483, 0.727132, 0.72711, 0.767074,
      0.850802, 0.831847, 0.805585, 0.817778, 0.804411, 0.807003,
      0.814046, 0.82229 },
    { 0.804692, 0.799141, 0.794776, 0.790953, 0.803637, 0.789965,
      0.826226, 0.834223, 0.755414, 0.709901, 0.70875, 0.756415,
      0.830037, 0.822851, 0.790829, 0.80359, 0.794334, 0.790506,
      0.801947, 0.803224 },
    { 0.782066, 0.776863, 0.777641, 0.779029, 0.789286, 0.775381,
      0.809164, 0.8014, 0.741867, 0.688316, 0.689865, 0.743061,
      0.801416, 0.805482, 0.773585, 0.789048, 0.781586, 0.777162,
      0.775198, 0.780561 },
    { 0.752632, 0.749598, 0.762307, 0.757796, 0.758955, 0.756511,
      0.782339, 0.75787, 0.721343, 0.666258, 0.668067, 0.720849,
      0.758386, 0.786926, 0.755133, 0.766845, 0.756555, 0.760879,
      0.748924, 0.750998 },
    { 0.711723, 0.711968, 0.738055, 0.724103, 0.741876, 0.72855,
      0.750054, 0.705974, 0.695288, 0.641971, 0.644456, 0.693081,
      0.704613, 0.749283, 0.730701, 0.735672, 0.724576, 0.735884,
      0.71195, 0.711646 },
    { 0.669569, 0.671043, 0.714502, 0.677259, 0.711067, 0.698983,
      0.701589, 0.636847, 0.662751, 0.6144, 0.615449, 0.664066,
      0.631455, 0.702509, 0.697745, 0.711669, 0.680096, 0.709219,
      0.671501, 0.667969 }
  };

  fFlukaFluxRatioNumu = new TH2D("hFlukaFluxRatioNumu","",EnergyBins,EnergyEdges,AngleBins,AngleEdges);
  for( Int_t nx=0; nx<EnergyBins; nx++ ){
    for( Int_t ny=0; ny<AngleBins; ny++ ){
      fFlukaFluxRatioNumu->SetBinContent(nx+1,ny+1,tempFluxRatioNumu[nx][ny]);
    }
  }

  Double_t tempFluxRatioNumuBar[30][20] = {
    { 0.813955, 0.880554, 0.953715, 0.963695, 0.979551, 0.977534,
      0.972933, 0.977053, 1.01003, 1.09478, 1.0922, 1.02583,
      0.958933, 0.91476, 0.877688, 0.842443, 0.81492, 0.812994,
      0.799772, 0.794956 },
    { 0.906648, 0.953322, 0.995764, 1.02309, 1.03268, 1.04295,
      1.03914, 1.03285, 1.06927, 1.16189, 1.16785, 1.0793,
      1.03057, 0.99315, 0.946973, 0.92694, 0.914156, 0.907307,
      0.898265, 0.892748 },
    { 0.967687, 0.989562, 1.00662, 1.03495, 1.04794, 1.05033,
      1.0605, 1.04846, 1.09523, 1.16135, 1.14993, 1.09518,
      1.06802, 1.01583, 1.00861, 0.985684, 0.973646, 0.981224,
      0.967247, 0.972892 },
    { 0.991858, 1.01151, 1.01583, 1.02905, 1.0576, 1.07458,
      1.04532, 1.04399, 1.07336, 1.12285, 1.18183, 1.09804,
      1.07076, 1.03873, 1.02046, 1.01448, 1.00555, 1.01874,
      1.01239, 1.01322 },
    { 1.02112, 1.01764, 0.992235, 1.01803, 1.04397, 1.04991,
      1.0389, 1.03678, 1.06082, 1.08782, 1.11203, 1.09425,
      1.05631, 1.04908, 1.02949, 1.01806, 1.0214, 1.0204,
      1.03596, 1.04356 },
    { 1.03206, 1.01305, 0.982932, 1.02285, 1.03961, 1.0437,
      1.02768, 1.01968, 1.06451, 1.0659, 1.10116, 1.12273,
      1.0602, 1.0577, 1.05059, 1.04353, 1.03921, 1.04453,
      1.05402, 1.0504 },
    { 1.04662, 0.997828, 0.974034, 0.990587, 1.01547, 1.03001,
      1.02316, 1.01201, 1.04113, 1.07662, 1.11572, 1.10016,
      1.05026, 1.02914, 1.03481, 1.01451, 1.04395, 1.04388,
      1.05508, 1.04726 },
    { 1.02127, 0.992092, 0.956757, 0.996096, 0.993806, 1.01754,
      1.00917, 0.987223, 1.03681, 1.03304, 1.0268, 1.0475,
      1.04834, 1.01599, 1.01068, 1.00583, 1.0196, 1.03302,
      1.03161, 1.01409 },
    { 0.984216, 0.982635, 0.954217, 0.975922, 0.982448, 0.993121,
      0.973113, 0.963795, 0.998541, 1.02415, 1.01221, 0.995614,
      1.00338, 0.993034, 1.00459, 1.00257, 1.00316, 1.02745,
      1.00104, 0.997168 },
    { 0.98761, 0.958462, 0.934194, 0.942014, 0.962714, 0.989743,
      0.96351, 0.937555, 0.975999, 0.969581, 0.9667, 1.01405,
      0.969301, 0.955468, 0.986127, 0.985209, 0.959008, 0.979765,
      0.988938, 0.98684 },
    { 0.950241, 0.917507, 0.934231, 0.937153, 0.936406, 0.922287,
      0.917466, 0.902756, 0.907597, 0.939338, 0.958714, 0.9671,
      0.9429, 0.929066, 0.93942, 0.937493, 0.922468, 0.941181,
      0.949811, 0.928142 },
    { 0.919697, 0.924763, 0.908541, 0.931705, 0.912169, 0.921021,
      0.931278, 0.909739, 0.894162, 0.891666, 0.903067, 0.906338,
      0.885935, 0.889556, 0.908473, 0.922189, 0.895008, 0.909152,
      0.914509, 0.913367 },
    { 0.861536, 0.896088, 0.880956, 0.886718, 0.904485, 0.906834,
      0.893743, 0.888632, 0.879259, 0.840086, 0.907926, 0.900853,
      0.887874, 0.895991, 0.897546, 0.891255, 0.900171, 0.873147,
      0.890441, 0.897142 },
    { 0.925017, 0.841357, 0.883133, 0.857694, 0.884814, 0.867021,
      0.86094, 0.840496, 0.862228, 0.840947, 0.855974, 0.847824,
      0.863373, 0.859504, 0.865201, 0.894354, 0.867702, 0.880532,
      0.862199, 0.906589 },
    { 0.883853, 0.876899, 0.869666, 0.867285, 0.868957, 0.851473,
      0.838532, 0.859709, 0.849183, 0.811147, 0.798071, 0.808903,
      0.8743, 0.867519, 0.84239, 0.857428, 0.874025, 0.857625,
      0.869227, 0.870517 },
    { 0.893944, 0.884452, 0.879849, 0.858588, 0.856081, 0.866517,
      0.881239, 0.887897, 0.846837, 0.923903, 0.911188, 0.856679,
      0.901609, 0.868033, 0.839657, 0.903914, 0.878702, 0.85994,
      0.878115, 0.874336 },
    { 0.876379, 0.903858, 0.875289, 0.936114, 0.892408, 0.885688,
      0.878715, 0.932544, 0.891969, 0.902956, 0.876699, 0.888734,
      0.891627, 0.863273, 0.897917, 0.854565, 0.842009, 0.916724,
      0.8826, 0.896507 },
    { 0.878963, 0.929851, 0.866767, 0.897388, 0.911348, 0.879581,
      0.897608, 0.882406, 0.866203, 0.938416, 0.85367, 0.875452,
      0.843557, 0.859121, 0.857549, 0.874313, 0.870554, 0.88027,
      0.896112, 0.890564 },
    { 0.897591, 0.908135, 0.904751, 0.869448, 0.865583, 0.880213,
      0.893247, 0.849962, 0.843268, 0.861255, 0.898265, 0.910786,
      0.866673, 0.879821, 0.88321, 0.868881, 0.844031, 0.85674,
      0.843284, 0.872949 },
    { 0.860163, 0.853109, 0.92984, 0.872262, 0.90232, 0.877299,
      0.866028, 0.855133, 0.853504, 0.901105, 0.94167, 0.829629,
      0.869636, 0.793252, 0.864147, 0.884708, 0.84564, 0.906451,
      0.892607, 0.86357 },
    { 0.867831, 0.877496, 0.875074, 0.870433, 0.866616, 0.861246,
      0.843616, 0.852774, 0.831459, 0.87798, 0.876609, 0.831759,
      0.853022, 0.842539, 0.864818, 0.868834, 0.871198, 0.871366,
      0.880759, 0.867405 },
    { 0.869339, 0.879671, 0.876674, 0.878389, 0.873719, 0.867182,
      0.847906, 0.852922, 0.829993, 0.881279, 0.882081, 0.827769,
      0.853044, 0.846821, 0.866685, 0.874014, 0.878079, 0.877752,
      0.880541, 0.868112 },
    { 0.877093, 0.8836, 0.880889, 0.883811, 0.881232, 0.869681,
      0.850304, 0.857053, 0.824268, 0.881685, 0.881795, 0.825359,
      0.856492, 0.852012, 0.86989, 0.877108, 0.886399, 0.880586,
      0.884923, 0.873408 },
    { 0.88204, 0.890896, 0.885366, 0.891302, 0.885061, 0.87298,
      0.855152, 0.862005, 0.821338, 0.873522, 0.877379, 0.822899,
      0.859078, 0.855205, 0.870471, 0.885673, 0.888663, 0.886682,
      0.889588, 0.88343 },
    { 0.887653, 0.897806, 0.890016, 0.887905, 0.888866, 0.871687,
      0.855953, 0.866537, 0.8189, 0.865546, 0.864905, 0.819642,
      0.866854, 0.858229, 0.868566, 0.892356, 0.888768, 0.890158,
      0.897351, 0.889105 },
    { 0.892832, 0.905664, 0.893565, 0.884684, 0.894954, 0.865424,
      0.854828, 0.872997, 0.815889, 0.843188, 0.846302, 0.81376,
      0.871165, 0.85856, 0.865709, 0.897249, 0.882239, 0.896822,
      0.910005, 0.892842 },
    { 0.894064, 0.923836, 0.900625, 0.874809, 0.895932, 0.858614,
      0.857579, 0.876727, 0.809252, 0.816875, 0.817175, 0.805324,
      0.878275, 0.855732, 0.857834, 0.896981, 0.869297, 0.89787,
      0.922487, 0.89789 },
    { 0.886917, 0.932699, 0.901197, 0.849552, 0.898207, 0.847619,
      0.853964, 0.886825, 0.796962, 0.779662, 0.77746, 0.796111,
      0.885423, 0.852868, 0.841331, 0.894837, 0.856341, 0.900216,
      0.931032, 0.89148 },
    { 0.870615, 0.948943, 0.900949, 0.828719, 0.893579, 0.827786,
      0.848829, 0.887204, 0.78355, 0.73544, 0.737145, 0.784997,
      0.889646, 0.848916, 0.830208, 0.891051, 0.831688, 0.901784,
      0.952889, 0.876058 },
    { 0.843181, 0.959424, 0.895298, 0.794099, 0.883033, 0.813929,
      0.851788, 0.896845, 0.77083, 0.683374, 0.681968, 0.769085,
      0.893734, 0.844724, 0.808426, 0.880979, 0.785459, 0.891535,
      0.960569, 0.841634 }
  };

  fFlukaFluxRatioNumuBar = new TH2D("hFlukaFluxRatioNumuBar","",EnergyBins,EnergyEdges,AngleBins,AngleEdges);
  for( Int_t nx=0; nx<EnergyBins; nx++ ){
    for( Int_t ny=0; ny<AngleBins; ny++ ){
      fFlukaFluxRatioNumuBar->SetBinContent(nx+1,ny+1,tempFluxRatioNumuBar[nx][ny]);
    }
  }

  return;
}

void AtNuFluxRW::LoadBartolFluxRatios()
{
  /*
  std::cout << "  AtNuFluxRW: Loading Bartol Flux Ratios " << std::endl 
            << "    Reweight: Bartol96 SolarMin to Bartol3D SolarMax " << std::endl;
  */
 
  Int_t EnergyBins = 50;
  Double_t EnergyEdges[51] = {
    0.1,
    0.125893,
    0.158489,
    0.199526,
    0.251189,
    0.316228,
    0.398107,
    0.501187,
    0.630957,
    0.794328,
    1,
    1.25893,
    1.58489,
    1.99526,
    2.51189,
    3.16228,
    3.98107,
    5.01187,
    6.30957,
    7.94328,
    10,
    12.5893,
    15.8489,
    19.9526,
    25.1189,
    31.6228,
    39.8107,
    50.1187,
    63.0957,
    79.4328,
    100,
    125.893,
    158.489,
    199.526,
    251.189,
    316.228,
    398.107,
    501.187,
    630.957,
    794.328,
    1000,
    1258.93,
    1584.89,
    1995.26,
    2511.89,
    3162.28,
    3981.07,
    5011.87,
    6309.57,
    7943.28,
    10000,
  };

  Int_t AngleBins = 12;
  Double_t AngleEdges[13] = {
    -1,
    -0.875,
    -0.625,
    -0.375,
    -0.2,
    -0.1,
    0,
    0.1,
    0.2,
    0.375,
    0.625,
    0.875,
    1,
  };

  Double_t tempFluxRatioNumu[50][12] = {
    { 0.407901, 0.87419, 0.616542, 0.688266, 0.915618, 1.58532,
      1.74734, 0.970917, 0.612829, 0.513663, 0.477047, 0.467546 },
    { 0.445901, 0.872082, 0.645515, 0.676216, 0.932318, 1.55963,
      1.73012, 0.955304, 0.616437, 0.532747, 0.501567, 0.504905 },
    { 0.492986, 0.848229, 0.681116, 0.694671, 0.92421, 1.52528,
      1.64609, 0.951425, 0.656118, 0.581295, 0.543958, 0.547079 },
    { 0.528652, 0.853431, 0.708676, 0.705465, 0.916687, 1.47216,
      1.63378, 0.972574, 0.68485, 0.621828, 0.593854, 0.601002 },
    { 0.561735, 0.831392, 0.722444, 0.722006, 0.893503, 1.34281,
      1.51554, 0.941985, 0.702086, 0.636388, 0.620098, 0.626743 },
    { 0.594873, 0.834399, 0.75474, 0.716433, 0.880281, 1.25631,
      1.33919, 0.919446, 0.714011, 0.669663, 0.648387, 0.657226 },
    { 0.620277, 0.829408, 0.780855, 0.719463, 0.832236, 1.21858,
      1.26482, 0.875715, 0.717299, 0.689813, 0.677408, 0.676059 },
    { 0.644834, 0.822422, 0.797919, 0.725557, 0.838948, 1.05533,
      1.16951, 0.861823, 0.727926, 0.700102, 0.684015, 0.693801 },
    { 0.691114, 0.840625, 0.816643, 0.737783, 0.799762, 1.03566,
      1.043, 0.816473, 0.729166, 0.717071, 0.730919, 0.717953 },
    { 0.719017, 0.850992, 0.836471, 0.752321, 0.816899, 0.973956,
      1.02467, 0.806696, 0.748375, 0.760005, 0.743555, 0.771294 },
    { 0.764811, 0.84433, 0.865975, 0.779598, 0.828076, 0.949407,
      0.974656, 0.829445, 0.778136, 0.799467, 0.78081, 0.817831 },
    { 0.801737, 0.871698, 0.890679, 0.817441, 0.874749, 0.969302,
      1.00796, 0.891016, 0.82688, 0.835656, 0.832187, 0.846983 },
    { 0.851308, 0.88483, 0.902341, 0.815974, 0.878693, 0.920148,
      0.918465, 0.87357, 0.834139, 0.857196, 0.86531, 0.895173 },
    { 0.904781, 0.933462, 0.895744, 0.864345, 0.902824, 0.960066,
      0.919705, 0.937197, 0.871655, 0.860518, 0.910428, 0.912769 },
    { 0.928794, 0.942555, 0.916024, 0.880167, 0.980136, 0.943384,
      0.950294, 0.912527, 0.908798, 0.893887, 0.927494, 0.937122 },
    { 0.975788, 0.906205, 0.954758, 0.929771, 0.944432, 0.923941,
      0.917028, 0.925558, 0.92427, 0.913486, 0.938226, 0.962792 },
    { 0.959312, 0.941718, 0.95912, 0.935943, 0.948323, 0.858838,
      0.893549, 0.958074, 0.909989, 0.976111, 0.928092, 0.965191 },
    { 0.950023, 0.982083, 0.954113, 0.936179, 0.947785, 0.899885,
      0.963736, 0.970708, 0.912453, 0.91452, 0.966331, 0.984249 },
    { 0.929443, 0.96308, 0.977871, 0.93561, 1.04203, 0.935262,
      0.971336, 0.989438, 0.934431, 0.991, 0.982813, 0.950919 },
    { 1.02293, 0.960648, 0.953228, 0.946973, 0.942216, 1.02765,
      0.945596, 0.935716, 0.941016, 0.928756, 0.922628, 1.06962 },
    { 0.975442, 0.94712, 0.963152, 0.909081, 0.961021, 0.944187,
      0.944215, 0.959489, 0.912573, 0.96302, 0.946721, 0.976709 },
    { 0.954521, 0.951504, 0.93789, 0.920042, 0.968064, 0.949156,
      0.94715, 0.969824, 0.919684, 0.937861, 0.951884, 0.954588 },
    { 0.980221, 0.980739, 0.933387, 0.895646, 0.947115, 0.944987,
      0.944963, 0.946951, 0.895734, 0.936308, 0.980007, 0.979505 },
    { 0.966654, 0.958888, 0.955399, 0.89338, 0.960503, 0.949842,
      0.947058, 0.961756, 0.893413, 0.953794, 0.958538, 0.966057 },
    { 0.967193, 0.935157, 0.943601, 0.925294, 0.945592, 0.954411,
      0.95444, 0.944864, 0.925321, 0.942915, 0.933895, 0.968163 },
    { 0.97509, 0.9379, 0.936289, 0.895395, 0.949033, 0.952325,
      0.953872, 0.947776, 0.900497, 0.936483, 0.937825, 0.974925 },
    { 0.974844, 0.949273, 0.951143, 0.888057, 0.950144, 0.930221,
      0.928132, 0.948617, 0.89038, 0.953059, 0.949053, 0.975188 },
    { 0.968934, 0.935001, 0.942496, 0.906561, 0.94847, 0.935474,
      0.932941, 0.949119, 0.904026, 0.939556, 0.936168, 0.968738 },
    { 0.973721, 0.9599, 0.945159, 0.890205, 0.966573, 0.940351,
      0.936725, 0.969651, 0.890462, 0.947307, 0.960404, 0.97169 },
    { 0.982156, 0.935489, 0.932213, 0.896294, 0.950343, 0.957662,
      0.95603, 0.948461, 0.897865, 0.931526, 0.936307, 0.982302 },
    { 0.977345, 0.970102, 0.963545, 0.896852, 0.964495, 0.942714,
      0.944971, 0.962793, 0.903185, 0.965066, 0.967328, 0.982911 },
    { 1.00439, 0.957662, 0.950131, 0.917232, 0.966896, 0.965724,
      0.96331, 0.971072, 0.918149, 0.948911, 0.962587, 1.00431 },
    { 0.967905, 1.00098, 0.967671, 0.89003, 0.971013, 0.961357,
      0.952684, 0.969011, 0.891553, 0.967983, 0.99899, 0.970126 },
    { 0.979338, 0.951973, 0.963948, 0.913014, 0.979198, 0.958417,
      0.969097, 0.987354, 0.915201, 0.969012, 0.951276, 0.983347 },
    { 1.04077, 0.942393, 0.958535, 0.900534, 0.985788, 0.968012,
      0.965025, 0.987405, 0.902251, 0.955214, 0.94414, 1.03586 },
    { 1.02074, 0.977487, 0.990401, 0.907411, 0.996908, 0.98743,
      0.988375, 0.986959, 0.89185, 0.99036, 0.978413, 1.02841 },
    { 1.02476, 1.01422, 0.95784, 0.898782, 0.994976, 1.01447,
      1.00977, 1.00243, 0.896561, 0.945832, 1.02195, 1.01565 },
    { 1.05813, 0.998507, 1.02375, 0.925399, 1.04316, 1.02827,
      1.02887, 1.03004, 0.914676, 1.00111, 1.01798, 1.05783 },
    { 1.05609, 1.00462, 0.992902, 0.916825, 1.01045, 1.01208,
      1.00721, 1.02104, 0.899057, 0.975705, 1.0157, 1.06022 },
    { 1.08674, 0.994028, 0.982287, 0.918874, 1.03912, 1.02558,
      1.04771, 1.06151, 0.920488, 0.982144, 0.955515, 1.07647 },
    { 1.04069, 1.0178, 1.02958, 0.981677, 1.06222, 1.05969,
      1.08046, 1.0426, 0.95391, 1.03907, 1.01133, 1.06558 },
    { 1.0953, 1.04574, 1.02943, 0.979644, 1.09715, 1.06773,
      1.07287, 1.11252, 0.988495, 1.04784, 1.03353, 1.13861 },
    { 1.05675, 1.02863, 1.03059, 0.95608, 1.04195, 1.04743,
      1.03884, 1.10631, 0.981417, 1.01581, 1.0475, 1.10792 },
    { 1.14721, 1.06448, 1.0757, 0.985877, 1.09084, 1.09804,
      1.09319, 1.09343, 0.994542, 1.03235, 1.03548, 1.02076 },
    { 1.12154, 1.08814, 1.05434, 0.958804, 1.06365, 1.05436,
      1.06191, 1.05345, 0.975105, 1.04579, 1.06025, 1.00413 },
    { 1.0512, 1.02273, 1.03443, 0.9821, 1.14604, 1.11062,
      1.19599, 1.10277, 0.970488, 1.01207, 0.982797, 1.07774 },
    { 1.09247, 0.983922, 0.995253, 0.997355, 1.11813, 1.06071,
      1.10208, 1.09786, 0.967704, 1.02564, 1.0107, 1.05864 },
    { 1.09168, 0.931565, 0.963358, 0.969301, 1.10096, 1.06495,
      1.06789, 1.13827, 0.952162, 1.01003, 0.890653, 0.957078 },
    { 1.09717, 1.03783, 1.11627, 0.987227, 1.26891, 1.15563,
      1.10412, 1.3178, 0.859991, 1.04903, 0.988425, 1.05713 },
    { 0.999478, 1.03798, 1.0156, 0.893401, 1.05337, 1.16683,
      1.06669, 1.07133, 0.954657, 0.874146, 1.14809, 1.23959 }
  };

  fBartolFluxRatioNumu = new TH2D("hBartolFluxRatioNumu","",EnergyBins,EnergyEdges,AngleBins,AngleEdges);
  for( Int_t nx=0; nx<EnergyBins; nx++ ){
    for( Int_t ny=0; ny<AngleBins; ny++ ){
      fBartolFluxRatioNumu->SetBinContent(nx+1,ny+1,tempFluxRatioNumu[nx][ny]);
    }
  }

  Double_t tempFluxRatioNumuBar[50][12] = {
    { 0.419909, 0.86696, 0.617805, 0.694347, 0.936305, 1.65452,
      1.80188, 0.977291, 0.621024, 0.521226, 0.486601, 0.478059 },
    { 0.457087, 0.860854, 0.640657, 0.685283, 0.904763, 1.5367,
      1.67528, 0.954441, 0.62423, 0.539869, 0.514669, 0.515865 },
    { 0.493547, 0.854435, 0.681477, 0.696709, 0.898733, 1.47652,
      1.6406, 0.964271, 0.647458, 0.574885, 0.550217, 0.549446 },
    { 0.540152, 0.841968, 0.706792, 0.702489, 0.906074, 1.45667,
      1.51981, 0.957375, 0.664831, 0.609647, 0.586539, 0.594554 },
    { 0.558002, 0.842949, 0.736069, 0.709655, 0.877461, 1.38931,
      1.50214, 0.925314, 0.686758, 0.639664, 0.618067, 0.612849 },
    { 0.599399, 0.834899, 0.746364, 0.71837, 0.844918, 1.30243,
      1.38521, 0.866158, 0.694154, 0.642933, 0.641561, 0.658374 },
    { 0.618166, 0.832418, 0.776161, 0.70901, 0.830014, 1.16974,
      1.23157, 0.845438, 0.702651, 0.686819, 0.670706, 0.682772 },
    { 0.655249, 0.822341, 0.786444, 0.717981, 0.805452, 1.1002,
      1.19497, 0.849365, 0.703929, 0.702722, 0.690835, 0.722096 },
    { 0.686463, 0.83479, 0.809993, 0.734854, 0.807546, 1.00753,
      1.09195, 0.852882, 0.727053, 0.718116, 0.719951, 0.738093 },
    { 0.698685, 0.841827, 0.832503, 0.755738, 0.793793, 1.00293,
      1.06289, 0.794962, 0.756311, 0.750372, 0.74752, 0.746921 },
    { 0.766557, 0.853241, 0.855347, 0.776714, 0.83853, 0.947186,
      0.968055, 0.809326, 0.762216, 0.77733, 0.786487, 0.819667 },
    { 0.79006, 0.842004, 0.857574, 0.776131, 0.836566, 0.946153,
      0.961994, 0.839372, 0.808635, 0.800766, 0.804363, 0.834603 },
    { 0.867501, 0.907733, 0.880754, 0.804795, 0.848948, 0.973121,
      0.921603, 0.834636, 0.808941, 0.842156, 0.872722, 0.86865 },
    { 0.835037, 0.907365, 0.937813, 0.859058, 0.87842, 0.966943,
      0.967174, 0.895318, 0.848295, 0.893034, 0.874765, 0.863815 },
    { 0.904532, 0.899746, 0.912705, 0.908259, 0.898149, 0.995856,
      1.02335, 0.94406, 0.886822, 0.888308, 0.887971, 0.928682 },
    { 0.902963, 0.890525, 0.917504, 0.88469, 0.924986, 0.876871,
      0.895394, 0.917661, 0.886617, 0.898717, 0.899884, 0.929002 },
    { 0.974293, 0.894307, 0.913058, 0.866757, 0.933625, 0.867731,
      0.893765, 0.937024, 0.897496, 0.9144, 0.907913, 0.961545 },
    { 0.933304, 0.913667, 0.875373, 0.875868, 0.943681, 0.849525,
      0.940958, 0.933712, 0.916019, 0.899732, 0.923929, 0.930549 },
    { 0.887209, 0.905562, 0.933587, 0.878906, 0.936487, 0.887606,
      0.851895, 0.867064, 0.87355, 0.931544, 0.949696, 0.921185 },
    { 0.958679, 0.879565, 0.930212, 0.904138, 0.923778, 0.903275,
      0.864285, 0.950362, 0.927507, 0.955342, 0.884841, 0.946768 },
    { 0.9197, 0.910948, 0.915674, 0.86876, 0.930058, 0.919897,
      0.92128, 0.929724, 0.868734, 0.91318, 0.911031, 0.919353 },
    { 0.90697, 0.918557, 0.911189, 0.871546, 0.919664, 0.927289,
      0.926447, 0.922136, 0.87174, 0.911396, 0.918127, 0.907802 },
    { 0.878294, 0.875183, 0.875293, 0.839296, 0.916957, 0.916563,
      0.916527, 0.915745, 0.838976, 0.87636, 0.874607, 0.880959 },
    { 0.895263, 0.853895, 0.858013, 0.835479, 0.927295, 0.910455,
      0.906606, 0.925536, 0.837389, 0.859093, 0.855112, 0.894417 },
    { 0.890479, 0.885531, 0.869181, 0.833928, 0.902262, 0.902183,
      0.903061, 0.901446, 0.83324, 0.86894, 0.885738, 0.88942 },
    { 0.86012, 0.846431, 0.860264, 0.820006, 0.89339, 0.895805,
      0.892723, 0.895727, 0.820076, 0.859227, 0.844828, 0.859259 },
    { 0.864054, 0.817477, 0.852692, 0.833401, 0.883778, 0.895342,
      0.895162, 0.888089, 0.833661, 0.85345, 0.820111, 0.8614 },
    { 0.830286, 0.822156, 0.819251, 0.799991, 0.881222, 0.901546,
      0.904078, 0.882164, 0.801415, 0.822587, 0.820294, 0.827235 },
    { 0.869974, 0.817073, 0.811808, 0.808597, 0.883987, 0.884445,
      0.882085, 0.882358, 0.807124, 0.811401, 0.81458, 0.864946 },
    { 0.807137, 0.80434, 0.831008, 0.797565, 0.878628, 0.868712,
      0.869756, 0.880621, 0.801195, 0.835496, 0.808747, 0.808106 },
    { 0.817639, 0.785971, 0.801099, 0.775249, 0.861652, 0.880509,
      0.883827, 0.859116, 0.777157, 0.803615, 0.784214, 0.819032 },
    { 0.793224, 0.791831, 0.778408, 0.770102, 0.864941, 0.860756,
      0.868271, 0.867974, 0.768037, 0.780019, 0.791509, 0.788658 },
    { 0.789749, 0.762709, 0.791717, 0.760806, 0.854404, 0.872839,
      0.871276, 0.856818, 0.759316, 0.791583, 0.759794, 0.789373 },
    { 0.777213, 0.755061, 0.768007, 0.776599, 0.87385, 0.887033,
      0.887403, 0.871721, 0.768832, 0.76777, 0.758978, 0.774709 },
    { 0.782047, 0.717187, 0.755737, 0.744113, 0.871312, 0.882519,
      0.881111, 0.867596, 0.739795, 0.760053, 0.72427, 0.794074 },
    { 0.750635, 0.731569, 0.758517, 0.693996, 0.856942, 0.898877,
      0.900438, 0.851152, 0.694862, 0.75975, 0.736655, 0.763372 },
    { 0.761118, 0.732715, 0.772226, 0.688345, 0.847226, 0.873483,
      0.873986, 0.851444, 0.697681, 0.770092, 0.728985, 0.769174 },
    { 0.774607, 0.73181, 0.736005, 0.680346, 0.838504, 0.877505,
      0.878862, 0.835941, 0.672439, 0.733923, 0.739178, 0.763102 },
    { 0.742569, 0.703755, 0.715088, 0.664591, 0.830326, 0.874947,
      0.880915, 0.833094, 0.674989, 0.716294, 0.708699, 0.738579 },
    { 0.737158, 0.689996, 0.730289, 0.701398, 0.836754, 0.881603,
      0.886054, 0.844662, 0.708866, 0.724037, 0.690787, 0.738451 },
    { 0.769987, 0.711012, 0.74531, 0.710789, 0.823083, 0.840787,
      0.836304, 0.808084, 0.712756, 0.722211, 0.714257, 0.745536 },
    { 0.751414, 0.71811, 0.718008, 0.707951, 0.811322, 0.874588,
      0.843233, 0.820841, 0.696452, 0.70208, 0.691169, 0.712002 },
    { 0.736834, 0.701795, 0.695325, 0.689812, 0.815717, 0.85955,
      0.844726, 0.780836, 0.689945, 0.705017, 0.680911, 0.74261 },
    { 0.701682, 0.675908, 0.682918, 0.657468, 0.773127, 0.830455,
      0.833518, 0.772601, 0.672326, 0.680118, 0.738788, 0.725297 },
    { 0.735046, 0.691412, 0.687795, 0.691496, 0.868924, 0.800936,
      0.815144, 0.848892, 0.678939, 0.69569, 0.700727, 0.730359 },
    { 0.749089, 0.68572, 0.673398, 0.656938, 0.867693, 0.745819,
      0.766935, 0.856235, 0.64287, 0.671985, 0.681261, 0.715962 },
    { 0.686553, 0.737141, 0.705175, 0.640149, 0.808821, 0.859627,
      0.835334, 0.827694, 0.654643, 0.694448, 0.641269, 0.75019 },
    { 0.833357, 0.668819, 0.704183, 0.617385, 0.78526, 0.911363,
      0.788458, 0.819574, 0.620585, 0.675242, 0.73846, 0.733908 },
    { 0.821739, 0.658415, 0.719287, 0.686069, 0.831316, 0.882383,
      0.799232, 0.799462, 0.632223, 0.686581, 0.771696, 0.80366 },
    { 0.786634, 0.66479, 0.750568, 0.718219, 0.813377, 0.879924,
      0.830972, 0.885847, 0.657697, 0.708878, 0.658219, 0.946438 }
  };

  fBartolFluxRatioNumuBar = new TH2D("hBartolFluxRatioNumuBar","",EnergyBins,EnergyEdges,AngleBins,AngleEdges);
  for( Int_t nx=0; nx<EnergyBins; nx++ ){
    for( Int_t ny=0; ny<AngleBins; ny++ ){
      fBartolFluxRatioNumuBar->SetBinContent(nx+1,ny+1,tempFluxRatioNumuBar[nx][ny]);
    }
  }

  return;
}

void AtNuFluxRW::LoadBartolFlavourRatios()
{
  /*
  std::cout << "  AtNuFluxRW: Loading Bartol Flavour Ratios " << std::endl 
            << "    Reweight: Bartol3D SolarMax Numu(Bar) to Nue(Bar)" << std::endl;
  */

  Int_t EnergyBins = 50;
  Double_t EnergyEdges[51] = {
    0.1,
    0.125893,
    0.158489,
    0.199526,
    0.251189,
    0.316228,
    0.398107,
    0.501187,
    0.630957,
    0.794328,
    1,
    1.25893,
    1.58489,
    1.99526,
    2.51189,
    3.16228,
    3.98107,
    5.01187,
    6.30957,
    7.94328,
    10,
    12.5893,
    15.8489,
    19.9526,
    25.1189,
    31.6228,
    39.8107,
    50.1187,
    63.0957,
    79.4328,
    100,
    125.893,
    158.489,
    199.526,
    251.189,
    316.228,
    398.107,
    501.187,
    630.957,
    794.328,
    1000,
    1258.93,
    1584.89,
    1995.26,
    2511.89,
    3162.28,
    3981.07,
    5011.87,
    6309.57,
    7943.28,
    10000,
  };

  Int_t AngleBins = 20;
  Double_t AngleEdges[21] = {
    -1,
    -0.9,
    -0.8,
    -0.7,
    -0.6,
    -0.5,
    -0.4,
    -0.3,
    -0.2,
    -0.1,
    0,
    0.1,
    0.2,
    0.3,
    0.4,
    0.5,
    0.6,
    0.7,
    0.8,
    0.9,
    1,
  };

  Double_t tempFlavourRatioNumu[50][20] = {
    { 0.483771, 0.482108, 0.468501, 0.467736, 0.487089,
      0.490678, 0.485871, 0.478179, 0.488448, 0.495532,
      0.487556, 0.493933, 0.497356, 0.490054, 0.49226,
      0.499126, 0.486597, 0.493977, 0.492738, 0.488128 },
    { 0.493145, 0.486648, 0.466921, 0.485041, 0.50345,
      0.49147, 0.517959, 0.515119, 0.507591, 0.536394,
      0.499488, 0.519409, 0.532653, 0.515699, 0.520559,
      0.51232, 0.528917, 0.511838, 0.505854, 0.503749 },
    { 0.50345, 0.488795, 0.47627, 0.489905, 0.49574,
      0.510479, 0.525588, 0.529513, 0.522196, 0.54615,
      0.558027, 0.533383, 0.530281, 0.526716, 0.531178,
      0.5302, 0.529989, 0.526617, 0.526447, 0.511257 },
    { 0.494412, 0.476955, 0.470938, 0.481523, 0.50026,
      0.500957, 0.50404, 0.522552, 0.51449, 0.514808,
      0.546352, 0.522192, 0.531277, 0.527599, 0.524156,
      0.533451, 0.522281, 0.520109, 0.510244, 0.497974 },
    { 0.485705, 0.480266, 0.477473, 0.484447, 0.50488,
      0.517759, 0.515405, 0.518466, 0.554715, 0.551937,
      0.543441, 0.534474, 0.548143, 0.534233, 0.538937,
      0.526618, 0.522721, 0.521193, 0.501988, 0.499907 },
    { 0.49551, 0.487358, 0.474522, 0.478369, 0.505141,
      0.524233, 0.540578, 0.547802, 0.549412, 0.572686,
      0.6097, 0.552332, 0.552638, 0.569476, 0.552663,
      0.533432, 0.539793, 0.518334, 0.524457, 0.502119 },
    { 0.485826, 0.481618, 0.467989, 0.48277, 0.496884,
      0.526026, 0.533962, 0.561812, 0.573883, 0.589635,
      0.599014, 0.593377, 0.574603, 0.562757, 0.565321,
      0.544403, 0.533308, 0.525284, 0.508585, 0.492155 },
    { 0.500213, 0.476787, 0.461245, 0.4819, 0.516985,
      0.520957, 0.537474, 0.576207, 0.581938, 0.591559,
      0.565733, 0.578725, 0.56797, 0.545511, 0.561957,
      0.551516, 0.523841, 0.51313, 0.502237, 0.495146 },
    { 0.45806, 0.464809, 0.452513, 0.479841, 0.491237,
      0.503463, 0.551025, 0.558151, 0.593469, 0.559609,
      0.590319, 0.588835, 0.582197, 0.565318, 0.54909,
      0.523059, 0.528094, 0.498051, 0.479855, 0.484975 },
    { 0.450572, 0.450117, 0.458697, 0.474582, 0.493752,
      0.534661, 0.560879, 0.567347, 0.593331, 0.644814,
      0.623036, 0.576097, 0.576903, 0.573431, 0.55282,
      0.52368, 0.51218, 0.493418, 0.49566, 0.454642 },
    { 0.428516, 0.433975, 0.453646, 0.457713, 0.490208,
      0.535391, 0.551684, 0.566821, 0.599828, 0.617741,
      0.62229, 0.5984, 0.592081, 0.552577, 0.546255,
      0.521942, 0.499915, 0.466678, 0.471527, 0.433905 },
    { 0.414213, 0.408831, 0.428097, 0.469648, 0.472421,
      0.512722, 0.53499, 0.594557, 0.591875, 0.605494,
      0.618075, 0.60744, 0.580938, 0.560856, 0.529507,
      0.511561, 0.488459, 0.447543, 0.431825, 0.417872 },
    { 0.387498, 0.395137, 0.422848, 0.439331, 0.469314,
      0.497287, 0.557325, 0.585158, 0.601038, 0.688688,
      0.638093, 0.587215, 0.568578, 0.55804, 0.525674,
      0.495038, 0.458604, 0.448664, 0.420186, 0.411906 },
    { 0.362212, 0.387219, 0.399575, 0.41824, 0.451681,
      0.468658, 0.538262, 0.550843, 0.602911, 0.606472,
      0.677418, 0.580113, 0.578858, 0.543293, 0.50707,
      0.465367, 0.435033, 0.426077, 0.403338, 0.381792 },
    { 0.335272, 0.335697, 0.380492, 0.366469, 0.433521,
      0.468824, 0.500071, 0.54058, 0.524332, 0.599026,
      0.598978, 0.573493, 0.553698, 0.49643, 0.470337,
      0.440934, 0.392765, 0.401862, 0.346422, 0.33904 },
    { 0.300426, 0.31122, 0.347986, 0.366456, 0.409579,
      0.475004, 0.473914, 0.522711, 0.605484, 0.595353,
      0.626054, 0.581413, 0.503464, 0.488057, 0.471954,
      0.397811, 0.363316, 0.371814, 0.315188, 0.303112 },
    { 0.281652, 0.307674, 0.293763, 0.320416, 0.359096,
      0.439461, 0.438454, 0.494095, 0.579221, 0.649487,
      0.632328, 0.559321, 0.521527, 0.456911, 0.375399,
      0.371075, 0.346508, 0.336367, 0.316047, 0.288439 },
    { 0.228568, 0.253805, 0.271087, 0.29186, 0.329319,
      0.407703, 0.411685, 0.467095, 0.552315, 0.692302,
      0.625927, 0.502497, 0.502373, 0.424052, 0.395593,
      0.364602, 0.322787, 0.297383, 0.276654, 0.240011 },
    { 0.246336, 0.249708, 0.235099, 0.295655, 0.281298,
      0.332034, 0.399846, 0.43911, 0.482904, 0.606731,
      0.538001, 0.472843, 0.442162, 0.390709, 0.340242,
      0.3173, 0.308261, 0.25273, 0.248806, 0.230625 },
    { 0.1631, 0.207982, 0.205912, 0.242623, 0.303877,
      0.288542, 0.338373, 0.437408, 0.521054, 0.597131,
      0.642662, 0.497307, 0.449432, 0.363256, 0.291221,
      0.321655, 0.234413, 0.240831, 0.229599, 0.207458 },
    { 0.169401, 0.184324, 0.200897, 0.222153, 0.249569,
      0.282551, 0.328894, 0.385506, 0.460316, 0.534085,
      0.534748, 0.462801, 0.385854, 0.327066, 0.282911,
      0.247291, 0.221632, 0.199857, 0.185385, 0.168774 },
    { 0.150167, 0.161729, 0.17881, 0.196742, 0.220957,
      0.252632, 0.293562, 0.352919, 0.431064, 0.512889,
      0.514071, 0.429936, 0.35205, 0.295427, 0.250656,
      0.219221, 0.195592, 0.177684, 0.161431, 0.150135 },
    { 0.132826, 0.144069, 0.155376, 0.174172, 0.194853,
      0.223914, 0.263341, 0.318823, 0.398018, 0.486413,
      0.490415, 0.396536, 0.318519, 0.261074, 0.223107,
      0.194106, 0.174984, 0.15696, 0.144289, 0.134147 },
    { 0.119721, 0.129023, 0.139478, 0.154179, 0.171923,
      0.19694, 0.234137, 0.285904, 0.363016, 0.458716,
      0.460095, 0.363648, 0.285564, 0.232981, 0.197246,
      0.173353, 0.154122, 0.14036, 0.127371, 0.118951 },
    { 0.108839, 0.113704, 0.123831, 0.137892, 0.152174,
      0.17414, 0.208053, 0.255379, 0.331325, 0.430478,
      0.431342, 0.329597, 0.253939, 0.206669, 0.17496,
      0.152544, 0.136938, 0.123935, 0.11439, 0.108862 },
    { 0.0973561, 0.103436, 0.111873, 0.121842, 0.134978,
      0.154286, 0.182153, 0.226978, 0.296657, 0.400109,
      0.400045, 0.297045, 0.224974, 0.180749, 0.153104,
      0.135703, 0.12208, 0.111167, 0.104326, 0.0995421 },
    { 0.0899228, 0.0953126, 0.101151, 0.110354, 0.121576,
      0.137762, 0.160898, 0.19966, 0.266258, 0.369011,
      0.369623, 0.266974, 0.199365, 0.161833, 0.136238,
      0.121946, 0.109746, 0.101529, 0.0953669, 0.091064 },
    { 0.0845741, 0.0881361, 0.0937866, 0.100352, 0.108946,
      0.123213, 0.143108, 0.176464, 0.237217, 0.337551,
      0.338259, 0.236706, 0.177374, 0.143615, 0.122609,
      0.109579, 0.0998899, 0.0943006, 0.0888237, 0.0846918 },
    { 0.0785875, 0.0810312, 0.0883395, 0.0917064, 0.100236,
      0.111409, 0.129469, 0.1586, 0.21045, 0.309442,
      0.310171, 0.208367, 0.155461, 0.129071, 0.109489,
      0.100451, 0.0940174, 0.0864788, 0.0830078, 0.0785639 },
    { 0.0744995, 0.0780492, 0.0813225, 0.0873548, 0.0934551,
      0.101645, 0.115324, 0.137047, 0.187218, 0.280671,
      0.279426, 0.186033, 0.139027, 0.115385, 0.0993897,
      0.0921002, 0.0859814, 0.0797118, 0.07803, 0.0766232 },
    { 0.0727529, 0.0733828, 0.0763106, 0.0804042, 0.0875126,
      0.0928214, 0.10437, 0.124952, 0.163905, 0.255326,
      0.254744, 0.165525, 0.12291, 0.104189, 0.0909762,
      0.085886, 0.0830063, 0.0789492, 0.0763853, 0.071196 },
    { 0.0669383, 0.0742784, 0.0745275, 0.076596, 0.0791858,
      0.0841418, 0.0954737, 0.111703, 0.148699, 0.229631,
      0.232904, 0.145686, 0.112031, 0.0978979, 0.0866845,
      0.0800304, 0.0759834, 0.0719356, 0.0731054, 0.0695679 },
    { 0.065199, 0.0675056, 0.0724757, 0.0699149, 0.0762228,
      0.079739, 0.0894119, 0.103028, 0.133076, 0.207396,
      0.212309, 0.129558, 0.10342, 0.0890477, 0.0818844,
      0.0775523, 0.0766142, 0.0696595, 0.0669001, 0.0666938 },
    { 0.064123, 0.0696213, 0.0655516, 0.0681806, 0.0738569,
      0.0781168, 0.0826077, 0.0945174, 0.114229, 0.194846,
      0.189752, 0.118109, 0.0922713, 0.083229, 0.0775152,
      0.0723676, 0.0686127, 0.0670341, 0.0640538, 0.0624709 },
    { 0.0591662, 0.0617666, 0.0615119, 0.0625595, 0.0693798,
      0.070784, 0.0763717, 0.0857307, 0.105883, 0.173673,
      0.174281, 0.10342, 0.0882207, 0.0774242, 0.0723012,
      0.0681079, 0.0638167, 0.0624689, 0.0604721, 0.0609709 },
    { 0.0558325, 0.057686, 0.0588582, 0.0665832, 0.0666299,
      0.0685581, 0.0737307, 0.0764614, 0.0964915, 0.16324,
      0.159981, 0.0995177, 0.0820128, 0.0743402, 0.0655594,
      0.0677487, 0.0670363, 0.0641988, 0.0596194, 0.0597144 },
    { 0.0587265, 0.0583574, 0.0607261, 0.0672015, 0.0623116,
      0.0654502, 0.0648915, 0.0760605, 0.0939325, 0.144998,
      0.143789, 0.088923, 0.0783105, 0.066246, 0.0710467,
      0.0619893, 0.0600036, 0.0534692, 0.0569893, 0.0561573 },
    { 0.0547225, 0.0596272, 0.0556315, 0.0568475, 0.0582189,
      0.0596421, 0.0703799, 0.0723212, 0.0803562, 0.129116,
      0.131531, 0.0858948, 0.0701371, 0.0620491, 0.0645727,
      0.0578932, 0.0585832, 0.0598886, 0.0559077, 0.0550951 },
    { 0.0552642, 0.0551161, 0.053484, 0.0591354, 0.0591139,
      0.059932, 0.0655763, 0.0678727, 0.0825035, 0.120444,
      0.118474, 0.0800536, 0.0659849, 0.0669378, 0.0606173,
      0.058444, 0.0576221, 0.0528608, 0.0556769, 0.0539962 },
    { 0.0489014, 0.0498292, 0.0475769, 0.0534821, 0.0549172,
      0.0559145, 0.0594824, 0.0625635, 0.0704311, 0.107317,
      0.109912, 0.0691772, 0.0630019, 0.0573904, 0.0587197,
      0.0556742, 0.0597354, 0.0520465, 0.0506518, 0.0496374 },
    { 0.0461725, 0.0528534, 0.0555651, 0.0557047, 0.0520061,
      0.060214, 0.0477867, 0.0589459, 0.0727363, 0.0973311,
      0.100299, 0.0693673, 0.0626349, 0.0622551, 0.0511697,
      0.0478783, 0.0493531, 0.0501689, 0.046332, 0.0582558 },
    { 0.0551047, 0.0498504, 0.0567048, 0.045668, 0.0483489,
      0.049971, 0.0519016, 0.0608866, 0.0646732, 0.0869788,
      0.0993125, 0.0616546, 0.0548603, 0.0621514, 0.0532224,
      0.0500748, 0.0494954, 0.0488039, 0.0493598, 0.0464522 },
    { 0.0578186, 0.0491193, 0.0521767, 0.0522645, 0.0500484,
      0.0458721, 0.062098, 0.0493422, 0.0583262, 0.0843645,
      0.0793672, 0.0646252, 0.0600312, 0.0530195, 0.051245,
      0.0598124, 0.0521268, 0.0508563, 0.0496217, 0.0590737 },
    { 0.0546895, 0.0470785, 0.0492463, 0.0401263, 0.0393656,
      0.0574015, 0.0425785, 0.0581796, 0.061947, 0.0794304,
      0.087147, 0.0586097, 0.0555622, 0.0532719, 0.050666,
      0.0613289, 0.0468641, 0.0475148, 0.0397562, 0.0490451 },
    { 0.0479479, 0.0525138, 0.0429685, 0.0499666, 0.0432402,
      0.0467873, 0.0499847, 0.050154, 0.0664756, 0.068267,
      0.0582761, 0.0602447, 0.051397, 0.0445861, 0.0556882,
      0.0453375, 0.0450909, 0.0404872, 0.057043, 0.0631275 },
    { 0.029464, 0.0388906, 0.0500278, 0.0442822, 0.0396983,
      0.0320882, 0.0536825, 0.0407348, 0.0453387, 0.0661033,
      0.0665525, 0.0531895, 0.0453006, 0.0705818, 0.0526684,
      0.0527037, 0.0451247, 0.0524748, 0.0545379, 0.070268 },
    { 0.0551849, 0.0458357, 0.0453332, 0.0490514, 0.0481249,
      0.052368, 0.0423143, 0.0575992, 0.0634309, 0.0660974,
      0.057854, 0.0495832, 0.0527541, 0.0430467, 0.0392717,
      0.0502788, 0.0372179, 0.0250528, 0.029136, 0.0571204 },
    { 0.0578826, 0.0250609, 0.0365661, 0.0539183, 0.0531626,
      0.085322, 0.057191, 0.0291048, 0.0593329, 0.0592411,
      0.0792843, 0.0364724, 0.0579378, 0.0600581, 0.0386804,
      0.0158845, 0.0504172, 0.0436397, 0.0410874, 0.0747608 },
    { 0.047278, 0.0108266, 0.0581433, 0.0435173, 0.0282802,
      0.0556926, 0.0608881, 0.0341592, 0.0468512, 0.0709938,
      0.0744972, 0.0344647, 0.0514944, 0.0616371, 0.0306792,
      0.0458638, 0.0402185, 0.0402181, 0.0544692, 0.0563191 },
    { 0.0488874, 0.016854, 0.0287852, 0.0486452, 0.0284928,
      0.0474454, 0.0162495, 0.0377561, 0.0634744, 0.0407885,
      0.060515, 0.0521133, 0.0392173, 0.0298258, 0.0295891,
      0.0556119, 0.0537151, 0.0231033, 0.00995481, 0.0147139 }
  };

  fBartolFlavourRatioNumu = new TH2D("hBartolFlavourRatioNumu","",EnergyBins,EnergyEdges,AngleBins,AngleEdges);
  for( Int_t nx=0; nx<EnergyBins; nx++ ){
    for( Int_t ny=0; ny<AngleBins; ny++ ){
      fBartolFlavourRatioNumu->SetBinContent(nx+1,ny+1,tempFlavourRatioNumu[nx][ny]);
    }
  }

  Double_t tempFlavourRatioNumuBar[50][20] = {
    { 0.412082, 0.446207, 0.465096, 0.460359, 0.453281,
      0.437887, 0.435633, 0.436573, 0.443728, 0.431594,
      0.42742, 0.424616, 0.426706, 0.422855, 0.426578,
      0.422978, 0.409966, 0.415512, 0.420574, 0.40902 },
    { 0.429984, 0.450796, 0.467453, 0.46128, 0.460112,
      0.471711, 0.456325, 0.461657, 0.443761, 0.46554,
      0.454107, 0.458411, 0.457236, 0.450232, 0.430269,
      0.443141, 0.428164, 0.435857, 0.428362, 0.424918 },
    { 0.433105, 0.446292, 0.466956, 0.472778, 0.459664,
      0.461006, 0.458654, 0.454493, 0.474067, 0.450731,
      0.432456, 0.443847, 0.464821, 0.443334, 0.451395,
      0.438165, 0.430345, 0.437035, 0.424072, 0.42411 },
    { 0.41547, 0.44927, 0.462523, 0.455846, 0.464306,
      0.463399, 0.437637, 0.452316, 0.449527, 0.443117,
      0.457429, 0.448325, 0.455673, 0.45393, 0.435496,
      0.441375, 0.428283, 0.427373, 0.419763, 0.422682 },
    { 0.423694, 0.436987, 0.448645, 0.458416, 0.45381,
      0.443066, 0.458565, 0.456408, 0.452544, 0.454817,
      0.453412, 0.452277, 0.451623, 0.449577, 0.444839,
      0.432767, 0.436776, 0.427372, 0.418028, 0.416475 },
    { 0.415695, 0.429557, 0.435888, 0.460423, 0.453169,
      0.453298, 0.452958, 0.461324, 0.474823, 0.481373,
      0.44374, 0.470241, 0.460486, 0.455262, 0.446434,
      0.432897, 0.431226, 0.425379, 0.423152, 0.414652 },
    { 0.414974, 0.41543, 0.433437, 0.451531, 0.442472,
      0.436401, 0.457848, 0.46187, 0.463416, 0.467512,
      0.468489, 0.472127, 0.468293, 0.448049, 0.440055,
      0.436197, 0.439518, 0.424458, 0.426239, 0.410102 },
    { 0.404801, 0.416695, 0.426806, 0.449237, 0.451902,
      0.449419, 0.458291, 0.453783, 0.475248, 0.461036,
      0.450015, 0.474451, 0.467286, 0.452803, 0.443891,
      0.439374, 0.429032, 0.416778, 0.415954, 0.3826 },
    { 0.386581, 0.417341, 0.419535, 0.447522, 0.451983,
      0.442891, 0.45073, 0.462137, 0.473213, 0.506575,
      0.472663, 0.470637, 0.456902, 0.452977, 0.448125,
      0.432245, 0.425759, 0.421255, 0.405408, 0.390853 },
    { 0.38914, 0.396454, 0.408412, 0.420317, 0.434366,
      0.448955, 0.469308, 0.463137, 0.506365, 0.473309,
      0.490495, 0.483602, 0.466677, 0.458701, 0.441111,
      0.441069, 0.416925, 0.401729, 0.392452, 0.386146 },
    { 0.364399, 0.376711, 0.395371, 0.415891, 0.414879,
      0.426647, 0.451585, 0.458355, 0.449111, 0.478984,
      0.480986, 0.492351, 0.469532, 0.44699, 0.445114,
      0.413849, 0.397488, 0.38728, 0.377112, 0.36034 },
    { 0.353007, 0.373926, 0.38936, 0.3944, 0.40981,
      0.43751, 0.456229, 0.461252, 0.491707, 0.451962,
      0.46186, 0.457081, 0.447181, 0.413908, 0.418922,
      0.420504, 0.390046, 0.380007, 0.365469, 0.348283 },
    { 0.303415, 0.363492, 0.358637, 0.376342, 0.4085,
      0.414435, 0.422241, 0.46363, 0.486646, 0.462166,
      0.448602, 0.494764, 0.487072, 0.447774, 0.415816,
      0.401511, 0.398458, 0.358477, 0.345809, 0.333924 },
    { 0.341528, 0.309434, 0.333661, 0.365176, 0.39094,
      0.407794, 0.420635, 0.419859, 0.470557, 0.470652,
      0.481248, 0.476548, 0.445981, 0.417633, 0.378496,
      0.388377, 0.362502, 0.355671, 0.310379, 0.305583 },
    { 0.29191, 0.314489, 0.318564, 0.350551, 0.336322,
      0.379333, 0.400711, 0.429679, 0.437953, 0.455759,
      0.401319, 0.454079, 0.448966, 0.405334, 0.374829,
      0.3625, 0.35585, 0.30944, 0.299335, 0.274751 },
    { 0.293613, 0.294652, 0.311178, 0.322521, 0.340589,
      0.367259, 0.379857, 0.423452, 0.466335, 0.470098,
      0.43711, 0.46497, 0.403181, 0.396796, 0.365467,
      0.355684, 0.321824, 0.288057, 0.275021, 0.281563 },
    { 0.239426, 0.273545, 0.313211, 0.308953, 0.315088,
      0.373761, 0.414211, 0.420854, 0.464901, 0.481912,
      0.487575, 0.458679, 0.418835, 0.389034, 0.357844,
      0.334657, 0.298711, 0.299819, 0.265162, 0.243532 },
    { 0.214475, 0.257205, 0.26995, 0.291824, 0.334321,
      0.328389, 0.375506, 0.402932, 0.445746, 0.500525,
      0.425476, 0.460156, 0.380677, 0.377773, 0.317102,
      0.312775, 0.279455, 0.280969, 0.264316, 0.221291 },
    { 0.202441, 0.225885, 0.255472, 0.267166, 0.287013,
      0.327647, 0.366529, 0.410451, 0.403483, 0.479243,
      0.419496, 0.482984, 0.402491, 0.382375, 0.348764,
      0.257041, 0.254401, 0.236084, 0.216587, 0.213084 },
    { 0.179954, 0.193844, 0.222862, 0.230639, 0.288656,
      0.286669, 0.336902, 0.362264, 0.45466, 0.49592,
      0.477066, 0.426637, 0.374788, 0.308645, 0.311289,
      0.276291, 0.228307, 0.20893, 0.238115, 0.23494 },
    { 0.174504, 0.185824, 0.202585, 0.218628, 0.239263,
      0.268425, 0.300748, 0.34331, 0.396976, 0.447518,
      0.446782, 0.39663, 0.345322, 0.302963, 0.267699,
      0.241377, 0.218861, 0.200956, 0.186883, 0.173887 },
    { 0.159006, 0.17012, 0.182721, 0.199817, 0.222308,
      0.24704, 0.279989, 0.322346, 0.379958, 0.434932,
      0.436291, 0.379098, 0.324053, 0.279305, 0.246665,
      0.220629, 0.200324, 0.182691, 0.169508, 0.15846 },
    { 0.145755, 0.154856, 0.166888, 0.181873, 0.200191,
      0.224519, 0.25634, 0.301256, 0.358876, 0.421369,
      0.421391, 0.361286, 0.299859, 0.257893, 0.225611,
      0.200406, 0.180594, 0.16616, 0.155017, 0.14533 },
    { 0.133866, 0.143334, 0.152407, 0.165488, 0.182148,
      0.205248, 0.23422, 0.279239, 0.337267, 0.404417,
      0.406198, 0.337136, 0.276908, 0.234852, 0.204838,
      0.182104, 0.165932, 0.152539, 0.143214, 0.133841 },
    { 0.123566, 0.131634, 0.139287, 0.151834, 0.16804,
      0.18802, 0.214236, 0.256979, 0.316113, 0.390613,
      0.386573, 0.317894, 0.25559, 0.215569, 0.187054,
      0.168076, 0.153555, 0.139994, 0.131555, 0.125902 },
    { 0.115897, 0.121767, 0.131581, 0.139697, 0.151212,
      0.171003, 0.19538, 0.233541, 0.293581, 0.370069,
      0.369491, 0.293331, 0.233459, 0.194641, 0.172387,
      0.153263, 0.139037, 0.129894, 0.123654, 0.116666 },
    { 0.111004, 0.116278, 0.121219, 0.129864, 0.140257,
      0.155804, 0.178266, 0.213658, 0.271505, 0.349669,
      0.351192, 0.269197, 0.213176, 0.179079, 0.157679,
      0.141758, 0.12908, 0.121258, 0.116909, 0.112008 },
    { 0.106056, 0.113147, 0.114718, 0.12155, 0.131086,
      0.144228, 0.162865, 0.19482, 0.246457, 0.330992,
      0.327864, 0.247094, 0.193329, 0.160031, 0.142282,
      0.130932, 0.122047, 0.115993, 0.110271, 0.107866 },
    { 0.105538, 0.106666, 0.110969, 0.116491, 0.121385,
      0.134591, 0.149833, 0.175646, 0.22705, 0.309264,
      0.310382, 0.226335, 0.177202, 0.14874, 0.134385,
      0.122798, 0.114333, 0.113435, 0.106166, 0.104194 },
    { 0.104183, 0.10512, 0.107658, 0.112914, 0.115858,
      0.123622, 0.141049, 0.165866, 0.20915, 0.290458,
      0.290746, 0.208201, 0.162713, 0.136834, 0.127399,
      0.118285, 0.109107, 0.104548, 0.102973, 0.104227 },
    { 0.0985275, 0.103386, 0.105375, 0.107694, 0.112493,
      0.119188, 0.132633, 0.150316, 0.189637, 0.271861,
      0.271673, 0.188251, 0.154303, 0.131769, 0.118001,
      0.114698, 0.106518, 0.106363, 0.102147, 0.0984905 },
    { 0.0986039, 0.0983164, 0.104228, 0.100105, 0.106353,
      0.113724, 0.122595, 0.138818, 0.177965, 0.252492,
      0.252904, 0.172283, 0.139934, 0.122001, 0.113189,
      0.10818, 0.105292, 0.0996621, 0.0995675, 0.0976611 },
    { 0.0952388, 0.0959157, 0.102642, 0.104566, 0.107753,
      0.111432, 0.117709, 0.131835, 0.15903, 0.239543,
      0.240279, 0.160642, 0.130822, 0.11649, 0.107555,
      0.105715, 0.103211, 0.0982097, 0.103424, 0.0970446 },
    { 0.0966725, 0.0984501, 0.101105, 0.100576, 0.0991159,
      0.109583, 0.11332, 0.12501, 0.147122, 0.224034,
      0.220114, 0.150081, 0.124032, 0.115834, 0.105434,
      0.103624, 0.105005, 0.101644, 0.0939662, 0.0933233 },
    { 0.099124, 0.0934924, 0.0917282, 0.101118, 0.104825,
      0.098954, 0.115263, 0.116816, 0.137998, 0.213482,
      0.207205, 0.138229, 0.119599, 0.110535, 0.104773,
      0.101785, 0.0996848, 0.0961159, 0.0906688, 0.0951277 },
    { 0.0922214, 0.0924557, 0.09404, 0.101346, 0.104847,
      0.0987602, 0.105733, 0.110381, 0.124784, 0.192419,
      0.195219, 0.127529, 0.1085, 0.107653, 0.108114,
      0.0984656, 0.101836, 0.0979752, 0.0880972, 0.0926429 },
    { 0.0869145, 0.0962932, 0.0990207, 0.0951642, 0.0909166,
      0.101027, 0.100303, 0.116116, 0.121938, 0.181358,
      0.187503, 0.12106, 0.112041, 0.101635, 0.106341,
      0.100281, 0.0836268, 0.0945236, 0.0882399, 0.084966 },
    { 0.0833919, 0.0833637, 0.0852624, 0.0929501, 0.0855587,
      0.0982595, 0.103193, 0.0989452, 0.120303, 0.171511,
      0.169597, 0.116081, 0.100552, 0.0975868, 0.102001,
      0.0992304, 0.0995348, 0.0894016, 0.0808943, 0.0871339 },
    { 0.0849254, 0.100271, 0.0852303, 0.0851574, 0.0922799,
      0.0914357, 0.103621, 0.101908, 0.110508, 0.161002,
      0.163773, 0.110799, 0.0911506, 0.101228, 0.0935324,
      0.0913888, 0.0829404, 0.0949507, 0.0862389, 0.0890564 },
    { 0.0866945, 0.0736328, 0.080485, 0.0895841, 0.0964735,
      0.0917, 0.0944815, 0.101367, 0.0972533, 0.152603,
      0.15325, 0.107244, 0.0944059, 0.095214, 0.0950256,
      0.0914785, 0.0831241, 0.0837255, 0.0873364, 0.0860993 },
    { 0.0779581, 0.0820391, 0.0916887, 0.0887483, 0.0827866,
      0.0827026, 0.0920711, 0.0967812, 0.103653, 0.143713,
      0.149082, 0.118462, 0.098462, 0.100916, 0.0892969,
      0.0852595, 0.0940911, 0.0742122, 0.0952267, 0.0936388 },
    { 0.0850528, 0.0791047, 0.0665815, 0.0837461, 0.0895629,
      0.0938629, 0.0976764, 0.102929, 0.102621, 0.137563,
      0.128113, 0.103627, 0.0927869, 0.102425, 0.0949018,
      0.0892192, 0.103132, 0.0932398, 0.0760799, 0.0901593 },
    { 0.0739873, 0.0657621, 0.0733051, 0.0810056, 0.0823362,
      0.080394, 0.082885, 0.0953372, 0.105809, 0.136333,
      0.142636, 0.101022, 0.0982706, 0.0851251, 0.0798855,
      0.0851645, 0.0799862, 0.087646, 0.0741014, 0.0812106 },
    { 0.0848307, 0.0738556, 0.0615958, 0.0647394, 0.0657774,
      0.0811213, 0.0891888, 0.0819757, 0.109486, 0.110607,
      0.117478, 0.104223, 0.0886376, 0.0773295, 0.0819245,
      0.0876723, 0.0725171, 0.0808546, 0.0813469, 0.095204 },
    { 0.0693136, 0.0975141, 0.0688994, 0.071388, 0.0883577,
      0.0654523, 0.0748456, 0.07202, 0.0828481, 0.115133,
      0.137037, 0.0930515, 0.0759633, 0.101175, 0.0873124,
      0.0629688, 0.0896685, 0.0692059, 0.0921279, 0.077337 },
    { 0.0686093, 0.0533264, 0.0434495, 0.0906161, 0.0764346,
      0.0809739, 0.0952387, 0.0845709, 0.0882732, 0.114535,
      0.0980442, 0.0891977, 0.0897545, 0.0827314, 0.0916396,
      0.127231, 0.0481388, 0.0897086, 0.0760457, 0.101561 },
    { 0.0740321, 0.0793183, 0.0488967, 0.0851262, 0.0764973,
      0.0579957, 0.121692, 0.0571608, 0.0703457, 0.113346,
      0.104348, 0.105721, 0.0784515, 0.0510833, 0.0689322,
      0.070756, 0.143634, 0.12002, 0.0406569, 0.0434415 },
    { 0.0325607, 0.170293, 0.0706702, 0.0635282, 0.0902266,
      0.118271, 0.0864529, 0.100342, 0.103515, 0.08514,
      0.113142, 0.0732865, 0.112173, 0.0644814, 0.0631694,
      0.0676042, 0.0594253, 0.102948, 0.0544197, 0.110599 },
    { 0.0396907, 0.0307508, 0.110389, 0.115979, 0.0810728,
      0.0655106, 0.042684, 0.0834926, 0.0970844, 0.0794593,
      0.114272, 0.0885057, 0.0870238, 0.121492, 0.0818473,
      0.108992, 0.0381046, 0.0274913, 0.0954279, 0.0635297 },
    { 0.015924, 0.0312376, 0.0517243, 0.0438412, 0.0503261,
      0.0501387, 0.105611, 0.0440046, 0.108186, 0.0893532,
      0.067842, 0.0789161, 0.0685462, 0.043805, 0.0707264,
      0.0298818, 0.10735, 0.0788188, 0.0314362, 0.0646983 }
  };

  fBartolFlavourRatioNumuBar = new TH2D("hBartolFlavourRatioNumuBar","",EnergyBins,EnergyEdges,AngleBins,AngleEdges);
  for( Int_t nx=0; nx<EnergyBins; nx++ ){
    for( Int_t ny=0; ny<AngleBins; ny++ ){
      fBartolFlavourRatioNumuBar->SetBinContent(nx+1,ny+1,tempFlavourRatioNumuBar[nx][ny]);
    }
  }

  return;
}

---