Public Member Functions
void	bookHistograms (DQMStore &dqmStore)

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

	tauDistributionExtra (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)

	~tauDistributionExtra ()

Public Attributes
StringCutObjectSelector < reco::Candidate > *	cutGen_

StringCutObjectSelector < pat::Tau > *	cutRec_

std::string	dqmDirectory_

double	dRmatch_

MonitorElement *	histogramGenTauDecayMode_

MonitorElement *	histogramRecTauDecayMode_

int	maxJets_

int	minJets_

edm::InputTag	srcGen_

edm::InputTag	srcRec_

Detailed Description

Definition at line 959 of file MCEmbeddingValidationAnalyzer.h.

Constructor & Destructor Documentation

MCEmbeddingValidationAnalyzer::tauDistributionExtra::tauDistributionExtra	(	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 961 of file MCEmbeddingValidationAnalyzer.h.

References cutGen_, cutRec_, dqmDirectory_, maxJets_, and 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<pat::Tau>(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_));
     }

MCEmbeddingValidationAnalyzer::tauDistributionExtra::~tauDistributionExtra ( )

inline

Definition at line 980 of file MCEmbeddingValidationAnalyzer.h.

References cutGen_, and cutRec_.

     {
       delete cutGen_;
       delete cutRec_;
     }

Member Function Documentation

void MCEmbeddingValidationAnalyzer::tauDistributionExtra::bookHistograms ( DQMStore & dqmStore )

inline

Definition at line 985 of file MCEmbeddingValidationAnalyzer.h.

References DQMStore::book1D(), dqmDirectory_, histogramGenTauDecayMode_, histogramRecTauDecayMode_, and DQMStore::setCurrentFolder().

     {
       dqmStore.setCurrentFolder(dqmDirectory_.data());
       histogramGenTauDecayMode_ = dqmStore.book1D("genTauDecayMode", "genTauDecayMode", 21, -1.5, +19.5);      
       histogramRecTauDecayMode_ = dqmStore.book1D("recTauDecayMode", "recTauDecayMode", 21, -1.5, +19.5);  
     }

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

inline

Definition at line 991 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);
     typedef std::vector<pat::Tau> recLeptonCollection;
     edm::Handle<recLeptonCollection> recLeptons;
     evt.getByLabel(srcRec_, recLeptons);
     for ( CandidateView::const_iterator genLepton = genLeptons->begin();
           genLepton != genLeptons->end(); ++genLepton ) {
       if ( cutGen_ && !(*cutGen_)(*genLepton) ) continue;
       const reco::CompositePtrCandidate* genLepton_composite = dynamic_cast<const reco::CompositePtrCandidate*>(&(*genLepton));
       for ( 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_ ) {
           std::string genTauDecayMode_string  = JetMCTagUtils::genTauDecayMode(*genLepton_composite);       
           //std::cout << "--> genTauDecayMode: = " << genTauDecayMode_string << std::endl;
           int genTauDecayMode = -1;
           if      ( genTauDecayMode_string == "electron"        ) genTauDecayMode = reco::PFTauDecayMode::tauDecaysElectron;
           else if ( genTauDecayMode_string == "muon"            ) genTauDecayMode = reco::PFTauDecayMode::tauDecayMuon;
           else if ( genTauDecayMode_string == "oneProng0Pi0"    ) genTauDecayMode = reco::PFTauDecayMode::tauDecay1ChargedPion0PiZero;
           else if ( genTauDecayMode_string == "oneProng1Pi0"    ) genTauDecayMode = reco::PFTauDecayMode::tauDecay1ChargedPion1PiZero;
           else if ( genTauDecayMode_string == "oneProng2Pi0"    ) genTauDecayMode = reco::PFTauDecayMode::tauDecay1ChargedPion2PiZero;
           else if ( genTauDecayMode_string == "oneProngOther"   ) genTauDecayMode = reco::PFTauDecayMode::tauDecay1ChargedPion3PiZero;
           else if ( genTauDecayMode_string == "threeProng0Pi0"  ) genTauDecayMode = reco::PFTauDecayMode::tauDecay3ChargedPion0PiZero;
           else if ( genTauDecayMode_string == "threeProng1Pi0"  ) genTauDecayMode = reco::PFTauDecayMode::tauDecay3ChargedPion1PiZero;
           else if ( genTauDecayMode_string == "threeProngOther" ) genTauDecayMode = reco::PFTauDecayMode::tauDecay3ChargedPion2PiZero;
           else if ( genTauDecayMode_string == "rare"            ) genTauDecayMode = reco::PFTauDecayMode::tauDecayOther;
           histogramGenTauDecayMode_->Fill(genTauDecayMode, evtWeight);
           int recTauDecayMode = recLepton->decayMode();
           histogramRecTauDecayMode_->Fill(recTauDecayMode, evtWeight);
         }
       } 
     }
     //edm::Handle<reco::GenParticleCollection> genParticles;
     //evt.getByLabel("genParticles", genParticles);     
     //for ( reco::GenParticleCollection::const_iterator genParticle = genParticles->begin();
     //      genParticle != genParticles->end(); ++genParticle ) {
     //  if ( TMath::Abs(genParticle->pdgId()) == 15 ) {
     //    for ( recLeptonCollection::const_iterator recLepton = recLeptons->begin();
     //          recLepton != recLeptons->end(); ++recLepton ) {
     //      if ( cutRec_ && !(*cutRec_)(*recLepton) ) continue;         
     //      double dR = deltaR(genParticle->p4(), recLepton->p4());
     //      if ( dR < dRmatch_ ) {
     //        std::string genTauDecayMode_string = getGenTauDecayMode(&(*genParticle));
     //        std::cout << "--> genTauDecayMode: = " << genTauDecayMode_string << std::endl;
     //      }
     //    }
     //  }
     //}
       }
     }