HEff4DHit Class Reference

A set of histograms for efficiency 4D RecHits. More...

#include <Histograms.h>

Public Member Functions
void	ComputeEfficiency ()
void	Fill (float etaSimSegm, float phiSimSegm, float xSimSegm, float ySimSegm, float alphaSimSegm, float betaSimSegm, bool fillRecHit, int nSeg)
	HEff4DHit (std::string name_, DQMStore *dbe_)
	~HEff4DHit ()

Public Attributes
MonitorElement *	hAlphaRecHit
MonitorElement *	hAlphaSimSegm
MonitorElement *	hBetaRecHit
MonitorElement *	hBetaSimSegm
TH1F *	hEffVsAlpha
TH1F *	hEffVsBeta
TH1F *	hEffVsEta
TH1F *	hEffVsPhi
TH1F *	hEffVsX
TH1F *	hEffVsY
MonitorElement *	hEtaRecHit
MonitorElement *	hEtaSimSegm
MonitorElement *	hNSeg
MonitorElement *	hPhiRecHit
MonitorElement *	hPhiSimSegm
MonitorElement *	hXRecHit
MonitorElement *	hXSimSegm
MonitorElement *	hYRecHit
MonitorElement *	hYSimSegm
TString	name

Detailed Description

A set of histograms for efficiency 4D RecHits.

Definition at line 1006 of file Histograms.h.

Constructor & Destructor Documentation

HEff4DHit::HEff4DHit ( std::string  name_, DQMStore *  dbe_  )

inline

Definition at line 1008 of file Histograms.h.

References DQMStore::book1D(), HRes1DHit::name, Pi, DQMStore::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                             {
      std::string pre ="4D_";
      pre += name_;
      name = pre;
      dbe_->setCurrentFolder("DT/4DSegments/");     
      hEtaSimSegm=0;hEtaSimSegm     = dbe_->book1D(pre+"_hEtaSimSegm", "Eta of SimHit segment", 100, -1.5, 1.5);
      hEtaRecHit=0;hEtaRecHit      = dbe_->book1D(pre+"_hEtaRecHit", "Eta distribution of SimHit segment with 4D RecHit",
                                 100, -1.5, 1.5);
      hEffVsEta       = 0;

      hPhiSimSegm=0;hPhiSimSegm     = dbe_->book1D(pre+"_hPhiSimSegm", "Phi of SimHit segment",
                                 100, -TMath::Pi(),TMath::Pi());
      hPhiRecHit=0;hPhiRecHit      = dbe_->book1D(pre+"_hPhiRecHit", "Phi distribution of SimHit segment with 4D RecHit",
                                 100, -TMath::Pi(),TMath::Pi());
      hEffVsPhi       = 0;


      hXSimSegm=0;hXSimSegm       = dbe_->book1D(pre+"_hXSimSegm", "X position in Chamber of SimHit segment (cm)",
                                 100, -200, 200);
      hXRecHit=0;hXRecHit        = dbe_->book1D(pre+"_hXRecHit", "X position in Chamber of SimHit segment with 4D RecHit (cm)",
                                 100, -200, 200);
      hEffVsX         = 0;

      hYSimSegm=0;hYSimSegm       = dbe_->book1D(pre+"_hYSimSegm", "Y position in Chamber of SimHit segment (cm)",
                                 100, -200, 200);
      hYRecHit=0;hYRecHit        = dbe_->book1D(pre+"_hYRecHit", "Y position in Chamber of SimHit segment with 4D RecHit (cm)",
                                 100, -200, 200);
      hEffVsY         = 0;

      hAlphaSimSegm=0;hAlphaSimSegm   = dbe_->book1D(pre+"_hAlphaSimSegm", "Alpha of SimHit segment (rad)",
                                 100, -1.5, 1.5);
      hAlphaRecHit=0;hAlphaRecHit    = dbe_->book1D(pre+"_hAlphaRecHit", "Alpha of SimHit segment with 4D RecHit (rad)",
                                 100, -1.5, 1.5);
      hEffVsAlpha     = 0;

      hBetaSimSegm=0;hBetaSimSegm   = dbe_->book1D(pre+"_hBetaSimSegm", "Beta of SimHit segment (rad)",
                                100, -2, 2);
      hBetaRecHit=0;hBetaRecHit    = dbe_->book1D(pre+"_hBetaRecHit", "Beta of SimHit segment with 4D RecHit (rad)",
                                100, -2, 2);
      hEffVsBeta     = 0;

      hNSeg =0; hNSeg            = dbe_->book1D(pre+"_hNSeg", "Number of rec segment per sim seg",
                                 20, 0, 20);

    }

HEff4DHit::~HEff4DHit ( )

inline

Definition at line 1054 of file Histograms.h.

                {

    }

Member Function Documentation

void HEff4DHit::ComputeEfficiency ( )

inline

Definition at line 1087 of file Histograms.h.

References stringResolutionProvider_cfi::bin, relativeConstraints::error, HRes1DHit::name, and mathSSE::sqrt().

                             {
      //hEffVsEta = (TH1F *) hEtaRecHit->Clone();
      hEffVsEta = (TH1F *) hEtaRecHit->getTH1();
      TH1F * hEffEtaSimSegm = (TH1F *) hEtaSimSegm->getTH1();
      hEffVsEta->SetName(name+"_hEffVsEta");
      hEffVsEta->SetTitle("4D RecHit Efficiency as a function of Eta");
      hEffVsEta->Divide(hEffEtaSimSegm);
      // Set the error accordingly to binomial statistics
      int nBinsEta = hEffVsEta->GetNbinsX();
      for(int bin = 1; bin <=  nBinsEta; bin++) {
        float nSimHit = hEffEtaSimSegm->GetBinContent(bin);
        float eff = hEffVsEta->GetBinContent(bin);
        float error = 0;
        if(nSimHit != 0) {
          error = sqrt(eff*(1-eff)/nSimHit);
        }
        hEffVsEta->SetBinError(bin, error);
      }

      hEffVsPhi = (TH1F *) hPhiRecHit->getTH1();
      TH1F * hEffPhiSimSegm = (TH1F *) hPhiSimSegm->getTH1();
      hEffVsPhi->SetName(name+"_hEffVsPhi");
      hEffVsPhi->SetTitle("4D RecHit Efficiency as a function of Phi");
      hEffVsPhi->Divide(hEffPhiSimSegm);
      // Set the error accordingly to binomial statistics
      int nBinsPhi = hEffVsPhi->GetNbinsX();
      for(int bin = 1; bin <=  nBinsPhi; bin++) {
        float nSimHit = hEffPhiSimSegm->GetBinContent(bin);
        float eff = hEffVsPhi->GetBinContent(bin);
        float error = 0;
        if(nSimHit != 0) {
          error = sqrt(eff*(1-eff)/nSimHit);
        }
        hEffVsPhi->SetBinError(bin, error);
      }

      hEffVsX = (TH1F *) hXRecHit->getTH1();
      TH1F * hEffXSimSegm = (TH1F *) hXSimSegm->getTH1();
      hEffVsX->SetName(name+"_hEffVsX");
      hEffVsX->SetTitle("4D RecHit Efficiency as a function of x position in Chamber");
      hEffVsX->Divide(hEffXSimSegm);
      // Set the error accordingly to binomial statistics
      int nBinsX = hEffVsX->GetNbinsX();
      for(int bin = 1; bin <=  nBinsX; bin++) {
        float nSimHit = hEffXSimSegm->GetBinContent(bin);
        float eff = hEffVsX->GetBinContent(bin);
        float error = 0;
        if(nSimHit != 0) {
          error = sqrt(eff*(1-eff)/nSimHit);
        }
        hEffVsX->SetBinError(bin, error);
      }


      hEffVsY = (TH1F *) hYRecHit->getTH1();
      TH1F * hEffYSimSegm = (TH1F *) hYSimSegm->getTH1();
      hEffVsY->SetName(name+"_hEffVsY");
      hEffVsY->SetTitle("4D RecHit Efficiency as a function of y position in Chamber");
      hEffVsY->Divide(hEffYSimSegm);
      // Set the error accordingly to binomial statistics
      int nBinsY = hEffVsY->GetNbinsX();
      for(int bin = 1; bin <=  nBinsY; bin++) {
        float nSimHit = hEffYSimSegm->GetBinContent(bin);
        float eff = hEffVsY->GetBinContent(bin);
        float error = 0;
        if(nSimHit != 0) {
          error = sqrt(eff*(1-eff)/nSimHit);
        }
        hEffVsY->SetBinError(bin, error);
      }

      hEffVsAlpha = (TH1F *) hAlphaRecHit->getTH1();
      TH1F * hEffAlphaSimSegm = (TH1F *) hAlphaSimSegm->getTH1();
      hEffVsAlpha->SetName(name+"_hEffVsAlpha");
      hEffVsAlpha->SetTitle("4D RecHit Efficiency as a function of alpha");
      hEffVsAlpha->Divide(hEffAlphaSimSegm);
      // Set the error accordingly to binomial statistics
      int nBinsAlpha = hEffVsAlpha->GetNbinsX();
      for(int bin = 1; bin <=  nBinsAlpha; bin++) {
        float nSimHit = hEffAlphaSimSegm->GetBinContent(bin);
        float eff = hEffVsAlpha->GetBinContent(bin);
        float error = 0;
        if(nSimHit != 0) {
          error = sqrt(eff*(1-eff)/nSimHit);
        }
        hEffVsAlpha->SetBinError(bin, error);
      }


      hEffVsBeta = (TH1F *) hBetaRecHit->getTH1();
      TH1F * hEffBetaSimSegm = (TH1F *) hBetaSimSegm->getTH1();
      hEffVsBeta->SetName(name+"_hEffVsBeta");
      hEffVsBeta->SetTitle("4D RecHit Efficiency as a function of beta");
      hEffVsBeta->Divide(hEffBetaSimSegm);
      // Set the error accordingly to binomial statistics
      int nBinsBeta = hEffVsBeta->GetNbinsX();
      for(int bin = 1; bin <=  nBinsBeta; bin++) {
        float nSimHit = hEffBetaSimSegm->GetBinContent(bin);
        float eff = hEffVsBeta->GetBinContent(bin);
        float error = 0;
        if(nSimHit != 0) {
          error = sqrt(eff*(1-eff)/nSimHit);
        }
        hEffVsBeta->SetBinError(bin, error);
      }
    }

void HEff4DHit::Fill ( float  etaSimSegm, float  phiSimSegm, float  xSimSegm, float  ySimSegm, float  alphaSimSegm, float  betaSimSegm, bool  fillRecHit, int  nSeg  )

inline

Definition at line 1058 of file Histograms.h.

Referenced by DTSegment4DQuality::analyze().

                        {

      hEtaSimSegm->Fill(etaSimSegm);
      hPhiSimSegm->Fill(phiSimSegm);
      hXSimSegm->Fill(xSimSegm);
      hYSimSegm->Fill(ySimSegm);
      hAlphaSimSegm->Fill(alphaSimSegm);
      hBetaSimSegm->Fill(betaSimSegm);
      hNSeg->Fill(nSeg);

      if(fillRecHit) {
        hEtaRecHit->Fill(etaSimSegm);  
        hPhiRecHit->Fill(phiSimSegm);
        hXRecHit->Fill(xSimSegm);
        hYRecHit->Fill(ySimSegm);
        hAlphaRecHit->Fill(alphaSimSegm);
        hBetaRecHit->Fill(betaSimSegm);
      }
    }

Member Data Documentation

MonitorElement* HEff4DHit::hAlphaRecHit

Definition at line 1208 of file Histograms.h.

MonitorElement* HEff4DHit::hAlphaSimSegm

Definition at line 1207 of file Histograms.h.

MonitorElement* HEff4DHit::hBetaRecHit

Definition at line 1211 of file Histograms.h.

MonitorElement* HEff4DHit::hBetaSimSegm

Definition at line 1210 of file Histograms.h.

TH1F* HEff4DHit::hEffVsAlpha

Definition at line 1209 of file Histograms.h.

TH1F* HEff4DHit::hEffVsBeta

Definition at line 1212 of file Histograms.h.

TH1F* HEff4DHit::hEffVsEta

Definition at line 1197 of file Histograms.h.

TH1F* HEff4DHit::hEffVsPhi

Definition at line 1200 of file Histograms.h.

TH1F* HEff4DHit::hEffVsX

Definition at line 1203 of file Histograms.h.

TH1F* HEff4DHit::hEffVsY

Definition at line 1206 of file Histograms.h.

MonitorElement* HEff4DHit::hEtaRecHit

Definition at line 1196 of file Histograms.h.

MonitorElement* HEff4DHit::hEtaSimSegm

Definition at line 1195 of file Histograms.h.

MonitorElement* HEff4DHit::hNSeg

Definition at line 1213 of file Histograms.h.

MonitorElement* HEff4DHit::hPhiRecHit

Definition at line 1199 of file Histograms.h.

MonitorElement* HEff4DHit::hPhiSimSegm

Definition at line 1198 of file Histograms.h.

MonitorElement* HEff4DHit::hXRecHit

Definition at line 1202 of file Histograms.h.

MonitorElement* HEff4DHit::hXSimSegm

Definition at line 1201 of file Histograms.h.

MonitorElement* HEff4DHit::hYRecHit

Definition at line 1205 of file Histograms.h.

MonitorElement* HEff4DHit::hYSimSegm

Definition at line 1204 of file Histograms.h.

TString HEff4DHit::name

Definition at line 1215 of file Histograms.h.

Referenced by config.CFG::__str__(), validation.Sample::digest(), VIDSelectorBase.VIDSelectorBase::initialize(), and Vispa.Views.PropertyView.Property::valueChanged().