make_rhcpred_2017.C
Go to the documentation of this file.
1 #include "CAFAna/Core/Spectrum.h"
3 
4 #include "CAFAna/Cuts/Cuts.h"
7 
8 #include "CAFAna/Vars/Vars.h"
13 
15 
16 #include "TFile.h"
17 
18 #include <iostream>
19 #include <cmath>
20 
21 using namespace ana;
22 
23 
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 
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 }
57 
Cuts and Vars for the 2020 FD DiF Study.
Definition: vars.h:6
const Var kPPFXFluxCVWgt
Definition: PPFXWeights.h:16
void SetSpillCut(const SpillCut &cut)
TFile * outFile
Definition: PlotXSec.C:135
#define pot
virtual void Go() override
Load all the registered spectra.
virtual void SaveTo(TDirectory *dir, const std::string &name) const override
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
void make_rhcpred_2017()
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
static Binning Simple(int n, double lo, double hi, const std::vector< std::string > &labels={})
Definition: Binning.cxx:107
const Cut kNue2017FDAllSamples
Our FD selection including all samples, for making predictions, etc.
Definition: NueCuts2017.h:155
enum BeamMode string