HCovarianceVSxy Class Reference

#include <Histograms.h>

Inherits Histograms.

Public Member Functions
void	Clear () override
void	Fill (const double &x, const double &y, const double &a, const double &b) override
double	Get (const double &x, const double &y) const
	HCovarianceVSxy (const TString &name, const TString &title, const int totBinsX, const double &xMin, const double &xMax, const int totBinsY, const double &yMin, const double &yMax, TDirectory *dir=nullptr, bool varianceCheck=false)
	HCovarianceVSxy (TFile *inputFile, const TString &name, const TString &dirName)
	Contructor to read histograms from file. More...
void	Write () override
	~HCovarianceVSxy () override

Protected Member Functions
int	getXindex (const double &x) const
int	getYindex (const double &y) const

Protected Attributes
Covariance **	covariances_
double	deltaX_
double	deltaY_
TH2D *	histoCovariance_
TH1D ***	histoVarianceCheck_
bool	readMode_
int	totBinsX_
int	totBinsY_
int	totBinsZ_
bool	varianceCheck_
double	xMin_
double	yMin_

Detailed Description

This class can be used to compute the covariance of two variables with respect to other two variables (to see e.g. how does the covariance of ptVSphi vary with respect to (pt,eta).

Definition at line 1906 of file Histograms.h.

Constructor & Destructor Documentation

HCovarianceVSxy() [1/2]

HCovarianceVSxy::HCovarianceVSxy ( const TString &  name, const TString &  title, const int  totBinsX, const double &  xMin, const double &  xMax, const int  totBinsY, const double &  yMin, const double &  yMax, TDirectory *  dir = nullptr, bool  varianceCheck = false  )

inline

Definition at line 1908 of file Histograms.h.

      : totBinsX_(totBinsX),
        totBinsY_(totBinsY),
        xMin_(xMin),
        deltaX_(xMax - xMin),
        yMin_(yMin),
        deltaY_(yMax - yMin),
        readMode_(false),
        varianceCheck_(varianceCheck) {
    name_ = name;
    histoDir_ = dir;
    histoCovariance_ = new TH2D(name + "Covariance", title + " covariance", totBinsX, xMin, xMax, totBinsY, yMin, yMax);
 
    covariances_ = new Covariance*[totBinsX];
    for (int i = 0; i < totBinsX; ++i) {
      covariances_[i] = new Covariance[totBinsY];
    }
    if (varianceCheck_) {
      histoVarianceCheck_ = new TH1D**[totBinsX_];
      for (int i = 0; i < totBinsX_; ++i) {
        histoVarianceCheck_[i] = new TH1D*[totBinsY_];
        for (int j = 0; j < totBinsY_; ++j) {
          std::stringstream namei;
          std::stringstream namej;
          namei << i;
          namej << j;
          histoVarianceCheck_[i][j] = new TH1D(name + "_" + namei.str() + "_" + namej.str(), name, 10000, -1, 1);
        }
      }
    }
  }

References covariances_, DeadROC_duringRun::dir, histoCovariance_, histoVarianceCheck_, mps_fire::i, dqmiolumiharvest::j, Skims_PA_cff::name, overlapproblemtsosanalyzer_cfi::title, totBinsX_, totBinsY_, varianceCheck_, multiplicitycorr_cfi::xMax, photonAnalyzer_cfi::xMin, multiplicitycorr_cfi::yMax, and photonAnalyzer_cfi::yMin.

HCovarianceVSxy() [2/2]

HCovarianceVSxy::HCovarianceVSxy ( TFile *  inputFile, const TString &  name, const TString &  dirName  )

inline

Contructor to read histograms from file.

Definition at line 1949 of file Histograms.h.

                                                                                 : readMode_(true) {
    histoDir_ = (TDirectory*)(inputFile->Get(dirName.Data()));
    if (histoDir_ == nullptr) {
      std::cout << "Error: directory not found" << std::endl;
      exit(0);
    }
    histoCovariance_ = (TH2D*)(histoDir_->Get(name));
    totBinsX_ = histoCovariance_->GetNbinsX();
    xMin_ = histoCovariance_->GetXaxis()->GetBinLowEdge(1);
    deltaX_ = histoCovariance_->GetXaxis()->GetBinUpEdge(totBinsX_) - xMin_;
    totBinsY_ = histoCovariance_->GetNbinsY();
    yMin_ = histoCovariance_->GetYaxis()->GetBinLowEdge(1);
    deltaY_ = histoCovariance_->GetYaxis()->GetBinUpEdge(totBinsY_) - yMin_;
  }

References gather_cfg::cout, deltaX_, deltaY_, TrackerOfflineValidation_Dqm_cff::dirName, beamvalidation::exit(), histoCovariance_, dtResolutionTest_cfi::inputFile, Skims_PA_cff::name, totBinsX_, totBinsY_, xMin_, and yMin_.

~HCovarianceVSxy()

HCovarianceVSxy::~HCovarianceVSxy ( )

inlineoverride

Definition at line 1964 of file Histograms.h.

                              {
    delete histoCovariance_;
    // Free covariances
    for (int i = 0; i < totBinsX_; ++i) {
      delete[] covariances_[i];
    }
    delete[] covariances_;
    // Free variance check histograms
    if (varianceCheck_) {
      for (int i = 0; i < totBinsX_; ++i) {
        for (int j = 0; j < totBinsY_; ++j) {
          delete histoVarianceCheck_[i][j];
        }
        delete[] histoVarianceCheck_[i];
      }
      delete[] histoVarianceCheck_;
    }
  }

References covariances_, histoCovariance_, histoVarianceCheck_, mps_fire::i, dqmiolumiharvest::j, totBinsX_, totBinsY_, and varianceCheck_.

Member Function Documentation

Clear()

void HCovarianceVSxy::Clear ( )

inlineoverride

Definition at line 2059 of file Histograms.h.

                        {
    histoCovariance_->Clear();
    if (varianceCheck_) {
      for (int i = 0; i < totBinsX_; ++i) {
        for (int j = 0; j < totBinsY_; ++j) {
          histoVarianceCheck_[i][j]->Clear();
        }
      }
    }
  }

References histoCovariance_, histoVarianceCheck_, mps_fire::i, dqmiolumiharvest::j, totBinsX_, totBinsY_, and varianceCheck_.

Fill()

void HCovarianceVSxy::Fill ( const double &  x, const double &  y, const double &  a, const double &  b  )

inlineoverride

x and y should be the variables VS which we are computing the covariance (pt and eta) a and b should be the variables OF which we are computing the covariance

Definition at line 1987 of file Histograms.h.

                                                                                         {
    // Need to convert the (x,y) values to the array indeces
    int xIndex = getXindex(x);
    int yIndex = getYindex(y);
    // Only fill values if they are in the selected range
    if (0 <= xIndex && xIndex < totBinsX_ && 0 <= yIndex && yIndex < totBinsY_) {
      // if( TString(histoCovariance_->GetName()).Contains("CovarianceCotgTheta_Covariance") )
      covariances_[xIndex][yIndex].fill(a, b);
      // Should be used with the variance, in which case a==b
      if (varianceCheck_)
        histoVarianceCheck_[xIndex][yIndex]->Fill(a);
    }
  }

References a, b, covariances_, Covariance::fill(), getXindex(), getYindex(), histoVarianceCheck_, totBinsX_, totBinsY_, varianceCheck_, x, and y.

Referenced by ResolutionAnalyzer::analyze().

Get()

double HCovarianceVSxy::Get ( const double &  x, const double &  y  ) const

inline

Definition at line 2002 of file Histograms.h.

                                                     {
    // Need to convert the (x,y) values to the array indeces
    int xIndex = getXindex(x);
    int yIndex = getYindex(y);
    // If the values exceed the limits of the histogram, return the border values
    if (xIndex < 0)
      xIndex = 0;
    if (xIndex >= totBinsX_)
      xIndex = totBinsX_ - 1;
    if (yIndex < 0)
      yIndex = 0;
    if (yIndex >= totBinsY_)
      yIndex = totBinsY_ - 1;
    return histoCovariance_->GetBinContent(xIndex + 1, yIndex + 1);
  }

References getXindex(), getYindex(), histoCovariance_, totBinsX_, totBinsY_, x, and y.

getXindex()

int HCovarianceVSxy::getXindex ( const double &  x ) const

inlineprotected

Definition at line 2071 of file Histograms.h.

{ return int((x - xMin_) / deltaX_ * totBinsX_); }

References deltaX_, createfilelist::int, totBinsX_, x, and xMin_.

Referenced by Fill(), and Get().

getYindex()

int HCovarianceVSxy::getYindex ( const double &  y ) const

inlineprotected

Definition at line 2072 of file Histograms.h.

{ return int((y - yMin_) / deltaY_ * totBinsY_); }

References deltaY_, createfilelist::int, totBinsY_, y, and yMin_.

Referenced by Fill(), and Get().

Write()

void HCovarianceVSxy::Write ( )

inlineoverride

Definition at line 2018 of file Histograms.h.

                        {
    if (!readMode_) {
      std::cout << "writing: " << histoCovariance_->GetName() << std::endl;
      for (int xBin = 0; xBin < totBinsX_; ++xBin) {
        for (int yBin = 0; yBin < totBinsY_; ++yBin) {
          double covariance = covariances_[xBin][yBin].covariance();
          // Histogram bins start from 1
          // std::cout << "covariance["<<xBin<<"]["<<yBin<<"] with N = "<<covariances_[xBin][yBin].getN()<<" is: " << covariance << std::endl;
          histoCovariance_->SetBinContent(xBin + 1, yBin + 1, covariance);
        }
      }
      if (histoDir_ != nullptr)
        histoDir_->cd();
      TCanvas canvas(TString(histoCovariance_->GetName()) + "_canvas",
                     TString(histoCovariance_->GetTitle()) + " canvas",
                     1000,
                     800);
      canvas.Divide(2);
      canvas.cd(1);
      histoCovariance_->Draw("lego");
      canvas.cd(2);
      histoCovariance_->Draw("surf5");
      canvas.Write();
      histoCovariance_->Write();
 
      TDirectory* binsDir = nullptr;
      if (varianceCheck_) {
        if (histoDir_ != nullptr) {
          histoDir_->cd();
          if (binsDir == nullptr)
            binsDir = histoDir_->mkdir(name_ + "Bins");
          binsDir->cd();
        }
        for (int xBin = 0; xBin < totBinsX_; ++xBin) {
          for (int yBin = 0; yBin < totBinsY_; ++yBin) {
            histoVarianceCheck_[xBin][yBin]->Write();
          }
        }
      }
    }
  }

References svgfig::canvas(), gather_cfg::cout, Covariance::covariance(), covariances_, histoCovariance_, histoVarianceCheck_, readMode_, totBinsX_, totBinsY_, varianceCheck_, photonAnalyzer_cfi::xBin, and photonAnalyzer_cfi::yBin.

Referenced by ResolutionAnalyzer::writeHistoMap().

Member Data Documentation

covariances_

Covariance** HCovarianceVSxy::covariances_

protected

Definition at line 2074 of file Histograms.h.

Referenced by Fill(), HCovarianceVSxy(), Write(), and ~HCovarianceVSxy().

deltaX_

double HCovarianceVSxy::deltaX_

protected

Definition at line 2076 of file Histograms.h.

Referenced by getXindex(), and HCovarianceVSxy().

deltaY_

double HCovarianceVSxy::deltaY_

protected

Definition at line 2076 of file Histograms.h.

Referenced by getYindex(), and HCovarianceVSxy().

histoCovariance_

TH2D* HCovarianceVSxy::histoCovariance_

protected

Definition at line 2073 of file Histograms.h.

Referenced by Clear(), Get(), HCovarianceVSxy(), Write(), and ~HCovarianceVSxy().

histoVarianceCheck_

TH1D*** HCovarianceVSxy::histoVarianceCheck_

protected

Definition at line 2079 of file Histograms.h.

Referenced by Clear(), Fill(), HCovarianceVSxy(), Write(), and ~HCovarianceVSxy().

readMode_

bool HCovarianceVSxy::readMode_

protected

Definition at line 2077 of file Histograms.h.

Referenced by Write().

totBinsX_

int HCovarianceVSxy::totBinsX_

protected

Definition at line 2075 of file Histograms.h.

Referenced by Clear(), Fill(), Get(), getXindex(), HCovarianceVSxy(), Write(), and ~HCovarianceVSxy().

totBinsY_

int HCovarianceVSxy::totBinsY_

protected

Definition at line 2075 of file Histograms.h.

Referenced by Clear(), Fill(), Get(), getYindex(), HCovarianceVSxy(), Write(), and ~HCovarianceVSxy().

totBinsZ_

int HCovarianceVSxy::totBinsZ_

protected

Definition at line 2075 of file Histograms.h.

varianceCheck_

bool HCovarianceVSxy::varianceCheck_

protected

Definition at line 2078 of file Histograms.h.

Referenced by Clear(), Fill(), HCovarianceVSxy(), Write(), and ~HCovarianceVSxy().

xMin_

double HCovarianceVSxy::xMin_

protected

Definition at line 2076 of file Histograms.h.

Referenced by getXindex(), and HCovarianceVSxy().

yMin_

double HCovarianceVSxy::yMin_

protected

Definition at line 2076 of file Histograms.h.

Referenced by getYindex(), and HCovarianceVSxy().