ana::ISurface Class Reference

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-08-04/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
void	SaveTo (TDirectory *dir) const
void	EnsureAxes () const

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

Protected Attributes
double	fBestLikelihood
double	fBestFitX
double	fBestFitY
TH2F *	fHist
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)
      {}

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.

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

{return fBestLikelihood;}

void ana::ISurface::Draw

(

)

const

Draw the surface itself.

Definition at line 20 of file ISurface.cxx.

References EnsureAxes(), and fHist.

Referenced by cc(), demo5(), ana::NueSurface::Draw(), ana::NumuSurface::Draw(), 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

Draw the best fit point.

Definition at line 33 of file ISurface.cxx.

References EnsureAxes(), fBestFitX, and fBestFitY.

Referenced by cc(), demo5(), demo_CPT(), ana::NumuSurface::DrawBestFit(), drawSensitivity(), 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
	)

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

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

Referenced by cc(), demo5(), Draw2DSurface(), ana::NueSurface::DrawContour(), ana::NumuSurface::DrawContour(), drawSensitivity(), 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

protected

Definition at line 64 of file ISurface.cxx.

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

Referenced by Draw(), DrawBestFit(), 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::ISurface::FillSurfObj ( ISurface &  surf, TDirectory *  dir  )

staticprotected

Definition at line 219 of file ISurface.cxx.

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

Referenced by GetBestFitY(), ana::FrequentistSurface::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]);

  }

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 plot_validation_datamc_2018::color, dir, Draw(), DrawBestFit(), DrawContour(), fBestFitY, FillSurfObj(), GetGraphs(), SaveTo(), SetTitle(), runNovaSAM::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 115 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::SaveTo ( TDirectory *  dir ) const

protected

Definition at line 198 of file ISurface.cxx.

References fBestFitX, fBestFitY, fBestLikelihood, fHist, fSeedValues, and registry_explorer::v.

Referenced by GetBestFitY(), ana::FrequentistSurface::SaveTo(), and ana::BayesianSurface::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 (oldDir)
      oldDir->cd();
  }

void ana::ISurface::SetTitle

(

const char *

str

)

Definition at line 192 of file ISurface.cxx.

References fHist.

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

  {
    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(), 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(), 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 60 of file ISurface.h.

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

double ana::ISurface::fBestFitY

protected

Definition at line 61 of file ISurface.h.

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

double ana::ISurface::fBestLikelihood

protected

Definition at line 59 of file ISurface.h.

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

TH2F* ana::ISurface::fHist

protected

Definition at line 62 of file ISurface.h.

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

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

protected

Definition at line 63 of file ISurface.h.

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

---

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

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