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

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

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

References cmsswSequenceInfo::tp.

Referenced by calculateBendAngle().

◆ calcCSCSpecificEta()

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

Definition at line 295 of file GeometryTranslator.cc.

295  {
296  return getCSCSpecificPoint(tp).eta();
297 }

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

Referenced by calculateGlobalEta().

◆ calcCSCSpecificPhi()

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

Definition at line 299 of file GeometryTranslator.cc.

299  {
300  return getCSCSpecificPoint(tp).phi();
301 }

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

Referenced by calculateGlobalPhi().

◆ calcDTSpecificBend()

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

Definition at line 351 of file GeometryTranslator.cc.

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

References cmsswSequenceInfo::tp.

Referenced by calculateBendAngle().

◆ calcDTSpecificEta()

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

◆ calcDTSpecificPhi()

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

◆ calcDTSpecificPoint()

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

Definition at line 314 of file GeometryTranslator.cc.

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 }

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().

◆ calcGEMSpecificBend()

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

Definition at line 205 of file GeometryTranslator.cc.

205 { return 0.0; }

Referenced by calculateBendAngle().

◆ calcGEMSpecificEta()

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

Definition at line 197 of file GeometryTranslator.cc.

197  {
198  return getGEMSpecificPoint(tp).eta();
199 }

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

Referenced by calculateGlobalEta().

◆ calcGEMSpecificPhi()

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

Definition at line 201 of file GeometryTranslator.cc.

201  {
202  return getGEMSpecificPoint(tp).phi();
203 }

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

Referenced by calculateGlobalPhi().

◆ calcME0SpecificBend()

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

Definition at line 179 of file GeometryTranslator.cc.

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

References cmsswSequenceInfo::tp.

Referenced by calculateBendAngle().

◆ calcME0SpecificEta()

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

Definition at line 171 of file GeometryTranslator.cc.

171  {
172  return getME0SpecificPoint(tp).eta();
173 }

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

Referenced by calculateGlobalEta().

◆ calcME0SpecificPhi()

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

Definition at line 175 of file GeometryTranslator.cc.

175  {
176  return getME0SpecificPoint(tp).phi();
177 }

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

Referenced by calculateGlobalPhi().

◆ calcRPCSpecificBend()

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

Definition at line 231 of file GeometryTranslator.cc.

231 { return 0.0; }

Referenced by calculateBendAngle().

◆ calcRPCSpecificEta()

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

Definition at line 221 of file GeometryTranslator.cc.

221  {
222  return getRPCSpecificPoint(tp).eta();
223 }

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

Referenced by calculateGlobalEta().

◆ calcRPCSpecificPhi()

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

Definition at line 225 of file GeometryTranslator.cc.

225  {
226  return getRPCSpecificPoint(tp).phi();
227 }

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

Referenced by calculateGlobalPhi().

◆ calculateBendAngle()

double GeometryTranslator::calculateBendAngle ( const TriggerPrimitive tp) const

Definition at line 83 of file GeometryTranslator.cc.

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 }

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

◆ calculateGlobalEta()

double GeometryTranslator::calculateGlobalEta ( const TriggerPrimitive tp) const

Definition at line 37 of file GeometryTranslator.cc.

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 }

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

◆ calculateGlobalPhi()

double GeometryTranslator::calculateGlobalPhi ( const TriggerPrimitive tp) const

Definition at line 60 of file GeometryTranslator.cc.

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 }

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

◆ checkAndUpdateGeometry()

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

Definition at line 131 of file GeometryTranslator.cc.

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 }

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 EMTFSetup::reload().

◆ getCSCGeometry()

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

Definition at line 58 of file GeometryTranslator.h.

58 { return *_geocsc; }

References _geocsc.

◆ getCSCSpecificPoint()

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

Definition at line 240 of file GeometryTranslator.cc.

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 }

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().

◆ getDTGeometry()

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

Definition at line 57 of file GeometryTranslator.h.

57 { return *_geodt; }

References _geodt.

◆ getGEMGeometry()

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

Definition at line 60 of file GeometryTranslator.h.

60 { return *_geogem; }

References _geogem.

◆ getGEMSpecificPoint()

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

Definition at line 185 of file GeometryTranslator.cc.

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 }

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().

◆ getGlobalPoint()

GlobalPoint GeometryTranslator::getGlobalPoint ( const TriggerPrimitive tp) const

Definition at line 106 of file GeometryTranslator.cc.

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 }

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 PrimitiveConversion::convert_csc(), PrimitiveConversion::convert_dt(), PrimitiveConversion::convert_gem(), PrimitiveConversion::convert_me0(), and PrimitiveConversion::convert_rpc().

◆ getMagneticField()

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

Definition at line 63 of file GeometryTranslator.h.

63 { return *_magfield; }

References _magfield.

◆ getME0Geometry()

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

Definition at line 61 of file GeometryTranslator.h.

61 { return *_geome0; }

References _geome0.

◆ getME0SpecificPoint()

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

Definition at line 153 of file GeometryTranslator.cc.

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 }

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().

◆ getRPCGeometry()

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

Definition at line 59 of file GeometryTranslator.h.

59 { return *_georpc; }

References _georpc.

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

◆ getRPCSpecificPoint()

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

Definition at line 209 of file GeometryTranslator.cc.

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 }

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().

◆ isCSCCounterClockwise()

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

Definition at line 305 of file GeometryTranslator.cc.

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 }

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

Referenced by getCSCSpecificPoint().

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().

runTheMatrix.ret
ret
prodAgent to be discontinued
Definition: runTheMatrix.py:543
RPCRoll
Definition: RPCRoll.h:12
Geom::Phi::value
T1 value() const
Explicit access to value in case implicit conversion not OK.
Definition: Phi.h:75
L1TMuon::GeometryTranslator::getME0SpecificPoint
GlobalPoint getME0SpecificPoint(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:153
L1TMuon::GeometryTranslator::_magfield_cache_id
unsigned long long _magfield_cache_id
Definition: GeometryTranslator.h:79
L1TMuon::GeometryTranslator::calcGEMSpecificPhi
double calcGEMSpecificPhi(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:201
L1TMuon::GeometryTranslator::calcME0SpecificPhi
double calcME0SpecificPhi(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:175
L1TMuon::GeometryTranslator::getCSCSpecificPoint
GlobalPoint getCSCSpecificPoint(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:240
digitizers_cfi.strip
strip
Definition: digitizers_cfi.py:19
L1TMuon::GeometryTranslator::calcCSCSpecificPhi
double calcCSCSpecificPhi(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:299
PV3DBase::theta
Geom::Theta< T > theta() const
Definition: PV3DBase.h:72
RPCDetId
Definition: RPCDetId.h:16
GEMEtaPartition
Definition: GEMEtaPartition.h:12
dqmiodatasetharvest.nan
nan
Definition: dqmiodatasetharvest.py:39
cms::cuda::assert
assert(be >=bs)
Geom::Spherical2Cartesian
Definition: CoordinateSets.h:58
CSCConstants::KEY_ALCT_LAYER
Definition: CSCConstants.h:103
RPCGeometry::roll
const RPCRoll * roll(RPCDetId id) const
Return a roll given its id.
Definition: RPCGeometry.cc:50
L1TMuon::GeometryTranslator::geodtToken_
edm::ESGetToken< DTGeometry, MuonGeometryRecord > geodtToken_
Definition: GeometryTranslator.h:73
L1TMuon::GeometryTranslator::_georpc
edm::ESHandle< RPCGeometry > _georpc
Definition: GeometryTranslator.h:69
L1TMuon::GeometryTranslator::_geome0
edm::ESHandle< ME0Geometry > _geome0
Definition: GeometryTranslator.h:71
L1TMuon::GeometryTranslator::calcRPCSpecificEta
double calcRPCSpecificEta(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:221
edm::ConsumesCollector::esConsumes
auto esConsumes()
Definition: ConsumesCollector.h:97
DTGeometry::chamber
const DTChamber * chamber(const DTChamberId &id) const
Return a DTChamber given its id.
Definition: DTGeometry.cc:90
ME0Layer::etaPartition
const ME0EtaPartition * etaPartition(ME0DetId id) const
Return the eta partition corresponding to the given id.
Definition: ME0Layer.cc:27
ME0Geometry::chamber
const ME0Chamber * chamber(ME0DetId id) const
Return a chamber given its id.
Definition: ME0Geometry.cc:43
IdealMagneticFieldRecord
Definition: IdealMagneticFieldRecord.h:11
GeomDet::surface
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:37
L1TMuon::kRPC
Definition: L1TMuonSubsystems.h:5
L1TMuon::GeometryTranslator::calcRPCSpecificPhi
double calcRPCSpecificPhi(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:225
L1TMuon::GeometryTranslator::_geogem
edm::ESHandle< GEMGeometry > _geogem
Definition: GeometryTranslator.h:70
L1TMuon::kME0
Definition: L1TMuonSubsystems.h:5
edm::EventSetup::get
T get() const
Definition: EventSetup.h:87
L1TMuon::GeometryTranslator::_geocsc
edm::ESHandle< CSCGeometry > _geocsc
Definition: GeometryTranslator.h:68
ME0EtaPartition::centreOfStrip
LocalPoint centreOfStrip(int strip) const
Definition: ME0EtaPartition.cc:26
Surface::toGlobal
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
Definition: Surface.h:79
CSCPatternBank::getLegacyPosition
static double getLegacyPosition(int pattern)
Definition: CSCPatternBank.cc:31
L1TMuon::GeometryTranslator::calcGEMSpecificBend
double calcGEMSpecificBend(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:205
L1TMuon::GeometryTranslator::calcCSCSpecificEta
double calcCSCSpecificEta(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:295
relativeConstraints.geom
geom
Definition: relativeConstraints.py:72
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
Point3DBase< float, GlobalTag >
L1TMuon::GeometryTranslator::calcDTSpecificBend
double calcDTSpecificBend(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:351
L1TMuon::GeometryTranslator::calcGEMSpecificEta
double calcGEMSpecificEta(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:197
L1TMuon::GeometryTranslator::_magfield
edm::ESHandle< MagneticField > _magfield
Definition: GeometryTranslator.h:80
cmsswSequenceInfo.tp
tp
Definition: cmsswSequenceInfo.py:17
L1TMuon::GeometryTranslator::calcME0SpecificEta
double calcME0SpecificEta(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:171
phase1PixelTopology::layer
constexpr std::array< uint8_t, layerIndexSize > layer
Definition: phase1PixelTopology.h:99
RPCRoll::centreOfStrip
LocalPoint centreOfStrip(int strip) const
Definition: RPCRoll.cc:26
runTauDisplay.gp
gp
Definition: runTauDisplay.py:431
L1TMuon::kCSC
Definition: L1TMuonSubsystems.h:5
L1TMuon::kDT
Definition: L1TMuonSubsystems.h:5
L1TMuon::GeometryTranslator::getRPCSpecificPoint
GlobalPoint getRPCSpecificPoint(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:209
L1TMuon::GeometryTranslator::calcCSCSpecificBend
double calcCSCSpecificBend(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:303
ME0Layer
Definition: ME0Layer.h:9
GEMDetId
Definition: GEMDetId.h:18
CSCDetId
Definition: CSCDetId.h:26
L1TMuon::kGEM
Definition: L1TMuonSubsystems.h:5
PV3DBase::eta
T eta() const
Definition: PV3DBase.h:73
ME0Chamber
Definition: ME0Chamber.h:10
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:49
topSingleLeptonDQM_PU_cfi.pattern
pattern
Definition: topSingleLeptonDQM_PU_cfi.py:39
ME0DetId
Definition: ME0DetId.h:16
L1TMuon::GeometryTranslator::georpcToken_
edm::ESGetToken< RPCGeometry, MuonGeometryRecord > georpcToken_
Definition: GeometryTranslator.h:75
L1TMuon::GeometryTranslator::magfieldToken_
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > magfieldToken_
Definition: GeometryTranslator.h:81
DTBtiId
Definition: DTBtiId.h:32
L1TMuon::GeometryTranslator::calcDTSpecificPhi
double calcDTSpecificPhi(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:348
me0TriggerPseudoDigis_cff.nStrips
nStrips
1.2 is to make the matching window safely the two nearest strips 0.35 is the size of an ME0 chamber i...
Definition: me0TriggerPseudoDigis_cff.py:26
GEMEtaPartition::centreOfPad
LocalPoint centreOfPad(int pad) const
Definition: GEMEtaPartition.cc:45
L1TMuon::GeometryTranslator::calcME0SpecificBend
double calcME0SpecificBend(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:179
DDAxes::phi
PV3DBase::mag
T mag() const
Definition: PV3DBase.h:64
L1TMuon::GeometryTranslator::geome0Token_
edm::ESGetToken< ME0Geometry, MuonGeometryRecord > geome0Token_
Definition: GeometryTranslator.h:77
L1TMuon::GeometryTranslator::isCSCCounterClockwise
bool isCSCCounterClockwise(const std::unique_ptr< const CSCLayer > &) const
Definition: GeometryTranslator.cc:305
triggerObjects_cff.id
id
Definition: triggerObjects_cff.py:29
DTTrigGeom
Definition: DTTrigGeom.h:41
relativeConstraints.chamber
chamber
Definition: relativeConstraints.py:53
L1TMuon::GeometryTranslator::calcDTSpecificPoint
GlobalPoint calcDTSpecificPoint(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:314
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
DTChamberId
Definition: DTChamberId.h:14
hltrates_dqm_sourceclient-live_cfg.offset
offset
Definition: hltrates_dqm_sourceclient-live_cfg.py:82
L1TMuon::GeometryTranslator::calcRPCSpecificBend
double calcRPCSpecificBend(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:231
L1TMuon::GeometryTranslator::geocscToken_
edm::ESGetToken< CSCGeometry, MuonGeometryRecord > geocscToken_
Definition: GeometryTranslator.h:74
volumeBasedMagneticField_160812_cfi.magfield
magfield
Definition: volumeBasedMagneticField_160812_cfi.py:11
MuonGeometryRecord
Definition: MuonGeometryRecord.h:34
L1TMuon::GeometryTranslator::calcDTSpecificEta
double calcDTSpecificEta(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:346
PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:66
CSCGeometry::chamber
const CSCChamber * chamber(CSCDetId id) const
Return the chamber corresponding to given DetId.
Definition: CSCGeometry.cc:100
L1TMuon::GeometryTranslator::getGEMSpecificPoint
GlobalPoint getGEMSpecificPoint(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:185
L1TMuon::GeometryTranslator::_geom_cache_id
unsigned long long _geom_cache_id
Definition: GeometryTranslator.h:66
L1TMuon::GeometryTranslator::_geodt
edm::ESHandle< DTGeometry > _geodt
Definition: GeometryTranslator.h:67
L1TMuon::GeometryTranslator::geogemToken_
edm::ESGetToken< GEMGeometry, MuonGeometryRecord > geogemToken_
Definition: GeometryTranslator.h:76
ME0EtaPartition
Definition: ME0EtaPartition.h:12
GEMGeometry::etaPartition
const GEMEtaPartition * etaPartition(GEMDetId id) const
Return a GEMEtaPartition given its id.
Definition: GEMGeometry.cc:77