Log-likelihood scan across two parameters. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-08-04/CAFAna/Fit/FrequentistSurface.h"

Inheritance diagram for ana::FrequentistSurface:

Public Member Functions
	FrequentistSurface (const IExperiment expt, osc::IOscCalculatorAdjustable calc, const IFitVar xvar, int nbinsx, double xmin, double xmax, const IFitVar yvar, int nbinsy, double ymin, double ymax, const std::vector< const IFitVar * > &profVars={}, const std::vector< const ISyst * > &profSysts={}, const SeedList &seedPts=SeedList(), const std::vector< SystShifts > &systSeedPts={}, bool parallel=false, MinuitFitter::FitOpts opts=MinuitFitter::kNormal)

virtual	~FrequentistSurface ()

std::vector< TH2 * >	GetProfiledHists ()
	Maps of the values taken on by the profiled parameters. More...

void	SaveTo (TDirectory *dir) const

double	BestLikelihood () const

double	GetBestFitX () const

double	GetBestFitY () const

void	Draw () const
	Draw the surface itself. More...

void	DrawBestFit (Color_t color, Int_t marker=kFullCircle) const
	Draw the best fit point. More...

void	DrawContour (TH2 *fc, Style_t style, Color_t color, double minchi=-1)

std::vector< TGraph * >	GetGraphs (TH2 *fc, double minchi=-1)
	For expert use, custom painting of contours. More...

TH2 *	ToTH2 (double minchi=-1) const

void	SetTitle (const char *str)

Static Public Member Functions
static std::unique_ptr< FrequentistSurface >	LoadFrom (TDirectory *dir)

Protected Member Functions
	FrequentistSurface ()

void	CreateHistograms (const IFitVar xvar, int nbinsx, double xmin, double xmax, const IFitVar yvar, int nbinsy, double ymin, double ymax, const std::vector< const IFitVar * > &profVars, const std::vector< const ISyst * > &profSysts)

std::string	ProgressBarTitle (const IFitVar xvar, const IFitVar yvar, const std::vector< const IFitVar * > &profVars, const std::vector< const ISyst * > &profSysts) const

virtual void	FillSurface (const IExperiment expt, osc::IOscCalculatorAdjustable calc, const IFitVar xvar, const IFitVar yvar, const std::vector< const IFitVar * > &profVars, const std::vector< const ISyst * > &profSysts, const SeedList &seedPts, const std::vector< SystShifts > &systSeedPts)

double	FillSurfacePoint (const IExperiment expt, osc::IOscCalculatorAdjustable calc, const IFitVar xvar, double x, const IFitVar yvar, double y, const std::vector< const IFitVar * > &profVars, const std::vector< const ISyst * > &profSysts, const SeedList &seedPts, const std::vector< SystShifts > &systSeedPts)

void	FindMinimum (const IExperiment expt, osc::IOscCalculatorAdjustable calc, const IFitVar xvar, const IFitVar yvar, const std::vector< const IFitVar * > &profVars, const std::vector< const ISyst * > &profSysts, const SeedList &seedPts, const std::vector< SystShifts > &systSeedPts)

void	EnsureAxes () const

Static Protected Member Functions
static void	FillSurfObj (ISurface &surf, TDirectory *dir)

Protected Attributes
bool	fParallel

MinuitFitter::FitOpts	fFitOpts

std::vector< TH2 * >	fProfHists

double	fBestLikelihood

double	fBestFitX

double	fBestFitY

TH2F *	fHist

std::vector< double >	fSeedValues

Friends
class	NumuSurface

class	NueSurface

Detailed Description

Log-likelihood scan across two parameters.

Definition at line 26 of file FrequentistSurface.h.

Constructor & Destructor Documentation

ana::FrequentistSurface::FrequentistSurface	(	const IExperiment *	expt,
		osc::IOscCalculatorAdjustable *	calc,
		const IFitVar *	xvar,
		int	nbinsx,
		double	xmin,
		double	xmax,
		const IFitVar *	yvar,
		int	nbinsy,
		double	ymin,
		double	ymax,
		const std::vector< const IFitVar * > &	profVars = `{}`,
		const std::vector< const ISyst * > &	profSysts = `{}`,
		const SeedList &	seedPts = `SeedList()`,
		const std::vector< SystShifts > &	systSeedPts = `{}`,
		bool	parallel = `false`,
		MinuitFitter::FitOpts	opts = `MinuitFitter::kNormal`
	)

Parameters

expt	The experiment object to draw $\chi^2$ values from
calc	Values for oscillation parameters to be held fixed
xvar	Oscillation parameter to place on the x axis
nbinsx	Number of bins along x axis
xmin	Minimum value of x axis
xmax	Maximum value of x axis
nbinsy	Number of bins along y axis
ymin	Minimum value of y axis
ymax	Maximum value of y axis
profVars	Oscillation parameters to profile over
profSysts	Systematic parameters to profile over
seedPts	Try all combinations of these params as seeds
systSeedPts	Try all of these systematic combinations as seeds
parallel	Use all the cores on this machine? Be careful...

Definition at line 31 of file FrequentistSurface.cxx.

References CreateHistograms(), FillSurface(), FindMinimum(), ana::ISurface::fSeedValues, and registry_explorer::v.

     : fParallel(parallel), fFitOpts(opts)
   {
     CreateHistograms(xvar, nbinsx, xmin, xmax,
                      yvar, nbinsy, ymin, ymax,
                      profVars, profSysts);
 
     for(const IFitVar* v: profVars) fSeedValues.push_back(v->GetValue( calc));
 
     FillSurface(expt, calc, xvar, yvar, profVars, profSysts, seedPts, systSeedPts);
 
     FindMinimum(expt, calc, xvar, yvar, profVars, profSysts, seedPts, systSeedPts);
   }

ana::FrequentistSurface::~FrequentistSurface ( )

virtual

Definition at line 55 of file FrequentistSurface.cxx.

References CreateHistograms().

56 {

57 }

ana::FrequentistSurface::FrequentistSurface ( )

inlineprotected

Definition at line 66 of file FrequentistSurface.h.

References CreateHistograms(), FillSurface(), FillSurfacePoint(), FindMinimum(), ProgressBarTitle(), submit_syst::x, and submit_syst::y.

66 {};

Member Function Documentation

double ana::ISurface::BestLikelihood ( ) const

inlineinherited

Definition at line 29 of file ISurface.h.

References ana::ISurface::fBestLikelihood.

Referenced by ana::NumuSurface::DrawBestFit(), ana::NueSurface::MinChi(), and ana::NumuSurface::MinChi().

29 {return fBestLikelihood;}

ana::ISurface::fBestLikelihood

double fBestLikelihood

Definition: ISurface.h:59

void ana::FrequentistSurface::CreateHistograms	(	const IFitVar *	xvar,
		int	nbinsx,
		double	xmin,
		double	xmax,
		const IFitVar *	yvar,
		int	nbinsy,
		double	ymin,
		double	ymax,
		const std::vector< const IFitVar * > &	profVars,
		const std::vector< const ISyst * > &	profSysts
	)

protected

Definition at line 61 of file FrequentistSurface.cxx.

References ana::ExpandedHistogram(), ana::ISurface::fHist, fProfHists, MECModelEnuComparisons::i, ana::IFitVar::LatexName(), nbinsx, ProgressBarTitle(), plotROC::title, xmax, make_mec_shifts_plots::xmin, ymax, and ymin.

Referenced by FrequentistSurface(), and ~FrequentistSurface().

   {
     fHist = ExpandedHistogram(";"+xvar->LatexName()+";"+yvar->LatexName(),
                               nbinsx, xmin, xmax,
                               nbinsy, ymin, ymax);
 
     for(unsigned int i = 0; i < profVars.size()+profSysts.size(); ++i){
       std::string title;
       if(i < profVars.size())
         title = profVars[i]->LatexName();
       else
         title = profSysts[i-profVars.size()]->LatexName();
 
       fProfHists.push_back(ExpandedHistogram(title+";"+xvar->LatexName()+";"+yvar->LatexName(),
                                              nbinsx, xmin, xmax,
                                              nbinsy, ymin, ymax));
     }
 
   }

void ana::ISurface::Draw ( ) const

inherited

Draw the surface itself.

Definition at line 20 of file ISurface.cxx.

References ana::ISurface::EnsureAxes(), and ana::ISurface::fHist.

Referenced by cc(), demo5(), ana::NueSurface::Draw(), ana::NumuSurface::Draw(), ana::ISurface::GetBestFitY(), modularextrap_demo_nue(), modularextrap_demo_numu(), Plot2DSlices(), test_ana(), test_external_constraints(), test_stanfit_statsonly(), and test_stanfit_withsysts().

   {
     EnsureAxes();
 
     fHist->Draw("colz same");
 
     // colz obliterated them
     gPad->RedrawAxis();
 
     gPad->Update();
   }

void ana::ISurface::DrawBestFit	(	Color_t	color,
		Int_t	marker = `kFullCircle`
	)		const

inherited

Draw the best fit point.

Definition at line 33 of file ISurface.cxx.

References ana::ISurface::EnsureAxes(), ana::ISurface::fBestFitX, and ana::ISurface::fBestFitY.

Referenced by cc(), demo5(), demo_CPT(), ana::NumuSurface::DrawBestFit(), drawSensitivity(), ana::ISurface::GetBestFitY(), getSensitivity(), make_contours(), modularextrap_demo_nue(), modularextrap_demo_numu(), Plot2DSlices(), starPlot(), syst_test(), template_basic(), template_GENIE_systs(), template_nonGENIE_systs(), test_ana(), test_external_constraints(), test_stanfit_statsonly(), and test_stanfit_withsysts().

   {
     EnsureAxes();
 
     TMarker *mark = new TMarker(fBestFitX, fBestFitY, marker);
     mark->SetMarkerSize(1.5);
     mark->SetMarkerColor(color);
     mark->Draw();
     gPad->Update();
   }

void ana::ISurface::DrawContour	(	TH2 *	fc,
		Style_t	style,
		Color_t	color,
		double	minchi = `-1`
	)

inherited

Parameters

fc	Surface to compare against for this significance level
style	Line style for TAttLine, solid, dotted, dashed etc
color	Line color for TAttLine
minchi	$\chi^2$ of best fit to compare against. Default: best fit from this surface.

Definition at line 45 of file ISurface.cxx.

References ana::ISurface::EnsureAxes(), MECModelEnuComparisons::g, ana::ISurface::GetGraphs(), and APDGainPoints::gs.

Referenced by cc(), demo5(), Draw2DSurface(), ana::NueSurface::DrawContour(), ana::NumuSurface::DrawContour(), drawSensitivity(), ana::ISurface::GetBestFitY(), getSensitivity(), make_contours(), make_nom_expt(), modularextrap_demo_nue(), modularextrap_demo_numu(), Plot2DSlice(), Plot2DSlices(), plotContProf(), plots(), starPlot(), syst_test(), template_basic(), template_GENIE_systs(), template_nonGENIE_systs(), test_ana(), test_external_constraints(), test_stanfit_statsonly(), and test_stanfit_withsysts().

   {
     EnsureAxes();
 
     std::vector<TGraph *> gs = GetGraphs(fc, minchi);
 
     for (TGraph *g: gs)
     {
       g->SetLineWidth(3);//2);
       g->SetLineStyle(style);
       g->SetLineColor(color);
       g->Draw("l");
     }
 
     gPad->Update();
   }

void ana::ISurface::EnsureAxes ( ) const

protectedinherited

Definition at line 64 of file ISurface.cxx.

References visualisationForPaperMasterPlot::ax, file_size_ana::axes, ana::ISurface::fHist, genie::utils::style::Format(), and ana::UniqueName().

Referenced by ana::ISurface::Draw(), ana::ISurface::DrawBestFit(), ana::ISurface::DrawContour(), and ana::NueSurface::DrawValley().

   {
     // Could have a file temporarily open
     DontAddDirectory guard;
 
     // If this pad has already been drawn in, already has axes
     if (gPad && !gPad->GetListOfPrimitives()->IsEmpty()) return;
 
     // Old, hackier solution
     /*
     std::cout << gPad->GetListOfPrimitives()->GetEntries() << std::endl;
     // Which pads have we already drawn axes in? Never draw axes in them
     // again. Unfortunately UniqueID() never seems to be set. If that's the
     // case, set it to a random value and hope...
     static std::set<UInt_t> already;
     if(already.count(gPad->GetUniqueID())) return;
     if(gPad->GetUniqueID() == 0) gPad->SetUniqueID(rand());
     already.insert(gPad->GetUniqueID());
     */
 
     const TAxis *ax = fHist->GetXaxis();
     const TAxis *ay = fHist->GetYaxis();
     const double Nx = ax->GetNbins();
     const double Ny = ay->GetNbins();
 
     // Axes with limits where the user originally requested them, which we
     // adjusted to be the centres of the first and last bins.
     TH2 *axes = new TH2C(UniqueName().c_str(),
                          TString::Format(";%s;%s",
                                          ax->GetTitle(), ay->GetTitle()),
                          Nx - 1, ax->GetBinCenter(1), ax->GetBinCenter(Nx),
                          Ny - 1, ay->GetBinCenter(1), ay->GetBinCenter(Ny));
     axes->Draw();
 
     if(fHist){
       // "colz same" will reuse axis's min and max, so set them helpfully here
       axes->SetMinimum(fHist->GetMinimum());
       axes->SetMaximum(fHist->GetMaximum());
     }
 
     axes->SetTitle(fHist->GetTitle());
     axes->GetXaxis()->SetLabelSize(ax->GetLabelSize());
     axes->GetYaxis()->SetLabelSize(ay->GetLabelSize());
     axes->GetXaxis()->SetLabelOffset(ax->GetLabelOffset());
     axes->GetYaxis()->SetLabelOffset(ay->GetLabelOffset());
     axes->GetXaxis()->CenterTitle();
     axes->GetYaxis()->CenterTitle();
     gPad->Update();
   }

void ana::FrequentistSurface::FillSurface	(	const IExperiment *	expt,
		osc::IOscCalculatorAdjustable *	calc,
		const IFitVar *	xvar,
		const IFitVar *	yvar,
		const std::vector< const IFitVar * > &	profVars,
		const std::vector< const ISyst * > &	profSysts,
		const SeedList &	seedPts,
		const std::vector< SystShifts > &	systSeedPts
	)

protectedvirtual

Reimplemented in ana::SurfaceKrige.

Definition at line 107 of file FrequentistSurface.cxx.

References ana::ThreadPool::AddMemberTask(), bin, om::cerr, ana::Progress::Done(), e, allTimeWatchdog::endl, ana::ISurface::fHist, FillSurfacePoint(), ana::ThreadPool::Finish(), fParallel, Get, ana::IFitVar::Penalty(), time_estimates::pool, cacheDefinitionData::prog, ProgressBarTitle(), ana::Progress::SetProgress(), ana::IFitVar::ShortName(), ana::ThreadPool::ShowProgress(), fillBadChanDBTables::step, submit_syst::x, and submit_syst::y.

Referenced by FrequentistSurface().

   {
     if (fParallel && !(profVars.empty() && profSysts.empty()))
     {
       // Minuit calls this at some point, and it creates a static. Which
       // doesn't like happening on multiple threads at once. So do it upfront
       // while we're still single-threaded.
       ROOT::Minuit2::StackAllocatorHolder::Get();
     }
 
     // Nothing created during surface filling belongs in a
     // directory. Unfortunately the local guards in Spectrum etc are racey when
     // run in parallel. But this should cover the whole lot safely.
     DontAddDirectory guard;
 
     const std::string progTitle = ProgressBarTitle(xvar, yvar, profVars, profSysts);
 
     Progress *prog = 0;
     // Difficult to keep a progress bar properly up to date when threaded
     if (!fParallel) prog = new Progress(progTitle);
     ThreadPool *pool = 0;
     if (fParallel)
     {
       pool = new ThreadPool;
       pool->ShowProgress(progTitle);
     }
 
     const int Nx = fHist->GetNbinsX();
     const int Ny = fHist->GetNbinsY();
 
     // Fill bins in "random" order so that the progress bar is accurate
     // straight away instead of possibly being misled by whatever atypical
     // points we start with. This step is a prime which guarantees we get every
     // cell.
     int step = 7919;
     // Very unlikely (Nx or Ny is a multiple of step), but just to be safe.
     if ((Nx * Ny) % step == 0) step = 1;
 
     int bin = 0;
     int neval = 0;
 
     do
     {
       const int x = bin % Nx + 1;
       const int y = bin / Nx + 1;
 
       const double xv = fHist->GetXaxis()->GetBinCenter(x);
       const double yv = fHist->GetYaxis()->GetBinCenter(y);
 
       if (xvar->Penalty(xv, calc) > 1e-10)
       {
         std::cerr << "Warning! " << xvar->ShortName() << " = " << xv
                   << " has penalty of " << xvar->Penalty(xv, calc)
                   << " that could have been applied in surface. "
                   << "This should never happen." << std::endl;
       }
       if (yvar->Penalty(yv, calc) > 1e-10)
       {
         std::cerr << "Warning! " << yvar->ShortName() << " = " << yv
                   << " has penalty of " << yvar->Penalty(yv, calc)
                   << " that could have been applied in surface. "
                   << "This should never happen." << std::endl;
       }
 
       if (fParallel)
       {
         pool->AddMemberTask(this, &FrequentistSurface::FillSurfacePoint,
                             expt, calc,
                             xvar, xv, yvar, yv,
                             profVars, profSysts, seedPts, systSeedPts);
       } else
       {
         FillSurfacePoint(expt, calc,
                           xvar, xv, yvar, yv,
                           profVars, profSysts, seedPts, systSeedPts);
         ++neval;
         prog->SetProgress(neval / double(Nx * Ny));
       }
 
       bin = (bin + step) % (Nx * Ny);
     } while (bin != 0);
 
 
     if (fParallel)
     {
       pool->Finish();
       delete pool;
     } else
     {
       prog->Done();
       delete prog;
     }
   }

double ana::FrequentistSurface::FillSurfacePoint	(	const IExperiment *	expt,
		osc::IOscCalculatorAdjustable *	calc,
		const IFitVar *	xvar,
		double	x,
		const IFitVar *	yvar,
		double	y,
		const std::vector< const IFitVar * > &	profVars,
		const std::vector< const ISyst * > &	profSysts,
		const SeedList &	seedPts,
		const std::vector< SystShifts > &	systSeedPts
	)

protected

Definition at line 245 of file FrequentistSurface.cxx.

References calc, ana::IExperiment::ChiSq(), osc::_IOscCalculatorAdjustable< T >::Copy(), fFitOpts, ana::ISurface::fHist, ana::IFitter::Fit(), fParallel, fProfHists, ana::ISurface::fSeedValues, ana::SystShifts::GetShift(), MECModelEnuComparisons::i, makeTrainCVSamples::int, calib::j, ana::IFitter::kQuiet, ana::MinuitFitter::SetFitOpts(), and ana::IFitVar::SetValue().

Referenced by FillSurface(), and FrequentistSurface().

   {
     osc::IOscCalculatorAdjustable* calcNoHash = 0; // specific to parallel mode
 
     if(fParallel){
       // Need to take our own copy so that we don't get overwritten by someone
       // else's changes.
       calc = calc->Copy();
     }
     
     xvar->SetValue(calc, x);
     yvar->SetValue(calc, y);
 
     if (fParallel){
       calcNoHash = new OscCalcNoHash(calc);
     }
 
     //Make sure that the profiled values of fitvars do not persist between steps.
     for(int i = 0; i < (int)fSeedValues.size(); ++i) profVars[i]->SetValue( calc, fSeedValues[i] );
 
     double chi;
     if(profVars.empty() && profSysts.empty()){
       chi = expt->ChiSq(fParallel ? calcNoHash : calc);
     }
     else{
       MinuitFitter fitter(expt, profVars, profSysts);
       fitter.SetFitOpts(fFitOpts);
       SystShifts bestSysts;
       chi = fitter.Fit(calc, bestSysts, seedPts, systSeedPts, MinuitFitter::kQuiet);
 
       for(unsigned int i = 0; i < profVars.size(); ++i){
         fProfHists[i]->Fill(x, y, profVars[i]->GetValue(calc));
       }
       for(unsigned int j = 0; j < profSysts.size(); ++j){
         fProfHists[j+profVars.size()]->Fill(x, y, bestSysts.GetShift(profSysts[j]));
       }
     }
 
     fHist->Fill(x, y, chi);
 
     if(fParallel){
       delete calc;
       delete calcNoHash;
     }
     return chi;
   }

void ana::ISurface::FillSurfObj	(	ISurface &	surf,
		TDirectory *	dir
	)

staticprotectedinherited

Definition at line 219 of file ISurface.cxx.

References ana::ISurface::fBestFitX, ana::ISurface::fBestFitY, ana::ISurface::fBestLikelihood, ana::ISurface::fHist, ana::ISurface::fSeedValues, compare_h5_caf::idx, and registry_explorer::v.

Referenced by ana::ISurface::GetBestFitY(), LoadFrom(), and ana::BayesianSurface::LoadFrom().

   {
 
     const TVectorD v = *(TVectorD *) dir->Get("minValues");
     const TVectorD s = *(TVectorD *) dir->Get("seeds");
 
     surf.fHist = (TH2F *) dir->Get("hist");
     surf.fBestLikelihood = v[0];
     surf.fBestFitX = v[1];
     surf.fBestFitY = v[2];
 
     for (int idx = 0; idx < s.GetNrows(); ++idx)
       surf.fSeedValues.push_back(s[idx]);
 
   }

void ana::FrequentistSurface::FindMinimum	(	const IExperiment *	expt,
		osc::IOscCalculatorAdjustable *	calc,
		const IFitVar *	xvar,
		const IFitVar *	yvar,
		const std::vector< const IFitVar * > &	profVars,
		const std::vector< const ISyst * > &	profSysts,
		const SeedList &	seedPts,
		const std::vector< SystShifts > &	systSeedPts
	)

protected

Definition at line 301 of file FrequentistSurface.cxx.

References ana::ISurface::fBestFitX, ana::ISurface::fBestFitY, ana::ISurface::fBestLikelihood, fFitOpts, ana::ISurface::fHist, ana::IFitter::Fit(), ana::ISurface::fSeedValues, ana::IFitVar::GetValue(), MECModelEnuComparisons::i, makeTrainCVSamples::int, ana::SystShifts::Nominal(), ana::MinuitFitter::SetFitOpts(), ana::IFitVar::SetValue(), submit_syst::x, and submit_syst::y.

Referenced by FrequentistSurface().

   {
     // Location of the best minimum found from filled surface
     double minchi = 1e10;
     int minx = fHist->GetNbinsX()/2;
     int miny = fHist->GetNbinsY()/2;
     for(int x = 1; x <= fHist->GetNbinsX(); ++x){
       for(int y = 1; y <= fHist->GetNbinsY(); ++y){
         const double chi = fHist->GetBinContent(x, y);
         if(chi < minchi){
           minchi = chi;
           minx = x;
           miny = y;
         }
       }
     }
 
     std::vector<const IFitVar*> allVars = {xvar, yvar};
     allVars.insert(allVars.end(), profVars.begin(), profVars.end());
     MinuitFitter fit(expt, allVars, profSysts);
     fit.SetFitOpts(fFitOpts);
     // Seed from best grid point
     xvar->SetValue(calc, fHist->GetXaxis()->GetBinCenter(minx));
     yvar->SetValue(calc, fHist->GetYaxis()->GetBinCenter(miny));
     for(int i = 0; i < (int)fSeedValues.size(); ++i) profVars[i]->SetValue( calc, fSeedValues[i] );
     SystShifts systSeed = SystShifts::Nominal();
     fBestLikelihood = fit.Fit(calc, systSeed, seedPts);
     fBestFitX = xvar->GetValue(calc);
     fBestFitY = yvar->GetValue(calc);
 
     for(int x = 0; x < fHist->GetNbinsX()+2; ++x){
       for(int y = 0; y < fHist->GetNbinsY()+2; ++y){
         fHist->SetBinContent(x, y, fHist->GetBinContent(x, y)-fBestLikelihood);
       }
     }
 
     fHist->SetMinimum(0);
   }

double ana::ISurface::GetBestFitX ( ) const

inlineinherited

Definition at line 30 of file ISurface.h.

References ana::ISurface::fBestFitX.

Referenced by test_stanfit_statsonly().

30 {return fBestFitX;}

ana::ISurface::fBestFitX

double fBestFitX

Definition: ISurface.h:60

double ana::ISurface::GetBestFitY ( ) const

inlineinherited

Definition at line 31 of file ISurface.h.

References plot_validation_datamc_2018::color, dir, ana::ISurface::Draw(), ana::ISurface::DrawBestFit(), ana::ISurface::DrawContour(), ana::ISurface::fBestFitY, ana::ISurface::FillSurfObj(), ana::ISurface::GetGraphs(), ana::ISurface::SaveTo(), ana::ISurface::SetTitle(), runNovaSAM::str, demo5::surf, and ana::ISurface::ToTH2().

Referenced by test_stanfit_statsonly().

31 {return fBestFitY;}

ana::ISurface::fBestFitY

double fBestFitY

Definition: ISurface.h:61

std::vector< TGraph * > ana::ISurface::GetGraphs	(	TH2 *	fc,
		double	minchi = `-1`
	)

inherited

For expert use, custom painting of contours.

Definition at line 115 of file ISurface.cxx.

References ana::ISurface::fBestLikelihood, ana::ISurface::fHist, runNovaSAM::ret, PandAna.Demos.tute_pid_validation::specs, demo5::surf, tmp, and ana::UniqueName().

Referenced by ana::ISurface::DrawContour(), and ana::ISurface::GetBestFitY().

   {
     std::vector<TGraph*> ret;
 
     if(minchi < 0) minchi = fBestLikelihood;
     std::unique_ptr<TH2> surf((TH2*)fHist->Clone(UniqueName().c_str()));
     surf->Add(fc, -1);
 
     TVirtualPad* bak = gPad;
 
     const bool wasbatch = gROOT->IsBatch();
     gROOT->SetBatch(); // User doesn't want to see our temporary canvas
     TCanvas tmp;
 
     gStyle->SetOptStat(0);
 
     const double level = minchi-fBestLikelihood;
     surf->SetContour(1, &level);
     surf->Draw("cont list");
 
     tmp.Update();
     tmp.Paint();
 
     gROOT->SetBatch(wasbatch);
     gPad = bak;
 
     // The graphs we need (contained inside TLists, contained inside
     // TObjArrays) are in the list of specials. But we need to be careful about
     // types, because other stuff can get in here too (TDatabasePDG for
     // example).
     TCollection* specs = gROOT->GetListOfSpecials();
 
     TIter nextSpec(specs);
     while(TObject* spec = nextSpec()){
       if(!spec->InheritsFrom(TObjArray::Class())) continue;
       TObjArray* conts = (TObjArray*)spec;
 
       if(conts->IsEmpty()) continue;
 
       if(!conts->At(0)->InheritsFrom(TList::Class())) continue;
       TList* cont = (TList*)conts->At(0);
 
       TIter nextObj(cont);
       // Contour could be split into multiple pieces
       while(TObject* obj = nextObj()){
         if(!obj->InheritsFrom(TGraph::Class())) continue;
 
         ret.push_back((TGraph*)obj->Clone(UniqueName().c_str()));
       } // end for obj
     } // end for spec
 
     return ret;
   }

std::vector<TH2*> ana::FrequentistSurface::GetProfiledHists ( )

inline

Maps of the values taken on by the profiled parameters.

Definition at line 60 of file FrequentistSurface.h.

References dir, fProfHists, LoadFrom(), and SaveTo().

Referenced by run_joint_fit_2020_contours().

60 {return fProfHists;}

ana::FrequentistSurface::fProfHists

std::vector< TH2 * > fProfHists

Definition: FrequentistSurface.h:107

std::unique_ptr< FrequentistSurface > ana::FrequentistSurface::LoadFrom ( TDirectory * dir )

static

Definition at line 367 of file FrequentistSurface.cxx.

References ana::assert(), ana::ISurface::FillSurfObj(), genie::utils::style::Format(), make_syst_table_plots::h, compare_h5_caf::idx, demo5::surf, and getGoodRuns4SAM::tag.

Referenced by compare_fits(), demoFitContours(), DrawExtrapSurface(), GetProfiledHists(), joint_fit_2018_contours(), joint_fit_2019_contours(), median_contours(), plot_joint_fit_2020_contours(), Plotting_CompareMultipleContours(), sensitivity2018(), and sensitivity2020().

   {
     DontAddDirectory guard;
 
     TObjString *tag = (TObjString *) dir->Get("type");
     assert(tag);
     assert(tag->GetString() == "FrequentistSurface" || tag->GetString() == "Surface");
 
     std::unique_ptr<FrequentistSurface> surf(new FrequentistSurface);
     ISurface::FillSurfObj(*surf, dir);
 
     for (std::size_t idx = 0; idx < surf->fSeedValues.size(); ++idx)
     {
       // Search for old "marg" name here too for backwards compatibility
       TH2 *h = (TH2 *) dir->Get(TString::Format("profHists/hist%lu", idx));
       if (h)
         surf->fProfHists.push_back(h);
       else
         surf->fProfHists.push_back((TH2 *) dir->Get(TString::Format("margHists/hist%lu", idx)));
     }
 
     return surf;
   }

std::string ana::FrequentistSurface::ProgressBarTitle	(	const IFitVar *	xvar,
		const IFitVar *	yvar,
		const std::vector< const IFitVar * > &	profVars,
		const std::vector< const ISyst * > &	profSysts
	)		const

protected

Definition at line 86 of file FrequentistSurface.cxx.

References ana::IFitVar::ShortName(), plotROC::title, and registry_explorer::v.

Referenced by CreateHistograms(), FillSurface(), and FrequentistSurface().

   {
     std::string title = "Filling surface for "+yvar->ShortName()+" vs "+xvar->ShortName();
 
     if(!profVars.empty() || !profSysts.empty()){
       title += " (profiling ";
 
       for(const IFitVar* v: profVars) title += v->ShortName() + ", ";
       for(const ISyst* s: profSysts)  title += s->ShortName() + ", ";
 
       // Always have one superfluous ", " at the end
       title.resize(title.size()-2);
       title += ")";
     }
 
     return title;
   }

void ana::FrequentistSurface::SaveTo ( TDirectory * dir ) const

Definition at line 347 of file FrequentistSurface.cxx.

References genie::utils::style::Format(), fProfHists, compare_h5_caf::idx, it, ana::ISurface::SaveTo(), and tmp.

Referenced by getContProf(), getContProf_Sensitivity(), GetProfiledHists(), joint_fit_future_contour_univ(), make_surfprof(), make_surfprof_sensitivity(), MakeExtrapSurface(), run_joint_fit_2020_contours(), sensitivity2018(), and sensitivity2020().

   {
     TDirectory *tmp = gDirectory;
     dir->cd();
     TObjString("FrequentistSurface").Write("type");
 
     ISurface::SaveTo(dir);
 
     TDirectory *profDir = dir->mkdir("profHists");
     int idx = 0;
     for (auto it: FrequentistSurface::fProfHists)
     {
       profDir->cd();
       it->Write(TString::Format("hist%d", idx++));
     }
 
     tmp->cd();
   }

void ana::ISurface::SetTitle ( const char * str )

inherited

Definition at line 192 of file ISurface.cxx.

References ana::ISurface::fHist.

Referenced by ana::ISurface::GetBestFitY(), plots(), ana::NumuSurface::SetTitle(), ana::NueSurface::SetTitle(), and test_external_constraints().

   {
     fHist->SetTitle(str);
   }

TH2 * ana::ISurface::ToTH2 ( double minchi = -1 ) const

inherited

Definition at line 170 of file ISurface.cxx.

References ana::ISurface::fBestLikelihood, ana::ISurface::fHist, runNovaSAM::ret, submit_syst::x, and submit_syst::y.

Referenced by CAF_makeCAFSensitivities_for_FNEX(), demoFitContours(), ana::Flat(), ana::ISurface::GetBestFitY(), joint_fit_2017_contours(), joint_fit_2018_contours(), joint_fit_2019_contours(), joint_fit_future_contour_univ(), MakeCAFSensitivities_for_FNEX(), Plot2DSlice(), Plot2DSlices(), run_joint_fit_2020_contours(), sensitivity2018(), sensitivity2020(), and test_stanfit_withsysts().

   {
     // Could have a file temporarily open
     DontAddDirectory guard;
 
     TH2 *ret = new TH2F(*fHist);
 
     if (minchi >= 0)
     {
       for (int x = 0; x < ret->GetNbinsX() + 2; ++x)
       {
         for (int y = 0; y < ret->GetNbinsY() + 2; ++y)
         {
           ret->SetBinContent(x, y, ret->GetBinContent(x, y) + fBestLikelihood - minchi);
         }
       }
     }
 
     return ret;
   }

Friends And Related Function Documentation

friend class NueSurface

friend

Definition at line 30 of file FrequentistSurface.h.

friend class NumuSurface

friend

Definition at line 29 of file FrequentistSurface.h.

Member Data Documentation

double ana::ISurface::fBestFitX

protectedinherited

Definition at line 60 of file ISurface.h.

Referenced by ana::BayesianSurface::BuildHist(), ana::ISurface::DrawBestFit(), ana::ISurface::FillSurfObj(), FindMinimum(), ana::ISurface::GetBestFitX(), and ana::ISurface::SaveTo().

double ana::ISurface::fBestFitY

protectedinherited

Definition at line 61 of file ISurface.h.

Referenced by ana::BayesianSurface::BuildHist(), ana::ISurface::DrawBestFit(), ana::ISurface::FillSurfObj(), FindMinimum(), ana::ISurface::GetBestFitY(), and ana::ISurface::SaveTo().

double ana::ISurface::fBestLikelihood

protectedinherited

Definition at line 59 of file ISurface.h.

Referenced by ana::ISurface::BestLikelihood(), ana::BayesianSurface::BuildHist(), ana::ISurface::FillSurfObj(), FindMinimum(), ana::ISurface::GetGraphs(), ana::ISurface::SaveTo(), and ana::ISurface::ToTH2().

MinuitFitter::FitOpts ana::FrequentistSurface::fFitOpts

protected

Definition at line 104 of file FrequentistSurface.h.

Referenced by FillSurfacePoint(), and FindMinimum().

TH2F* ana::ISurface::fHist

protectedinherited

Definition at line 62 of file ISurface.h.

Referenced by ana::BayesianSurface::BuildHist(), CreateHistograms(), ana::ISurface::Draw(), ana::ISurface::EnsureAxes(), ana::NumuSurface::EnsurePad(), ana::NueSurface::EnsurePad(), FillSurface(), FillSurfacePoint(), ana::ISurface::FillSurfObj(), FindMinimum(), ana::ISurface::GetGraphs(), ana::ISurface::SaveTo(), ana::ISurface::SetTitle(), and ana::ISurface::ToTH2().

bool ana::FrequentistSurface::fParallel

protected

Definition at line 102 of file FrequentistSurface.h.

Referenced by FillSurface(), and FillSurfacePoint().

std::vector<TH2*> ana::FrequentistSurface::fProfHists

protected

Definition at line 107 of file FrequentistSurface.h.

Referenced by CreateHistograms(), FillSurfacePoint(), GetProfiledHists(), and SaveTo().

std::vector<double> ana::ISurface::fSeedValues

protectedinherited

Definition at line 63 of file ISurface.h.

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

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

/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-08-04/CAFAna/Fit/FrequentistSurface.h
/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-08-04/CAFAna/Fit/FrequentistSurface.cxx

Public Member Functions

Static Public Member Functions

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Friends And Related Function Documentation

Member Data Documentation