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/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::EDAnalyzer {

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")))
{
}

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);