#include <MuonTrackingRegionByPtBuilder.h>

Inheritance diagram for MuonTrackingRegionByPtBuilder:

Public Member Functions
	MuonTrackingRegionByPtBuilder (const edm::ParameterSet &par, edm::ConsumesCollector &iC)
	Constructor. More...

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

std::unique_ptr< RectangularEtaPhiTrackingRegion >	region (const reco::TrackRef &) const
	Define tracking region. More...

std::unique_ptr< RectangularEtaPhiTrackingRegion >	region (const reco::Track &t) const

std::unique_ptr< RectangularEtaPhiTrackingRegion >	region (const reco::Track &, const edm::Event &, const edm::EventSetup &) const

std::vector< std::unique_ptr< TrackingRegion > >	regions (const edm::Event &, const edm::EventSetup &) const override
	Create Region of Interest. More...

void	setEvent (const edm::Event &, const edm::EventSetup &)
	Pass the Event to the algo at each event. More...

	~MuonTrackingRegionByPtBuilder () override=default
	Destructor. More...

Public Member Functions inherited from TrackingRegionProducer
virtual	~TrackingRegionProducer ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
	Add Fill Descriptions. More...

Private Member Functions
void	build (const edm::ParameterSet &, edm::ConsumesCollector &)

Private Attributes
edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	bfieldToken

std::vector< double >	deltaEtas_

std::vector< double >	deltaPhis_

edm::EDGetTokenT< reco::TrackCollection >	inputCollectionToken

edm::ESGetToken< MultipleScatteringParametrisationMaker, TrackerMultipleScatteringRecord >	msmakerToken

std::vector< double >	ptRanges_

double	theDeltaR

const edm::Event *	theEvent

const edm::EventSetup *	theEventSetup

double	theHalfZ

int	theMaxRegions

edm::EDGetTokenT< MeasurementTrackerEvent >	theMeasurementTrackerToken

double	theNsigmaDz

RectangularEtaPhiTrackingRegion::UseMeasurementTracker	theOnDemand

bool	thePrecise

double	thePtMin

bool	useFixedPt

bool	useFixedZ

bool	useVertex

edm::EDGetTokenT< reco::VertexCollection >	vertexCollectionToken

Detailed Description

Definition at line 26 of file MuonTrackingRegionByPtBuilder.h.

Constructor & Destructor Documentation

◆ MuonTrackingRegionByPtBuilder() [1/2]

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

inlineexplicit

Constructor.

Definition at line 29 of file MuonTrackingRegionByPtBuilder.h.

References build().

29 { build(par, iC); }

MuonTrackingRegionByPtBuilder::build

void build(const edm::ParameterSet &, edm::ConsumesCollector &)

Definition: MuonTrackingRegionByPtBuilder.cc:28

◆ MuonTrackingRegionByPtBuilder() [2/2]

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

inlineexplicit

Definition at line 30 of file MuonTrackingRegionByPtBuilder.h.

References build().

30 { build(par, iC); }

MuonTrackingRegionByPtBuilder::build

void build(const edm::ParameterSet &, edm::ConsumesCollector &)

Definition: MuonTrackingRegionByPtBuilder.cc:28

◆ ~MuonTrackingRegionByPtBuilder()

MuonTrackingRegionByPtBuilder::~MuonTrackingRegionByPtBuilder ( )

overridedefault

Destructor.

Member Function Documentation

◆ build()

void MuonTrackingRegionByPtBuilder::build	(	const edm::ParameterSet &	par,
		edm::ConsumesCollector &	iC
	)

private

Definition at line 28 of file MuonTrackingRegionByPtBuilder.cc.

References beamSpotToken, bfieldToken, edm::ConsumesCollector::consumes(), deltaEtas_, deltaPhis_, edm::ConsumesCollector::esConsumes(), edm::ParameterSet::getParameter(), inputCollectionToken, RectangularEtaPhiTrackingRegion::intToUseMeasurementTracker(), RectangularEtaPhiTrackingRegion::kNever, msmakerToken, ptRanges_, theDeltaR, theHalfZ, theMaxRegions, theMeasurementTrackerToken, theNsigmaDz, theOnDemand, thePrecise, thePtMin, useFixedPt, useFixedZ, useVertex, and vertexCollectionToken.

Referenced by MuonTrackingRegionByPtBuilder().

                                                                                               {
   // Adjust errors on dz
   theNsigmaDz = par.getParameter<double>("Rescale_Dz");
 
   // Flag to switch to use Vertices instead of BeamSpot
   useVertex = par.getParameter<bool>("UseVertex");
 
   // Flag to use fixed limits for Eta, Phi, Z, pT
   useFixedZ = par.getParameter<bool>("Z_fixed");
   useFixedPt = par.getParameter<bool>("Pt_fixed");
 
   // Minimum value for pT
   thePtMin = par.getParameter<double>("Pt_min");
 
   // The static region size along the Z direction
   theHalfZ = par.getParameter<double>("DeltaZ");
 
   // The transverse distance of the region from the BS/PV
   theDeltaR = par.getParameter<double>("DeltaR");
 
   // Parameterized ROIs depending on the pt
   ptRanges_ = par.getParameter<std::vector<double>>("ptRanges");
   if (ptRanges_.size() < 2) {
     edm::LogError("MuonTrackingRegionByPtBuilder") << "Size of ptRanges does not be less than 2." << std::endl;
   }
 
   deltaEtas_ = par.getParameter<std::vector<double>>("deltaEtas");
   if (deltaEtas_.size() != ptRanges_.size() - 1) {
     edm::LogError("MuonTrackingRegionByPtBuilder")
         << "Size of deltaEtas does not match number of pt bins." << std::endl;
   }
 
   deltaPhis_ = par.getParameter<std::vector<double>>("deltaPhis");
   if (deltaPhis_.size() != ptRanges_.size() - 1) {
     edm::LogError("MuonTrackingRegionByPtBuilder")
         << "Size of deltaPhis does not match number of pt bins." << std::endl;
   }
 
   // Maximum number of regions to build when looping over Muons
   theMaxRegions = par.getParameter<int>("maxRegions");
 
   // Flag to use precise??
   thePrecise = par.getParameter<bool>("precise");
 
   // perigee reference point
   theOnDemand = RectangularEtaPhiTrackingRegion::intToUseMeasurementTracker(par.getParameter<int>("OnDemand"));
   if (theOnDemand != RectangularEtaPhiTrackingRegion::UseMeasurementTracker::kNever) {
     theMeasurementTrackerToken =
         iC.consumes<MeasurementTrackerEvent>(par.getParameter<edm::InputTag>("MeasurementTrackerName"));
   }
 
   // Vertex collection and Beam Spot
   beamSpotToken = iC.consumes<reco::BeamSpot>(par.getParameter<edm::InputTag>("beamSpot"));
   vertexCollectionToken = iC.consumes<reco::VertexCollection>(par.getParameter<edm::InputTag>("vertexCollection"));
 
   // Input muon collection
   inputCollectionToken = iC.consumes<reco::TrackCollection>(par.getParameter<edm::InputTag>("input"));
 
   bfieldToken = iC.esConsumes();
   if (thePrecise) {
     msmakerToken = iC.esConsumes();
   }
 }

◆ fillDescriptions()

void MuonTrackingRegionByPtBuilder::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Add Fill Descriptions.

Definition at line 219 of file MuonTrackingRegionByPtBuilder.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, HLT_2022v12_cff::InputTag, and edm::ConfigurationDescriptions::setComment().

                                                                                                {
   edm::ParameterSetDescription desc;
   desc.add<double>("DeltaR", 0.2);
   desc.add<edm::InputTag>("beamSpot", edm::InputTag(""));
   desc.add<int>("OnDemand", -1);
   desc.add<edm::InputTag>("vertexCollection", edm::InputTag(""));
   desc.add<edm::InputTag>("MeasurementTrackerName", edm::InputTag(""));
   desc.add<bool>("UseVertex", false);
   desc.add<double>("Rescale_Dz", 3.0);
   desc.add<bool>("Pt_fixed", false);
   desc.add<bool>("Z_fixed", true);
   desc.add<double>("Pt_min", 1.5);
   desc.add<double>("DeltaZ", 15.9);
   desc.add<std::vector<double>>("ptRanges", {0., 1.e9});
   desc.add<std::vector<double>>("deltaEtas", {0.2});
   desc.add<std::vector<double>>("deltaPhis", {0.15});
   desc.add<int>("maxRegions", 1);
   desc.add<bool>("precise", true);
   desc.add<edm::InputTag>("input", edm::InputTag(""));
   descriptions.add("MuonTrackingRegionByPtBuilder", desc);
 
   descriptions.setComment(
       "Build a TrackingRegion around a standalone muon. Options to define region around beamspot or primary vertex and "
       "dynamic regions are included.");
 }

◆ region() [1/3]

std::unique_ptr< RectangularEtaPhiTrackingRegion > MuonTrackingRegionByPtBuilder::region ( const reco::TrackRef & track ) const

Define tracking region.

Definition at line 114 of file MuonTrackingRegionByPtBuilder.cc.

References HLT_2022v12_cff::track.

Referenced by region(), and regions().

                                      {
   return region(*track);
 }

◆ region() [2/3]

std::unique_ptr<RectangularEtaPhiTrackingRegion> MuonTrackingRegionByPtBuilder::region ( const reco::Track & t ) const

inline

Definition at line 40 of file MuonTrackingRegionByPtBuilder.h.

References region(), submitPVValidationJobs::t, theEvent, and theEventSetup.

                                                                                   {
     return region(t, *theEvent, *theEventSetup);
   }

◆ region() [3/3]

std::unique_ptr< RectangularEtaPhiTrackingRegion > MuonTrackingRegionByPtBuilder::region	(	const reco::Track &	staTrack,
		const edm::Event &	ev,
		const edm::EventSetup &	es
	)		const

Definition at line 130 of file MuonTrackingRegionByPtBuilder.cc.

References funct::abs(), beamSpotToken, bfieldToken, cms::cuda::bs, deltaEtas_, deltaPhis_, PbPb_ZMuSkimMuonDPG_cff::deltaR, L1TkHTMissProducer_cfi::deltaZ, makeMEIFBenchmarkPlots::ev, edm::EventSetup::getData(), reco::Track::innerMomentum(), edm::EDGetTokenT< T >::isUninitialized(), LogDebug, SiStripPI::max, HLTSiStripMonitoring_cff::measurementTracker, beam_dqm_sourceclient-live_cfg::minPt, reco::TrackBase::momentum(), msmakerToken, reco::TrackBase::p(), edm::Handle< T >::product(), DiDispStaMuonMonitor_cfi::pt, reco::TrackBase::pt(), ptRanges_, region(), theDeltaR, theHalfZ, theMeasurementTrackerToken, theNsigmaDz, theOnDemand, thePrecise, thePtMin, pfDeepBoostedJetPreprocessParams_cfi::upper_bound, useFixedPt, useFixedZ, useVertex, spclusmultinvestigator_cfi::vertexCollection, vertexCollectionToken, and extraflags_cff::vtx.

                                                                                   {
   // get track momentum/direction at vertex
   const math::XYZVector& mom = staTrack.momentum();
   GlobalVector dirVector(mom.x(), mom.y(), mom.z());
   double pt = staTrack.pt();
 
   // Fix for StandAlone tracks with low momentum
   const math::XYZVector& innerMomentum = staTrack.innerMomentum();
   GlobalVector forSmallMomentum(innerMomentum.x(), innerMomentum.y(), innerMomentum.z());
   if (staTrack.p() <= 1.5) {
     pt = std::abs(forSmallMomentum.perp());
   }
 
   // initial vertex position - in the following it is replaced with beamspot/vertexing
   GlobalPoint vertexPos(0.0, 0.0, 0.0);
   // standard 15.9, if useVertex than use error from  vertex
   double deltaZ = theHalfZ;
 
   // retrieve beam spot information
   edm::Handle<reco::BeamSpot> bs;
   bool bsHandleFlag = ev.getByToken(beamSpotToken, bs);
 
   // check the validity, otherwise vertexing
   if (bsHandleFlag && bs.isValid() && !useVertex) {
     vertexPos = GlobalPoint(bs->x0(), bs->y0(), bs->z0());
     deltaZ = useFixedZ ? theHalfZ : bs->sigmaZ() * theNsigmaDz;
   } else {
     // get originZPos from list of reconstructed vertices (first or all)
     edm::Handle<reco::VertexCollection> vertexCollection;
     bool vtxHandleFlag = ev.getByToken(vertexCollectionToken, vertexCollection);
     // check if there exists at least one reconstructed vertex
     if (vtxHandleFlag && !vertexCollection->empty()) {
       // use the first vertex in the collection and assume it is the primary event vertex
       reco::VertexCollection::const_iterator vtx = vertexCollection->begin();
       if (!vtx->isFake() && vtx->isValid()) {
         vertexPos = GlobalPoint(vtx->x(), vtx->y(), vtx->z());
         deltaZ = useFixedZ ? theHalfZ : vtx->zError() * theNsigmaDz;
       }
     }
   }
 
   // set delta eta and delta phi depending on the pt
   auto region_dEta = deltaEtas_.at(0);
   auto region_dPhi = deltaPhis_.at(0);
   if (pt < ptRanges_.back()) {
     auto lowEdge = std::upper_bound(ptRanges_.begin(), ptRanges_.end(), pt);
     region_dEta = deltaEtas_.at(lowEdge - ptRanges_.begin() - 1);
     region_dPhi = deltaPhis_.at(lowEdge - ptRanges_.begin() - 1);
   }
 
   float deltaR = theDeltaR;
   double minPt = useFixedPt ? thePtMin : std::max(thePtMin, pt * 0.6);
 
   const MeasurementTrackerEvent* measurementTracker = nullptr;
   if (!theMeasurementTrackerToken.isUninitialized()) {
     edm::Handle<MeasurementTrackerEvent> hmte;
     ev.getByToken(theMeasurementTrackerToken, hmte);
     measurementTracker = hmte.product();
   }
 
   const auto& bfield = es.getData(bfieldToken);
   const MultipleScatteringParametrisationMaker* msmaker = nullptr;
   if (thePrecise) {
     msmaker = &es.getData(msmakerToken);
   }
 
   auto region = std::make_unique<RectangularEtaPhiTrackingRegion>(dirVector,
                                                                   vertexPos,
                                                                   minPt,
                                                                   deltaR,
                                                                   deltaZ,
                                                                   region_dEta,
                                                                   region_dPhi,
                                                                   bfield,
                                                                   msmaker,
                                                                   thePrecise,
                                                                   theOnDemand,
                                                                   measurementTracker);
 
   LogDebug("MuonTrackingRegionByPtBuilder")
       << "the region parameters are:\n"
       << "\n dirVector: " << dirVector << "\n vertexPos: " << vertexPos << "\n minPt: " << minPt
       << "\n deltaR:" << deltaR << "\n deltaZ:" << deltaZ << "\n region_dEta:" << region_dEta
       << "\n region_dPhi:" << region_dPhi << "\n on demand parameter: " << static_cast<int>(theOnDemand);
 
   return region;
 }

◆ regions()

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

overridevirtual

Create Region of Interest.

Implements TrackingRegionProducer.

Definition at line 95 of file MuonTrackingRegionByPtBuilder.cc.

References makeMEIFBenchmarkPlots::ev, inputCollectionToken, region(), mps_fire::result, theMaxRegions, and tracks.

                                                                                                                    {
   std::vector<std::unique_ptr<TrackingRegion>> result;
 
   edm::Handle<reco::TrackCollection> tracks;
   ev.getByToken(inputCollectionToken, tracks);
 
   int nRegions = 0;
   for (auto it = tracks->cbegin(), ed = tracks->cend(); it != ed && nRegions < theMaxRegions; ++it) {
     result.push_back(region(*it, ev, es));
     nRegions++;
   }
 
   return result;
 }

◆ setEvent()

void MuonTrackingRegionByPtBuilder::setEvent	(	const edm::Event &	event,
		const edm::EventSetup &	es
	)

Pass the Event to the algo at each event.

Definition at line 122 of file MuonTrackingRegionByPtBuilder.cc.

References edmPickEvents::event, theEvent, and theEventSetup.

                                                                                            {
   theEvent = &event;
   theEventSetup = &es;
 }

Member Data Documentation

◆ beamSpotToken

edm::EDGetTokenT<reco::BeamSpot> MuonTrackingRegionByPtBuilder::beamSpotToken

private

Definition at line 78 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and region().

◆ bfieldToken

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> MuonTrackingRegionByPtBuilder::bfieldToken

private

Definition at line 81 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and region().

◆ deltaEtas_

std::vector<double> MuonTrackingRegionByPtBuilder::deltaEtas_

private

Definition at line 73 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and region().

◆ deltaPhis_

std::vector<double> MuonTrackingRegionByPtBuilder::deltaPhis_

private

Definition at line 74 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and region().

◆ inputCollectionToken

edm::EDGetTokenT<reco::TrackCollection> MuonTrackingRegionByPtBuilder::inputCollectionToken

private

Definition at line 80 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and regions().

◆ msmakerToken

edm::ESGetToken<MultipleScatteringParametrisationMaker, TrackerMultipleScatteringRecord> MuonTrackingRegionByPtBuilder::msmakerToken

private

Definition at line 82 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and region().

◆ ptRanges_

std::vector<double> MuonTrackingRegionByPtBuilder::ptRanges_

private

Definition at line 72 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and region().

◆ theDeltaR

double MuonTrackingRegionByPtBuilder::theDeltaR

private

Definition at line 69 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and region().

◆ theEvent

const edm::Event* MuonTrackingRegionByPtBuilder::theEvent

private

Definition at line 56 of file MuonTrackingRegionByPtBuilder.h.

Referenced by region(), and setEvent().

◆ theEventSetup

const edm::EventSetup* MuonTrackingRegionByPtBuilder::theEventSetup

private

Definition at line 57 of file MuonTrackingRegionByPtBuilder.h.

Referenced by region(), and setEvent().

◆ theHalfZ

double MuonTrackingRegionByPtBuilder::theHalfZ

private

Definition at line 70 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and region().

◆ theMaxRegions

int MuonTrackingRegionByPtBuilder::theMaxRegions

private

Definition at line 64 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and regions().

◆ theMeasurementTrackerToken

edm::EDGetTokenT<MeasurementTrackerEvent> MuonTrackingRegionByPtBuilder::theMeasurementTrackerToken

private

Definition at line 77 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and region().

◆ theNsigmaDz

double MuonTrackingRegionByPtBuilder::theNsigmaDz

private

Definition at line 66 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and region().

◆ theOnDemand

RectangularEtaPhiTrackingRegion::UseMeasurementTracker MuonTrackingRegionByPtBuilder::theOnDemand

private

Definition at line 76 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and region().

◆ thePrecise

bool MuonTrackingRegionByPtBuilder::thePrecise

private

Definition at line 62 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and region().

◆ thePtMin

double MuonTrackingRegionByPtBuilder::thePtMin

private

Definition at line 68 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and region().

◆ useFixedPt

bool MuonTrackingRegionByPtBuilder::useFixedPt

private

Definition at line 61 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and region().

◆ useFixedZ

bool MuonTrackingRegionByPtBuilder::useFixedZ

private

Definition at line 60 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and region().

◆ useVertex

bool MuonTrackingRegionByPtBuilder::useVertex

private

Definition at line 59 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and region().

◆ vertexCollectionToken

edm::EDGetTokenT<reco::VertexCollection> MuonTrackingRegionByPtBuilder::vertexCollectionToken

private

Definition at line 79 of file MuonTrackingRegionByPtBuilder.h.

Referenced by build(), and region().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ MuonTrackingRegionByPtBuilder() [1/2]

◆ MuonTrackingRegionByPtBuilder() [2/2]

◆ ~MuonTrackingRegionByPtBuilder()

Member Function Documentation

◆ build()

◆ fillDescriptions()

◆ region() [1/3]

◆ region() [2/3]

◆ region() [3/3]

◆ regions()

◆ setEvent()

Member Data Documentation

◆ beamSpotToken

◆ bfieldToken

◆ deltaEtas_

◆ deltaPhis_

◆ inputCollectionToken

◆ msmakerToken

◆ ptRanges_

◆ theDeltaR

◆ theEvent

◆ theEventSetup

◆ theHalfZ

◆ theMaxRegions

◆ theMeasurementTrackerToken

◆ theNsigmaDz

◆ theOnDemand

◆ thePrecise

◆ thePtMin

◆ useFixedPt

◆ useFixedZ

◆ useVertex

◆ vertexCollectionToken