Inheritance diagram for JetIdSelector< T >:

Public Types
typedef std::vector< T >	JetCollection

Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
void	endJob () override

	JetIdSelector (const edm::ParameterSet &iConfig)

void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override

	~JetIdSelector () override

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDProducer () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector < edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector < edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) 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

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > 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

bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
JetIDSelectionFunctor *	jetIDFunctor

edm::InputTag	jetIDMap_

std::string	moduleLabel_

unsigned int	nJetsPassed_

unsigned int	nJetsTot_

std::string	qualityStr

edm::InputTag	src_

bool	use_pfloose

bool	use_pfmedium

bool	use_pftight

Additional Inherited Members
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

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 ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

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 () noexcept

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

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

template<typename T>
class JetIdSelector< T >

Definition at line 59 of file JetIdSelector.cc.

Member Typedef Documentation

template<typename T >

typedef std::vector<T> JetIdSelector< T >::JetCollection

Definition at line 61 of file JetIdSelector.cc.

Constructor & Destructor Documentation

template<typename T >

JetIdSelector< T >::JetIdSelector ( const edm::ParameterSet & iConfig )

explicit

Definition at line 92 of file JetIdSelector.cc.

References Exception, JetIdSelector< T >::jetIDFunctor, JetIDSelectionFunctor::LOOSE, JetIDSelectionFunctor::LOOSE_AOD, JetIDSelectionFunctor::MINIMAL, JetIDSelectionFunctor::PURE09, JetIdSelector< T >::qualityStr, JetIDSelectionFunctor::TIGHT, JetIdSelector< T >::use_pfloose, JetIdSelector< T >::use_pfmedium, and JetIdSelector< T >::use_pftight.

     : src_(iConfig.getParameter<edm::InputTag>("src")),
       qualityStr(iConfig.getParameter<string>("idLevel")),
       jetIDMap_(iConfig.getUntrackedParameter<edm::InputTag>("jetIDMap", edm::InputTag("ak5JetID"))),
       moduleLabel_(iConfig.getParameter<string>("@module_label")),
       nJetsTot_(0),
       nJetsPassed_(0) {
   produces<JetCollection>();
 
   use_pfloose = false;
   use_pfmedium = false;
   use_pftight = false;
 
   if (qualityStr == "MINIMAL") {
     jetIDFunctor = new JetIDSelectionFunctor(JetIDSelectionFunctor::PURE09, JetIDSelectionFunctor::MINIMAL);
     use_pfloose = true;
   } else if (qualityStr == "LOOSE_AOD") {
     jetIDFunctor = new JetIDSelectionFunctor(JetIDSelectionFunctor::PURE09, JetIDSelectionFunctor::LOOSE_AOD);
     use_pfloose = true;
   } else if (qualityStr == "LOOSE") {
     jetIDFunctor = new JetIDSelectionFunctor(JetIDSelectionFunctor::PURE09, JetIDSelectionFunctor::LOOSE);
     use_pfloose = true;
   } else if (qualityStr == "MEDIUM") {
     jetIDFunctor = new JetIDSelectionFunctor(JetIDSelectionFunctor::PURE09, JetIDSelectionFunctor::LOOSE);
     // There is no medium quality for CaloJet !!
     use_pfmedium = true;
   } else if (qualityStr == "TIGHT") {
     jetIDFunctor = new JetIDSelectionFunctor(JetIDSelectionFunctor::PURE09, JetIDSelectionFunctor::TIGHT);
     use_pftight = true;
   } else
     throw cms::Exception("InvalidInput") << "Expect quality to be one of MINIMAL, LOOSE_AOD, LOOSE, MEDIUM, TIGHT"
                                          << std::endl;
 }

template<typename T >

JetIdSelector< T >::~JetIdSelector ( )

override

Definition at line 128 of file JetIdSelector.cc.

                                  {
   if (jetIDFunctor)
     delete jetIDFunctor;
 }

Member Function Documentation

template<typename T >

void JetIdSelector< T >::endJob ( void )

overridevirtual

Reimplemented from edm::EDProducer.

Definition at line 226 of file JetIdSelector.cc.

References gather_cfg::cout, moduleLabel_(), and contentValuesCheck::ss.

                               {
   stringstream ss;
   ss << "nJetsTot=" << nJetsTot_ << " nJetsPassed=" << nJetsPassed_
      << " fJetsPassed=" << 100. * (nJetsPassed_ / (double)nJetsTot_) << "%\n";
   cout << "++++++++++++++++++++++++++++++++++++++++++++++++++"
        << "\n"
        << moduleLabel_ << "(JetIdSelector) SUMMARY:\n"
        << ss.str() << "++++++++++++++++++++++++++++++++++++++++++++++++++" << endl;
 }

template<typename T >

void JetIdSelector< T >::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Implements edm::EDProducer.

Definition at line 139 of file JetIdSelector.cc.

References reco::PFJet::chargedEmEnergy(), reco::PFJet::chargedHadronEnergy(), reco::PFJet::chargedMultiplicity(), reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), edm::Event::getByLabel(), edm::HandleBase::isValid(), metsig::jet, fwrapper::jets, eostools::move(), reco::Jet::nConstituents(), reco::PFJet::neutralEmEnergy(), reco::PFJet::neutralHadronEnergy(), edm::Event::put(), and patRefSel_refMuJets_cfi::selectedJets.

                                                                             {
   auto selectedJets = std::make_unique<JetCollection>();
   edm::Handle<reco::JetView> hjets;  // uncorrected jets!
   iEvent.getByLabel(src_, hjets);
   auto const& jets = *hjets;
 
   // handle to the jet ID variables, filled in the loop
   edm::Handle<reco::JetIDValueMap> hJetIDMap;
   unsigned int idx = 0;
   bool passed = false;
 
   for (edm::View<reco::Jet>::const_iterator ibegin = jets.begin(), iend = jets.end(), iJet = ibegin; iJet != iend;
        ++iJet) {
     // initialize the boolean flag to false
     passed = false;
 
     //calculate the Calo jetID
     auto const& jet = *iJet;
     const std::type_info& type = typeid(jet);
     if (type == typeid(reco::CaloJet)) {
       const reco::CaloJet& calojet = static_cast<const reco::CaloJet&>(jet);
       edm::RefToBase<reco::Jet> jetRef = jets.refAt(idx);
       if (not hJetIDMap.isValid()) {
         iEvent.getByLabel(jetIDMap_, hJetIDMap);
       }
       reco::JetID const& jetId = (*hJetIDMap)[jetRef];
       passed = (*jetIDFunctor)(calojet, jetId);
     }
 
     //calculate the PF jetID
     if (type == typeid(reco::PFJet)) {
       const reco::PFJet& pfjet = static_cast<const reco::PFJet&>(jet);
       bool passingLoose = false;
       bool passingMedium = false;
       bool passingTight = false;
       bool ThisIsClean = true;
       //apply following only if |eta|<2.4: CHF>0, CEMF<0.99, chargedMultiplicity>0
       if ((pfjet.chargedHadronEnergy() / pfjet.energy()) <= 0.0 && fabs(pfjet.eta()) < 2.4)
         ThisIsClean = false;
       if ((pfjet.chargedEmEnergy() / pfjet.energy()) >= 0.99 && fabs(pfjet.eta()) < 2.4)
         ThisIsClean = false;
       if (pfjet.chargedMultiplicity() <= 0 && fabs(pfjet.eta()) < 2.4)
         ThisIsClean = false;
 
       // always require #Constituents > 1
       if (pfjet.nConstituents() <= 1)
         ThisIsClean = false;
 
       if (ThisIsClean && (pfjet.neutralHadronEnergy() / pfjet.energy()) < 0.99 &&
           (pfjet.neutralEmEnergy() / pfjet.energy()) < 0.99)
         passingLoose = true;
 
       if (ThisIsClean && (pfjet.neutralHadronEnergy() / pfjet.energy()) < 0.95 &&
           (pfjet.neutralEmEnergy() / pfjet.energy()) < 0.95)
         passingMedium = true;
 
       if (ThisIsClean && (pfjet.neutralHadronEnergy() / pfjet.energy()) < 0.90 &&
           (pfjet.neutralEmEnergy() / pfjet.energy()) < 0.90)
         passingTight = true;
 
       if (use_pfloose && passingLoose)
         passed = true;
       if (use_pfmedium && passingMedium)
         passed = true;
       if (use_pftight && passingTight)
         passed = true;
     }
 
     if (type == typeid(reco::GenJet) || type == typeid(reco::JPTJet)) {
       edm::LogWarning("JetId") << "Criteria for jets other than CaloJets and PFJets are not yet implemented";
       passed = true;
     }  // close GenJet, JPT jet
 
     const T& goodJet = static_cast<const T&>(jet);
     if (passed)
       selectedJets->push_back(goodJet);
 
     idx++;
   }  // close jet iterator
 
   nJetsTot_ += jets.size();
   nJetsPassed_ += selectedJets->size();
   iEvent.put(std::move(selectedJets));
 }