#include <EgammaTrackExtractor.h>

Inheritance diagram for egammaisolation::EgammaTrackExtractor:

Public Member Functions
reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Candidate &muon) const override

reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &muon) const override

	EgammaTrackExtractor ()

	EgammaTrackExtractor (const edm::ParameterSet &par, edm::ConsumesCollector &&iC)

	EgammaTrackExtractor (const edm::ParameterSet &par, edm::ConsumesCollector &iC)

void	fillVetos (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackCollection &track) override

virtual reco::IsoDeposit::Vetos	vetos (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &track) const

	~EgammaTrackExtractor () override

Public Member Functions inherited from reco::isodeposit::IsoDepositExtractor
virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::CandidateBaseRef &track) const
	make single IsoDeposit based on a CandidateBaseRef as input More...

virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackBaseRef &track) const
	make single IsoDeposit based on trackRef as input More...

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Candidate &track) const

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::CandidateBaseRef &track) const

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &track) const

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackBaseRef &track) const

virtual void	initEvent (const edm::Event &, const edm::EventSetup &)
	perform some constly initialization More...

virtual	~IsoDepositExtractor ()
	Destructor. More...

Private Member Functions
reco::IsoDeposit::Veto	veto (const reco::IsoDeposit::Direction &dir) const

Private Attributes
edm::InputTag	barrelEcalHitsTag_
	"NONE", "BeamSpotFromEvent" More...

int	dzOption

std::string	dzOptionString
	Endcap requirements to determine if isolated for selective filling. More...

edm::InputTag	endcapEcalHitsTag_

double	minCandEt_
	name for deposit More...

std::vector< double >	paramForIsolBarrel_
	min track pt to include into iso deposit More...

std::vector< double >	paramForIsolEndcap_
	Barrel requirements to determine if isolated for selective filling. More...

std::string	theBeamlineOption
	Veto cone angle. More...

edm::EDGetTokenT< reco::BeamSpot >	theBeamSpotToken

double	theChi2Ndof_Max
	trk.numberOfValidHits >= theNHits_Min More...

double	theChi2Prob_Min
	trk.normalizedChi2 < theChi2Ndof_Max More...

std::string	theDepositLabel
	Track Collection Label. More...

double	theDiff_r
	minimum candidate et More...

double	theDiff_z
	transverse distance to vertex More...

double	theDR_Max
	z distance to vertex More...

double	theDR_Veto
	Maximum cone angle for deposits. More...

unsigned int	theNHits_Min
	BeamSpot name. More...

double	thePt_Min
	ChiSquaredProbability(trk.chi2,trk.ndof) > theChi2Prob_Min. More...

edm::EDGetTokenT< edm::View< reco::Track > >	theTrackCollectionToken

Detailed Description

Definition at line 18 of file EgammaTrackExtractor.h.

Constructor & Destructor Documentation

◆ EgammaTrackExtractor() [1/3]

egammaisolation::EgammaTrackExtractor::EgammaTrackExtractor ( )

inline

Definition at line 20 of file EgammaTrackExtractor.h.

20 {};

◆ EgammaTrackExtractor() [2/3]

egammaisolation::EgammaTrackExtractor::EgammaTrackExtractor	(	const edm::ParameterSet &	par,
		edm::ConsumesCollector &&	iC
	)

inline

Definition at line 21 of file EgammaTrackExtractor.h.

21 : EgammaTrackExtractor(par, iC) {}

◆ EgammaTrackExtractor() [3/3]

EgammaTrackExtractor::EgammaTrackExtractor	(	const edm::ParameterSet &	par,
		edm::ConsumesCollector &	iC
	)

Definition at line 28 of file EgammaTrackExtractor.cc.

     : theTrackCollectionToken(iC.consumes<View<Track> >(par.getParameter<edm::InputTag>("inputTrackCollection"))),
       theDepositLabel(par.getUntrackedParameter<std::string>("DepositLabel")),
       theDiff_r(par.getParameter<double>("Diff_r")),
       theDiff_z(par.getParameter<double>("Diff_z")),
       theDR_Max(par.getParameter<double>("DR_Max")),
       theDR_Veto(par.getParameter<double>("DR_Veto")),
       theBeamlineOption(par.getParameter<std::string>("BeamlineOption")),
       theBeamSpotToken(iC.mayConsume<reco::BeamSpot>(par.getParameter<edm::InputTag>("BeamSpotLabel"))),
       theNHits_Min(par.getParameter<unsigned int>("NHits_Min")),
       theChi2Ndof_Max(par.getParameter<double>("Chi2Ndof_Max")),
       theChi2Prob_Min(par.getParameter<double>("Chi2Prob_Min")),
       thePt_Min(par.getParameter<double>("Pt_Min")),
       dzOptionString(par.getParameter<std::string>("dzOption")) {
   if (!dzOptionString.compare("dz"))
     dzOption = EgammaTrackSelector::dz;
   else if (!dzOptionString.compare("vz"))
     dzOption = EgammaTrackSelector::vz;
   else if (!dzOptionString.compare("bs"))
     dzOption = EgammaTrackSelector::bs;
   else if (!dzOptionString.compare("vtx"))
     dzOption = EgammaTrackSelector::vtx;
   else
     dzOption = EgammaTrackSelector::dz;
 }

References egammaisolation::EgammaTrackSelector::bs, egammaisolation::EgammaTrackSelector::dz, dzOption, dzOptionString, egammaisolation::EgammaTrackSelector::vtx, and egammaisolation::EgammaTrackSelector::vz.

◆ ~EgammaTrackExtractor()

egammaisolation::EgammaTrackExtractor::~EgammaTrackExtractor ( )

inlineoverride

Definition at line 24 of file EgammaTrackExtractor.h.

24 {}

Member Function Documentation

◆ deposit() [1/2]

IsoDeposit EgammaTrackExtractor::deposit	(	const edm::Event &	ev,
		const edm::EventSetup &	evSetup,
		const reco::Candidate &	track
	)		const

overridevirtual

make single IsoDeposit based on a candidate as input purely virtual: have to implement in concrete implementations

Reimplemented from reco::isodeposit::IsoDepositExtractor.

Definition at line 68 of file EgammaTrackExtractor.cc.

                                                                         {
   static const std::string metname = "EgammaIsolationAlgos|EgammaTrackExtractor";
  
   reco::isodeposit::Direction candDir;
   double dzCut = 0;
  
   reco::TrackBase::Point beamPoint(0, 0, 0);
   if (theBeamlineOption == "BeamSpotFromEvent") {
     //pick beamSpot
     reco::BeamSpot beamSpot;
     edm::Handle<reco::BeamSpot> beamSpotH;
  
     event.getByToken(theBeamSpotToken, beamSpotH);
  
     if (beamSpotH.isValid()) {
       beamPoint = beamSpotH->position();
     }
   }
  
   Handle<View<Track> > tracksH;
   event.getByToken(theTrackCollectionToken, tracksH);
  
   if (candTk.isElectron()) {
     const reco::GsfElectron* elec = dynamic_cast<const reco::GsfElectron*>(&candTk);
     candDir = reco::isodeposit::Direction(elec->gsfTrack()->eta(), elec->gsfTrack()->phi());
   } else {
     candDir = reco::isodeposit::Direction(candTk.eta(), candTk.phi());
   }
  
   IsoDeposit deposit(candDir);
   deposit.setVeto(veto(candDir));
   deposit.addCandEnergy(candTk.et());
  
   View<Track>::const_iterator itrTr = tracksH->begin();
   View<Track>::const_iterator trEnd = tracksH->end();
   for (itrTr = tracksH->begin(); itrTr != trEnd; ++itrTr) {
     if (candDir.deltaR(reco::isodeposit::Direction(itrTr->eta(), itrTr->phi())) > theDR_Max)
       continue;
  
     if (itrTr->normalizedChi2() > theChi2Ndof_Max)
       continue;
  
     if (itrTr->pt() < thePt_Min)
       continue;
  
     if (theChi2Prob_Min > 0 && ChiSquaredProbability(itrTr->chi2(), itrTr->ndof()) < theChi2Prob_Min)
       continue;
  
     if (theNHits_Min > 0 && itrTr->numberOfValidHits() < theNHits_Min)
       continue;
  
     if (candTk.isElectron()) {
       const reco::GsfElectron* elec = dynamic_cast<const reco::GsfElectron*>(&candTk);
       switch (dzOption) {
         case EgammaTrackSelector::dz:
           dzCut = elec->gsfTrack()->dz() - itrTr->dz();
           break;
         case EgammaTrackSelector::vz:
           dzCut = elec->gsfTrack()->vz() - itrTr->vz();
           break;
         case EgammaTrackSelector::bs:
           dzCut = elec->gsfTrack()->dz(beamPoint) - itrTr->dz(beamPoint);
           break;
         case EgammaTrackSelector::vtx:
           dzCut = itrTr->dz(elec->gsfTrack()->vertex());
           break;
         default:
           dzCut = elec->gsfTrack()->vz() - itrTr->vz();
           break;
       }
     } else {
       switch (dzOption) {
         case EgammaTrackSelector::dz:
           dzCut = (*itrTr).dz() - candTk.vertex().z();
           break;
         case EgammaTrackSelector::vz:
           dzCut = (*itrTr).vz() - candTk.vertex().z();
           break;
         case EgammaTrackSelector::bs:
           dzCut = (*itrTr).dz(beamPoint) - candTk.vertex().z();
           break;
         case EgammaTrackSelector::vtx:
           dzCut = (*itrTr).dz(candTk.vertex());
           break;
         default:
           dzCut = (*itrTr).vz() - candTk.vertex().z();
           break;
       }
     }
  
     if (fabs(dzCut) > theDiff_z)
       continue;
  
     if (fabs(itrTr->dxy(beamPoint)) > theDiff_r)
       continue;
  
     deposit.addDeposit(reco::isodeposit::Direction(itrTr->eta(), itrTr->phi()), itrTr->pt());
   }
  
   return deposit;
 }

References reco::IsoDeposit::addCandEnergy(), reco::IsoDeposit::addDeposit(), pwdgSkimBPark_cfi::beamSpot, egammaisolation::EgammaTrackSelector::bs, ChiSquaredProbability(), reco::isodeposit::Direction::deltaR(), deposit(), egammaisolation::EgammaTrackSelector::dz, TrackSplittingMonitor_cfi::dzCut, dzOption, reco::Candidate::et(), reco::Candidate::eta(), reco::GsfElectron::gsfTrack(), reco::Candidate::isElectron(), edm::HandleBase::isValid(), metname, reco::Candidate::phi(), reco::BeamSpot::position(), reco::IsoDeposit::setVeto(), AlCaHLTBitMon_QueryRunRegistry::string, theBeamlineOption, theBeamSpotToken, theChi2Ndof_Max, theChi2Prob_Min, theDiff_r, theDiff_z, theDR_Max, theNHits_Min, thePt_Min, theTrackCollectionToken, reco::Candidate::vertex(), veto(), egammaisolation::EgammaTrackSelector::vtx, and egammaisolation::EgammaTrackSelector::vz.

◆ deposit() [2/2]

reco::IsoDeposit egammaisolation::EgammaTrackExtractor::deposit	(	const edm::Event &	ev,
		const edm::EventSetup &	evSetup,
		const reco::Track &	track
	)		const

inlineoverridevirtual

make single IsoDeposit based on track as input purely virtual: have to implement in concrete implementations

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 32 of file EgammaTrackExtractor.h.

                                                                    {
       edm::LogWarning("EgammaIsolationAlgos|EgammaTrackExtractor")
           << "This Function is not implemented, bad IsoDeposit Returned";
       return reco::IsoDeposit(reco::isodeposit::Direction(1, 1));
     }

Referenced by deposit().

◆ fillVetos()

void egammaisolation::EgammaTrackExtractor::fillVetos	(	const edm::Event &	ev,
		const edm::EventSetup &	evSetup,
		const reco::TrackCollection &	tracks
	)

inlineoverridevirtual

fill vetoes: to exclude deposits at IsoDeposit creation stage check concrete extractors if it's no-op !

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 26 of file EgammaTrackExtractor.h.

26 {}

◆ veto()

reco::IsoDeposit::Veto EgammaTrackExtractor::veto ( const reco::IsoDeposit::Direction & dir ) const

private

Definition at line 61 of file EgammaTrackExtractor.cc.

                                                                                       {
   reco::IsoDeposit::Veto result;
   result.vetoDir = dir;
   result.dR = theDR_Veto;
   return result;
 }

References DeadROC_duringRun::dir, mps_fire::result, and theDR_Veto.

Referenced by deposit(), and vetos().

◆ vetos()

reco::IsoDeposit::Vetos EgammaTrackExtractor::vetos	(	const edm::Event &	ev,
		const edm::EventSetup &	evSetup,
		const reco::Track &	track
	)		const

virtual

Definition at line 54 of file EgammaTrackExtractor.cc.

                                                                                   {
   reco::isodeposit::Direction dir(track.eta(), track.phi());
   return reco::IsoDeposit::Vetos(1, veto(dir));
 }

References DeadROC_duringRun::dir, HLT_2018_cff::track, and veto().

Member Data Documentation

◆ barrelEcalHitsTag_

edm::InputTag egammaisolation::EgammaTrackExtractor::barrelEcalHitsTag_

private

"NONE", "BeamSpotFromEvent"

Definition at line 57 of file EgammaTrackExtractor.h.

◆ dzOption

int egammaisolation::EgammaTrackExtractor::dzOption

private

Definition at line 67 of file EgammaTrackExtractor.h.

Referenced by deposit(), and EgammaTrackExtractor().

◆ dzOptionString

std::string egammaisolation::EgammaTrackExtractor::dzOptionString

private

Endcap requirements to determine if isolated for selective filling.

Definition at line 66 of file EgammaTrackExtractor.h.

Referenced by EgammaTrackExtractor().

◆ endcapEcalHitsTag_

edm::InputTag egammaisolation::EgammaTrackExtractor::endcapEcalHitsTag_

private

Definition at line 58 of file EgammaTrackExtractor.h.

◆ minCandEt_

double egammaisolation::EgammaTrackExtractor::minCandEt_

private

name for deposit

Definition at line 51 of file EgammaTrackExtractor.h.

◆ paramForIsolBarrel_

std::vector<double> egammaisolation::EgammaTrackExtractor::paramForIsolBarrel_

private

min track pt to include into iso deposit

Definition at line 64 of file EgammaTrackExtractor.h.

◆ paramForIsolEndcap_

std::vector<double> egammaisolation::EgammaTrackExtractor::paramForIsolEndcap_

private

Barrel requirements to determine if isolated for selective filling.

Definition at line 65 of file EgammaTrackExtractor.h.

◆ theBeamlineOption

std::string egammaisolation::EgammaTrackExtractor::theBeamlineOption

private

Veto cone angle.

Definition at line 56 of file EgammaTrackExtractor.h.

Referenced by deposit().

◆ theBeamSpotToken

edm::EDGetTokenT<reco::BeamSpot> egammaisolation::EgammaTrackExtractor::theBeamSpotToken

private

Definition at line 59 of file EgammaTrackExtractor.h.

Referenced by deposit().

◆ theChi2Ndof_Max

double egammaisolation::EgammaTrackExtractor::theChi2Ndof_Max

private

trk.numberOfValidHits >= theNHits_Min

Definition at line 61 of file EgammaTrackExtractor.h.

Referenced by deposit().

◆ theChi2Prob_Min

double egammaisolation::EgammaTrackExtractor::theChi2Prob_Min

private

trk.normalizedChi2 < theChi2Ndof_Max

Definition at line 62 of file EgammaTrackExtractor.h.

Referenced by deposit().

◆ theDepositLabel

std::string egammaisolation::EgammaTrackExtractor::theDepositLabel

private

Track Collection Label.

Definition at line 50 of file EgammaTrackExtractor.h.

◆ theDiff_r

double egammaisolation::EgammaTrackExtractor::theDiff_r

private

minimum candidate et

Definition at line 52 of file EgammaTrackExtractor.h.

Referenced by deposit().

◆ theDiff_z

double egammaisolation::EgammaTrackExtractor::theDiff_z

private

transverse distance to vertex

Definition at line 53 of file EgammaTrackExtractor.h.

Referenced by deposit().

◆ theDR_Max

double egammaisolation::EgammaTrackExtractor::theDR_Max

private

z distance to vertex

Definition at line 54 of file EgammaTrackExtractor.h.

Referenced by deposit().

◆ theDR_Veto

double egammaisolation::EgammaTrackExtractor::theDR_Veto

private

Maximum cone angle for deposits.

Definition at line 55 of file EgammaTrackExtractor.h.

Referenced by veto().

◆ theNHits_Min

unsigned int egammaisolation::EgammaTrackExtractor::theNHits_Min

private

BeamSpot name.

Definition at line 60 of file EgammaTrackExtractor.h.

Referenced by deposit().

◆ thePt_Min

double egammaisolation::EgammaTrackExtractor::thePt_Min

private

ChiSquaredProbability(trk.chi2,trk.ndof) > theChi2Prob_Min.

Definition at line 63 of file EgammaTrackExtractor.h.

Referenced by deposit().

◆ theTrackCollectionToken

edm::EDGetTokenT<edm::View<reco::Track> > egammaisolation::EgammaTrackExtractor::theTrackCollectionToken

private

Definition at line 49 of file EgammaTrackExtractor.h.

Referenced by deposit().

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ EgammaTrackExtractor() [1/3]

◆ EgammaTrackExtractor() [2/3]

◆ EgammaTrackExtractor() [3/3]

◆ ~EgammaTrackExtractor()

Member Function Documentation

◆ deposit() [1/2]

◆ deposit() [2/2]

◆ fillVetos()

◆ veto()

◆ vetos()

Member Data Documentation

◆ barrelEcalHitsTag_

◆ dzOption

◆ dzOptionString

◆ endcapEcalHitsTag_

◆ minCandEt_

◆ paramForIsolBarrel_

◆ paramForIsolEndcap_

◆ theBeamlineOption

◆ theBeamSpotToken

◆ theChi2Ndof_Max

◆ theChi2Prob_Min

◆ theDepositLabel

◆ theDiff_r

◆ theDiff_z

◆ theDR_Max

◆ theDR_Veto

◆ theNHits_Min

◆ thePt_Min

◆ theTrackCollectionToken