#include <Histograms.h>

Public Member Functions
void	ComputeEfficiency ()

void	Fill (float etaSimSegm, float phiSimSegm, float posSimSegm, float angleSimSegm, bool fillRecHit)

	HEff2DHit (std::string name_)

	HEff2DHit (TString name_, TFile *file)

void	Write ()

	~HEff2DHit ()

Public Attributes
TH1F *	hAngleRecHit

TH1F *	hAngleSimSegm

TH1F *	hEffVsAngle

TH1F *	hEffVsEta

TH1F *	hEffVsPhi

TH1F *	hEffVsPos

TH1F *	hEtaRecHit

TH1F *	hEtaSimSegm

TH1F *	hPhiRecHit

TH1F *	hPhiSimSegm

TH1F *	hPosRecHit

TH1F *	hPosSimSegm

TString	name

Detailed Description

Definition at line 427 of file Histograms.h.

Constructor & Destructor Documentation

HEff2DHit::HEff2DHit ( std::string name_ )

inline

Definition at line 429 of file Histograms.h.

References hAngleRecHit, hAngleSimSegm, hEffVsAngle, hEffVsEta, hEffVsPhi, hEffVsPos, hEtaRecHit, hEtaSimSegm, hPhiRecHit, hPhiSimSegm, hPosRecHit, hPosSimSegm, MultiGaussianStateTransform::N, name, and Pi.

                               {
       TString N = name_.c_str();
       name=N;
 
       hEtaSimSegm     = new TH1F("2D_"+N+"_hEtaSimSegm", "Eta of SimHit segment", 100, -1.5, 1.5);
       hEtaRecHit      = new TH1F("2D_"+N+"_hEtaRecHit", "Eta distribution of SimHit segment with 2D RecHit",
                                  100, -1.5, 1.5);
       hEffVsEta       = 0;
 
       hPhiSimSegm     = new TH1F("2D_"+N+"_hPhiSimSegm", "Phi of SimHit segment",
                                  100, -TMath::Pi(),TMath::Pi());
       hPhiRecHit      = new TH1F("2D_"+N+"_hPhiRecHit", "Phi distribution of SimHit segment with 2D RecHit",
                                  100, -TMath::Pi(),TMath::Pi());
       hEffVsPhi       = 0;
 
 
       hPosSimSegm     = new TH1F("2D_"+N+"_hPosSimSegm", "Position in SL of SimHit segment (cm)",
                                  100, -250, 250);
       hPosRecHit      = new TH1F("2D_"+N+"_hPosRecHit", "Position in SL of SimHit segment with 2D RecHit (cm)",
                                  100, -250, 250);
       hEffVsPos       = 0;
 
 
       hAngleSimSegm   = new TH1F("2D_"+N+"_hAngleSimSegm", "Angle of SimHit segment (rad)",
                                  100, -2, 2);
       hAngleRecHit    = new TH1F("2D_"+N+"_hAngleRecHit", "Angle of SimHit segment with 2D RecHit (rad)",
                                  100, -2, 2);
       hEffVsAngle     = 0;
 
     }

HEff2DHit::HEff2DHit	(	TString	name_,
		TFile *	file
	)

inline

Definition at line 460 of file Histograms.h.

References hAngleRecHit, hAngleSimSegm, hEffVsAngle, hEffVsEta, hEffVsPhi, hEffVsPos, hEtaRecHit, hEtaSimSegm, hPhiRecHit, hPhiSimSegm, hPosRecHit, hPosSimSegm, and name.

                                           {
       name=name_;
       hEtaSimSegm = (TH1F *) file->Get("2D_"+name+"_hEtaSimSegm");
       hEtaRecHit = (TH1F *) file->Get("2D_"+name+"_hEtaRecHit");
       hEffVsEta = (TH1F *) file->Get("2D_"+name+"_hEffVsEta");
 
       hPhiSimSegm = (TH1F *) file->Get("2D_"+name+"_hPhiSimSegm");
       hPhiRecHit = (TH1F *) file->Get("2D_"+name+"_hPhiRecHit");
       hEffVsPhi = (TH1F *) file->Get("2D_"+name+"_hEffVsPhi");
 
       hPosSimSegm = (TH1F *) file->Get("2D_"+name+"_hPosSimSegm");
       hPosRecHit = (TH1F *) file->Get("2D_"+name+"_hPosRecHit");
       hEffVsPos = (TH1F *) file->Get("2D_"+name+"_hEffVsPos");
 
       hAngleSimSegm = (TH1F *) file->Get("2D_"+name+"_hAngleSimSegm");
       hAngleRecHit = (TH1F *) file->Get("2D_"+name+"_hAngleRecHit");
       hEffVsAngle = (TH1F *) file->Get("2D_"+name+"_hEffVsAngle");
     }

HEff2DHit::~HEff2DHit ( )

inline

Definition at line 480 of file Histograms.h.

                 {
       //delete hEtaSimSegm;    
       //delete hEtaRecHit;     
       //if(hEffVsEta != 0)
       // delete hEffVsEta;   
 
       //delete hPhiSimSegm;    
       //delete hPhiRecHit;
       //if(hEffVsPhi != 0)
       //  delete hEffVsPhi;      
 
       //delete hPosSimSegm;    
       //delete hPosRecHit;
       //if(hEffVsPos != 0)
       //  delete hEffVsPos;      
 
       //delete hAngleSimSegm;  
       //delete hAngleRecHit;
       //if(hEffVsAngle != 0)
       //  delete hEffVsAngle;    
     }

Member Function Documentation

void HEff2DHit::ComputeEfficiency ( )

inline

Definition at line 523 of file Histograms.h.

References newFWLiteAna::bin, error, hAngleRecHit, hAngleSimSegm, hEffVsAngle, hEffVsEta, hEffVsPhi, hEffVsPos, hEtaRecHit, hEtaSimSegm, hPhiRecHit, hPhiSimSegm, hPosRecHit, hPosSimSegm, name, and mathSSE::sqrt().

                              {
 
       hEffVsEta = (TH1F *) hEtaRecHit->Clone();
       hEffVsEta->SetName("2D_"+name+"_hEffVsEta");
       hEffVsEta->SetTitle("2D RecHit Efficiency as a function of Eta");
       hEffVsEta->Divide(hEtaSimSegm);
       // Set the error accordingly to binomial statistics
       int nBinsEta = hEffVsEta->GetNbinsX();
       for(int bin = 1; bin <=  nBinsEta; bin++) {
         float nSimHit = hEtaSimSegm->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("2D_"+name+"_hEffVsPhi");
       hEffVsPhi->SetTitle("2D RecHit Efficiency as a function of Phi");
       hEffVsPhi->Divide(hPhiSimSegm);
       // Set the error accordingly to binomial statistics
       int nBinsPhi = hEffVsPhi->GetNbinsX();
       for(int bin = 1; bin <=  nBinsPhi; bin++) {
         float nSimHit = hPhiSimSegm->GetBinContent(bin);
         float eff = hEffVsPhi->GetBinContent(bin);
         float error = 0;
         if(nSimHit != 0) {
           error = sqrt(eff*(1-eff)/nSimHit);
         }
         hEffVsPhi->SetBinError(bin, error);
       }
 
       hEffVsPos = (TH1F *) hPosRecHit->Clone();
       hEffVsPos->SetName("2D_"+name+"_hEffVsPos");
       hEffVsPos->SetTitle("2D RecHit Efficiency as a function of position in SL");
       hEffVsPos->Divide(hPosSimSegm);
       // Set the error accordingly to binomial statistics
       int nBinsPos = hEffVsPos->GetNbinsX();
       for(int bin = 1; bin <=  nBinsPos; bin++) {
         float nSimHit = hPosSimSegm->GetBinContent(bin);
         float eff = hEffVsPos->GetBinContent(bin);
         float error = 0;
         if(nSimHit != 0) {
           error = sqrt(eff*(1-eff)/nSimHit);
         }
         hEffVsPos->SetBinError(bin, error);
       }
 
       hEffVsAngle = (TH1F *) hAngleRecHit->Clone();
       hEffVsAngle->SetName("2D_"+name+"_hEffVsAngle");
       hEffVsAngle->SetTitle("2D RecHit Efficiency as a function of angle");
       hEffVsAngle->Divide(hAngleSimSegm);
       // Set the error accordingly to binomial statistics
       int nBinsAngle = hEffVsAngle->GetNbinsX();
       for(int bin = 1; bin <=  nBinsAngle; bin++) {
         float nSimHit = hAngleSimSegm->GetBinContent(bin);
         float eff = hEffVsAngle->GetBinContent(bin);
         float error = 0;
         if(nSimHit != 0) {
           error = sqrt(eff*(1-eff)/nSimHit);
         }
         hEffVsAngle->SetBinError(bin, error);
       }
 
     }

void HEff2DHit::Fill	(	float	etaSimSegm,
		float	phiSimSegm,
		float	posSimSegm,
		float	angleSimSegm,
		bool	fillRecHit
	)

inline

Definition at line 502 of file Histograms.h.

References hAngleRecHit, hAngleSimSegm, hEtaRecHit, hEtaSimSegm, hPhiRecHit, hPhiSimSegm, hPosRecHit, and hPosSimSegm.

Referenced by DTSegment2DQuality::analyze().

                                {
 
       hEtaSimSegm->Fill(etaSimSegm);       
       hPhiSimSegm->Fill(phiSimSegm);   
       hPosSimSegm->Fill(posSimSegm);   
       hAngleSimSegm->Fill(angleSimSegm); 
 
       if(fillRecHit) {
         hEtaRecHit->Fill(etaSimSegm);    
         hPhiRecHit->Fill(phiSimSegm);    
         hPosRecHit->Fill(posSimSegm);    
         hAngleRecHit->Fill(angleSimSegm);  
       }
     }

void HEff2DHit::Write ( )

inline

Definition at line 591 of file Histograms.h.

References hAngleRecHit, hAngleSimSegm, hEffVsAngle, hEffVsEta, hEffVsPhi, hEffVsPos, hEtaRecHit, hEtaSimSegm, hPhiRecHit, hPhiSimSegm, hPosRecHit, and hPosSimSegm.

                  {
       hEtaSimSegm->Write();   
       hEtaRecHit->Write();    
       if(hEffVsEta != 0)
         hEffVsEta->Write();     
       hPhiSimSegm->Write();   
       hPhiRecHit->Write();    
       if(hEffVsPhi != 0)
         hEffVsPhi->Write();     
       hPosSimSegm->Write();   
       hPosRecHit->Write();    
       if(hEffVsPos != 0)
         hEffVsPos->Write();     
       hAngleSimSegm->Write(); 
       hAngleRecHit->Write();  
       if(hEffVsAngle != 0)
         hEffVsAngle->Write();   
 
     }