#include <PixelTrackExtractor.h>

Inheritance diagram for muonisolation::PixelTrackExtractor:

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

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

	PixelTrackExtractor ()

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

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

	~PixelTrackExtractor () override

Public Member Functions inherited from reco::isodeposit::IsoDepositExtractor
virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Candidate &track) const

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::Direction	directionAtPresetRadius (const reco::Track &tk, double bz) const

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

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

edm::EDGetTokenT< reco::BeamSpot >	theBeamSpotToken
	"NONE", "BeamSpotFromEvent" More...

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

double	theChi2Prob_Min
	trk.normalizedChi2 < theChi2Ndof_Max More...

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

double	theDiff_r
	name for deposit 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...

double	theDR_VetoPt
	.. it is above this threshold More...

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	theFieldToken

unsigned int	theNHits_Min
	BeamSpot name. More...

bool	thePropagateTracksToRadius
	min track pt to include into iso deposit More...

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

double	thePtVeto_Min
	will veto leading track if More...

double	theReferenceRadius
	If set to true will compare track eta-phi at ... More...

edm::EDGetTokenT< reco::TrackCollection >	theTrackCollectionToken

bool	theVetoLeadingTrack
	... this radius More...

Detailed Description

Definition at line 20 of file PixelTrackExtractor.h.

Constructor & Destructor Documentation

◆ PixelTrackExtractor() [1/2]

muonisolation::PixelTrackExtractor::PixelTrackExtractor ( )

inline

Definition at line 22 of file PixelTrackExtractor.h.

22 {};

◆ PixelTrackExtractor() [2/2]

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

Definition at line 19 of file PixelTrackExtractor.cc.

     : theFieldToken(iC.esConsumes()),
       theTrackCollectionToken(iC.consumes<TrackCollection>(par.getParameter<edm::InputTag>("inputTrackCollection"))),
       theDepositLabel(par.getUntrackedParameter<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<string>("BeamlineOption")),
       theBeamSpotToken(iC.mayConsume<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")),
       thePropagateTracksToRadius(par.getParameter<bool>("PropagateTracksToRadius")),
       theReferenceRadius(par.getParameter<double>("ReferenceRadius")),
       theVetoLeadingTrack(par.getParameter<bool>("VetoLeadingTrack")),  
       thePtVeto_Min(par.getParameter<double>("PtVeto_Min")),            
       theDR_VetoPt(par.getParameter<double>("DR_VetoPt"))               
 {}

◆ ~PixelTrackExtractor()

muonisolation::PixelTrackExtractor::~PixelTrackExtractor ( )

inlineoverride

Definition at line 25 of file PixelTrackExtractor.h.

25 {}

Member Function Documentation

◆ deposit()

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

overridevirtual

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

Note, this can be reset below if theVetoLeadingTrack is set and conditions are met

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 70 of file PixelTrackExtractor.cc.

                                                                                                                  {
   static const std::string metname = "MuonIsolation|PixelTrackExtractor";
  
   auto const& bField = eventSetup.getData(theFieldToken);
   double bz = bField.inInverseGeV(GlobalPoint(0., 0., 0.)).z();
  
   Direction muonDir(directionAtPresetRadius(muon, bz));
   IsoDeposit deposit(muonDir);
   deposit.setVeto(veto(muonDir));
  
   deposit.addCandEnergy(muon.pt());
  
   Handle<TrackCollection> tracksH;
   event.getByToken(theTrackCollectionToken, tracksH);
   //  const TrackCollection tracks = *(tracksH.product());
   LogTrace(metname) << "***** TRACK COLLECTION SIZE: " << tracksH->size();
  
   double vtx_z = muon.vz();
   LogTrace(metname) << "***** Muon vz: " << vtx_z;
   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();
       LogTrace(metname) << "Extracted beam point at " << beamPoint << std::endl;
     }
   }
  
   LogTrace(metname) << "Using beam point at " << beamPoint << std::endl;
  
   TrackSelector::Parameters pars(
       TrackSelector::Range(vtx_z - theDiff_z, vtx_z + theDiff_z), theDiff_r, muonDir, theDR_Max, beamPoint);
  
   pars.nHitsMin = theNHits_Min;
   pars.chi2NdofMax = theChi2Ndof_Max;
   pars.chi2ProbMin = theChi2Prob_Min;
   pars.ptMin = thePt_Min;
  
   TrackSelector selection(pars);
   TrackSelector::result_type sel_tracks = selection(*tracksH);
   LogTrace(metname) << "all tracks: " << tracksH->size() << " selected: " << sel_tracks.size();
  
   double maxPt = -1;
   Direction maxPtDir;
   TrackSelector::result_type::const_iterator tkI = sel_tracks.begin();
   for (; tkI != sel_tracks.end(); ++tkI) {
     const reco::Track* tk = *tkI;
     LogTrace(metname) << "This track has: pt= " << tk->pt() << ", eta= " << tk->eta() << ", phi= " << tk->phi();
     Direction dirTrk(directionAtPresetRadius(*tk, bz));
     deposit.addDeposit(dirTrk, tk->pt());
     double tkDr = (muonDir - dirTrk).deltaR;
     double tkPt = tk->pt();
     if (theVetoLeadingTrack && tkPt > thePtVeto_Min && tkDr < theDR_VetoPt && maxPt < tkPt) {
       maxPt = tkPt;
       maxPtDir = dirTrk;
     }
   }
   if (maxPt > 0) {
     deposit.setVeto(veto(maxPtDir));
     LogTrace(metname) << " Set track veto the leading track with pt " << maxPt << " in direction  (eta,phi) "
                       << maxPtDir.eta() << ", " << maxPtDir.phi();
   }
  
   return deposit;
 }

References reco::IsoDeposit::addCandEnergy(), reco::IsoDeposit::addDeposit(), pwdgSkimBPark_cfi::beamSpot, Calorimetry_cff::bField, PbPb_ZMuSkimMuonDPG_cff::deltaR, directionAtPresetRadius(), reco::isodeposit::Direction::eta(), reco::TrackBase::eta(), edm::EventSetup::getData(), edm::HandleBase::isValid(), LogTrace, L1TrackerEtMissProducer_cfi::maxPt, metname, reco::isodeposit::Direction::phi(), reco::TrackBase::phi(), reco::BeamSpot::position(), reco::TrackBase::pt(), corrVsCorr::selection, reco::IsoDeposit::setVeto(), AlCaHLTBitMon_QueryRunRegistry::string, theBeamlineOption, theBeamSpotToken, theChi2Ndof_Max, theChi2Prob_Min, theDiff_r, theDiff_z, theDR_Max, theDR_VetoPt, theFieldToken, theNHits_Min, thePt_Min, thePtVeto_Min, theTrackCollectionToken, theVetoLeadingTrack, and veto().

◆ directionAtPresetRadius()

Direction PixelTrackExtractor::directionAtPresetRadius	(	const reco::Track &	tk,
		double	bz
	)		const

private

Definition at line 54 of file PixelTrackExtractor.cc.

                                                                                        {
   if (!thePropagateTracksToRadius) {
     return Direction(tk.eta(), tk.phi());
   }
  
   // this should represent a cylinder in global frame at R=refRadius cm, roughly where mid-layer of pixels is
   double psRadius = theReferenceRadius;
   double tkDxy = tk.dxy();
   double s2D = fabs(tk.dxy()) < psRadius ? sqrt(psRadius * psRadius - tkDxy * tkDxy) : 0;
  
   // the field we get from the caller is already in units of GeV/cm
   double dPhi = -s2D * tk.charge() * bz / tk.pt();
  
   return Direction(tk.eta(), tk.phi() + dPhi);
 }

References reco::TrackBase::charge(), HLT_FULL_cff::dPhi, reco::TrackBase::dxy(), reco::TrackBase::eta(), reco::TrackBase::phi(), reco::TrackBase::pt(), mathSSE::sqrt(), thePropagateTracksToRadius, and theReferenceRadius.

Referenced by deposit().

◆ fillVetos()

void muonisolation::PixelTrackExtractor::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 27 of file PixelTrackExtractor.h.

27 {}

◆ veto()

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

private

Definition at line 47 of file PixelTrackExtractor.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 PixelTrackExtractor::vetos	(	const edm::Event &	ev,
		const edm::EventSetup &	evSetup,
		const reco::Track &	track
	)		const

virtual

Definition at line 40 of file PixelTrackExtractor.cc.

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

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

Member Data Documentation

◆ theBeamlineOption

std::string muonisolation::PixelTrackExtractor::theBeamlineOption

private

Veto cone angle.

Definition at line 51 of file PixelTrackExtractor.h.

Referenced by deposit().

◆ theBeamSpotToken

edm::EDGetTokenT<reco::BeamSpot> muonisolation::PixelTrackExtractor::theBeamSpotToken

private

"NONE", "BeamSpotFromEvent"

Definition at line 52 of file PixelTrackExtractor.h.

Referenced by deposit().

◆ theChi2Ndof_Max

double muonisolation::PixelTrackExtractor::theChi2Ndof_Max

private

trk.numberOfValidHits >= theNHits_Min

Definition at line 54 of file PixelTrackExtractor.h.

Referenced by deposit().

◆ theChi2Prob_Min

double muonisolation::PixelTrackExtractor::theChi2Prob_Min

private

trk.normalizedChi2 < theChi2Ndof_Max

Definition at line 55 of file PixelTrackExtractor.h.

Referenced by deposit().

◆ theDepositLabel

std::string muonisolation::PixelTrackExtractor::theDepositLabel

private

Track Collection Token.

Definition at line 46 of file PixelTrackExtractor.h.

◆ theDiff_r

double muonisolation::PixelTrackExtractor::theDiff_r

private

name for deposit

Definition at line 47 of file PixelTrackExtractor.h.

Referenced by deposit().

◆ theDiff_z

double muonisolation::PixelTrackExtractor::theDiff_z

private

transverse distance to vertex

Definition at line 48 of file PixelTrackExtractor.h.

Referenced by deposit().

◆ theDR_Max

double muonisolation::PixelTrackExtractor::theDR_Max

private

z distance to vertex

Definition at line 49 of file PixelTrackExtractor.h.

Referenced by deposit().

◆ theDR_Veto

double muonisolation::PixelTrackExtractor::theDR_Veto

private

Maximum cone angle for deposits.

Definition at line 50 of file PixelTrackExtractor.h.

Referenced by veto().

◆ theDR_VetoPt

double muonisolation::PixelTrackExtractor::theDR_VetoPt

private

.. it is above this threshold

Definition at line 63 of file PixelTrackExtractor.h.

Referenced by deposit().

◆ theFieldToken

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> muonisolation::PixelTrackExtractor::theFieldToken

private

Definition at line 44 of file PixelTrackExtractor.h.

Referenced by deposit().

◆ theNHits_Min

unsigned int muonisolation::PixelTrackExtractor::theNHits_Min

private

BeamSpot name.

Definition at line 53 of file PixelTrackExtractor.h.

Referenced by deposit().

◆ thePropagateTracksToRadius

bool muonisolation::PixelTrackExtractor::thePropagateTracksToRadius

private

min track pt to include into iso deposit

Definition at line 58 of file PixelTrackExtractor.h.

Referenced by directionAtPresetRadius().

◆ thePt_Min

double muonisolation::PixelTrackExtractor::thePt_Min

private

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

Definition at line 56 of file PixelTrackExtractor.h.

Referenced by deposit().

◆ thePtVeto_Min

double muonisolation::PixelTrackExtractor::thePtVeto_Min

private

will veto leading track if

Definition at line 62 of file PixelTrackExtractor.h.

Referenced by deposit().

◆ theReferenceRadius

double muonisolation::PixelTrackExtractor::theReferenceRadius

private

If set to true will compare track eta-phi at ...

Definition at line 59 of file PixelTrackExtractor.h.

Referenced by directionAtPresetRadius().

◆ theTrackCollectionToken

edm::EDGetTokenT<reco::TrackCollection> muonisolation::PixelTrackExtractor::theTrackCollectionToken

private

Definition at line 45 of file PixelTrackExtractor.h.

Referenced by deposit().

◆ theVetoLeadingTrack

bool muonisolation::PixelTrackExtractor::theVetoLeadingTrack

private

... this radius

Definition at line 61 of file PixelTrackExtractor.h.

Referenced by deposit().

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PixelTrackExtractor() [1/2]

◆ PixelTrackExtractor() [2/2]

◆ ~PixelTrackExtractor()

Member Function Documentation

◆ deposit()

◆ directionAtPresetRadius()

◆ fillVetos()

◆ veto()

◆ vetos()

Member Data Documentation

◆ theBeamlineOption

◆ theBeamSpotToken

◆ theChi2Ndof_Max

◆ theChi2Prob_Min

◆ theDepositLabel

◆ theDiff_r

◆ theDiff_z

◆ theDR_Max

◆ theDR_Veto

◆ theDR_VetoPt

◆ theFieldToken

◆ theNHits_Min

◆ thePropagateTracksToRadius

◆ thePt_Min

◆ thePtVeto_Min

◆ theReferenceRadius

◆ theTrackCollectionToken

◆ theVetoLeadingTrack