#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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

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<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	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.

       : 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");
   }

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

◆ ~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.

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

References should_type_MuonCandidateValuesUsed().

Referenced by produce().

◆ 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.

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

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

Referenced by produce().

◆ produce()

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

overrideprivate

Definition at line 80 of file MuonMETValueMapProducer.cc.

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

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

◆ 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.

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