#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_

Additional Inherited Members
Public Attributes inherited from Histograms
dqm::reco::MonitorElement *	lumiVsLS

dqm::reco::MonitorElement *	numberOfPixelClustersVsLS

dqm::reco::MonitorElement *	numberOfPixelClustersVsLumi

dqm::reco::MonitorElement *	pixelLumiVsLS

dqm::reco::MonitorElement *	pixelLumiVsLumi

dqm::reco::MonitorElement *	puVsLS

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 1908 of file Histograms.h.

Constructor & Destructor Documentation

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 1910 of file Histograms.h.

References covariances_, DeadROC_duringRun::dir, histoCovariance_, histoVarianceCheck_, mps_fire::i, dqmiolumiharvest::j, mergeVDriftHistosByStation::name, totBinsX_, totBinsY_, and varianceCheck_.

       : 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);
         }
       }
     }
   }

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

inline

Contructor to read histograms from file.

Definition at line 1951 of file Histograms.h.

References gather_cfg::cout, deltaX_, deltaY_, beamvalidation::exit(), histoCovariance_, totBinsX_, totBinsY_, xMin_, and yMin_.

                                                                                  : 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_;
   }

HCovarianceVSxy::~HCovarianceVSxy ( )

inlineoverride

Definition at line 1966 of file Histograms.h.

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

                               {
     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_;
     }
   }

Member Function Documentation

void HCovarianceVSxy::Clear ( )

inlineoverride

Definition at line 2061 of file Histograms.h.

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

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

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 1989 of file Histograms.h.

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

Referenced by ResolutionAnalyzer::analyze().

                                                                                          {
     // 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);
     }
   }

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

inline

Definition at line 2004 of file Histograms.h.

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

                                                      {
     // 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);
   }

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

inlineprotected

Definition at line 2073 of file Histograms.h.

References deltaX_, totBinsX_, and xMin_.

Referenced by Fill(), and Get().

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

HCovarianceVSxy::xMin_

double xMin_

Definition: Histograms.h:2078

HCovarianceVSxy::deltaX_

double deltaX_

Definition: Histograms.h:2078

HCovarianceVSxy::totBinsX_

int totBinsX_

Definition: Histograms.h:2077

DDAxes::x

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

inlineprotected

Definition at line 2074 of file Histograms.h.

References deltaY_, totBinsY_, and yMin_.

Referenced by Fill(), and Get().

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

HCovarianceVSxy::deltaY_

double deltaY_

Definition: Histograms.h:2078

HCovarianceVSxy::totBinsY_

int totBinsY_

Definition: Histograms.h:2077

HCovarianceVSxy::yMin_

double yMin_

Definition: Histograms.h:2078

DDAxes::y

void HCovarianceVSxy::Write ( )

inlineoverride

Definition at line 2020 of file Histograms.h.

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

Referenced by ResolutionAnalyzer::writeHistoMap().

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