CaloMCTruthTreeProducer Class Reference

#include <CaloMCTruthTreeProducer.h>

Inheritance diagram for CaloMCTruthTreeProducer:

Public Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &iSetup) override
void	beginJob () override
	CaloMCTruthTreeProducer (edm::ParameterSet const &cfg)
void	endJob () override
	~CaloMCTruthTreeProducer () override
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
	~EDAnalyzer () 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
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::vector< ModuleDescription const * > &modules, 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
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Private Attributes
float	dR_
float	emfJet_
float	etaGen_
float	etaJet_
TFile *	file_
edm::EDGetTokenT< GenEventInfoProduct >	gen_
edm::EDGetTokenT< reco::GenJetCollection >	genjets_
std::string	histogramFile_
edm::EDGetTokenT< reco::CaloJetCollection >	jets_
TTree *	mcTruthTree_
float	phiGen_
float	phiJet_
float	ptGen_
float	ptHat_
float	ptJet_
int	rank_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
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 ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
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)

Detailed Description

Definition at line 13 of file CaloMCTruthTreeProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 26 of file CaloMCTruthTreeProducer.cc.

References edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

{
  jets_          = consumes<CaloJetCollection>(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");
}

CaloMCTruthTreeProducer::~CaloMCTruthTreeProducer ( )

override

Definition at line 103 of file CaloMCTruthTreeProducer.cc.

{
  delete file_;
  delete mcTruthTree_;
}

Member Function Documentation

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

overridevirtual

Implements edm::EDAnalyzer.

Definition at line 61 of file CaloMCTruthTreeProducer.cc.

References GenEventInfoProduct::binningValues(), fwrapper::jets, njet, TCMET_cfi::radius, and findQualityFiles::rr.

{
  edm::Handle<GenJetCollection> genjets;
  edm::Handle<CaloJetCollection> jets;
  edm::Handle<GenEventInfoProduct> hEventInfo;
  CaloJetCollection::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();
         emfJet_  = i_matched->emEnergyFraction();
         dR_      = rmin;
         rank_    = njet;
         mcTruthTree_->Fill();
         njet++;
       }
    }
}

void CaloMCTruthTreeProducer::beginJob ( void  )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 34 of file CaloMCTruthTreeProducer.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("emfJet",     &emfJet_,     "emfJet_/F");
  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");
}

void CaloMCTruthTreeProducer::endJob ( void  )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 51 of file CaloMCTruthTreeProducer.cc.

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

Member Data Documentation

float CaloMCTruthTreeProducer::dR_

private

Definition at line 29 of file CaloMCTruthTreeProducer.h.

float CaloMCTruthTreeProducer::emfJet_

private

Definition at line 29 of file CaloMCTruthTreeProducer.h.

float CaloMCTruthTreeProducer::etaGen_

private

Definition at line 29 of file CaloMCTruthTreeProducer.h.

float CaloMCTruthTreeProducer::etaJet_

private

Definition at line 29 of file CaloMCTruthTreeProducer.h.

TFile* CaloMCTruthTreeProducer::file_

private

Definition at line 27 of file CaloMCTruthTreeProducer.h.

edm::EDGetTokenT<GenEventInfoProduct> CaloMCTruthTreeProducer::gen_

private

Definition at line 26 of file CaloMCTruthTreeProducer.h.

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

private

Definition at line 25 of file CaloMCTruthTreeProducer.h.

std::string CaloMCTruthTreeProducer::histogramFile_

private

Definition at line 23 of file CaloMCTruthTreeProducer.h.

edm::EDGetTokenT<reco::CaloJetCollection> CaloMCTruthTreeProducer::jets_

private

Definition at line 24 of file CaloMCTruthTreeProducer.h.

TTree* CaloMCTruthTreeProducer::mcTruthTree_

private

Definition at line 28 of file CaloMCTruthTreeProducer.h.

float CaloMCTruthTreeProducer::phiGen_

private

Definition at line 29 of file CaloMCTruthTreeProducer.h.

float CaloMCTruthTreeProducer::phiJet_

private

Definition at line 29 of file CaloMCTruthTreeProducer.h.

float CaloMCTruthTreeProducer::ptGen_

private

Definition at line 29 of file CaloMCTruthTreeProducer.h.

float CaloMCTruthTreeProducer::ptHat_

private

Definition at line 29 of file CaloMCTruthTreeProducer.h.

float CaloMCTruthTreeProducer::ptJet_

private

Definition at line 29 of file CaloMCTruthTreeProducer.h.

int CaloMCTruthTreeProducer::rank_

private

Definition at line 30 of file CaloMCTruthTreeProducer.h.