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
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

virtual	~JetIdSelector ()

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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

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, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

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

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

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)

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

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 62 of file JetIdSelector.cc.

Constructor & Destructor Documentation

template<typename T >

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

explicit

Definition at line 96 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 ( )

virtual

Definition at line 144 of file JetIdSelector.cc.

References JetIdSelector< T >::jetIDFunctor.

                                 {
   if(jetIDFunctor) delete jetIDFunctor;
 }

Member Function Documentation

template<typename T >

void JetIdSelector< T >::endJob ( void )

overridevirtual

Reimplemented from edm::EDProducer.

Definition at line 251 of file JetIdSelector.cc.

References gather_cfg::cout, JetIdSelector< T >::moduleLabel_, JetIdSelector< T >::nJetsPassed_, and JetIdSelector< T >::nJetsTot_.

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

override

Definition at line 156 of file JetIdSelector.cc.

References reco::PFJet::chargedEmEnergy(), reco::PFJet::chargedHadronEnergy(), reco::PFJet::chargedMultiplicity(), reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), edm::Event::getByLabel(), training_settings::idx, JetIdSelector< T >::jetIDMap_, fwrapper::jets, eostools::move(), reco::Jet::nConstituents(), reco::PFJet::neutralEmEnergy(), reco::PFJet::neutralHadronEnergy(), JetIdSelector< T >::nJetsPassed_, JetIdSelector< T >::nJetsTot_, TriggerAnalyzer::passed, edm::Event::put(), objectSelection_cff::selectedJets, JetIdSelector< T >::src_, JetIdSelector< T >::use_pfloose, JetIdSelector< T >::use_pfmedium, and JetIdSelector< T >::use_pftight.

 {
   auto selectedJets = std::make_unique<JetCollection>();
   edm::Handle<reco::JetView> jets;  // uncorrected jets!
   iEvent.getByLabel(src_,jets);
 
 
  // handle to the jet ID variables
    edm::Handle<reco::JetIDValueMap> hJetIDMap;
 
    if(typeid((*jets)[0]) == typeid(reco::CaloJet)) 
      iEvent.getByLabel( jetIDMap_, 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
      const std::type_info & type = typeid((*jets)[idx]);
      if( type == typeid(reco::CaloJet) ) {
        const reco::CaloJet calojet = static_cast<const reco::CaloJet &>((*jets)[idx]);
        edm::RefToBase<reco::Jet> jetRef = jets->refAt(idx);
        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 &>((*jets)[idx]);
        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&>((*jets)[idx]);
      if(passed) selectedJets->push_back( goodJet );
 
      idx++;
    } // close jet iterator
 
 
 
   nJetsTot_  +=jets->size();
   nJetsPassed_+=selectedJets->size();  
   iEvent.put(std::move(selectedJets));
 }