#include <MuonMETValueMapProducer.h>

Inheritance diagram for cms::MuonMETValueMapProducer:

Public Member Functions
	MuonMETValueMapProducer (const edm::ParameterSet &)

	~MuonMETValueMapProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Private Member Functions
reco::MuonMETCorrectionData::Type	decide_correction_type (const reco::Muon &muon, const math::XYZPoint &beamSpotPosition)

void	determine_deltax_deltay (double &deltax, double &deltay, const reco::Muon &muon, double bfield, const MagneticField &magneticField, edm::Event &iEvent, const edm::EventSetup &iSetup)

void	produce (edm::Event &, const edm::EventSetup &) override

bool	should_type_MuonCandidateValuesUsed (const reco::Muon &muon, const math::XYZPoint &beamSpotPosition)

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

bool	isAlsoTkMu_

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magFieldToken_

double	maxd0_

double	maxEta_

double	maxNormChi2_

int	minnHits_

int	minnValidStaHits_

double	minPt_

edm::EDGetTokenT< edm::View< reco::Muon > >	muonToken_

double	towerEtThreshold_

TrackDetectorAssociator	trackAssociator_

TrackAssociatorParameters	trackAssociatorParameters_

bool	useHO_

bool	useRecHits_

bool	useTrackAssociatorPositions_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 40 of file MuonMETValueMapProducer.h.

Constructor & Destructor Documentation

◆ MuonMETValueMapProducer()

MuonMETValueMapProducer::MuonMETValueMapProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 55 of file MuonMETValueMapProducer.cc.

References beamSpotToken_, edm::ParameterSet::getParameter(), TrackAssociatorParameters::loadParameters(), magFieldToken_, muonToken_, trackAssociator_, trackAssociatorParameters_, and TrackDetectorAssociator::useDefaultPropagator().

       : minPt_(iConfig.getParameter<double>("minPt")),
         maxEta_(iConfig.getParameter<double>("maxEta")),
         isAlsoTkMu_(iConfig.getParameter<bool>("isAlsoTkMu")),
         maxNormChi2_(iConfig.getParameter<double>("maxNormChi2")),
         maxd0_(iConfig.getParameter<double>("maxd0")),
         minnHits_(iConfig.getParameter<int>("minnHits")),
         minnValidStaHits_(iConfig.getParameter<int>("minnValidStaHits")),
         useTrackAssociatorPositions_(iConfig.getParameter<bool>("useTrackAssociatorPositions")),
         useHO_(iConfig.getParameter<bool>("useHO")),
         towerEtThreshold_(iConfig.getParameter<double>("towerEtThreshold")),
         useRecHits_(iConfig.getParameter<bool>("useRecHits")) {
     muonToken_ = consumes<edm::View<reco::Muon>>(iConfig.getParameter<edm::InputTag>("muonInputTag"));
     beamSpotToken_ = consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("beamSpotInputTag"));
     magFieldToken_ = esConsumes<MagneticField, IdealMagneticFieldRecord>();
 
     edm::ParameterSet trackAssociatorParams = iConfig.getParameter<edm::ParameterSet>("TrackAssociatorParameters");
     edm::ConsumesCollector iC = consumesCollector();
     trackAssociatorParameters_.loadParameters(trackAssociatorParams, iC);
     trackAssociator_.useDefaultPropagator();
 
     produces<edm::ValueMap<reco::MuonMETCorrectionData>>("muCorrData");
   }

◆ ~MuonMETValueMapProducer()

cms::MuonMETValueMapProducer::~MuonMETValueMapProducer ( )

inlineoverride

Definition at line 43 of file MuonMETValueMapProducer.h.

43 {}

Member Function Documentation

◆ decide_correction_type()

reco::MuonMETCorrectionData::Type MuonMETValueMapProducer::decide_correction_type	(	const reco::Muon &	muon,
		const math::XYZPoint &	beamSpotPosition
	)

private

Definition at line 139 of file MuonMETValueMapProducer.cc.

References should_type_MuonCandidateValuesUsed().

Referenced by produce().

                                                                   {
     if (should_type_MuonCandidateValuesUsed(muon, beamSpotPosition))
       return reco::MuonMETCorrectionData::Type::MuonCandidateValuesUsed;
 
     return reco::MuonMETCorrectionData::Type::NotUsed;
   }

◆ determine_deltax_deltay()

void MuonMETValueMapProducer::determine_deltax_deltay	(	double &	deltax,
		double &	deltay,
		const reco::Muon &	muon,
		double	bfield,
		const MagneticField &	magneticField,
		edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

private

Definition at line 115 of file MuonMETValueMapProducer.cc.

References TrackDetectorAssociator::associate(), TrackDetectorAssociator::getFreeTrajectoryState(), MuonMETAlgo::GetMuDepDeltas(), iEvent, info(), HLT_2023v11_cff::magneticField, towerEtThreshold_, trackAssociator_, trackAssociatorParameters_, useHO_, useRecHits_, and useTrackAssociatorPositions_.

Referenced by produce().

                                                                                      {
     reco::TrackRef mu_track;
     if (muon.isGlobalMuon())
       mu_track = muon.globalTrack();
     else if (muon.isTrackerMuon() || muon.isRPCMuon() || muon.isGEMMuon() || muon.isME0Muon())
       mu_track = muon.innerTrack();
     else
       mu_track = muon.outerTrack();
 
     TrackDetMatchInfo info = trackAssociator_.associate(
         iEvent, iSetup, trackAssociator_.getFreeTrajectoryState(&magneticField, *mu_track), trackAssociatorParameters_);
 
     MuonMETAlgo alg;
     alg.GetMuDepDeltas(
         &muon, info, useTrackAssociatorPositions_, useRecHits_, useHO_, towerEtThreshold_, deltax, deltay, bfield);
   }

◆ produce()

void MuonMETValueMapProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 80 of file MuonMETValueMapProducer.cc.

References pwdgSkimBPark_cfi::beamSpot, beamSpotToken_, decide_correction_type(), determine_deltax_deltay(), edm::EventSetup::getData(), iEvent, magFieldToken_, HLT_2023v11_cff::magneticField, eostools::move(), HLT_2023v11_cff::muon, PDWG_BPHSkim_cff::muons, and muonToken_.

                                                                                      {
     edm::Handle<edm::View<reco::Muon>> muons;
     iEvent.getByToken(muonToken_, muons);
 
     edm::Handle<reco::BeamSpot> beamSpot;
     iEvent.getByToken(beamSpotToken_, beamSpot);
 
     const MagneticField& magneticField = iSetup.getData(magFieldToken_);
 
     double bfield = magneticField.inTesla(GlobalPoint(0., 0., 0.)).z();
 
     std::vector<reco::MuonMETCorrectionData> muCorrDataList;
 
     for (edm::View<reco::Muon>::const_iterator muon = muons->begin(); muon != muons->end(); ++muon) {
       double deltax = 0.0;
       double deltay = 0.0;
       determine_deltax_deltay(deltax, deltay, *muon, bfield, magneticField, iEvent, iSetup);
 
       reco::MuonMETCorrectionData::Type muCorrType = decide_correction_type(*muon, beamSpot->position());
 
       reco::MuonMETCorrectionData muMETCorrData(muCorrType, deltax, deltay);
       muCorrDataList.push_back(muMETCorrData);
     }
 
     auto valueMapMuCorrData = std::make_unique<edm::ValueMap<reco::MuonMETCorrectionData>>();
 
     edm::ValueMap<reco::MuonMETCorrectionData>::Filler dataFiller(*valueMapMuCorrData);
 
     dataFiller.insert(muons, muCorrDataList.begin(), muCorrDataList.end());
     dataFiller.fill();
 
     iEvent.put(std::move(valueMapMuCorrData), "muCorrData");
   }

◆ should_type_MuonCandidateValuesUsed()

bool MuonMETValueMapProducer::should_type_MuonCandidateValuesUsed	(	const reco::Muon &	muon,
		const math::XYZPoint &	beamSpotPosition
	)

private

Definition at line 148 of file MuonMETValueMapProducer.cc.

References isAlsoTkMu_, maxd0_, maxEta_, maxNormChi2_, minnHits_, minnValidStaHits_, and minPt_.

Referenced by decide_correction_type().

                                                                                                           {
     if (!muon.isGlobalMuon())
       return false;
     if (!muon.isTrackerMuon() && isAlsoTkMu_)
       return false;
     reco::TrackRef globTk = muon.globalTrack();
     reco::TrackRef siTk = muon.innerTrack();
 
     if (muon.pt() < minPt_ || fabs(muon.eta()) > maxEta_)
       return false;
     if (globTk->chi2() / globTk->ndof() > maxNormChi2_)
       return false;
     if (fabs(globTk->dxy(beamSpotPosition)) > fabs(maxd0_))
       return false;
     if (siTk->numberOfValidHits() < minnHits_)
       return false;
     if (globTk->hitPattern().numberOfValidMuonHits() < minnValidStaHits_)
       return false;
     return true;
   }