nus17_fd_cut_tables2D.C

Go to the documentation of this file.

//  This macro is used to make the nue energy blessed plots. (see doc 11789)
//  Oscillation parameters are same as those in nue_ana_basic macro

#include "CAFAna/Cuts/Cuts.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "CAFAna/Analysis/Exposures.h"
#include "CAFAna/Core/Loaders.h"
#include "CAFAna/Core/Utilities.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Vars/Vars.h"
#include "NuXAna/Vars/NusVars.h"
#include "CAFAna/Vars/GenieWeights.h"
#include "CAFAna/Cuts/TimingCuts.h"
#include "OscLib/OscCalcDumb.h"
#include "CAFAna/Prediction/PredictionNoExtrap.h"
#include "CAFAna/Vars/PPFXWeights.h"
#include "NuXAna/macros/Nus17/nus17_fd_arrays.h"
#include "NuXAna/macros/Nus17/NusLoadProd3.h"

#include <string>

using namespace ana;


void nus17_fd_cut_tables2D(std::string outfile)
{

  std::string fFDmcDECAF_mcv1 = "/nova/ana/users/gsdavies/nus/prod3/gsdavies_ncid_decaf_nus17_20-06-17_1k_fd_genie_nonswap.root";

  std::string fFDcosmicdataDECAF_mcv1 = "/nova/ana/users/gsdavies/nus/prod3/gsdavies_ncid_decaf_nus17_20-06-17_11k_concat.root";
  Loaders loadermc;
  loadermc.SetLoaderPath(fFDmcDECAF_mcv1, caf::kFARDET, Loaders::kMC, kBeam, Loaders::kNonSwap);
  //SpectrumLoader loadermc(fFDmcCAF_test);
  SpectrumLoader loadercos(fFDcosmicdataDECAF_mcv1, kCosmic);

  loadermc.SetSpillCut(kStandardSpillCuts);
  loadercos.SetSpillCut(kStandardSpillCuts);
  
  //------------------------------------------------------------

  Spectrum* spec2D_nminus1[ncuts_nminus1][2];
  IPrediction* preds2D_nminus1[ncuts_nminus1][2];


  const Var kReweight = kPPFXFluxCVWgt*kXSecCVWgt2017;
  //------------------------------------------------------------


  for(int ic = 0; ic < ncuts_nminus1; ++ic){
    preds2D_nminus1[ic][0] = new PredictionNoExtrap(loadermc,
                                                    maxyptpaxis,
                                                    cutsarray_nminus1[ic],
                                                    kNoShift, kReweight);
    spec2D_nminus1[ic][0] = new Spectrum(loadercos, maxyptpaxis,
                                         cutsarray_nminus1[ic] && kInCosmicTimingWindow);
    preds2D_nminus1[ic][1] = new PredictionNoExtrap(loadermc,
                                                    timedistaxis,
                                                    cutsarray_nminus1[ic],
                                                    kNoShift, kReweight);
    spec2D_nminus1[ic][1] = new Spectrum(loadercos, timedistaxis,
                                         cutsarray_nminus1[ic] && kInCosmicTimingWindow);
  }

  loadercos.Go();
  loadermc.Go();
    
  std::cout<<"making output file \n";
  TFile fout(outfile.c_str(), "RECREATE");


  for(int ic = 7; ic < ncuts_nminus1; ++ic){
    char name[50];
    sprintf(name, "spec_nminus1_%s_%s", cutnames_nminus1[ic].c_str(), 
            "2Dptp");
    spec2D_nminus1[ic][0]->SaveTo(&fout, name);    
      
    sprintf(name, "pred_nminus1_%s_%s", cutnames_nminus1[ic].c_str(), 
            "2Dptp");
    preds2D_nminus1[ic][0]->SaveTo(&fout, name);

    
    sprintf(name, "spec_nminus1_%s_%s", cutnames_nminus1[ic].c_str(), 
            "2Dtime");
    spec2D_nminus1[ic][1]->SaveTo(&fout, name);    
    
    sprintf(name, "pred_nminus1_%s_%s", cutnames_nminus1[ic].c_str(), 
            "2Dtime");
    preds2D_nminus1[ic][1]->SaveTo(&fout, name);
  }  
}