CMS 3D CMS Logo

List of all members | Public Member Functions | Private Member Functions | Private Attributes
L1TMuon::GeometryTranslator Class Reference

#include <GeometryTranslator.h>

Public Member Functions

double calculateBendAngle (const TriggerPrimitive &) const
 
double calculateGlobalEta (const TriggerPrimitive &) const
 
double calculateGlobalPhi (const TriggerPrimitive &) const
 
void checkAndUpdateGeometry (const edm::EventSetup &)
 
 GeometryTranslator (edm::ConsumesCollector)
 
const CSCGeometrygetCSCGeometry () const
 
const DTGeometrygetDTGeometry () const
 
const GEMGeometrygetGEMGeometry () const
 
GlobalPoint getGlobalPoint (const TriggerPrimitive &) const
 
const MagneticFieldgetMagneticField () const
 
const ME0GeometrygetME0Geometry () const
 
const RPCGeometrygetRPCGeometry () const
 
 ~GeometryTranslator ()
 

Private Member Functions

double calcCSCSpecificBend (const TriggerPrimitive &) const
 
double calcCSCSpecificEta (const TriggerPrimitive &) const
 
double calcCSCSpecificPhi (const TriggerPrimitive &) const
 
double calcDTSpecificBend (const TriggerPrimitive &) const
 
double calcDTSpecificEta (const TriggerPrimitive &) const
 
double calcDTSpecificPhi (const TriggerPrimitive &) const
 
GlobalPoint calcDTSpecificPoint (const TriggerPrimitive &) const
 
double calcGEMSpecificBend (const TriggerPrimitive &) const
 
double calcGEMSpecificEta (const TriggerPrimitive &) const
 
double calcGEMSpecificPhi (const TriggerPrimitive &) const
 
double calcME0SpecificBend (const TriggerPrimitive &) const
 
double calcME0SpecificEta (const TriggerPrimitive &) const
 
double calcME0SpecificPhi (const TriggerPrimitive &) const
 
double calcRPCSpecificBend (const TriggerPrimitive &) const
 
double calcRPCSpecificEta (const TriggerPrimitive &) const
 
double calcRPCSpecificPhi (const TriggerPrimitive &) const
 
GlobalPoint getCSCSpecificPoint (const TriggerPrimitive &) const
 
GlobalPoint getGEMSpecificPoint (const TriggerPrimitive &) const
 
GlobalPoint getME0SpecificPoint (const TriggerPrimitive &) const
 
GlobalPoint getRPCSpecificPoint (const TriggerPrimitive &) const
 
bool isCSCCounterClockwise (const std::unique_ptr< const CSCLayer > &) const
 

Private Attributes

edm::ESHandle< CSCGeometry_geocsc
 
edm::ESHandle< DTGeometry_geodt
 
edm::ESHandle< GEMGeometry_geogem
 
unsigned long long _geom_cache_id
 
edm::ESHandle< ME0Geometry_geome0
 
edm::ESHandle< RPCGeometry_georpc
 
edm::ESHandle< MagneticField_magfield
 
unsigned long long _magfield_cache_id
 
edm::ESGetToken< CSCGeometry, MuonGeometryRecordgeocscToken_
 
edm::ESGetToken< DTGeometry, MuonGeometryRecordgeodtToken_
 
edm::ESGetToken< GEMGeometry, MuonGeometryRecordgeogemToken_
 
edm::ESGetToken< ME0Geometry, MuonGeometryRecordgeome0Token_
 
edm::ESGetToken< RPCGeometry, MuonGeometryRecordgeorpcToken_
 
edm::ESGetToken< MagneticField, IdealMagneticFieldRecordmagfieldToken_
 

Detailed Description

Definition at line 44 of file GeometryTranslator.h.

Constructor & Destructor Documentation

◆ GeometryTranslator()

GeometryTranslator::GeometryTranslator ( edm::ConsumesCollector  iC)

Definition at line 25 of file GeometryTranslator.cc.

26  : _geom_cache_id(0ULL),
27  geodtToken_(iC.esConsumes()),
32  _magfield_cache_id(0ULL),
33  magfieldToken_(iC.esConsumes()) {}
unsigned long long _geom_cache_id
edm::ESGetToken< DTGeometry, MuonGeometryRecord > geodtToken_
edm::ESGetToken< ME0Geometry, MuonGeometryRecord > geome0Token_
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > magfieldToken_
edm::ESGetToken< RPCGeometry, MuonGeometryRecord > georpcToken_
edm::ESGetToken< CSCGeometry, MuonGeometryRecord > geocscToken_
edm::ESGetToken< GEMGeometry, MuonGeometryRecord > geogemToken_
unsigned long long _magfield_cache_id

◆ ~GeometryTranslator()

GeometryTranslator::~GeometryTranslator ( )

Definition at line 35 of file GeometryTranslator.cc.

35 {}

Member Function Documentation

◆ calcCSCSpecificBend()

double GeometryTranslator::calcCSCSpecificBend ( const TriggerPrimitive tp) const
private

Definition at line 303 of file GeometryTranslator.cc.

References cmsswSequenceInfo::tp.

Referenced by calculateBendAngle().

303 { return tp.getCSCData().bend; }

◆ calcCSCSpecificEta()

double GeometryTranslator::calcCSCSpecificEta ( const TriggerPrimitive tp) const
private

Definition at line 295 of file GeometryTranslator.cc.

References PV3DBase< T, PVType, FrameType >::eta(), getCSCSpecificPoint(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalEta().

295  {
296  return getCSCSpecificPoint(tp).eta();
297 }
GlobalPoint getCSCSpecificPoint(const TriggerPrimitive &) const
T eta() const
Definition: PV3DBase.h:73

◆ calcCSCSpecificPhi()

double GeometryTranslator::calcCSCSpecificPhi ( const TriggerPrimitive tp) const
private

Definition at line 299 of file GeometryTranslator.cc.

References getCSCSpecificPoint(), PV3DBase< T, PVType, FrameType >::phi(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalPhi().

299  {
300  return getCSCSpecificPoint(tp).phi();
301 }
GlobalPoint getCSCSpecificPoint(const TriggerPrimitive &) const
Geom::Phi< T > phi() const
Definition: PV3DBase.h:66

◆ calcDTSpecificBend()

double GeometryTranslator::calcDTSpecificBend ( const TriggerPrimitive tp) const
private

Definition at line 351 of file GeometryTranslator.cc.

References cmsswSequenceInfo::tp.

Referenced by calculateBendAngle().

351  {
352  int bend = tp.getDTData().bendingAngle;
353  double bendf = bend / 512.0;
354  return bendf;
355 }

◆ calcDTSpecificEta()

double GeometryTranslator::calcDTSpecificEta ( const TriggerPrimitive tp) const
private

Definition at line 346 of file GeometryTranslator.cc.

References calcDTSpecificPoint(), PV3DBase< T, PVType, FrameType >::eta(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalEta().

346 { return calcDTSpecificPoint(tp).eta(); }
T eta() const
Definition: PV3DBase.h:73
GlobalPoint calcDTSpecificPoint(const TriggerPrimitive &) const

◆ calcDTSpecificPhi()

double GeometryTranslator::calcDTSpecificPhi ( const TriggerPrimitive tp) const
private

Definition at line 348 of file GeometryTranslator.cc.

References calcDTSpecificPoint(), PV3DBase< T, PVType, FrameType >::phi(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalPhi().

348 { return calcDTSpecificPoint(tp).phi(); }
Geom::Phi< T > phi() const
Definition: PV3DBase.h:66
GlobalPoint calcDTSpecificPoint(const TriggerPrimitive &) const

◆ calcDTSpecificPoint()

GlobalPoint GeometryTranslator::calcDTSpecificPoint ( const TriggerPrimitive tp) const
private

Definition at line 314 of file GeometryTranslator.cc.

References _geodt, DTGeometry::chamber(), M_PI, PV3DBase< T, PVType, FrameType >::mag(), DTChamberId::station(), PV3DBase< T, PVType, FrameType >::theta(), and cmsswSequenceInfo::tp.

Referenced by calcDTSpecificEta(), calcDTSpecificPhi(), and getGlobalPoint().

314  {
315  const DTChamberId baseid(tp.detId<DTChamberId>());
316  // do not use this pointer for anything other than creating a trig geom
317  std::unique_ptr<DTChamber> chamb(const_cast<DTChamber*>(_geodt->chamber(baseid)));
318  std::unique_ptr<DTTrigGeom> trig_geom(new DTTrigGeom(chamb.get(), false));
319  chamb.release(); // release it here so no one gets funny ideas
320  // super layer one is the theta superlayer in a DT chamber
321  // station 4 does not have a theta super layer
322  // the BTI index from the theta trigger is an OR of some BTI outputs
323  // so, we choose the BTI that's in the middle of the group
324  // as the BTI that we get theta from
325  // TODO:::::>>> need to make sure this ordering doesn't flip under wheel sign
326  const int NBTI_theta = ((baseid.station() != 4) ? trig_geom->nCell(2) : trig_geom->nCell(3));
327  const int bti_group = tp.getDTData().theta_bti_group;
328  const unsigned bti_actual = bti_group * NBTI_theta / 7 + NBTI_theta / 14 + 1;
329  DTBtiId thetaBTI;
330  if (baseid.station() != 4 && bti_group != -1) {
331  thetaBTI = DTBtiId(baseid, 2, bti_actual);
332  } else {
333  // since this is phi oriented it'll give us theta in the middle
334  // of the chamber
335  thetaBTI = DTBtiId(baseid, 3, 1);
336  }
337  const GlobalPoint& theta_gp = trig_geom->CMSPosition(thetaBTI);
338 
339  // local phi in sector -> global phi
340  double phi = static_cast<double>(tp.getDTData().radialAngle) / 4096.0; // 12 bits for 1 radian
341  phi += tp.getDTData().sector * M_PI / 6.0; // add sector offset, sector is [0,11]
342 
343  return GlobalPoint(GlobalPoint::Polar(theta_gp.theta(), phi, theta_gp.mag()));
344 }
edm::ESHandle< DTGeometry > _geodt
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
T mag() const
Definition: PV3DBase.h:64
#define M_PI
const DTChamber * chamber(const DTChamberId &id) const
Return a DTChamber given its id.
Definition: DTGeometry.cc:90
Geom::Theta< T > theta() const
Definition: PV3DBase.h:72

◆ calcGEMSpecificBend()

double GeometryTranslator::calcGEMSpecificBend ( const TriggerPrimitive tp) const
private

Definition at line 205 of file GeometryTranslator.cc.

Referenced by calculateBendAngle().

205 { return 0.0; }

◆ calcGEMSpecificEta()

double GeometryTranslator::calcGEMSpecificEta ( const TriggerPrimitive tp) const
private

Definition at line 197 of file GeometryTranslator.cc.

References PV3DBase< T, PVType, FrameType >::eta(), getGEMSpecificPoint(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalEta().

197  {
198  return getGEMSpecificPoint(tp).eta();
199 }
T eta() const
Definition: PV3DBase.h:73
GlobalPoint getGEMSpecificPoint(const TriggerPrimitive &) const

◆ calcGEMSpecificPhi()

double GeometryTranslator::calcGEMSpecificPhi ( const TriggerPrimitive tp) const
private

Definition at line 201 of file GeometryTranslator.cc.

References getGEMSpecificPoint(), PV3DBase< T, PVType, FrameType >::phi(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalPhi().

201  {
202  return getGEMSpecificPoint(tp).phi();
203 }
Geom::Phi< T > phi() const
Definition: PV3DBase.h:66
GlobalPoint getGEMSpecificPoint(const TriggerPrimitive &) const

◆ calcME0SpecificBend()

double GeometryTranslator::calcME0SpecificBend ( const TriggerPrimitive tp) const
private

Definition at line 179 of file GeometryTranslator.cc.

References cmsswSequenceInfo::tp.

Referenced by calculateBendAngle().

179  {
180  return tp.getME0Data().deltaphi * (tp.getME0Data().bend == 0 ? 1 : -1);
181 }

◆ calcME0SpecificEta()

double GeometryTranslator::calcME0SpecificEta ( const TriggerPrimitive tp) const
private

Definition at line 171 of file GeometryTranslator.cc.

References PV3DBase< T, PVType, FrameType >::eta(), getME0SpecificPoint(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalEta().

171  {
172  return getME0SpecificPoint(tp).eta();
173 }
T eta() const
Definition: PV3DBase.h:73
GlobalPoint getME0SpecificPoint(const TriggerPrimitive &) const

◆ calcME0SpecificPhi()

double GeometryTranslator::calcME0SpecificPhi ( const TriggerPrimitive tp) const
private

Definition at line 175 of file GeometryTranslator.cc.

References getME0SpecificPoint(), PV3DBase< T, PVType, FrameType >::phi(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalPhi().

175  {
176  return getME0SpecificPoint(tp).phi();
177 }
Geom::Phi< T > phi() const
Definition: PV3DBase.h:66
GlobalPoint getME0SpecificPoint(const TriggerPrimitive &) const

◆ calcRPCSpecificBend()

double GeometryTranslator::calcRPCSpecificBend ( const TriggerPrimitive tp) const
private

Definition at line 231 of file GeometryTranslator.cc.

Referenced by calculateBendAngle().

231 { return 0.0; }

◆ calcRPCSpecificEta()

double GeometryTranslator::calcRPCSpecificEta ( const TriggerPrimitive tp) const
private

Definition at line 221 of file GeometryTranslator.cc.

References PV3DBase< T, PVType, FrameType >::eta(), getRPCSpecificPoint(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalEta().

221  {
222  return getRPCSpecificPoint(tp).eta();
223 }
T eta() const
Definition: PV3DBase.h:73
GlobalPoint getRPCSpecificPoint(const TriggerPrimitive &) const

◆ calcRPCSpecificPhi()

double GeometryTranslator::calcRPCSpecificPhi ( const TriggerPrimitive tp) const
private

Definition at line 225 of file GeometryTranslator.cc.

References getRPCSpecificPoint(), PV3DBase< T, PVType, FrameType >::phi(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalPhi().

225  {
226  return getRPCSpecificPoint(tp).phi();
227 }
Geom::Phi< T > phi() const
Definition: PV3DBase.h:66
GlobalPoint getRPCSpecificPoint(const TriggerPrimitive &) const

◆ calculateBendAngle()

double GeometryTranslator::calculateBendAngle ( const TriggerPrimitive tp) const

Definition at line 83 of file GeometryTranslator.cc.

References calcCSCSpecificBend(), calcDTSpecificBend(), calcGEMSpecificBend(), calcME0SpecificBend(), calcRPCSpecificBend(), L1TMuon::kCSC, L1TMuon::kDT, L1TMuon::kGEM, L1TMuon::kME0, L1TMuon::kRPC, dqmiodatasetharvest::nan, and cmsswSequenceInfo::tp.

83  {
84  switch (tp.subsystem()) {
85  case L1TMuon::kDT:
86  return calcDTSpecificBend(tp);
87  break;
88  case L1TMuon::kCSC:
89  return calcCSCSpecificBend(tp);
90  break;
91  case L1TMuon::kRPC:
92  return calcRPCSpecificBend(tp);
93  break;
94  case L1TMuon::kGEM:
95  return calcGEMSpecificBend(tp);
96  break;
97  case L1TMuon::kME0:
98  return calcME0SpecificBend(tp);
99  break;
100  default:
101  return std::nan("Invalid TP type!");
102  break;
103  }
104 }
double calcDTSpecificBend(const TriggerPrimitive &) const
double calcME0SpecificBend(const TriggerPrimitive &) const
double calcGEMSpecificBend(const TriggerPrimitive &) const
double calcCSCSpecificBend(const TriggerPrimitive &) const
double calcRPCSpecificBend(const TriggerPrimitive &) const

◆ calculateGlobalEta()

double GeometryTranslator::calculateGlobalEta ( const TriggerPrimitive tp) const

Definition at line 37 of file GeometryTranslator.cc.

References calcCSCSpecificEta(), calcDTSpecificEta(), calcGEMSpecificEta(), calcME0SpecificEta(), calcRPCSpecificEta(), L1TMuon::kCSC, L1TMuon::kDT, L1TMuon::kGEM, L1TMuon::kME0, L1TMuon::kRPC, dqmiodatasetharvest::nan, and cmsswSequenceInfo::tp.

37  {
38  switch (tp.subsystem()) {
39  case L1TMuon::kDT:
40  return calcDTSpecificEta(tp);
41  break;
42  case L1TMuon::kCSC:
43  return calcCSCSpecificEta(tp);
44  break;
45  case L1TMuon::kRPC:
46  return calcRPCSpecificEta(tp);
47  break;
48  case L1TMuon::kGEM:
49  return calcGEMSpecificEta(tp);
50  break;
51  case L1TMuon::kME0:
52  return calcME0SpecificEta(tp);
53  break;
54  default:
55  return std::nan("Invalid TP type!");
56  break;
57  }
58 }
double calcGEMSpecificEta(const TriggerPrimitive &) const
double calcRPCSpecificEta(const TriggerPrimitive &) const
double calcCSCSpecificEta(const TriggerPrimitive &) const
double calcME0SpecificEta(const TriggerPrimitive &) const
double calcDTSpecificEta(const TriggerPrimitive &) const

◆ calculateGlobalPhi()

double GeometryTranslator::calculateGlobalPhi ( const TriggerPrimitive tp) const

Definition at line 60 of file GeometryTranslator.cc.

References calcCSCSpecificPhi(), calcDTSpecificPhi(), calcGEMSpecificPhi(), calcME0SpecificPhi(), calcRPCSpecificPhi(), L1TMuon::kCSC, L1TMuon::kDT, L1TMuon::kGEM, L1TMuon::kME0, L1TMuon::kRPC, dqmiodatasetharvest::nan, and cmsswSequenceInfo::tp.

60  {
61  switch (tp.subsystem()) {
62  case L1TMuon::kDT:
63  return calcDTSpecificPhi(tp);
64  break;
65  case L1TMuon::kCSC:
66  return calcCSCSpecificPhi(tp);
67  break;
68  case L1TMuon::kRPC:
69  return calcRPCSpecificPhi(tp);
70  break;
71  case L1TMuon::kGEM:
72  return calcGEMSpecificPhi(tp);
73  break;
74  case L1TMuon::kME0:
75  return calcME0SpecificPhi(tp);
76  break;
77  default:
78  return std::nan("Invalid TP type!");
79  break;
80  }
81 }
double calcGEMSpecificPhi(const TriggerPrimitive &) const
double calcCSCSpecificPhi(const TriggerPrimitive &) const
double calcRPCSpecificPhi(const TriggerPrimitive &) const
double calcDTSpecificPhi(const TriggerPrimitive &) const
double calcME0SpecificPhi(const TriggerPrimitive &) const

◆ checkAndUpdateGeometry()

void GeometryTranslator::checkAndUpdateGeometry ( const edm::EventSetup es)

Definition at line 131 of file GeometryTranslator.cc.

References _geocsc, _geodt, _geogem, _geom_cache_id, _geome0, _georpc, _magfield, _magfield_cache_id, geocscToken_, geodtToken_, geogemToken_, relativeConstraints::geom, geome0Token_, georpcToken_, edm::EventSetup::get(), volumeBasedMagneticField_160812_cfi::magfield, and magfieldToken_.

Referenced by Phase2L1TGMTStubProducer::produce(), and EMTFSetup::reload().

131  {
133  unsigned long long geomid = geom.cacheIdentifier();
134  if (_geom_cache_id != geomid) {
135  _geodt = geom.getHandle(geodtToken_);
136  _geocsc = geom.getHandle(geocscToken_);
137  _georpc = geom.getHandle(georpcToken_);
138  _geogem = geom.getHandle(geogemToken_);
139  _geome0 = geom.getHandle(geome0Token_);
140  _geom_cache_id = geomid;
141  }
142 
144  unsigned long long magfieldid = magfield.cacheIdentifier();
145  if (_magfield_cache_id != magfieldid) {
146  _magfield = magfield.getHandle(magfieldToken_);
147  _magfield_cache_id = magfieldid;
148  }
149 }
edm::ESHandle< DTGeometry > _geodt
unsigned long long _geom_cache_id
edm::ESHandle< RPCGeometry > _georpc
edm::ESHandle< GEMGeometry > _geogem
edm::ESHandle< ME0Geometry > _geome0
edm::ESHandle< MagneticField > _magfield
edm::ESGetToken< DTGeometry, MuonGeometryRecord > geodtToken_
edm::ESHandle< CSCGeometry > _geocsc
T get() const
Definition: EventSetup.h:82
edm::ESGetToken< ME0Geometry, MuonGeometryRecord > geome0Token_
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > magfieldToken_
edm::ESGetToken< RPCGeometry, MuonGeometryRecord > georpcToken_
edm::ESGetToken< CSCGeometry, MuonGeometryRecord > geocscToken_
edm::ESGetToken< GEMGeometry, MuonGeometryRecord > geogemToken_
unsigned long long _magfield_cache_id

◆ getCSCGeometry()

const CSCGeometry& L1TMuon::GeometryTranslator::getCSCGeometry ( ) const
inline

Definition at line 58 of file GeometryTranslator.h.

References _geocsc.

58 { return *_geocsc; }
edm::ESHandle< CSCGeometry > _geocsc

◆ getCSCSpecificPoint()

GlobalPoint GeometryTranslator::getCSCSpecificPoint ( const TriggerPrimitive tp) const
private

Definition at line 240 of file GeometryTranslator.cc.

References _geocsc, cms::cuda::assert(), CSCGeometry::chamber(), CSCPatternBank::getLegacyPosition(), triggerObjects_cff::id, isCSCCounterClockwise(), CSCConstants::KEY_ALCT_LAYER, phase1PixelTopology::layer, PV3DBase< T, PVType, FrameType >::mag(), hltrates_dqm_sourceclient-live_cfg::offset, topSingleLeptonDQM_PU_cfi::pattern, PV3DBase< T, PVType, FrameType >::phi(), digitizers_cfi::strip, PV3DBase< T, PVType, FrameType >::theta(), cmsswSequenceInfo::tp, and Geom::Phi< T1, Range >::value().

Referenced by calcCSCSpecificEta(), calcCSCSpecificPhi(), and getGlobalPoint().

240  {
241  const CSCDetId id(tp.detId<CSCDetId>());
242  // we should change this to weak_ptrs at some point
243  // requires introducing std::shared_ptrs to geometry
244  std::unique_ptr<const CSCChamber> chamb(_geocsc->chamber(id));
245  assert(chamb != nullptr); // failed to get CSC chamber
246  std::unique_ptr<const CSCLayerGeometry> layer_geom(chamb->layer(CSCConstants::KEY_ALCT_LAYER)->geometry());
247  std::unique_ptr<const CSCLayer> layer(chamb->layer(CSCConstants::KEY_ALCT_LAYER));
248 
249  const uint16_t halfstrip = tp.getCSCData().strip;
250  const uint16_t pattern = tp.getCSCData().pattern;
251  const uint16_t keyWG = tp.getCSCData().keywire;
252  //const unsigned maxStrips = layer_geom->numberOfStrips();
253 
254  // so we can extend this later
255  // assume TMB2007 half-strips only as baseline
256  double offset = 0.0;
257  switch (1) {
258  case 1:
260  }
261  const unsigned halfstrip_offs = static_cast<unsigned>(0.5 + halfstrip + offset);
262  const unsigned strip = halfstrip_offs / 2 + 1; // geom starts from 1
263 
264  // the rough location of the hit at the ALCT key layer
265  // we will refine this using the half strip information
266  const LocalPoint& coarse_lp = layer_geom->stripWireGroupIntersection(strip, keyWG);
267  const GlobalPoint& coarse_gp = layer->surface().toGlobal(coarse_lp);
268 
269  // the strip width/4.0 gives the offset of the half-strip
270  // center with respect to the strip center
271  const double hs_offset = layer_geom->stripPhiPitch() / 4.0;
272 
273  // determine handedness of the chamber
274  const bool ccw = isCSCCounterClockwise(layer);
275  // we need to subtract the offset of even half strips and add the odd ones
276  const double phi_offset = ((halfstrip_offs % 2 ? 1 : -1) * (ccw ? -hs_offset : hs_offset));
277 
278  // the global eta calculation uses the middle of the strip
279  // so no need to increment it
280  const GlobalPoint final_gp(
281  GlobalPoint::Polar(coarse_gp.theta(), (coarse_gp.phi().value() + phi_offset), coarse_gp.mag()));
282 
283  // We need to add in some notion of the 'error' on trigger primitives
284  // like the width of the wire group by the width of the strip
285  // or something similar
286 
287  // release ownership of the pointers
288  chamb.release();
289  layer_geom.release();
290  layer.release();
291 
292  return final_gp;
293 }
const CSCChamber * chamber(CSCDetId id) const
Return the chamber corresponding to given DetId.
Definition: CSCGeometry.cc:100
Geom::Phi< T > phi() const
Definition: PV3DBase.h:66
assert(be >=bs)
constexpr std::array< uint8_t, layerIndexSize > layer
edm::ESHandle< CSCGeometry > _geocsc
T mag() const
Definition: PV3DBase.h:64
static double getLegacyPosition(int pattern)
T1 value() const
Explicit access to value in case implicit conversion not OK.
Definition: Phi.h:75
bool isCSCCounterClockwise(const std::unique_ptr< const CSCLayer > &) const
Geom::Theta< T > theta() const
Definition: PV3DBase.h:72

◆ getDTGeometry()

const DTGeometry& L1TMuon::GeometryTranslator::getDTGeometry ( ) const
inline

Definition at line 57 of file GeometryTranslator.h.

References _geodt.

57 { return *_geodt; }
edm::ESHandle< DTGeometry > _geodt

◆ getGEMGeometry()

const GEMGeometry& L1TMuon::GeometryTranslator::getGEMGeometry ( ) const
inline

Definition at line 60 of file GeometryTranslator.h.

References _geogem.

60 { return *_geogem; }
edm::ESHandle< GEMGeometry > _geogem

◆ getGEMSpecificPoint()

GlobalPoint GeometryTranslator::getGEMSpecificPoint ( const TriggerPrimitive tp) const
private

Definition at line 185 of file GeometryTranslator.cc.

References _geogem, cms::cuda::assert(), GEMEtaPartition::centreOfPad(), GEMGeometry::etaPartition(), runTauDisplay::gp, triggerObjects_cff::id, GeomDet::surface(), Surface::toGlobal(), and cmsswSequenceInfo::tp.

Referenced by calcGEMSpecificEta(), calcGEMSpecificPhi(), and getGlobalPoint().

185  {
186  const GEMDetId id(tp.detId<GEMDetId>());
187  const GEMEtaPartition* roll = _geogem->etaPartition(id);
188  assert(roll != nullptr); // failed to get GEM roll
189  //const uint16_t pad = tp.getGEMData().pad;
190  // Use half-pad precision, + 0.5 at the end to get the center of the pad (pad starts from 0)
191  const float pad = (0.5 * static_cast<float>(tp.getGEMData().pad_low + tp.getGEMData().pad_hi)) + 0.5f;
192  const LocalPoint& lp = roll->centreOfPad(pad);
193  const GlobalPoint& gp = roll->surface().toGlobal(lp);
194  return gp;
195 }
const GEMEtaPartition * etaPartition(GEMDetId id) const
Return a GEMEtaPartition given its id.
Definition: GEMGeometry.cc:77
edm::ESHandle< GEMGeometry > _geogem
assert(be >=bs)
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
Definition: Surface.h:79
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:37
LocalPoint centreOfPad(int pad) const

◆ getGlobalPoint()

GlobalPoint GeometryTranslator::getGlobalPoint ( const TriggerPrimitive tp) const

Definition at line 106 of file GeometryTranslator.cc.

References calcDTSpecificPoint(), getCSCSpecificPoint(), getGEMSpecificPoint(), getME0SpecificPoint(), getRPCSpecificPoint(), L1TMuon::kCSC, L1TMuon::kDT, L1TMuon::kGEM, L1TMuon::kME0, L1TMuon::kRPC, dqmiodatasetharvest::nan, runTheMatrix::ret, and cmsswSequenceInfo::tp.

Referenced by L1TPhase2GMTEndcapStubProcessor::buildCSCOnlyStub(), L1TPhase2GMTEndcapStubProcessor::buildRPCOnlyStub(), PrimitiveConversion::convert_csc(), PrimitiveConversion::convert_dt(), PrimitiveConversion::convert_gem(), PrimitiveConversion::convert_me0(), and PrimitiveConversion::convert_rpc().

106  {
107  switch (tp.subsystem()) {
108  case L1TMuon::kDT:
109  return calcDTSpecificPoint(tp);
110  break;
111  case L1TMuon::kCSC:
112  return getCSCSpecificPoint(tp);
113  break;
114  case L1TMuon::kRPC:
115  return getRPCSpecificPoint(tp);
116  break;
117  case L1TMuon::kGEM:
118  return getGEMSpecificPoint(tp);
119  break;
120  case L1TMuon::kME0:
121  return getME0SpecificPoint(tp);
122  break;
123  default:
125  GlobalPoint::Polar(std::nan("Invalid TP type!"), std::nan("Invalid TP type!"), std::nan("Invalid TP type!")));
126  return ret;
127  break;
128  }
129 }
GlobalPoint getCSCSpecificPoint(const TriggerPrimitive &) const
ret
prodAgent to be discontinued
GlobalPoint getRPCSpecificPoint(const TriggerPrimitive &) const
GlobalPoint getGEMSpecificPoint(const TriggerPrimitive &) const
GlobalPoint getME0SpecificPoint(const TriggerPrimitive &) const
GlobalPoint calcDTSpecificPoint(const TriggerPrimitive &) const

◆ getMagneticField()

const MagneticField& L1TMuon::GeometryTranslator::getMagneticField ( ) const
inline

Definition at line 63 of file GeometryTranslator.h.

References _magfield.

63 { return *_magfield; }
edm::ESHandle< MagneticField > _magfield

◆ getME0Geometry()

const ME0Geometry& L1TMuon::GeometryTranslator::getME0Geometry ( ) const
inline

Definition at line 61 of file GeometryTranslator.h.

References _geome0.

61 { return *_geome0; }
edm::ESHandle< ME0Geometry > _geome0

◆ getME0SpecificPoint()

GlobalPoint GeometryTranslator::getME0SpecificPoint ( const TriggerPrimitive tp) const
private

Definition at line 153 of file GeometryTranslator.cc.

References _geome0, cms::cuda::assert(), ME0EtaPartition::centreOfStrip(), ME0Geometry::chamber(), relativeConstraints::chamber, ME0Layer::etaPartition(), runTauDisplay::gp, triggerObjects_cff::id, GeomDet::surface(), Surface::toGlobal(), and cmsswSequenceInfo::tp.

Referenced by calcME0SpecificEta(), calcME0SpecificPhi(), and getGlobalPoint().

153  {
154  const ME0DetId id(tp.detId<ME0DetId>());
155  const ME0Chamber* chamber = _geome0->chamber(id);
156  const ME0Layer* keylayer = chamber->layer(3); // ME0 key layer is layer 3
157  int partition = tp.getME0Data().partition; // 'partition' is in half-roll unit
158  int iroll = (partition >> 1) + 1;
159  const ME0EtaPartition* roll = keylayer->etaPartition(iroll);
160  assert(roll != nullptr); // failed to get ME0 roll
161  // See L1Trigger/ME0Trigger/src/ME0TriggerPseudoBuilder.cc
162  int phiposition = tp.getME0Data().phiposition; // 'phiposition' is in half-strip unit
163  int istrip = (phiposition >> 1);
164  int phiposition2 = (phiposition & 0x1);
165  float centreOfStrip = istrip + 0.25 + phiposition2 * 0.5;
166  const LocalPoint& lp = roll->centreOfStrip(centreOfStrip);
167  const GlobalPoint& gp = roll->surface().toGlobal(lp);
168  return gp;
169 }
edm::ESHandle< ME0Geometry > _geome0
assert(be >=bs)
const ME0Chamber * chamber(ME0DetId id) const
Return a chamber given its id.
Definition: ME0Geometry.cc:43
LocalPoint centreOfStrip(int strip) const
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
Definition: Surface.h:79
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:37
const ME0EtaPartition * etaPartition(ME0DetId id) const
Return the eta partition corresponding to the given id.
Definition: ME0Layer.cc:27

◆ getRPCGeometry()

const RPCGeometry& L1TMuon::GeometryTranslator::getRPCGeometry ( ) const
inline

Definition at line 59 of file GeometryTranslator.h.

References _georpc.

Referenced by PrimitiveConversion::convert_rpc(), and EMTFSubsystemCollector::extractPrimitives().

59 { return *_georpc; }
edm::ESHandle< RPCGeometry > _georpc

◆ getRPCSpecificPoint()

GlobalPoint GeometryTranslator::getRPCSpecificPoint ( const TriggerPrimitive tp) const
private

Definition at line 209 of file GeometryTranslator.cc.

References _georpc, cms::cuda::assert(), RPCRoll::centreOfStrip(), runTauDisplay::gp, triggerObjects_cff::id, RPCGeometry::roll(), digitizers_cfi::strip, GeomDet::surface(), Surface::toGlobal(), and cmsswSequenceInfo::tp.

Referenced by calcRPCSpecificEta(), calcRPCSpecificPhi(), and getGlobalPoint().

209  {
210  const RPCDetId id(tp.detId<RPCDetId>());
211  const RPCRoll* roll = _georpc->roll(id);
212  assert(roll != nullptr); // failed to get RPC roll
213  //const int strip = static_cast<int>(tp.getRPCData().strip);
214  // Use half-strip precision, - 0.5 at the end to get the center of the strip (strip starts from 1)
215  const float strip = (0.5 * static_cast<float>(tp.getRPCData().strip_low + tp.getRPCData().strip_hi)) - 0.5f;
216  const LocalPoint& lp = roll->centreOfStrip(strip);
217  const GlobalPoint& gp = roll->surface().toGlobal(lp);
218  return gp;
219 }
edm::ESHandle< RPCGeometry > _georpc
const RPCRoll * roll(RPCDetId id) const
Return a roll given its id.
Definition: RPCGeometry.cc:50
assert(be >=bs)
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
Definition: Surface.h:79
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:37
LocalPoint centreOfStrip(int strip) const
Definition: RPCRoll.cc:26

◆ isCSCCounterClockwise()

bool GeometryTranslator::isCSCCounterClockwise ( const std::unique_ptr< const CSCLayer > &  layer) const
private

Definition at line 305 of file GeometryTranslator.cc.

References funct::abs(), phase1PixelTopology::layer, M_PI, and me0TriggerPseudoDigis_cff::nStrips.

Referenced by getCSCSpecificPoint().

305  {
306  const int nStrips = layer->geometry()->numberOfStrips();
307  const double phi1 = layer->centerOfStrip(1).phi();
308  const double phiN = layer->centerOfStrip(nStrips).phi();
309  return ((std::abs(phi1 - phiN) < M_PI && phi1 >= phiN) || (std::abs(phi1 - phiN) >= M_PI && phi1 < phiN));
310 }
nStrips
1.2 is to make the matching window safely the two nearest strips 0.35 is the size of an ME0 chamber i...
constexpr std::array< uint8_t, layerIndexSize > layer
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
#define M_PI

Member Data Documentation

◆ _geocsc

edm::ESHandle<CSCGeometry> L1TMuon::GeometryTranslator::_geocsc
private

◆ _geodt

edm::ESHandle<DTGeometry> L1TMuon::GeometryTranslator::_geodt
private

Definition at line 67 of file GeometryTranslator.h.

Referenced by calcDTSpecificPoint(), checkAndUpdateGeometry(), and getDTGeometry().

◆ _geogem

edm::ESHandle<GEMGeometry> L1TMuon::GeometryTranslator::_geogem
private

◆ _geom_cache_id

unsigned long long L1TMuon::GeometryTranslator::_geom_cache_id
private

Definition at line 66 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().

◆ _geome0

edm::ESHandle<ME0Geometry> L1TMuon::GeometryTranslator::_geome0
private

◆ _georpc

edm::ESHandle<RPCGeometry> L1TMuon::GeometryTranslator::_georpc
private

◆ _magfield

edm::ESHandle<MagneticField> L1TMuon::GeometryTranslator::_magfield
private

Definition at line 80 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry(), and getMagneticField().

◆ _magfield_cache_id

unsigned long long L1TMuon::GeometryTranslator::_magfield_cache_id
private

Definition at line 79 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().

◆ geocscToken_

edm::ESGetToken<CSCGeometry, MuonGeometryRecord> L1TMuon::GeometryTranslator::geocscToken_
private

Definition at line 74 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().

◆ geodtToken_

edm::ESGetToken<DTGeometry, MuonGeometryRecord> L1TMuon::GeometryTranslator::geodtToken_
private

Definition at line 73 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().

◆ geogemToken_

edm::ESGetToken<GEMGeometry, MuonGeometryRecord> L1TMuon::GeometryTranslator::geogemToken_
private

Definition at line 76 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().

◆ geome0Token_

edm::ESGetToken<ME0Geometry, MuonGeometryRecord> L1TMuon::GeometryTranslator::geome0Token_
private

Definition at line 77 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().

◆ georpcToken_

edm::ESGetToken<RPCGeometry, MuonGeometryRecord> L1TMuon::GeometryTranslator::georpcToken_
private

Definition at line 75 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().

◆ magfieldToken_

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> L1TMuon::GeometryTranslator::magfieldToken_
private

Definition at line 81 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().