Public Member Functions
void	bookHistograms (DQMStore &dqmStore)

void	fillHistograms (int numJets, const edm::Event &evt, double evtWeight)

	leptonEfficiencyT (int minJets, int maxJets, const edm::InputTag &srcGen, const std::string &cutGen, const edm::InputTag &srcRec, const std::string &cutRec, double dRmatch, const std::string &dqmDirectory)

	~leptonEfficiencyT ()

Public Attributes
StringCutObjectSelector < reco::Candidate > *	cutGen_

StringCutObjectSelector< T > *	cutRec_

std::string	dqmDirectory_

double	dRmatch_

MonitorElement *	histogramDenominatorEta_

MonitorElement *	histogramDenominatorPhi_

MonitorElement *	histogramDenominatorPt_

MonitorElement *	histogramNumeratorEta_

MonitorElement *	histogramNumeratorPhi_

MonitorElement *	histogramNumeratorPt_

MonitorElement *	histogramNumGenLeptons_

MonitorElement *	histogramNumRecLeptons_

int	maxJets_

int	minJets_

edm::InputTag	srcGen_

edm::InputTag	srcRec_

Detailed Description

template<typename T>
struct MCEmbeddingValidationAnalyzer::leptonEfficiencyT< T >

Definition at line 1072 of file MCEmbeddingValidationAnalyzer.h.

Constructor & Destructor Documentation

template<typename T>

MCEmbeddingValidationAnalyzer::leptonEfficiencyT< T >::leptonEfficiencyT	(	int	minJets,
		int	maxJets,
		const edm::InputTag &	srcGen,
		const std::string &	cutGen,
		const edm::InputTag &	srcRec,
		const std::string &	cutRec,
		double	dRmatch,
		const std::string &	dqmDirectory
	)

inline

Definition at line 1074 of file MCEmbeddingValidationAnalyzer.h.

References MCEmbeddingValidationAnalyzer::leptonEfficiencyT< T >::cutGen_, MCEmbeddingValidationAnalyzer::leptonEfficiencyT< T >::cutRec_, MCEmbeddingValidationAnalyzer::leptonEfficiencyT< T >::dqmDirectory_, MCEmbeddingValidationAnalyzer::leptonEfficiencyT< T >::maxJets_, and MCEmbeddingValidationAnalyzer::leptonEfficiencyT< T >::minJets_.

       : minJets_(minJets),
     maxJets_(maxJets),
     srcGen_(srcGen),
     cutGen_(0),
     srcRec_(srcRec),
     cutRec_(0),
     dRmatch_(dRmatch)
     {
       if ( cutGen != "" ) cutGen_ = new StringCutObjectSelector<reco::Candidate>(cutGen);
       if ( cutRec != "" ) cutRec_ = new StringCutObjectSelector<T>(cutRec);
       dqmDirectory_ = dqmDirectory;
       if      ( minJets_ < 0 && maxJets_ < 0 ) dqmDirectory_.append("");
       else if (                 maxJets_ < 0 ) dqmDirectory_.append(Form("_numJetsGe%i", minJets_));
       else if ( minJets_ < 0                 ) dqmDirectory_.append(Form("_numJetsLe%i", maxJets_));
       else if ( maxJets_     == minJets_     ) dqmDirectory_.append(Form("_numJetsEq%i", minJets_));
       else                                     dqmDirectory_.append(Form("_numJets%ito%i", minJets_, maxJets_));
     }

template<typename T>

MCEmbeddingValidationAnalyzer::leptonEfficiencyT< T >::~leptonEfficiencyT ( )

inline

Definition at line 1093 of file MCEmbeddingValidationAnalyzer.h.

References MCEmbeddingValidationAnalyzer::leptonEfficiencyT< T >::cutGen_, and MCEmbeddingValidationAnalyzer::leptonEfficiencyT< T >::cutRec_.

     {
       delete cutGen_;
       delete cutRec_;
     }

Member Function Documentation

template<typename T>

void MCEmbeddingValidationAnalyzer::leptonEfficiencyT< T >::bookHistograms ( DQMStore & dqmStore )

inline

Definition at line 1098 of file MCEmbeddingValidationAnalyzer.h.

     {
       dqmStore.setCurrentFolder(dqmDirectory_.data());
       histogramNumGenLeptons_ = dqmStore.book1D("numGenLeptons", "numGenLeptons", 10, -0.5, 9.5);
       histogramNumRecLeptons_ = dqmStore.book1D("numRecLeptons", "numRecLeptons", 10, -0.5, 9.5);
       histogramNumeratorPt_ = dqmStore.book1D("numeratorPt", "numeratorPt", 250, 0., 250.);
       histogramDenominatorPt_ = dqmStore.book1D("denominatorPt", "denominatorPt", 250, 0., 250.);
       histogramNumeratorEta_ = dqmStore.book1D("numeratorEta", "numeratorEta", 198, -9.9, +9.9);
       histogramDenominatorEta_ = dqmStore.book1D("denominatorEta", "denominatorEta", 198, -9.9, +9.9);
       histogramNumeratorPhi_ = dqmStore.book1D("numeratorPhi", "numeratorPhi", 72, -TMath::Pi(), +TMath::Pi());
       histogramDenominatorPhi_ = dqmStore.book1D("denominatorPhi", "denominatorPhi", 72, -TMath::Pi(), +TMath::Pi());
     }

template<typename T>

void MCEmbeddingValidationAnalyzer::leptonEfficiencyT< T >::fillHistograms	(	int	numJets,
		const edm::Event &	evt,
		double	evtWeight
	)

inline

Definition at line 1110 of file MCEmbeddingValidationAnalyzer.h.

     {
       if ( (minJets_ == -1 || numJets >= minJets_) &&
        (maxJets_ == -1 || numJets <= maxJets_) ) {
     typedef edm::View<reco::Candidate> CandidateView;
     edm::Handle<CandidateView> genLeptons;
     evt.getByLabel(srcGen_, genLeptons);
     histogramNumGenLeptons_->Fill(genLeptons->size(), evtWeight);
     typedef std::vector<T> recLeptonCollection;
     edm::Handle<recLeptonCollection> recLeptons;
     evt.getByLabel(srcRec_, recLeptons);
     histogramNumRecLeptons_->Fill(recLeptons->size(), evtWeight);
     for ( CandidateView::const_iterator genLepton = genLeptons->begin();
           genLepton != genLeptons->end(); ++genLepton ) {
       if ( cutGen_ && !(*cutGen_)(*genLepton) ) continue;   
       bool isMatched = false;
       for ( typename recLeptonCollection::const_iterator recLepton = recLeptons->begin();
         recLepton != recLeptons->end(); ++recLepton ) {
         if ( cutRec_ && !(*cutRec_)(*recLepton) ) continue;
         double dR = deltaR(genLepton->p4(), recLepton->p4());
         if ( dR < dRmatch_ ) isMatched = true;
       }
       histogramDenominatorPt_->Fill(genLepton->pt(), evtWeight);
       histogramDenominatorEta_->Fill(genLepton->eta(), evtWeight);
       histogramDenominatorPhi_->Fill(genLepton->phi(), evtWeight);
       if ( isMatched ) {
         histogramNumeratorPt_->Fill(genLepton->pt(), evtWeight);
         histogramNumeratorEta_->Fill(genLepton->eta(), evtWeight);
         histogramNumeratorPhi_->Fill(genLepton->phi(), evtWeight);
       }
     }
       }
     }