20 #include "Utilities/func/MathUtil.h" 55 std::string anaDir =
"/nova/ana/nu_e_ana/Ana2018/Predictions/cosmic/cosmic_prediction_real_data.root";
67 MinuitFitter::FitOpts prec[4] = {MinuitFitter::FitOpts::kFast, MinuitFitter::FitOpts::kNormal, MinuitFitter::FitOpts::kCareful, MinuitFitter::FitOpts::kGradDesc};
73 clock_t begin = clock();
75 chi = fitter.Fit(calc, bestSysts, {}, {})->EvalMetricVal();
76 clock_t
end = clock();
77 counter = (double)(end - begin)/CLOCKS_PER_SEC;
Cuts and Vars for the 2020 FD DiF Study.
void ResetOscCalcToDefault(osc::IOscCalcAdjustable *calc)
Reset calculator to default assumptions for all parameters.
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.
Loads shifted spectra from files.
const FitDmSq32 kFitDmSq32
Simple record of shifts applied to systematic parameters.
void SetValue(osc::IOscCalcAdjustable *osc, double val) const override
virtual T GetTh23() const
void SetValue(osc::IOscCalcAdjustable *osc, double val) const override
Forward to wrapped Var's SetValue()
osc::IOscCalcAdjustable * DefaultOscCalc()
Create a new calculator with default assumptions for all parameters.
Representation of a spectrum in any variable, with associated POT.
const XML_Char const XML_Char * data
Sum up livetimes from individual cosmic triggers.
void test_nue2018_fitter()
void SetValue(osc::IOscCalcAdjustable *osc, double val) const override
fvar< T > exp(const fvar< T > &x)
virtual T GetDmsq32() const
std::vector< const ISyst * > getAna2018LargeXsecSysts(const EAnaType2018 ana)
const ConstrainedFitVarWithPrior fitSsqTh23_UniformPriorSsqTh23 & kFitSinSqTheta23
const FitVarWithPrior fitDeltaInPiUnits_UniformPriordCP & kFitDeltaInPiUnits
Compare a single data spectrum to the MC + cosmics expectation.
Perform MINUIT fits in one or two dimensions.