EleTkIsolFromCands Class Reference

#include <EleTkIsolFromCands.h>

Classes
struct	TrkCuts

Public Member Functions
std::pair< int, double >	calIsol (const reco::TrackBase &trk, const pat::PackedCandidateCollection &cands, const edm::View< reco::GsfElectron > &eles)
std::pair< int, double >	calIsol (const double eleEta, const double elePhi, const double eleVZ, const pat::PackedCandidateCollection &cands, const edm::View< reco::GsfElectron > &eles)
double	calIsolPt (const reco::TrackBase &trk, const pat::PackedCandidateCollection &cands, const edm::View< reco::GsfElectron > &eles)
double	calIsolPt (const double eleEta, const double elePhi, const double eleVZ, const pat::PackedCandidateCollection &cands, const edm::View< reco::GsfElectron > &eles)
	EleTkIsolFromCands (const edm::ParameterSet &para)
	EleTkIsolFromCands (const EleTkIsolFromCands &)=default
EleTkIsolFromCands &	operator= (const EleTkIsolFromCands &)=default
	~EleTkIsolFromCands ()=default

Static Public Member Functions
static bool	passTrkSel (const reco::Track &trk, const double trkPt, const TrkCuts &cuts, const double eleEta, const double elePhi, const double eleVZ)
static edm::ParameterSetDescription	pSetDescript ()

Private Member Functions
double	getTrkPt (const reco::TrackBase &trk, const edm::View< reco::GsfElectron > &eles)

Static Private Member Functions
static bool	passAlgo (const reco::TrackBase &trk, const std::vector< reco::TrackBase::TrackAlgorithm > &algosToRej)
static bool	passQual (const reco::TrackBase &trk, const std::vector< reco::TrackBase::TrackQuality > &quals)

Private Attributes
TrkCuts	barrelCuts_
TrkCuts	endcapCuts_

Detailed Description

Definition at line 11 of file EleTkIsolFromCands.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 47 of file EleTkIsolFromCands.cc.

                                                                 :
  barrelCuts_(para.getParameter<edm::ParameterSet>("barrelCuts")),
  endcapCuts_(para.getParameter<edm::ParameterSet>("endcapCuts"))          
{
  

}

EleTkIsolFromCands::EleTkIsolFromCands ( const EleTkIsolFromCands &  )

default

EleTkIsolFromCands::~EleTkIsolFromCands ( )

default

Member Function Documentation

std::pair< int, double > EleTkIsolFromCands::calIsol	(	const reco::TrackBase &	trk,
		const pat::PackedCandidateCollection &	cands,
		const edm::View< reco::GsfElectron > &	eles
	)

Definition at line 64 of file EleTkIsolFromCands.cc.

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

Referenced by calIsolPt().

{
  return calIsol(eleTrk.eta(),eleTrk.phi(),eleTrk.vz(),cands,eles);
}

std::pair< int, double > EleTkIsolFromCands::calIsol	(	const double	eleEta,
		const double	elePhi,
		const double	eleVZ,
		const pat::PackedCandidateCollection &	cands,
		const edm::View< reco::GsfElectron > &	eles
	)

Definition at line 72 of file EleTkIsolFromCands.cc.

References funct::abs(), barrelCuts_, particleFlowClusterECALTimeSelected_cfi::cuts, endcapCuts_, getTrkPt(), passTrkSel(), and reco::TrackBase::pt().

{

  double ptSum=0.;
  int nrTrks=0;

  const TrkCuts& cuts = std::abs(eleEta)<1.5 ? barrelCuts_ : endcapCuts_;
  
  for(auto& cand  : cands){
    if(cand.charge()!=0){
      const reco::Track& trk = cand.pseudoTrack(); 
      double trkPt = std::abs(cand.pdgId())!=11 ? trk.pt() : getTrkPt(trk,eles);
      if(passTrkSel(trk,trkPt,cuts,eleEta,elePhi,eleVZ)){   
    ptSum+=trkPt;
    nrTrks++;
      }
    }
  }
  return {nrTrks,ptSum};    
}

double EleTkIsolFromCands::calIsolPt ( const reco::TrackBase &  trk, const pat::PackedCandidateCollection &  cands, const edm::View< reco::GsfElectron > &  eles  )

inline

Definition at line 45 of file EleTkIsolFromCands.h.

References calIsol().

Referenced by ElectronHEEPIDValueMapProducer::calTrkIso().

                                                {
    return calIsol(trk,cands,eles).second;
  }

double EleTkIsolFromCands::calIsolPt ( const double  eleEta, const double  elePhi, const double  eleVZ, const pat::PackedCandidateCollection &  cands, const edm::View< reco::GsfElectron > &  eles  )

inline

Definition at line 50 of file EleTkIsolFromCands.h.

References calIsol(), particleFlowClusterECALTimeSelected_cfi::cuts, getTrkPt(), passAlgo(), passQual(), and passTrkSel().

                                                {
    return calIsol(eleEta,elePhi,eleVZ,cands,eles).second;
  }

double EleTkIsolFromCands::getTrkPt ( const reco::TrackBase &  trk, const edm::View< reco::GsfElectron > &  eles  )

private

Definition at line 143 of file EleTkIsolFromCands.cc.

References funct::abs(), reco::deltaPhi(), reco::TrackBase::eta(), match(), reco::TrackBase::phi(), and reco::TrackBase::pt().

Referenced by calIsol(), calIsolPt(), and passAlgo().

{
  //note, the trk.eta(),trk.phi() should be identical to the gsf track eta,phi
  //although this may not be the case due to roundings after packing
  auto match=[](const reco::TrackBase& trk,const reco::GsfElectron& ele){
    return std::abs(trk.eta()-ele.gsfTrack()->eta())<0.001 &&
    reco::deltaPhi(trk.phi(),ele.gsfTrack()->phi())<0.001;// && 
  };
  for(auto& ele : eles){
    if(ele.gsfTrack().isNonnull()){
      if(match(trk,ele)){
    return ele.gsfTrack()->pt();
      }
    }
  }
  return trk.pt();

}

EleTkIsolFromCands& EleTkIsolFromCands::operator= ( const EleTkIsolFromCands &  )

default

bool EleTkIsolFromCands::passAlgo ( const reco::TrackBase &  trk, const std::vector< reco::TrackBase::TrackAlgorithm > &  algosToRej  )

staticprivate

Definition at line 133 of file EleTkIsolFromCands.cc.

References reco::TrackBase::algo(), and getTrkPt().

Referenced by calIsolPt(), passQual(), and passTrkSel().

{
  return algosToRej.empty() || !std::binary_search(algosToRej.begin(),algosToRej.end(),trk.algo());
}

bool EleTkIsolFromCands::passQual ( const reco::TrackBase &  trk, const std::vector< reco::TrackBase::TrackQuality > &  quals  )

staticprivate

Definition at line 120 of file EleTkIsolFromCands.cc.

References passAlgo(), and reco::TrackBase::quality().

Referenced by calIsolPt(), and passTrkSel().

{
  if(quals.empty()) return true;

  for(auto qual : quals) {
    if(trk.quality(qual)) return true;
  }

  return false;  
}

bool EleTkIsolFromCands::passTrkSel ( const reco::Track &  trk, const double  trkPt, const TrkCuts &  cuts, const double  eleEta, const double  elePhi, const double  eleVZ  )

static

Definition at line 97 of file EleTkIsolFromCands.cc.

References funct::abs(), EleTkIsolFromCands::TrkCuts::algosToReject, EleTkIsolFromCands::TrkCuts::allowedQualities, reco::deltaR2(), reco::TrackBase::eta(), reco::TrackBase::hitPattern(), EleTkIsolFromCands::TrkCuts::maxDPtPt, EleTkIsolFromCands::TrkCuts::maxDR2, EleTkIsolFromCands::TrkCuts::maxDZ, EleTkIsolFromCands::TrkCuts::minDEta, EleTkIsolFromCands::TrkCuts::minDR2, EleTkIsolFromCands::TrkCuts::minHits, EleTkIsolFromCands::TrkCuts::minPixelHits, EleTkIsolFromCands::TrkCuts::minPt, reco::HitPattern::numberOfValidHits(), reco::HitPattern::numberOfValidPixelHits(), passAlgo(), passQual(), reco::TrackBase::phi(), reco::TrackBase::ptError(), and reco::TrackBase::vz().

Referenced by calIsol(), and calIsolPt().

{
  const float dR2 = reco::deltaR2(eleEta,elePhi,trk.eta(),trk.phi());
  const float dEta = trk.eta()-eleEta;
  const float dZ = eleVZ - trk.vz();

  return dR2>=cuts.minDR2 && dR2<=cuts.maxDR2 && 
    std::abs(dEta)>=cuts.minDEta && 
    std::abs(dZ)<cuts.maxDZ &&
    trk.hitPattern().numberOfValidHits() >= cuts.minHits &&
    trk.hitPattern().numberOfValidPixelHits() >=cuts.minPixelHits &&
    (trk.ptError()/trkPt < cuts.maxDPtPt || cuts.maxDPtPt<0) && 
    passQual(trk,cuts.allowedQualities) &&
    passAlgo(trk,cuts.algosToReject) &&
    trkPt > cuts.minPt;
}

edm::ParameterSetDescription EleTkIsolFromCands::pSetDescript ( )

static

Definition at line 55 of file EleTkIsolFromCands.cc.

References edm::ParameterSetDescription::add(), and EleTkIsolFromCands::TrkCuts::pSetDescript().

Referenced by ElectronHEEPIDValueMapProducer::fillDescriptions().

{
  edm::ParameterSetDescription desc;
  desc.add("barrelCuts",TrkCuts::pSetDescript());
  desc.add("endcapCuts",TrkCuts::pSetDescript());
  return desc;
}

Member Data Documentation

TrkCuts EleTkIsolFromCands::barrelCuts_

private

Definition at line 29 of file EleTkIsolFromCands.h.

Referenced by calIsol().

TrkCuts EleTkIsolFromCands::endcapCuts_

private

Definition at line 29 of file EleTkIsolFromCands.h.

Referenced by calIsol().