Inheritance diagram for GenJetFlavourTableProducer:

Public Member Functions
	GenJetFlavourTableProducer (const edm::ParameterSet &iConfig)

	~GenJetFlavourTableProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
void	produce (edm::Event &, edm::EventSetup const &) override

Private Attributes
const StringCutObjectSelector< reco::GenJet >	cut_

const double	deltaR_

edm::EDGetTokenT< reco::JetFlavourInfoMatchingCollection >	jetFlavourInfosToken_

std::string	name_

edm::EDGetTokenT< std::vector< reco::GenJet > >	src_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Definition at line 24 of file GenJetFlavourTableProducer.cc.

Constructor & Destructor Documentation

GenJetFlavourTableProducer::GenJetFlavourTableProducer ( const edm::ParameterSet & iConfig )

inlineexplicit

Definition at line 26 of file GenJetFlavourTableProducer.cc.

                                                                             :
             name_(iConfig.getParameter<std::string>("name")),
             src_(consumes<std::vector<reco::GenJet> >(iConfig.getParameter<edm::InputTag>("src"))),
             cut_(iConfig.getParameter<std::string>("cut"), true),
             deltaR_(iConfig.getParameter<double>("deltaR")),
             jetFlavourInfosToken_(consumes<reco::JetFlavourInfoMatchingCollection>(iConfig.getParameter<edm::InputTag>("jetFlavourInfos")))
         {
             produces<nanoaod::FlatTable>();
         }

GenJetFlavourTableProducer::~GenJetFlavourTableProducer ( )

inlineoverride

Definition at line 36 of file GenJetFlavourTableProducer.cc.

36 {};

Member Function Documentation

static void GenJetFlavourTableProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

inlinestatic

Definition at line 38 of file GenJetFlavourTableProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), produce(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                   {
             edm::ParameterSetDescription desc;
             desc.add<edm::InputTag>("src")->setComment("input genJet collection");
             desc.add<edm::InputTag>("jetFlavourInfos")->setComment("input flavour info collection");
             desc.add<std::string>("name")->setComment("name of the genJet FlatTable we are extending with flavour information");
             desc.add<std::string>("cut")->setComment("cut on input genJet collection");
             desc.add<double>("deltaR")->setComment("deltaR to match genjets");
             descriptions.add("genJetFlavourTable", desc);
         }

void GenJetFlavourTableProducer::produce	(	edm::Event &	,
		edm::EventSetup const &
	)

overrideprivate

Definition at line 61 of file GenJetFlavourTableProducer.cc.

References cut_, DEFINE_FWK_MODULE, deltaR(), deltaR_, edm::Event::getByToken(), jets_cff::hadronFlavour, nanoaod::FlatTable::IntColumn, metsig::jet, GenHFHadronMatcher_cfi::jetFlavourInfos, jetFlavourInfosToken_, fwrapper::jets, electrons_cff::matched, eostools::move(), name_, jets_cff::partonFlavour, edm::Event::put(), src_, and nanoaod::FlatTable::UInt8Column.

Referenced by fillDescriptions().

 {
     edm::Handle<reco::GenJetCollection> jets;
     iEvent.getByToken(src_, jets);
     
     edm::Handle<reco::JetFlavourInfoMatchingCollection> jetFlavourInfos;
     iEvent.getByToken(jetFlavourInfosToken_, jetFlavourInfos);
 
     unsigned int ncand = 0;
     std::vector<int> partonFlavour;
     std::vector<uint8_t> hadronFlavour;
     
     for (const reco::GenJet & jet : *jets) {
       if (!cut_(jet)) continue;
       ++ncand;
       bool matched = false;
       for (const reco::JetFlavourInfoMatching & jetFlavourInfoMatching : *jetFlavourInfos) {
         if (deltaR(jet.p4(), jetFlavourInfoMatching.first->p4()) < deltaR_) {
           partonFlavour.push_back(jetFlavourInfoMatching.second.getPartonFlavour());
           hadronFlavour.push_back(jetFlavourInfoMatching.second.getHadronFlavour());
           matched = true;
           break;
         }
       }
       if (!matched) {
         partonFlavour.push_back(0);
         hadronFlavour.push_back(0);
       }
     }
 
     auto tab  = std::make_unique<nanoaod::FlatTable>(ncand, name_, false, true);
     tab->addColumn<int>("partonFlavour", partonFlavour, "flavour from parton matching", nanoaod::FlatTable::IntColumn);
     tab->addColumn<uint8_t>("hadronFlavour", hadronFlavour, "flavour from hadron ghost clustering", nanoaod::FlatTable::UInt8Column);
 
     iEvent.put(std::move(tab));
 }