test_saloaders.C
Go to the documentation of this file.
2 #include "CAFAna/Analysis/SALoaders.h"
5 #include "CAFAna/Vars/Vars.h"
6 using namespace ana;
7 
8 #include "OscLib/OscCalcDumb.h"
9 
10 #include "TPad.h"
11 
13 {
14  SANueDecafLoaders loaders;
15 
16  HistAxis axis("LEM", Binning::Simple(50, 0, 1), kLEM);
17 
18  PredictionNoExtrap pred(loaders, axis, kNoCut);
19 
20  loaders.Go();
21 
23 
24  DataMCComparisonComponents(pred.Predict(&calc).MockData(6e20), &pred, &calc);
25  pred.PredictComponent(&calc, Flavors::kNuMuToNuE, Current::kCC, Sign::kBoth).ToTH1(6e20, kNueSignalColor)->Draw("hist same");
26 
27  gPad->SetLogy();
28 }
Spectrum PredictComponent(osc::IOscCalc *calc, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const override
const Var kLEM
PID
Definition: Vars.cxx:26
Cuts and Vars for the 2020 FD DiF Study.
Definition: vars.h:6
(' appearance')
Definition: IPrediction.h:18
TH1D * ToTH1(double exposure, Color_t col=kBlack, Style_t style=kSolid, EExposureType expotype=kPOT, EBinType bintype=kBinContent) const
Histogram made from this Spectrum, scaled to some exposure.
Definition: Spectrum.cxx:176
osc::OscCalcDumb calc
Spectrum MockData(double pot, int seed=0) const
Mock data is FakeData with Poisson fluctuations applied.
Definition: Spectrum.cxx:328
const Color_t kNueSignalColor
Definition: Style.h:19
Charged-current interactions.
Definition: IPrediction.h:39
Spectrum Predict(osc::IOscCalc *calc) const override
Simple oscillation probability calculator that has no solar term or mass hierarchy or delta so it's s...
Definition: OscCalcDumb.h:16
TH1 * DataMCComparisonComponents(const Spectrum &data, const IPrediction *mc, osc::IOscCalc *calc)
Plot MC broken down into flavour components, overlayed with data.
Definition: Plots.cxx:114
std::vector< Loaders * > loaders
Definition: syst_header.h:386
void test_saloaders()
Both neutrinos and antineutrinos.
Definition: IPrediction.h:52
Prediction that just uses FD MC, with no extrapolation.
static Binning Simple(int n, double lo, double hi, const std::vector< std::string > &labels={})
Definition: Binning.cxx:107
const Cut kNoCut
Definition: Cut.h:13