ana::ISurface Class Reference

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-10-20/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

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

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.

{}

Member Function Documentation

double ana::ISurface::BestLikelihood ( ) const

inline

Definition at line 29 of file ISurface.h.

References fBestLikelihood.

{return fBestLikelihood;}

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

protected

Definition at line 244 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 251 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::Draw

(

)

const

Draw the surface itself.

Definition at line 19 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(), 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

Draw the best fit point.

Definition at line 32 of file ISurface.cxx.

References 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().

  {
    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	of best fit to compare against. Default: best fit from this surface.

Definition at line 44 of file ISurface.cxx.

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

Referenced by cc(), demo5(), Draw2DSurface(), 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().

  {
    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 63 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 214 of file ISurface.cxx.

References fBestFitX, fBestFitY, fBestLikelihood, fHist, fSeedValues, compare_h5_caf::idx, 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");

    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]);

  }

double ana::ISurface::GetBestFitX ( ) const

inline

Definition at line 30 of file ISurface.h.

References fBestFitX.

Referenced by test_stanfit_statsonly().

{return fBestFitX;}

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(), GetGraphs(), SaveToHelper(), SetTitle(), submit_syst::str, demo5::surf, and ToTH2().

Referenced by test_stanfit_statsonly().

{return fBestFitY;}

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

For expert use, custom painting of contours.

Definition at line 106 of file ISurface.cxx.

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

Referenced by DrawContour(), and 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;
  }

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

protected

dir should already be the appropriate sub-directory

Definition at line 189 of file ISurface.cxx.

References fBestFitX, fBestFitY, fBestLikelihood, fHist, fLogX, fLogY, fSeedValues, 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");

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

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

Referenced by CAF_makeCAFSensitivities_for_FNEX(), demoFitContours(), ana::Flat(), GetBestFitY(), joint_fit_2017_contours(), joint_fit_2018_contours(), joint_fit_2019_contours(), joint_fit_future_contour_univ(), 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;
  }

Member Data Documentation

double ana::ISurface::fBestFitX

protected

Definition at line 64 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 65 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 63 of file ISurface.h.

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

TH2F* ana::ISurface::fHist

protected

Definition at line 66 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(), GetGraphs(), ana::BayesianSurface::QuantileSurface(), SaveToHelper(), SetTitle(), and ToTH2().

bool ana::ISurface::fLogX

protected

Definition at line 67 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 67 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 68 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/N20-10-20/CAFAna/Fit/ISurface.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-10-20/CAFAna/Fit/ISurface.cxx