#include <PatTauAnalyzer.h>

Inheritance diagram for PatTauAnalyzer:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &)

void	beginJob ()

void	endJob ()

	PatTauAnalyzer (const edm::ParameterSet &)

	~PatTauAnalyzer ()

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Private Attributes
std::string	discrByIso_

std::string	discrByLeadTrack_

std::string	discrByTaNC_

TH1 *	hGenTauEnergy_

TH1 *	hGenTauEta_

TH1 *	hGenTauPhi_

TH1 *	hGenTauPt_

TH1 *	hNumTauJets_

TH1 *	hTauDiscrAgainstElectrons_

TH1 *	hTauDiscrAgainstMuons_

TH1 *	hTauDiscrByIso_

TH1 *	hTauDiscrByTaNC_

TH1 *	hTauJetEnergy_

TH1 *	hTauJetEnergyIsoPassed_

TH1 *	hTauJetEnergyTaNCpassed_

TH1 *	hTauJetEta_

TH1 *	hTauJetEtaIsoPassed_

TH1 *	hTauJetEtaTaNCpassed_

TH1 *	hTauJetPhi_

TH1 *	hTauJetPhiIsoPassed_

TH1 *	hTauJetPhiTaNCpassed_

TH1 *	hTauJetPt_

TH1 *	hTauJetPtIsoPassed_

TH1 *	hTauJetPtTaNCpassed_

TH1 *	hTauLeadTrackPt_

TH1 *	hTauNumIsoConeTracks_

TH1 *	hTauNumSigConeTracks_

bool	requireGenTauMatch_

edm::InputTag	src_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 14 of file PatTauAnalyzer.h.

Constructor & Destructor Documentation

PatTauAnalyzer::PatTauAnalyzer ( const edm::ParameterSet & cfg )

explicit

Definition at line 25 of file PatTauAnalyzer.cc.

References discrByIso_, discrByLeadTrack_, discrByTaNC_, edm::ParameterSet::getParameter(), requireGenTauMatch_, and src_.

 {
   //std::cout << "<PatTauAnalyzer::PatTauAnalyzer>:" << std::endl;
 
 //--- read name of pat::Tau collection
   src_ = cfg.getParameter<edm::InputTag>("src");
   //std::cout << " src = " << src_ << std::endl;
 
 //--- fill histograms for all tau-jet candidates or for "real" taus only ?
   requireGenTauMatch_ = cfg.getParameter<bool>("requireGenTauMatch");
   //std::cout << " requireGenTauMatch = " << requireGenTauMatch_ << std::endl;
 
 //--- read names of tau id. discriminators
   discrByLeadTrack_ = cfg.getParameter<std::string>("discrByLeadTrack");
   //std::cout << " discrByLeadTrack = " << discrByLeadTrack_ << std::endl;
 
   discrByIso_ = cfg.getParameter<std::string>("discrByIso");
   //std::cout << " discrByIso = " << discrByIso_ << std::endl;
 
   discrByTaNC_ = cfg.getParameter<std::string>("discrByTaNC");
   //std::cout << " discrByTaNC = " << discrByTaNC_ << std::endl;
 }

PatTauAnalyzer::~PatTauAnalyzer ( )

Definition at line 48 of file PatTauAnalyzer.cc.

 {
   //std::cout << "<PatTauAnalyzer::~PatTauAnalyzer>:" << std::endl;
 
 //--- clean-up memory;
 //    delete all histograms
 /*
   deletion of histograms taken care of by TFileService;
   do not delete them here (if the histograms are deleted here,
   they will not appear in the ROOT file written by TFileService)
 
   delete hGenTauEnergy_;
   delete hGenTauPt_;
   delete hGenTauEta_;
   delete hGenTauPhi_;
   delete hTauJetEnergy_;
   delete hTauJetPt_;
   delete hTauJetEta_;
   delete hTauJetPhi_;
   delete hNumTauJets_;
   delete hTauLeadTrackPt_;
   delete hTauNumSigConeTracks_;
   delete hTauNumIsoConeTracks_;
   delete hTauDiscrByIso_;
   delete hTauDiscrByTaNC_;
   delete hTauDiscrAgainstElectrons_;
   delete hTauDiscrAgainstMuons_;
   delete hTauJetEnergyIsoPassed_;
   delete hTauJetPtIsoPassed_;
   delete hTauJetEtaIsoPassed_;
   delete hTauJetPhiIsoPassed_;
   delete hTauJetEnergyTaNCpassed_;
   delete hTauJetPtTaNCpassed_;
   delete hTauJetEtaTaNCpassed_;
   delete hTauJetPhiTaNCpassed_;
  */
 }

Member Function Documentation

void PatTauAnalyzer::analyze	(	const edm::Event &	evt,
		const edm::EventSetup &	es
	)

virtual

Implements edm::EDAnalyzer.

Definition at line 152 of file PatTauAnalyzer.cc.

References discrByIso_, discrByLeadTrack_, discrByTaNC_, reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), edm::Event::getByLabel(), getGenTau(), hGenTauEnergy_, hGenTauEta_, hGenTauPhi_, hGenTauPt_, hNumTauJets_, hTauDiscrAgainstMuons_, hTauDiscrByIso_, hTauDiscrByTaNC_, hTauJetEnergy_, hTauJetEnergyIsoPassed_, hTauJetEnergyTaNCpassed_, hTauJetEtaIsoPassed_, hTauJetEtaTaNCpassed_, hTauJetPhiIsoPassed_, hTauJetPhiTaNCpassed_, hTauJetPt_, hTauJetPtIsoPassed_, hTauJetPtTaNCpassed_, hTauNumIsoConeTracks_, hTauNumSigConeTracks_, tauProducer_cfi::patTaus, reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), requireGenTauMatch_, and src_.

 {  
   //std::cout << "<PatTauAnalyzer::analyze>:" << std::endl; 
 
   edm::Handle<pat::TauCollection> patTaus;
   evt.getByLabel(src_, patTaus);
 
   hNumTauJets_->Fill(patTaus->size());
 
   for ( pat::TauCollection::const_iterator patTau = patTaus->begin(); 
     patTau != patTaus->end(); ++patTau ) {
 
 //--- skip fake taus in case configuration parameters set to do so...
     const reco::GenParticle* genTau = getGenTau(*patTau);
     if ( requireGenTauMatch_ && !genTau ) continue;
 
 //--- fill generator level histograms    
     if ( genTau ) {
       hGenTauEnergy_->Fill(genTau->energy());
       hGenTauPt_->Fill(genTau->pt());
       hGenTauEta_->Fill(genTau->eta());
       hGenTauPhi_->Fill(genTau->phi());
     }
 
 //--- fill reconstruction level histograms
 //    for Pt of highest Pt track within signal cone tau-jet...
     hTauJetEnergy_->Fill(patTau->energy());
     hTauJetPt_->Fill(patTau->pt());
 // 
 // TO-DO: 
 //  1.) fill histograms 
 //       hTauJetEta_ and hTauJetPhi_
 //      with the pseudo-rapidity and azimuthal angle
 //      of the tau-jet candidate respectively
 //  hTauJetEta_->...
 //  hTauJetPhi_->...
 //
 //  2.) fill histogram
 //       hTauLeadTrackPt_
 //      with the transverse momentum of the highest Pt ("leading") track within the tau-jet
 //       
 // NOTE: 
 //  1.) please have a look at
 //       http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/CMSSW/DataFormats/Candidate/interface/Particle.h?revision=1.28&view=markup
 //      to find the methods for accessing eta and phi of the tau-jet
 //
 //  2.) please have a look at
 //       http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/CMSSW/DataFormats/PatCandidates/interface/Tau.h?revision=1.25&view=markup
 //      to find the method for accessing the leading track
 //
 //  3.) the method pat::Tau::leadTrack returns a reference (reco::TrackRef) to a reco::Track object
 //      this reference can be null (in case no high Pt track has been reconstructed within the tau-jet),
 //      so a check if the leadTrack exists is needed before dereferencing the reco::TrackRef via operator->
 //
 //  if ( patTau->leadTrack().isAvailable() ) hTauLeadTrackPt_->Fill(patTau->leadTrack()->pt());
   
 //... for total number of tracks within signal/isolation cones
     hTauNumSigConeTracks_->Fill(patTau->signalTracks().size());
     hTauNumIsoConeTracks_->Fill(patTau->isolationTracks().size());
 
 //... for values of tau id. discriminators based on track isolation cut/
 //    neural network-based tau id.
 //    (combine with requirement of at least one "leading" track of Pt > 5. GeV 
 //     within the signal cone of the tau-jet)
     float discrByIso = ( patTau->tauID(discrByLeadTrack_.data()) > 0.5 ) ? patTau->tauID(discrByIso_.data()) : 0.;
     hTauDiscrByIso_->Fill(discrByIso);
     float discrByTaNC = ( patTau->tauID(discrByLeadTrack_.data()) > 0.5 ) ? patTau->tauID(discrByTaNC_.data()) : 0.;
     hTauDiscrByTaNC_->Fill(discrByTaNC);
 
 //... for values of tau id. discriminators against (unidentified) electrons and muons
 //
 // TO-DO: fill histogram
 //         hTauDiscrAgainstElectrons_
 //        with the value of the discriminatorAgainstElectrons
 //
 // NOTE:
 //  1.) please have a look at
 //       http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/CMSSW/DataFormats/PatCandidates/interface/Tau.h?revision=1.25&view=markup
 //      to find the method for accessing the tau id. information
 //  
 //  2.) please have a look at
 //       http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/CMSSW/PhysicsTools/PatAlgos/python/tools/tauTools.py?revision=1.18&view=markup
 //      and convince yourself that the string "againstElectron" needs to be passed as argument 
 //      of the pat::Tau::tauID method
 //
 //  hTauDiscrAgainstElectrons_->Fill...
     hTauDiscrAgainstMuons_->Fill(patTau->tauID("againstMuon"));
 
 //... for Energy, Pt, Eta, Phi of tau-jets passing the discriminatorByIsolation selection
     if ( discrByIso > 0.5 ) {
       hTauJetEnergyIsoPassed_->Fill(patTau->energy());
       hTauJetPtIsoPassed_->Fill(patTau->pt());
       hTauJetEtaIsoPassed_->Fill(patTau->eta());
       hTauJetPhiIsoPassed_->Fill(patTau->phi());
     }
 
 //... for Energy, Pt, Eta, Phi of tau-jets passing the discriminatorByTaNC ("Tau Neural Classifier") selection
     if ( discrByTaNC > 0.5 ) {
       hTauJetEnergyTaNCpassed_->Fill(patTau->energy());
       hTauJetPtTaNCpassed_->Fill(patTau->pt());
       hTauJetEtaTaNCpassed_->Fill(patTau->eta());
       hTauJetPhiTaNCpassed_->Fill(patTau->phi());
     }
   }
 }

void PatTauAnalyzer::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 86 of file PatTauAnalyzer.cc.

References hGenTauEnergy_, hGenTauEta_, hGenTauPhi_, hGenTauPt_, hNumTauJets_, hTauDiscrAgainstElectrons_, hTauDiscrAgainstMuons_, hTauDiscrByIso_, hTauDiscrByTaNC_, hTauJetEnergy_, hTauJetEnergyIsoPassed_, hTauJetEnergyTaNCpassed_, hTauJetEtaIsoPassed_, hTauJetEtaTaNCpassed_, hTauJetPhiIsoPassed_, hTauJetPhiTaNCpassed_, hTauJetPt_, hTauJetPtIsoPassed_, hTauJetPtTaNCpassed_, hTauLeadTrackPt_, hTauNumIsoConeTracks_, hTauNumSigConeTracks_, TFileDirectory::make(), and Pi.

 {
 //--- retrieve handle to auxiliary service
 //    used for storing histograms into ROOT file
   edm::Service<TFileService> fs;
 
 //--- book generator level histograms
   hGenTauEnergy_ = fs->make<TH1F>("GenTauEnergy", "GenTauEnergy", 30, 0., 150.);
   hGenTauPt_ = fs->make<TH1F>("GenTauPt", "GenTauPt", 30, 0., 150.);
   hGenTauEta_ = fs->make<TH1F>("GenTauEta", "GenTauEta", 24, -3., +3.);
   hGenTauPhi_ = fs->make<TH1F>("GenTauPhi", "GenTauPhi", 18, -TMath::Pi(), +TMath::Pi());
   
 //--- book reconstruction level histograms
 //    for tau-jet Energy, Pt, Eta, Phi
   hTauJetEnergy_ = fs->make<TH1F>("TauJetEnergy", "TauJetEnergy", 30, 0., 150.);
   hTauJetPt_ = fs->make<TH1F>("TauJetPt", "TauJetPt", 30, 0., 150.);
 // 
 // TO-DO: add histograms for eta and phi of the tau-jet candidate
 //      
 // NOTE:
 //  1.) please use 
 //       "TauJetEta" and "TauJetPhi" 
 //      for the names of the histograms and choose the exact same binning 
 //      as is used for the histograms 
 //       "TauJetEtaIsoPassed" and "TauJetPhiIsoPassed" 
 //      below
 //
 //  2.) please check the histograms
 //       hTauJetEta_ and hTauJetPt_
 //      have already been defined in PatTauAnalyzer.h
 // 
 //hTauJetEta_ =...
 //hTauJetPt_ =...
 
 //... for number of tau-jet candidates
   hNumTauJets_ = fs->make<TH1F>("NumTauJets", "NumTauJets", 10, -0.5, 9.5);
 
 //... for Pt of highest Pt track within signal cone tau-jet...
   hTauLeadTrackPt_ = fs->make<TH1F>("TauLeadTrackPt", "TauLeadTrackPt", 40, 0., 100.);
   
 //... for total number of tracks within signal/isolation cones
   hTauNumSigConeTracks_ = fs->make<TH1F>("TauNumSigConeTracks", "TauNumSigConeTracks", 10, -0.5,  9.5);
   hTauNumIsoConeTracks_ = fs->make<TH1F>("TauNumIsoConeTracks", "TauNumIsoConeTracks", 20, -0.5, 19.5);
 
 //... for values of tau id. discriminators based on track isolation cut/
 //    neural network-based tau id.
   hTauDiscrByIso_ = fs->make<TH1F>("TauDiscrByIso", "TauDiscrByIso", 103, -0.015, 1.015);
   hTauDiscrByTaNC_ = fs->make<TH1F>("TauDiscrByTaNC", "TauDiscrByTaNC", 103, -0.015, 1.015);
   
 //... for values of tau id. discriminators against (unidentified) electrons and muons
   hTauDiscrAgainstElectrons_ = fs->make<TH1F>("TauDiscrAgainstElectrons", "TauDiscrAgainstElectrons", 103, -0.015, 1.015);
   hTauDiscrAgainstMuons_ = fs->make<TH1F>("TauDiscrAgainstMuons", "TauDiscrAgainstMuons", 103, -0.015, 1.015);
 
 //... for Energy, Pt, Eta, Phi of tau-jets passing the discriminatorByIsolation selection
   hTauJetEnergyIsoPassed_ = fs->make<TH1F>("TauJetEnergyIsoPassed", "TauJetEnergyIsoPassed", 30, 0., 150.);
   hTauJetPtIsoPassed_ = fs->make<TH1F>("TauJetPtIsoPassed", "TauJetPtIsoPassed", 30, 0., 150.);
   hTauJetEtaIsoPassed_ = fs->make<TH1F>("TauJetEtaIsoPassed", "TauJetEtaIsoPassed", 24, -3., +3.);
   hTauJetPhiIsoPassed_ = fs->make<TH1F>("TauJetPhiIsoPassed", "TauJetPhiIsoPassed", 18, -TMath::Pi(), +TMath::Pi());
 
 //... for Energy, Pt, Eta, Phi of tau-jets passing the discriminatorByTaNC ("Tau Neural Classifier") selection
   hTauJetEnergyTaNCpassed_ = fs->make<TH1F>("TauJetEnergyTaNCpassed", "TauJetEnergyTaNCpassed", 30, 0., 150.);
   hTauJetPtTaNCpassed_ = fs->make<TH1F>("TauJetPtTaNCpassed", "TauJetPtTaNCpassed", 30, 0., 150.);
   hTauJetEtaTaNCpassed_ = fs->make<TH1F>("TauJetEtaTaNCpassed", "TauJetEtaTaNCpassed", 24, -3., +3.);
   hTauJetPhiTaNCpassed_ = fs->make<TH1F>("TauJetPhiTaNCpassed", "TauJetPhiTaNCpassed", 18, -TMath::Pi(), +TMath::Pi());
 }

void PatTauAnalyzer::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 258 of file PatTauAnalyzer.cc.

 {
 //--- nothing to be done yet...
 }

Member Data Documentation

std::string PatTauAnalyzer::discrByIso_

private

Definition at line 32 of file PatTauAnalyzer.h.

Referenced by analyze(), and PatTauAnalyzer().

std::string PatTauAnalyzer::discrByLeadTrack_

private

Definition at line 31 of file PatTauAnalyzer.h.

Referenced by analyze(), and PatTauAnalyzer().

std::string PatTauAnalyzer::discrByTaNC_

private

Definition at line 33 of file PatTauAnalyzer.h.

Referenced by analyze(), and PatTauAnalyzer().

TH1* PatTauAnalyzer::hGenTauEnergy_

private

Definition at line 36 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hGenTauEta_

private

Definition at line 38 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hGenTauPhi_

private

Definition at line 39 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hGenTauPt_

private

Definition at line 37 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hNumTauJets_

private

Definition at line 47 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hTauDiscrAgainstElectrons_

private

Definition at line 56 of file PatTauAnalyzer.h.

Referenced by beginJob().

TH1* PatTauAnalyzer::hTauDiscrAgainstMuons_

private

Definition at line 57 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hTauDiscrByIso_

private

Definition at line 54 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hTauDiscrByTaNC_

private

Definition at line 55 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hTauJetEnergy_

private

Definition at line 42 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hTauJetEnergyIsoPassed_

private

Definition at line 59 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hTauJetEnergyTaNCpassed_

private

Definition at line 64 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hTauJetEta_

private

Definition at line 44 of file PatTauAnalyzer.h.

TH1* PatTauAnalyzer::hTauJetEtaIsoPassed_

private

Definition at line 61 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hTauJetEtaTaNCpassed_

private

Definition at line 66 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hTauJetPhi_

private

Definition at line 45 of file PatTauAnalyzer.h.

TH1* PatTauAnalyzer::hTauJetPhiIsoPassed_

private

Definition at line 62 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hTauJetPhiTaNCpassed_

private

Definition at line 67 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hTauJetPt_

private

Definition at line 43 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hTauJetPtIsoPassed_

private

Definition at line 60 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hTauJetPtTaNCpassed_

private

Definition at line 65 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hTauLeadTrackPt_

private

Definition at line 49 of file PatTauAnalyzer.h.

Referenced by beginJob().

TH1* PatTauAnalyzer::hTauNumIsoConeTracks_

private

Definition at line 52 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1* PatTauAnalyzer::hTauNumSigConeTracks_

private

Definition at line 51 of file PatTauAnalyzer.h.

Referenced by analyze(), and beginJob().

bool PatTauAnalyzer::requireGenTauMatch_

private

Definition at line 29 of file PatTauAnalyzer.h.

Referenced by analyze(), and PatTauAnalyzer().

edm::InputTag PatTauAnalyzer::src_

private

Definition at line 27 of file PatTauAnalyzer.h.

Referenced by analyze(), and PatTauAnalyzer().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation