---

muonisolation::PixelTrackExtractor Class Reference

#include <RecoMuon/MuonIsolation/plugins/PixelTrackExtractor.h>

Inheritance diagram for muonisolation::PixelTrackExtractor:

List of all members.

Public Member Functions
virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &muon) const
	make single IsoDeposit based on track as input purely virtual: have to implement in concrete implementations
virtual void	fillVetos (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackCollection &track)
	fill vetoes: to exclude deposits at IsoDeposit creation stage check concrete extractors if it's no-op !
	PixelTrackExtractor (const edm::ParameterSet &par)
	PixelTrackExtractor ()
virtual reco::IsoDeposit::Vetos	vetos (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &track) const
virtual	~PixelTrackExtractor ()
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.
edm::InputTag	theBeamSpotLabel
	"NONE", "BeamSpotFromEvent"
double	theChi2Ndof_Max
	trk.numberOfValidHits >= theNHits_Min
double	theChi2Prob_Min
	trk.normalizedChi2 < theChi2Ndof_Max
std::string	theDepositLabel
	Track Collection Label.
double	theDiff_r
	name for deposit
double	theDiff_z
	transverse distance to vertex
double	theDR_Max
	z distance to vertex
double	theDR_Veto
	Maximum cone angle for deposits.
double	theDR_VetoPt
	.. it is above this threshold
uint	theNHits_Min
	BeamSpot name.
bool	thePropagateTracksToRadius
	min track pt to include into iso deposit
double	thePt_Min
	ChiSquaredProbability(trk.chi2,trk.ndof) > theChi2Prob_Min.
double	thePtVeto_Min
	will veto leading track if
double	theReferenceRadius
	If set to true will compare track eta-phi at ...
edm::InputTag	theTrackCollectionTag
bool	theVetoLeadingTrack
	... this radius

---

Detailed Description

Definition at line 15 of file PixelTrackExtractor.h.

---

Constructor & Destructor Documentation

muonisolation::PixelTrackExtractor::PixelTrackExtractor

(

)

[inline]

Definition at line 19 of file PixelTrackExtractor.h.

{};

PixelTrackExtractor::PixelTrackExtractor

(

const edm::ParameterSet &

par

)

Definition at line 24 of file PixelTrackExtractor.cc.

                                                                  :
  theTrackCollectionTag(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")),
  theBeamSpotLabel(par.getParameter<edm::InputTag>("BeamSpotLabel")),
  theNHits_Min(par.getParameter<uint>("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"))              
{
}

virtual muonisolation::PixelTrackExtractor::~PixelTrackExtractor

(

)

[inline, virtual]

Definition at line 22 of file PixelTrackExtractor.h.

{}

---

Member Function Documentation

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

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

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 76 of file PixelTrackExtractor.cc.

References reco::IsoDeposit::addCandEnergy(), reco::IsoDeposit::addDeposit(), deltaR(), directionAtPresetRadius(), lat::endl(), reco::TrackBase::eta(), edm::EventSetup::get(), edm::Handle< T >::isValid(), LogTrace, pars, reco::TrackBase::phi(), reco::TrackBase::pt(), reco::IsoDeposit::setVeto(), theBeamlineOption, theBeamSpotLabel, theChi2Ndof_Max, theChi2Prob_Min, theDiff_r, theDiff_z, theDR_Max, theDR_VetoPt, theNHits_Min, thePt_Min, thePtVeto_Min, theTrackCollectionTag, theVetoLeadingTrack, veto(), and reco::TrackBase::vz().

{
  static std::string metname = "MuonIsolation|PixelTrackExtractor";

  edm::ESHandle<MagneticField> bField;
  eventSetup.get<IdealMagneticFieldRecord>().get(bField);
  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<View<Track> > tracksH;
  event.getByLabel(theTrackCollectionTag, 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.compare("BeamSpotFromEvent") == 0){
    //pick beamSpot
    reco::BeamSpot beamSpot;
    edm::Handle<reco::BeamSpot> beamSpotH;

    event.getByLabel(theBeamSpotLabel,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;
}

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

Definition at line 60 of file PixelTrackExtractor.cc.

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

Referenced by deposit().

                                                                                       {
  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);
}

virtual void muonisolation::PixelTrackExtractor::fillVetos	(	const edm::Event &	ev,
		const edm::EventSetup &	evSetup,
		const reco::TrackCollection &	tracks
	)			[inline, virtual]

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

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 24 of file PixelTrackExtractor.h.

                                                                        {}

reco::IsoDeposit::Veto muonisolation::PixelTrackExtractor::veto

(

const reco::IsoDeposit::Direction &

dir

)

const [private]

Referenced by deposit(), and vetos().

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

Definition at line 45 of file PixelTrackExtractor.cc.

References dir, reco::TrackBase::eta(), reco::TrackBase::phi(), and veto().

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

---

Member Data Documentation

std::string muonisolation::PixelTrackExtractor::theBeamlineOption [private]

Veto cone angle.

Definition at line 45 of file PixelTrackExtractor.h.

Referenced by deposit().

edm::InputTag muonisolation::PixelTrackExtractor::theBeamSpotLabel [private]

"NONE", "BeamSpotFromEvent"

Definition at line 46 of file PixelTrackExtractor.h.

Referenced by deposit().

double muonisolation::PixelTrackExtractor::theChi2Ndof_Max [private]

trk.numberOfValidHits >= theNHits_Min

Definition at line 48 of file PixelTrackExtractor.h.

Referenced by deposit().

double muonisolation::PixelTrackExtractor::theChi2Prob_Min [private]

trk.normalizedChi2 < theChi2Ndof_Max

Definition at line 49 of file PixelTrackExtractor.h.

Referenced by deposit().

std::string muonisolation::PixelTrackExtractor::theDepositLabel [private]

Track Collection Label.

Definition at line 40 of file PixelTrackExtractor.h.

double muonisolation::PixelTrackExtractor::theDiff_r [private]

name for deposit

Definition at line 41 of file PixelTrackExtractor.h.

Referenced by deposit().

double muonisolation::PixelTrackExtractor::theDiff_z [private]

transverse distance to vertex

Definition at line 42 of file PixelTrackExtractor.h.

Referenced by deposit().

double muonisolation::PixelTrackExtractor::theDR_Max [private]

z distance to vertex

Definition at line 43 of file PixelTrackExtractor.h.

Referenced by deposit().

double muonisolation::PixelTrackExtractor::theDR_Veto [private]

Maximum cone angle for deposits.

Definition at line 44 of file PixelTrackExtractor.h.

double muonisolation::PixelTrackExtractor::theDR_VetoPt [private]

.. it is above this threshold

Definition at line 57 of file PixelTrackExtractor.h.

Referenced by deposit().

uint muonisolation::PixelTrackExtractor::theNHits_Min [private]

BeamSpot name.

Definition at line 47 of file PixelTrackExtractor.h.

Referenced by deposit().

bool muonisolation::PixelTrackExtractor::thePropagateTracksToRadius [private]

min track pt to include into iso deposit

Definition at line 52 of file PixelTrackExtractor.h.

Referenced by directionAtPresetRadius().

double muonisolation::PixelTrackExtractor::thePt_Min [private]

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

Definition at line 50 of file PixelTrackExtractor.h.

Referenced by deposit().

double muonisolation::PixelTrackExtractor::thePtVeto_Min [private]

will veto leading track if

Definition at line 56 of file PixelTrackExtractor.h.

Referenced by deposit().

double muonisolation::PixelTrackExtractor::theReferenceRadius [private]

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

Definition at line 53 of file PixelTrackExtractor.h.

Referenced by directionAtPresetRadius().

edm::InputTag muonisolation::PixelTrackExtractor::theTrackCollectionTag [private]

Definition at line 39 of file PixelTrackExtractor.h.

Referenced by deposit().

bool muonisolation::PixelTrackExtractor::theVetoLeadingTrack [private]

... this radius

Definition at line 55 of file PixelTrackExtractor.h.

Referenced by deposit().

---

The documentation for this class was generated from the following files:

• RecoMuon/MuonIsolation/plugins/PixelTrackExtractor.h
• RecoMuon/MuonIsolation/plugins/PixelTrackExtractor.cc

---