HResolution Class Reference

#include <Histograms.h>

Inheritance diagram for HResolution:

Public Types
typedef dqm::legacy::DQMStore	DQMStore
typedef dqm::legacy::MonitorElement	MonitorElement

Public Member Functions
void	Fill (double p, double pt, double eta, double phi, double rp, double rpt)
void	Fill (double p, double pt, double eta, double phi, double rp, double rpt, double reta, double rphi, double rcharge)
void	Fill (double rp, double rpt, double reta, double rphi, double rcharge)
Int_t	Fill (Double_t x, Double_t y) override
	HResolution (const TString &name, const TString &title, const int totBins, const double &xMin, const double &xMax, const double &yMin, const double &yMax, TDirectory *dir=nullptr)
	HResolution (DQMStore::IBooker &ibooker, std::string dirName_, std::string name, std::string whereIs)
	HResolution (DQMStore::IBooker &ibooker, std::string name, TFile *file)
Int_t	Write (const char *name=nullptr, Int_t option=0, Int_t bufsize=0) override
	~HResolution ()
	~HResolution () override

Protected Attributes
TDirectory *	diffDir_
TProfile *	diffHisto_
TDirectory *	dir2D_
TDirectory *	dir_
TH2F *	histo2D_
TH1F *	resoHisto_

Private Attributes
MonitorElement *	h2Eta
MonitorElement *	h2EtaVsPhi
MonitorElement *	h2EtaVsPt
MonitorElement *	h2P
MonitorElement *	h2Phi
MonitorElement *	h2PhiVsEta
MonitorElement *	h2PhiVsPt
MonitorElement *	h2Pt
MonitorElement *	h2PtVsEta
MonitorElement *	h2PtVsPhi
MonitorElement *	hCharge
MonitorElement *	hEta
MonitorElement *	hP
MonitorElement *	hPhi
MonitorElement *	hPt
std::string	theName
std::string	where

Detailed Description

This histogram class can be used to evaluate the resolution of a variable. It has a TProfile, a TH2F and a TH1F. The TProfile is used to compute the rms of the distribution which is filled in the TH1F (the resolution histogram) in the Write method. If a TDirectory is passed to the constructor, the different histograms are placed in subdirectories.

Definition at line 1800 of file Histograms.h.

Member Typedef Documentation

DQMStore

typedef dqm::legacy::DQMStore HResolution::DQMStore

Definition at line 151 of file Histograms.h.

MonitorElement

typedef dqm::legacy::MonitorElement HResolution::MonitorElement

Definition at line 152 of file Histograms.h.

Constructor & Destructor Documentation

HResolution() [1/3]

HResolution::HResolution ( const TString &  name, const TString &  title, const int  totBins, const double &  xMin, const double &  xMax, const double &  yMin, const double &  yMax, TDirectory *  dir = nullptr  )

inline

Definition at line 1802 of file Histograms.h.

      : dir_(dir), dir2D_(nullptr), diffDir_(nullptr) {
    if (dir_ != nullptr) {
      dir2D_ = (TDirectory*)dir_->Get("2D");
      if (dir2D_ == nullptr)
        dir2D_ = dir_->mkdir("2D");
      diffDir_ = (TDirectory*)dir_->Get("deltaXoverX");
      if (diffDir_ == nullptr)
        diffDir_ = dir->mkdir("deltaXoverX");
    }
    diffHisto_ = new TProfile(name + "_prof", title + " profile", totBins, xMin, xMax, yMin, yMax);
    histo2D_ = new TH2F(name + "2D", title, totBins, xMin, xMax, 4000, yMin, yMax);
    resoHisto_ = new TH1F(name, title, totBins, xMin, xMax);
  }

References diffDir_, diffHisto_, DeadROC_duringRun::dir, dir2D_, dir_, histo2D_, Skims_PA_cff::name, resoHisto_, overlapproblemtsosanalyzer_cfi::title, multiplicitycorr_cfi::xMax, photonAnalyzer_cfi::xMin, multiplicitycorr_cfi::yMax, and photonAnalyzer_cfi::yMin.

~HResolution() [1/2]

HResolution::~HResolution ( )

inlineoverride

Definition at line 1823 of file Histograms.h.

                          {
    delete diffHisto_;
    delete histo2D_;
    delete resoHisto_;
  }

References diffHisto_, histo2D_, and resoHisto_.

HResolution() [2/3]

HResolution::HResolution ( DQMStore::IBooker &  ibooker, std::string  dirName_, std::string  name, std::string  whereIs  )

inline

Definition at line 154 of file Histograms.h.

      : theName(name), where(whereIs) {
    ibooker.cd();
    const std::string &dirName = dirName_;
    //dirName+="/";
    //dirName+=name.c_str();
 
    ibooker.setCurrentFolder(dirName);
 
    double eta = 15.;
    int neta = 800;
    double phi = 12.;
    int nphi = 400;
    double pt = 60.;
    int npt = 2000;
 
    hEta = ibooker.book1D(theName + "_Eta_" + where, "#eta " + theName, neta, -eta, eta);  // 400
    hPhi = ibooker.book1D(theName + "_Phi_" + where, "#phi " + theName, nphi, -phi, phi);  // 100
 
    hP = ibooker.book1D(theName + "_P_" + where, "P " + theName, 400, -4, 4);          // 200
    hPt = ibooker.book1D(theName + "_Pt_" + where, "P_{t} " + theName, npt, -pt, pt);  // 200
 
    hCharge = ibooker.book1D(theName + "_charge_" + where, "Charge " + theName, 4, -2., 2.);
 
    h2Eta = ibooker.book2D(
        theName + "_Eta_vs_Eta" + where, "#eta " + theName + " as a function of #eta", 200, -2.5, 2.5, neta, -eta, eta);
    h2Phi = ibooker.book2D(theName + "_Phi_vs_Phi" + where,
                           "#phi " + theName + " as a function of #phi",
                           100,
                           -Geom::pi(),
                           Geom::pi(),
                           nphi,
                           -phi,
                           phi);
 
    h2P =
        ibooker.book2D(theName + "_P_vs_P" + where, "P " + theName + " as a function of P", 1000, 0, 2000, 400, -4, 4);
    h2Pt = ibooker.book2D(
        theName + "_Pt_vs_Pt" + where, "P_{t} " + theName + " as a function of P_{t}", 1000, 0, 2000, npt, -pt, pt);
 
    h2PtVsEta = ibooker.book2D(
        theName + "_Pt_vs_Eta" + where, "P_{t} " + theName + " as a function of #eta", 200, -2.5, 2.5, npt, -pt, pt);
    h2PtVsPhi = ibooker.book2D(theName + "_Pt_vs_Phi" + where,
                               "P_{t} " + theName + " as a function of #phi",
                               100,
                               -Geom::pi(),
                               Geom::pi(),
                               npt,
                               -pt,
                               pt);
 
    h2EtaVsPt = ibooker.book2D(
        theName + "_Eta_vs_Pt" + where, "#eta " + theName + " as a function of P_{t}", 1000, 0, 2000, neta, -eta, eta);
    h2EtaVsPhi = ibooker.book2D(theName + "_Eta_vs_Phi" + where,
                                "#eta " + theName + " as a function of #phi",
                                100,
                                -Geom::pi(),
                                Geom::pi(),
                                neta,
                                -eta,
                                eta);
 
    h2PhiVsPt = ibooker.book2D(
        theName + "_Phi_vs_Pt" + where, "#phi " + theName + " as a function of P_{t}", 1000, 0, 2000, nphi, -phi, phi);
    h2PhiVsEta = ibooker.book2D(
        theName + "_Phi_vs_Eta" + where, "#phi " + theName + " as a function of #eta", 200, -2.5, 2.5, nphi, -phi, phi);
  }

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::NavigatorBase::cd(), TrackerOfflineValidation_Dqm_cff::dirName, PVValHelper::eta, h2Eta, h2EtaVsPhi, h2EtaVsPt, h2P, h2Phi, h2PhiVsEta, h2PhiVsPt, h2Pt, h2PtVsEta, h2PtVsPhi, hCharge, hEta, hP, hPhi, hPt, neta, nphi, phi, Geom::pi(), DiDispStaMuonMonitor_cfi::pt, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, theName, and where.

HResolution() [3/3]

HResolution::HResolution ( DQMStore::IBooker &  ibooker, std::string  name, TFile *  file  )

inline

Definition at line 222 of file Histograms.h.

                                                                     : theName(name) {
    //    dynamic_cast<TH1F*>( file->Get(theName+"") );
  }

~HResolution() [2/2]

HResolution::~HResolution ( )

inline

Definition at line 226 of file Histograms.h.

{}

Member Function Documentation

Fill() [1/4]

void HResolution::Fill ( double  p, double  pt, double  eta, double  phi, double  rp, double  rpt  )

inline

Definition at line 248 of file Histograms.h.

                                                                                {
    hP->Fill(rp);
    hPt->Fill(rpt);
 
    h2P->Fill(p, rp);
    // h2PVsEta->Fill(eta,rp);
    // h2PVsPhi->Fill(phi,rp);
 
    h2Pt->Fill(pt, rpt);
    h2PtVsEta->Fill(eta, rpt);
    h2PtVsPhi->Fill(phi, rpt);
  }

References PVValHelper::eta, dqm::impl::MonitorElement::Fill(), h2P, h2Pt, h2PtVsEta, h2PtVsPhi, hP, hPt, AlCaHLTBitMon_ParallelJobs::p, phi, and DiDispStaMuonMonitor_cfi::pt.

Fill() [2/4]

void HResolution::Fill ( double  p, double  pt, double  eta, double  phi, double  rp, double  rpt, double  reta, double  rphi, double  rcharge  )

inline

Definition at line 228 of file Histograms.h.

                                                                                                                    {
    Fill(rp, rpt, reta, rphi, rcharge);
 
    h2Eta->Fill(eta, reta);
    h2Phi->Fill(phi, rphi);
 
    h2P->Fill(p, rp);
    h2Pt->Fill(pt, rpt);
 
    h2PtVsEta->Fill(eta, rpt);
    h2PtVsPhi->Fill(phi, rpt);
 
    h2EtaVsPt->Fill(pt, reta);
    h2EtaVsPhi->Fill(phi, reta);
 
    h2PhiVsPt->Fill(pt, rphi);
    h2PhiVsEta->Fill(eta, rphi);
  }

References PVValHelper::eta, dqm::impl::MonitorElement::Fill(), Fill(), h2Eta, h2EtaVsPhi, h2EtaVsPt, h2P, h2Phi, h2PhiVsEta, h2PhiVsPt, h2Pt, h2PtVsEta, h2PtVsPhi, AlCaHLTBitMon_ParallelJobs::p, phi, and DiDispStaMuonMonitor_cfi::pt.

Fill() [3/4]

void HResolution::Fill ( double  rp, double  rpt, double  reta, double  rphi, double  rcharge  )

inline

Definition at line 261 of file Histograms.h.

                                                                             {
    hEta->Fill(reta);
    hPhi->Fill(rphi);
 
    hP->Fill(rp);
    hPt->Fill(rpt);
 
    hCharge->Fill(rcharge);
  }

References dqm::impl::MonitorElement::Fill(), hCharge, hEta, hP, hPhi, and hPt.

Fill() [4/4]

Int_t HResolution::Fill ( Double_t  x, Double_t  y  )

inlineoverride

Definition at line 1828 of file Histograms.h.

                                              {
    diffHisto_->Fill(x, y);
    histo2D_->Fill(x, y);
    return 0;
  }

References diffHisto_, histo2D_, x, and y.

Referenced by HTrack::computePull(), HTrack::computeResolution(), HTrack::computeTDRResolution(), and Fill().

Write()

Int_t HResolution::Write ( const char *  name = nullptr, Int_t  option = 0, Int_t  bufsize = 0  )

inlineoverride

Definition at line 1833 of file Histograms.h.

                                                                                        {
    // Loop on all the bins and take the rms.
    // The TProfile bin error is by default the standard error on the mean, that is
    // rms/sqrt(N). If it is created with the "S" option (as we did NOT do), it would
    // already be the rms. Thus we take the error and multiply it by the sqrt of the
    // bin entries to get the rms.
    // bin 0 is the underflow, bin totBins+1 is the overflow.
    unsigned int totBins = diffHisto_->GetNbinsX();
    // std::cout << "totBins = " << totBins << std::endl;
    for (unsigned int iBin = 1; iBin <= totBins; ++iBin) {
      // std::cout << "iBin = " << iBin << ", " << diffHisto_->GetBinError(iBin)*sqrt(diffHisto_->GetBinEntries(iBin)) << std::endl;
      resoHisto_->SetBinContent(iBin, diffHisto_->GetBinError(iBin) * sqrt(diffHisto_->GetBinEntries(iBin)));
    }
    if (dir_ != nullptr)
      dir_->cd();
    resoHisto_->Write();
    if (diffDir_ != nullptr)
      diffDir_->cd();
    diffHisto_->Write();
    if (dir2D_ != nullptr)
      dir2D_->cd();
    histo2D_->Write();
 
    return 0;
  }

References diffDir_, diffHisto_, dir2D_, dir_, histo2D_, resoHisto_, and mathSSE::sqrt().

Member Data Documentation

diffDir_

TDirectory* HResolution::diffDir_

protected

Definition at line 1862 of file Histograms.h.

Referenced by HResolution(), and Write().

diffHisto_

TProfile* HResolution::diffHisto_

protected

Definition at line 1863 of file Histograms.h.

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

dir2D_

TDirectory* HResolution::dir2D_

protected

Definition at line 1861 of file Histograms.h.

Referenced by HResolution(), and Write().

dir_

TDirectory* HResolution::dir_

protected

Definition at line 1860 of file Histograms.h.

Referenced by HResolution(), and Write().

h2Eta

MonitorElement* HResolution::h2Eta

private

Definition at line 283 of file Histograms.h.

Referenced by Fill(), and HResolution().

h2EtaVsPhi

MonitorElement* HResolution::h2EtaVsPhi

private

Definition at line 293 of file Histograms.h.

Referenced by Fill(), and HResolution().

h2EtaVsPt

MonitorElement* HResolution::h2EtaVsPt

private

Definition at line 292 of file Histograms.h.

Referenced by Fill(), and HResolution().

h2P

MonitorElement* HResolution::h2P

private

Definition at line 286 of file Histograms.h.

Referenced by Fill(), and HResolution().

h2Phi

MonitorElement* HResolution::h2Phi

private

Definition at line 284 of file Histograms.h.

Referenced by Fill(), and HResolution().

h2PhiVsEta

MonitorElement* HResolution::h2PhiVsEta

private

Definition at line 296 of file Histograms.h.

Referenced by Fill(), and HResolution().

h2PhiVsPt

MonitorElement* HResolution::h2PhiVsPt

private

Definition at line 295 of file Histograms.h.

Referenced by Fill(), and HResolution().

h2Pt

MonitorElement* HResolution::h2Pt

private

Definition at line 287 of file Histograms.h.

Referenced by Fill(), and HResolution().

h2PtVsEta

MonitorElement* HResolution::h2PtVsEta

private

Definition at line 289 of file Histograms.h.

Referenced by Fill(), and HResolution().

h2PtVsPhi

MonitorElement* HResolution::h2PtVsPhi

private

Definition at line 290 of file Histograms.h.

Referenced by Fill(), and HResolution().

hCharge

MonitorElement* HResolution::hCharge

private

Definition at line 281 of file Histograms.h.

Referenced by Fill(), and HResolution().

hEta

MonitorElement* HResolution::hEta

private

Definition at line 275 of file Histograms.h.

Referenced by Fill(), and HResolution().

histo2D_

TH2F* HResolution::histo2D_

protected

Definition at line 1864 of file Histograms.h.

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

hP

MonitorElement* HResolution::hP

private

Definition at line 278 of file Histograms.h.

Referenced by Fill(), and HResolution().

hPhi

MonitorElement* HResolution::hPhi

private

Definition at line 276 of file Histograms.h.

Referenced by Fill(), and HResolution().

hPt

MonitorElement* HResolution::hPt

private

Definition at line 279 of file Histograms.h.

Referenced by Fill(), and HResolution().

resoHisto_

TH1F* HResolution::resoHisto_

protected

Definition at line 1865 of file Histograms.h.

Referenced by HResolution(), Write(), and ~HResolution().

theName

std::string HResolution::theName

private

Definition at line 272 of file Histograms.h.

Referenced by HResolution().

where

std::string HResolution::where

private

Definition at line 273 of file Histograms.h.

Referenced by HResolution().