TopTauAnalyzer.cc

Go to the documentation of this file.

#include "TopQuarkAnalysis/Examples/plugins/TopTauAnalyzer.h"


TopTauAnalyzer::TopTauAnalyzer(const edm::ParameterSet& cfg):
  inputToken_(consumes<std::vector<pat::Tau> >(cfg.getParameter<edm::InputTag>("input")))
{
  edm::Service<TFileService> fs;

  mult_ = fs->make<TH1F>("mult", "multiplicity (taus)", 30,  0 ,   30);
  en_   = fs->make<TH1F>("en"  , "energy (taus)",       60,  0., 300.);
  pt_   = fs->make<TH1F>("pt"  , "pt (taus}",           60,  0., 300.);
  eta_  = fs->make<TH1F>("eta" , "eta (taus)",          30, -3.,   3.);
  phi_  = fs->make<TH1F>("phi" , "phi (taus)",          40, -4.,   4.);
}

TopTauAnalyzer::~TopTauAnalyzer()
{
}

void
TopTauAnalyzer::analyze(const edm::Event& evt, const edm::EventSetup& setup)
{
  edm::Handle<std::vector<pat::Tau> > taus;
  evt.getByToken(inputToken_, taus);

  // fill histograms

  mult_->Fill( taus->size() );
  for(std::vector<pat::Tau>::const_iterator tau=taus->begin(); tau!=taus->end(); ++tau){
    en_ ->Fill( tau->energy() );
    pt_ ->Fill( tau->pt()     );
    eta_->Fill( tau->eta()    );
    phi_->Fill( tau->phi()    );
  }
}

void TopTauAnalyzer::beginJob()
{
}

void TopTauAnalyzer::endJob()
{
}