SelectionStep< Object > Class Template Reference

Templated helper class to allow a selection on a certain object collection. More...

#include <TopDQMHelpers.h>

Public Member Functions
bool	select (const edm::Event &event)
	apply selection More...
bool	select (const edm::Event &event, const std::string &type)
	apply selection with special treatment for PFCandidates More...
bool	select (const edm::Event &event, const edm::EventSetup &setup)
	apply selection override for jets More...
	SelectionStep (const edm::ParameterSet &cfg, edm::ConsumesCollector &&iC)
	default constructor More...
bool	selectVertex (const edm::Event &event)
	~SelectionStep ()
	default destructor More...

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...
double	eidCutValue_
edm::EDGetTokenT< edm::ValueMap< float > >	electronId_
	electronId label as extra selection type More...
edm::EDGetTokenT< edm::View< reco::GsfElectron > >	gsfEs_
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 SelectionStep< Object >

Templated helper class to allow a selection on a certain object collection.

"DQM/Physics/interface/TopDQMHelpers.h"

Templated helper class to allow a selection on a certain object collection, which may be monitored by a separate class afterwards. The class wraps and slightly extends the features of the StringCutParser to allow also to apply event based selections, according to a minimal or maximal number of elements in the collection after the object selection has been applied. It takes an edm::ParameterSet in the constructor, which should contain the following elements:

• src : the input collection (mandatory).
• select : the selection string (mandatory).
• min : whether there is a min value on which to reject the whole event after the selection (optional).
• max : whether there is a max value on which to reject the whole event after the selection (optional).
• electronId : input tag of an electronId association map and cut value (optional).
• jetCorrector : label of jet corrector (optional).
• jetBTagger : parameters defining the btag algorithm and working point of choice (optional).
• jetID : parameters defining the jetID value map and selection (optional).

The parameters src and select are mandatory. The parameters min and max are optional. The parameters electronId and jetCorrector are optional. They are added to keep the possibility to apply selections on id'ed electrons or on corrected jets. They may be omitted in the PSet for simplification reasons if not needed at any time. They are not effiective for other object collections but electrons or jets. If none of the two parameters min or max is found in the event the select function returns true if at least one object fullfilled the requirements.

The class has one template value, which is the object collection to apply the selection on. This has to be parsed to the StringCutParser class. The function select is overrided for jets to circumvent problems with the template specialisation. Note that for MET not type1 or muon corrections are supported on reco candidates.

Definition at line 237 of file TopDQMHelpers.h.

Constructor & Destructor Documentation

template<typename Object >

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

default constructor

Definition at line 289 of file TopDQMHelpers.h.

References SelectionStep< Object >::btagLabel_, SelectionStep< Object >::btagWorkingPoint_, SelectionStep< Object >::eidCutValue_, SelectionStep< Object >::electronId_, edm::ParameterSet::exists(), edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), SelectionStep< Object >::gsfEs_, SelectionStep< Object >::jetCorrector_, singleTopDQM_cfi::jetID, SelectionStep< Object >::jetIDLabel_, SelectionStep< Object >::jetIDSelect_, SelectionStep< Object >::max_, SelectionStep< Object >::min_, SelectionStep< Object >::src_, and AlCaHLTBitMon_QueryRunRegistry::string.

    : select_(cfg.getParameter<std::string>("select")), jetIDSelect_(nullptr) {
  src_ = iC.consumes<edm::View<Object> >(cfg.getParameter<edm::InputTag>("src"));
  // 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;
  std::string mygSF = "gedGsfElectrons";
  gsfEs_ = iC.consumes<edm::View<reco::GsfElectron> >(cfg.getUntrackedParameter<edm::InputTag>("myGSF", mygSF));
  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"));
    eidCutValue_ = elecId.getParameter<double>("cutValue");
  }
  if (cfg.exists("jetCorrector")) {
    jetCorrector_ = cfg.getParameter<std::string>("jetCorrector");
  }
  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");
  }
  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 >

SelectionStep< Object >::~SelectionStep ( )

inline

default destructor

Definition at line 242 of file TopDQMHelpers.h.

References singleTopDQM_cfi::select, singleTopDQM_cfi::setup, and AlCaHLTBitMon_QueryRunRegistry::string.

{};

Member Function Documentation

template<typename Object >

bool SelectionStep< Object >::select

(

const edm::Event &

event

)

apply selection

Definition at line 319 of file TopDQMHelpers.h.

References accept(), SelectionStep< Object >::eidCutValue_, topSingleLeptonDQM_PU_cfi::electronId, SelectionStep< Object >::electronId_, edm::Event::getByToken(), training_settings::idx, edm::EDGetTokenT< T >::isUninitialized(), SelectionStep< Object >::max_, SelectionStep< Object >::min_, dqmiodumpmetadata::n, getGTfromDQMFile::obj, SelectionStep< Object >::select_, TrackRefitter_38T_cff::src, and SelectionStep< 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 : ((double)(*electronId)[src->refAt(idx)] >= eidCutValue_)) {
        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 SelectionStep< Object >::select	(	const edm::Event &	event,
		const std::string &	type
	)

apply selection with special treatment for PFCandidates

Definition at line 355 of file TopDQMHelpers.h.

References accept(), SelectionStep< Object >::eidCutValue_, topSingleLeptonDQM_PU_cfi::electronId, SelectionStep< Object >::electronId_, edm::Event::getByToken(), reco::PFCandidate::gsfElectronRef(), training_settings::idx, edm::Ref< C, T, F >::isNonnull(), edm::EDGetTokenT< T >::isUninitialized(), SelectionStep< Object >::max_, SelectionStep< Object >::min_, reco::PFCandidate::muonRef(), dqmiodumpmetadata::n, getGTfromDQMFile::obj, SelectionStep< Object >::select_, TrackRefitter_38T_cff::src, and SelectionStep< 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;

  // special for gsfElectron
  edm::Handle<edm::View<reco::GsfElectron> > elecs_gsf;

  // 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 PF candidates
    if (dynamic_cast<const reco::PFCandidate*>(&*obj)) {
      reco::PFCandidate objtmp = dynamic_cast<const reco::PFCandidate&>(*obj);

      if (objtmp.muonRef().isNonnull() && type == "muon") {
        if (select_(*obj)) {
          ++n;
        }
      } else if (objtmp.gsfElectronRef().isNonnull() && type == "electron") {
        if (select_(*obj)) {
          if (electronId_.isUninitialized()) {
            ++n;
          } else if (((double)(*electronId)[obj->gsfElectronRef()] >= eidCutValue_)) {
            ++n;
          }
        }
        //        idx_gsf++;
      }
    }

    // special treatment for electrons
    else if (dynamic_cast<const reco::GsfElectron*>(&*obj)) {
      unsigned int idx = obj - src->begin();
      if (electronId_.isUninitialized() ? true : ((double)(*electronId)[src->refAt(idx)] >= eidCutValue_)) {
        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 SelectionStep< Object >::select	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

apply selection override for jets

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

Definition at line 439 of file TopDQMHelpers.h.

References accept(), SelectionStep< Object >::btagLabel_, SelectionStep< Object >::btagWorkingPoint_, JetCorrector::correction(), mitigatedMETSequence_cff::corrector, edm::EventSetup::find(), edm::Event::getByToken(), JetCorrector::getJetCorrector(), training_settings::idx, edm::EDGetTokenT< T >::isUninitialized(), metsig::jet, SelectionStep< Object >::jetCorrector_, singleTopDQM_cfi::jetID, SelectionStep< Object >::jetIDLabel_, SelectionStep< Object >::jetIDSelect_, SelectionStep< Object >::max_, SelectionStep< Object >::min_, dqmiodumpmetadata::n, getGTfromDQMFile::obj, SelectionStep< Object >::pvs_, SelectionStep< Object >::select_, TrackRefitter_38T_cff::src, and SelectionStep< 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 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("TopDQMHelpers") << "\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";
    }
  }
  // determine multiplicity of selected objects
  int n = 0;
  for (typename edm::View<Object>::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
        Object 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 SelectionStep< Object >::selectVertex

(

const edm::Event &

event

)

Definition at line 414 of file TopDQMHelpers.h.

References accept(), topSingleLeptonDQM_PU_cfi::electronId, SelectionStep< Object >::electronId_, edm::Event::getByToken(), edm::EDGetTokenT< T >::isUninitialized(), SelectionStep< Object >::max_, SelectionStep< Object >::min_, dqmiodumpmetadata::n, getGTfromDQMFile::obj, SelectionStep< Object >::select_, TrackRefitter_38T_cff::src, and SelectionStep< 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> SelectionStep< Object >::btagLabel_

private

choice for b-tag as extra selection type

Definition at line 273 of file TopDQMHelpers.h.

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

template<typename Object >

double SelectionStep< Object >::btagWorkingPoint_

private

choice of b-tag working point as extra selection type

Definition at line 275 of file TopDQMHelpers.h.

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

template<typename Object >

double SelectionStep< Object >::eidCutValue_

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 269 of file TopDQMHelpers.h.

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

template<typename Object >

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

private

electronId label as extra selection type

Definition at line 257 of file TopDQMHelpers.h.

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

template<typename Object >

edm::EDGetTokenT<edm::View<reco::GsfElectron> > SelectionStep< Object >::gsfEs_

private

Definition at line 280 of file TopDQMHelpers.h.

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

template<typename Object >

std::string SelectionStep< Object >::jetCorrector_

private

jet corrector as extra selection type

Definition at line 271 of file TopDQMHelpers.h.

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

template<typename Object >

edm::EDGetTokenT<reco::JetIDValueMap> SelectionStep< Object >::jetIDLabel_

private

jetID as an extra selection type

Definition at line 277 of file TopDQMHelpers.h.

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

template<typename Object >

StringCutObjectSelector<reco::JetID>* SelectionStep< Object >::jetIDSelect_

private

selection string on the jetID

Definition at line 284 of file TopDQMHelpers.h.

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

template<typename Object >

int SelectionStep< Object >::max_

private

Definition at line 255 of file TopDQMHelpers.h.

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

template<typename Object >

int SelectionStep< Object >::min_

private

min/max for object multiplicity

Definition at line 255 of file TopDQMHelpers.h.

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

template<typename Object >

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

private

Definition at line 279 of file TopDQMHelpers.h.

Referenced by SelectionStep< Object >::select().

template<typename Object >

StringCutObjectSelector<Object> SelectionStep< Object >::select_

private

string cut selector

Definition at line 282 of file TopDQMHelpers.h.

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

template<typename Object >

edm::EDGetTokenT<edm::View<Object> > SelectionStep< Object >::src_

private

input collection

Definition at line 253 of file TopDQMHelpers.h.

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