24 #include "Utilities/func/MathUtil.h" 37 "Calorimetric energy", bins,
66 TH1* uoh = unosc.
ToTH1(18e20);
67 uoh->SetLineColor(
kBlue);
78 const double chi = fit.
Fit(&calc)->EvalMetricVal();
126 curve.
ToTH1()->Draw();
130 curve2.
ToTH1()->Draw();
157 hExp->SetLineColor(
kRed);
159 hExp->SetLineStyle(7);
163 hExp2->SetLineColor(
kGreen+2);
174 const double chi = fit.
Fit(&calc)->EvalMetricVal();
Far Detector at Ash River.
Represent the binning of a Spectrum's x-axis.
Cuts and Vars for the 2020 FD DiF Study.
void SetTh13(const T &th13) override
const Cut kCrudeCCSel([](const caf::SRProxy *sr){if(sr->trk.kalman.ntracks==0) return false;return sr->trk.kalman.tracks[0].len > 200;})
Extremely crude CC selection, only for use in some old demos.
void SetL(double L) override
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.
const FitDmSq32 kFitDmSq32
const FitSinSq2Theta23 kFitSinSq2Theta23
Collection of SpectrumLoaders for many configurations.
TH2 * Gaussian68Percent2D(const FrequentistSurface &s)
Up-value surface for 68% confidence in 2D in gaussian approximation.
void SetValue(osc::IOscCalcAdjustable *osc, double val) const override
virtual Spectrum PredictUnoscillated() const
void DisableLoader(caf::Det_t det, DataMC datamc, DataSource src=kBeam, SwappingConfig swap=kNonSwap)
virtual T GetTh23() const
void SetValue(osc::IOscCalcAdjustable *osc, double val) const override
T sqr(T x)
More efficient square function than pow(x,2)
Representation of a spectrum in any variable, with associated POT.
void Go()
Call Go() on all the loaders.
Log-likelihood scan across two parameters.
void DrawBestFit(Color_t color, Int_t marker=kFullCircle) const
Draw the best fit point.
const Var kCaloE
Summed calorimetric energy of all hits in slice, uncorrected.
Spectrum FakeData(double pot) const
Synonymous with AsimovData(). Retained for compatibility.
void Draw() const
Draw the surface itself.
void DrawContour(TH2 *fc, Style_t style, Color_t color, double minchi=-1)
virtual std::unique_ptr< IFitSummary > Fit(osc::IOscCalcAdjustable *seed, SystShifts &bestSysts=junkShifts, const SeedList &seedPts=SeedList(), const std::vector< SystShifts > &systSeedPts={}, Verbosity verb=kVerbose) const
Master fitting method. Depends on FitHelper and FitHelperSeeded.
TH1D * ToTH1(bool title=false) const
Optimized version of OscCalcPMNS.
const Cut kPassesPresel(PassesPreselFunc)
Does this event pass the nue preselection?
void SetTh23(const T &th23) override
virtual T GetDmsq32() const
void SetDmsq21(const T &dmsq21) override
void SetDmsq32(const T &dmsq32) override
virtual T GetTh13() const
Transition probability for any one channel as a function of energy.
TH1 * DataMCComparisonComponents(const Spectrum &data, const IPrediction *mc, osc::IOscCalc *calc)
Plot MC broken down into flavour components, overlayed with data.
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
double SimpleFOM(const Spectrum &obs, const Spectrum &unosc, double pot)
Figure-of-merit with no systematics, for binned data.
const FitSinSq2Theta13 kFitSinSq2Theta13
const FitVarWithPrior fitDeltaInPiUnits_UniformPriordCP & kFitDeltaInPiUnits
static Binning Simple(int n, double lo, double hi, const std::vector< std::string > &labels={})
TCanvas * RatioPlot(std::vector< TH1 * > topHistos, std::vector< TString > topOption, std::vector< TH1 * > bottomHistos, std::vector< TString > bottomOption, TString pidtag, bool pidaxis=false)
void SetRho(double rho) override
Compare a single data spectrum to the MC + cosmics expectation.
Perform MINUIT fits in one or two dimensions.