ana::MinuitFitter Class Reference

Perform MINUIT fits in one or two dimensions. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-04-14/CAFAna/Fit/MinuitFitter.h"

Inheritance diagram for ana::MinuitFitter:

Classes
class	MinuitFitSummary
	Concrete instance of IFitSummary for use in MinuitFitter. More...

Public Types
enum	FitOpts {   kNormal = 0, kCareful = 1, kFast = 2, kGradDesc = 3,   kPrecisionMask = 3, kPrefitSysts = 4 }
enum	Verbosity { kQuiet, kVerbose }

Public Member Functions
void	SetFitOpts (FitOpts opts)
	MinuitFitter (const IExperiment *expt, std::vector< const IFitVar * > vars, std::vector< const ISyst * > systs={}, FitOpts opts=kNormal)
	~MinuitFitter ()
virtual double	operator() (const std::vector< double > &pars) const override
	Evaluate the log-likelihood, as required by MINUT interface. More...
std::vector< double >	Gradient (const std::vector< double > &pars) const override
bool	CheckGradient () const override
TMatrixDSym	GetCovariance () const
virtual double	Up () const override
	Definition of one-sigma, required by MINUIT. More...
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. More...
virtual std::unique_ptr< IFitSummary >	Fit (osc::IOscCalcAdjustable *seed, SystShifts &systSeed, Verbosity verb) const
	Variant with no seedPts. More...
virtual std::unique_ptr< IFitSummary >	Fit (osc::IOscCalcAdjustable *seed, Verbosity verb) const
	Variant with no seedPts and no systematics result returned. More...
virtual std::unique_ptr< IFitSummary >	Fit (SystShifts &systSeed, Verbosity verb=kVerbose) const
	Variant with no oscillations - useful for ND fits. More...
std::unique_ptr< SystShifts >	GetSystShifts () const

Protected Member Functions
int	GetMINUITPrecision () const
void	DecodePars (const std::vector< double > &pars, osc::IOscCalcAdjustable *calc) const
	Stuff the parameters into the calculator and/or syst shifts object. More...
std::unique_ptr< IFitSummary >	FitHelperSeeded (osc::IOscCalcAdjustable *seed, SystShifts &systSeed, Verbosity verb) const override
	Helper for FitHelper. More...
bool	SupportsDerivatives () const
std::vector< SeedPt >	ExpandSeeds (const SeedList &seedPts, const std::vector< SystShifts > &systSeedPts) const
virtual std::unique_ptr< IFitSummary >	FitHelper (osc::IOscCalcAdjustable *seed, SystShifts &bestSysts, const SeedList &seedPts, const std::vector< SystShifts > &systSeedPts, Verbosity verb) const
	Helper for Fit. More...
void	ValidateSeeds (osc::IOscCalcAdjustable *seed, const SeedList &seedPts, const std::vector< SystShifts > &systSeedPts) const
	Check that the seeds that were specified are compatible with the vars being fitted. More...

Protected Attributes
osc::IOscCalcAdjustable *	fCalc
const IExperiment *	fExpt
FitOpts	fFitOpts
bool	fSupportsDerivatives
int	fNEval = 0
int	fNEvalGrad = 0
int	fNEvalFiniteDiff = 0
TMatrixDSym	fCovariance
std::vector< const IFitVar * >	fVars
std::vector< const ISyst * >	fSysts
std::unique_ptr< SystShifts >	fShifts
Verbosity	fVerb

Static Protected Attributes
static SystShifts	junkShifts = SystShifts()

Detailed Description

Perform MINUIT fits in one or two dimensions.

Definition at line 17 of file MinuitFitter.h.

Member Enumeration Documentation

enum ana::MinuitFitter::FitOpts

Enumerator
kNormal	default if un-specified
kCareful
kFast
kGradDesc
kPrecisionMask
kPrefitSysts

Definition at line 21 of file MinuitFitter.h.

      {
        // These first three correspond to MIGRAD options (in a different
        // order)
        kNormal = 0, ///< default if un-specified
        kCareful = 1,
        kFast = 2,
        // A simple alternate fitter, exclusive with the options above
        kGradDesc = 3,
        kPrecisionMask = 3,
        // This option is independent, you may xor it in with one of the above
        kPrefitSysts = 4
      };

enum ana::IFitter::Verbosity

inherited

Enumerator
kQuiet
kVerbose

Definition at line 22 of file IFitter.h.

{kQuiet, kVerbose};

Constructor & Destructor Documentation

ana::MinuitFitter::MinuitFitter	(	const IExperiment *	expt,
		std::vector< const IFitVar * >	vars,
		std::vector< const ISyst * >	systs = {},
		FitOpts	opts = kNormal
	)

Definition at line 64 of file MinuitFitter.cxx.

    : IFitter(vars, systs),
      fExpt(expt),
      fFitOpts(opts),
      fSupportsDerivatives(SupportsDerivatives())
  {
  }

ana::MinuitFitter::~MinuitFitter

(

)

Definition at line 76 of file MinuitFitter.cxx.

  {
  }

Member Function Documentation

bool ana::MinuitFitter::CheckGradient ( ) const

inlineoverride

Definition at line 49 of file MinuitFitter.h.

References fFitOpts, kFast, and kPrecisionMask.

{ return (fFitOpts & kPrecisionMask) != kFast; }      

void ana::MinuitFitter::DecodePars ( const std::vector< double > &  pars, osc::IOscCalcAdjustable *  calc  ) const

protected

Stuff the parameters into the calculator and/or syst shifts object.

Definition at line 284 of file MinuitFitter.cxx.

References ana::assert(), ana::IFitter::fShifts, ana::IFitter::fSysts, ana::IFitter::fVars, MECModelEnuComparisons::i, calib::j, and febshutoff_auto::val.

Referenced by GetMINUITPrecision(), and operator()().

  {
    assert(pars.size() == fVars.size()+fSysts.size());

    if (fVars.size() > 0)
    {
      assert(calc);
      for(unsigned int i = 0; i < fVars.size(); ++i){
        auto val = pars[i];
        fVars[i]->SetValue(calc, val);
      }
    }

    for(unsigned int j = 0; j < fSysts.size(); ++j){
      auto val = pars[fVars.size()+j];
      fShifts->SetShift(fSysts[j], val);
    }
  }

std::vector< IFitter::SeedPt > ana::IFitter::ExpandSeeds ( const SeedList &  seedPts, const std::vector< SystShifts > &  systSeedPts  ) const

protectedinherited

Definition at line 46 of file IFitter.cxx.

References ana::SeedList::GetSeeds(), ana::SystShifts::Nominal(), gen_hdf5record::pt, runNovaSAM::ret, and seed.

Referenced by ana::IFitter::FitHelper().

  {
    std::vector<SeedPt> ret;
    for(Seed seed: seedPts.GetSeeds()) ret.emplace_back(seed, SystShifts::Nominal());

    // Now duplicate as many times as required for the syst seeds
    if(!systSeedPts.empty()){
      std::vector<SeedPt> newret;
      for(const SystShifts& s: systSeedPts){
        for(SeedPt pt: ret){
          pt.shift = s;
          newret.push_back(pt);
        }
      }
      ret = newret;
    }

    return ret;
  }

std::unique_ptr< IFitter::IFitSummary > ana::IFitter::Fit ( osc::IOscCalcAdjustable *  seed, SystShifts &  bestSysts = junkShifts, const SeedList &  seedPts = SeedList(), const std::vector< SystShifts > &  systSeedPts = {}, Verbosity  verb = kVerbose  ) const

virtualinherited

Master fitting method. Depends on FitHelper and FitHelperSeeded.

Parameters

[out]	seed	Seed parameter and output best-fit point
[out]	bestSysts	Best systematics result returned here
	seedPts	List of oscillation parameter seeds
	systSeedPts	If non-empty, try fit starting at each of these
	verb	If quiet, no printout

Returns

-2x the log-likelihood of the best-fit point

Reimplemented in ana::StanFitter.

Definition at line 69 of file IFitter.cxx.

References om::cout, allTimeWatchdog::endl, ana::IFitter::FitHelper(), ana::IFitter::fSysts, ana::IFitter::fVars, ana::IFitter::fVerb, ana::SystShifts::GetShift(), ana::IFitter::kVerbose, registry_explorer::v, and ana::IFitter::ValidateSeeds().

Referenced by BlessedPlotsAna(), BlessedPlotsAnaByPeriod(), cc(), fake_future_data(), FCContour(), FCTutorial2020(), fill_col(), ana::FrequentistSurface::FillSurfacePoint(), ana::FrequentistSurface::FindMinimum(), ana::FindValley(), ana::IFitter::Fit(), ana::StanFitter::Fit(), FitHelperSeeded(), getBestFit(), joint_fit_2017_contours(), joint_fit_2017_make_fc_slice(), joint_fit_2017_make_fc_surf(), joint_fit_2017_slices(), joint_fit_2018_contours(), joint_fit_2018_slices(), joint_fit_future_bestfit_univ(), joint_fit_future_contour_univ(), make_fc_mass_and_oct_nersc_2018(), make_fc_mh_nersc_2018(), make_fc_nus_surfs_nersc_2018(), make_fc_nus_surfs_nersc_2019(), make_fc_oct_nersc_2018(), make_fc_slices_nersc_2018(), make_fc_slices_nersc_2018_stats(), make_fc_surfaces_2020(), make_fc_surfaces_2020_validation(), make_fc_surfaces_nersc_2018(), make_fc_surfaces_nersc_2018_stats(), make_nus17_fc_surfs(), make_nus_fc_surfs(), make_starplots(), MakeCAFSensitivities_for_FNEX(), numu_sig_nonmax(), plot_shifts(), PlotNus17Prediction(), ana::Profile(), reach_2018_dCPfractions(), ana::RefineSeeds(), run_joint_fit_2020_bestfit(), run_joint_fit_2020_contours(), run_joint_fit_2020_slices(), RunFitter(), sensitivity2018(), sensitivity2020(), starPlot(), syst_test(), and test_ana().

  {
    fVerb = verb;
    ValidateSeeds(seed, seedPts, systSeedPts);

    if (fVerb == kVerbose)
    {
      std::cout << "Finding best fit for";
      for (const auto *v: fVars) std::cout << " " << v->ShortName();
      for (const auto *s: fSysts) std::cout << " " << s->ShortName();
//      if (fSupportsDerivatives) std::cout << " using analytic derivatives";
      std::cout << "..." << std::endl;
    }

    // Do all the actual work. This wrapper function is just so we can have
    // better control of printing.
    auto fitSummary = FitHelper(seed, bestSysts, seedPts, systSeedPts, verb);

    if (fVerb == kVerbose)
    {
      std::cout << "Best fit";
      for (const auto *v: fVars)
      {
        std::cout << ", " << v->ShortName() << " = " << v->GetValue(seed);
      }
      for (const auto *s: fSysts)
      {
        std::cout << ", " << s->ShortName() << " = " << bestSysts.GetShift(s);
      }
      std::cout << ", eval metric = " << fitSummary->EvalMetricVal() << std::endl;

//      std::cout << "  found with " << fNEval << " evaluations of the likelihood";
//      if (fNEvalFiniteDiff > 0)
//        std::cout << " (" << fNEvalFiniteDiff << " to evaluate gradients numerically)";
//      if (fNEvalGrad > 0)
//        std::cout << " and " << fNEvalGrad << " evaluations of the gradient";
//      std::cout << std::endl;
    }

    return fitSummary;
  }

virtual std::unique_ptr<IFitSummary> ana::IFitter::Fit ( osc::IOscCalcAdjustable *  seed, SystShifts &  systSeed, Verbosity  verb  ) const

inlinevirtualinherited

Variant with no seedPts.

Definition at line 61 of file IFitter.h.

References ana::IFitter::Fit().

      {
        return Fit(seed, systSeed, {}, std::vector<SystShifts>(1, systSeed), verb);
      }

virtual std::unique_ptr<IFitSummary> ana::IFitter::Fit ( osc::IOscCalcAdjustable *  seed, Verbosity  verb  ) const

inlinevirtualinherited

Variant with no seedPts and no systematics result returned.

Definition at line 69 of file IFitter.h.

References ana::IFitter::Fit().

      {
        return Fit(seed, junkShifts, {}, {}, verb);
      }

virtual std::unique_ptr<IFitSummary> ana::IFitter::Fit ( SystShifts &  systSeed, Verbosity  verb = kVerbose  ) const

inlinevirtualinherited

Variant with no oscillations - useful for ND fits.

Definition at line 76 of file IFitter.h.

References ana::IFitter::Fit().

      {
        return Fit(nullptr, systSeed, {}, std::vector<SystShifts>(1, systSeed), verb);
      }

std::unique_ptr< IFitter::IFitSummary > ana::IFitter::FitHelper ( osc::IOscCalcAdjustable *  seed, SystShifts &  bestSysts, const SeedList &  seedPts, const std::vector< SystShifts > &  systSeedPts, Verbosity  verb  ) const

protectedvirtualinherited

Helper for Fit.

Definition at line 117 of file IFitter.cxx.

References ana::assert(), ana::SystShifts::Copy(), osc::_IOscCalcAdjustable< T >::Copy(), ana::IFitter::ExpandSeeds(), ana::IFitter::FitHelperSeeded(), ana::IFitter::fShifts, ana::IFitter::fSysts, ana::IFitter::fVars, MECModelEnuComparisons::i, gen_hdf5record::pt, seed, ana::SystShifts::SetShift(), and registry_explorer::v.

Referenced by ana::IFitter::Fit().

  {
    // if user passed a derived kind of SystShifts, this preserves it
    fShifts = bestSysts.Copy();
    fShifts->ResetToNominal();

    const std::vector<SeedPt> pts = ExpandSeeds(seedPts, systSeedPts);

    std::unique_ptr<IFitSummary> bestFitSummary;
    std::vector<double> bestFitPars, bestSystPars;

    for (const SeedPt &pt: pts)
    {
      osc::IOscCalcAdjustable* seed = nullptr;
      if(initseed) seed = initseed->Copy();

      pt.fitvars.ResetCalc(seed);

      // be sure to keep any derived class stuff around
      auto shift = fShifts->Copy();
      *shift = pt.shift;
      auto fitSummary = FitHelperSeeded(seed, *shift, verb);
      if (fitSummary->IsBetterThan(bestFitSummary.get()))
      {
         bestFitSummary = std::move(fitSummary);
        // Store the best fit values of all the parameters we know are being
        // varied.
        bestFitPars.clear();
        bestSystPars.clear();
        for (const auto *v: fVars)
          bestFitPars.push_back(v->GetValue(seed));
        for (const auto *s: fSysts)
          bestSystPars.push_back(shift->GetShift<double>(s));
      }

      delete seed;
    } // end for pt

    assert(bestFitSummary);
    // Stuff the results of the actual best fit back into the seeds
    for (unsigned int i = 0; i < fVars.size(); ++i)
      fVars[i]->SetValue(initseed, bestFitPars[i]);
    for (unsigned int i = 0; i < fSysts.size(); ++i)
      bestSysts.SetShift(fSysts[i], bestSystPars[i]);

    return std::move(bestFitSummary);
  }

std::unique_ptr< IFitter::IFitSummary > ana::MinuitFitter::FitHelperSeeded ( osc::IOscCalcAdjustable *  seed, SystShifts &  systSeed, Verbosity  verb  ) const

overrideprotectedvirtual

Helper for FitHelper.

Implements ana::IFitter.

Definition at line 103 of file MinuitFitter.cxx.

References ana::assert(), om::cout, allTimeWatchdog::endl, fCalc, fCovariance, fExpt, fFitOpts, ana::IFitter::Fit(), ana::IFitter::fShifts, fSupportsDerivatives, ana::IFitter::fSysts, ana::IFitter::fVars, ana::IFitter::fVerb, runNovaSAM::gErrorIgnoreLevel, GetMINUITPrecision(), ana::SystShifts::GetShift(), MECModelEnuComparisons::i, calib::j, kGradDesc, kPrecisionMask, kPrefitSysts, elec2geo::pos, ana::IExperiment::Reset(), seed, ana::SystShifts::SetShift(), registry_explorer::v, and febshutoff_auto::val.

  {
    fVerb = verb;
    if((fFitOpts & kPrefitSysts) && !fVars.empty()){
      // Update systSeed to the best syst values, while holding the oscillation
      // parameters fixed
      MinuitFitter prefit(fExpt, {}, fSysts, fFitOpts);
      prefit.Fit(seed, systSeed, verb);
      // Then continue with the full fit as usual
    }

    fCalc = seed;
    *fShifts = systSeed;
    fExpt->Reset();

    ROOT::Minuit2::MnUserParameters mnPars;

    for (const IFitVar *v: fVars)
    {
      const double val = v->GetValue(seed);
      // name, value, error
      mnPars.Add(v->ShortName(), val, val ? val / 2 : .1);
    }
    // One way this can go wrong is if two variables have the same ShortName
    assert(mnPars.Params().size() == fVars.size());
    for (const ISyst *s: fSysts)
    {
      const double val = systSeed.GetShift(s);
      // name, value, error
      mnPars.Add(s->ShortName(), val, 1);
    }
    // One way this can go wrong is if two variables have the same ShortName
    assert(mnPars.Params().size() == fVars.size() + fSysts.size());

    ROOT::Minuit2::MnApplication *mnApp = 0;

    if ((fFitOpts & kPrecisionMask) == kGradDesc)
    {
      mnApp = new GradientDescent(*this, mnPars);
    } else
    {
      if (fSupportsDerivatives)
      {
        mnApp = new ROOT::Minuit2::MnMigrad(*this, mnPars, GetMINUITPrecision());
      } else
      {
        mnApp = new ROOT::Minuit2::MnMigrad(*((ROOT::Minuit2::FCNBase *) this), mnPars, GetMINUITPrecision());
      }
    }

    // Minuit2 doesn't give a good way to control verbosity...
    const int olderr = gErrorIgnoreLevel;
    if (fVerb <= Verbosity::kQuiet) gErrorIgnoreLevel = 1001; // Ignore warnings
    auto minpt = std::make_unique<ROOT::Minuit2::FunctionMinimum>((*mnApp)());

    gErrorIgnoreLevel = olderr;

    if(minpt->IsValid()){
      const std::vector<double> cov = minpt->UserCovariance().Data();
      const unsigned int nPars = mnPars.Params().size();

      assert(cov.size() == (nPars*(nPars+1))/2);

      fCovariance.ResizeTo(nPars, nPars);
      int pos = 0;
      for(unsigned int iRow = 0; iRow < nPars; iRow++){
        for(unsigned int iCol = 0; iCol < iRow+1; iCol++){
          fCovariance(iRow,iCol) = fCovariance(iCol, iRow) = cov[pos];
          pos++;
        }
      }
    }
    else{
      fCovariance.ResizeTo(0, 0);
      if(fVerb > Verbosity::kQuiet)
        std::cout << "*** ERROR: minimum is not valid ***" << std::endl;
    }

    const std::vector<double> minvec = minpt->UserParameters().Params();

    // Store results back to the "seed" variable
    for (unsigned int i = 0; i < fVars.size(); ++i)
    {
      fVars[i]->SetValue(seed, minvec[i]);
    }
    // Store systematic results back into "systSeed"
    for (unsigned int j = 0; j < fSysts.size(); ++j)
    {
      systSeed.SetShift(fSysts[j], minvec[fVars.size() + j]);
    }

    delete mnApp;

    return std::make_unique<MinuitFitSummary>(std::move(minpt));
  }

TMatrixDSym ana::MinuitFitter::GetCovariance ( ) const

inline

Definition at line 50 of file MinuitFitter.h.

References fCovariance.

{return fCovariance;}

int ana::MinuitFitter::GetMINUITPrecision ( ) const

inlineprotected

Convert from our enum order to MINUIT's codes (0=fast, 1=normal, 2=careful)

Definition at line 58 of file MinuitFitter.h.

References calc, DecodePars(), fFitOpts, kPrecisionMask, and pars().

Referenced by FitHelperSeeded().

{return ((fFitOpts & kPrecisionMask)+1)%3;}

std::unique_ptr<SystShifts> ana::IFitter::GetSystShifts ( ) const

inlineinherited

Definition at line 82 of file IFitter.h.

References ana::IFitter::fShifts.

Referenced by plot_shifts().

{return fShifts->Copy();}

std::vector< double > ana::MinuitFitter::Gradient ( const std::vector< double > &  pars ) const

override

Definition at line 238 of file MinuitFitter.cxx.

  {
    abort();
    /*
    ++fNEvalGrad;

    std::vector<double> ret(pars.size());

    // TODO handling of FitVars including penalty terms

    if (!fVars.empty())
    {
      // Have to use finite differences to calculate these derivatives
      const double dx = 1e-9;
      const double nom = (*this)(pars);
      ++fNEvalFiniteDiff;
      std::vector<double> parsCopy = pars;
      for (unsigned int i = 0; i < fVars.size(); ++i)
      {
        parsCopy[i] += dx;
        ret[i] = ((*this)(parsCopy) - nom) / dx;
        ++fNEvalFiniteDiff;
        parsCopy[i] = pars[i];
      }
    }

    // Get systematic parts analytically

    DecodePars(pars, fCalc); // Updates fCalc and fShifts

    std::unordered_map<const ISyst *, double> dchi;
    for (const ISyst *s: fSysts) dchi[s] = 0;
    fExpt->Derivative(fCalc, *fShifts, dchi);

    for (unsigned int i = 0; i < fSysts.size(); ++i)
    {
      // Include the derivative of the penalty terms too. TODO this is only
      // right for quadratic penalties (ie all the currently existing ones)
      ret[fVars.size() + i] = dchi[fSysts[i]] + 2 * fShifts->GetShift(fSysts[i]);
    }

    return ret;
    */
  }

double ana::MinuitFitter::operator() ( const std::vector< double > &  pars ) const

overridevirtual

Evaluate the log-likelihood, as required by MINUT interface.

Definition at line 218 of file MinuitFitter.cxx.

References ana::assert(), ana::IExperiment::ChiSq(), DecodePars(), fCalc, fExpt, fNEval, ana::IFitter::fShifts, ana::IFitter::fSysts, ana::IFitter::fVars, and MECModelEnuComparisons::i.

  {
    ++fNEval;

    assert(pars.size() == fVars.size() + fSysts.size());

    DecodePars(pars, fCalc); // Updates fCalc and fShifts

    // Have to re-fetch the FitVar values because DecodePars() will have
    // truncated values to the physical range where necessary.
    double penalty = 0;
    for (unsigned int i = 0; i < fVars.size(); ++i)
    {
      penalty += fVars[i]->Penalty(pars[i], fCalc);
    }

    return fExpt->ChiSq(fCalc, *fShifts) + penalty + fShifts->Penalty();
  }

void ana::MinuitFitter::SetFitOpts

(

FitOpts

opts

)

Definition at line 202 of file MinuitFitter.cxx.

References om::cout, allTimeWatchdog::endl, fFitOpts, fSupportsDerivatives, ana::IFitter::fSysts, kGradDesc, kPrecisionMask, and plot_validation_datamc::opts.

Referenced by ana::FrequentistSurface::FillSurfacePoint(), and ana::FrequentistSurface::FindMinimum().

  {
    if ((opts & kPrecisionMask) == kGradDesc &&
        !fSysts.empty() &&
        !fSupportsDerivatives)
    {
      std::cout
        << "Warning - not setting precision to kGradDesc, since analytic gradients are not supported by this experiment"
        << std::endl;
      return;
    }

    fFitOpts = opts;
  }

bool ana::MinuitFitter::SupportsDerivatives ( ) const

protected

Intended to be called only once (from constructor) to initialize fSupportsDerivatives

Definition at line 81 of file MinuitFitter.cxx.

  {
    return false;
    /*
    // Make completely opt-in for now
    if (getenv("CAFANA_ANALYTIC_DERIVATIVES") == 0) return false;

    // No point using derivatives for FitVars only, we do finite differences,
    // probably worse than MINUIT would.
    if (fSysts.empty()) return false;

    // Otherwise, do the minimal trial to see if the experiment will return a
    // gradient.
    std::unordered_map<const ISyst *, double> dchi = {{fSysts[0], 0}};
    osc::NoOscillations calc;
    fExpt->Derivative(&calc, SystShifts::Nominal(), dchi);
    return !dchi.empty();
    */
  }

virtual double ana::MinuitFitter::Up ( ) const

inlineoverridevirtual

Definition of one-sigma, required by MINUIT.

Definition at line 53 of file MinuitFitter.h.

{ return 1; }

void ana::IFitter::ValidateSeeds ( osc::IOscCalcAdjustable *  seed, const SeedList &  seedPts, const std::vector< SystShifts > &  systSeedPts  ) const

protectedinherited

Check that the seeds that were specified are compatible with the vars being fitted.

Definition at line 170 of file IFitter.cxx.

References ana::SeedList::ActiveFitVars(), om::cout, allTimeWatchdog::endl, ana::IFitter::fSysts, ana::IFitter::fVars, gen_hdf5record::pt, and registry_explorer::v.

Referenced by ana::IFitter::Fit().

  {
    if(!seed && !fVars.empty()){
      std::cout << "ERROR: MinuitFitter::Fit() trying to fit oscillation parameters without an oscillation calculator" << std::endl;
      abort();
    }

    for(const IFitVar* v: seedPts.ActiveFitVars()){
      if (std::find(fVars.begin(), fVars.end(), v) == fVars.end()){
        std::cout << "ERROR MinuitFitter::Fit() trying to seed '"
                  << v->ShortName()
                  << "' which is not part of the fit." << std::endl;
        abort();
      }
    }

    for(const SystShifts& pt: systSeedPts){
      for(const ISyst* s: pt.ActiveSysts()){
        if(std::find(fSysts.begin(), fSysts.end(), s) == fSysts.end()){
          std::cout << "ERROR MinuitFitter::Fit() trying to seed '"
                    << s->ShortName()
                    << "' which is not part of the fit." << std::endl;
          abort();
        }
      }
    }
  }

Member Data Documentation

osc::IOscCalcAdjustable* ana::MinuitFitter::fCalc

mutableprotected

Definition at line 88 of file MinuitFitter.h.

Referenced by FitHelperSeeded(), and operator()().

TMatrixDSym ana::MinuitFitter::fCovariance

mutableprotected

Definition at line 98 of file MinuitFitter.h.

Referenced by FitHelperSeeded(), and GetCovariance().

const IExperiment* ana::MinuitFitter::fExpt

protected

Definition at line 89 of file MinuitFitter.h.

Referenced by FitHelperSeeded(), and operator()().

FitOpts ana::MinuitFitter::fFitOpts

protected

Definition at line 91 of file MinuitFitter.h.

Referenced by CheckGradient(), FitHelperSeeded(), GetMINUITPrecision(), and SetFitOpts().

int ana::MinuitFitter::fNEval = 0

mutableprotected

Definition at line 95 of file MinuitFitter.h.

Referenced by operator()().

int ana::MinuitFitter::fNEvalFiniteDiff = 0

mutableprotected

Definition at line 97 of file MinuitFitter.h.

int ana::MinuitFitter::fNEvalGrad = 0

mutableprotected

Definition at line 96 of file MinuitFitter.h.

std::unique_ptr<SystShifts> ana::IFitter::fShifts

mutableprotectedinherited

Definition at line 116 of file IFitter.h.

Referenced by DecodePars(), ana::IFitter::FitHelper(), FitHelperSeeded(), ana::StanFitter::FitHelperSeeded(), ana::IFitter::GetSystShifts(), ana::StanFitter::log_prob(), operator()(), and ana::StanFitter::TestGradients().

bool ana::MinuitFitter::fSupportsDerivatives

protected

Definition at line 93 of file MinuitFitter.h.

Referenced by FitHelperSeeded(), and SetFitOpts().

std::vector<const ISyst*> ana::IFitter::fSysts

protectedinherited

Definition at line 115 of file IFitter.h.

Referenced by ana::StanFitter::BuildInitContext(), DecodePars(), ana::IFitter::Fit(), ana::IFitter::FitHelper(), FitHelperSeeded(), ana::StanFitter::FitHelperSeeded(), ana::StanFitter::get_param_names(), ana::StanFitter::log_prob(), operator()(), SetFitOpts(), ana::StanFitter::transform_inits(), and ana::IFitter::ValidateSeeds().

std::vector<const IFitVar*> ana::IFitter::fVars

protectedinherited

Definition at line 114 of file IFitter.h.

Referenced by ana::StanFitter::BuildInitContext(), DecodePars(), ana::IFitter::Fit(), ana::IFitter::FitHelper(), FitHelperSeeded(), ana::StanFitter::FitHelperSeeded(), ana::StanFitter::get_param_names(), ana::StanFitter::log_prob(), operator()(), ana::StanFitter::transform_inits(), ana::IFitter::ValidateSeeds(), and ana::StanFitter::write_array().

Verbosity ana::IFitter::fVerb

mutableprotectedinherited

Definition at line 117 of file IFitter.h.

Referenced by ana::IFitter::Fit(), and FitHelperSeeded().

SystShifts ana::IFitter::junkShifts = SystShifts()

staticprotectedinherited

Definition at line 19 of file IFitter.h.

---

The documentation for this class was generated from the following files:

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-04-14/CAFAna/Fit/MinuitFitter.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-04-14/CAFAna/Fit/MinuitFitter.cxx