PFMCTruthTreeProducer Class Reference

#include <PFMCTruthTreeProducer.h>

Inheritance diagram for PFMCTruthTreeProducer:

Public Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &iSetup) override
void	beginJob () override
void	endJob () override
	PFMCTruthTreeProducer (edm::ParameterSet const &cfg)
	~PFMCTruthTreeProducer () override
Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () const final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Attributes
float	cemfJet_
float	chfJet_
int	cmultiJet_
float	dR_
float	etaGen_
float	etaJet_
TFile *	file_
edm::EDGetTokenT< GenEventInfoProduct >	gen_
edm::EDGetTokenT< reco::GenJetCollection >	genjets_
std::string	histogramFile_
edm::EDGetTokenT< reco::PFJetCollection >	jets_
TTree *	mcTruthTree_
float	nemfJet_
float	nhfJet_
int	nmultiJet_
float	phiGen_
float	phiJet_
float	ptGen_
float	ptHat_
float	ptJet_
int	rank_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 13 of file PFMCTruthTreeProducer.h.

Constructor & Destructor Documentation

PFMCTruthTreeProducer()

PFMCTruthTreeProducer::PFMCTruthTreeProducer ( edm::ParameterSet const &  cfg )

explicit

Definition at line 26 of file PFMCTruthTreeProducer.cc.

References looper::cfg, ProducerED_cfi::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                       {
  jets_ = consumes<PFJetCollection>(edm::InputTag(cfg.getParameter<std::string>("jets")));
  genjets_ = consumes<GenJetCollection>(edm::InputTag(cfg.getParameter<std::string>("genjets")));
  gen_ = consumes<GenEventInfoProduct>(edm::InputTag("generator"));
  histogramFile_ = cfg.getParameter<std::string>("histogramFile");
}

~PFMCTruthTreeProducer()

PFMCTruthTreeProducer::~PFMCTruthTreeProducer ( )

override

Definition at line 104 of file PFMCTruthTreeProducer.cc.

                                              {
  delete file_;
  delete mcTruthTree_;
}

Member Function Documentation

analyze()

void PFMCTruthTreeProducer::analyze ( edm::Event const &  e, edm::EventSetup const &  iSetup  )

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 62 of file PFMCTruthTreeProducer.cc.

References GenEventInfoProduct::binningValues(), PDWG_EXODelayedJetMET_cff::jets, njet, CosmicsPD_Skims::radius, and findQualityFiles::rr.

                                                                                      {
  edm::Handle<GenJetCollection> genjets;
  edm::Handle<PFJetCollection> jets;
  edm::Handle<GenEventInfoProduct> hEventInfo;
  PFJetCollection::const_iterator i_jet, i_matched;
  GenJetCollection::const_iterator i_genjet;
  event.getByToken(genjets_, genjets);
  event.getByToken(jets_, jets);
  event.getByToken(gen_, hEventInfo);
  ptHat_ = hEventInfo->binningValues()[0];
  float rr;
  int njet(0);
  if (!jets->empty() && !genjets->empty()) {
    for (i_genjet = genjets->begin(); i_genjet != genjets->end(); i_genjet++) {
      float rmin(99);
      for (i_jet = jets->begin(); i_jet != jets->end(); i_jet++) {
        rr = radius(i_genjet, i_jet);
        if (rr < rmin) {
          rmin = rr;
          i_matched = i_jet;
        }
      }
      ptGen_ = i_genjet->pt();
      etaGen_ = i_genjet->eta();
      phiGen_ = i_genjet->phi();
      ptJet_ = i_matched->pt();
      etaJet_ = i_matched->eta();
      phiJet_ = i_matched->phi();
      chfJet_ = i_matched->chargedHadronEnergyFraction();
      nhfJet_ = i_matched->neutralHadronEnergyFraction();
      cemfJet_ = i_matched->chargedEmEnergyFraction();
      nemfJet_ = i_matched->neutralEmEnergyFraction();
      cmultiJet_ = i_matched->chargedMultiplicity();
      nmultiJet_ = i_matched->neutralMultiplicity();
      dR_ = rmin;
      rank_ = njet;
      mcTruthTree_->Fill();
      njet++;
    }
  }
}

beginJob()

void PFMCTruthTreeProducer::beginJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 33 of file PFMCTruthTreeProducer.cc.

                                     {
  file_ = new TFile(histogramFile_.c_str(), "RECREATE");
  mcTruthTree_ = new TTree("mcTruthTree", "mcTruthTree");

  mcTruthTree_->Branch("ptJet", &ptJet_, "ptJet_/F");
  mcTruthTree_->Branch("ptGen", &ptGen_, "ptGen_/F");
  mcTruthTree_->Branch("ptHat", &ptHat_, "ptHat_/F");
  mcTruthTree_->Branch("chfJet", &chfJet_, "chfJet_/F");
  mcTruthTree_->Branch("nhfJet", &nhfJet_, "nhfJet_/F");
  mcTruthTree_->Branch("cemfJet", &cemfJet_, "cemfJet_/F");
  mcTruthTree_->Branch("nemfJet", &nemfJet_, "nemfJet_/F");
  mcTruthTree_->Branch("cmultiJet", &cmultiJet_, "cmultiJet_/I");
  mcTruthTree_->Branch("nmultiJet", &nmultiJet_, "nmultiJet_/I");
  mcTruthTree_->Branch("etaJet", &etaJet_, "etaJet_/F");
  mcTruthTree_->Branch("etaGen", &etaGen_, "etaGen_/F");
  mcTruthTree_->Branch("phiJet", &phiJet_, "phiJet_/F");
  mcTruthTree_->Branch("phiGen", &phiGen_, "phiGen_/F");
  mcTruthTree_->Branch("dR", &dR_, "dR_/F");
  mcTruthTree_->Branch("rank", &rank_, "rank_/I");
}

endJob()

void PFMCTruthTreeProducer::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 54 of file PFMCTruthTreeProducer.cc.

                                   {
  if (file_ != nullptr) {
    file_->cd();
    mcTruthTree_->Write();
  }
  file_ = nullptr;
}

Member Data Documentation

cemfJet_

float PFMCTruthTreeProducer::cemfJet_

private

Definition at line 28 of file PFMCTruthTreeProducer.h.

chfJet_

float PFMCTruthTreeProducer::chfJet_

private

Definition at line 28 of file PFMCTruthTreeProducer.h.

cmultiJet_

int PFMCTruthTreeProducer::cmultiJet_

private

Definition at line 29 of file PFMCTruthTreeProducer.h.

dR_

float PFMCTruthTreeProducer::dR_

private

Definition at line 28 of file PFMCTruthTreeProducer.h.

etaGen_

float PFMCTruthTreeProducer::etaGen_

private

Definition at line 28 of file PFMCTruthTreeProducer.h.

etaJet_

float PFMCTruthTreeProducer::etaJet_

private

Definition at line 28 of file PFMCTruthTreeProducer.h.

file_

TFile* PFMCTruthTreeProducer::file_

private

Definition at line 26 of file PFMCTruthTreeProducer.h.

gen_

edm::EDGetTokenT<GenEventInfoProduct> PFMCTruthTreeProducer::gen_

private

Definition at line 25 of file PFMCTruthTreeProducer.h.

genjets_

edm::EDGetTokenT<reco::GenJetCollection> PFMCTruthTreeProducer::genjets_

private

Definition at line 24 of file PFMCTruthTreeProducer.h.

histogramFile_

std::string PFMCTruthTreeProducer::histogramFile_

private

Definition at line 22 of file PFMCTruthTreeProducer.h.

jets_

edm::EDGetTokenT<reco::PFJetCollection> PFMCTruthTreeProducer::jets_

private

Definition at line 23 of file PFMCTruthTreeProducer.h.

mcTruthTree_

TTree* PFMCTruthTreeProducer::mcTruthTree_

private

Definition at line 27 of file PFMCTruthTreeProducer.h.

nemfJet_

float PFMCTruthTreeProducer::nemfJet_

private

Definition at line 28 of file PFMCTruthTreeProducer.h.

nhfJet_

float PFMCTruthTreeProducer::nhfJet_

private

Definition at line 28 of file PFMCTruthTreeProducer.h.

nmultiJet_

int PFMCTruthTreeProducer::nmultiJet_

private

Definition at line 29 of file PFMCTruthTreeProducer.h.

phiGen_

float PFMCTruthTreeProducer::phiGen_

private

Definition at line 28 of file PFMCTruthTreeProducer.h.

phiJet_

float PFMCTruthTreeProducer::phiJet_

private

Definition at line 28 of file PFMCTruthTreeProducer.h.

ptGen_

float PFMCTruthTreeProducer::ptGen_

private

Definition at line 28 of file PFMCTruthTreeProducer.h.

ptHat_

float PFMCTruthTreeProducer::ptHat_

private

Definition at line 28 of file PFMCTruthTreeProducer.h.

ptJet_

float PFMCTruthTreeProducer::ptJet_

private

Definition at line 28 of file PFMCTruthTreeProducer.h.

rank_

int PFMCTruthTreeProducer::rank_

private

Definition at line 29 of file PFMCTruthTreeProducer.h.