TTbar_GenJetAnalyzer.cc

Go to the documentation of this file.

#include "Validation/EventGenerator/interface/TTbar_GenJetAnalyzer.h"
#include "Validation/EventGenerator/interface/DQMHelper.h"

TTbar_GenJetAnalyzer::TTbar_GenJetAnalyzer(const edm::ParameterSet& iConfig)
    : jets_(iConfig.getParameter<edm::InputTag>("jets")),
      genEventInfoProductTag_(iConfig.getParameter<edm::InputTag>("genEventInfoProductTag")) {
  genEventInfoProductTagToken_ = consumes<GenEventInfoProduct>(genEventInfoProductTag_);
  jetsToken_ = consumes<std::vector<reco::GenJet> >(jets_);
}

TTbar_GenJetAnalyzer::~TTbar_GenJetAnalyzer() {}

void TTbar_GenJetAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
  using namespace edm;

  // --- the MC weights ---
  Handle<GenEventInfoProduct> evt_info;
  iEvent.getByToken(genEventInfoProductTagToken_, evt_info);
  if (!evt_info.isValid())
    return;
  weight = evt_info->weight();

  // Gather information in the GenJet collection
  edm::Handle<std::vector<reco::GenJet> > jets;
  iEvent.getByToken(jetsToken_, jets);

  if (!jets.isValid())
    return;
  // loop Jet collection and fill histograms
  int njets = 0;
  for (std::vector<reco::GenJet>::const_iterator jet_it = jets->begin(); jet_it != jets->end(); ++jet_it) {
    ++njets;

    hists_["jetPtAll"]->Fill(jet_it->pt(), weight);
    hists_["jetEtaAll"]->Fill(jet_it->eta(), weight);

    if (njets == 1) {
      hists_["jetPt1"]->Fill(jet_it->pt(), weight);
      hists_["jetEta1"]->Fill(jet_it->eta(), weight);
    }
    if (njets == 2) {
      hists_["jetPt2"]->Fill(jet_it->pt(), weight);
      hists_["jetEta2"]->Fill(jet_it->eta(), weight);
    }
    if (njets == 3) {
      hists_["jetPt3"]->Fill(jet_it->pt(), weight);
      hists_["jetEta3"]->Fill(jet_it->eta(), weight);
    }
    if (njets == 4) {
      hists_["jetPt4"]->Fill(jet_it->pt(), weight);
      hists_["jetEta4"]->Fill(jet_it->eta(), weight);
    }
  }
}

void TTbar_GenJetAnalyzer::bookHistograms(DQMStore::IBooker& i, edm::Run const&, edm::EventSetup const&) {
  DQMHelper dqm(&i);
  i.setCurrentFolder("Generator/TTbar");
  hists_["jetPtAll"] =
      dqm.book1dHisto("TTbar_jetPtAll", "pt", 1000, 0., 1000., "P_{t}^{All-Jets} (GeV)", "Number of Events");
  hists_["jetPt1"] =
      dqm.book1dHisto("TTbar_jetPt1", "pt", 1000, 0., 1000., "P_{t}^{1st-Jet} (GeV)", "Number of Events");
  hists_["jetPt2"] =
      dqm.book1dHisto("TTbar_jetPt2", "pt", 1000, 0., 1000., "P_{t}^{2nd-Jet} (GeV)", "Number of Events");
  hists_["jetPt3"] =
      dqm.book1dHisto("TTbar_jetPt3", "pt", 1000, 0., 1000., "P_{t}^{3rd-Jet} (GeV)", "Number of Events");
  hists_["jetPt4"] =
      dqm.book1dHisto("TTbar_jetPt4", "pt", 1000, 0., 1000., "P_{t}^{4th-Jet} (GeV)", "Number of Events");

  hists_["jetEtaAll"] = dqm.book1dHisto("TTbar_jetEtaAll", "eta", 100, -5., 5., "#eta^{All-Jets}", "Number of Events");
  hists_["jetEta1"] = dqm.book1dHisto("TTbar_jetEta1", "eta", 100, -5., 5., "#eta^{1st-Jet}", "Number of Events");
  hists_["jetEta2"] = dqm.book1dHisto("TTbar_jetEta2", "eta", 100, -5., 5., "#eta^{2nd-Jet}", "Number of Events");
  hists_["jetEta3"] = dqm.book1dHisto("TTbar_jetEta3", "eta", 100, -5., 5., "#eta^{3rd-Jet}", "Number of Events");
  hists_["jetEta4"] = dqm.book1dHisto("TTbar_jetEta4", "eta", 100, -5., 5., "#eta^{4th-Jet}", "Number of Events");
}