#include <Histograms.h>

Inheritance diagram for HCovarianceVSxy:

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=0, bool varianceCheck=false)

	HCovarianceVSxy (TFile *inputFile, const TString &name, const TString &dirName)
	Contructor to read histograms from file. More...

void	Write () override

	~HCovarianceVSxy () override

Public Member Functions inherited from Histograms
virtual void	Fill (const reco::Particle::LorentzVector &p1, const reco::Particle::LorentzVector &p2)

virtual void	Fill (const reco::Particle::LorentzVector &p1, const reco::Particle::LorentzVector &p2, const int charge, const double &weight=1.)

virtual void	Fill (const CLHEP::HepLorentzVector &momentum1, const CLHEP::HepLorentzVector &momentum2)

virtual void	Fill (const CLHEP::HepLorentzVector &momentum1, const CLHEP::HepLorentzVector &momentum2, const int charge, const double &weight=1.)

virtual void	Fill (const CLHEP::HepLorentzVector &p1, const reco::Particle::LorentzVector &p2)

virtual void	Fill (const reco::Particle::LorentzVector &p4, const double &weight=1.)

virtual void	Fill (const reco::Particle::LorentzVector &p4, const int charge, const double &weight=1.)

virtual void	Fill (const CLHEP::HepLorentzVector &momentum, const int charge, const double &weight=1.)

virtual void	Fill (const reco::Particle::LorentzVector &p4, const double &resValue, const int charge)

virtual void	Fill (const reco::Particle::LorentzVector &p4, const double &genValue, const double recValue, const int charge)

virtual void	Fill (const CLHEP::HepLorentzVector &p, const double &likeValue)

virtual void	Fill (const unsigned int number)

virtual void	Fill (const reco::Particle::LorentzVector &recoP1, const int charge1, const reco::Particle::LorentzVector &genP1, const reco::Particle::LorentzVector &recoP2, const int charge2, const reco::Particle::LorentzVector &genP2, const double &recoMass, const double &genMass)

virtual void	Fill (const reco::Particle::LorentzVector &recoP1, const int charge1, const reco::Particle::LorentzVector &recoP2, const int charge2, const double &recoMass, const double &genMass)

virtual void	Fill (const reco::Particle::LorentzVector &recoP1, const reco::Particle::LorentzVector &genP1, const reco::Particle::LorentzVector &recoP2, const reco::Particle::LorentzVector &genP2)

virtual void	Fill (const double &x, const double &y)

virtual void	Fill (const reco::Particle::LorentzVector &p41, const reco::Particle::LorentzVector &p42, const reco::Particle::LorentzVector &p4Res, const double &weight=1.)

virtual void	Fill (const CLHEP::HepLorentzVector &momentum1, const CLHEP::HepLorentzVector &momentum2, const CLHEP::HepLorentzVector &momentumRes, const double &weight=1.)

virtual double	Get (const reco::Particle::LorentzVector &recoP1, const TString &covarianceName)

virtual TString	GetName ()

	Histograms ()

	Histograms (const TString &name)

	Histograms (TFile *outputFile, const TString &name)

virtual void	SetWeight (double weight)

virtual	~Histograms ()

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_

Protected Attributes inherited from Histograms
TDirectory *	histoDir_

TString	name_

TFile *	outputFile_

double	theWeight_

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 1818 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 = `0`,
		bool	varianceCheck = `false`
	)

inline

Definition at line 1821 of file Histograms.h.

References dir, Histograms::histoDir_, mps_fire::i, dataset::name, and Histograms::name_.

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

References gather_cfg::cout, cmsRelvalreport::exit, and Histograms::histoDir_.

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

References mps_fire::i.

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

inlineoverridevirtual

Implements Histograms.

Definition at line 1954 of file Histograms.h.

References mps_fire::i.

                         {
     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
	)

inlineoverridevirtual

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

Reimplemented from Histograms.

Definition at line 1893 of file Histograms.h.

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

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

inlineprotected

Definition at line 1965 of file Histograms.h.

References createfilelist::int.

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

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

inlineprotected

Definition at line 1968 of file Histograms.h.

References createfilelist::int.

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

void HCovarianceVSxy::Write ( )

inlineoverridevirtual

Implements Histograms.

Definition at line 1918 of file Histograms.h.

References svgfig::canvas(), gather_cfg::cout, Histograms::histoDir_, Histograms::name_, 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();
           }
         }
       }
     }
   }

Member Data Documentation

Covariance** HCovarianceVSxy::covariances_

protected

Definition at line 1972 of file Histograms.h.

double HCovarianceVSxy::deltaX_

protected

Definition at line 1974 of file Histograms.h.

double HCovarianceVSxy::deltaY_

protected

Definition at line 1974 of file Histograms.h.

TH2D* HCovarianceVSxy::histoCovariance_

protected

Definition at line 1971 of file Histograms.h.

TH1D*** HCovarianceVSxy::histoVarianceCheck_

protected

Definition at line 1977 of file Histograms.h.

bool HCovarianceVSxy::readMode_

protected

Definition at line 1975 of file Histograms.h.

int HCovarianceVSxy::totBinsX_

protected

Definition at line 1973 of file Histograms.h.

int HCovarianceVSxy::totBinsY_

protected

Definition at line 1973 of file Histograms.h.

int HCovarianceVSxy::totBinsZ_

protected

Definition at line 1973 of file Histograms.h.

bool HCovarianceVSxy::varianceCheck_

protected

Definition at line 1976 of file Histograms.h.

double HCovarianceVSxy::xMin_

protected

Definition at line 1974 of file Histograms.h.

double HCovarianceVSxy::yMin_

protected

Definition at line 1974 of file Histograms.h.

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation