#include <IsolationRegionAroundL3Muon.h>

Inheritance diagram for IsolationRegionAroundL3Muon:

Public Member Functions
	IsolationRegionAroundL3Muon (const edm::ParameterSet &cfg, edm::ConsumesCollector &&iC)

virtual std::vector < TrackingRegion * >	regions (const edm::Event &ev, const edm::EventSetup &es) const

virtual	~IsolationRegionAroundL3Muon ()

Public Member Functions inherited from TrackingRegionProducer
virtual	~TrackingRegionProducer ()

Private Attributes
double	theDeltaEta

double	theDeltaPhi

edm::EDGetTokenT < reco::TrackCollection >	theInputTrkToken

edm::EDGetTokenT < MeasurementTrackerEvent >	theMeasurementTrackerToken

double	theOriginHalfLength

double	theOriginRadius

double	theOriginZPos

double	thePtMin

edm::InputTag	theVertexSrc

edm::EDGetTokenT < reco::VertexCollection >	theVertexToken

bool	theVertexZconstrained

Detailed Description

Definition at line 16 of file IsolationRegionAroundL3Muon.h.

Constructor & Destructor Documentation

IsolationRegionAroundL3Muon::IsolationRegionAroundL3Muon	(	const edm::ParameterSet &	cfg,
		edm::ConsumesCollector &&	iC
	)

inline

Definition at line 20 of file IsolationRegionAroundL3Muon.h.

References edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, theDeltaEta, theDeltaPhi, theInputTrkToken, theMeasurementTrackerToken, theOriginHalfLength, theOriginRadius, theOriginZPos, thePtMin, theVertexSrc, theVertexToken, and theVertexZconstrained.

                                 { 
 
     edm::ParameterSet regionPSet = cfg.getParameter<edm::ParameterSet>("RegionPSet");
 
     theVertexSrc   = regionPSet.getParameter<edm::InputTag>("vertexSrc");
     if (theVertexSrc.label().length()>1) theVertexToken   = iC.consumes<reco::VertexCollection>(theVertexSrc);
     theInputTrkToken = iC.consumes<reco::TrackCollection>(regionPSet.getParameter<edm::InputTag>("TrkSrc"));
 
     thePtMin              = regionPSet.getParameter<double>("ptMin");
     theOriginRadius       = regionPSet.getParameter<double>("originRadius");
     theOriginHalfLength   = regionPSet.getParameter<double>("originHalfLength");
     theVertexZconstrained = regionPSet.getParameter<bool>("vertexZConstrained");
     theOriginZPos  = regionPSet.getParameter<double>("vertexZDefault");
 
     theDeltaEta = regionPSet.getParameter<double>("deltaEtaRegion");
     theDeltaPhi =  regionPSet.getParameter<double>("deltaPhiRegion");
     theMeasurementTrackerToken = iC.consumes<MeasurementTrackerEvent>(regionPSet.getParameter<std::string>("measurementTrackerName"));
   }   

virtual IsolationRegionAroundL3Muon::~IsolationRegionAroundL3Muon ( )

inlinevirtual

Definition at line 40 of file IsolationRegionAroundL3Muon.h.

40 {}

Member Function Documentation

virtual std::vector<TrackingRegion* > IsolationRegionAroundL3Muon::regions	(	const edm::Event &	ev,
		const edm::EventSetup &	es
	)		const

inlinevirtual

Implements TrackingRegionProducer.

Definition at line 42 of file IsolationRegionAroundL3Muon.h.

References edm::Event::getByToken(), RectangularEtaPhiTrackingRegion::kForSiStrips, edm::InputTag::label(), edm::Handle< T >::product(), query::result, theDeltaEta, theDeltaPhi, theInputTrkToken, theMeasurementTrackerToken, theOriginHalfLength, theOriginRadius, theOriginZPos, thePtMin, theVertexSrc, theVertexToken, and theVertexZconstrained.

                                      {
 
     std::vector<TrackingRegion* > result;
 
     // optional constraint for vertex
     // get highest Pt pixel vertex (if existing)
     double deltaZVertex =  theOriginHalfLength;
     double originz = theOriginZPos;
     if (theVertexSrc.label().length()>1) {
       edm::Handle<reco::VertexCollection> vertices;
       ev.getByToken(theVertexToken,vertices);
       const reco::VertexCollection vertCollection = *(vertices.product());
       reco::VertexCollection::const_iterator ci = vertCollection.begin();
       if (vertCollection.size()>0) {
             originz = ci->z();
       } else {
             originz = theOriginZPos;
             deltaZVertex = 15.;
       }
     }
 
     edm::Handle<reco::TrackCollection> trks;
     ev.getByToken(theInputTrkToken, trks);
 
     edm::Handle<MeasurementTrackerEvent> hmte;
     ev.getByToken(theMeasurementTrackerToken, hmte);
     const MeasurementTrackerEvent *measurementTrackerEvent = hmte.product();
 
     for(reco::TrackCollection::const_iterator iTrk = trks->begin();iTrk != trks->end();iTrk++) {
       double vz = (theVertexZconstrained) ? iTrk->dz() : originz;
       GlobalVector dirVector((iTrk)->px(),(iTrk)->py(),(iTrk)->pz());
       result.push_back( 
           new RectangularEtaPhiTrackingRegion( dirVector, GlobalPoint(0,0,float(vz)), 
                            thePtMin, theOriginRadius, deltaZVertex, theDeltaEta, theDeltaPhi,
                            RectangularEtaPhiTrackingRegion::UseMeasurementTracker::kForSiStrips,
                                                true,measurementTrackerEvent) );
     }
 
     return result;
   }