#include <Histograms.h>

Inheritance diagram for HMassResolutionVSPart:

Public Member Functions
virtual void	Clear ()

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)

	HMassResolutionVSPart (TFile *outputFile, const TString &name)

virtual void	Write ()

	~HMassResolutionVSPart ()

Public Member Functions inherited from Histograms
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 int charge, const double &weight=1.)

virtual void	Fill (const CLHEP::HepLorentzVector &momentum, const int charge, 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 unsigned int number)

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.)

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)

virtual void	SetWeight (double weight)

virtual	~Histograms ()

Protected Attributes
std::map< TString, TH1 * >	mapHisto_

std::auto_ptr< HDelta >	muMinus

std::auto_ptr< HDelta >	muPlus

TString	nameSuffix_ [2]

Protected Attributes inherited from Histograms
TDirectory *	histoDir_

TString	name_

TFile *	outputFile_

double	theWeight_

Detailed Description

A set of histograms for resolution. The fill method requires the two selected muons, their charges and the reconstructed and generated masses. It evaluates the resolution on the mass vs: Pt of the pair DeltaEta of the pair DeltaPhi of the pair pt, eta and phi of the plus and minus muon separately

Definition at line 2144 of file Histograms.h.

Constructor & Destructor Documentation

HMassResolutionVSPart::HMassResolutionVSPart	(	TFile *	outputFile,
		const TString &	name
	)

inline

Definition at line 2147 of file Histograms.h.

References Histograms::histoDir_, i, mapHisto_, muMinus, muPlus, mergeVDriftHistosByStation::name, and nameSuffix_.

                                                                   : Histograms( outputFile, name ) {
     // Kinematical variables
     nameSuffix_[0] = "Plus";
     nameSuffix_[1] = "Minus";
     TString titleSuffix[] = {" for mu+", " for mu-"};
 
     mapHisto_[name]                  = new TH1F (name, "#Delta M/M", 4000, -1, 1);
     mapHisto_[name+"VSPairPt"]       = new HResolution (name+"VSPairPt", "resolution VS pt of the pair", 100, 0, 200, -1, 1, histoDir_);
     mapHisto_[name+"VSPairDeltaEta"] = new HResolution (name+"VSPairDeltaEta", "resolution VS #Delta#eta of the pair", 100, -0.1, 6.2, -1, 1, histoDir_);
     mapHisto_[name+"VSPairDeltaPhi"] = new HResolution (name+"VSPairDeltaPhi", "resolution VS #Delta#phi of the pair", 100, -0.1, 3.2, -1, 1, histoDir_);
 
     for( int i=0; i<2; ++i ) {
       mapHisto_[name+"VSPt"+nameSuffix_[i]]  = new HResolution (name+"VSPt"+nameSuffix_[i], "resolution VS pt"+titleSuffix[i], 100, 0, 200, -1, 1, histoDir_);
       mapHisto_[name+"VSEta"+nameSuffix_[i]] = new HResolution (name+"VSEta"+nameSuffix_[i], "resolution VS #eta"+titleSuffix[i], 100, -3, 3, -1, 1, histoDir_);
       mapHisto_[name+"VSPhi"+nameSuffix_[i]] = new HResolution (name+"VSPhi"+nameSuffix_[i], "resolution VS #phi"+titleSuffix[i], 100, -3.2, 3.2, -1, 1, histoDir_);
     }
 
     // single particles histograms
     muMinus.reset( new HDelta("muMinus") );
     muPlus.reset( new HDelta("muPlus") );
   }

HMassResolutionVSPart::~HMassResolutionVSPart ( )

inline

Definition at line 2169 of file Histograms.h.

References timingPdfMaker::histo, and mapHisto_.

                           {
     for (std::map<TString, TH1*>::const_iterator histo=mapHisto_.begin(); 
          histo!=mapHisto_.end(); histo++) {
       delete (*histo).second;
     }
   }

Member Function Documentation

virtual void HMassResolutionVSPart::Clear ( )

inlinevirtual

Implements Histograms.

Definition at line 2247 of file Histograms.h.

References timingPdfMaker::histo, mapHisto_, muMinus, and muPlus.

                        {
     for (std::map<TString, TH1*>::const_iterator histo=mapHisto_.begin(); 
          histo!=mapHisto_.end(); histo++) {
       (*histo).second->Clear();
     }
     muMinus->Clear();
     muPlus->Clear();
   }

virtual void HMassResolutionVSPart::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
	)

inlinevirtual

Reimplemented from Histograms.

Definition at line 2176 of file Histograms.h.

References muMinus, and muPlus.

                                                                        {
     muMinus->Fill(recoP1, genP1);
     muPlus->Fill(recoP2, genP2);
 
     Fill( recoP1, charge1, recoP2, charge2, recoMass, genMass );
   }

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

inlinevirtual

Reimplemented from Histograms.

Definition at line 2187 of file Histograms.h.

References gather_cfg::cout, MuScleFitUtils::deltaPhi(), i, cmsHarvester::index, mapHisto_, Histograms::name_, and nameSuffix_.

                                                                        {
 
     if ( charge1 == charge2 ) std::cout << "Error: must get two opposite charge particles" << std::endl;
 
     double massRes = (recoMass - genMass)/genMass;
 
     reco::Particle::LorentzVector recoPair( recoP1 + recoP2 );
     double pairPt = recoPair.Pt();
 
     double recoPt[2]  = {recoP1.Pt(),  recoP2.Pt()};
     double recoEta[2] = {recoP1.Eta(), recoP2.Eta()};
     double recoPhi[2] = {recoP1.Phi(), recoP2.Phi()};
 
     // std::cout << "pairPt = " << pairPt << ", massRes = ("<<recoMass<<" - "<<genMass<<")/"<<genMass<<" = " << massRes
     //      << ", recoPt[0] = " << recoPt[0] << ", recoPt[1] = " << recoPt[1] << std::endl;
 
     // Index of the histogram. If the muons have charge1 = -1 and charge2 = 1, they already have the
     // correct histogram indeces. Otherwise, swap the indeces.
     // In any case the negative muon has index = 0  and the positive muon has index = 1.
     int id[2] = {0,1};
     if ( charge1 > 0 ) {
       id[0] = 1;
       id[1] = 0;
     }
 
     double pairDeltaEta = fabs(recoEta[0] - recoEta[1]);
     double pairDeltaPhi = MuScleFitUtils::deltaPhi(recoPhi[0], recoPhi[1]);
 
     mapHisto_[name_]->Fill(massRes);
     mapHisto_[name_+"VSPairPt"]->Fill(pairPt, massRes);
     mapHisto_[name_+"VSPairDeltaEta"]->Fill(pairDeltaEta, massRes);
     mapHisto_[name_+"VSPairDeltaPhi"]->Fill(pairDeltaPhi, massRes);
 
     // Resolution vs single muon quantities
     // ------------------------------------
     int index = 0;
     for( int i=0; i<2; ++i ) {
       index = id[i];
       mapHisto_[name_+"VSPt"+nameSuffix_[i]]->Fill(recoPt[index], massRes); // EM [index] or [i]???
       mapHisto_[name_+"VSEta"+nameSuffix_[i]]->Fill(recoEta[index], massRes);
       mapHisto_[name_+"VSPhi"+nameSuffix_[i]]->Fill(recoPhi[index], massRes);
     }
   } 

virtual void HMassResolutionVSPart::Write ( )

inlinevirtual

Implements Histograms.

Definition at line 2235 of file Histograms.h.

References timingPdfMaker::histo, Histograms::histoDir_, mapHisto_, muMinus, and muPlus.

                        {
     histoDir_->cd();
     for (std::map<TString, TH1*>::const_iterator histo=mapHisto_.begin(); 
          histo!=mapHisto_.end(); histo++) {
       (*histo).second->Write();
     }
     // Create the new dir and cd into it
     (histoDir_->mkdir("singleMuonsVSgen"))->cd();
     muMinus->Write();
     muPlus->Write();
   }