JetIdSelector< T > Class Template Reference

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
std::vector< bool > const &	recordProvenanceList () 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)
TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...
	~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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
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
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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...
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

template<typename T>
class JetIdSelector< T >

Definition at line 59 of file JetIdSelector.cc.

Member Typedef Documentation

JetCollection

template<typename T >

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

Definition at line 61 of file JetIdSelector.cc.

Constructor & Destructor Documentation

JetIdSelector()

template<typename T >

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

explicit

Definition at line 92 of file JetIdSelector.cc.

    : 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;
}

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.

~JetIdSelector()

template<typename T >

JetIdSelector< T >::~JetIdSelector ( )

override

Definition at line 128 of file JetIdSelector.cc.

                                 {
  if (jetIDFunctor)
    delete jetIDFunctor;
}

Member Function Documentation

endJob()

template<typename T >

void JetIdSelector< T >::endJob ( void  )

overridevirtual

Reimplemented from edm::EDProducer.

Definition at line 226 of file JetIdSelector.cc.

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

References gather_cfg::cout, and contentValuesCheck::ss.

produce()

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.

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

References reco::PFJet::chargedEmEnergy(), reco::PFJet::chargedHadronEnergy(), reco::PFJet::chargedMultiplicity(), reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), heavyIonCSV_trainingSettings::idx, iEvent, edm::HandleBase::isValid(), metsig::jet, PhotonMonitor_cff::jetId, singleTopDQM_cfi::jets, eostools::move(), reco::Jet::nConstituents(), reco::PFJet::neutralEmEnergy(), reco::PFJet::neutralHadronEnergy(), TriggerAnalyzer::passed, and objectSelection_cff::selectedJets.

Member Data Documentation

jetIDFunctor

template<typename T >

JetIDSelectionFunctor* JetIdSelector< T >::jetIDFunctor

private

Definition at line 79 of file JetIdSelector.cc.

Referenced by JetIdSelector< T >::JetIdSelector().

jetIDMap_

template<typename T >

edm::InputTag JetIdSelector< T >::jetIDMap_

private

Definition at line 74 of file JetIdSelector.cc.

moduleLabel_

template<typename T >

std::string JetIdSelector< T >::moduleLabel_

private

Definition at line 75 of file JetIdSelector.cc.

Referenced by ExternalGeneratorFilter.ExternalGeneratorFilter::insertInto(), Modules.SwitchProducer::insertInto(), Mixins._TypedParameterizable::insertInto(), Modules.ESSource::nameInProcessDesc_(), and Modules.ESProducer::nameInProcessDesc_().

nJetsPassed_

template<typename T >

unsigned int JetIdSelector< T >::nJetsPassed_

private

Definition at line 78 of file JetIdSelector.cc.

nJetsTot_

template<typename T >

unsigned int JetIdSelector< T >::nJetsTot_

private

Definition at line 77 of file JetIdSelector.cc.

qualityStr

template<typename T >

std::string JetIdSelector< T >::qualityStr

private

Definition at line 73 of file JetIdSelector.cc.

Referenced by JetIdSelector< T >::JetIdSelector().

src_

template<typename T >

edm::InputTag JetIdSelector< T >::src_

private

Definition at line 72 of file JetIdSelector.cc.

use_pfloose

template<typename T >

bool JetIdSelector< T >::use_pfloose

private

Definition at line 80 of file JetIdSelector.cc.

Referenced by JetIdSelector< T >::JetIdSelector().

use_pfmedium

template<typename T >

bool JetIdSelector< T >::use_pfmedium

private

Definition at line 81 of file JetIdSelector.cc.

Referenced by JetIdSelector< T >::JetIdSelector().

use_pftight

template<typename T >

bool JetIdSelector< T >::use_pftight

private

Definition at line 82 of file JetIdSelector.cc.

Referenced by JetIdSelector< T >::JetIdSelector().