A set of histograms for efficiency 1D DT RecHits. More...

#include <Histograms.h>

Public Member Functions
void	ComputeEfficiency ()
void	Fill (float distSimHit, float etaSimHit, float phiSimHit, bool fillRecHit)
	HEff1DHit (std::string name_)
	HEff1DHit (TString name_, TFile *file)
void	Write ()
	~HEff1DHit ()
Public Attributes
TH1F *	hDistMuSimHit
TH1F *	hDistRecHit
TH1F *	hEffVsDist
TH1F *	hEffVsEta
TH1F *	hEffVsPhi
TH1F *	hEtaMuSimHit
TH1F *	hEtaRecHit
TH1F *	hPhiMuSimHit
TH1F *	hPhiRecHit
TString	name

Detailed Description

A set of histograms for efficiency 1D DT RecHits.

Definition at line 139 of file Histograms.h.

Constructor & Destructor Documentation

HEff1DHit::HEff1DHit ( std::string name_ ) [inline]

Definition at line 141 of file Histograms.h.

References hDistMuSimHit, hDistRecHit, hEffVsDist, hEffVsEta, hEffVsPhi, hEtaMuSimHit, hEtaRecHit, hPhiMuSimHit, hPhiRecHit, MultiGaussianStateTransform::N, name, and Pi.

                              {
      TString N = name_.c_str();
      name=N;

      hEtaMuSimHit       = new TH1F("1D_"+N+"_hEtaMuSimHit", "SimHit Eta distribution",
                                    100, -1.5, 1.5);
      hEtaRecHit         = new TH1F("1D_"+N+"_hEtaRecHit", "SimHit Eta distribution with 1D RecHit",
                                    100, -1.5, 1.5);
      hEffVsEta = 0;


      hPhiMuSimHit       = new TH1F("1D_"+N+"_hPhiMuSimHit", "SimHit Phi distribution",
                                    100, -TMath::Pi(),TMath::Pi());
      hPhiRecHit         = new TH1F("1D_"+N+"_hPhiRecHit", "SimHit Phi distribution with 1D RecHit",
                                    100, -TMath::Pi(),TMath::Pi());
      hEffVsPhi = 0;


      hDistMuSimHit       = new TH1F("1D_"+N+"_hDistMuSimHit", "SimHit Distance from wire distribution",
                                     100, 0, 2.5);
      hDistRecHit         = new TH1F("1D_"+N+"_hDistRecHit", "SimHit Distance from wire distribution with 1D RecHit",
                                     100, 0, 2.5);
      hEffVsDist = 0;

    }

HEff1DHit::HEff1DHit	(	TString	name_,
		TFile *	file
	)		`[inline]`

Definition at line 167 of file Histograms.h.

References hDistMuSimHit, hDistRecHit, hEffVsDist, hEffVsEta, hEffVsPhi, hEtaMuSimHit, hEtaRecHit, hPhiMuSimHit, hPhiRecHit, and name.

                                          {
      name=name_;
      hEtaMuSimHit        = (TH1F *) file->Get("1D_"+name+"_hEtaMuSimHit");
      hEtaRecHit          = (TH1F *) file->Get("1D_"+name+"_hEtaRecHit");
      hEffVsEta           = (TH1F *) file->Get("1D_"+name+"_hEffVsEta");

      hPhiMuSimHit        = (TH1F *) file->Get("1D_"+name+"_hPhiMuSimHit");
      hPhiRecHit          = (TH1F *) file->Get("1D_"+name+"_hPhiRecHit");
      hEffVsPhi           = (TH1F *) file->Get("1D_"+name+"_hEffVsPhi");

      hDistMuSimHit       = (TH1F *) file->Get("1D_"+name+"_hDistMuSimHit");
      hDistRecHit         = (TH1F *) file->Get("1D_"+name+"_hDistRecHit");
      hEffVsDist          = (TH1F *) file->Get("1D_"+name+"_hEffVsDist");
    }

HEff1DHit::~HEff1DHit ( ) [inline]

Definition at line 183 of file Histograms.h.

                {

      //     delete hEtaMuSimHit;
      //     delete hEtaRecHit;
      //     if(hEffVsEta != 0)
      //       delete hEffVsEta;

      //     delete hPhiMuSimHit;
      //     delete hPhiRecHit;
      //     if(hEffVsPhi != 0)
      //       delete hEffVsPhi;

      //     delete hDistMuSimHit;
      //     delete hDistRecHit;
      //     if(hEffVsDist != 0)
      //       delete hEffVsDist;

    }

Member Function Documentation

void HEff1DHit::ComputeEfficiency ( ) [inline]

Definition at line 219 of file Histograms.h.

References newFWLiteAna::bin, error, hDistMuSimHit, hDistRecHit, hEffVsDist, hEffVsEta, hEffVsPhi, hEtaMuSimHit, hEtaRecHit, hPhiMuSimHit, hPhiRecHit, name, and mathSSE::sqrt().

                             {

      hEffVsEta = (TH1F *) hEtaRecHit->Clone();
      hEffVsEta->SetName("1D_"+name+"_hEffVsEta");
      hEffVsEta->SetTitle("1D RecHit Efficiency as a function of Eta");
      hEffVsEta->Divide(hEtaMuSimHit);
      // Set the error accordingly to binomial statistics
      int nBinsEta = hEffVsEta->GetNbinsX();
      for(int bin = 1; bin <=  nBinsEta; bin++) {
        float nSimHit = hEtaMuSimHit->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->Clone();
      hEffVsPhi->SetName("1D_"+name+"_hEffVsPhi");
      hEffVsPhi->SetTitle("1D RecHit Efficiency as a function of Phi");
      hEffVsPhi->Divide(hPhiMuSimHit);
      // Set the error accordingly to binomial statistics
      int nBinsPhi = hEffVsPhi->GetNbinsX();
      for(int bin = 1; bin <=  nBinsPhi; bin++) {
        float nSimHit = hPhiMuSimHit->GetBinContent(bin);
        float eff = hEffVsPhi->GetBinContent(bin);
        float error = 0;
        if(nSimHit != 0) {
          error = sqrt(eff*(1-eff)/nSimHit);
        }
        hEffVsPhi->SetBinError(bin, error);
      }

      hEffVsDist = (TH1F *) hDistRecHit->Clone();
      hEffVsDist->SetName("1D_"+name+"_hEffVsDist");
      hEffVsDist->SetTitle("1D RecHit Efficiency as a function of Dist");
      hEffVsDist->Divide(hDistMuSimHit);
      // Set the error accordingly to binomial statistics
      int nBinsDist = hEffVsDist->GetNbinsX();
      for(int bin = 1; bin <=  nBinsDist; bin++) {
        float nSimHit = hDistMuSimHit->GetBinContent(bin);
        float eff = hEffVsDist->GetBinContent(bin);
        float error = 0;
        if(nSimHit != 0) {
          error = sqrt(eff*(1-eff)/nSimHit);
        }
        hEffVsDist->SetBinError(bin, error);
      }


    }

void HEff1DHit::Fill	(	float	distSimHit,
		float	etaSimHit,
		float	phiSimHit,
		bool	fillRecHit
	)		`[inline]`

Definition at line 202 of file Histograms.h.

References hDistMuSimHit, hDistRecHit, hEtaMuSimHit, hEtaRecHit, hPhiMuSimHit, and hPhiRecHit.

Referenced by DTRecHitQuality::compute().

                               {

      hEtaMuSimHit->Fill(etaSimHit);
      hPhiMuSimHit->Fill(phiSimHit);
      hDistMuSimHit->Fill(distSimHit);
      if(fillRecHit) {
        hEtaRecHit->Fill(etaSimHit);
        hPhiRecHit->Fill(phiSimHit);
        hDistRecHit->Fill(distSimHit);
      }
    }

void HEff1DHit::Write ( ) [inline]

Definition at line 272 of file Histograms.h.

References hDistMuSimHit, hDistRecHit, hEffVsDist, hEffVsEta, hEffVsPhi, hEtaMuSimHit, hEtaRecHit, hPhiMuSimHit, and hPhiRecHit.

                 {
      hEtaMuSimHit->Write();                         
      hEtaRecHit->Write();
      if(hEffVsEta != 0)
        hEffVsEta->Write();                    
      hPhiMuSimHit->Write();                         
      hPhiRecHit->Write();                         
      if(hEffVsPhi != 0)
        hEffVsPhi->Write();               
      hDistMuSimHit->Write();                         
      hDistRecHit->Write();                         
      if(hEffVsDist != 0)
        hEffVsDist->Write();
    }