CMS 3D CMS Logo

Functions
L1TkElectronTrackMatchAlgo Namespace Reference

Functions

GlobalPoint calorimeterPosition (double phi, double eta, double e)
 
double deltaEta (const GlobalPoint &epos, const edm::Ptr< L1TTTrackType > &pTrk)
 
double deltaPhi (const GlobalPoint &epos, const edm::Ptr< L1TTTrackType > &pTrk)
 
double deltaPhiClusterET (BXVector< l1t::EGamma >::const_iterator egIter, const edm::Ptr< L1TTTrackType > &pTrk)
 
double deltaPhiClusterET (l1t::EGammaBxCollection::const_iterator egIter, const edm::Ptr< L1TTTrackType > &pTrk)
 
double deltaR (const GlobalPoint &epos, const edm::Ptr< L1TTTrackType > &pTrk)
 
void doMatch (BXVector< l1t::EGamma >::const_iterator egIter, const edm::Ptr< L1TTTrackType > &pTrk, double &dph, double &dr, double &deta)
 
void doMatch (const GlobalPoint &epos, const edm::Ptr< L1TTTrackType > &pTrk, double &dph, double &dr, double &deta)
 
void doMatch (l1t::EGammaBxCollection::const_iterator egIter, const edm::Ptr< L1TTTrackType > &pTrk, double &dph, double &dr, double &deta)
 
void doMatchClusterET (BXVector< l1t::EGamma >::const_iterator egIter, const edm::Ptr< L1TTTrackType > &pTrk, double &dph, double &dr, double &deta)
 
void doMatchClusterET (l1t::EGammaBxCollection::const_iterator egIter, const edm::Ptr< L1TTTrackType > &pTrk, double &dph, double &dr, double &deta)
 

Function Documentation

◆ calorimeterPosition()

GlobalPoint L1TkElectronTrackMatchAlgo::calorimeterPosition ( double  phi,
double  eta,
double  e 
)

Definition at line 99 of file L1TkElectronTrackMatchAlgo.cc.

99  {
100  double x = 0.;
101  double y = 0.;
102  double z = 0.;
103  double depth = 0.89 * (7.7 + log(e));
104  double theta = 2 * atan(exp(-1 * eta));
105  double r = 0;
106  if (fabs(eta) > max_eb_eta) {
107  double ecalZ = max_eb_z * fabs(eta) / eta;
108 
109  r = ecalZ / cos(2 * atan(exp(-1 * eta))) + depth;
110  x = r * cos(phi) * sin(theta);
111  y = r * sin(phi) * sin(theta);
112  z = r * cos(theta);
113  } else {
114  double zface = sqrt(cos(theta) * cos(theta) / (1 - cos(theta) * cos(theta)) * eb_rperp * eb_rperp);
115  r = sqrt(eb_rperp * eb_rperp + zface * zface) + depth;
116  x = r * cos(phi) * sin(theta);
117  y = r * sin(phi) * sin(theta);
118  z = r * cos(theta);
119  }
120  GlobalPoint pos(x, y, z);
121  return pos;
122  }

References funct::cos(), LEDCalibrationChannels::depth, MillePedeFileConverter_cfg::e, eb_rperp, PVValHelper::eta, JetChargeProducer_cfi::exp, dqm-mbProfile::log, max_eb_eta, max_eb_z, alignCSCRings::r, funct::sin(), mathSSE::sqrt(), and theta().

Referenced by deltaPhiClusterET(), doMatch(), and doMatchClusterET().

◆ deltaEta()

double L1TkElectronTrackMatchAlgo::deltaEta ( const GlobalPoint epos,
const edm::Ptr< L1TTTrackType > &  pTrk 
)

Definition at line 84 of file L1TkElectronTrackMatchAlgo.cc.

84  {
85  double corr_eta = 999.0;
86  double er = epos.perp();
87  double ez = epos.z();
88  double z0 = pTrk->POCA().z();
89  double theta = 0.0;
90  if (ez - z0 >= 0)
91  theta = atan(er / fabs(ez - z0));
92  else
93  theta = M_PI - atan(er / fabs(ez - z0));
94  corr_eta = -1.0 * log(tan(theta / 2.0));
95  double deleta = (corr_eta - pTrk->momentum().eta());
96  return deleta;
97  }

References dqm-mbProfile::log, M_PI, PV3DBase< T, PVType, FrameType >::perp(), funct::tan(), theta(), PV3DBase< T, PVType, FrameType >::z(), and HLTMuonOfflineAnalyzer_cfi::z0.

Referenced by deltaR(), doMatch(), and doMatchClusterET().

◆ deltaPhi()

double L1TkElectronTrackMatchAlgo::deltaPhi ( const GlobalPoint epos,
const edm::Ptr< L1TTTrackType > &  pTrk 
)

Definition at line 52 of file L1TkElectronTrackMatchAlgo.cc.

52  {
53  double er = epos.perp();
54  double curv = pTrk->rInv();
55 
56  double dphi_curv = (asin(er * curv / (2.0)));
57  double trk_phi_ecal = reco::deltaPhi(pTrk->momentum().phi(), dphi_curv);
58 
59  double dphi = reco::deltaPhi(trk_phi_ecal, epos.phi());
60  return dphi;
61  }

References reco::deltaPhi(), PV3DBase< T, PVType, FrameType >::perp(), and PV3DBase< T, PVType, FrameType >::phi().

Referenced by deltaR(), and doMatch().

◆ deltaPhiClusterET() [1/2]

double L1TkElectronTrackMatchAlgo::deltaPhiClusterET ( BXVector< l1t::EGamma >::const_iterator  egIter,
const edm::Ptr< L1TTTrackType > &  pTrk 
)

Referenced by doMatchClusterET().

◆ deltaPhiClusterET() [2/2]

double L1TkElectronTrackMatchAlgo::deltaPhiClusterET ( l1t::EGammaBxCollection::const_iterator  egIter,
const edm::Ptr< L1TTTrackType > &  pTrk 
)

Definition at line 63 of file L1TkElectronTrackMatchAlgo.cc.

63  {
64  GlobalPoint epos = L1TkElectronTrackMatchAlgo::calorimeterPosition(egIter->phi(), egIter->eta(), egIter->energy());
65  double er = epos.perp();
66  double et = egIter->et();
67  double pt = pTrk->momentum().perp();
68  double curv = pTrk->rInv();
69 
70  double dphi_curv = (asin(er * curv * pt / (2.0 * et)));
71  double trk_phi_ecal = reco::deltaPhi(pTrk->momentum().phi(), dphi_curv);
72 
73  double dphi = reco::deltaPhi(trk_phi_ecal, epos.phi());
74  return dphi;
75  }

References calorimeterPosition(), reco::deltaPhi(), EgHLTOffHistBins_cfi::et, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), and DiDispStaMuonMonitor_cfi::pt.

◆ deltaR()

double L1TkElectronTrackMatchAlgo::deltaR ( const GlobalPoint epos,
const edm::Ptr< L1TTTrackType > &  pTrk 
)

Definition at line 77 of file L1TkElectronTrackMatchAlgo.cc.

77  {
78  //double dPhi = fabs(reco::deltaPhi(epos.phi(), pTrk->momentum().phi()));
79  double dPhi = L1TkElectronTrackMatchAlgo::deltaPhi(epos, pTrk);
80  double dEta = deltaEta(epos, pTrk);
81  return sqrt(dPhi * dPhi + dEta * dEta);
82  }

References deltaEta(), deltaPhi(), HLT_FULL_cff::dEta, HLT_FULL_cff::dPhi, and mathSSE::sqrt().

Referenced by doMatch(), and doMatchClusterET().

◆ doMatch() [1/3]

void L1TkElectronTrackMatchAlgo::doMatch ( BXVector< l1t::EGamma >::const_iterator  egIter,
const edm::Ptr< L1TTTrackType > &  pTrk,
double &  dph,
double &  dr,
double &  deta 
)

◆ doMatch() [2/3]

void L1TkElectronTrackMatchAlgo::doMatch ( const GlobalPoint epos,
const edm::Ptr< L1TTTrackType > &  pTrk,
double &  dph,
double &  dr,
double &  deta 
)

Definition at line 46 of file L1TkElectronTrackMatchAlgo.cc.

46  {
47  dph = L1TkElectronTrackMatchAlgo::deltaPhi(epos, pTrk);
49  deta = L1TkElectronTrackMatchAlgo::deltaEta(epos, pTrk);
50  }

References deltaEta(), deltaPhi(), deltaR(), and flavorHistoryFilter_cfi::dr.

◆ doMatch() [3/3]

void L1TkElectronTrackMatchAlgo::doMatch ( l1t::EGammaBxCollection::const_iterator  egIter,
const edm::Ptr< L1TTTrackType > &  pTrk,
double &  dph,
double &  dr,
double &  deta 
)

Definition at line 24 of file L1TkElectronTrackMatchAlgo.cc.

28  {
29  GlobalPoint egPos = L1TkElectronTrackMatchAlgo::calorimeterPosition(egIter->phi(), egIter->eta(), egIter->energy());
30  dph = L1TkElectronTrackMatchAlgo::deltaPhi(egPos, pTrk);
32  deta = L1TkElectronTrackMatchAlgo::deltaEta(egPos, pTrk);
33  }

References calorimeterPosition(), deltaEta(), deltaPhi(), deltaR(), and flavorHistoryFilter_cfi::dr.

◆ doMatchClusterET() [1/2]

void L1TkElectronTrackMatchAlgo::doMatchClusterET ( BXVector< l1t::EGamma >::const_iterator  egIter,
const edm::Ptr< L1TTTrackType > &  pTrk,
double &  dph,
double &  dr,
double &  deta 
)

◆ doMatchClusterET() [2/2]

void L1TkElectronTrackMatchAlgo::doMatchClusterET ( l1t::EGammaBxCollection::const_iterator  egIter,
const edm::Ptr< L1TTTrackType > &  pTrk,
double &  dph,
double &  dr,
double &  deta 
)

Definition at line 35 of file L1TkElectronTrackMatchAlgo.cc.

39  {
40  GlobalPoint egPos = L1TkElectronTrackMatchAlgo::calorimeterPosition(egIter->phi(), egIter->eta(), egIter->energy());
43  deta = L1TkElectronTrackMatchAlgo::deltaEta(egPos, pTrk);
44  }

References calorimeterPosition(), deltaEta(), deltaPhiClusterET(), deltaR(), and flavorHistoryFilter_cfi::dr.

L1TkElectronTrackMatchAlgo::deltaPhi
double deltaPhi(const GlobalPoint &epos, const edm::Ptr< L1TTTrackType > &pTrk)
Definition: L1TkElectronTrackMatchAlgo.cc:52
detailsBasic3DVector::z
float float float z
Definition: extBasic3DVector.h:14
DiDispStaMuonMonitor_cfi.pt
pt
Definition: DiDispStaMuonMonitor_cfi.py:39
reco::deltaPhi
constexpr double deltaPhi(double phi1, double phi2)
Definition: deltaPhi.h:26
pos
Definition: PixelAliasList.h:18
L1TkElectronTrackMatchAlgo::deltaR
double deltaR(const GlobalPoint &epos, const edm::Ptr< L1TTTrackType > &pTrk)
Definition: L1TkElectronTrackMatchAlgo.cc:77
HLT_FULL_cff.dPhi
dPhi
Definition: HLT_FULL_cff.py:13702
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
PV3DBase::z
T z() const
Definition: PV3DBase.h:61
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
spr::deltaEta
static const double deltaEta
Definition: CaloConstants.h:8
PVValHelper::eta
Definition: PVValidationHelpers.h:69
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
theta
Geom::Theta< T > theta() const
Definition: Basic3DVectorLD.h:150
HLTMuonOfflineAnalyzer_cfi.z0
z0
Definition: HLTMuonOfflineAnalyzer_cfi.py:98
Point3DBase< float, GlobalTag >
LEDCalibrationChannels.depth
depth
Definition: LEDCalibrationChannels.py:65
L1TkElectronTrackMatchAlgo::deltaEta
double deltaEta(const GlobalPoint &epos, const edm::Ptr< L1TTTrackType > &pTrk)
Definition: L1TkElectronTrackMatchAlgo.cc:84
PVValHelper::phi
Definition: PVValidationHelpers.h:68
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:49
EgHLTOffHistBins_cfi.et
et
Definition: EgHLTOffHistBins_cfi.py:8
alignCSCRings.r
r
Definition: alignCSCRings.py:93
HLT_FULL_cff.dEta
dEta
Definition: HLT_FULL_cff.py:13701
flavorHistoryFilter_cfi.dr
dr
Definition: flavorHistoryFilter_cfi.py:37
genVertex_cff.x
x
Definition: genVertex_cff.py:12
L1TkElectronTrackMatchAlgo::calorimeterPosition
GlobalPoint calorimeterPosition(double phi, double eta, double e)
Definition: L1TkElectronTrackMatchAlgo.cc:99
L1TkElectronTrackMatchAlgo::deltaPhiClusterET
double deltaPhiClusterET(BXVector< l1t::EGamma >::const_iterator egIter, const edm::Ptr< L1TTTrackType > &pTrk)
detailsBasic3DVector::y
float float y
Definition: extBasic3DVector.h:14
dqm-mbProfile.log
log
Definition: dqm-mbProfile.py:17
JetChargeProducer_cfi.exp
exp
Definition: JetChargeProducer_cfi.py:6
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:69
PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:66
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37