MuonGmtPair Class Reference

#include <L1TMuonDQMOffline.h>

Public Member Functions
int	charge () const
double	dR ()
double	eta () const
L1TMuonDQMOffline::EtaRegion	etaRegion () const
double	getDeltaVar (const L1TMuonDQMOffline::ResType) const
double	getVar (const L1TMuonDQMOffline::EffType) const
int	gmtCharge () const
double	gmtEta () const
double	gmtPhi () const
double	gmtPt () const
int	gmtQual () const
	MuonGmtPair (const reco::Muon *muon, const l1t::Muon *regMu, bool useAtVtxCoord)
	MuonGmtPair (const MuonGmtPair &muonGmtPair)
double	phi () const
void	propagate (edm::ESHandle< MagneticField > bField, edm::ESHandle< Propagator > propagatorAlong, edm::ESHandle< Propagator > propagatorOpposite)
double	pt () const
	~MuonGmtPair ()

Private Member Functions
TrajectoryStateOnSurface	cylExtrapTrkSam (reco::TrackRef track, double rho)
FreeTrajectoryState	freeTrajStateMuon (reco::TrackRef track)
TrajectoryStateOnSurface	surfExtrapTrkSam (reco::TrackRef track, double z)

Private Attributes
edm::ESHandle< MagneticField >	m_BField
double	m_eta
double	m_gmtEta
double	m_gmtPhi
const reco::Muon *	m_muon
double	m_phi_bar
double	m_phi_end
edm::ESHandle< Propagator >	m_propagatorAlong
edm::ESHandle< Propagator >	m_propagatorOpposite
const l1t::Muon *	m_regMu

Detailed Description

Definition at line 177 of file L1TMuonDQMOffline.h.

Constructor & Destructor Documentation

MuonGmtPair::MuonGmtPair	(	const reco::Muon *	muon,
		const l1t::Muon *	regMu,
		bool	useAtVtxCoord
	)

Definition at line 25 of file L1TMuonDQMOffline.cc.

References reco::LeafCandidate::eta(), l1t::Muon::etaAtVtx(), m_gmtEta, m_gmtPhi, m_regMu, reco::LeafCandidate::phi(), and l1t::Muon::phiAtVtx().

                                                                                       :
        m_muon(muon), m_regMu(regMu), m_eta(999.), m_phi_bar(999.), m_phi_end(999.)
{
    if (m_regMu) {
        if (useAtVtxCoord) {
            m_gmtEta = m_regMu->etaAtVtx();
            m_gmtPhi = m_regMu->phiAtVtx();
        } else {
            m_gmtEta = m_regMu->eta();
            m_gmtPhi = m_regMu->phi();
        }
    } else {
        m_gmtEta = -5.;
        m_gmtPhi = -5.;
    }
};

MuonGmtPair::MuonGmtPair

(

const MuonGmtPair &

muonGmtPair

)

Definition at line 42 of file L1TMuonDQMOffline.cc.

References m_eta, m_gmtEta, m_gmtPhi, m_muon, m_phi_bar, m_phi_end, and m_regMu.

                                                       {
    m_muon    = muonGmtPair.m_muon;
    m_regMu   = muonGmtPair.m_regMu;

    m_gmtEta  = muonGmtPair.m_gmtEta;
    m_gmtPhi  = muonGmtPair.m_gmtPhi;

    m_eta     = muonGmtPair.m_eta;
    m_phi_bar = muonGmtPair.m_phi_bar;
    m_phi_end = muonGmtPair.m_phi_end;
}

MuonGmtPair::~MuonGmtPair ( )

inline

Definition at line 181 of file L1TMuonDQMOffline.h.

References PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR.

{ };

Member Function Documentation

int MuonGmtPair::charge ( ) const

inline

Definition at line 187 of file L1TMuonDQMOffline.h.

Referenced by getDeltaVar().

{ return m_muon->charge(); };

TrajectoryStateOnSurface MuonGmtPair::cylExtrapTrkSam ( reco::TrackRef  track, double  rho  )

private

Definition at line 109 of file L1TMuonDQMOffline.cc.

References Cylinder::build(), freeTrajStateMuon(), TrajectoryStateOnSurface::isValid(), m_propagatorAlong, m_propagatorOpposite, Propagator::propagate(), and makeMuonMisalignmentScenario::rot.

Referenced by propagate().

{
    Cylinder::PositionType pos(0, 0, 0);
    Cylinder::RotationType rot;
    Cylinder::CylinderPointer myCylinder = Cylinder::build(pos, rot, rho);

    FreeTrajectoryState recoStart = freeTrajStateMuon(track);
    TrajectoryStateOnSurface recoProp;
    recoProp = m_propagatorAlong->propagate(recoStart, *myCylinder);
    if (!recoProp.isValid()) {
      recoProp = m_propagatorOpposite->propagate(recoStart, *myCylinder);
    }
    return recoProp;
}

double MuonGmtPair::dR

(

)

Definition at line 54 of file L1TMuonDQMOffline.cc.

References reco::deltaPhi(), particleFlow_cfi::dEta, particleFlow_cfi::dPhi, eta(), m_gmtEta, m_gmtPhi, m_regMu, phi(), and mathSSE::sqrt().

Referenced by L1TMuonDQMOffline::getMuonGmtPairs().

                       {
    float dEta = m_regMu ? (m_gmtEta - eta()) : 999.;
    float dPhi = m_regMu ? reco::deltaPhi(m_gmtPhi, phi()) : 999.;
    return sqrt(dEta*dEta + dPhi*dPhi);
}

double MuonGmtPair::eta ( ) const

inline

Definition at line 185 of file L1TMuonDQMOffline.h.

Referenced by Particle.Particle::__str__(), dR(), etaRegion(), getDeltaVar(), getVar(), Jet.Jet::jetID(), and Jet.Jet::puJetId().

{ return m_muon->eta(); };

L1TMuonDQMOffline::EtaRegion MuonGmtPair::etaRegion

(

)

const

Definition at line 85 of file L1TMuonDQMOffline.cc.

References funct::abs(), eta(), L1TMuonDQMOffline::kEtaRegionBmtf, L1TMuonDQMOffline::kEtaRegionEmtf, L1TMuonDQMOffline::kEtaRegionOmtf, and L1TMuonDQMOffline::kEtaRegionOut.

                                                        {
    if (std::abs(eta()) < 0.83) return L1TMuonDQMOffline::kEtaRegionBmtf;
    if (std::abs(eta()) < 1.24) return L1TMuonDQMOffline::kEtaRegionOmtf;
    if (std::abs(eta()) < 2.4)  return L1TMuonDQMOffline::kEtaRegionEmtf;
    return L1TMuonDQMOffline::kEtaRegionOut;
}

FreeTrajectoryState MuonGmtPair::freeTrajStateMuon ( reco::TrackRef  track )

private

Definition at line 139 of file L1TMuonDQMOffline.cc.

References m_BField.

Referenced by cylExtrapTrkSam(), and surfExtrapTrkSam().

{
    GlobalPoint  innerPoint(track->innerPosition().x(), track->innerPosition().y(),  track->innerPosition().z());
    GlobalVector innerVec  (track->innerMomentum().x(),  track->innerMomentum().y(),  track->innerMomentum().z());
    FreeTrajectoryState recoStart(innerPoint, innerVec, track->charge(), &*m_BField);
    return recoStart;
}

double MuonGmtPair::getDeltaVar

(

const L1TMuonDQMOffline::ResType

type

)

const

Definition at line 92 of file L1TMuonDQMOffline.cc.

References charge(), reco::deltaPhi(), eta(), gmtCharge(), gmtEta(), gmtPhi(), gmtPt(), L1TMuonDQMOffline::kRes1OverPt, L1TMuonDQMOffline::kResCh, L1TMuonDQMOffline::kResEta, L1TMuonDQMOffline::kResPhi, L1TMuonDQMOffline::kResPt, L1TMuonDQMOffline::kResQOverPt, phi(), and pt().

                                                                         {
    if (type == L1TMuonDQMOffline::kResPt)      return (gmtPt() - pt()) / pt();
    if (type == L1TMuonDQMOffline::kRes1OverPt) return (pt() - gmtPt()) / gmtPt(); // (1/gmtPt - 1/pt) / (1/pt)
    if (type == L1TMuonDQMOffline::kResQOverPt) return (gmtCharge()*charge()*pt() - gmtPt()) / gmtPt(); // (gmtCharge/gmtPt - charge/pt) / (charge/pt) with gmtCharge/charge = gmtCharge*charge
    if (type == L1TMuonDQMOffline::kResPhi)     return reco::deltaPhi(gmtPhi(), phi());
    if (type == L1TMuonDQMOffline::kResEta)     return gmtEta() - eta();
    if (type == L1TMuonDQMOffline::kResCh)      return gmtCharge() - charge();
    return -999.;
}

double MuonGmtPair::getVar

(

const L1TMuonDQMOffline::EffType

type

)

const

Definition at line 102 of file L1TMuonDQMOffline.cc.

References eta(), L1TMuonDQMOffline::kEffEta, L1TMuonDQMOffline::kEffPhi, L1TMuonDQMOffline::kEffPt, phi(), and pt().

                                                                    {
    if (type == L1TMuonDQMOffline::kEffPt)  return pt();
    if (type == L1TMuonDQMOffline::kEffPhi) return phi();
    if (type == L1TMuonDQMOffline::kEffEta) return eta();
    return -999.;
}

int MuonGmtPair::gmtCharge ( ) const

inline

Definition at line 191 of file L1TMuonDQMOffline.h.

Referenced by getDeltaVar().

{return m_regMu ? m_regMu->charge() : -5; };

double MuonGmtPair::gmtEta ( ) const

inline

Definition at line 189 of file L1TMuonDQMOffline.h.

Referenced by getDeltaVar().

{ return m_regMu ? m_gmtEta : -5.; };

double MuonGmtPair::gmtPhi ( ) const

inline

Definition at line 190 of file L1TMuonDQMOffline.h.

Referenced by getDeltaVar().

{ return m_regMu ? m_gmtPhi : -5.; };

double MuonGmtPair::gmtPt ( ) const

inline

Definition at line 188 of file L1TMuonDQMOffline.h.

Referenced by getDeltaVar().

{ return m_regMu ? m_regMu->pt() : -1.; };

int MuonGmtPair::gmtQual ( ) const

inline

Definition at line 192 of file L1TMuonDQMOffline.h.

References ecalTB2006H4_GenSimDigiReco_cfg::bField, gsfElectronCkfTrackCandidateMaker_cff::propagatorAlong, gsfElectronCkfTrackCandidateMaker_cff::propagatorOpposite, rho, HiIsolationCommonParameters_cff::track, and z.

{ return m_regMu ? m_regMu->hwQual() : -1; };

double MuonGmtPair::phi ( void  ) const

inline

Definition at line 186 of file L1TMuonDQMOffline.h.

Referenced by Particle.Particle::__str__(), dR(), getDeltaVar(), getVar(), and ntupleDataFormat.Track::phiPull().

{ return m_muon->phi(); };

void MuonGmtPair::propagate	(	edm::ESHandle< MagneticField >	bField,
		edm::ESHandle< Propagator >	propagatorAlong,
		edm::ESHandle< Propagator >	propagatorOpposite
	)

Definition at line 60 of file L1TMuonDQMOffline.cc.

References ecalTB2006H4_GenSimDigiReco_cfg::bField, cylExtrapTrkSam(), PV3DBase< T, PVType, FrameType >::eta(), TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::isValid(), m_BField, m_eta, m_muon, m_phi_bar, m_phi_end, m_propagatorAlong, m_propagatorOpposite, reco::Muon::outerTrack(), PV3DBase< T, PVType, FrameType >::phi(), gsfElectronCkfTrackCandidateMaker_cff::propagatorAlong, gsfElectronCkfTrackCandidateMaker_cff::propagatorOpposite, and surfExtrapTrkSam().

                                                         {
    m_BField = bField;
    m_propagatorAlong = propagatorAlong;
    m_propagatorOpposite = propagatorOpposite;
    TrackRef standaloneMuon = m_muon->outerTrack();
    TrajectoryStateOnSurface trajectory;
    trajectory = cylExtrapTrkSam(standaloneMuon, 500);  // track at MB2 radius - extrapolation
    if (trajectory.isValid()) {
        m_eta     = trajectory.globalPosition().eta();
        m_phi_bar = trajectory.globalPosition().phi();
    }
    trajectory = surfExtrapTrkSam(standaloneMuon, 790);   // track at ME2+ plane - extrapolation
    if (trajectory.isValid()) {
        m_eta     = trajectory.globalPosition().eta();
        m_phi_end = trajectory.globalPosition().phi();
    }
    trajectory = surfExtrapTrkSam(standaloneMuon, -790); // track at ME2- disk - extrapolation
    if (trajectory.isValid()) {
        m_eta     = trajectory.globalPosition().eta();
        m_phi_end = trajectory.globalPosition().phi();
    }
}

double MuonGmtPair::pt ( ) const

inline

Definition at line 184 of file L1TMuonDQMOffline.h.

Referenced by Particle.Particle::__str__(), DiObject.DiMuon::__str__(), Photon.Photon::calScaledIsoValueExp(), Photon.Photon::calScaledIsoValueLin(), Photon.Photon::calScaledIsoValueQuadr(), getDeltaVar(), getVar(), Electron.Electron::mvaIDLoose(), Electron.Electron::mvaIDRun2(), Electron.Electron::mvaIDTight(), Electron.Electron::ptErr(), ntupleDataFormat.Track::ptPull(), Lepton.Lepton::relIso(), Lepton.Lepton::relIsoFromEA(), Lepton.Lepton::relIsoR(), and Jet.Jet::setCorrP4().

{ return m_muon->pt(); };

TrajectoryStateOnSurface MuonGmtPair::surfExtrapTrkSam ( reco::TrackRef  track, double  z  )

private

Definition at line 124 of file L1TMuonDQMOffline.cc.

References Plane::build(), freeTrajStateMuon(), TrajectoryStateOnSurface::isValid(), m_propagatorAlong, m_propagatorOpposite, Propagator::propagate(), and makeMuonMisalignmentScenario::rot.

Referenced by propagate().

{
    Plane::PositionType pos(0, 0, z);
    Plane::RotationType rot;
    Plane::PlanePointer myPlane = Plane::build(pos, rot);

    FreeTrajectoryState recoStart = freeTrajStateMuon(track);
    TrajectoryStateOnSurface recoProp;
    recoProp = m_propagatorAlong->propagate(recoStart, *myPlane);
    if (!recoProp.isValid()) {
      recoProp = m_propagatorOpposite->propagate(recoStart, *myPlane);
    }
    return recoProp;
}

Member Data Documentation

edm::ESHandle<MagneticField> MuonGmtPair::m_BField

private

Definition at line 212 of file L1TMuonDQMOffline.h.

Referenced by freeTrajStateMuon(), and propagate().

double MuonGmtPair::m_eta

private

Definition at line 221 of file L1TMuonDQMOffline.h.

Referenced by MuonGmtPair(), and propagate().

double MuonGmtPair::m_gmtEta

private

Definition at line 218 of file L1TMuonDQMOffline.h.

Referenced by dR(), and MuonGmtPair().

double MuonGmtPair::m_gmtPhi

private

Definition at line 219 of file L1TMuonDQMOffline.h.

Referenced by dR(), and MuonGmtPair().

const reco::Muon* MuonGmtPair::m_muon

private

Definition at line 209 of file L1TMuonDQMOffline.h.

Referenced by MuonGmtPair(), and propagate().

double MuonGmtPair::m_phi_bar

private

Definition at line 222 of file L1TMuonDQMOffline.h.

Referenced by MuonGmtPair(), and propagate().

double MuonGmtPair::m_phi_end

private

Definition at line 223 of file L1TMuonDQMOffline.h.

Referenced by MuonGmtPair(), and propagate().

edm::ESHandle<Propagator> MuonGmtPair::m_propagatorAlong

private

Definition at line 213 of file L1TMuonDQMOffline.h.

Referenced by cylExtrapTrkSam(), propagate(), and surfExtrapTrkSam().

edm::ESHandle<Propagator> MuonGmtPair::m_propagatorOpposite

private

Definition at line 214 of file L1TMuonDQMOffline.h.

Referenced by cylExtrapTrkSam(), propagate(), and surfExtrapTrkSam().

const l1t::Muon* MuonGmtPair::m_regMu

private

Definition at line 210 of file L1TMuonDQMOffline.h.

Referenced by dR(), and MuonGmtPair().