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

Detailed Description

Definition at line 41 of file GeometryTranslator.h.

Constructor & Destructor Documentation

◆ GeometryTranslator()

GeometryTranslator::GeometryTranslator ( )

Definition at line 24 of file GeometryTranslator.cc.

24 : _geom_cache_id(0ULL), _magfield_cache_id(0ULL) {}

◆ ~GeometryTranslator()

GeometryTranslator::~GeometryTranslator ( )

Definition at line 26 of file GeometryTranslator.cc.

26 {}

Member Function Documentation

◆ calcCSCSpecificBend()

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

Definition at line 294 of file GeometryTranslator.cc.

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

References cmsswSequenceInfo::tp.

Referenced by calculateBendAngle().

◆ calcCSCSpecificEta()

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

Definition at line 286 of file GeometryTranslator.cc.

286  {
287  return getCSCSpecificPoint(tp).eta();
288 }

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 290 of file GeometryTranslator.cc.

290  {
291  return getCSCSpecificPoint(tp).phi();
292 }

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 342 of file GeometryTranslator.cc.

342  {
343  int bend = tp.getDTData().bendingAngle;
344  double bendf = bend / 512.0;
345  return bendf;
346 }

References trklet::bend(), and 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 305 of file GeometryTranslator.cc.

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

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 196 of file GeometryTranslator.cc.

196 { return 0.0; }

Referenced by calculateBendAngle().

◆ calcGEMSpecificEta()

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

Definition at line 188 of file GeometryTranslator.cc.

188  {
189  return getGEMSpecificPoint(tp).eta();
190 }

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 192 of file GeometryTranslator.cc.

192  {
193  return getGEMSpecificPoint(tp).phi();
194 }

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 170 of file GeometryTranslator.cc.

170  {
171  return tp.getME0Data().deltaphi * (tp.getME0Data().bend == 0 ? 1 : -1);
172 }

References cmsswSequenceInfo::tp.

Referenced by calculateBendAngle().

◆ calcME0SpecificEta()

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

Definition at line 162 of file GeometryTranslator.cc.

162  {
163  return getME0SpecificPoint(tp).eta();
164 }

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 166 of file GeometryTranslator.cc.

166  {
167  return getME0SpecificPoint(tp).phi();
168 }

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 222 of file GeometryTranslator.cc.

222 { return 0.0; }

Referenced by calculateBendAngle().

◆ calcRPCSpecificEta()

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

Definition at line 212 of file GeometryTranslator.cc.

212  {
213  return getRPCSpecificPoint(tp).eta();
214 }

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 216 of file GeometryTranslator.cc.

216  {
217  return getRPCSpecificPoint(tp).phi();
218 }

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

Referenced by calculateGlobalPhi().

◆ calculateBendAngle()

double GeometryTranslator::calculateBendAngle ( const TriggerPrimitive tp) const

Definition at line 74 of file GeometryTranslator.cc.

74  {
75  switch (tp.subsystem()) {
76  case L1TMuon::kDT:
77  return calcDTSpecificBend(tp);
78  break;
79  case L1TMuon::kCSC:
80  return calcCSCSpecificBend(tp);
81  break;
82  case L1TMuon::kRPC:
83  return calcRPCSpecificBend(tp);
84  break;
85  case L1TMuon::kGEM:
86  return calcGEMSpecificBend(tp);
87  break;
88  case L1TMuon::kME0:
89  return calcME0SpecificBend(tp);
90  break;
91  default:
92  return std::nan("Invalid TP type!");
93  break;
94  }
95 }

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 28 of file GeometryTranslator.cc.

28  {
29  switch (tp.subsystem()) {
30  case L1TMuon::kDT:
31  return calcDTSpecificEta(tp);
32  break;
33  case L1TMuon::kCSC:
34  return calcCSCSpecificEta(tp);
35  break;
36  case L1TMuon::kRPC:
37  return calcRPCSpecificEta(tp);
38  break;
39  case L1TMuon::kGEM:
40  return calcGEMSpecificEta(tp);
41  break;
42  case L1TMuon::kME0:
43  return calcME0SpecificEta(tp);
44  break;
45  default:
46  return std::nan("Invalid TP type!");
47  break;
48  }
49 }

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 51 of file GeometryTranslator.cc.

51  {
52  switch (tp.subsystem()) {
53  case L1TMuon::kDT:
54  return calcDTSpecificPhi(tp);
55  break;
56  case L1TMuon::kCSC:
57  return calcCSCSpecificPhi(tp);
58  break;
59  case L1TMuon::kRPC:
60  return calcRPCSpecificPhi(tp);
61  break;
62  case L1TMuon::kGEM:
63  return calcGEMSpecificPhi(tp);
64  break;
65  case L1TMuon::kME0:
66  return calcME0SpecificPhi(tp);
67  break;
68  default:
69  return std::nan("Invalid TP type!");
70  break;
71  }
72 }

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 122 of file GeometryTranslator.cc.

122  {
124  unsigned long long geomid = geom.cacheIdentifier();
125  if (_geom_cache_id != geomid) {
126  geom.get(_geodt);
127  geom.get(_geocsc);
128  geom.get(_georpc);
129  geom.get(_geogem);
130  geom.get(_geome0);
131  _geom_cache_id = geomid;
132  }
133 
135  unsigned long long magfieldid = magfield.cacheIdentifier();
136  if (_magfield_cache_id != magfieldid) {
137  magfield.get(_magfield);
138  _magfield_cache_id = magfieldid;
139  }
140 }

References _geocsc, _geodt, _geogem, _geom_cache_id, _geome0, _georpc, _magfield, _magfield_cache_id, relativeConstraints::geom, edm::EventSetup::get(), and volumeBasedMagneticField_160812_cfi::magfield.

Referenced by EMTFSetup::reload().

◆ getCSCGeometry()

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

Definition at line 55 of file GeometryTranslator.h.

55 { return *_geocsc; }

References _geocsc.

◆ getCSCSpecificPoint()

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

Definition at line 231 of file GeometryTranslator.cc.

231  {
232  const CSCDetId id(tp.detId<CSCDetId>());
233  // we should change this to weak_ptrs at some point
234  // requires introducing std::shared_ptrs to geometry
235  std::unique_ptr<const CSCChamber> chamb(_geocsc->chamber(id));
236  assert(chamb != nullptr); // failed to get CSC chamber
237  std::unique_ptr<const CSCLayerGeometry> layer_geom(chamb->layer(CSCConstants::KEY_ALCT_LAYER)->geometry());
238  std::unique_ptr<const CSCLayer> layer(chamb->layer(CSCConstants::KEY_ALCT_LAYER));
239 
240  const uint16_t halfstrip = tp.getCSCData().strip;
241  const uint16_t pattern = tp.getCSCData().pattern;
242  const uint16_t keyWG = tp.getCSCData().keywire;
243  //const unsigned maxStrips = layer_geom->numberOfStrips();
244 
245  // so we can extend this later
246  // assume TMB2007 half-strips only as baseline
247  double offset = 0.0;
248  switch (1) {
249  case 1:
251  }
252  const unsigned halfstrip_offs = static_cast<unsigned>(0.5 + halfstrip + offset);
253  const unsigned strip = halfstrip_offs / 2 + 1; // geom starts from 1
254 
255  // the rough location of the hit at the ALCT key layer
256  // we will refine this using the half strip information
257  const LocalPoint& coarse_lp = layer_geom->stripWireGroupIntersection(strip, keyWG);
258  const GlobalPoint& coarse_gp = layer->surface().toGlobal(coarse_lp);
259 
260  // the strip width/4.0 gives the offset of the half-strip
261  // center with respect to the strip center
262  const double hs_offset = layer_geom->stripPhiPitch() / 4.0;
263 
264  // determine handedness of the chamber
265  const bool ccw = isCSCCounterClockwise(layer);
266  // we need to subtract the offset of even half strips and add the odd ones
267  const double phi_offset = ((halfstrip_offs % 2 ? 1 : -1) * (ccw ? -hs_offset : hs_offset));
268 
269  // the global eta calculation uses the middle of the strip
270  // so no need to increment it
271  const GlobalPoint final_gp(
272  GlobalPoint::Polar(coarse_gp.theta(), (coarse_gp.phi().value() + phi_offset), coarse_gp.mag()));
273 
274  // We need to add in some notion of the 'error' on trigger primitives
275  // like the width of the wire group by the width of the strip
276  // or something similar
277 
278  // release ownership of the pointers
279  chamb.release();
280  layer_geom.release();
281  layer.release();
282 
283  return final_gp;
284 }

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 54 of file GeometryTranslator.h.

54 { return *_geodt; }

References _geodt.

◆ getGEMGeometry()

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

Definition at line 57 of file GeometryTranslator.h.

57 { return *_geogem; }

References _geogem.

◆ getGEMSpecificPoint()

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

Definition at line 176 of file GeometryTranslator.cc.

176  {
177  const GEMDetId id(tp.detId<GEMDetId>());
178  const GEMEtaPartition* roll = _geogem->etaPartition(id);
179  assert(roll != nullptr); // failed to get GEM roll
180  //const uint16_t pad = tp.getGEMData().pad;
181  // Use half-pad precision, + 0.5 at the end to get the center of the pad (pad starts from 0)
182  const float pad = (0.5 * static_cast<float>(tp.getGEMData().pad_low + tp.getGEMData().pad_hi)) + 0.5f;
183  const LocalPoint& lp = roll->centreOfPad(pad);
184  const GlobalPoint& gp = roll->surface().toGlobal(lp);
185  return gp;
186 }

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 97 of file GeometryTranslator.cc.

97  {
98  switch (tp.subsystem()) {
99  case L1TMuon::kDT:
100  return calcDTSpecificPoint(tp);
101  break;
102  case L1TMuon::kCSC:
103  return getCSCSpecificPoint(tp);
104  break;
105  case L1TMuon::kRPC:
106  return getRPCSpecificPoint(tp);
107  break;
108  case L1TMuon::kGEM:
109  return getGEMSpecificPoint(tp);
110  break;
111  case L1TMuon::kME0:
112  return getME0SpecificPoint(tp);
113  break;
114  default:
116  GlobalPoint::Polar(std::nan("Invalid TP type!"), std::nan("Invalid TP type!"), std::nan("Invalid TP type!")));
117  return ret;
118  break;
119  }
120 }

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 60 of file GeometryTranslator.h.

60 { return *_magfield; }

References _magfield.

◆ getME0Geometry()

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

Definition at line 58 of file GeometryTranslator.h.

58 { return *_geome0; }

References _geome0.

◆ getME0SpecificPoint()

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

Definition at line 144 of file GeometryTranslator.cc.

144  {
145  const ME0DetId id(tp.detId<ME0DetId>());
146  const ME0Chamber* chamber = _geome0->chamber(id);
147  const ME0Layer* keylayer = chamber->layer(3); // ME0 key layer is layer 3
148  int partition = tp.getME0Data().partition; // 'partition' is in half-roll unit
149  int iroll = (partition >> 1) + 1;
150  const ME0EtaPartition* roll = keylayer->etaPartition(iroll);
151  assert(roll != nullptr); // failed to get ME0 roll
152  // See L1Trigger/ME0Trigger/src/ME0TriggerPseudoBuilder.cc
153  int phiposition = tp.getME0Data().phiposition; // 'phiposition' is in half-strip unit
154  int istrip = (phiposition >> 1);
155  int phiposition2 = (phiposition & 0x1);
156  float centreOfStrip = istrip + 0.25 + phiposition2 * 0.5;
157  const LocalPoint& lp = roll->centreOfStrip(centreOfStrip);
158  const GlobalPoint& gp = roll->surface().toGlobal(lp);
159  return gp;
160 }

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 56 of file GeometryTranslator.h.

56 { return *_georpc; }

References _georpc.

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

◆ getRPCSpecificPoint()

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

Definition at line 200 of file GeometryTranslator.cc.

200  {
201  const RPCDetId id(tp.detId<RPCDetId>());
202  const RPCRoll* roll = _georpc->roll(id);
203  assert(roll != nullptr); // failed to get RPC roll
204  //const int strip = static_cast<int>(tp.getRPCData().strip);
205  // Use half-strip precision, - 0.5 at the end to get the center of the strip (strip starts from 1)
206  const float strip = (0.5 * static_cast<float>(tp.getRPCData().strip_low + tp.getRPCData().strip_hi)) - 0.5f;
207  const LocalPoint& lp = roll->centreOfStrip(strip);
208  const GlobalPoint& gp = roll->surface().toGlobal(lp);
209  return gp;
210 }

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 296 of file GeometryTranslator.cc.

296  {
297  const int nStrips = layer->geometry()->numberOfStrips();
298  const double phi1 = layer->centerOfStrip(1).phi();
299  const double phiN = layer->centerOfStrip(nStrips).phi();
300  return ((std::abs(phi1 - phiN) < M_PI && phi1 >= phiN) || (std::abs(phi1 - phiN) >= M_PI && phi1 < phiN));
301 }

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 64 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 63 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 71 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry(), and getMagneticField().

◆ _magfield_cache_id

unsigned long long L1TMuon::GeometryTranslator::_magfield_cache_id
private

Definition at line 70 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().

runTheMatrix.ret
ret
prodAgent to be discontinued
Definition: runTheMatrix.py:542
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:144
L1TMuon::GeometryTranslator::_magfield_cache_id
unsigned long long _magfield_cache_id
Definition: GeometryTranslator.h:70
L1TMuon::GeometryTranslator::calcGEMSpecificPhi
double calcGEMSpecificPhi(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:192
L1TMuon::GeometryTranslator::calcME0SpecificPhi
double calcME0SpecificPhi(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:166
L1TMuon::GeometryTranslator::getCSCSpecificPoint
GlobalPoint getCSCSpecificPoint(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:231
digitizers_cfi.strip
strip
Definition: digitizers_cfi.py:19
L1TMuon::GeometryTranslator::calcCSCSpecificPhi
double calcCSCSpecificPhi(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:290
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:46
RPCGeometry::roll
const RPCRoll * roll(RPCDetId id) const
Return a roll given its id.
Definition: RPCGeometry.cc:50
L1TMuon::GeometryTranslator::_georpc
edm::ESHandle< RPCGeometry > _georpc
Definition: GeometryTranslator.h:66
L1TMuon::GeometryTranslator::_geome0
edm::ESHandle< ME0Geometry > _geome0
Definition: GeometryTranslator.h:68
L1TMuon::GeometryTranslator::calcRPCSpecificEta
double calcRPCSpecificEta(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:212
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:216
L1TMuon::GeometryTranslator::_geogem
edm::ESHandle< GEMGeometry > _geogem
Definition: GeometryTranslator.h:67
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:65
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:196
L1TMuon::GeometryTranslator::calcCSCSpecificEta
double calcCSCSpecificEta(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:286
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:342
L1TMuon::GeometryTranslator::calcGEMSpecificEta
double calcGEMSpecificEta(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:188
L1TMuon::GeometryTranslator::_magfield
edm::ESHandle< MagneticField > _magfield
Definition: GeometryTranslator.h:71
cmsswSequenceInfo.tp
tp
Definition: cmsswSequenceInfo.py:17
L1TMuon::GeometryTranslator::calcME0SpecificEta
double calcME0SpecificEta(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:162
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:200
L1TMuon::GeometryTranslator::calcCSCSpecificBend
double calcCSCSpecificBend(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:294
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
DTBtiId
Definition: DTBtiId.h:32
L1TMuon::GeometryTranslator::calcDTSpecificPhi
double calcDTSpecificPhi(const TriggerPrimitive &) const
Definition: GeometryTranslator.cc:339
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:170
DDAxes::phi
PV3DBase::mag
T mag() const
Definition: PV3DBase.h:64
L1TMuon::GeometryTranslator::isCSCCounterClockwise
bool isCSCCounterClockwise(const std::unique_ptr< const CSCLayer > &) const
Definition: GeometryTranslator.cc:296
trklet::bend
double bend(double r, double rinv, double stripPitch)
Definition: Util.h:166
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:305
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:222
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:337
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:176
L1TMuon::GeometryTranslator::_geom_cache_id
unsigned long long _geom_cache_id
Definition: GeometryTranslator.h:63
L1TMuon::GeometryTranslator::_geodt
edm::ESHandle< DTGeometry > _geodt
Definition: GeometryTranslator.h:64
ME0EtaPartition
Definition: ME0EtaPartition.h:12
GEMGeometry::etaPartition
const GEMEtaPartition * etaPartition(GEMDetId id) const
Return a GEMEtaPartition given its id.
Definition: GEMGeometry.cc:77