egHLT::OffEgSel Class Reference

#include <EgHLTOffEgSel.h>

Public Member Functions
const EgCutValues &	ebCuts () const
const EgCutValues &	eeCuts () const
int	getCutCode (const OffEle &ele, int cutMask=~0x0) const
int	getCutCode (const OffPho &pho, int cutMask=~0x0) const
	OffEgSel ()=default
	OffEgSel (const edm::ParameterSet &config)
bool	passCuts (const OffEle &ele, int cutMask=~0x0) const
bool	passCuts (const OffPho &pho, int cutMask=~0x0) const
void	setEBCuts (const EgCutValues &cuts)
void	setEECuts (const EgCutValues &cuts)
void	setup (const edm::ParameterSet &)
	~OffEgSel ()=default

Static Public Member Functions
static int	getCutCode (const OffEle &ele, const EgCutValues &cuts, int cutMask=~0x0)
static int	getCutCode (const OffPho &pho, const EgCutValues &cuts, int cutMask=~0x0)

Private Attributes
EgCutValues	ebCutValues_
EgCutValues	eeCutValues_

Detailed Description

Definition at line 24 of file EgHLTOffEgSel.h.

Constructor & Destructor Documentation

OffEgSel() [1/2]

egHLT::OffEgSel::OffEgSel ( )

default

OffEgSel() [2/2]

egHLT::OffEgSel::OffEgSel ( const edm::ParameterSet &  config )

inlineexplicit

Definition at line 31 of file EgHLTOffEgSel.h.

References setup().

{ setup(config); }

~OffEgSel()

egHLT::OffEgSel::~OffEgSel ( )

default

Member Function Documentation

ebCuts()

const EgCutValues& egHLT::OffEgSel::ebCuts ( ) const

inline

Definition at line 45 of file EgHLTOffEgSel.h.

References ebCutValues_.

{ return ebCutValues_; }

eeCuts()

const EgCutValues& egHLT::OffEgSel::eeCuts ( ) const

inline

Definition at line 46 of file EgHLTOffEgSel.h.

References eeCutValues_.

{ return eeCutValues_; }

getCutCode() [1/4]

int OffEgSel::getCutCode	(	const OffEle &	ele,
		int	cutMask = ~0x0
	)		const

Definition at line 18 of file EgHLTOffEgSel.cc.

References egHLT::OffEle::detEta(), ebCutValues_, and eeCutValues_.

Referenced by egHLT::OffHelper::fillOffEleVec(), egHLT::OffHelper::fillOffPhoVec(), getCutCode(), and passCuts().

                                                             {
  if (std::fabs(ele.detEta()) < 1.5)
    return getCutCode(ele, ebCutValues_, cutMask);
  else
    return getCutCode(ele, eeCutValues_, cutMask);
}

getCutCode() [2/4]

int OffEgSel::getCutCode ( const OffEle &  ele, const EgCutValues &  cuts, int  cutMask = ~0x0  )

static

Definition at line 25 of file EgHLTOffEgSel.cc.

References egHLT::OffEle::caloEnergy(), egHLT::EgCutCodes::CRACK, egHLT::EgCutCodes::CTFTRACK, egHLT::OffEle::ctfTrack(), egHLT::OffEle::ctfTrkHitsFound(), egHLT::OffEle::ctfTrkHitsLost(), egHLT::OffEle::ctfTrkInnerRadius(), egHLT::OffEle::ctfTrkOuterRadius(), egHLT::EgCutCodes::DETAIN, egHLT::OffEle::dEtaIn(), egHLT::EgCutCodes::DETETA, egHLT::EgCutCodes::DPHIIN, egHLT::OffEle::dPhiIn(), egHLT::EgCutCodes::ET, egHLT::OffEle::et(), egHLT::OffEle::etaSC(), egHLT::EgCutCodes::HADEM, egHLT::EgCutCodes::HLTDETAIN, egHLT::OffEle::hltDEtaIn(), egHLT::EgCutCodes::HLTDPHIIN, egHLT::OffEle::hltDPhiIn(), egHLT::EgCutCodes::HLTINVEINVP, egHLT::OffEle::hltInvEInvP(), egHLT::EgCutCodes::HLTISOLEM, egHLT::OffEle::hltIsolEm(), egHLT::EgCutCodes::HLTISOLHAD, egHLT::OffEle::hltIsolHad(), egHLT::EgCutCodes::HLTISOLTRKSELE, egHLT::OffEle::hltIsolTrksEle(), egHLT::EgCutCodes::HLTISOLTRKSPHO, egHLT::OffEle::hltIsolTrksPho(), egHLT::OffEle::hOverE(), egHLT::EgCutCodes::INVEINVP, egHLT::OffEle::invEInvP(), egHLT::OffEle::isGap(), egHLT::EgCutCodes::ISOLEM, egHLT::OffEle::isolEm(), egHLT::EgCutCodes::ISOLHAD, egHLT::OffEle::isolHad(), egHLT::EgCutCodes::ISOLPTTRKS, egHLT::OffEle::isolPtTrks(), egHLT::EgCutCodes::MAXR9, egHLT::EgCutCodes::MINR9, egHLT::OffEle::r9(), egHLT::EgCutCodes::SIGMAETAETA, egHLT::OffEle::sigmaEtaEta(), egHLT::EgCutCodes::SIGMAIETAIETA, egHLT::OffEle::sigmaIEtaIEta(), and egHLT::OffEle::validCTFTrack().

                                                                                {
  int cutCode = 0x0;
  //kinematic cuts
  if (ele.et() < cuts.minEt)
    cutCode |= EgCutCodes::ET;
  if (std::fabs(ele.etaSC()) < cuts.minEta || std::fabs(ele.etaSC()) > cuts.maxEta)
    cutCode |= EgCutCodes::DETETA;
  if (ele.isGap())
    cutCode |= EgCutCodes::CRACK;
  //track cuts
  if (std::fabs(ele.dEtaIn()) > cuts.maxDEtaIn)
    cutCode |= EgCutCodes::DETAIN;
  if (std::fabs(ele.dPhiIn()) > cuts.maxDPhiIn)
    cutCode |= EgCutCodes::DPHIIN;
  if (ele.invEInvP() > cuts.maxInvEInvP)
    cutCode |= EgCutCodes::INVEINVP;
  //supercluster cuts
  if (ele.hOverE() > cuts.maxHadem && ele.hOverE() * ele.caloEnergy() > cuts.maxHadEnergy)
    cutCode |= EgCutCodes::HADEM;
  if (ele.sigmaIEtaIEta() > cuts.maxSigmaIEtaIEta)
    cutCode |= EgCutCodes::SIGMAIETAIETA;
  if (ele.sigmaEtaEta() > cuts.maxSigmaEtaEta)
    cutCode |= EgCutCodes::SIGMAETAETA;
  //---Morse-------
  //if(ele.r9()<cuts.minR9) cutCode |= EgCutCodes::R9;
  if (ele.r9() < cuts.minR9)
    cutCode |= EgCutCodes::MINR9;
  if (ele.r9() > cuts.maxR9)
    cutCode |= EgCutCodes::MAXR9;
  //----------------

  //std isolation cuts
  if (ele.isolEm() > (cuts.isolEmConstTerm +
                      cuts.isolEmGradTerm * (ele.et() < cuts.isolEmGradStart ? 0. : (ele.et() - cuts.isolEmGradStart))))
    cutCode |= EgCutCodes::ISOLEM;
  if (ele.isolHad() >
      (cuts.isolHadConstTerm +
       cuts.isolHadGradTerm * (ele.et() < cuts.isolHadGradStart ? 0. : (ele.et() - cuts.isolHadGradStart))))
    cutCode |= EgCutCodes::ISOLHAD;
  if (ele.isolPtTrks() >
      (cuts.isolPtTrksConstTerm +
       cuts.isolPtTrksGradTerm * (ele.et() < cuts.isolPtTrksGradStart ? 0. : (ele.et() - cuts.isolPtTrksGradStart))))
    cutCode |= EgCutCodes::ISOLPTTRKS;
  //ele Nr trks not defined, assume it passes
  //hlt isolation cuts
  if (ele.et() <=
      0.) {  //even it if et<=0, we give it a shot at passing isolation. Note this should be an impossible case
    if (ele.hltIsolTrksEle() > cuts.maxHLTIsolTrksEle)
      cutCode |= EgCutCodes::HLTISOLTRKSELE;
    if (ele.hltIsolTrksPho() > cuts.maxHLTIsolTrksPho)
      cutCode |= EgCutCodes::HLTISOLTRKSPHO;
    if (ele.hltIsolHad() > cuts.maxHLTIsolHad)
      cutCode |= EgCutCodes::HLTISOLHAD;
    if (ele.hltIsolEm() > cuts.maxHLTIsolEm)
      cutCode |= EgCutCodes::HLTISOLEM;
  } else {
    if (ele.hltIsolTrksEle() > cuts.maxHLTIsolTrksEle &&
        ele.hltIsolTrksEle() / ele.et() > cuts.maxHLTIsolTrksEleOverPt &&
        ele.hltIsolTrksEle() / ele.et() / ele.et() > cuts.maxHLTIsolTrksEleOverPt2)
      cutCode |= EgCutCodes::HLTISOLTRKSELE;
    if (ele.hltIsolTrksPho() > cuts.maxHLTIsolTrksPho &&
        ele.hltIsolTrksPho() / ele.et() > cuts.maxHLTIsolTrksPhoOverPt &&
        ele.hltIsolTrksPho() / ele.et() / ele.et() > cuts.maxHLTIsolTrksPhoOverPt2)
      cutCode |= EgCutCodes::HLTISOLTRKSPHO;
    if (ele.hltIsolHad() > cuts.maxHLTIsolHad && ele.hltIsolHad() / ele.et() > cuts.maxHLTIsolHadOverEt &&
        ele.hltIsolHad() / ele.et() / ele.et() > cuts.maxHLTIsolHadOverEt2)
      cutCode |= EgCutCodes::HLTISOLHAD;
    if (ele.hltIsolEm() > cuts.maxHLTIsolEm && ele.hltIsolEm() / ele.et() > cuts.maxHLTIsolEmOverEt &&
        ele.hltIsolEm() / ele.et() / ele.et() > cuts.maxHLTIsolEmOverEt2)
      cutCode |= EgCutCodes::HLTISOLEM;
  }

  //cuts on CTF track, HLT tracking doesnt handle poor quaility tracks
  if (ele.validCTFTrack()) {
    if (!(ele.ctfTrkOuterRadius() >= cuts.minCTFTrkOuterRadius &&  //note I'm NOTing the result of the AND
          ele.ctfTrkInnerRadius() <= cuts.maxCTFTrkInnerRadius &&

          ele.ctfTrkHitsFound() >= cuts.minNrCTFTrkHits && ele.ctfTrkHitsLost() <= cuts.maxNrCTFTrkHitsLost))
      cutCode |= EgCutCodes::CTFTRACK;  //the next line can also set this bit
    if (cuts.requirePixelHitsIfOuterInOuter) {
      DetId innerDetId(ele.ctfTrack()->extra()->innerDetId());
      DetId outerDetId(ele.ctfTrack()->extra()->outerDetId());

      if (outerDetId.subdetId() >= 5 && innerDetId.subdetId() >= 3)
        cutCode |= EgCutCodes::CTFTRACK;  //1,2 = pixel, 3,4,5,6 sistrip
    }
    // std::cout <<"eta "<<ele.detEta()<<" max inner "<<cuts.maxCTFTrkInnerRadius<<" inner "<<ele.ctfTrkInnerRadius()<<std::endl;
  } else
    cutCode |= EgCutCodes::CTFTRACK;

  if (std::fabs(ele.hltDEtaIn()) > cuts.maxHLTDEtaIn)
    cutCode |= EgCutCodes::HLTDETAIN;
  if (std::fabs(ele.hltDPhiIn()) > cuts.maxHLTDPhiIn)
    cutCode |= EgCutCodes::HLTDPHIIN;
  if (std::fabs(ele.hltInvEInvP()) > cuts.maxHLTInvEInvP)
    cutCode |= EgCutCodes::HLTINVEINVP;

  return (cutCode & cuts.cutMask & cutMask);
}

getCutCode() [3/4]

int OffEgSel::getCutCode	(	const OffPho &	pho,
		int	cutMask = ~0x0
	)		const

Definition at line 125 of file EgHLTOffEgSel.cc.

References egHLT::OffPho::detEta(), ebCutValues_, eeCutValues_, and getCutCode().

                                                             {
  if (std::fabs(pho.detEta()) < 1.5)
    return getCutCode(pho, ebCutValues_, cutMask);
  else
    return getCutCode(pho, eeCutValues_, cutMask);
}

getCutCode() [4/4]

int OffEgSel::getCutCode ( const OffPho &  pho, const EgCutValues &  cuts, int  cutMask = ~0x0  )

static

Definition at line 133 of file EgHLTOffEgSel.cc.

References egHLT::EgCutCodes::CRACK, egHLT::EgCutCodes::CTFTRACK, egHLT::EgCutCodes::DETAIN, egHLT::EgCutCodes::DETETA, egHLT::EgCutCodes::DPHIIN, egHLT::OffPho::energy(), egHLT::EgCutCodes::ET, egHLT::OffPho::et(), egHLT::OffPho::etaSC(), egHLT::EgCutCodes::HADEM, egHLT::EgCutCodes::HLTDETAIN, egHLT::EgCutCodes::HLTDPHIIN, egHLT::EgCutCodes::HLTINVEINVP, egHLT::EgCutCodes::HLTISOLEM, egHLT::OffPho::hltIsolEm(), egHLT::EgCutCodes::HLTISOLHAD, egHLT::OffPho::hltIsolHad(), egHLT::OffPho::hltIsolTrks(), egHLT::EgCutCodes::HLTISOLTRKSELE, egHLT::EgCutCodes::HLTISOLTRKSPHO, egHLT::OffPho::hOverE(), egHLT::EgCutCodes::INVEINVP, egHLT::OffPho::isGap(), egHLT::EgCutCodes::ISOLEM, egHLT::OffPho::isolEm(), egHLT::EgCutCodes::ISOLHAD, egHLT::OffPho::isolHad(), egHLT::EgCutCodes::ISOLNRTRKS, egHLT::OffPho::isolNrTrks(), egHLT::EgCutCodes::ISOLPTTRKS, egHLT::OffPho::isolPtTrks(), egHLT::EgCutCodes::MAXR9, egHLT::EgCutCodes::MINR9, egHLT::OffPho::r9(), egHLT::EgCutCodes::SIGMAETAETA, egHLT::OffPho::sigmaEtaEta(), egHLT::EgCutCodes::SIGMAIETAIETA, and egHLT::OffPho::sigmaIEtaIEta().

                                                                                {
  int cutCode = 0x0;
  //kinematic cuts
  if (pho.et() < cuts.minEt)
    cutCode |= EgCutCodes::ET;
  if (std::fabs(pho.etaSC()) < cuts.minEta || std::fabs(pho.etaSC()) > cuts.maxEta)
    cutCode |= EgCutCodes::DETETA;
  if (pho.isGap())
    cutCode |= EgCutCodes::CRACK;
  //track cuts (all fail)
  cutCode |= EgCutCodes::DETAIN;
  cutCode |= EgCutCodes::DPHIIN;
  cutCode |= EgCutCodes::INVEINVP;
  //supercluster cuts
  if (pho.hOverE() > cuts.maxHadem && pho.hOverE() * pho.energy() > cuts.maxHadEnergy)
    cutCode |= EgCutCodes::HADEM;
  if (pho.sigmaIEtaIEta() > cuts.maxSigmaIEtaIEta)
    cutCode |= EgCutCodes::SIGMAIETAIETA;
  if (pho.sigmaEtaEta() > cuts.maxSigmaEtaEta)
    cutCode |= EgCutCodes::SIGMAETAETA;
  //----Morse--------
  //if(pho.r9()<cuts.minR9) cutCode |= EgCutCodes::R9;
  if (pho.r9() < cuts.minR9)
    cutCode |= EgCutCodes::MINR9;
  if (pho.r9() > cuts.maxR9)
    cutCode |= EgCutCodes::MAXR9;
  //----------------
  //std isolation cuts
  if (pho.isolEm() > (cuts.isolEmConstTerm +
                      cuts.isolEmGradTerm * (pho.et() < cuts.isolEmGradStart ? 0. : (pho.et() - cuts.isolEmGradStart))))
    cutCode |= EgCutCodes::ISOLEM;
  if (pho.isolHad() >
      (cuts.isolHadConstTerm +
       cuts.isolHadGradTerm * (pho.et() < cuts.isolHadGradStart ? 0. : (pho.et() - cuts.isolHadGradStart))))
    cutCode |= EgCutCodes::ISOLHAD;
  if (pho.isolPtTrks() >
      (cuts.isolPtTrksConstTerm +
       cuts.isolPtTrksGradTerm * (pho.et() < cuts.isolPtTrksGradStart ? 0. : (pho.et() - cuts.isolPtTrksGradStart))))
    cutCode |= EgCutCodes::ISOLPTTRKS;
  if (pho.isolNrTrks() > cuts.isolNrTrksConstTerm)
    cutCode |= EgCutCodes::ISOLNRTRKS;

  //hlt isolation cuts
  cutCode |= EgCutCodes::HLTISOLTRKSELE;  //automatically fails ele track isolation
  if (pho.et() <=
      0.) {  //even it if et<=0, we give it a shot at passing isolation. Note this should be an impossible case
    if (pho.hltIsolTrks() > cuts.maxHLTIsolTrksPho)
      cutCode |= EgCutCodes::HLTISOLTRKSPHO;
    if (pho.hltIsolHad() > cuts.maxHLTIsolHad)
      cutCode |= EgCutCodes::HLTISOLHAD;
    if (pho.hltIsolEm() > cuts.maxHLTIsolEm)
      cutCode |= EgCutCodes::HLTISOLEM;
  } else {
    if (pho.hltIsolTrks() > cuts.maxHLTIsolTrksPho && pho.hltIsolTrks() / pho.et() > cuts.maxHLTIsolTrksPhoOverPt &&
        pho.hltIsolTrks() / pho.et() / pho.et() > cuts.maxHLTIsolTrksPhoOverPt2)
      cutCode |= EgCutCodes::HLTISOLTRKSPHO;
    if (pho.hltIsolHad() > cuts.maxHLTIsolHad && pho.hltIsolHad() / pho.et() > cuts.maxHLTIsolHadOverEt &&
        pho.hltIsolHad() / pho.et() / pho.et() > cuts.maxHLTIsolHadOverEt2)
      cutCode |= EgCutCodes::HLTISOLHAD;
    if (pho.hltIsolEm() > cuts.maxHLTIsolEm && pho.hltIsolEm() / pho.et() > cuts.maxHLTIsolEmOverEt &&
        pho.hltIsolEm() / pho.et() / pho.et() > cuts.maxHLTIsolEmOverEt2)
      cutCode |= EgCutCodes::HLTISOLEM;
  }

  //track cuts, photon will automatically fail them (for now)
  cutCode |= EgCutCodes::CTFTRACK;
  cutCode |= EgCutCodes::HLTDETAIN;
  cutCode |= EgCutCodes::HLTDPHIIN;
  cutCode |= EgCutCodes::HLTINVEINVP;

  return (cutCode & cuts.cutMask & cutMask);
}

passCuts() [1/2]

bool egHLT::OffEgSel::passCuts ( const OffEle &  ele, int  cutMask = ~0x0  ) const

inline

Definition at line 34 of file EgHLTOffEgSel.h.

References getCutCode().

{ return getCutCode(ele, cutMask) == 0x0; }

passCuts() [2/2]

bool egHLT::OffEgSel::passCuts ( const OffPho &  pho, int  cutMask = ~0x0  ) const

inline

Definition at line 38 of file EgHLTOffEgSel.h.

References getCutCode().

{ return getCutCode(pho, cutMask) == 0x0; }

setEBCuts()

void egHLT::OffEgSel::setEBCuts ( const EgCutValues &  cuts )

inline

Definition at line 42 of file EgHLTOffEgSel.h.

References ebCutValues_.

{ ebCutValues_ = cuts; }

setEECuts()

void egHLT::OffEgSel::setEECuts ( const EgCutValues &  cuts )

inline

Definition at line 43 of file EgHLTOffEgSel.h.

References eeCutValues_.

{ eeCutValues_ = cuts; }

setup()

void OffEgSel::setup

(

const edm::ParameterSet &

iConfig

)

Definition at line 13 of file EgHLTOffEgSel.cc.

References ebCutValues_, eeCutValues_, edm::ParameterSet::getParameter(), and egHLT::EgCutValues::setup().

Referenced by OffEgSel(), and egHLT::OffHelper::setup().

                                                   {
  ebCutValues_.setup(iConfig.getParameter<edm::ParameterSet>("barrel"));
  eeCutValues_.setup(iConfig.getParameter<edm::ParameterSet>("endcap"));
}

Member Data Documentation

ebCutValues_

EgCutValues egHLT::OffEgSel::ebCutValues_

private

Definition at line 26 of file EgHLTOffEgSel.h.

Referenced by ebCuts(), getCutCode(), setEBCuts(), and setup().

eeCutValues_

EgCutValues egHLT::OffEgSel::eeCutValues_

private

Definition at line 27 of file EgHLTOffEgSel.h.

Referenced by eeCuts(), getCutCode(), setEECuts(), and setup().