PatZjetsJetAnalyzer.cc

Go to the documentation of this file.

#include <map>
#include <string>

#include "TH1.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "CommonTools/UtilAlgos/interface/TFileService.h"

#include "DataFormats/PatCandidates/interface/Jet.h"

class PatZjetsJetAnalyzer : public edm::one::EDAnalyzer<edm::one::SharedResources> {
public:
  explicit PatZjetsJetAnalyzer(const edm::ParameterSet&);
  ~PatZjetsJetAnalyzer() override;

private:
  void beginJob() override;
  void analyze(const edm::Event&, const edm::EventSetup&) override;
  void endJob() override;

  // simple map to contain all histograms;
  // histograms are booked in the beginJob()
  // method
  std::map<std::string, TH1F*> histContainer_;

  // input tags
  edm::EDGetTokenT<edm::View<pat::Jet> > srcToken_;
};

PatZjetsJetAnalyzer::PatZjetsJetAnalyzer(const edm::ParameterSet& iConfig)
    : histContainer_(), srcToken_(consumes<edm::View<pat::Jet> >(iConfig.getUntrackedParameter<edm::InputTag>("src"))) {
  usesResource(TFileService::kSharedResource);
}

PatZjetsJetAnalyzer::~PatZjetsJetAnalyzer() {}

void PatZjetsJetAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
  // get electron collection
  edm::Handle<edm::View<pat::Jet> > jets;
  iEvent.getByToken(srcToken_, jets);

  // loop jets
  for (edm::View<pat::Jet>::const_iterator ijet = jets->begin(); ijet != jets->end(); ++ijet) {
    // fill simple histograms
    pat::Jet jet = ijet->correctedJet("had", "uds");
    histContainer_["pt"]->Fill(jet.pt());
    histContainer_["eta"]->Fill(jet.eta());
    histContainer_["phi"]->Fill(jet.phi());
    histContainer_["emf"]->Fill(jet.emEnergyFraction());
    for (unsigned int i = 0; i < jet.getCaloConstituents().size(); ++i) {
      histContainer_["dEta"]->Fill(jet.getCaloConstituent(i)->eta() - jet.eta());
    }
  }
}

void PatZjetsJetAnalyzer::beginJob() {
  // register to the TFileService
  edm::Service<TFileService> fs;

  // book histograms:
  histContainer_["pt"] = fs->make<TH1F>("pt", "pt", 150, 0., 150.);
  histContainer_["eta"] = fs->make<TH1F>("eta", "eta", 50, 0., 5.);
  histContainer_["phi"] = fs->make<TH1F>("phi", "phi", 60, 3.14, 3.14);
  histContainer_["emf"] = fs->make<TH1F>("emf", "emf", 40, 0., 1.);
  histContainer_["dEta"] = fs->make<TH1F>("dEta", "dEta", 40, 0., 1.);
}

void PatZjetsJetAnalyzer::endJob() {}

#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(PatZjetsJetAnalyzer);