#include <Histograms.h>

Inheritance diagram for HParticle:

Public Member Functions
virtual void	Clear ()

virtual void	Fill (const reco::Particle::LorentzVector &p4, const int charge, const double &weight=1.)

virtual void	Fill (const CLHEP::HepLorentzVector &momentum, const int charge, const double &weight=1.)

virtual void	Fill (unsigned int number)

	HParticle (const TString &name, const double &minMass=0., const double &maxMass=200., const double &maxPt=100.)

	HParticle (TFile *outputFile, const TString &name, const double &minMass=0., const double &maxMass=200., const double &maxPt=100.)
	Constructor that puts the histograms inside a TDirectory. More...

	HParticle (const TString &name, TFile *file)

virtual void	Write ()

	~HParticle ()

Public Member Functions inherited from Histograms
void	declareHistograms ()

virtual void	Fill (const reco::Particle::LorentzVector &p1, const reco::Particle::LorentzVector &p2)

virtual void	Fill (const reco::Particle::LorentzVector &p1, const reco::Particle::LorentzVector &p2, const int charge, const double &weight=1.)

virtual void	Fill (const CLHEP::HepLorentzVector &momentum1, const CLHEP::HepLorentzVector &momentum2)

virtual void	Fill (const CLHEP::HepLorentzVector &momentum1, const CLHEP::HepLorentzVector &momentum2, const int charge, const double &weight=1.)

virtual void	Fill (const CLHEP::HepLorentzVector &p1, const reco::Particle::LorentzVector &p2)

virtual void	Fill (const reco::Particle::LorentzVector &p4, const double &weight=1.)

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

virtual void	Fill (const reco::Particle::LorentzVector &p4, const double &genValue, const double recValue, const int charge)

virtual void	Fill (const CLHEP::HepLorentzVector &p, const double &likeValue)

virtual void	Fill (const reco::Particle::LorentzVector &recoP1, const int charge1, const reco::Particle::LorentzVector &genP1, const reco::Particle::LorentzVector &recoP2, const int charge2, const reco::Particle::LorentzVector &genP2, const double &recoMass, const double &genMass)

virtual void	Fill (const reco::Particle::LorentzVector &recoP1, const int charge1, const reco::Particle::LorentzVector &recoP2, const int charge2, const double &recoMass, const double &genMass)

virtual void	Fill (const reco::Particle::LorentzVector &recoP1, const reco::Particle::LorentzVector &genP1, const reco::Particle::LorentzVector &recoP2, const reco::Particle::LorentzVector &genP2)

virtual void	Fill (const double &x, const double &y)

virtual void	Fill (const double &x, const double &y, const double &a, const double &b)

virtual void	Fill (const reco::Particle::LorentzVector &p41, const reco::Particle::LorentzVector &p42, const reco::Particle::LorentzVector &p4Res, const double &weight=1.)

virtual void	Fill (const CLHEP::HepLorentzVector &momentum1, const CLHEP::HepLorentzVector &momentum2, const CLHEP::HepLorentzVector &momentumRes, const double &weight=1.)

void	fillEventInfo (int proc, int strk, int ntrkr)

void	fillRecHistograms (const RecTrack_t &r)

void	fillSimHistograms (const SimTrack_t &s)

void	fillVzeroHistograms (const RecVzero_t &r, int part)

virtual double	Get (const reco::Particle::LorentzVector &recoP1, const TString &covarianceName)

virtual TString	GetName ()

	Histograms ()

	Histograms (const TString &name)

	Histograms (TFile *outputFile, const TString &name)

	Histograms (const edm::ParameterSet &pset)

virtual void	SetWeight (double weight)

void	writeHistograms ()

virtual	~Histograms ()

	~Histograms ()

Protected Attributes
TProfile *	hCurvVsEtaNeg_

TProfile *	hCurvVsEtaPos_

TProfile *	hCurvVsPhiNeg_

TProfile *	hCurvVsPhiPos_

TH1F *	hEta_

TH1F *	hMass_

TH1F *	hNumber_

TH1F *	hPhi_

TH1F *	hPt_

TH2F *	hPtVsEta_

TProfile *	hPtVsPhiNeg_

TProfile *	hPtVsPhiPos_

Protected Attributes inherited from Histograms
TDirectory *	histoDir_

TString	name_

TFile *	outputFile_

double	theWeight_

Detailed Description

Definition at line 238 of file Histograms.h.

Constructor & Destructor Documentation

HParticle::HParticle	(	const TString &	name,
		const double &	minMass = `0.`,
		const double &	maxMass = `200.`,
		const double &	maxPt = `100.`
	)

inline

Definition at line 240 of file Histograms.h.

                                                                                                                              :
     Histograms(name),
     // Kinematical variables
     hPt_(      new TH1F (name+"_Pt",      "transverse momentum", 100, 0, maxPt) ),
     hPtVsEta_( new TH2F (name+"_PtVsEta", "transverse momentum vs #eta", 100, 0, maxPt, 100, -3.0, 3.0) ),
    
     hCurvVsEtaNeg_( new TProfile(name+"_CurvVsEtaNeg", "q/pT vs #eta neg.", 64, -3.2, 3.2, -1., 0.) ),
     hCurvVsEtaPos_( new TProfile(name+"_CurvVsEtaPos", "q/pT vs #eta pos.", 64, -3.2, 3.2,  0., 1.) ),
     hCurvVsPhiNeg_( new TProfile(name+"_CurvVsPhiNeg", "q/pT vs #phi neg.", 32, -3.2, 3.2, -1., 0.) ),
     hCurvVsPhiPos_( new TProfile(name+"_CurvVsPhiPos", "q/pT vs #phi pos.", 32, -3.2, 3.2,  0., 1.) ),
 
     hPtVsPhiNeg_( new TProfile(name+"_PtVsPhiNeg", "pT vs #phi neg.", 32, -3.2, 3.2, 0.,100) ),
     hPtVsPhiPos_( new TProfile(name+"_PtVsPhiPos", "pT vs #phi pos.", 32, -3.2, 3.2, 0.,100) ),
 
     
     hEta_(     new TH1F (name+"_Eta",     "pseudorapidity", 64, -3.2, 3.2) ),
     hPhi_(     new TH1F (name+"_Phi",     "phi angle",      64, -3.2, 3.2) ),
     hMass_(    new TH1F (name+"_Mass",    "mass", 10000, minMass, maxMass) ),
     hNumber_( new TH1F (name+"_Number", "number", 20, -0.5, 19.5) )
   {}

HParticle::HParticle	(	TFile *	outputFile,
		const TString &	name,
		const double &	minMass = `0.`,
		const double &	maxMass = `200.`,
		const double &	maxPt = `100.`
	)

inline

Constructor that puts the histograms inside a TDirectory.

Definition at line 262 of file Histograms.h.

References hCurvVsEtaNeg_, hCurvVsEtaPos_, hCurvVsPhiNeg_, hCurvVsPhiPos_, hEta_, hMass_, hNumber_, hPhi_, hPt_, hPtVsEta_, hPtVsPhiNeg_, and hPtVsPhiPos_.

                                                                                                                                                 :
     Histograms(outputFile, name)
   {
     // Kinematical variables
     hPt_ =      new TH1F (name+"_Pt", "transverse momentum", 100, 0, maxPt);
     hPtVsEta_ = new TH2F (name+"_PtVsEta", "transverse momentum vs #eta", 100, 0, maxPt, 100, -3.0, 3.0);
     
     hPtVsEta_ = new TH2F (name+"_PtVsEta", "transverse momentum vs #eta", 100, 0, maxPt, 100, -3.0, 3.0);
    
     hCurvVsEtaNeg_ = new TProfile(name+"_CurvVsEtaNeg", "q/pT vs #eta neg.", 100, -3.0, 3.0, -1. ,0.);
     hCurvVsEtaPos_ = new TProfile(name+"_CurvVsEtaPos", "q/pT vs #eta pos.", 100, -3.0, 3.0, 0., 1.);
     hCurvVsPhiNeg_ = new TProfile(name+"_CurvVsPhiNeg", "q/pT vs #phi neg.", 32, -3.2, 3.2, -1. ,0.);
     hCurvVsPhiPos_ = new TProfile(name+"_CurvVsPhiPos", "q/pT vs #phi pos.", 32, -3.2, 3.2, 0., 1.);
 
     hPtVsPhiNeg_ = new TProfile(name+"_PtVsPhiNeg", "pT vs #phi neg.", 32, -3.2, 3.2, 0.,100);
     hPtVsPhiPos_ = new TProfile(name+"_PtVsPhiPos", "pT vs #phi pos.", 32, -3.2, 3.2, 0.,100);
 
 
     //hPtVSPhi_prof_ = new TProfile (name+"_PtVSPhi_prof", "pt vs phi angle",12, -3.2, 3.2, 0, 200);
 
     hEta_ =     new TH1F (name+"_Eta", "pseudorapidity", 64, -3.2, 3.2);
     hPhi_ =     new TH1F (name+"_Phi", "phi angle",      64, -3.2, 3.2);
     hMass_ =    new TH1F (name+"_Mass", "mass", 40000, minMass, maxMass);
     hNumber_ = new TH1F (name+"_Number", "number", 20, -0.5, 19.5);
   }

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

inline

Definition at line 288 of file Histograms.h.

                                                  :
     Histograms(name),
     hPt_(      (TH1F *) file->Get(name_+"_Pt") ),
     hPtVsEta_( (TH2F *) file->Get(name_+"_PtVsEta") ),
 
    
     hCurvVsEtaNeg_( (TProfile *) file->Get(name_+"_CurvVsEtaNeg") ),   
     hCurvVsEtaPos_( (TProfile *) file->Get(name_+"_CurvVsEtaPos") ),
     hCurvVsPhiNeg_( (TProfile *) file->Get(name_+"_CurvVsPhiNeg") ),   
     hCurvVsPhiPos_( (TProfile *) file->Get(name_+"_CurvVsPhiPos") ),
 
     hPtVsPhiNeg_( (TProfile *) file->Get(name_+"_PtVsPhiNeg") ),   
     hPtVsPhiPos_( (TProfile *) file->Get(name_+"_PtVsPhiPos") ),
 
     hEta_(     (TH1F *) file->Get(name_+"_Eta") ),
     hPhi_(     (TH1F *) file->Get(name_+"_Phi") ),
     hMass_(    (TH1F *) file->Get(name_+"_Mass") ),
     //hMass_fine_ = (TH1F *) file->Get(name_+"_Mass_fine");
     hNumber_( (TH1F *) file->Get(name_+"_Number") )
   {}

HParticle::~HParticle ( )

inline

Definition at line 309 of file Histograms.h.

References hCurvVsEtaNeg_, hCurvVsEtaPos_, hCurvVsPhiNeg_, hCurvVsPhiPos_, hEta_, hMass_, hNumber_, hPhi_, hPt_, hPtVsEta_, hPtVsPhiNeg_, and hPtVsPhiPos_.

   {
     delete hPt_;
     delete hPtVsEta_;
 
     delete hCurvVsEtaNeg_; 
     delete hCurvVsEtaPos_;
     delete hCurvVsPhiNeg_; 
     delete hCurvVsPhiPos_;
 
     delete hPtVsPhiNeg_; 
     delete hPtVsPhiPos_;
 
     delete hEta_;
     delete hPhi_;
     delete hMass_;
     // delete hMass_fine_;
     delete hNumber_;
   }

Member Function Documentation

virtual void HParticle::Clear ( )

inlinevirtual

Implements Histograms.

Definition at line 382 of file Histograms.h.

References hCurvVsEtaNeg_, hCurvVsEtaPos_, hCurvVsPhiNeg_, hCurvVsPhiPos_, hEta_, hMass_, hNumber_, hPhi_, hPt_, hPtVsEta_, hPtVsPhiNeg_, and hPtVsPhiPos_.

   {
     hPt_->Clear();
     hPtVsEta_->Clear();
     
     hCurvVsEtaNeg_->Clear();
     hCurvVsEtaPos_->Clear();
     hCurvVsPhiNeg_->Clear();
     hCurvVsPhiPos_->Clear();
     
     hPtVsPhiNeg_->Clear();
     hPtVsPhiPos_->Clear();
 
     hEta_->Clear();    
     hPhi_->Clear();
     hMass_->Clear();
     //hMass_fine_->Clear();
     hNumber_->Clear();
   }

virtual void HParticle::Fill	(	const reco::Particle::LorentzVector &	p4,
		const int	charge,
		const double &	weight = `1.`
	)

inlinevirtual

Reimplemented from Histograms.

Definition at line 329 of file Histograms.h.

References CommonMethods::weight().

   {
     Fill(CLHEP::HepLorentzVector(p4.x(),p4.y(),p4.z(),p4.t()),charge, weight);
   }

virtual void HParticle::Fill	(	const CLHEP::HepLorentzVector &	momentum,
		const int	charge,
		const double &	weight = `1.`
	)

inlinevirtual

Reimplemented from Histograms.

Definition at line 334 of file Histograms.h.

References hCurvVsEtaNeg_, hCurvVsEtaPos_, hCurvVsPhiNeg_, hCurvVsPhiPos_, hEta_, hMass_, hPhi_, hPt_, hPtVsEta_, hPtVsPhiNeg_, hPtVsPhiPos_, and CommonMethods::weight().

   {
     hPt_->Fill(momentum.perp(), weight);
     hPtVsEta_->Fill(momentum.perp(), momentum.eta(), weight);
    
     //   std::cout<< "charge-> " <<charge<<std::endl;
     if(charge<0)hCurvVsEtaNeg_->Fill( momentum.eta(),charge/(momentum.perp()),weight);     
     if(charge>0)hCurvVsEtaPos_->Fill( momentum.eta(),charge/(momentum.perp()),weight);
     if(charge<0)hCurvVsPhiNeg_->Fill( momentum.phi(),charge/(momentum.perp()),weight);     
     if(charge>0)hCurvVsPhiPos_->Fill( momentum.phi(),charge/(momentum.perp()),weight);
 
     if(charge<0)hPtVsPhiNeg_->Fill( momentum.phi(),momentum.perp(),weight);     
     if(charge>0)hPtVsPhiPos_->Fill( momentum.phi(),momentum.perp(),weight);
 
     hEta_->Fill(momentum.eta(), weight);
     hPhi_->Fill(momentum.phi(), weight);
     hMass_->Fill(momentum.m(), weight);
     //hMass_fine_->Fill(momentum.m(), weight);
   }

virtual void HParticle::Fill ( unsigned int number )

inlinevirtual

Reimplemented from Histograms.

Definition at line 356 of file Histograms.h.

References hNumber_.

   {
     hNumber_->Fill(number);
   }

virtual void HParticle::Write ( )

inlinevirtual

Implements Histograms.

Definition at line 361 of file Histograms.h.

References hCurvVsEtaNeg_, hCurvVsEtaPos_, hCurvVsPhiNeg_, hCurvVsPhiPos_, hEta_, Histograms::histoDir_, hMass_, hNumber_, hPhi_, hPt_, hPtVsEta_, hPtVsPhiNeg_, and hPtVsPhiPos_.

   {
     if(histoDir_ != 0) histoDir_->cd();
     hPt_->Write();
     hPtVsEta_->Write();
    
     hCurvVsEtaNeg_->Write();
     hCurvVsEtaPos_->Write();
     hCurvVsPhiNeg_->Write();
     hCurvVsPhiPos_->Write();
 
     hPtVsPhiNeg_->Write();
     hPtVsPhiPos_->Write();
 
     hEta_->Write();    
     hPhi_->Write();
     hMass_->Write();
     //hMass_fine_->Write();
     hNumber_->Write();
   }