3 #include "CAFAna/Core/HistAxis.h" 5 #include "CAFAna/Core/Binning.h" 6 #include "CAFAna/Core/Cut.h" 28 const double pot(1e20);
31 std::string haddDir(
"/nova/ana/users/tamsett/data/sam_definitions/");
32 SpectrumLoader loaderNDData(haddDir +
"hadd_prod_decaf_S15-05-22a_" 33 "nd_numi_numu_contain_goodruns.root");
35 "nd_genie_fhc_nonswap_numu_contain_s00-04_without_zero_entry_files.root");
36 SpectrumLoader loaderFDNonSwap(haddDir +
"hadd_prod_decaf_S15-05-22a_" 37 "fd_genie_fhc_nonswap_numu_contain_matchedruns.root");
38 SpectrumLoader loaderFDSwap(haddDir +
"hadd_prod_decaf_S15-05-22a_" 39 "fd_genie_fhc_fluxswap_numu_contain_matchedruns.root");
44 const HistAxis axis(
"Quasielastic Energy Estimator (GeV)",
56 loaderNDMC, loaderFDSwap, loaderFDNonSwap, loaderFDTauSwap,
57 decomp, axis, cutFD, cutND );
63 loaderFDNonSwap, loaderFDSwap,
64 axis.label, axis.bins, axis.var, cutFD );
75 TFile
file(
"modularextrap_demo_numu.root",
"RECREATE");
77 extrap.
SaveTo(& file,
"numu_extrap") ;
79 pred.
SaveTo(& file,
"numu_pred") ;
80 predFDfake.
SaveTo(& file,
"numu_pred_noextrap") ;
95 Spectrum fakeFDData( fakeFDMC.FakeData(pot) );
Cuts and Vars for the 2020 FD DiF Study.
void SetTh13(const T &th13) override
void SetL(double L) override
const FitSinSq2Theta23 kFitSinSq2Theta23
TH2 * Gaussian68Percent2D(const FrequentistSurface &s)
Up-value surface for 68% confidence in 2D in gaussian approximation.
virtual void Go() override
Load all the registered spectra.
Representation of a spectrum in any variable, with associated POT.
Log-likelihood scan across two parameters.
void DrawBestFit(Color_t color, Int_t marker=kFullCircle) const
Draw the best fit point.
Uses MC for NC and CC components, assigns remainder of data to CC.
void Draw() const
Draw the surface itself.
void DrawContour(TH2 *fc, Style_t style, Color_t color, double minchi=-1)
Optimized version of OscCalcPMNS.
virtual void Go() override
Load all the registered spectra.
void SetTh23(const T &th23) override
const ConstrainedFitVarWithPrior fitDmSq32Scaled_UniformPrior & kFitDmSq32Scaled
void SetDmsq21(const T &dmsq21) override
void SetDmsq32(const T &dmsq32) override
Collaborates with Spectrum and OscillatableSpectrum to fill spectra from CAF files.
TH2 * Gaussian90Percent2D(const FrequentistSurface &s)
Up-value surface for 90% confidence in 2D in gaussian approximation.
void SetdCP(const T &dCP) override
void SetTh12(const T &th12) override
const Var kQEE
Energy estimator for quasielastic CC events.
Dummy loader that doesn't load any files.
static Binning Simple(int n, double lo, double hi, const std::vector< std::string > &labels={})
void SetRho(double rho) override
Compare a single data spectrum to the MC + cosmics expectation.