HResolutionVSPart Class Reference

A set of histograms for resolution. More...

#include <Histograms.h>

Inherits Histograms.

Public Member Functions
void	Clear () override
void	Fill (const reco::Particle::LorentzVector &p4, const double &resValue, const int charge) override
	HResolutionVSPart (TFile *outputFile, const TString &name, const double maxPt=100, const double &yMinEta=0., const double &yMaxEta=2., const double &yMinPt=0., const double &yMaxPt=2., const bool doProfiles=false)
	HResolutionVSPart (const TString &name, TFile *file)
void	Write () override
	~HResolutionVSPart () override

Public Attributes
bool	doProfiles_
TH1F *	hAbsReso_
TH2F *	hAbsResoVSEta_
TH2F *	hAbsResoVSPhi_
TH2F *	hAbsResoVSPt_
TH1F *	hReso_
TH2F *	hResoVSEta_
TProfile *	hResoVSEta_prof_
TProfile *	hResoVSPhi_prof_
TH2F *	hResoVSPhiMinus_
TH2F *	hResoVSPhiPlus_
TH2F *	hResoVSPt_
TH2F *	hResoVSPt_Bar_
TProfile *	hResoVSPt_Bar_prof_
TH2F *	hResoVSPt_Endc_17_
TProfile *	hResoVSPt_Endc_17_prof_
TH2F *	hResoVSPt_Endc_20_
TProfile *	hResoVSPt_Endc_20_prof_
TH2F *	hResoVSPt_Endc_24_
TProfile *	hResoVSPt_Endc_24_prof_
TH2F *	hResoVSPt_Ovlap_
TProfile *	hResoVSPt_Ovlap_prof_
TProfile *	hResoVSPt_prof_
TH2F *	hResoVSPtEta_
TH2F *	hResoVSTheta_

Detailed Description

A set of histograms for resolution.

Definition at line 1139 of file Histograms.h.

Constructor & Destructor Documentation

HResolutionVSPart::HResolutionVSPart ( TFile *  outputFile, const TString &  name, const double  maxPt = 100, const double &  yMinEta = 0., const double &  yMaxEta = 2., const double &  yMinPt = 0., const double &  yMaxPt = 2., const bool  doProfiles = false  )

inline

Definition at line 1142 of file Histograms.h.

References MuonErrorMatrixAnalyzer_cfi::maxPt, and Pi.

                                                   : Histograms(outputFile, name), doProfiles_(doProfiles)
  {
    // Kinematical variables

    // hReso           = new TH1F (name+"_Reso", "resolution", 4000, -1, 1);
    // hResoVSPtEta    = new TH2F (name+"_ResoVSPtEta", "resolution VS pt and #eta", 200, 0, 200, 60, -3, 3);
    // hResoVSPt       = new TH2F (name+"_ResoVSPt", "resolution VS pt", 200, 0, 200, 8000, -1, 1);
    // //hResoVSPt_prof  = new TProfile (name+"_ResoVSPt_prof", "resolution VS pt", 100, 0, 200, -1, 1);
    // hResoVSEta      = new TH2F (name+"_ResoVSEta", "resolution VS eta", 60, -3, 3, 8000, yMinEta, yMaxEta);
    // hResoVSTheta    = new TH2F (name+"_ResoVSTheta", "resolution VS theta", 30, 0, TMath::Pi(), 8000, -1, 1);
    // //hResoVSEta_prof = new TProfile (name+"_ResoVSEta_prof", "resolution VS eta", 10, -2.5, 2.5, -1, 1);
    // hResoVSPhiPlus  = new TH2F (name+"_ResoVSPhiPlus", "resolution VS phi mu+", 14, -3.2, 3.2, 8000, -1, 1);
    // hResoVSPhiMinus = new TH2F (name+"_ResoVSPhiMinus", "resolution VS phi mu-", 14, -3.2, 3.2, 8000, -1, 1);
    // //hResoVSPhi_prof = new TProfile (name+"_ResoVSPhi_prof", "resolution VS phi", 14, -3.2, 3.2, -1, 1);
    // hAbsReso        = new TH1F (name+"_AbsReso", "resolution", 100, 0, 1);
    // hAbsResoVSPt    = new TH2F (name+"_AbsResoVSPt", "Abs resolution VS pt", 200, 0, 500, 100, 0, 1);
    // hAbsResoVSEta   = new TH2F (name+"_AbsResoVSEta", "Abs resolution VS eta", 60, -3, 3, 100, 0, 1);
    // hAbsResoVSPhi   = new TH2F (name+"_AbsResoVSPhi", "Abs resolution VS phi", 14, -3.2, 3.2, 100, 0, 1);

    hReso_           = new TH1F( name+"_Reso", "resolution", 200, -1, 1 );
    hResoVSPtEta_    = new TH2F( name+"_ResoVSPtEta", "resolution VS pt and #eta", 200, 0, maxPt, 60, -3, 3 );
    hResoVSPt_       = new TH2F( name+"_ResoVSPt", "resolution VS pt", 200, 0, maxPt, 200, yMinPt, yMaxPt );
    hResoVSPt_Bar_   = new TH2F( name+"_ResoVSPt_Bar", "resolution VS pt Barrel", 200, 0, maxPt, 200, yMinPt, yMaxPt );
    hResoVSPt_Endc_17_  = new TH2F( name+"_ResoVSPt_Endc_1.7", "resolution VS pt Endcap (1.4<eta<1.7)", 200, 0, maxPt, 200, yMinPt, yMaxPt );
    hResoVSPt_Endc_20_  = new TH2F( name+"_ResoVSPt_Endc_2.0", "resolution VS pt Endcap (1.7<eta<2.0)", 200, 0, maxPt, 200, yMinPt, yMaxPt );
    hResoVSPt_Endc_24_  = new TH2F( name+"_ResoVSPt_Endc_2.4", "resolution VS pt Endcap (2.0<eta<2.4)", 200, 0, maxPt, 200, yMinPt, yMaxPt );
    hResoVSPt_Ovlap_ = new TH2F( name+"_ResoVSPt_Ovlap", "resolution VS pt overlap", 200, 0, maxPt, 200, yMinPt, yMaxPt );
    hResoVSEta_      = new TH2F( name+"_ResoVSEta", "resolution VS eta", 200, -3, 3, 200, yMinEta, yMaxEta );
    hResoVSTheta_    = new TH2F( name+"_ResoVSTheta", "resolution VS theta", 30, 0, TMath::Pi(), 200, yMinEta, yMaxEta );
    hResoVSPhiPlus_  = new TH2F( name+"_ResoVSPhiPlus", "resolution VS phi mu+",  16, -3.2, 3.2, 200, -1, 1 );
    hResoVSPhiMinus_ = new TH2F( name+"_ResoVSPhiMinus", "resolution VS phi mu-", 16, -3.2, 3.2, 200, -1, 1 );
    hAbsReso_        = new TH1F( name+"_AbsReso", "resolution", 100, 0, 1 );
    hAbsResoVSPt_    = new TH2F( name+"_AbsResoVSPt", "Abs resolution VS pt", 200, 0, maxPt, 100, 0, 1 );
    hAbsResoVSEta_   = new TH2F( name+"_AbsResoVSEta", "Abs resolution VS eta", 64, -3.2, 3.2, 100, 0, 1 );
    hAbsResoVSPhi_   = new TH2F( name+"_AbsResoVSPhi", "Abs resolution VS phi", 16, -3.2, 3.2, 100, 0, 1 );
    if( doProfiles_ ) {
      hResoVSPt_prof_    = new TProfile(name+"_ResoVSPt_prof", "resolution VS pt", 100, 0, maxPt, yMinPt, yMaxPt );
      hResoVSPt_Bar_prof_    = new TProfile(name+"_ResoVSPt_Bar_prof", "resolution VS pt Barrel", 100, 0, maxPt, yMinPt, yMaxPt );
      hResoVSPt_Endc_17_prof_   = new TProfile(name+"_ResoVSPt_Endc_1.7_prof", "resolution VS pt Endcap (1.4<eta<1.7)", 100, 0, maxPt, yMinPt, yMaxPt );
      hResoVSPt_Endc_20_prof_   = new TProfile(name+"_ResoVSPt_Endc_2.0_prof", "resolution VS pt Endcap (1.7<eta<2.0)", 100, 0, maxPt, yMinPt, yMaxPt );
      hResoVSPt_Endc_24_prof_   = new TProfile(name+"_ResoVSPt_Endc_2.4_prof", "resolution VS pt Endcap (2.0<eta<2.4)", 100, 0, maxPt, yMinPt, yMaxPt );
      hResoVSPt_Ovlap_prof_  = new TProfile(name+"_ResoVSPt_Ovlap_prof", "resolution VS pt Overlap", 100, 0, maxPt, yMinPt, yMaxPt );
      hResoVSEta_prof_       = new TProfile(name+"_ResoVSEta_prof", "resolution VS eta", 200, -3.0, 3.0, yMinEta, yMaxEta );
      hResoVSPhi_prof_       = new TProfile(name+"_ResoVSPhi_prof", "resolution VS phi", 16, -3.2, 3.2, -1, 1 );
    }
  }

HResolutionVSPart::HResolutionVSPart ( const TString &  name, TFile *  file  )

inline

Definition at line 1192 of file Histograms.h.

References dataset::name.

                                                       {
    name_=name;
    hReso_           = (TH1F *) file->Get(name+"_Reso");
    hResoVSPtEta_    = (TH2F *) file->Get(name+"_ResoVSPtEta");
    hResoVSPt_       = (TH2F *) file->Get(name+"_ResoVSPt");
    hResoVSPt_Bar_   = (TH2F *) file->Get(name+"_ResoVSPt");
    hResoVSPt_Endc_17_  = (TH2F *) file->Get(name+"_ResoVSPt");
    hResoVSPt_Endc_20_  = (TH2F *) file->Get(name+"_ResoVSPt");
    hResoVSPt_Endc_24_  = (TH2F *) file->Get(name+"_ResoVSPt");
    hResoVSPt_Ovlap_ = (TH2F *) file->Get(name+"_ResoVSPt");
    hResoVSEta_      = (TH2F *) file->Get(name+"_ResoVSEta");
    hResoVSTheta_    = (TH2F *) file->Get(name+"_ResoVSTheta");
    hResoVSPhiPlus_  = (TH2F *) file->Get(name+"_ResoVSPhiPlus");
    hResoVSPhiMinus_ = (TH2F *) file->Get(name+"_ResoVSPhiMinus");
    hAbsReso_        = (TH1F *) file->Get(name+"_AbsReso");
    hAbsResoVSPt_    = (TH2F *) file->Get(name+"_AbsResoVSPt");
    hAbsResoVSEta_   = (TH2F *) file->Get(name+"_AbsResoVSEta");
    hAbsResoVSPhi_   = (TH2F *) file->Get(name+"_AbsResoVSPhi");
    if( doProfiles_ ) {
      hResoVSPt_prof_  = (TProfile *) file->Get(name+"_ResoVSPt_prof");
      hResoVSPt_Bar_prof_ = (TProfile *) file->Get(name+"_ResoVSPt_prof"); 
      hResoVSPt_Endc_17_prof_ = (TProfile *) file->Get(name+"_ResoVSPt_1.7_prof"); 
      hResoVSPt_Endc_20_prof_ = (TProfile *) file->Get(name+"_ResoVSPt_2.0_prof"); 
      hResoVSPt_Endc_24_prof_ = (TProfile *) file->Get(name+"_ResoVSPt_2.4_prof"); 
      hResoVSPt_Ovlap_prof_= (TProfile *) file->Get(name+"_ResoVSPt_prof"); 
      hResoVSEta_prof_ = (TProfile *) file->Get(name+"_ResoVSEta_prof");
      hResoVSPhi_prof_ = (TProfile *) file->Get(name+"_ResoVSPhi_prof");
    }
  }

HResolutionVSPart::~HResolutionVSPart ( )

inlineoverride

Definition at line 1222 of file Histograms.h.

                                {
    delete hReso_;
    delete hResoVSPtEta_;
    delete hResoVSPt_;
    delete hResoVSPt_Bar_;
    delete hResoVSPt_Endc_17_;
    delete hResoVSPt_Endc_20_;
    delete hResoVSPt_Endc_24_;
    delete hResoVSPt_Ovlap_;
    delete hResoVSEta_;
    delete hResoVSTheta_;
    delete hResoVSPhiMinus_;
    delete hResoVSPhiPlus_;
    delete hAbsReso_;
    delete hAbsResoVSPt_;
    delete hAbsResoVSEta_;
    delete hAbsResoVSPhi_;
    if( doProfiles_ ) {
      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 hResoVSPhi_prof_;
    }
  }

Member Function Documentation

void HResolutionVSPart::Clear ( )

inlineoverride

Definition at line 1327 of file Histograms.h.

                        {
    hReso_->Clear();
    hResoVSPtEta_->Clear();
    hResoVSPt_->Clear();
    hResoVSPt_Bar_->Clear();
    hResoVSPt_Endc_17_->Clear();
    hResoVSPt_Endc_20_->Clear();
    hResoVSPt_Endc_24_->Clear();
    hResoVSPt_Ovlap_->Clear();
    hResoVSEta_->Clear();
    hResoVSTheta_->Clear();
    hResoVSPhiPlus_->Clear();
    hResoVSPhiMinus_->Clear();
    hAbsReso_->Clear();
    hAbsResoVSPt_->Clear();
    hAbsResoVSEta_->Clear();
    hAbsResoVSPhi_->Clear();
    if( doProfiles_ ) {
      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();
      hResoVSPhi_prof_->Clear();
    }
  }

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

inlineoverride

Definition at line 1251 of file Histograms.h.

References PVValHelper::eta, and EnergyCorrector::pt.

                                                                                                      { 
    double pt = p4.Pt();
    double eta = p4.Eta();
    hReso_->Fill(resValue);
    hResoVSPtEta_->Fill(pt, eta,resValue); 
    hResoVSPt_->Fill(pt,resValue); 
    if(fabs(eta)<=0.9)
      hResoVSPt_Bar_->Fill(pt,resValue);
    else if(fabs(eta)>0.9 && fabs(eta)<=1.4)
      hResoVSPt_Ovlap_->Fill(pt,resValue);
    else if(fabs(eta)>1.4 && fabs(eta)<=1.7)
      hResoVSPt_Endc_17_->Fill(pt,resValue);
    else if(fabs(eta)>1.7 && fabs(eta)<=2.0)
      hResoVSPt_Endc_20_->Fill(pt,resValue);
    else
      hResoVSPt_Endc_24_->Fill(pt,resValue);
     
    hResoVSEta_->Fill(eta,resValue);
    hResoVSTheta_->Fill(p4.Theta(),resValue);
    if(charge>0)
      hResoVSPhiPlus_->Fill(p4.Phi(),resValue); 
    else if(charge<0)
      hResoVSPhiMinus_->Fill(p4.Phi(),resValue); 
    hAbsReso_->Fill(fabs(resValue)); 
    hAbsResoVSPt_->Fill(pt,fabs(resValue)); 
    hAbsResoVSEta_->Fill(eta,fabs(resValue)); 
    hAbsResoVSPhi_->Fill(p4.Phi(),fabs(resValue));     
    if( doProfiles_ ) {
      hResoVSPt_prof_->Fill(p4.Pt(),resValue); 
      if(fabs(eta)<=0.9)
    hResoVSPt_Bar_prof_->Fill(p4.Pt(),resValue); 
      else if(fabs(eta)>0.9 && fabs(eta)<=1.4)
    hResoVSPt_Ovlap_prof_->Fill(pt,resValue);
      else if(fabs(eta)>1.4 && fabs(eta)<=1.7)
    hResoVSPt_Endc_17_prof_->Fill(pt,resValue);
      else if(fabs(eta)>1.7 && fabs(eta)<=2.0)
    hResoVSPt_Endc_20_prof_->Fill(pt,resValue);
      else
    hResoVSPt_Endc_24_prof_->Fill(pt,resValue);
      
      hResoVSEta_prof_->Fill(p4.Eta(),resValue); 
      hResoVSPhi_prof_->Fill(p4.Phi(),resValue); 
    }
  }

void HResolutionVSPart::Write ( )

inlineoverride

Definition at line 1296 of file Histograms.h.

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

    hReso_->Write();
    hResoVSPtEta_->Write();
    hResoVSPt_->Write();
    hResoVSPt_Bar_->Write();
    hResoVSPt_Endc_17_->Write();
    hResoVSPt_Endc_20_->Write();
    hResoVSPt_Endc_24_->Write();
    hResoVSPt_Ovlap_->Write();
    hResoVSEta_->Write();
    hResoVSTheta_->Write();
    hResoVSPhiMinus_->Write();
    hResoVSPhiPlus_->Write();
    hAbsReso_->Write();
    hAbsResoVSPt_->Write();
    hAbsResoVSEta_->Write();
    hAbsResoVSPhi_->Write();
    if( doProfiles_ ) {
      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();
      hResoVSPhi_prof_->Write();
    }
  }

Member Data Documentation

bool HResolutionVSPart::doProfiles_

Definition at line 1381 of file Histograms.h.

TH1F* HResolutionVSPart::hAbsReso_

Definition at line 1377 of file Histograms.h.

TH2F* HResolutionVSPart::hAbsResoVSEta_

Definition at line 1379 of file Histograms.h.

TH2F* HResolutionVSPart::hAbsResoVSPhi_

Definition at line 1380 of file Histograms.h.

TH2F* HResolutionVSPart::hAbsResoVSPt_

Definition at line 1378 of file Histograms.h.

TH1F* HResolutionVSPart::hReso_

Definition at line 1357 of file Histograms.h.

TH2F* HResolutionVSPart::hResoVSEta_

Definition at line 1371 of file Histograms.h.

TProfile* HResolutionVSPart::hResoVSEta_prof_

Definition at line 1373 of file Histograms.h.

TProfile* HResolutionVSPart::hResoVSPhi_prof_

Definition at line 1376 of file Histograms.h.

TH2F* HResolutionVSPart::hResoVSPhiMinus_

Definition at line 1374 of file Histograms.h.

TH2F* HResolutionVSPart::hResoVSPhiPlus_

Definition at line 1375 of file Histograms.h.

TH2F* HResolutionVSPart::hResoVSPt_

Definition at line 1359 of file Histograms.h.

TH2F* HResolutionVSPart::hResoVSPt_Bar_

Definition at line 1360 of file Histograms.h.

TProfile* HResolutionVSPart::hResoVSPt_Bar_prof_

Definition at line 1366 of file Histograms.h.

TH2F* HResolutionVSPart::hResoVSPt_Endc_17_

Definition at line 1361 of file Histograms.h.

TProfile* HResolutionVSPart::hResoVSPt_Endc_17_prof_

Definition at line 1367 of file Histograms.h.

TH2F* HResolutionVSPart::hResoVSPt_Endc_20_

Definition at line 1362 of file Histograms.h.

TProfile* HResolutionVSPart::hResoVSPt_Endc_20_prof_

Definition at line 1368 of file Histograms.h.

TH2F* HResolutionVSPart::hResoVSPt_Endc_24_

Definition at line 1363 of file Histograms.h.

TProfile* HResolutionVSPart::hResoVSPt_Endc_24_prof_

Definition at line 1369 of file Histograms.h.

TH2F* HResolutionVSPart::hResoVSPt_Ovlap_

Definition at line 1364 of file Histograms.h.

TProfile* HResolutionVSPart::hResoVSPt_Ovlap_prof_

Definition at line 1370 of file Histograms.h.

TProfile* HResolutionVSPart::hResoVSPt_prof_

Definition at line 1365 of file Histograms.h.

TH2F* HResolutionVSPart::hResoVSPtEta_

Definition at line 1358 of file Histograms.h.

TH2F* HResolutionVSPart::hResoVSTheta_

Definition at line 1372 of file Histograms.h.