TopElecAnalyzer.cc

Go to the documentation of this file.

#include "TopQuarkAnalysis/Examples/plugins/TopElecAnalyzer.h"
 
TopElecAnalyzer::TopElecAnalyzer(const edm::ParameterSet& cfg)
    : inputToken_(consumes<std::vector<pat::Electron> >(cfg.getParameter<edm::InputTag>("input"))),
      verbose_(cfg.getParameter<bool>("verbose")) {
  edm::Service<TFileService> fs;
 
  mult_ = fs->make<TH1F>("mult", "multiplicity (electrons)", 10, 0, 10);
  en_ = fs->make<TH1F>("en", "energy (electrons)", 60, 0., 300.);
  pt_ = fs->make<TH1F>("pt", "pt (electrons)", 60, 0., 300.);
  eta_ = fs->make<TH1F>("eta", "eta (electrons)", 30, -3., 3.);
  phi_ = fs->make<TH1F>("phi", "phi (electrons)", 40, -4., 4.);
}
 
TopElecAnalyzer::~TopElecAnalyzer() {}
 
void TopElecAnalyzer::analyze(const edm::Event& evt, const edm::EventSetup& setup) {
  edm::Handle<std::vector<pat::Electron> > elecs;
  evt.getByToken(inputToken_, elecs);
 
  // fill histograms
 
  mult_->Fill(elecs->size());
  for (std::vector<pat::Electron>::const_iterator elec = elecs->begin(); elec != elecs->end(); ++elec) {
    en_->Fill(elec->energy());
    pt_->Fill(elec->pt());
    eta_->Fill(elec->eta());
    phi_->Fill(elec->phi());
  }
 
  // produce printout if desired
 
  if (elecs->empty() || !verbose_)
    return;
 
  unsigned i = 0;
 
  std::cout << "=======================================================" << std::endl;
  std::cout << std::setw(5) << "ele :" << std::setw(13) << "et :" << std::setw(13) << "eta :" << std::setw(13)
            << "phi :" << std::setw(11) << "relIso" << std::endl;
  std::cout << "-------------------------------------------------------" << std::endl;
  for (std::vector<pat::Electron>::const_iterator elec = elecs->begin(); elec != elecs->end(); ++elec) {
    std::cout << std::setw(3) << i << " : " << std::setw(10) << elec->pt() << " : " << std::setw(10) << elec->eta()
              << " : " << std::setw(10) << elec->phi() << " : " << std::setw(10)
              << (elec->dr03TkSumPt() + elec->dr03EcalRecHitSumEt() + elec->dr03HcalTowerSumEt()) / elec->et()
              << std::endl;
    i++;
  }
  std::cout << "=======================================================" << std::endl;
}
 
void TopElecAnalyzer::beginJob() {}
 
void TopElecAnalyzer::endJob() {}