HFunctionResolution Class Reference

#include <Histograms.h>

Inheritance diagram for HFunctionResolution:

Public Member Functions
void	Clear () override
void	Fill (const reco::Particle::LorentzVector &p4, const double &resValue, const int charge) override
	HFunctionResolution (TFile *outputFile, const TString &name, const double &ptMax=100, const int totBinsY=300)
void	Write () override
	~HFunctionResolution () override

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

Protected Attributes
double	deltaX_
double	deltaY_
TH1F *	hReso_
TProfile *	hResoVSEta_prof_
TProfile *	hResoVSPhi_prof_
TProfile *	hResoVSPhiMinus_prof_
TProfile *	hResoVSPhiPlus_prof_
TProfile *	hResoVSPt_Bar_prof_
TProfile *	hResoVSPt_Endc_17_prof_
TProfile *	hResoVSPt_Endc_20_prof_
TProfile *	hResoVSPt_Endc_24_prof_
TProfile *	hResoVSPt_Ovlap_prof_
TProfile *	hResoVSPt_prof_
TH2F *	hResoVSPtEta_
int **	resoCount_
double **	resoVsPtEta_
int	totBinsX_
int	totBinsY_
double	xMin_
double	yMin_

Detailed Description

This histogram class fills a TProfile with the resolution evaluated from the resolution functions for single muon quantities. The resolution functions are used by MuScleFit to evaluate the mass resolution, which is the value seen by minuit and through it, corrections are evaluated.
In the end we will compare the histograms filled by this class (from the resolution function, reflecting the parameters changes done by minuit) with those filled comparing recoMuons with genMuons (the real resolutions).

Definition at line 1471 of file Histograms.h.

Constructor & Destructor Documentation

HFunctionResolution::HFunctionResolution ( TFile *  outputFile, const TString &  name, const double &  ptMax = 100, const int  totBinsY = 300  )

inline

Definition at line 1474 of file Histograms.h.

References mps_fire::i, dataset::name, AlignmentTrackSelector_cfi::ptMax, anotherprimaryvertexanalyzer_cfi::xMax, and CMSBoostedTauSeedingParameters_cfi::yMax.

                                                                                                                      : Histograms(outputFile, name) {
    name_ = name;
    totBinsX_ = 300;
    totBinsY_ = totBinsY;
    xMin_ = 0.;
    yMin_ = -3.0;
    double xMax = ptMax;
    double yMax = 3.0;
    deltaX_ = xMax - xMin_;
    deltaY_ = yMax - yMin_;
    hReso_        = new TH1F( name+"_Reso", "resolution", 1000, 0, 1 );
    hResoVSPtEta_ = new TH2F( name+"_ResoVSPtEta", "resolution vs pt and #eta", totBinsX_, xMin_, xMax, totBinsY_, yMin_, yMax );
    // Create and initialize the resolution arrays
    resoVsPtEta_  = new double*[totBinsX_];
    resoCount_    = new int*[totBinsX_];
    for( int i=0; i<totBinsX_; ++i ) {
      resoVsPtEta_[i] = new double[totBinsY_];
      resoCount_[i]   = new int[totBinsY_];
      for( int j=0; j<totBinsY_; ++j ) {
        resoVsPtEta_[i][j] = 0;
        resoCount_[i][j] = 0;
      }
    }
    hResoVSPt_prof_       = new TProfile( name+"_ResoVSPt_prof", "resolution VS pt", totBinsX_, xMin_, xMax, yMin_, yMax);
    hResoVSPt_Bar_prof_   = new TProfile( name+"_ResoVSPt_Bar_prof", "resolution VS pt Barrel", totBinsX_, xMin_, xMax, yMin_, yMax);
    hResoVSPt_Endc_17_prof_  = new TProfile( name+"_ResoVSPt_Endc_1.7_prof", "resolution VS pt Endcap (1.4<eta<1.7)", totBinsX_, xMin_, xMax, yMin_, yMax);
    hResoVSPt_Endc_20_prof_  = new TProfile( name+"_ResoVSPt_Endc_2.0_prof", "resolution VS pt Endcap (1.7<eta<2.0)", totBinsX_, xMin_, xMax, yMin_, yMax);
    hResoVSPt_Endc_24_prof_  = new TProfile( name+"_ResoVSPt_Endc_2.4_prof", "resolution VS pt Endcap (2.0<eta<2.4)", totBinsX_, xMin_, xMax, yMin_, yMax);
    hResoVSPt_Ovlap_prof_ = new TProfile( name+"_ResoVSPt_Ovlap_prof", "resolution VS pt Overlap", totBinsX_, xMin_, xMax, yMin_, yMax);
    hResoVSEta_prof_      = new TProfile( name+"_ResoVSEta_prof", "resolution VS eta", totBinsY_, yMin_, yMax, 0, 1);
    //hResoVSTheta_prof_    = new TProfile( name+"_ResoVSTheta_prof", "resolution VS theta", 30, 0, TMath::Pi(), 0, 1);
    hResoVSPhiPlus_prof_  = new TProfile( name+"_ResoVSPhiPlus_prof", "resolution VS phi mu+", 16, -3.2, 3.2, 0, 1);
    hResoVSPhiMinus_prof_ = new TProfile( name+"_ResoVSPhiMinus_prof", "resolution VS phi mu-", 16, -3.2, 3.2, 0, 1);
    hResoVSPhi_prof_      = new TProfile( name+"_ResoVSPhi_prof", "resolution VS phi", 16, -3.2, 3.2, -1, 1);
  }

HFunctionResolution::~HFunctionResolution ( )

inlineoverride

Definition at line 1509 of file Histograms.h.

References mps_fire::i.

                                  {
    delete hReso_;
    delete hResoVSPtEta_;
    // Free the resolution arrays
    for( int i=0; i<totBinsX_; ++i ) {
      delete[] resoVsPtEta_[i];
      delete[] resoCount_[i];
    }
    delete[] resoVsPtEta_;
    delete[] resoCount_;
    // Free the profile histograms
    delete hResoVSPt_prof_;
    delete hResoVSPt_Bar_prof_;
    delete hResoVSPt_Endc_17_prof_;
    delete hResoVSPt_Endc_20_prof_;
    delete hResoVSPt_Endc_24_prof_;
    delete hResoVSPt_Ovlap_prof_;
    delete hResoVSEta_prof_;
    //delete hResoVSTheta_prof_;
    delete hResoVSPhiPlus_prof_;
    delete hResoVSPhiMinus_prof_;
    delete hResoVSPhi_prof_;
  }

Member Function Documentation

void HFunctionResolution::Clear ( )

inlineoverride

Definition at line 1612 of file Histograms.h.

                        {
    hReso_->Clear();
    hResoVSPtEta_->Clear();
    hResoVSPt_prof_->Clear();
    hResoVSPt_Bar_prof_->Clear();
    hResoVSPt_Endc_17_prof_->Clear();
    hResoVSPt_Endc_20_prof_->Clear();
    hResoVSPt_Endc_24_prof_->Clear();
    hResoVSPt_Ovlap_prof_->Clear();
    hResoVSEta_prof_->Clear();
    //hResoVSTheta_prof_->Clear();
    hResoVSPhiPlus_prof_->Clear();
    hResoVSPhiMinus_prof_->Clear();
    hResoVSPhi_prof_->Clear();
  }

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

inlineoverride

Definition at line 1532 of file Histograms.h.

Referenced by HFunctionResolutionVarianceCheck::Fill().

                                                                                                      {
    if( resValue != resValue ) return;
    hReso_->Fill(resValue);

    // Fill the arrays with the resolution value and count
    int xIndex = getXindex(p4.Pt());
    int yIndex = getYindex(p4.Eta());
    if ( 0 <= xIndex && xIndex < totBinsX_ && 0 <= yIndex && yIndex < totBinsY_ ) {
      resoVsPtEta_[xIndex][yIndex] += resValue;
      // ATTENTION: we count only for positive muons because we are considering di-muon resonances
      // and we use this counter to compute the mean in the end. The resoVsPtEta value is filled with each muon,
      // but they must be considered independently (analogous to a factor 2) so in the end we would have
      // to divide by N/2, that is why we only increase the counter for half the muons.
      // if( charge > 0 )
      // No more. Changing it here influences also other uses of this class. The macro FunctionTerms.cc
      // multiplies the terms by the 2 factor.

      resoCount_[xIndex][yIndex] += 1;

      // hResoVSPtEta->Fill(p4.Pt(), p4.Eta(), resValue);
      hResoVSPt_prof_->Fill(p4.Pt(),resValue);
      if(fabs(p4.Eta())<=0.9)
    hResoVSPt_Bar_prof_->Fill(p4.Pt(),resValue);
      else if(fabs(p4.Eta())>0.9 && fabs(p4.Eta())<=1.4 )
    hResoVSPt_Ovlap_prof_->Fill(p4.Pt(),resValue);
      else if(fabs(p4.Eta())>1.4 && fabs(p4.Eta())<=1.7 )
    hResoVSPt_Endc_17_prof_->Fill(p4.Pt(),resValue);
      else if(fabs(p4.Eta())>1.7 && fabs(p4.Eta())<=2.0 )
    hResoVSPt_Endc_20_prof_->Fill(p4.Pt(),resValue);
      else
        hResoVSPt_Endc_24_prof_->Fill(p4.Pt(),resValue);
      hResoVSEta_prof_->Fill(p4.Eta(),resValue);
      //hResoVSTheta_prof_->Fill(p4.Theta(),resValue);
      if(charge>0)
        hResoVSPhiPlus_prof_->Fill(p4.Phi(),resValue);
      else if(charge<0)
        hResoVSPhiMinus_prof_->Fill(p4.Phi(),resValue);
      hResoVSPhi_prof_->Fill(p4.Phi(),resValue);
    }
  }

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

inlineprotected

Definition at line 1629 of file Histograms.h.

References createfilelist::int.

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

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

inlineprotected

Definition at line 1632 of file Histograms.h.

References createfilelist::int.

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

void HFunctionResolution::Write ( )

inlineoverride

Definition at line 1573 of file Histograms.h.

References svgfig::canvas(), mps_fire::i, and N.

Referenced by HFunctionResolutionVarianceCheck::Write().

                        {
    if(histoDir_ != nullptr) histoDir_->cd();

    hReso_->Write();

    for( int i=0; i<totBinsX_; ++i ) {
      for( int j=0; j<totBinsY_; ++j ) {
        int N = resoCount_[i][j];
        // Fill with the mean value
        if( N != 0 ) hResoVSPtEta_->SetBinContent( i+1, j+1, resoVsPtEta_[i][j]/N );
        else hResoVSPtEta_->SetBinContent( i+1, j+1, 0 );
      }
    }
    hResoVSPtEta_->Write();

    hResoVSPt_prof_->Write();
    hResoVSPt_Bar_prof_->Write();
    hResoVSPt_Endc_17_prof_->Write();
    hResoVSPt_Endc_20_prof_->Write();
    hResoVSPt_Endc_24_prof_->Write();
    hResoVSPt_Ovlap_prof_->Write();
    hResoVSEta_prof_->Write();
    //hResoVSTheta_prof_->Write();
    hResoVSPhiMinus_prof_->Write();
    hResoVSPhiPlus_prof_->Write();
    hResoVSPhi_prof_->Write();

    TCanvas canvas(TString(hResoVSPtEta_->GetName())+"_canvas", TString(hResoVSPtEta_->GetTitle())+" canvas", 1000, 800);
    canvas.Divide(2);
    canvas.cd(1);
    hResoVSPtEta_->Draw("lego");
    canvas.cd(2);
    hResoVSPtEta_->Draw("surf5");
    canvas.Write();
    hResoVSPtEta_->Write();

    outputFile_->cd();
  }

Member Data Documentation

double HFunctionResolution::deltaX_

protected

Definition at line 1652 of file Histograms.h.

double HFunctionResolution::deltaY_

protected

Definition at line 1652 of file Histograms.h.

TH1F* HFunctionResolution::hReso_

protected

Definition at line 1635 of file Histograms.h.

TProfile* HFunctionResolution::hResoVSEta_prof_

protected

Definition at line 1645 of file Histograms.h.

TProfile* HFunctionResolution::hResoVSPhi_prof_

protected

Definition at line 1649 of file Histograms.h.

TProfile* HFunctionResolution::hResoVSPhiMinus_prof_

protected

Definition at line 1647 of file Histograms.h.

TProfile* HFunctionResolution::hResoVSPhiPlus_prof_

protected

Definition at line 1648 of file Histograms.h.

TProfile* HFunctionResolution::hResoVSPt_Bar_prof_

protected

Definition at line 1640 of file Histograms.h.

TProfile* HFunctionResolution::hResoVSPt_Endc_17_prof_

protected

Definition at line 1641 of file Histograms.h.

TProfile* HFunctionResolution::hResoVSPt_Endc_20_prof_

protected

Definition at line 1642 of file Histograms.h.

TProfile* HFunctionResolution::hResoVSPt_Endc_24_prof_

protected

Definition at line 1643 of file Histograms.h.

TProfile* HFunctionResolution::hResoVSPt_Ovlap_prof_

protected

Definition at line 1644 of file Histograms.h.

TProfile* HFunctionResolution::hResoVSPt_prof_

protected

Definition at line 1639 of file Histograms.h.

TH2F* HFunctionResolution::hResoVSPtEta_

protected

Definition at line 1636 of file Histograms.h.

int** HFunctionResolution::resoCount_

protected

Definition at line 1638 of file Histograms.h.

double** HFunctionResolution::resoVsPtEta_

protected

Definition at line 1637 of file Histograms.h.

int HFunctionResolution::totBinsX_

protected

Definition at line 1650 of file Histograms.h.

int HFunctionResolution::totBinsY_

protected

Definition at line 1650 of file Histograms.h.

double HFunctionResolution::xMin_

protected

Definition at line 1651 of file Histograms.h.

double HFunctionResolution::yMin_

protected

Definition at line 1651 of file Histograms.h.