make_rhcpred_2017.C

Go to the documentation of this file.

#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/SpectrumLoader.h"

#include "CAFAna/Cuts/Cuts.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "3FlavorAna/Cuts/NueCuts2017.h"

#include "CAFAna/Vars/Vars.h"
#include "3FlavorAna/Vars/NumuVars.h"
#include "3FlavorAna/Vars/NueVars.h"
#include "CAFAna/Vars/PPFXWeights.h"
#include "CAFAna/Vars/GenieWeights.h"

#include "CAFAna/Prediction/PredictionNoExtrap.h"

#include "TFile.h"

#include <iostream>
#include <cmath>

using namespace ana;


void make_rhcpred_2017()
{
  double pot = 9.48e20; // Pick something


  std::string nNonS = "prod_decaf_R17-03-01-prod3reco.l_fd_genie_nonswap_rhc_nova_v08_full_nue_or_numu_or_nus_contain_v1";
  std::string nSwap = "prod_decaf_R17-03-01-prod3reco.l_fd_genie_fluxswap_rhc_nova_v08_full_nue_or_numu_or_nus_contain_v1";
  std::string nTauS = "prod_decaf_R17-03-01-prod3reco.l_fd_genie_tau_rhc_nova_v08_full_nue_or_numu_or_nus_contain_v1";

  SpectrumLoader lNonS(nNonS);
  SpectrumLoader lSwap(nSwap);
  SpectrumLoader lTauS(nTauS);

  lNonS.SetSpillCut(kStandardSpillCuts);
  lSwap.SetSpillCut(kStandardSpillCuts);
  lTauS.SetSpillCut(kStandardSpillCuts);

  Var wei = kXSecCVWgt2017 * kPPFXFluxCVWgt;
  PredictionNoExtrap pred(lNonS,lSwap,lTauS,"AnaBin",Binning::Simple(32,0,32),
                          kNue2017AnaBinMergedPeripheral,kNue2017FDAllSamples,
                          kNoShift,wei);

  lNonS.Go();
  lSwap.Go();
  lTauS.Go();

  TFile* outFile = new TFile("fdrhc_nue2017.root","RECREATE");
  outFile->cd();

  pred.SaveTo(outFile, "pred_rhc");

  outFile->Close();
}