make_rhcpred_2017.C File Reference

#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>

Go to the source code of this file.

Functions
void	make_rhcpred_2017 ()

Function Documentation

void make_rhcpred_2017

(

)

Definition at line 24 of file make_rhcpred_2017.C.

References ana::SpectrumLoader::Go(), ana::kNoShift, ana::kNue2017AnaBinMergedPeripheral, ana::kNue2017FDAllSamples, ana::kPPFXFluxCVWgt, ana::kStandardSpillCuts, ana::kXSecCVWgt2017, outFile, pot, plot_validation_datamc::pred, ana::PredictionNoExtrap::SaveTo(), ana::SpectrumLoaderBase::SetSpillCut(), ana::Binning::Simple(), and string.

{
  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();
}