ana::MedianSurface Class Reference

#include <MedianSurface.h>

Inheritance diagram for ana::MedianSurface:

Public Member Functions
	MedianSurface (const std::vector< Surface > &throws)
void	DrawEnsemble (TH2 *fc, Color_t color=kGray)
void	DrawBand (TH2 *fc)
void	SaveTo (TDirectory *dir) const
Public Member Functions inherited from ana::Surface
	Surface (const IExperiment *expt, osc::IOscCalcAdjustable *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 std::map< const IFitVar *, std::vector< double >> &seedPts={}, const std::vector< SystShifts > &systSeedPts={}, bool parallel=false, Fitter::Precision prec=Fitter::kNormal)
	Surface (const IExperiment *expt, osc::IOscCalcAdjustable *calc, const FitAxis &xax, const FitAxis &yax, const std::vector< const IFitVar * > &profVars={}, const std::vector< const ISyst * > &profSysts={}, const std::map< const IFitVar *, std::vector< double >> &seedPts={}, const std::vector< SystShifts > &systSeedPts={}, bool parallel=false, Fitter::Precision prec=Fitter::kNormal)
void	Draw () const
	Draw the surface itself. More...
void	DrawBestFit (Color_t color, Int_t marker=kFullCircle) const
	Draw the best fit point. More...
double	MinChi () const
double	GetMinX () const
double	GetMinY () const
void	DrawContour (TH2 *fc, Style_t style, Color_t color, double minchi=-1)
TH2F *	ToTH2 (double minchi=-1) const
void	SetTitle (const char *str)
std::vector< TH2 * >	GetProfiledHists ()
	Maps of the values taken on by the profiled parameters. More...
std::vector< TH2 * >	GetMarginalizedHists ()
	Deprecated. Retained for backwards compatibility. More...
std::vector< TGraph * >	GetGraphs (TH2 *fc, double minchi=-1)
	For expert use, custom painting of contours. More...
void	SaveTo (TDirectory *dir) const

Static Public Member Functions
static std::unique_ptr < MedianSurface >	LoadFrom (TDirectory *dir)
Static Public Member Functions inherited from ana::Surface
static std::unique_ptr< Surface >	LoadFrom (TDirectory *dir)
static std::unique_ptr< Surface >	LoadFromMulti (const std::vector< TFile * > &files, const std::string &label)
static std::unique_ptr< Surface >	LoadFromMulti (const std::string &wildcard, const std::string &label)

Protected Attributes
std::vector< Surface >	fThrows
TH2F *	fHistUp1
TH2F *	fHistDn1
TH2F *	fHistUp2
TH2F *	fHistDn2
Protected Attributes inherited from ana::Surface
bool	fParallel
Fitter::Precision	fPrec
double	fMinChi
double	fMinX
double	fMinY
TH2F *	fHist
bool	fLogX
bool	fLogY
std::vector< TH2 * >	fProfHists
std::vector< double >	fSeedValues
std::vector< int >	fBinMask

Additional Inherited Members
Protected Member Functions inherited from ana::Surface
	Surface ()
void	EnsureAxes () const
void	CheckMask (const std::string &func) const
void	FillSurface (const std::string &progTitle, const IExperiment *expt, osc::IOscCalcAdjustable *calc, const FitAxis &xax, const FitAxis &yax, const std::vector< const IFitVar * > &profVars, const std::vector< const ISyst * > &profSysts, const std::map< const IFitVar *, std::vector< double >> &seedPts, const std::vector< SystShifts > &systSeedPts)
void	FillSurfacePoint (const IExperiment *expt, osc::IOscCalcAdjustable *calc, const FitAxis &xax, double x, const FitAxis &yax, double y, const std::vector< const IFitVar * > &profVars, const std::vector< const ISyst * > &profSysts, const std::map< const IFitVar *, std::vector< double >> &seedPts, const std::vector< SystShifts > &systSeedPts)

Detailed Description

Definition at line 9 of file MedianSurface.h.

Constructor & Destructor Documentation

ana::MedianSurface::MedianSurface

(

const std::vector< Surface > &

throws

)

Definition at line 15 of file MedianSurface.cxx.

    : fThrows(throws)
  {
    assert(!throws.empty());

    for(const Surface& s: throws){
      assert(s.fHist->GetNbinsX() == throws[0].fHist->GetNbinsX());
      assert(s.fHist->GetNbinsY() == throws[0].fHist->GetNbinsY());
      // TODO check min and max match at least
    }

    fLogX = throws[0].fLogX;
    fLogY = throws[0].fLogY;

    // Think about what to do with these
    fMinX = -1;
    fMinY = -1;
    // Though this is right
    fMinChi = 0;

    fHist = new TH2F(*throws[0].fHist);
    fHist->Reset();
    fHistUp1 = new TH2F(*fHist);
    fHistDn1 = new TH2F(*fHist);
    fHistUp2 = new TH2F(*fHist);
    fHistDn2 = new TH2F(*fHist);

    for(int ix = 0; ix < fHist->GetNbinsX()+2; ++ix){
      for(int iy = 0; iy < fHist->GetNbinsY()+2; ++iy){
        std::vector<float> chis;
        chis.reserve(throws.size());
        for(const Surface& s: throws){
          chis.push_back(s.fHist->GetBinContent(ix, iy));
        }
        std::sort(chis.begin(), chis.end());

        // TODO think about off-by-one errors and interpolation here
        // Median
        fHist->SetBinContent(ix, iy, chis[throws.size()/2]);
        // One sigma bounds
        const double tail1 = (1-0.6827)/2;
        fHistUp1->SetBinContent(ix, iy, chis[throws.size()*tail1]);
        fHistDn1->SetBinContent(ix, iy, chis[throws.size()*(1-tail1)]);

        // Two sigma bounds
        const double tail2 = (1-0.9545)/2;
        fHistUp2->SetBinContent(ix, iy, chis[throws.size()*tail2]);
        fHistDn2->SetBinContent(ix, iy, chis[(throws.size()-1)*(1-tail2)]);
      } // end for iy
    } // end for ix
  }

Member Function Documentation

void ana::MedianSurface::DrawBand

(

TH2 *

fc

)

Definition at line 86 of file MedianSurface.cxx.

  {
    EnsureAxes();

    TH2F surf1(*fHist);
    TH2F surf2(*fHist);

    for(int ix = 0; ix < fHist->GetNbinsX()+2; ++ix){
      for(int iy = 0; iy < fHist->GetNbinsY()+2; ++iy){
        const double c = fc->GetBinContent(ix, iy);

        const double u1 = fHistUp1->GetBinContent(ix, iy);
        const double d1 = fHistDn1->GetBinContent(ix, iy);
        surf1.SetBinContent(ix, iy, std::max(u1-c, c-d1));

        const double u2 = fHistUp2->GetBinContent(ix, iy);
        const double d2 = fHistDn2->GetBinContent(ix, iy);
        surf2.SetBinContent(ix, iy, std::max(u2-c, c-d2));
      }
    }

    const double level = 0;

    // I haven't been able to figure out how to draw these filled properly. cont0 does it, but can't handle countours that reach the edge of the space
    surf2.SetContour(1, &level);
    surf2.SetLineColor(kYellow);
    surf2.DrawCopy("cont3 same");

    surf1.SetContour(1, &level);
    surf1.SetLineColor(kGreen);
    surf1.SetFillColor(kGreen+2);
    surf1.DrawCopy("cont3 same");

    gPad->Update();
  }

void ana::MedianSurface::DrawEnsemble	(	TH2 *	fc,
		Color_t	color = kGray
	)

Definition at line 68 of file MedianSurface.cxx.

  {
    EnsureAxes();

    for(Surface& s: fThrows){
      std::vector<TGraph*> gs = s.GetGraphs(fc, -1);

      for(TGraph* g: gs){
        g->SetLineWidth(1);
        g->SetLineColor(color);
        g->Draw("l");
      }
    }

    gPad->Update();
  }

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

static

Definition at line 137 of file MedianSurface.cxx.

  {
    DontAddDirectory guard;

    TObjString* tag = (TObjString*)dir->Get("type");
    assert(tag);
    assert(tag->GetString() == "MedianSurface");

    std::vector<Surface> surfs;
    for(unsigned int i = 0; ; ++i){
      TDirectory* surfdir = dir->GetDirectory(TString::Format("surf%d", i));
      if(!surfdir) break; // we got all of them
      surfs.push_back(*ana::LoadFrom<Surface>(surfdir));
    }

    return std::make_unique<MedianSurface>(surfs);
  }

void ana::MedianSurface::SaveTo

(

TDirectory *

dir

)

const

Definition at line 123 of file MedianSurface.cxx.

  {
    TDirectory* tmp = gDirectory;
    dir->cd();
    TObjString("MedianSurface").Write("type");

    for(unsigned int i = 0; i < fThrows.size(); ++i){
      fThrows[i].SaveTo(dir->mkdir(TString::Format("surf%d", i)));
    }

    tmp->cd();
  }

Member Data Documentation

TH2F * ana::MedianSurface::fHistDn1

protected

Definition at line 22 of file MedianSurface.h.

TH2F * ana::MedianSurface::fHistDn2

protected

Definition at line 22 of file MedianSurface.h.

TH2F* ana::MedianSurface::fHistUp1

protected

Definition at line 22 of file MedianSurface.h.

TH2F * ana::MedianSurface::fHistUp2

protected

Definition at line 22 of file MedianSurface.h.

std::vector<Surface> ana::MedianSurface::fThrows

protected

Definition at line 20 of file MedianSurface.h.

---

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

• MedianSurface.h
• MedianSurface.cxx