#include <TopHLTOfflineDQMHelper.h>

Inheritance diagram for SelectionStepHLT< Object >:

Public Member Functions
template<>
bool	select (const edm::Event &event, const edm::EventSetup &setup)
	apply selection (w/o using the template class Object), override for jets More...

bool	select (const edm::Event &event) override
	apply selection More...

bool	select (const edm::Event &event, const edm::EventSetup &setup) override
	apply selection override for jets More...

	SelectionStepHLT (const edm::ParameterSet &cfg, edm::ConsumesCollector &&iC)
	default constructor More...

bool	selectVertex (const edm::Event &event)

	~SelectionStepHLT () override=default
	default destructor More...

Public Member Functions inherited from SelectionStepHLTBase
virtual	~SelectionStepHLTBase ()=default

Private Attributes
edm::EDGetTokenT< reco::JetTagCollection >	btagLabel_
	choice for b-tag as extra selection type More...

double	btagWorkingPoint_
	choice of b-tag working point as extra selection type More...

int	eidPattern_

edm::EDGetTokenT< edm::ValueMap< float > >	electronId_
	electronId label as extra selection type More...

std::string	jetCorrector_
	jet corrector as extra selection type More...

edm::EDGetTokenT< reco::JetIDValueMap >	jetIDLabel_
	jetID as an extra selection type More...

StringCutObjectSelector< reco::JetID > *	jetIDSelect_
	selection string on the jetID More...

int	max_

int	min_
	min/max for object multiplicity More...

edm::EDGetTokenT< edm::View< reco::Vertex > >	pvs_

StringCutObjectSelector< Object >	select_
	string cut selector More...

edm::EDGetTokenT< edm::View< Object > >	src_
	input collection More...

Detailed Description

template<typename Object>
class SelectionStepHLT< Object >

Definition at line 190 of file TopHLTOfflineDQMHelper.h.

Constructor & Destructor Documentation

template<typename Object >

SelectionStepHLT< Object >::SelectionStepHLT	(	const edm::ParameterSet &	cfg,
		edm::ConsumesCollector &&	iC
	)

default constructor

Definition at line 239 of file TopHLTOfflineDQMHelper.h.

                                                                                                 :
   src_( iC.consumes< edm::View<Object> >(cfg.getParameter<edm::InputTag>( "src"   ))),
   select_( cfg.getParameter<std::string>("select")),
   jetIDSelect_(nullptr)
 {
   // construct min/max if the corresponding params
   // exist otherwise they are initialized with -1
   cfg.exists("min") ? min_= cfg.getParameter<int>("min") : min_= -1;
   cfg.exists("max") ? max_= cfg.getParameter<int>("max") : max_= -1;
   // read electron extras if they exist
   if(cfg.existsAs<edm::ParameterSet>("electronId")){ 
     edm::ParameterSet elecId=cfg.getParameter<edm::ParameterSet>("electronId");
     electronId_= iC.consumes< edm::ValueMap<float> >(elecId.getParameter<edm::InputTag>("src"));
     eidPattern_= elecId.getParameter<int>("pattern");
   }
   // read jet corrector label if it exists
   if(cfg.exists("jetCorrector")){ jetCorrector_= cfg.getParameter<std::string>("jetCorrector"); }
   // read btag information if it exists
   if(cfg.existsAs<edm::ParameterSet>("jetBTagger")){
     edm::ParameterSet jetBTagger=cfg.getParameter<edm::ParameterSet>("jetBTagger");
     btagLabel_= iC.consumes< reco::JetTagCollection >(jetBTagger.getParameter<edm::InputTag>("label"));
     btagWorkingPoint_=jetBTagger.getParameter<double>("workingPoint");
   }
   // read jetID information if it exists
   if(cfg.existsAs<edm::ParameterSet>("jetID")){
     edm::ParameterSet jetID=cfg.getParameter<edm::ParameterSet>("jetID");
     jetIDLabel_ = iC.consumes< reco::JetIDValueMap >(jetID.getParameter<edm::InputTag>("label"));
     jetIDSelect_= new StringCutObjectSelector<reco::JetID>(jetID.getParameter<std::string>("select"));
   }
 }

template<typename Object >

SelectionStepHLT< Object >::~SelectionStepHLT ( )

overridedefault

default destructor

Member Function Documentation

template<>

bool SelectionStepHLT< reco::Jet >::select	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

virtual

apply selection (w/o using the template class Object), override for jets

Reimplemented from SelectionStepHLTBase.

Definition at line 9 of file TopHLTOfflineDQMHelper.cc.

References accept(), JetCorrector::correction(), mitigatedMETSequence_cff::corrector, edm::EventSetup::find(), edm::Event::getByToken(), JetCorrector::getJetCorrector(), training_settings::idx, metsig::jet, gen::n, hgcalPlots::obj, reco::Jet::scaleEnergy(), and TrackRefitter_38T_cff::src.

 {
   // fetch input collection
   edm::Handle<edm::View<reco::Jet> > src; 
   if( !event.getByToken(src_, src) ) return false;
 
   // load btag collection if configured such
   // NOTE that the JetTagCollection needs an
   // edm::View to reco::Jets; we have to add
   // another Handle bjets for this purpose
   edm::Handle<edm::View<reco::Jet> > bjets; 
   edm::Handle<reco::JetTagCollection> btagger;
   edm::Handle<edm::View<reco::Vertex> > pvertex; 
   if(!btagLabel_.isUninitialized()){ 
     if( !event.getByToken(src_, bjets) ) return false;
     if( !event.getByToken(btagLabel_, btagger) ) return false;
     if( !event.getByToken(pvs_, pvertex) ) return false;
   }
 
   // load jetID value map if configured such 
   edm::Handle<reco::JetIDValueMap> jetID;
   if(jetIDSelect_){
     if( !event.getByToken(jetIDLabel_, jetID) ) return false;
 
   }
 
   // load jet corrector if configured such
   const JetCorrector* corrector=nullptr;
   if(!jetCorrector_.empty()){
     // check whether a jet correcto is in the event setup or not
     if(setup.find( edm::eventsetup::EventSetupRecordKey::makeKey<JetCorrectionsRecord>() )){
       corrector = JetCorrector::getJetCorrector(jetCorrector_, setup);
     }
     else{
       edm::LogVerbatim( "TopHLTOfflineDQMHelper" ) 
         << "\n"
         << "------------------------------------------------------------------------------------- \n"
         << " No JetCorrectionsRecord available from EventSetup:                                   \n" 
         << "  - Jets will not be corrected.                                                       \n"
         << "  - If you want to change this add the following lines to your cfg file               \n"
         << "                                                                                      \n"
         << "  ## load jet corrections                                                             \n"
         << "  process.load(\"JetMETCorrections.Configuration.JetCorrectionServicesAllAlgos_cff\") \n"
         << "  process.prefer(\"ak5CaloL2L3\")                                                     \n"
         << "                                                                                      \n"
         << "------------------------------------------------------------------------------------- \n";
     }
   }
   // determine multiplicity of selected objects
   int n=0;
   for(typename edm::View<reco::Jet>::const_iterator obj=src->begin(); obj!=src->end(); ++obj){
     // check for chosen btag discriminator to be above the 
     // corresponding working point if configured such 
     unsigned int idx = obj-src->begin();
     if( btagLabel_.isUninitialized() ? true : (*btagger)[bjets->refAt(idx)]>btagWorkingPoint_ ){   
       bool passedJetID=true;
       // check jetID for calo jets
       if( jetIDSelect_ && dynamic_cast<const reco::CaloJet*>(src->refAt(idx).get())){
     passedJetID=(*jetIDSelect_)((*jetID)[src->refAt(idx)]);
       }
       if(passedJetID){
     // scale jet energy if configured such
         reco::Jet jet=*obj; jet.scaleEnergy(corrector ? corrector->correction(*obj) : 1.);
     if(select_(jet))++n;
       }
     }
   }
   bool accept=(min_>=0 ? n>=min_:true) && (max_>=0 ? n<=max_:true);
   return (min_<0 && max_<0) ? (n>0):accept;
 }

template<typename Object >

bool SelectionStepHLT< Object >::select ( const edm::Event & event )

overridevirtual

apply selection

Reimplemented from SelectionStepHLTBase.

Definition at line 272 of file TopHLTOfflineDQMHelper.h.

References accept(), SelectionStepHLT< Object >::eidPattern_, SelectionStepHLT< Object >::electronId_, edm::Event::getByToken(), training_settings::idx, createfilelist::int, edm::EDGetTokenT< T >::isUninitialized(), SelectionStepHLT< Object >::max_, SelectionStepHLT< Object >::min_, gen::n, hgcalPlots::obj, SelectionStepHLT< Object >::select_, TrackRefitter_38T_cff::src, and SelectionStepHLT< Object >::src_.

Referenced by Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), Vispa.Views.LineDecayView.LineDecayContainer::mousePressEvent(), Vispa.Gui.PortConnection.PointToPointConnection::mousePressEvent(), Vispa.Gui.VispaWidget.VispaWidget::mousePressEvent(), and Vispa.Views.AbstractView.AbstractView::restoreSelection().

 {
   // fetch input collection
   edm::Handle<edm::View<Object> > src; 
   if( !event.getByToken(src_, src) ) return false;
 
   // load electronId value map if configured such
   edm::Handle<edm::ValueMap<float> > electronId;
   if(!electronId_.isUninitialized()) {
     if( !event.getByToken(electronId_, electronId) ) return false;
   }
 
   // determine multiplicity of selected objects
   int n=0;
   for(typename edm::View<Object>::const_iterator obj=src->begin(); obj!=src->end(); ++obj){
     // special treatment for electrons
     if(dynamic_cast<const reco::GsfElectron*>(&*obj)){
       unsigned int idx = obj-src->begin();
       if( electronId_.isUninitialized() ? true : ((int)(*electronId)[src->refAt(idx)] & eidPattern_) ){   
     if(select_(*obj))++n;
       }
     }
     // normal treatment
     else{
       if(select_(*obj))++n;
     }
   }
   bool accept=(min_>=0 ? n>=min_:true) && (max_>=0 ? n<=max_:true);
   return (min_<0 && max_<0) ? (n>0):accept;
 }

template<typename Object >

bool SelectionStepHLT< Object >::select	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

overridevirtual

apply selection override for jets

Reimplemented from SelectionStepHLTBase.

Definition at line 326 of file TopHLTOfflineDQMHelper.h.

References Exception.

Referenced by Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), Vispa.Views.LineDecayView.LineDecayContainer::mousePressEvent(), Vispa.Gui.PortConnection.PointToPointConnection::mousePressEvent(), Vispa.Gui.VispaWidget.VispaWidget::mousePressEvent(), and Vispa.Views.AbstractView.AbstractView::restoreSelection().

 {
   throw cms::Exception("SelectionStepHLT") << "you fail" << std::endl;
   return false;
 }

template<typename Object >

bool SelectionStepHLT< Object >::selectVertex ( const edm::Event & event )

Definition at line 303 of file TopHLTOfflineDQMHelper.h.

References accept(), SelectionStepHLT< Object >::electronId_, edm::Event::getByToken(), edm::EDGetTokenT< T >::isUninitialized(), SelectionStepHLT< Object >::max_, SelectionStepHLT< Object >::min_, gen::n, hgcalPlots::obj, SelectionStepHLT< Object >::select_, TrackRefitter_38T_cff::src, and SelectionStepHLT< Object >::src_.

 {
   // fetch input collection
   edm::Handle<edm::View<Object> > src; 
   if( !event.getByToken(src_, src) ) return false;
 
   // load electronId value map if configured such
   edm::Handle<edm::ValueMap<float> > electronId;
   if(!electronId_.isUninitialized()) {
     if( !event.getByToken(electronId_, electronId) ) return false;
   }
 
   // determine multiplicity of selected objects
   int n=0;
   for(typename edm::View<Object>::const_iterator obj=src->begin(); obj!=src->end(); ++obj){
    
       if(select_(*obj))++n;
   }
   bool accept=(min_>=0 ? n>=min_:true) && (max_>=0 ? n<=max_:true);
   return (min_<0 && max_<0) ? (n>0):accept;
 }

Member Data Documentation

template<typename Object >

edm::EDGetTokenT< reco::JetTagCollection > SelectionStepHLT< Object >::btagLabel_

private

choice for b-tag as extra selection type

Definition at line 223 of file TopHLTOfflineDQMHelper.h.

Referenced by SelectionStepHLT< Object >::SelectionStepHLT().

template<typename Object >

double SelectionStepHLT< Object >::btagWorkingPoint_

private

choice of b-tag working point as extra selection type

Definition at line 225 of file TopHLTOfflineDQMHelper.h.

Referenced by SelectionStepHLT< Object >::SelectionStepHLT().

template<typename Object >

int SelectionStepHLT< Object >::eidPattern_

private

electronId pattern we expect the following pattern: 0: fails 1: passes electron ID only 2: passes electron Isolation only 3: passes electron ID and Isolation only 4: passes conversion rejection 5: passes conversion rejection and ID 6: passes conversion rejection and Isolation 7: passes the whole selection As described on https://twiki.cern.ch/twiki/bin/view/CMS/SimpleCutBasedEleID

Definition at line 219 of file TopHLTOfflineDQMHelper.h.

Referenced by SelectionStepHLT< Object >::select(), and SelectionStepHLT< Object >::SelectionStepHLT().

template<typename Object >

edm::EDGetTokenT< edm::ValueMap<float> > SelectionStepHLT< Object >::electronId_

private

electronId label as extra selection type

Definition at line 208 of file TopHLTOfflineDQMHelper.h.

Referenced by SelectionStepHLT< Object >::select(), SelectionStepHLT< Object >::SelectionStepHLT(), and SelectionStepHLT< Object >::selectVertex().

template<typename Object >

std::string SelectionStepHLT< Object >::jetCorrector_

private

jet corrector as extra selection type