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

Inheritance diagram for ana::ISurface:

Public Member Functions
	ISurface ()

virtual	~ISurface ()

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)

Protected Member Functions
double	BinCenterX (int bin) const

double	BinCenterY (int bin) const

void	SaveToHelper (TDirectory *dir) const
	dir should already be the appropriate sub-directory More...

void	EnsureAxes () const

void	CheckMask (const std::string &func) const

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

Protected Attributes
double	fBestLikelihood

double	fBestFitX

double	fBestFitY

TH2F *	fHist

bool	fLogX

bool	fLogY

std::vector< double >	fSeedValues

std::vector< int >	fBinMask

Friends
TH2 *	Flat (double, const ISurface &)
	Helper function for the gaussian approximation surfaces. More...

Detailed Description

Definition at line 17 of file ISurface.h.

Constructor & Destructor Documentation

ana::ISurface::ISurface ( )

inline

Definition at line 20 of file ISurface.h.

         : fBestLikelihood(std::numeric_limits<double>::signaling_NaN()),
           fBestFitX(std::numeric_limits<double>::signaling_NaN()),
           fBestFitY(std::numeric_limits<double>::signaling_NaN()),
           fHist(nullptr), fLogX(false), fLogY(false)
       {}

virtual ana::ISurface::~ISurface ( )

inlinevirtual

Definition at line 27 of file ISurface.h.

27 {}

Member Function Documentation

double ana::ISurface::BestLikelihood ( ) const

inline

Definition at line 29 of file ISurface.h.

References fBestLikelihood.

29 {return fBestLikelihood;}

ana::ISurface::fBestLikelihood

double fBestLikelihood

Definition: ISurface.h:65

double ana::ISurface::BinCenterX ( int bin ) const

protected

Definition at line 265 of file ISurface.cxx.

References visualisationForPaperMasterPlot::ax, fHist, and fLogX.

Referenced by ana::BayesianSurface::BuildHist(), EnsureAxes(), ana::FrequentistSurface::FillSurface(), ana::FrequentistSurface::FindMinimum(), and GetBestFitY().

   {
     const TAxis* ax = fHist->GetXaxis();
     return fLogX ? ax->GetBinCenterLog(bin) : ax->GetBinCenter(bin);
   }

double ana::ISurface::BinCenterY ( int bin ) const

protected

Definition at line 272 of file ISurface.cxx.

References visualisationForPaperMasterPlot::ax, fHist, and fLogY.

Referenced by ana::BayesianSurface::BuildHist(), EnsureAxes(), ana::FrequentistSurface::FillSurface(), ana::FrequentistSurface::FindMinimum(), and GetBestFitY().

   {
     const TAxis* ax = fHist->GetYaxis();
     return fLogY ? ax->GetBinCenterLog(bin) : ax->GetBinCenter(bin);
   }

void ana::ISurface::CheckMask ( const std::string & func ) const

protected

Definition at line 279 of file ISurface.cxx.

References om::cout, allTimeWatchdog::endl, and fBinMask.

Referenced by DrawBestFit(), DrawContour(), and GetGraphs().

   {
     if(!fBinMask.empty()) {
       std::cout << "Cannot call " << func << "() on a partial surface!" << std::endl;
       abort();
     }
   }

void ana::ISurface::Draw ( ) const

Draw the surface itself.

Definition at line 21 of file ISurface.cxx.

References EnsureAxes(), and fHist.

Referenced by cc(), demo5(), GetBestFitY(), modularextrap_demo_nue(), modularextrap_demo_numu(), Plot2DSlices(), plot_3flavor_withsysts(), test_ana(), test_stanfit_statsonly(), test_stanfit_withsysts(), and test_surf_stride().

   {
     // Can be useful to draw a partial surface for debugging
     //    CheckMask();
 
     EnsureAxes();
 
     fHist->Draw("colz same");
 
     // colz obliterated them
     gPad->RedrawAxis();
 
     gPad->Update();
   }

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

Draw the best fit point.

Definition at line 37 of file ISurface.cxx.

References CheckMask(), EnsureAxes(), fBestFitX, and fBestFitY.

Referenced by cc(), demo5(), demo_CPT(), drawSensitivity(), GetBestFitY(), getSensitivity(), make_contours(), modularextrap_demo_nue(), modularextrap_demo_numu(), Plot2DSlices(), plot_3flavor_withsysts(), starPlot(), syst_test(), template_basic(), template_GENIE_systs(), template_nonGENIE_systs(), test_ana(), test_stanfit_statsonly(), and test_stanfit_withsysts().

   {
     CheckMask("DrawBestFit");
     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`
	)

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 50 of file ISurface.cxx.

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

Referenced by cc(), demo5(), Draw2DSurface(), DrawContours(), drawSensitivity(), GetBestFitY(), getSensitivity(), make_contours(), make_nom_expt(), modularextrap_demo_nue(), modularextrap_demo_numu(), Plot2DSlice(), Plot2DSlices(), plot_3flavor_withsysts(), plotContProf(), plots(), starPlot(), syst_test(), template_basic(), template_GENIE_systs(), template_nonGENIE_systs(), test_ana(), test_stanfit_statsonly(), and test_stanfit_withsysts().

   {
     CheckMask("DrawContour");
     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

protected

Definition at line 70 of file ISurface.cxx.

References visualisationForPaperMasterPlot::ax, file_size_ana::axes, BinCenterX(), BinCenterY(), fHist, fLogX, fLogY, genie::utils::style::Format(), and ana::UniqueName().

Referenced by Draw(), DrawBestFit(), and DrawContour().

   {
     // 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;
 
     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, BinCenterX(1), BinCenterX(Nx),
                          Ny - 1, BinCenterY(1), BinCenterY(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();
 
     if(fLogX) gPad->SetLogx();
     if(fLogY) gPad->SetLogy();
 
     gPad->Update();
   }

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

staticprotected

Definition at line 229 of file ISurface.cxx.

References fBestFitX, fBestFitY, fBestLikelihood, fBinMask, fHist, fSeedValues, compare_h5_caf::idx, m, and registry_explorer::v.

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

   {
 
     const TVectorD v = *(TVectorD *) dir->Get("minValues");
     const TVectorD s = *(TVectorD *) dir->Get("seeds");
     TVectorD* m = (TVectorD *) dir->Get("mask");
 
     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]);
 
     if (m) {
       for (int idx = 0; idx < m->GetNrows(); ++idx)
         surf.fBinMask.push_back((*m)[idx]);
     }
 
   }

double ana::ISurface::GetBestFitX ( ) const

inline

Definition at line 30 of file ISurface.h.

References fBestFitX.

Referenced by test_stanfit_statsonly().

30 {return fBestFitX;}

ana::ISurface::fBestFitX

double fBestFitX

Definition: ISurface.h:66

double ana::ISurface::GetBestFitY ( ) const

inline

Definition at line 31 of file ISurface.h.

References bin, BinCenterX(), BinCenterY(), plot_validation_datamc_2018::color, dir, Draw(), DrawBestFit(), DrawContour(), fBestFitY, FillSurfObj(), Flat, GetGraphs(), SaveToHelper(), SetTitle(), submit_syst::str, demo5::surf, and ToTH2().

Referenced by test_stanfit_statsonly().

31 {return fBestFitY;}

ana::ISurface::fBestFitY

double fBestFitY

Definition: ISurface.h:67

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

For expert use, custom painting of contours.

Definition at line 113 of file ISurface.cxx.

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

Referenced by DrawContour(), and GetBestFitY().

   {
     CheckMask("GetGraphs");
 
     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;
   }

void ana::ISurface::SaveToHelper ( TDirectory * dir ) const

protected

dir should already be the appropriate sub-directory

Definition at line 198 of file ISurface.cxx.

References fBestFitX, fBestFitY, fBestLikelihood, fBinMask, fHist, fLogX, fLogY, fSeedValues, m, tmp, and registry_explorer::v.

Referenced by GetBestFitY(), ana::BayesianSurface::SaveTo(), and ana::FrequentistSurface::SaveTo().

   {
     TDirectory *oldDir = gDirectory;
 
     dir->cd();
 
     TVectorD v(3);
     v[0] = fBestLikelihood;
     v[1] = fBestFitX;
     v[2] = fBestFitY;
     v.Write("minValues");
 
     fHist->Write("hist");
 
     TVectorD s(fSeedValues.size(), &fSeedValues[0]);
     s.Write("seeds");
 
     if (!fBinMask.empty()) {
       std::vector<double> tmp(fBinMask.begin(), fBinMask.end());
       TVectorD m(tmp.size(), &tmp[0]);
       m.Write("mask");
     }
 
     TObjString(fLogX ? "yes" : "no").Write("logx");
     TObjString(fLogY ? "yes" : "no").Write("logy");
 
     if (oldDir)
       oldDir->cd();
   }

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

Definition at line 192 of file ISurface.cxx.

References fHist.

Referenced by GetBestFitY(), and plots().

   {
     fHist->SetTitle(str);
   }

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

Definition at line 170 of file ISurface.cxx.

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

Referenced by CAF_makeCAFSensitivities_for_FNEX(), demoFitContours(), GetBestFitY(), joint_fit_2017_contours(), joint_fit_2018_contours(), joint_fit_2019_contours(), joint_fit_future_contour_univ(), ana::FrequentistSurface::LoadFromMulti(), MakeCAFSensitivities_for_FNEX(), Plot2DSlice(), Plot2DSlices(), plot_3flavor_withsysts(), 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

TH2* Flat	(	double	,
		const ISurface &
	)

friend

Helper function for the gaussian approximation surfaces.

Definition at line 253 of file ISurface.cxx.

Referenced by ana::Gaussian2Sigma1D(), ana::Gaussian2Sigma2D(), ana::Gaussian3Sigma1D(), ana::Gaussian3Sigma2D(), ana::Gaussian68Percent1D(), ana::Gaussian68Percent2D(), ana::Gaussian90Percent1D(), ana::Gaussian90Percent2D(), ana::Gaussian95Percent1D(), ana::Gaussian95Percent2D(), ana::Gaussian99Percent1D(), ana::Gaussian99Percent2D(), GetBestFitY(), and ana::BayesianSurface::QuantileSurface().

   {
     TH2* h = new TH2F(*s.fHist);
 
     for(int x = 0; x < h->GetNbinsX()+2; ++x)
       for(int y = 0; y < h->GetNbinsY()+2; ++y)
         h->SetBinContent(x, y, level);
 
     return h;
   }

Member Data Documentation

double ana::ISurface::fBestFitX

protected

Definition at line 66 of file ISurface.h.

Referenced by ana::BayesianSurface::BuildHist(), DrawBestFit(), FillSurfObj(), ana::FrequentistSurface::FindMinimum(), GetBestFitX(), and SaveToHelper().

double ana::ISurface::fBestFitY

protected

Definition at line 67 of file ISurface.h.

Referenced by ana::BayesianSurface::BuildHist(), DrawBestFit(), FillSurfObj(), ana::FrequentistSurface::FindMinimum(), GetBestFitY(), and SaveToHelper().

double ana::ISurface::fBestLikelihood

protected

Definition at line 65 of file ISurface.h.

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

std::vector<int> ana::ISurface::fBinMask

protected

Definition at line 71 of file ISurface.h.

Referenced by CheckMask(), ana::FrequentistSurface::FillSurface(), FillSurfObj(), ana::FrequentistSurface::FindMinimum(), ana::FrequentistSurface::LoadFromMulti(), and SaveToHelper().

TH2F* ana::ISurface::fHist

protected

Definition at line 68 of file ISurface.h.

Referenced by BinCenterX(), BinCenterY(), ana::BayesianSurface::BuildHist(), ana::FrequentistSurface::CreateHistograms(), Draw(), EnsureAxes(), ana::FrequentistSurface::FillSurface(), ana::FrequentistSurface::FillSurfacePoint(), FillSurfObj(), ana::FrequentistSurface::FindMinimum(), ana::Flat(), GetGraphs(), ana::BayesianSurface::QuantileSurface(), SaveToHelper(), SetTitle(), and ToTH2().

bool ana::ISurface::fLogX

protected

Definition at line 69 of file ISurface.h.

Referenced by ana::BayesianSurface::BayesianSurface(), BinCenterX(), ana::BayesianSurface::BuildHist(), EnsureAxes(), ana::FrequentistSurface::FrequentistSurface(), and SaveToHelper().

bool ana::ISurface::fLogY

protected

Definition at line 69 of file ISurface.h.

Referenced by ana::BayesianSurface::BayesianSurface(), BinCenterY(), ana::BayesianSurface::BuildHist(), EnsureAxes(), ana::FrequentistSurface::FrequentistSurface(), and SaveToHelper().

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

protected

Definition at line 70 of file ISurface.h.

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

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

/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-04-15/CAFAna/Fit/ISurface.h
/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-04-15/CAFAna/Fit/ISurface.cxx

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