#include <MuonTCMETValueMapProducer.h>

Inheritance diagram for cms::MuonTCMETValueMapProducer:

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

	~MuonTCMETValueMapProducer () 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
bool	isGoodCaloMuon (const reco::Muon *, const unsigned int)

bool	isGoodMuon (const reco::Muon *)

bool	isGoodTrack (const reco::Muon *)

bool	isValidVertex ()

int	nLayers (const reco::TrackRef)

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

class TVector3	propagateTrack (const reco::Muon *)

Private Attributes
edm::Handle< reco::BeamSpot >	beamSpot_

edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_

const class MagneticField *	bField

double	d0cuta_

double	d0cutb_

bool	hasValidVertex

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magFieldToken_

double	maxchi2_

double	maxchi2_tight_

double	maxd0cut_

double	maxeta_

double	maxpt_

double	maxpt_eta20_

double	maxpt_eta25_

double	maxPtErr_

double	maxPtErr_tight_

double	minhits_

double	minhits_tight_

double	minpt_

double	muon_dptrel_

double	muonchi2_

double	muond0_

double	muonDeltaR_

double	muoneta_

bool	muonGlobal_

double	muonhits_

int	muonMinValidStaHits_

double	muonpt_

edm::Handle< reco::MuonCollection >	muons_

edm::EDGetTokenT< reco::MuonCollection >	muonToken_

bool	muonTracker_

int	nLayers_

int	nLayersTight_

class TH2D *	response_function

int	rfType_

TCMETAlgo *	tcmetAlgo_

std::vector< reco::TrackBase::TrackAlgorithm >	trackAlgos_

std::vector< reco::TrackBase::TrackAlgorithm >	trkAlgos_

std::vector< int >	trkQuality_

bool	useCaloMuons_

bool	usePvtxd0_

edm::Handle< reco::VertexCollection >	vertexHandle_

double	vertexMaxDZ_

int	vertexNdof_

double	vertexRho_

edm::EDGetTokenT< reco::VertexCollection >	vertexToken_

double	vertexZ_

const reco::VertexCollection *	vertices_

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 43 of file MuonTCMETValueMapProducer.h.

Constructor & Destructor Documentation

◆ MuonTCMETValueMapProducer()

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

explicit

Definition at line 46 of file MuonTCMETValueMapProducer.cc.

References a, reco::TrackBase::algoByName(), hltPFPuppi_cfi::algos, beamSpotToken_, d0cuta_, d0cutb_, edm::ParameterSet::getParameter(), TCMETAlgo::getResponseFunction_fit(), TCMETAlgo::getResponseFunction_mode(), magFieldToken_, maxchi2_, maxchi2_tight_, maxd0cut_, maxeta_, maxpt_, maxpt_eta20_, maxpt_eta25_, maxPtErr_, maxPtErr_tight_, minhits_, minhits_tight_, minpt_, muon_dptrel_, muonchi2_, muond0_, muonDeltaR_, muoneta_, muonGlobal_, muonhits_, muonMinValidStaHits_, muonpt_, muonToken_, muonTracker_, nLayers_, nLayersTight_, response_function, rfType_, AlCaHLTBitMon_QueryRunRegistry::string, tcmetAlgo_, HcalDetIdTransform::transform(), trkAlgos_, trkQuality_, useCaloMuons_, usePvtxd0_, vertexMaxDZ_, vertexNdof_, vertexRho_, vertexToken_, and vertexZ_.

                                                                                      {
     produces<edm::ValueMap<reco::MuonMETCorrectionData>>("muCorrData");
 
     rfType_ = iConfig.getParameter<int>("rf_type");
 
     nLayers_ = iConfig.getParameter<int>("nLayers");
     nLayersTight_ = iConfig.getParameter<int>("nLayersTight");
     vertexNdof_ = iConfig.getParameter<int>("vertexNdof");
     vertexZ_ = iConfig.getParameter<double>("vertexZ");
     vertexRho_ = iConfig.getParameter<double>("vertexRho");
     vertexMaxDZ_ = iConfig.getParameter<double>("vertexMaxDZ");
     maxpt_eta20_ = iConfig.getParameter<double>("maxpt_eta20");
     maxpt_eta25_ = iConfig.getParameter<double>("maxpt_eta25");
 
     // get configuration parameters
     std::vector<std::string> algos = iConfig.getParameter<std::vector<std::string>>("trackAlgos");
     std::transform(algos.begin(), algos.end(), std::back_inserter(trkAlgos_), [](const std::string& a) {
       return reco::TrackBase::algoByName(a);
     });
     maxd0cut_ = iConfig.getParameter<double>("d0_max");
     minpt_ = iConfig.getParameter<double>("pt_min");
     maxpt_ = iConfig.getParameter<double>("pt_max");
     maxeta_ = iConfig.getParameter<double>("eta_max");
     maxchi2_ = iConfig.getParameter<double>("chi2_max");
     minhits_ = iConfig.getParameter<double>("nhits_min");
     maxPtErr_ = iConfig.getParameter<double>("ptErr_max");
 
     trkQuality_ = iConfig.getParameter<std::vector<int>>("track_quality");
     algos = iConfig.getParameter<std::vector<std::string>>("track_algos");
     std::transform(algos.begin(), algos.end(), std::back_inserter(trkAlgos_), [](const std::string& a) {
       return reco::TrackBase::algoByName(a);
     });
     maxchi2_tight_ = iConfig.getParameter<double>("chi2_max_tight");
     minhits_tight_ = iConfig.getParameter<double>("nhits_min_tight");
     maxPtErr_tight_ = iConfig.getParameter<double>("ptErr_max_tight");
     usePvtxd0_ = iConfig.getParameter<bool>("usePvtxd0");
     d0cuta_ = iConfig.getParameter<double>("d0cuta");
     d0cutb_ = iConfig.getParameter<double>("d0cutb");
 
     muon_dptrel_ = iConfig.getParameter<double>("muon_dptrel");
     muond0_ = iConfig.getParameter<double>("d0_muon");
     muonpt_ = iConfig.getParameter<double>("pt_muon");
     muoneta_ = iConfig.getParameter<double>("eta_muon");
     muonchi2_ = iConfig.getParameter<double>("chi2_muon");
     muonhits_ = iConfig.getParameter<double>("nhits_muon");
     muonGlobal_ = iConfig.getParameter<bool>("global_muon");
     muonTracker_ = iConfig.getParameter<bool>("tracker_muon");
     muonDeltaR_ = iConfig.getParameter<double>("deltaR_muon");
     useCaloMuons_ = iConfig.getParameter<bool>("useCaloMuons");
     muonMinValidStaHits_ = iConfig.getParameter<int>("muonMinValidStaHits");
 
     response_function = nullptr;
     tcmetAlgo_ = new TCMETAlgo();
 
     if (rfType_ == 1)
       response_function = tcmetAlgo_->getResponseFunction_fit();
     else if (rfType_ == 2)
       response_function = tcmetAlgo_->getResponseFunction_mode();
 
     muonToken_ = consumes<reco::MuonCollection>(iConfig.getParameter<edm::InputTag>("muonInputTag"));
     beamSpotToken_ = consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("beamSpotInputTag"));
     vertexToken_ = consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("vertexInputTag"));
     magFieldToken_ = esConsumes<MagneticField, IdealMagneticFieldRecord>();
   }

◆ ~MuonTCMETValueMapProducer()

MuonTCMETValueMapProducer::~MuonTCMETValueMapProducer ( )

override

Definition at line 112 of file MuonTCMETValueMapProducer.cc.

References tcmetAlgo_.

112 { delete tcmetAlgo_; }

cms::MuonTCMETValueMapProducer::tcmetAlgo_

TCMETAlgo * tcmetAlgo_

Definition: MuonTCMETValueMapProducer.h:105

Member Function Documentation

◆ isGoodCaloMuon()

bool MuonTCMETValueMapProducer::isGoodCaloMuon	(	const reco::Muon *	muon,
		const unsigned int	index
	)

private

Definition at line 243 of file MuonTCMETValueMapProducer.cc.

References funct::cos(), spr::deltaEta, SiPixelRawToDigiRegional_cfi::deltaPhi, PbPb_ZMuSkimMuonDPG_cff::deltaR, edm::Ref< C, T, F >::get(), isGoodTrack(), edm::Ref< C, T, F >::isNonnull(), amptDefaultParameters_cff::mu, muonDeltaR_, muons_, and custom_jme_cff::nMuons.

Referenced by produce().

                                                                                                {
     if (muon->pt() < 10)
       return false;
 
     if (!isGoodTrack(muon))
       return false;
 
     const reco::TrackRef inputSiliconTrack = muon->innerTrack();
     if (!inputSiliconTrack.isNonnull())
       return false;
 
     //check if it is in the vicinity of a global or tracker muon
     unsigned int nMuons = muons_->size();
     for (unsigned int iMu = 0; iMu < nMuons; iMu++) {
       if (iMu == index)
         continue;
 
       const reco::Muon* mu = &(*muons_)[iMu];
 
       const reco::TrackRef testSiliconTrack = mu->innerTrack();
       if (!testSiliconTrack.isNonnull())
         continue;
 
       double deltaEta = inputSiliconTrack.get()->eta() - testSiliconTrack.get()->eta();
       double deltaPhi = acos(cos(inputSiliconTrack.get()->phi() - testSiliconTrack.get()->phi()));
       double deltaR = TMath::Sqrt(deltaEta * deltaEta + deltaPhi * deltaPhi);
 
       if (deltaR < muonDeltaR_)
         return false;
     }
 
     return true;
   }

◆ isGoodMuon()

bool MuonTCMETValueMapProducer::isGoodMuon ( const reco::Muon * muon )

private

Definition at line 193 of file MuonTCMETValueMapProducer.cc.

References beamSpot_, hltPixelTracks_cff::chi2, d0, cutBasedMuonId_MuonPOG_V0_cff::globalTrack, edm::Ref< C, T, F >::isNonnull(), edm::HandleBase::isValid(), muon_dptrel_, muonchi2_, muond0_, muoneta_, muonGlobal_, muonhits_, muonMinValidStaHits_, muonpt_, muonTracker_, nhits, and reco::BeamSpot::position().

Referenced by produce().

                                                                  {
     double d0 = -999;
     double nhits = 0;
     double chi2 = 999;
 
     // get d0 corrected for beam spot
     bool haveBeamSpot = true;
     if (!beamSpot_.isValid())
       haveBeamSpot = false;
 
     if (muonGlobal_ && !muon->isGlobalMuon())
       return false;
     if (muonTracker_ && !muon->isTrackerMuon())
       return false;
 
     const reco::TrackRef siTrack = muon->innerTrack();
     const reco::TrackRef globalTrack = muon->globalTrack();
 
     Point bspot = haveBeamSpot ? beamSpot_->position() : Point(0, 0, 0);
     if (siTrack.isNonnull())
       nhits = siTrack->numberOfValidHits();
     if (globalTrack.isNonnull()) {
       d0 = -1 * globalTrack->dxy(bspot);
       chi2 = globalTrack->normalizedChi2();
     }
 
     if (fabs(d0) > muond0_)
       return false;
     if (muon->pt() < muonpt_)
       return false;
     if (fabs(muon->eta()) > muoneta_)
       return false;
     if (nhits < muonhits_)
       return false;
     if (chi2 > muonchi2_)
       return false;
     if (globalTrack->hitPattern().numberOfValidMuonHits() < muonMinValidStaHits_)
       return false;
 
     //reject muons with tracker dpt/pt > X
     if (!siTrack.isNonnull())
       return false;
     if (siTrack->ptError() / siTrack->pt() > muon_dptrel_)
       return false;
 
     else
       return true;
   }

◆ isGoodTrack()

bool MuonTCMETValueMapProducer::isGoodTrack ( const reco::Muon * muon )

private

Definition at line 278 of file MuonTCMETValueMapProducer.cc.

References beamSpot_, DMR_cfg::cut, d0, d0cuta_, d0cutb_, PVValHelper::dz, spr::find(), hasValidVertex, mps_fire::i, edm::Ref< C, T, F >::isNonnull(), edm::HandleBase::isValid(), maxchi2_, maxchi2_tight_, maxd0cut_, maxeta_, maxpt_, maxpt_eta20_, maxpt_eta25_, maxPtErr_, maxPtErr_tight_, minhits_, minhits_tight_, nLayers(), nLayers_, nLayersTight_, reco::BeamSpot::position(), funct::pow(), mathSSE::sqrt(), trackAlgos_, trkAlgos_, trkQuality_, vertexMaxDZ_, and vertices_.

Referenced by isGoodCaloMuon(), and produce().

                                                                   {
     double d0;
 
     const reco::TrackRef siTrack = muon->innerTrack();
     if (!siTrack.isNonnull())
       return false;
 
     if (hasValidVertex) {
       //get d0 corrected for primary vertex
 
       const Point pvtx = Point(vertices_->begin()->x(), vertices_->begin()->y(), vertices_->begin()->z());
 
       double dz = siTrack->dz(pvtx);
 
       if (fabs(dz) < vertexMaxDZ_) {
         //get d0 corrected for pvtx
         d0 = -1 * siTrack->dxy(pvtx);
 
       } else {
         // get d0 corrected for beam spot
         bool haveBeamSpot = true;
         if (!beamSpot_.isValid())
           haveBeamSpot = false;
 
         Point bspot = haveBeamSpot ? beamSpot_->position() : Point(0, 0, 0);
         d0 = -1 * siTrack->dxy(bspot);
       }
     } else {
       // get d0 corrected for beam spot
       bool haveBeamSpot = true;
       if (!beamSpot_.isValid())
         haveBeamSpot = false;
 
       Point bspot = haveBeamSpot ? beamSpot_->position() : Point(0, 0, 0);
       d0 = -1 * siTrack->dxy(bspot);
     }
 
     if (std::find(trackAlgos_.begin(), trackAlgos_.end(), siTrack->algo()) != trackAlgos_.end()) {
       //1st 4 tracking iterations (pT-dependent d0 cut)
 
       float d0cut = sqrt(std::pow(d0cuta_, 2) + std::pow(d0cutb_ / siTrack->pt(), 2));
       if (d0cut > maxd0cut_)
         d0cut = maxd0cut_;
 
       if (fabs(d0) > d0cut)
         return false;
       if (nLayers(siTrack) < nLayers_)
         return false;
     } else {
       //last 2 tracking iterations (tighten chi2, nhits, pt error cuts)
 
       if (siTrack->normalizedChi2() > maxchi2_tight_)
         return false;
       if (siTrack->numberOfValidHits() < minhits_tight_)
         return false;
       if ((siTrack->ptError() / siTrack->pt()) > maxPtErr_tight_)
         return false;
       if (nLayers(siTrack) < nLayersTight_)
         return false;
     }
 
     if (siTrack->numberOfValidHits() < minhits_)
       return false;
     if (siTrack->normalizedChi2() > maxchi2_)
       return false;
     if (fabs(siTrack->eta()) > maxeta_)
       return false;
     if (siTrack->pt() > maxpt_)
       return false;
     if ((siTrack->ptError() / siTrack->pt()) > maxPtErr_)
       return false;
     if (fabs(siTrack->eta()) > 2.5 && siTrack->pt() > maxpt_eta25_)
       return false;
     if (fabs(siTrack->eta()) > 2.0 && siTrack->pt() > maxpt_eta20_)
       return false;
 
     int cut = 0;
     for (unsigned int i = 0; i < trkQuality_.size(); i++) {
       cut |= (1 << trkQuality_.at(i));
     }
 
     if (!((siTrack->qualityMask() & cut) == cut))
       return false;
 
     bool isGoodAlgo = false;
     if (trkAlgos_.empty())
       isGoodAlgo = true;
     for (unsigned int i = 0; i < trkAlgos_.size(); i++) {
       if (siTrack->algo() == trkAlgos_.at(i))
         isGoodAlgo = true;
     }
 
     if (!isGoodAlgo)
       return false;
 
     return true;
   }

◆ isValidVertex()

bool MuonTCMETValueMapProducer::isValidVertex ( )

private

Definition at line 434 of file MuonTCMETValueMapProducer.cc.

References funct::pow(), mathSSE::sqrt(), vertexNdof_, vertexRho_, vertexZ_, and vertices_.

Referenced by produce().

                                                 {
     if (vertices_->begin()->isFake())
       return false;
     if (vertices_->begin()->ndof() < vertexNdof_)
       return false;
     if (fabs(vertices_->begin()->z()) > vertexZ_)
       return false;
     if (sqrt(std::pow(vertices_->begin()->x(), 2) + std::pow(vertices_->begin()->y(), 2)) > vertexRho_)
       return false;
 
     return true;
   }

◆ nLayers()

int MuonTCMETValueMapProducer::nLayers ( const reco::TrackRef track )

private

Definition at line 428 of file MuonTCMETValueMapProducer.cc.

References AlCaHLTBitMon_ParallelJobs::p, and HLT_2023v12_cff::track.

Referenced by isGoodTrack().

                                                                  {
     const reco::HitPattern& p = track->hitPattern();
     return p.trackerLayersWithMeasurement();
   }

◆ produce()

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

overrideprivate

Definition at line 115 of file MuonTCMETValueMapProducer.cc.

References beamSpot_, beamSpotToken_, bField, funct::cos(), edm::helper::Filler< Map >::fill(), edm::EventSetup::getHandle(), hasValidVertex, iEvent, edm::helper::Filler< Map >::insert(), isGoodCaloMuon(), isGoodMuon(), isGoodTrack(), edm::HandleBase::isValid(), isValidVertex(), magFieldToken_, minpt_, eostools::move(), amptDefaultParameters_cff::mu, reco::MuonMETCorrectionData::MuonCandidateValuesUsed, muons_, muonToken_, custom_jme_cff::nMuons, reco::MuonMETCorrectionData::NotUsed, edm::Handle< T >::product(), edm::ESHandle< T >::product(), propagateTrack(), response_function, funct::sin(), reco::MuonMETCorrectionData::TreatedAsPion, useCaloMuons_, usePvtxd0_, vertexHandle_, vertexToken_, and vertices_.

                                                                                        {
     iEvent.getByToken(muonToken_, muons_);
     iEvent.getByToken(beamSpotToken_, beamSpot_);
 
     hasValidVertex = false;
     if (usePvtxd0_) {
       iEvent.getByToken(vertexToken_, vertexHandle_);
 
       if (vertexHandle_.isValid()) {
         vertices_ = vertexHandle_.product();
         hasValidVertex = isValidVertex();
       }
     }
 
     edm::ESHandle<MagneticField> theMagField = iSetup.getHandle(magFieldToken_);
     bField = theMagField.product();
 
     auto vm_muCorrData = std::make_unique<edm::ValueMap<reco::MuonMETCorrectionData>>();
 
     std::vector<reco::MuonMETCorrectionData> v_muCorrData;
 
     unsigned int nMuons = muons_->size();
 
     for (unsigned int iMu = 0; iMu < nMuons; iMu++) {
       const reco::Muon* mu = &(*muons_)[iMu];
       double deltax = 0.0;
       double deltay = 0.0;
 
       reco::MuonMETCorrectionData muMETCorrData(reco::MuonMETCorrectionData::NotUsed, deltax, deltay);
 
       reco::TrackRef mu_track;
       if (mu->isGlobalMuon() || mu->isTrackerMuon() || mu->isCaloMuon())
         mu_track = mu->innerTrack();
       else {
         v_muCorrData.push_back(muMETCorrData);
         continue;
       }
 
       // figure out depositions muons would make if they were treated as pions
       if (isGoodTrack(mu)) {
         if (mu_track->pt() < minpt_)
           muMETCorrData = reco::MuonMETCorrectionData(reco::MuonMETCorrectionData::TreatedAsPion, deltax, deltay);
 
         else {
           int bin_index = response_function->FindBin(mu_track->eta(), mu_track->pt());
           double response = response_function->GetBinContent(bin_index);
 
           TVector3 outerTrkPosition = propagateTrack(mu);
 
           deltax = response * mu_track->p() * sin(outerTrkPosition.Theta()) * cos(outerTrkPosition.Phi());
           deltay = response * mu_track->p() * sin(outerTrkPosition.Theta()) * sin(outerTrkPosition.Phi());
 
           muMETCorrData = reco::MuonMETCorrectionData(reco::MuonMETCorrectionData::TreatedAsPion, deltax, deltay);
         }
       }
 
       // figure out muon flag
       if (isGoodMuon(mu))
         v_muCorrData.push_back(
             reco::MuonMETCorrectionData(reco::MuonMETCorrectionData::MuonCandidateValuesUsed, deltax, deltay));
 
       else if (useCaloMuons_ && isGoodCaloMuon(mu, iMu))
         v_muCorrData.push_back(
             reco::MuonMETCorrectionData(reco::MuonMETCorrectionData::MuonCandidateValuesUsed, deltax, deltay));
 
       else
         v_muCorrData.push_back(muMETCorrData);
     }
 
     edm::ValueMap<reco::MuonMETCorrectionData>::Filler dataFiller(*vm_muCorrData);
 
     dataFiller.insert(muons_, v_muCorrData.begin(), v_muCorrData.end());
     dataFiller.fill();
 
     iEvent.put(std::move(vm_muCorrData), "muCorrData");
   }

◆ propagateTrack()

TVector3 MuonTCMETValueMapProducer::propagateTrack ( const reco::Muon * muon )

private

Definition at line 377 of file MuonTCMETValueMapProducer.cc.

References alongMomentum, Reference_intrackfit_cff::barrel, bField, Plane::build(), Cylinder::build(), TCMET_cfi::corner, TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::isValid(), M_PI, Pi, AnalyticalPropagator::propagate(), CosmicsPD_Skims::radius, HLT_2023v12_cff::track, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::z(), and TCMET_cfi::zdist.

Referenced by produce().

                                                                          {
     TVector3 outerTrkPosition;
 
     outerTrkPosition.SetPtEtaPhi(999., -10., 2 * TMath::Pi());
 
     const reco::TrackRef track = muon->innerTrack();
 
     if (!track.isNonnull()) {
       return outerTrkPosition;
     }
 
     GlobalPoint tpVertex(track->vx(), track->vy(), track->vz());
     GlobalVector tpMomentum(track.get()->px(), track.get()->py(), track.get()->pz());
     int tpCharge(track->charge());
 
     FreeTrajectoryState fts(tpVertex, tpMomentum, tpCharge, bField);
 
     const double zdist = 314.;
 
     const double radius = 130.;
 
     const double corner = 1.479;
 
     Plane::PlanePointer lendcap = Plane::build(Plane::PositionType(0, 0, -zdist), Plane::RotationType());
     Plane::PlanePointer rendcap = Plane::build(Plane::PositionType(0, 0, zdist), Plane::RotationType());
 
     Cylinder::CylinderPointer barrel =
         Cylinder::build(Cylinder::PositionType(0, 0, 0), Cylinder::RotationType(), radius);
 
     AnalyticalPropagator myAP(bField, alongMomentum, 2 * M_PI);
 
     TrajectoryStateOnSurface tsos;
 
     if (track.get()->eta() < -corner) {
       tsos = myAP.propagate(fts, *lendcap);
     } else if (fabs(track.get()->eta()) < corner) {
       tsos = myAP.propagate(fts, *barrel);
     } else if (track.get()->eta() > corner) {
       tsos = myAP.propagate(fts, *rendcap);
     }
 
     if (tsos.isValid())
       outerTrkPosition.SetXYZ(tsos.globalPosition().x(), tsos.globalPosition().y(), tsos.globalPosition().z());
 
     else
       outerTrkPosition.SetPtEtaPhi(999., -10., 2 * TMath::Pi());
 
     return outerTrkPosition;
   }

Member Data Documentation

◆ beamSpot_

edm::Handle<reco::BeamSpot> cms::MuonTCMETValueMapProducer::beamSpot_

private

Definition at line 52 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), isGoodTrack(), and produce().

◆ beamSpotToken_

edm::EDGetTokenT<reco::BeamSpot> cms::MuonTCMETValueMapProducer::beamSpotToken_

private

Definition at line 56 of file MuonTCMETValueMapProducer.h.

Referenced by MuonTCMETValueMapProducer(), and produce().

◆ bField

const class MagneticField* cms::MuonTCMETValueMapProducer::bField

private

Definition at line 60 of file MuonTCMETValueMapProducer.h.

Referenced by produce(), and propagateTrack().

◆ d0cuta_

double cms::MuonTCMETValueMapProducer::d0cuta_

private

Definition at line 91 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

◆ d0cutb_

double cms::MuonTCMETValueMapProducer::d0cutb_

private

Definition at line 92 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

◆ hasValidVertex

bool cms::MuonTCMETValueMapProducer::hasValidVertex

private

Definition at line 69 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and produce().

◆ magFieldToken_

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> cms::MuonTCMETValueMapProducer::magFieldToken_

private

Definition at line 58 of file MuonTCMETValueMapProducer.h.

Referenced by MuonTCMETValueMapProducer(), and produce().

◆ maxchi2_

double cms::MuonTCMETValueMapProducer::maxchi2_

private

Definition at line 84 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

◆ maxchi2_tight_

double cms::MuonTCMETValueMapProducer::maxchi2_tight_

private

Definition at line 88 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

◆ maxd0cut_

double cms::MuonTCMETValueMapProducer::maxd0cut_

private

Definition at line 87 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

◆ maxeta_

double cms::MuonTCMETValueMapProducer::maxeta_

private

Definition at line 83 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

◆ maxpt_

double cms::MuonTCMETValueMapProducer::maxpt_

private

Definition at line 82 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

◆ maxpt_eta20_

double cms::MuonTCMETValueMapProducer::maxpt_eta20_

private

Definition at line 79 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

◆ maxpt_eta25_

double cms::MuonTCMETValueMapProducer::maxpt_eta25_

private

Definition at line 78 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

◆ maxPtErr_

double cms::MuonTCMETValueMapProducer::maxPtErr_

private

Definition at line 86 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

◆ maxPtErr_tight_

double cms::MuonTCMETValueMapProducer::maxPtErr_tight_

private

Definition at line 90 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

◆ minhits_

double cms::MuonTCMETValueMapProducer::minhits_

private

Definition at line 85 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

◆ minhits_tight_

double cms::MuonTCMETValueMapProducer::minhits_tight_

private

Definition at line 89 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

◆ minpt_

double cms::MuonTCMETValueMapProducer::minpt_

private

Definition at line 81 of file MuonTCMETValueMapProducer.h.

Referenced by MuonTCMETValueMapProducer(), and produce().

◆ muon_dptrel_

double cms::MuonTCMETValueMapProducer::muon_dptrel_

private

Definition at line 104 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), and MuonTCMETValueMapProducer().

◆ muonchi2_

double cms::MuonTCMETValueMapProducer::muonchi2_

private

Definition at line 100 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), and MuonTCMETValueMapProducer().

◆ muond0_

double cms::MuonTCMETValueMapProducer::muond0_

private

Definition at line 102 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), and MuonTCMETValueMapProducer().

◆ muonDeltaR_

double cms::MuonTCMETValueMapProducer::muonDeltaR_

private

Definition at line 103 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodCaloMuon(), and MuonTCMETValueMapProducer().

◆ muoneta_

double cms::MuonTCMETValueMapProducer::muoneta_

private

Definition at line 99 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), and MuonTCMETValueMapProducer().

◆ muonGlobal_

bool cms::MuonTCMETValueMapProducer::muonGlobal_

private

Definition at line 66 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), and MuonTCMETValueMapProducer().

◆ muonhits_

double cms::MuonTCMETValueMapProducer::muonhits_

private

Definition at line 101 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), and MuonTCMETValueMapProducer().

◆ muonMinValidStaHits_

int cms::MuonTCMETValueMapProducer::muonMinValidStaHits_

private

Definition at line 97 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), and MuonTCMETValueMapProducer().

◆ muonpt_

double cms::MuonTCMETValueMapProducer::muonpt_

private

Definition at line 98 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), and MuonTCMETValueMapProducer().

◆ muons_

edm::Handle<reco::MuonCollection> cms::MuonTCMETValueMapProducer::muons_

private

Definition at line 51 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodCaloMuon(), and produce().

◆ muonToken_

edm::EDGetTokenT<reco::MuonCollection> cms::MuonTCMETValueMapProducer::muonToken_

private

Definition at line 55 of file MuonTCMETValueMapProducer.h.

Referenced by MuonTCMETValueMapProducer(), and produce().

◆ muonTracker_

bool cms::MuonTCMETValueMapProducer::muonTracker_

private

Definition at line 67 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), and MuonTCMETValueMapProducer().

◆ nLayers_

int cms::MuonTCMETValueMapProducer::nLayers_

private

Definition at line 72 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

◆ nLayersTight_

int cms::MuonTCMETValueMapProducer::nLayersTight_

private

Definition at line 73 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

◆ response_function

class TH2D* cms::MuonTCMETValueMapProducer::response_function

private

Definition at line 64 of file MuonTCMETValueMapProducer.h.

Referenced by MuonTCMETValueMapProducer(), and produce().

◆ rfType_

int cms::MuonTCMETValueMapProducer::rfType_

private

Definition at line 71 of file MuonTCMETValueMapProducer.h.

Referenced by MuonTCMETValueMapProducer().

◆ tcmetAlgo_

TCMETAlgo* cms::MuonTCMETValueMapProducer::tcmetAlgo_

private

Definition at line 105 of file MuonTCMETValueMapProducer.h.

Referenced by MuonTCMETValueMapProducer(), and ~MuonTCMETValueMapProducer().

◆ trackAlgos_

std::vector<reco::TrackBase::TrackAlgorithm> cms::MuonTCMETValueMapProducer::trackAlgos_

private

Definition at line 80 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack().

◆ trkAlgos_

std::vector<reco::TrackBase::TrackAlgorithm> cms::MuonTCMETValueMapProducer::trkAlgos_

private

Definition at line 95 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

◆ trkQuality_

std::vector<int> cms::MuonTCMETValueMapProducer::trkQuality_

private

Definition at line 94 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

◆ useCaloMuons_

bool cms::MuonTCMETValueMapProducer::useCaloMuons_

private

Definition at line 68 of file MuonTCMETValueMapProducer.h.

Referenced by MuonTCMETValueMapProducer(), and produce().

◆ usePvtxd0_

bool cms::MuonTCMETValueMapProducer::usePvtxd0_

private

Definition at line 93 of file MuonTCMETValueMapProducer.h.

Referenced by MuonTCMETValueMapProducer(), and produce().

◆ vertexHandle_

edm::Handle<reco::VertexCollection> cms::MuonTCMETValueMapProducer::vertexHandle_

private

Definition at line 53 of file MuonTCMETValueMapProducer.h.

Referenced by produce().

◆ vertexMaxDZ_

double cms::MuonTCMETValueMapProducer::vertexMaxDZ_

private

Definition at line 77 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

◆ vertexNdof_

int cms::MuonTCMETValueMapProducer::vertexNdof_

private

Definition at line 74 of file MuonTCMETValueMapProducer.h.

Referenced by isValidVertex(), and MuonTCMETValueMapProducer().

◆ vertexRho_

double cms::MuonTCMETValueMapProducer::vertexRho_

private

Definition at line 76 of file MuonTCMETValueMapProducer.h.

Referenced by isValidVertex(), and MuonTCMETValueMapProducer().

◆ vertexToken_

edm::EDGetTokenT<reco::VertexCollection> cms::MuonTCMETValueMapProducer::vertexToken_

private

Definition at line 57 of file MuonTCMETValueMapProducer.h.

Referenced by MuonTCMETValueMapProducer(), and produce().

◆ vertexZ_

double cms::MuonTCMETValueMapProducer::vertexZ_

private

Definition at line 75 of file MuonTCMETValueMapProducer.h.

Referenced by isValidVertex(), and MuonTCMETValueMapProducer().

◆ vertices_

const reco::VertexCollection* cms::MuonTCMETValueMapProducer::vertices_

private

Definition at line 62 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), isValidVertex(), and produce().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MuonTCMETValueMapProducer()

◆ ~MuonTCMETValueMapProducer()

Member Function Documentation

◆ isGoodCaloMuon()

◆ isGoodMuon()

◆ isGoodTrack()

◆ isValidVertex()

◆ nLayers()

◆ produce()

◆ propagateTrack()

Member Data Documentation

◆ beamSpot_

◆ beamSpotToken_

◆ bField

◆ d0cuta_

◆ d0cutb_

◆ hasValidVertex

◆ magFieldToken_

◆ maxchi2_

◆ maxchi2_tight_

◆ maxd0cut_

◆ maxeta_

◆ maxpt_

◆ maxpt_eta20_

◆ maxpt_eta25_

◆ maxPtErr_

◆ maxPtErr_tight_

◆ minhits_

◆ minhits_tight_

◆ minpt_

◆ muon_dptrel_

◆ muonchi2_

◆ muond0_

◆ muonDeltaR_

◆ muoneta_

◆ muonGlobal_

◆ muonhits_

◆ muonMinValidStaHits_

◆ muonpt_

◆ muons_

◆ muonToken_

◆ muonTracker_

◆ nLayers_

◆ nLayersTight_

◆ response_function

◆ rfType_

◆ tcmetAlgo_

◆ trackAlgos_

◆ trkAlgos_

◆ trkQuality_

◆ useCaloMuons_

◆ usePvtxd0_

◆ vertexHandle_

◆ vertexMaxDZ_

◆ vertexNdof_

◆ vertexRho_

◆ vertexToken_

◆ vertexZ_

◆ vertices_