Functions
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.

25 {
26  double pot = 9.48e20; // Pick something
27 
28 
29  std::string nNonS = "prod_decaf_R17-03-01-prod3reco.l_fd_genie_nonswap_rhc_nova_v08_full_nue_or_numu_or_nus_contain_v1";
30  std::string nSwap = "prod_decaf_R17-03-01-prod3reco.l_fd_genie_fluxswap_rhc_nova_v08_full_nue_or_numu_or_nus_contain_v1";
31  std::string nTauS = "prod_decaf_R17-03-01-prod3reco.l_fd_genie_tau_rhc_nova_v08_full_nue_or_numu_or_nus_contain_v1";
32 
33  SpectrumLoader lNonS(nNonS);
34  SpectrumLoader lSwap(nSwap);
35  SpectrumLoader lTauS(nTauS);
36 
37  lNonS.SetSpillCut(kStandardSpillCuts);
38  lSwap.SetSpillCut(kStandardSpillCuts);
39  lTauS.SetSpillCut(kStandardSpillCuts);
40 
42  PredictionNoExtrap pred(lNonS,lSwap,lTauS,"AnaBin",Binning::Simple(32,0,32),
44  kNoShift,wei);
45 
46  lNonS.Go();
47  lSwap.Go();
48  lTauS.Go();
49 
50  TFile* outFile = new TFile("fdrhc_nue2017.root","RECREATE");
51  outFile->cd();
52 
53  pred.SaveTo(outFile, "pred_rhc");
54 
55  outFile->Close();
56 }
const Var kPPFXFluxCVWgt
Definition: PPFXWeights.h:16
TFile * outFile
Definition: PlotXSec.C:135
#define pot
const SystShifts kNoShift
Definition: SystShifts.cxx:22
const Var kNue2017AnaBinMergedPeripheral([](const caf::SRProxy *sr){int anaBin=kNue2017AnaBin(sr);anaBin=anaBin > 27?29:anaBin;return anaBin;})
Definition: NueCuts2017.h:319
Collaborates with Spectrum and OscillatableSpectrum to fill spectra from CAF files.
const SpillCut kStandardSpillCuts
Apply this unless you&#39;re doing something special.
Definition: SpillCuts.h:49
Prediction that just uses FD MC, with no extrapolation.
const Var kXSecCVWgt2017
Definition: XsecTunes.h:36
const Cut kNue2017FDAllSamples
Our FD selection including all samples, for making predictions, etc.
Definition: NueCuts2017.h:155
enum BeamMode string