loon: BeamMatrix.cxx Source File

00001 #include "NCUtils/Extrapolation/BeamMatrix.h"
00002 
00003 //#include "NCUtils/Extrapolation/NCBeam.h" // for binning scheme
00004 
00005 //#include "TCanvas.h"
00006 //#include "TLatex.h"
00007 //#include "TMath.h"
00008 //#include "TPad.h"
00009 //#include "TStyle.h"
00010 //#include "TLegend.h"
00011 #include "MessageService/MsgService.h"
00012 
00013 //#include <cmath>
00014 //#include <iostream>
00015 //#include <algorithm>
00016 #include "TH1D.h"
00017 #include "TH2D.h"
00018 
00019 using namespace std;
00020 
00021 ClassImp(BeamMatrix)
00022 
00023   const int kNumTrueEnergyBins = 1000;
00024 //const int kNumTrueEnergyBins = 910;
00025 //const int kNumTrueEnergyBins = 165;
00026 
00027 BeamMatrix::BeamMatrix(string bmname, string title, string name)
00028   : TNamed(bmname, title)
00029 {
00030 
00031   int numBins = kNumEnergyBinsFar;
00032   double bins[kNumEnergyBinsFar+1];
00033 
00034   //use NCType to fill array of bins as it has the right binning already
00035   for( int i = 0; i < numBins+1; ++i ){
00036     bins[i] = kEnergyBinsFar[i];
00037   }
00038 
00039   // Create a true energy binning scheme for the
00040   // F/N matrix
00041 
00042   //fTrue_bin_width = bins[numBins]/double(kNumTrueEnergyBins);
00043   //vector<double> true_bins(kNumTrueEnergyBins+1, 0);
00044   
00045   //for( int i = 0; i < kNumTrueEnergyBins+1; ++i )
00046   //true_bins[i] = i*fTrue_bin_width;
00047   
00048   std::vector<double> true_bins;  
00049   /*
00050   // making new far bin scheme
00051   //165
00052   for( int i = 0; i < 100; ++i )
00053   true_bins.push_back(i*0.01);
00054 
00055   for( int i = 10; i < 61; ++i )
00056   true_bins.push_back(i*0.1);
00057 
00058   // including 20to120bin
00059   for( int i = 7; i < 21; ++i )
00060   true_bins.push_back(i*1);
00061   
00062   true_bins.push_back(120) ;
00063   */
00064   
00065   // 910 
00066   //for( int i = 0; i < 100; ++i )
00067   //true_bins.push_back(i*0.01);
00068   
00069   //for( int i = 40; i < 801; ++i )
00070   //true_bins.push_back(i*0.025);
00071   
00072   //for( int i = 22; i < 120; i+=2 )
00073   //true_bins.push_back(i);
00074   
00075   //1000
00076   for( int i = 0; i < 200; ++i )
00077     true_bins.push_back(i*0.01);
00078   for( int i = 80; i < 801; ++i )
00079     true_bins.push_back(i*0.025);
00080   for( int i = 41; i < 81; ++i )
00081     true_bins.push_back(i*0.5);
00082   for( int i = 42; i < 120; i+=2 )
00083     true_bins.push_back(i);
00084   true_bins.push_back(120) ;
00085   
00086 
00087   // including 20to120bin
00088   //for( int i = 7; i < 21; ++i )
00089   // true_bins.push_back(i*1);
00090 
00091   true_bins.push_back(120) ;
00092 
00093 
00094 
00095 
00096 
00097 
00098   fBMHistsMC.push_back(new TH1D( ("NCRecoEnergy_ND"+name).c_str() , ("NCRecoEnergy_ND"+name).c_str() , numBins, bins) );
00099   fBMHistsMC.push_back(new TH1D( ("CCRecoEnergy_ND"+name).c_str() , ("CCRecoEnergy_ND"+name).c_str() , numBins, bins) );
00100   fBMHistsMC.push_back(new TH1D( ("NCRecoBkg_ND"+name).c_str() , ("NCRecoBkg_ND"+name).c_str() , numBins, bins) );
00101   fBMHistsMC.push_back(new TH1D( ("CCRecoBkg_ND"+name).c_str() , ("CCRecoBkg_ND"+name).c_str() , numBins, bins) );
00102   fBMHistsMC.push_back(new TH1D( ("NCPurity_ND"+name).c_str() , ("NCPurity_ND"+name).c_str() , numBins, bins) );
00103   fBMHistsMC.push_back(new TH1D( ("CCPurity_ND"+name).c_str() , ("CCPurity_ND"+name).c_str() , numBins, bins) );
00104   fBMHistsMC.push_back(new TH1D( ("NCRecoEnergyPurCorr_ND"+name).c_str() , ("NCRecoEnergyPurCorr_ND"+name).c_str() , numBins, bins) );
00105   fBMHistsMC.push_back(new TH1D( ("CCRecoEnergyPurCorr_ND"+name).c_str() , ("CCRecoEnergyPurCorr_ND"+name).c_str() , numBins, bins) );
00106   fBMHistsMC.push_back(new TH1D( ("NCTrueEnergyPurCorr_ND"+name).c_str() , ("NCTrueEnergyPurCorr_ND"+name).c_str() , numBins, bins) );
00107   fBMHistsMC.push_back(new TH1D( ("CCTrueEnergyPurCorr_ND"+name).c_str() , ("CCTrueEnergyPurCorr_ND"+name).c_str() , numBins, bins) );
00108   fBMHistsMC.push_back(new TH1D( ("NCTrueSelEfficiency_ND"+name).c_str() , ("NCTrueSelEfficiency_ND"+name).c_str() , numBins, bins) );
00109   fBMHistsMC.push_back(new TH1D( ("CCTrueSelEfficiency_ND"+name).c_str() , ("CCTrueSelEfficiency_ND"+name).c_str() , numBins, bins) );
00110   fBMHistsMC.push_back(new TH1D( ("NCTrueEnergySelEffCorr_ND"+name).c_str() , ("NCTrueEnergySelEffCorr_ND"+name).c_str() , numBins, bins) );
00111   fBMHistsMC.push_back(new TH1D( ("CCTrueEnergySelEffCorr_ND"+name).c_str() , ("CCTrueEnergySelEffCorr_ND"+name).c_str() , numBins, bins) );
00112   fBMHistsMC.push_back(new TH1D( ("NCTrueRecoEfficiency_ND"+name).c_str() , ("NCTrueRecoEfficiencyCorr_ND"+name).c_str() , numBins, bins) );
00113   fBMHistsMC.push_back(new TH1D( ("CCTrueRecoEfficiency_ND"+name).c_str() , ("CCTrueRecoEfficiencyCorr_ND"+name).c_str() , numBins, bins) );
00114   fBMHistsMC.push_back(new TH1D( ("NCTrueEnergyEffCorr_ND"+name).c_str() , ("NCTrueEnergyEffCorr_ND"+name).c_str() , numBins, bins) );
00115   fBMHistsMC.push_back(new TH1D( ("CCTrueEnergyEffCorr_ND"+name).c_str() , ("CCTrueEnergyEffCorr_ND"+name).c_str() , numBins, bins) );
00116   fBMHistsMC.push_back(new TH1D( ("NCTrueEnergy_ND"+name).c_str() , ("NCTrueEnergy_ND"+name).c_str() , numBins, bins) );
00117   fBMHistsMC.push_back(new TH1D( ("CCTrueEnergy_ND"+name).c_str() , ("CCTrueEnergy_ND"+name).c_str() , numBins, bins) );
00118   fBMHistsMC.push_back(new TH1D( ("NCTrueEnergyBkg_ND"+name).c_str() , ("NCTrueEnergyBkg_ND"+name).c_str() , numBins, bins) );
00119   fBMHistsMC.push_back(new TH1D( ("CCTrueEnergyBkg_ND"+name).c_str() , ("CCTrueEnergyBkg_ND"+name).c_str() , numBins, bins) );
00120   fBMHistsMC.push_back(new TH1D( ("NCTrueEnergyFlux_ND"+name).c_str() , ("NCTrueEnergyFlux_ND"+name).c_str() , numBins, bins) );
00121   fBMHistsMC.push_back(new TH1D( ("CCTrueEnergyFlux_ND"+name).c_str() , ("CCTrueEnergyFlux_ND"+name).c_str() , numBins, bins) );
00122   fBMHistsMC.push_back(new TH1D( ("NCTrueEnergyMatrix_FD"+name).c_str() , ("NCTrueEnergyMatrix_FD"+name).c_str() , numBins, bins) );
00123   fBMHistsMC.push_back(new TH1D( ("CCTrueEnergyMatrix_FD"+name).c_str() , ("CCTrueEnergyMatrix_FD"+name).c_str() , numBins, bins) );
00124   fBMHistsMC.push_back(new TH1D( ("NCTrueEnergyFlux_FD"+name).c_str() , ("NCTrueEnergyFlux_FD"+name).c_str() , numBins, bins) );
00125   fBMHistsMC.push_back(new TH1D( ("CCTrueEnergyFlux_FD"+name).c_str() , ("CCTrueEnergyFlux_FD"+name).c_str() , numBins, bins) );
00126   fBMHistsMC.push_back(new TH1D( ("NCTrueSelEfficiency_FD"+name).c_str() , ("NCTrueSelEfficiency_FD"+name).c_str() , numBins, bins) );
00127   fBMHistsMC.push_back(new TH1D( ("CCTrueSelEfficiency_FD"+name).c_str() , ("CCTrueSelEfficiency_FD"+name).c_str() , numBins, bins) );
00128   fBMHistsMC.push_back(new TH1D( ("NCTrueEnergySelEffCorr_FD"+name).c_str() , ("NCTrueEnergySelEffCorr_FD"+name).c_str() , numBins, bins) );
00129   fBMHistsMC.push_back(new TH1D( ("CCTrueEnergySelEffCorr_FD"+name).c_str() , ("CCTrueEnergySelEffCorr_FD"+name).c_str() , numBins, bins) );
00130   fBMHistsMC.push_back(new TH1D( ("NCTrueRecoEfficiency_FD"+name).c_str() , ("NCTrueRecoEfficiency_FD"+name).c_str() , numBins, bins) );
00131   fBMHistsMC.push_back(new TH1D( ("CCTrueRecoEfficiency_FD"+name).c_str() , ("CCTrueRecoEfficiency_FD"+name).c_str() , numBins, bins) );
00132   fBMHistsMC.push_back(new TH1D( ("NCTrueEnergyEffCorr_FD"+name).c_str() , ("NCTrueEnergyEffCorr_FD"+name).c_str() , numBins, bins) );
00133   fBMHistsMC.push_back(new TH1D( ("CCTrueEnergyEffCorr_FD"+name).c_str() , ("CCTrueEnergyEffCorr_FD"+name).c_str() , numBins, bins) );
00134   fBMHistsMC.push_back(new TH1D( ("NCTrueEnergy_FD"+name).c_str() , ("NCTrueEnergy_FD"+name).c_str() , numBins, bins) );
00135   fBMHistsMC.push_back(new TH1D( ("CCTrueEnergy_FD"+name).c_str() , ("CCTrueEnergy_FD"+name).c_str() , numBins, bins) );
00136   fBMHistsMC.push_back(new TH1D( ("NCTrueEnergyXsecFit_FD"+name).c_str() , ("NCTrueEnergyXsecFit_FD"+name).c_str() , numBins, bins) );
00137   fBMHistsMC.push_back(new TH1D( ("CCTrueEnergyXsecFit_FD"+name).c_str() , ("CCTrueEnergyXsecFit_FD"+name).c_str() , numBins, bins) );
00138   fBMHistsMC.push_back(new TH1D( ("NCTrueEnergyBkg_FD"+name).c_str() , ("NCTrueEnergyBkg_FD"+name).c_str() , numBins, bins) );
00139   fBMHistsMC.push_back(new TH1D( ("CCTrueEnergyBkg_FD"+name).c_str() , ("CCTrueEnergyBkg_FD"+name).c_str() , numBins, bins) );
00140   fBMHistsMC.push_back(new TH1D( ("NCRecoEnergy_FD"+name).c_str() , ("NCRecoEnergy_FD"+name).c_str() , numBins, bins) );
00141   fBMHistsMC.push_back(new TH1D( ("CCRecoEnergy_FD"+name).c_str() , ("CCRecoEnergy_FD"+name).c_str() , numBins, bins) );
00142   fBMHistsMC.push_back(new TH1D( ("NCRecoEnergyPurCorr_FD"+name).c_str() , ("NCRecoEnergyPurCorr_FD"+name).c_str() , numBins, bins) );
00143   fBMHistsMC.push_back(new TH1D( ("CCRecoEnergyPurCorr_FD"+name).c_str() , ("CCRecoEnergyPurCorr_FD"+name).c_str() , numBins, bins) );
00144   fBMHistsMC.push_back(new TH1D( ("NCPurity_FD"+name).c_str() , ("NCPurity_FD"+name).c_str() , numBins, bins) );
00145   fBMHistsMC.push_back(new TH1D( ("CCPurity_FD"+name).c_str() , ("CCPurity_FD"+name).c_str() , numBins, bins) );
00146   fBMHistsMC.push_back(new TH1D( ("NCRecoBkg_FD"+name).c_str() , ("NCRecoBkg_FD"+name).c_str() , numBins, bins) );
00147   fBMHistsMC.push_back(new TH1D( ("CCRecoBkg_FD"+name).c_str() , ("CCRecoBkg_FD"+name).c_str() , numBins, bins) );
00148   fBMHistsMC.push_back(new TH1D( ("NCRecoEnergyPred_FD"+name).c_str() , ("NCRecoEnergyPred_FD"+name).c_str() , numBins, bins) );
00149   fBMHistsMC.push_back(new TH1D( ("CCRecoEnergyPred_FD"+name).c_str() , ("CCRecoEnergyPRed_FD"+name).c_str() , numBins, bins) );
00150   fBMHistsMC.push_back(new TH1D( ("NCRecoEnergyBeamNueBkg_FD"+name).c_str() , ("NCRecoEnergyBeamNueBkg_FD"+name).c_str() , numBins, bins) );
00151   fBMHistsMC.push_back(new TH1D( ("CCRecoEnergyBeamNueBkg_FD"+name).c_str() , ("CCRecoEnergyBeamNueBkg_FD"+name).c_str() , numBins, bins) );
00152   fBMHistsMC.push_back(new TH1D( ("NCRecoEnergyBeamNueBkg_ND"+name).c_str() , ("NCRecoEnergyBeamNueBkg_ND"+name).c_str() , numBins, bins) );
00153   fBMHistsMC.push_back(new TH1D( ("CCRecoEnergyBeamNueBkg_ND"+name).c_str() , ("CCRecoEnergyBeamNueBkg_ND"+name).c_str() , numBins, bins) );
00154   fBMHistsMC.push_back(new TH1D( ("NCRecoEnergyTau_ND"+name).c_str() , ("NCRecoEnergyTau_ND"+name).c_str() , numBins, bins) );
00155   fBMHistsMC.push_back(new TH1D( ("CCRecoEnergyTau_ND"+name).c_str() , ("CCRecoEnergyTau_ND"+name).c_str() , numBins, bins) );
00156   fBMHistsMC.push_back(new TH1D( ("NCRecoEnergyOscNue_ND"+name).c_str() , ("NCRecoEnergyOscNue_ND"+name).c_str() , numBins, bins) );
00157   fBMHistsMC.push_back(new TH1D( ("CCRecoEnergyOscNue_ND"+name).c_str() , ("CCRecoEnergyOscNue_ND"+name).c_str() , numBins, bins) );
00158   fBMHistsMC.push_back(new TH1D( ("NCSelRecoEnergy_ND"+name).c_str() , ("NCSelRecoEnergy_ND"+name).c_str() , numBins, bins) );
00159   fBMHistsMC.push_back(new TH1D( ("CCSelRecoEnergy_ND"+name).c_str() , ("CCSelRecoEnergy_ND"+name).c_str() , numBins, bins) );
00160   fBMHistsMC.push_back(new TH1D( ("NCSelRecoEnergy_FD"+name).c_str() , ("NCSelRecoEnergy_FD"+name).c_str() , numBins, bins) );
00161   fBMHistsMC.push_back(new TH1D( ("CCSelRecoEnergy_FD"+name).c_str() , ("CCSelRecoEnergy_FD"+name).c_str() , numBins, bins) );
00162   //fBMHistsMC.push_back(new TH1D( ("NCRecoEnergySig_FD"+name).c_str() , ("NCRecoEnergySig_FD"+name).c_str() , numBins, bins) );
00163   //fBMHistsMC.push_back(new TH1D( ("CCRecoEnergySig_FD"+name).c_str() , ("CCRecoEnergySig_FD"+name).c_str() , numBins, bins) );
00164   //fBMHistsMC.push_back(new TH1D( ("NCRecoEnergyWrongSignBkg_FD"+name).c_str() , ("NCRecoEnergyWrongSignBkg_FD"+name).c_str() , numBins, bins) );
00165   //fBMHistsMC.push_back(new TH1D( ("CCRecoEnergyWrongSignBkg_FD"+name).c_str() , ("CCRecoEnergyWrongSignBkg_FD"+name).c_str() , numBins, bins) );
00166   //fBMHistsMC.push_back(new TH1D( ("NCTrueEnergyMatrixNDCorr_FD"+name).c_str() , ("NCTrueEnergyMatrixNDCorr_FD"+name).c_str() , numBins, bins) );
00167   //fBMHistsMC.push_back(new TH1D( ("CCTrueEnergyMatrixNDCorr_FD"+name).c_str() , ("CCTrueEnergyMatrixNDCorr_FD"+name).c_str() , numBins, bins) );
00168 
00169 
00170   fBMHistsData.push_back(new TH1D( ("NCDataRecoEnergy_ND"+name).c_str() , ("NCDataRecoEnergy_ND"+name).c_str() , numBins, bins) );
00171   fBMHistsData.push_back(new TH1D( ("CCDataRecoEnergy_ND"+name).c_str() , ("CCDataRecoEnergy_ND"+name).c_str() , numBins, bins) );
00172   fBMHistsData.push_back(new TH1D( ("NCDataRecoEnergyPurCorr_ND"+name).c_str() , ("NCDataRecoEnergyPurCorr_ND"+name).c_str() , numBins, bins) );
00173   fBMHistsData.push_back(new TH1D( ("CCDataRecoEnergyPurCorr_ND"+name).c_str() , ("CCDataRecoEnergyPurCorr_ND"+name).c_str() , numBins, bins) );
00174   fBMHistsData.push_back(new TH1D( ("NCDataTrueEnergyPurCorr_ND"+name).c_str() , ("NCDataTrueEnergyPurCorr_ND"+name).c_str() , numBins, bins) );
00175   fBMHistsData.push_back(new TH1D( ("CCDataTrueEnergyPurCorr_ND"+name).c_str() , ("CCDataTrueEnergyPurCorr_ND"+name).c_str() , numBins, bins) );
00176   fBMHistsData.push_back(new TH1D( ("NCDataTrueEnergySelEffCorr_ND"+name).c_str() , ("NCDataTrueEnergySelEffCorr_ND"+name).c_str() , numBins, bins) );
00177   fBMHistsData.push_back(new TH1D( ("CCDataTrueEnergySelEffCorr_ND"+name).c_str() , ("CCDataTrueEnergySelEffCorr_ND"+name).c_str() , numBins, bins) );
00178   fBMHistsData.push_back(new TH1D( ("NCDataTrueEnergyEffCorr_ND"+name).c_str() , ("NCDataTrueEnergyEffCorr_ND"+name).c_str() , numBins, bins) );
00179   fBMHistsData.push_back(new TH1D( ("CCDataTrueEnergyEffCorr_ND"+name).c_str() , ("CCDataTrueEnergyEffCorr_ND"+name).c_str() , numBins, bins) );
00180   fBMHistsData.push_back(new TH1D( ("NCDataTrueEnergyFlux_ND"+name).c_str() , ("NCDataTrueEnergyFlux_ND"+name).c_str() , numBins, bins) );
00181   fBMHistsData.push_back(new TH1D( ("CCDataTrueEnergyFlux_ND"+name).c_str() , ("CCDataTrueEnergyFlux_ND"+name).c_str() , numBins, bins) );
00182   fBMHistsData.push_back(new TH1D( ("NCDataTrueEnergyMatrix_FD"+name).c_str() , ("NCDataTrueEnergyMatrix_FD"+name).c_str() , numBins, bins) );
00183   fBMHistsData.push_back(new TH1D( ("CCDataTrueEnergyMatrix_FD"+name).c_str() , ("CCDataTrueEnergyMatrix_FD"+name).c_str() , numBins, bins) );
00184   fBMHistsData.push_back(new TH1D( ("NCDataTrueEnergyFlux_FD"+name).c_str() , ("NCDataTrueEnergyFlux_FD"+name).c_str() , numBins, bins) );
00185   fBMHistsData.push_back(new TH1D( ("CCDataTrueEnergyFlux_FD"+name).c_str() , ("CCDataTrueEnergyFlux_FD"+name).c_str() , numBins, bins) );
00186   fBMHistsData.push_back(new TH1D( ("NCDataTrueEnergySelEffCorr_FD"+name).c_str() , ("NCDataTrueEnergySelEffCorr_FD"+name).c_str() , numBins, bins) );
00187   fBMHistsData.push_back(new TH1D( ("CCDataTrueEnergySelEffCorr_FD"+name).c_str() , ("CCDataTrueEnergySelEffCorr_FD"+name).c_str() , numBins, bins) );
00188   fBMHistsData.push_back(new TH1D( ("NCDataTrueEnergyEffCorr_FD"+name).c_str() , ("NCDataTrueEnergyEffCorr_FD"+name).c_str() , numBins, bins) );
00189   fBMHistsData.push_back(new TH1D( ("CCDataTrueEnergyEffCorr_FD"+name).c_str() , ("CCDataTrueEnergyEffCorr_FD"+name).c_str() , numBins, bins) );
00190   fBMHistsData.push_back(new TH1D( ("NCDataRecoEnergyPurCorr_FD"+name).c_str() , ("NCDataRecoEnergyPurCorr_FD"+name).c_str() , numBins, bins) );
00191   fBMHistsData.push_back(new TH1D( ("CCDataRecoEnergyPurCorr_FD"+name).c_str() , ("CCDataRecoEnergyPurCorr_FD"+name).c_str() , numBins, bins) );
00192   fBMHistsData.push_back(new TH1D( ("NCDataRecoEnergyPred_FD"+name).c_str() , ("NCDataRecoEnergyPred_FD"+name).c_str() , numBins, bins) );
00193   fBMHistsData.push_back(new TH1D( ("CCDataRecoEnergyPred_FD"+name).c_str() , ("CCDataRecoEnergyPred_FD"+name).c_str() , numBins, bins) );
00194   fBMHistsData.push_back(new TH1D( ("NCDataRecoEnergySig_FD"+name).c_str() , ("NCDataRecoEnergySig_FD"+name).c_str() , numBins, bins) );
00195   fBMHistsData.push_back(new TH1D( ("CCDataRecoEnergySig_FD"+name).c_str() , ("CCDataRecoEnergySig_FD"+name).c_str() , numBins, bins) );
00196   fBMHistsData.push_back(new TH1D( ("NCDataRecoEnergyWrongSignBkg_FD"+name).c_str() , ("NCDataRecoEnergyWrongSignBkg_FD"+name).c_str() , numBins, bins) );
00197   fBMHistsData.push_back(new TH1D( ("CCDataRecoEnergyWrongSignBkg_FD"+name).c_str() , ("CCDataRecoEnergyWrongSignBkg_FD"+name).c_str() , numBins, bins) );
00198   
00199  //fBMHistsData.push_back(new TH1D( ("NCDataTrueEnergyMatrixNDCorr_FD"+name).c_str() , ("NCDataTrueEnergyMatrixNDCorr_FD"+name).c_str() , numBins, bins) );
00200   //fBMHistsData.push_back(new TH1D( ("CCDataTrueEnergyMatrixNDCorr_FD"+name).c_str() , ("CCDataTrueEnergyMatrixNDCorr_FD"+name).c_str() , numBins, bins) );
00201 
00202   //kNumTrueEnergyBins, &true_bins.at(0) , numBins, bins ));
00203   fBMHists2D.push_back(new TH2D ( ("NCNomrecovstrueND"+name).c_str(), ("NC Reco vs. True (Nom) in ND"+name).c_str(),
00204                                   numBins,bins , numBins, bins ));
00205   fBMHists2D.push_back(new TH2D ( ("CCNomrecovstrueND"+name).c_str(), ("CC Reco vs. True (Nom) in ND"+name).c_str(),
00206                                   numBins,bins , numBins, bins ));
00207   fBMHists2D.push_back(new TH2D ( ("NCNomrecovstrueND_bkg"+name).c_str(), ("NC Reco vs. True (Nom)bkg in ND"+name).c_str(),
00208                                   numBins,bins , numBins, bins ));
00209   fBMHists2D.push_back(new TH2D ( ("CCNomrecovstrueND_bkg"+name).c_str(), ("CC Reco vs. True (Nom)bkg in ND"+name).c_str(),
00210                                   numBins,bins , numBins, bins ));
00211   fBMHists2D.push_back(new TH2D ( ("NCNomrecovstrueNDNorm"+name).c_str(), ("NC Reco vs. True (Nom) Normalised in ND"+name).c_str(),
00212                                   numBins,bins , numBins, bins ));
00213   fBMHists2D.push_back(new TH2D ( ("CCNomrecovstrueNDNorm"+name).c_str(), ("CC Reco vs. True (Nom) Normalised in ND"+name).c_str(),
00214                                   numBins,bins , numBins, bins ));
00215 
00216   
00217   
00218   fBMHists2D.push_back(new TH2D ( ("NCNomrecovstrueFD"+name).c_str(), ("NC Reco vs. True (Nom) in FD"+name).c_str(),
00219                                   kNumTrueEnergyBins, &true_bins.at(0) , numBins, bins ));
00220   fBMHists2D.push_back(new TH2D ( ("CCNomrecovstrueFD"+name).c_str(), ("CC Reco vs. True (Nom) in FD"+name).c_str(),
00221                                   kNumTrueEnergyBins, &true_bins.at(0) , numBins, bins ));
00222   fBMHists2D.push_back(new TH2D ( ("NCNomrecovstrueFD_bkg"+name).c_str(), ("NC Reco vs. True (Nom)bkg in FD"+name).c_str(),
00223                                   kNumTrueEnergyBins, &true_bins.at(0) , numBins, bins ));
00224   fBMHists2D.push_back(new TH2D ( ("CCNomrecovstrueFD_bkg"+name).c_str(), ("CC Reco vs. True (Nom)bkg in FD"+name).c_str(),
00225                                   kNumTrueEnergyBins, &true_bins.at(0) , numBins, bins ));
00226   fBMHists2D.push_back(new TH2D ( ("NCNomrecovstrueFDNorm"+name).c_str(), ("NC Reco vs. True (Nom) Normalised in FD"+name).c_str(),
00227                                   kNumTrueEnergyBins, &true_bins.at(0) , numBins, bins ));
00228   fBMHists2D.push_back(new TH2D ( ("CCNomrecovstrueFDNorm"+name).c_str(), ("CC Reco vs. True (Nom) Normalised in FD"+name).c_str(),
00229                                   kNumTrueEnergyBins, &true_bins.at(0) , numBins, bins ));
00230 
00231 
00232 
00233 
00234 
00235 
00236 
00237 
00238 
00239 
00240   fBMHists2DOsc.push_back(new TH2D ( ("NCNomtruevsreco"+name).c_str(), ("NC True vs. Reco (Nom)"+name).c_str(),
00241                                      numBins, bins, kNumTrueEnergyBins, &true_bins.at(0) ));
00242   fBMHists2DOsc.push_back(new TH2D ( ("CCNomtruevsreco"+name).c_str(), ("CC True vs. Reco (Nom)"+name).c_str(),
00243                                      numBins, bins, kNumTrueEnergyBins, &true_bins.at(0) ));
00244   fBMHists2DOsc.push_back(new TH2D ( ("NCNomtruevsreco_bkg"+name).c_str(), ("NC True vs. Reco (Nom)bkg"+name).c_str(),
00245                                      numBins, bins, kNumTrueEnergyBins, &true_bins.at(0) ));
00246   fBMHists2DOsc.push_back(new TH2D ( ("CCNomtruevsreco_bkg"+name).c_str(), ("CC True vs. Reco (Nom)bkg"+name).c_str(),
00247                                      numBins, bins, kNumTrueEnergyBins, &true_bins.at(0) ));
00248   fBMHists2DOsc.push_back(new TH2D ( ("NCNomtruevsreco_tau"+name).c_str(), ("NC True vs. Reco (Nom)tau"+name).c_str(),
00249                                      numBins, bins, kNumTrueEnergyBins, &true_bins.at(0) ));
00250   fBMHists2DOsc.push_back(new TH2D ( ("CCNomtruevsreco_tau"+name).c_str(), ("CC True vs. Reco (Nom)tau"+name).c_str(),
00251                                      numBins, bins, kNumTrueEnergyBins, &true_bins.at(0) ));
00252   fBMHists2DOsc.push_back(new TH2D ( ("NCNomtruevsreco_beamnue"+name).c_str(), ("NC True vs. Reco (Nom)beamnue"+name).c_str(),
00253                                      numBins, bins, kNumTrueEnergyBins, &true_bins.at(0) ));
00254   fBMHists2DOsc.push_back(new TH2D ( ("CCNomtruevsreco_beamnue"+name).c_str(), ("CC True vs. Reco (Nom)beamnue"+name).c_str(),
00255                                      numBins, bins, kNumTrueEnergyBins, &true_bins.at(0) ));
00256   fBMHists2DOsc.push_back(new TH2D ( ("NCNomtruevsreco_oscnue"+name).c_str(), ("NC True vs. Reco (Nom)oscnue"+name).c_str(),
00257                                      numBins, bins, kNumTrueEnergyBins, &true_bins.at(0) ));
00258   fBMHists2DOsc.push_back(new TH2D ( ("CCNomtruevsreco_oscnue"+name).c_str(), ("CC True vs. Reco (Nom)oscnue"+name).c_str(),
00259                                      numBins, bins, kNumTrueEnergyBins, &true_bins.at(0) ));
00260   fBMHists2DOsc.push_back(new TH2D ( ("NCNomtruevsreco_wrongsignbkg"+name).c_str(), ("NC True vs. Reco (Nom)wrong sign bkg"+name).c_str(),
00261                                      numBins, bins, kNumTrueEnergyBins, &true_bins.at(0) ));
00262   fBMHists2DOsc.push_back(new TH2D ( ("CCNomtruevsreco_wrongsignbkg"+name).c_str(), ("CC True vs. Reco (Nom)wrong sign bkg"+name).c_str(),
00263                                      numBins, bins, kNumTrueEnergyBins, &true_bins.at(0) ));
00264 
00265 
00266   fFD_BeamMatrixfit.push_back( new TH1D( ("FD_BeamMatrixfit_nc"+name).c_str() , ("FD_BeamMatrixfit NC"+name).c_str(), numBins, bins ) );
00267   fFD_BeamMatrixfit.push_back( new TH1D( ("FD_BeamMatrixfit_cc"+name).c_str() , ("FD_BeamMatrixfit CC"+name).c_str(), numBins, bins ) );
00268 
00269   fFD_data.push_back( new TH1D( ("FD_data_NC"+name).c_str(), ("FD data NC"+name).c_str(), kNumEnergyBinsFar, kEnergyBinsFar) );
00270   fFD_data.push_back( new TH1D( ("FD_data_CC"+name).c_str(), ("FD data CC"+name).c_str(), kNumEnergyBinsFar, kEnergyBinsFar) );
00271 
00272   mScaleNCbm.push_back(1.0);
00273   mScaleNCbm.push_back(1.0);
00274   mScaleNCbm.push_back(1.0);
00275   mScaleNCbm.push_back(1.0);
00276 
00277   /*
00278   std::vector< TH1D* > fNDByEnergyForXsecLE; 
00279   fNDByEnergyForXsecLE.push_back(new TH1D( ("NC0to4LE"+name).c_str() , ("NC0to4LE"+name).c_str() , numBins, bins) );
00280   fNDByEnergyForXsecLE.push_back(new TH1D( ("NC5to8LE"+name).c_str() , ("NC5to8LE"+name).c_str() , numBins, bins) );
00281   fNDByEnergyForXsecLE.push_back(new TH1D( ("NC8to15LE"+name).c_str() , ("NC8to15LE"+name).c_str() , numBins, bins) );
00282   fNDByEnergyForXsecLE.push_back(new TH1D( ("NCover15LE"+name).c_str() , ("NCover15LE"+name).c_str() , numBins, bins) );
00283   fNDByEnergyForXsecLE.push_back(new TH1D( ("NCnotMuLE"+name).c_str() , ("NCnotMuLE"+name).c_str() , numBins, bins) );
00284 
00285   std::vector< TH1D* > fNDByEnergyForXsecME; 
00286   fNDByEnergyForXsecME.push_back(new TH1D( ("NC0to4ME"+name).c_str() , ("NC0to4ME"+name).c_str() , numBins, bins) );
00287   fNDByEnergyForXsecME.push_back(new TH1D( ("NC5to8ME"+name).c_str() , ("NC5to8ME"+name).c_str() , numBins, bins) );
00288   fNDByEnergyForXsecME.push_back(new TH1D( ("NC8to15ME"+name).c_str() , ("NC8to15ME"+name).c_str() , numBins, bins) );
00289   fNDByEnergyForXsecME.push_back(new TH1D( ("NCover15ME"+name).c_str() , ("NCover15ME"+name).c_str() , numBins, bins) );
00290   fNDByEnergyForXsecME.push_back(new TH1D( ("NCnotMuME"+name).c_str() , ("NCnotMuME"+name).c_str() , numBins, bins) );
00291   
00292   std::vector< TH1D* > fNDByEnergyForXsecHE; 
00293   fNDByEnergyForXsecHE.push_back(new TH1D( ("NC0to4HE"+name).c_str() , ("NC0to4HE"+name).c_str() , numBins, bins) );
00294   fNDByEnergyForXsecHE.push_back(new TH1D( ("NC5to8HE"+name).c_str() , ("NC5to8HE"+name).c_str() , numBins, bins) );
00295   fNDByEnergyForXsecHE.push_back(new TH1D( ("NC8to15HE"+name).c_str() , ("NC8to15HE"+name).c_str() , numBins, bins) );
00296   fNDByEnergyForXsecHE.push_back(new TH1D( ("NCover15HE"+name).c_str() , ("NCover15HE"+name).c_str() , numBins, bins) );
00297   fNDByEnergyForXsecHE.push_back(new TH1D( ("NCnotMuHE"+name).c_str() , ("NCnotMuHE"+name).c_str() , numBins, bins) );
00298  
00299 
00300 
00301   fND_XsecByParams.push_back(fNDByEnergyForXsecLE);
00302   fND_XsecByParams.push_back(fNDByEnergyForXsecME);
00303   fND_XsecByParams.push_back(fNDByEnergyForXsecHE);
00304 
00305 
00306   fND_XSectionFitRew.push_back(new TH1D("NCDataLE"+name).c_str() , "NCDataLE"+name).c_str() , numBins, bins) );
00307   fND_XSectionFitRew.push_back(new TH1D("NCMCLE"+name).c_str() , "NCMCLE"+name).c_str() , numBins, bins) );
00308   fND_XSectionFitRew.push_back(new TH1D("NCDataME"+name).c_str() , "NCDataME"+name).c_str() , numBins, bins) );
00309   fND_XSectionFitRew.push_back(new TH1D("NCMCME"+name).c_str() , "NCMCME"+name).c_str() , numBins, bins) );
00310   fND_XSectionFitRew.push_back(new TH1D("NCDataHE"+name).c_str() , "NCDataHE"+name).c_str() , numBins, bins) );
00311   fND_XSectionFitRew.push_back(new TH1D("NCMCHE"+name).c_str() , "NCMCHE"+name).c_str() , numBins, bins) );
00312   fND_XSectionFitRew.push_back(new TH1D("CCDataLE"+name).c_str() , "CCDataLE"+name).c_str() , numBins, bins) );
00313   fND_XSectionFitRew.push_back(new TH1D("CCMCLE"+name).c_str() , "CCMCLE"+name).c_str() , numBins, bins) );
00314   fND_XSectionFitRew.push_back(new TH1D("CCDataME"+name).c_str() , "CCDataME"+name).c_str() , numBins, bins) );
00315   fND_XSectionFitRew.push_back(new TH1D("CCMCME"+name).c_str() , "CCMCME"+name).c_str() , numBins, bins) );
00316   fND_XSectionFitRew.push_back(new TH1D("CCDataHE"+name).c_str() , "CCDataHE"+name).c_str() , numBins, bins) );
00317   fND_XSectionFitRew.push_back(new TH1D("CCMCHE"+name).c_str() , "CCMCHE"+name).c_str() , numBins, bins) );
00318 
00319 
00320   fND_XSectionFit.push_back(new TH1D( ("NCDataLEOrig"+name).c_str() , ("NCDataLEOrig"+name).c_str() , numBins, bins) );
00321   fND_XSectionFit.push_back(new TH1D( ("NCMCLEOrig"+name).c_str() , ("NCMCLEOrig"+name).c_str() , numBins, bins) );
00322   fND_XSectionFit.push_back(new TH1D( ("NCDataMEOrig"+name).c_str() , ("NCDataMEOrig"+name).c_str() , numBins, bins) );
00323   fND_XSectionFit.push_back(new TH1D( ("NCMCMEOrig"+name).c_str() , ("NCMCMEOrig"+name).c_str() , numBins, bins) );
00324   fND_XSectionFit.push_back(new TH1D( ("NCDataHEOrig"+name).c_str() , ("NCDataHEOrig"+name).c_str() , numBins, bins) );
00325   fND_XSectionFit.push_back(new TH1D( ("NCMCHEOrig"+name).c_str() , ("NCMCHEOrig"+name).c_str() , numBins, bins) );
00326   fND_XSectionFit.push_back(new TH1D( ("CCDataLEOrig"+name).c_str() , ("CCDataLEOrig"+name).c_str() , numBins, bins) );
00327   fND_XSectionFit.push_back(new TH1D( ("CCMCLEOrig"+name).c_str() , ("CCMCLEOrig"+name).c_str() , numBins, bins) );
00328   fND_XSectionFit.push_back(new TH1D( ("CCDataMEOrig"+name).c_str() , ("CCDataMEOrig"+name).c_str() , numBins, bins) );
00329   fND_XSectionFit.push_back(new TH1D( ("CCMCMEOrig"+name).c_str() , ("CCMCMEOrig"+name).c_str() , numBins, bins) );
00330   fND_XSectionFit.push_back(new TH1D( ("CCDataHEOrig"+name).c_str() , ("CCDataHEOrig"+name).c_str() , numBins, bins) );
00331   fND_XSectionFit.push_back(new TH1D( ("CCMCHEOrig"+name).c_str() , ("CCMCHEOrig"+name).c_str() , numBins, bins) );
00332 
00333   */ 
00334   //fFD_data[0] = new TH1D( ("FD_data_NC"+name).c_str(), ("FD data NC"+name).c_str(), kNumEnergyBinsFar, kEnergyBinsFar);
00335   //fFD_data[1] = new TH1D( ("FD_data_CC"+name).c_str(), ("FD data CC"+name).c_str(), kNumEnergyBinsFar, kEnergyBinsFar);
00336   
00337   /*
00338   std::vector<TH1D*>  readInHists;
00339   readinHists.push_back(selEffNDHist);
00340   readinHists.push_back(selCCEfFDHist);
00341   readinHists.push_back(selNCEffFDHist);
00342   readinHists.push_back(recoEffNDHist);
00343   readinHists.push_back(recoCCEffFDHist);
00344   readinHists.push_back(recoNCEffFDHist);
00345   readinHists.push_back(extractNCFromCCFlux);
00346   */
00347 
00348   //std::vector<TH2D*>  beamMatrix;
00349   //beamMatrix.push_back(fBeamMatrix);
00350 
00351  
00352  //how do I get 2d and 1d in here? 
00353  //std::vector<std::vector<TH1* > > fBMMatricesForABeam ;
00354  //fBMMatricesForABeam.push_back(fBMHistsMC);
00355  //fBMMatricesForABeam.push_back(fBMHistsData);
00356  //fBMMatricesForABeam.push_back(fBMHists2D);
00357  //fBMMatricesForABeam.push_back(fBMHists2DOsc);
00358  //fBMMatricesForABeam.push_back(fFD_BeamMatrixfit);
00359  //fBMMatricesForABeam.push_back(fFD_data);
00360  //fBMMatricesForABeam.push_back(readInHists);
00361  //fBMMatricesForABeam.push_back(beamMatrix);
00362  
00363  //return fBMMatricesForABeam ;
00364 
00365 }
00366 void BeamMatrix::
00367 SetmScaleNCparam(int param , double val){ 
00368   mScaleNCbm.at(param) = val ; 
00369 return;
00370 }