EgammaL1TkIsolation Class Reference

#include <EgammaL1TkIsolation.h>

Classes
struct	TrkCuts

Public Member Functions
std::pair< int, double >	calIsol (const reco::TrackBase &trk, const L1TrackCollection &l1Tks) const
std::pair< int, double >	calIsol (const double objEta, const double objPhi, const double objZ, const L1TrackCollection &l1Tks) const
template<typename... Args>
double	calIsolPt (Args &&... args) const
	EgammaL1TkIsolation (const edm::ParameterSet &para)

Static Public Member Functions
static void	fillPSetDescription (edm::ParameterSetDescription &desc)
static edm::ParameterSetDescription	makePSetDescription ()

Private Member Functions
size_t	etaBinNr (double eta) const

Static Private Member Functions
static bool	passTrkSel (const L1Track &trk, const double trkPt, const TrkCuts &cuts, const double objEta, const double objPhi, const double objZ)

Private Attributes
std::vector< double >	etaBoundaries_
std::vector< TrkCuts >	trkCuts_
bool	useAbsEta_

Detailed Description

Definition at line 13 of file EgammaL1TkIsolation.h.

Constructor & Destructor Documentation

EgammaL1TkIsolation()

EgammaL1TkIsolation::EgammaL1TkIsolation ( const edm::ParameterSet &  para )

explicit

Definition at line 7 of file EgammaL1TkIsolation.cc.

References etaBoundaries_, Exception, edm::ParameterSet::getParameter(), submitPVValidationJobs::params, and trkCuts_.

    : useAbsEta_(para.getParameter<bool>("useAbsEta")),
      etaBoundaries_(para.getParameter<std::vector<double>>("etaBoundaries")) {
  const auto& trkCutParams = para.getParameter<std::vector<edm::ParameterSet>>("trkCuts");
  for (const auto& params : trkCutParams) {
    trkCuts_.emplace_back(TrkCuts(params));
  }
  if (etaBoundaries_.size() + 1 != trkCuts_.size()) {
    throw cms::Exception("ConfigError") << "EgammaL1TkIsolation: etaBoundaries parameters size ("
                                        << etaBoundaries_.size()
                                        << ") should be one less than the size of trkCuts VPSet (" << trkCuts_.size()
                                        << ")";
  }
  if (!std::is_sorted(etaBoundaries_.begin(), etaBoundaries_.end())) {
    throw cms::Exception("ConfigError")
        << "EgammaL1TkIsolation: etaBoundaries parameter's entries should be in increasing value";
  }
}

Member Function Documentation

calIsol() [1/2]

std::pair< int, double > EgammaL1TkIsolation::calIsol	(	const reco::TrackBase &	trk,
		const L1TrackCollection &	l1Tks
	)		const

Definition at line 32 of file EgammaL1TkIsolation.cc.

References reco::TrackBase::eta(), reco::TrackBase::phi(), pwdgSkimBPark_cfi::tracks, and reco::TrackBase::vz().

Referenced by calIsolPt(), and EgammaHLTEleL1TrackIsolProducer::produce().

                                                                                                                 {
  return calIsol(trk.eta(), trk.phi(), trk.vz(), tracks);
}

calIsol() [2/2]

std::pair< int, double > EgammaL1TkIsolation::calIsol	(	const double	objEta,
		const double	objPhi,
		const double	objZ,
		const L1TrackCollection &	l1Tks
	)		const

Definition at line 36 of file EgammaL1TkIsolation.cc.

References etaBinNr(), passTrkSel(), pwdgSkimBPark_cfi::tracks, and trkCuts_.

                                                                                           {
  double ptSum = 0.;
  int nrTrks = 0;

  const TrkCuts& cuts = trkCuts_[etaBinNr(objEta)];

  for (const auto& trk : tracks) {
    const float trkPt = trk.momentum().perp();
    if (passTrkSel(trk, trkPt, cuts, objEta, objPhi, objZ)) {
      ptSum += trkPt;
      nrTrks++;
    }
  }
  return {nrTrks, ptSum};
}

calIsolPt()

template<typename... Args>

double EgammaL1TkIsolation::calIsolPt ( Args &&...  args ) const

inline

Definition at line 33 of file EgammaL1TkIsolation.h.

References writedatasetfile::args, and calIsol().

                                         {
    return calIsol(std::forward<Args>(args)...).second;
  }

etaBinNr()

size_t EgammaL1TkIsolation::etaBinNr ( double  eta ) const

private

Definition at line 76 of file EgammaL1TkIsolation.cc.

References funct::abs(), HLT_2023v12_cff::distance, PVValHelper::eta, etaBoundaries_, pfDeepBoostedJetPreprocessParams_cfi::upper_bound, and useAbsEta_.

Referenced by calIsol().

                                                     {
  if (useAbsEta_) {
    eta = std::abs(eta);
  }
  auto res = std::upper_bound(etaBoundaries_.begin(), etaBoundaries_.end(), eta);
  size_t binNr = std::distance(etaBoundaries_.begin(), res);
  return binNr;
}

fillPSetDescription()

void EgammaL1TkIsolation::fillPSetDescription ( edm::ParameterSetDescription &  desc )

static

Definition at line 26 of file EgammaL1TkIsolation.cc.

References submitPVResolutionJobs::desc, EgammaL1TkIsolation::TrkCuts::makePSetDescription(), and trackerHitRTTI::vector.

Referenced by makePSetDescription().

                                                                              {
  desc.add("useAbsEta", true);
  desc.add("etaBoundaries", std::vector{1.5});
  desc.addVPSet("trkCuts", TrkCuts::makePSetDescription(), {edm::ParameterSet(), edm::ParameterSet()});
}

makePSetDescription()

static edm::ParameterSetDescription EgammaL1TkIsolation::makePSetDescription ( )

inlinestatic

Definition at line 18 of file EgammaL1TkIsolation.h.

References submitPVResolutionJobs::desc, and fillPSetDescription().

Referenced by EgammaHLTEleL1TrackIsolProducer::fillDescriptions().

                                                          {
    edm::ParameterSetDescription desc;
    fillPSetDescription(desc);
    return desc;
  }

passTrkSel()

bool EgammaL1TkIsolation::passTrkSel ( const L1Track &  trk, const double  trkPt, const TrkCuts &  cuts, const double  objEta, const double  objPhi, const double  objZ  )

staticprivate

Definition at line 85 of file EgammaL1TkIsolation.cc.

References funct::abs(), reco::deltaR2(), and HLT_2023v12_cff::dEta.

Referenced by calIsol().

                                                        {
  if (trkPt > cuts.minPt && std::abs(objZ - trk.z0()) < cuts.maxDZ) {
    const float trkEta = trk.eta();
    const float dEta = trkEta - objEta;
    const float dR2 = reco::deltaR2(objEta, objPhi, trkEta, trk.phi());
    return dR2 >= cuts.minDR2 && dR2 <= cuts.maxDR2 && std::abs(dEta) >= cuts.minDEta;
  }

  return false;
}

Member Data Documentation

etaBoundaries_

std::vector<double> EgammaL1TkIsolation::etaBoundaries_

private

Definition at line 57 of file EgammaL1TkIsolation.h.

Referenced by EgammaL1TkIsolation(), and etaBinNr().

trkCuts_

std::vector<TrkCuts> EgammaL1TkIsolation::trkCuts_

private

Definition at line 58 of file EgammaL1TkIsolation.h.

Referenced by calIsol(), and EgammaL1TkIsolation().

useAbsEta_

bool EgammaL1TkIsolation::useAbsEta_

private

Definition at line 56 of file EgammaL1TkIsolation.h.

Referenced by etaBinNr().