CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
FWRecoGeometryESProducer.cc
Go to the documentation of this file.
4 
29 
30 #include "TNamed.h"
31 # define ADD_PIXEL_TOPOLOGY( rawid, detUnit ) \
32  const PixelGeomDetUnit* det = dynamic_cast<const PixelGeomDetUnit*>( detUnit ); \
33  if( det ) \
34  { \
35  const PixelTopology* topo = &det->specificTopology(); \
36  m_fwGeometry->idToName[rawid].topology[0] = topo->nrows(); \
37  m_fwGeometry->idToName[rawid].topology[1] = topo->ncolumns(); \
38  } \
39 
40 # define ADD_SISTRIP_TOPOLOGY( rawid, detUnit ) \
41  const StripGeomDetUnit* det = dynamic_cast<const StripGeomDetUnit*>( detUnit ); \
42  if( det ) \
43  { \
44  const StripTopology* topo = dynamic_cast<const StripTopology*>( &det->specificTopology() ); \
45  m_fwGeometry->idToName[rawid].topology[0] = 0; \
46  m_fwGeometry->idToName[rawid].topology[1] = topo->nstrips(); \
47  m_fwGeometry->idToName[rawid].topology[2] = topo->stripLength(); \
48  if( const RadialStripTopology* rtop = dynamic_cast<const RadialStripTopology*>( &(det->specificType().specificTopology()) ) ) \
49  { \
50  m_fwGeometry->idToName[rawid].topology[0] = 1; \
51  m_fwGeometry->idToName[rawid].topology[3] = rtop->yAxisOrientation(); \
52  m_fwGeometry->idToName[rawid].topology[4] = rtop->originToIntersection(); \
53  m_fwGeometry->idToName[rawid].topology[5] = rtop->phiOfOneEdge(); \
54  m_fwGeometry->idToName[rawid].topology[6] = rtop->angularWidth(); \
55  } \
56  else if( dynamic_cast<const RectangularStripTopology*>( &(det->specificType().specificTopology()) ) ) \
57  { \
58  m_fwGeometry->idToName[rawid].topology[0] = 2; \
59  m_fwGeometry->idToName[rawid].topology[3] = topo->pitch(); \
60  } \
61  else if( dynamic_cast<const TrapezoidalStripTopology*>( &(det->specificType().specificTopology()) ) ) \
62  { \
63  m_fwGeometry->idToName[rawid].topology[0] = 3; \
64  m_fwGeometry->idToName[rawid].topology[3] = topo->pitch(); \
65  } \
66  } \
67 
68 namespace {
69  static const std::array<std::string,3> hgcal_geom_names = { { "HGCalEESensitive",
70  "HGCalHESiliconSensitive",
71  "HGCalHEScintillatorSensitive" } };
72 }
73 
75  : m_current( -1 )
76 {
77  setWhatProduced( this );
78 }
79 
81 {}
82 
83 boost::shared_ptr<FWRecoGeometry>
85 {
86  using namespace edm;
87 
88  m_fwGeometry = boost::shared_ptr<FWRecoGeometry>( new FWRecoGeometry );
89 
91 
92  DetId detId( DetId::Tracker, 0 );
93  m_trackerGeom = (const TrackerGeometry*) m_geomRecord->slaveGeometry( detId );
94 
95  record.getRecord<CaloGeometryRecord>().get( m_caloGeom );
97  for( const auto& name : hgcal_geom_names ) {
98  const auto& calogr = record.getRecord<CaloGeometryRecord>();
99  calogr.getRecord<IdealGeometryRecord>().get( name , test );
100  if( test.isValid() ) {
101  m_hgcalGeoms.push_back(test);
102  }
103  }
104 
107  addTIBGeometry();
108  addTIDGeometry();
109  addTOBGeometry();
110  addTECGeometry();
111  addDTGeometry();
112  addCSCGeometry();
113  addRPCGeometry();
114  addGEMGeometry();
115  addME0Geometry();
116  addCaloGeometry();
117 
118  m_fwGeometry->idToName.resize( m_current + 1 );
119  std::vector<FWRecoGeom::Info>( m_fwGeometry->idToName ).swap( m_fwGeometry->idToName );
120  std::sort( m_fwGeometry->idToName.begin(), m_fwGeometry->idToName.end());
121 
122  return m_fwGeometry;
123 }
124 
125 void
127 {
128  DetId detId( DetId::Muon, 2 );
129  const CSCGeometry* cscGeometry = (const CSCGeometry*) m_geomRecord->slaveGeometry( detId );
130  for( auto it = cscGeometry->chambers().begin(),
131  end = cscGeometry->chambers().end();
132  it != end; ++it )
133  {
134  const CSCChamber *chamber = *it;
135 
136  if( chamber )
137  {
138  unsigned int rawid = chamber->geographicalId();
139  unsigned int current = insert_id( rawid );
140  fillShapeAndPlacement( current, chamber );
141  //
142  // CSC layers geometry
143  //
144  for( std::vector< const CSCLayer* >::const_iterator lit = chamber->layers().begin(),
145  lend = chamber->layers().end();
146  lit != lend; ++lit )
147  {
148  const CSCLayer* layer = *lit;
149 
150  if( layer )
151  {
152  unsigned int rawid = layer->geographicalId();
153  unsigned int current = insert_id( rawid );
154  fillShapeAndPlacement( current, layer );
155 
156  const CSCStripTopology* stripTopology = layer->geometry()->topology();
157  m_fwGeometry->idToName[current].topology[0] = stripTopology->yAxisOrientation();
158  m_fwGeometry->idToName[current].topology[1] = stripTopology->centreToIntersection();
159  m_fwGeometry->idToName[current].topology[2] = stripTopology->yCentreOfStripPlane();
160  m_fwGeometry->idToName[current].topology[3] = stripTopology->phiOfOneEdge();
161  m_fwGeometry->idToName[current].topology[4] = stripTopology->stripOffset();
162  m_fwGeometry->idToName[current].topology[5] = stripTopology->angularWidth();
163 
164  const CSCWireTopology* wireTopology = layer->geometry()->wireTopology();
165  m_fwGeometry->idToName[current].topology[6] = wireTopology->wireSpacing();
166  m_fwGeometry->idToName[current].topology[7] = wireTopology->wireAngle();
167  }
168  }
169  }
170  }
171 }
172 
173 void
175 {
176  DetId detId( DetId::Muon, 1 );
177  const DTGeometry* dtGeometry = (const DTGeometry*) m_geomRecord->slaveGeometry( detId );
178 
179  //
180  // DT chambers geometry
181  //
182  for( auto it = dtGeometry->chambers().begin(),
183  end = dtGeometry->chambers().end();
184  it != end; ++it )
185  {
186  const DTChamber *chamber = *it;
187 
188  if( chamber )
189  {
190  unsigned int rawid = chamber->geographicalId().rawId();
191  unsigned int current = insert_id( rawid );
192  fillShapeAndPlacement( current, chamber );
193  }
194  }
195 
196  // Fill in DT layer parameters
197  for( auto it = dtGeometry->layers().begin(),
198  end = dtGeometry->layers().end();
199  it != end; ++it )
200  {
201  const DTLayer* layer = *it;
202 
203  if( layer )
204  {
205  unsigned int rawid = layer->id().rawId();
206  unsigned int current = insert_id( rawid );
207  fillShapeAndPlacement( current, layer );
208 
209  const DTTopology& topo = layer->specificTopology();
210  const BoundPlane& surf = layer->surface();
211  // Topology W/H/L:
212  m_fwGeometry->idToName[current].topology[0] = topo.cellWidth();
213  m_fwGeometry->idToName[current].topology[1] = topo.cellHeight();
214  m_fwGeometry->idToName[current].topology[2] = topo.cellLenght();
215  m_fwGeometry->idToName[current].topology[3] = topo.firstChannel();
216  m_fwGeometry->idToName[current].topology[4] = topo.lastChannel();
217  m_fwGeometry->idToName[current].topology[5] = topo.channels();
218 
219  // Bounds W/H/L:
220  m_fwGeometry->idToName[current].topology[6] = surf.bounds().width();
221  m_fwGeometry->idToName[current].topology[7] = surf.bounds().thickness();
222  m_fwGeometry->idToName[current].topology[8] = surf.bounds().length();
223  }
224  }
225 }
226 
227 void
229 {
230  //
231  // RPC rolls geometry
232  //
233  DetId detId( DetId::Muon, 3 );
234  const RPCGeometry* rpcGeom = (const RPCGeometry*) m_geomRecord->slaveGeometry( detId );
235  for( auto it = rpcGeom->rolls().begin(),
236  end = rpcGeom->rolls().end();
237  it != end; ++it )
238  {
239  const RPCRoll* roll = (*it);
240  if( roll )
241  {
242  unsigned int rawid = roll->geographicalId().rawId();
243  unsigned int current = insert_id( rawid );
244  fillShapeAndPlacement( current, roll );
245 
246  const StripTopology& topo = roll->specificTopology();
247  m_fwGeometry->idToName[current].topology[0] = topo.nstrips();
248  m_fwGeometry->idToName[current].topology[1] = topo.stripLength();
249  m_fwGeometry->idToName[current].topology[2] = topo.pitch();
250  }
251  }
252 
253 
254  try {
255  RPCDetId id(1, 1, 4, 1, 1, 1, 1 );
256  m_geomRecord->slaveGeometry( detId );
257  m_fwGeometry->extraDet.Add(new TNamed("RE4", "RPC endcap station 4"));
258  }
259  catch (...) {}
260 }
261 
262 void
264 {
265  //
266  // GEM geometry
267  //
268  DetId detId( DetId::Muon, 4 );
269 
270  try
271  {
272  const GEMGeometry* gemGeom = (const GEMGeometry*) m_geomRecord->slaveGeometry( detId );
273  for(auto roll : gemGeom->etaPartitions())
274  {
275  if( roll )
276  {
277  unsigned int rawid = roll->geographicalId().rawId();
278  unsigned int current = insert_id( rawid );
279  fillShapeAndPlacement( current, roll );
280 
281  const StripTopology& topo = roll->specificTopology();
282  m_fwGeometry->idToName[current].topology[0] = topo.nstrips();
283  m_fwGeometry->idToName[current].topology[1] = topo.stripLength();
284  m_fwGeometry->idToName[current].topology[2] = topo.pitch();
285 
286  float height = topo.stripLength()/2;
287  LocalPoint lTop( 0., height, 0.);
288  LocalPoint lBottom( 0., -height, 0.);
289  m_fwGeometry->idToName[current].topology[3] = roll->localPitch(lTop);
290  m_fwGeometry->idToName[current].topology[4] = roll->localPitch(lBottom);
291  m_fwGeometry->idToName[current].topology[5] = roll->npads();
292  }
293  }
294 
295  m_fwGeometry->extraDet.Add(new TNamed("GEM", "GEM muon detector"));
296  try {
297  GEMDetId id(1, 1, 2, 1, 1, 1 );
298  m_geomRecord->slaveGeometry( detId );
299  m_fwGeometry->extraDet.Add(new TNamed("GE2", "GEM endcap station 2"));
300  }
301  catch (...) {}
302 
303  }
304  catch( cms::Exception &exception )
305  {
306  edm::LogInfo("FWRecoGeometry") << "failed to produce GEM geometry " << exception.what() << std::endl;
307  }
308 }
309 
310 void
312 {
313  //
314  // ME0 geometry
315  //
316 
317  DetId detId( DetId::Muon, 5 );
318  try
319  {
320  const ME0Geometry* me0Geom = (const ME0Geometry*) m_geomRecord->slaveGeometry( detId );
321  for(auto roll : me0Geom->etaPartitions())
322  {
323  if( roll )
324  {
325  unsigned int rawid = roll->geographicalId().rawId();
326  unsigned int current = insert_id( rawid );
327  fillShapeAndPlacement( current, roll );
328 
329  const StripTopology& topo = roll->specificTopology();
330  m_fwGeometry->idToName[current].topology[0] = topo.nstrips();
331  m_fwGeometry->idToName[current].topology[1] = topo.stripLength();
332  m_fwGeometry->idToName[current].topology[2] = topo.pitch();
333 
334  float height = topo.stripLength()/2;
335  LocalPoint lTop( 0., height, 0.);
336  LocalPoint lBottom( 0., -height, 0.);
337  m_fwGeometry->idToName[current].topology[3] = roll->localPitch(lTop);
338  m_fwGeometry->idToName[current].topology[4] = roll->localPitch(lBottom);
339  m_fwGeometry->idToName[current].topology[5] = roll->npads();
340  }
341  }
342  m_fwGeometry->extraDet.Add(new TNamed("ME0", "ME0 muon detector"));
343  }
344  catch( cms::Exception &exception )
345  {
346  edm::LogInfo("FWRecoGeometry") << "failed to produce ME0 geometry " << exception.what() << std::endl;
347  }
348 }
349 
350 void
352 {
353  for( TrackerGeometry::DetContainer::const_iterator it = m_trackerGeom->detsPXB().begin(),
354  end = m_trackerGeom->detsPXB().end();
355  it != end; ++it)
356  {
357  const GeomDet *det = *it;
358 
359  if( det )
360  {
361  DetId detid = det->geographicalId();
362  unsigned int rawid = detid.rawId();
363  unsigned int current = insert_id( rawid );
364  fillShapeAndPlacement( current, det );
365 
366  ADD_PIXEL_TOPOLOGY( current, m_trackerGeom->idToDetUnit( detid ));
367  }
368  }
369 }
370 
371 void
373 {
374  for( TrackerGeometry::DetContainer::const_iterator it = m_trackerGeom->detsPXF().begin(),
375  end = m_trackerGeom->detsPXF().end();
376  it != end; ++it )
377  {
378  const GeomDet *det = *it;
379 
380  if( det )
381  {
382  DetId detid = det->geographicalId();
383  unsigned int rawid = detid.rawId();
384  unsigned int current = insert_id( rawid );
385  fillShapeAndPlacement( current, det );
386 
387  ADD_PIXEL_TOPOLOGY( current, m_trackerGeom->idToDetUnit( detid ));
388  }
389  }
390 }
391 
392 void
394 {
395  for( TrackerGeometry::DetContainer::const_iterator it = m_trackerGeom->detsTIB().begin(),
396  end = m_trackerGeom->detsTIB().end();
397  it != end; ++it )
398  {
399  const GeomDet *det = *it;
400 
401  if( det )
402  {
403  DetId detid = det->geographicalId();
404  unsigned int rawid = detid.rawId();
405  unsigned int current = insert_id( rawid );
406  fillShapeAndPlacement( current, det );
407 
408  ADD_SISTRIP_TOPOLOGY( current, m_trackerGeom->idToDet( detid ));
409  }
410  }
411 }
412 
413 void
415 {
416  for( TrackerGeometry::DetContainer::const_iterator it = m_trackerGeom->detsTOB().begin(),
417  end = m_trackerGeom->detsTOB().end();
418  it != end; ++it )
419  {
420  const GeomDet *det = *it;
421 
422  if( det )
423  {
424  DetId detid = det->geographicalId();
425  unsigned int rawid = detid.rawId();
426  unsigned int current = insert_id( rawid );
427  fillShapeAndPlacement( current, det );
428 
429  ADD_SISTRIP_TOPOLOGY( current, m_trackerGeom->idToDet( detid ));
430  }
431  }
432 }
433 
434 void
436 {
437  for( TrackerGeometry::DetContainer::const_iterator it = m_trackerGeom->detsTID().begin(),
438  end = m_trackerGeom->detsTID().end();
439  it != end; ++it)
440  {
441  const GeomDet *det = *it;
442 
443  if( det )
444  {
445  DetId detid = det->geographicalId();
446  unsigned int rawid = detid.rawId();
447  unsigned int current = insert_id( rawid );
448  fillShapeAndPlacement( current, det );
449 
450  ADD_SISTRIP_TOPOLOGY( current, m_trackerGeom->idToDet( detid ));
451  }
452  }
453 }
454 
455 void
457 {
458  for( TrackerGeometry::DetContainer::const_iterator it = m_trackerGeom->detsTEC().begin(),
459  end = m_trackerGeom->detsTEC().end();
460  it != end; ++it )
461  {
462  const GeomDet *det = *it;
463 
464  if( det )
465  {
466  DetId detid = det->geographicalId();
467  unsigned int rawid = detid.rawId();
468  unsigned int current = insert_id( rawid );
469  fillShapeAndPlacement( current, det );
470 
471  ADD_SISTRIP_TOPOLOGY( current, m_trackerGeom->idToDet( detid ));
472  }
473  }
474 }
475 
476 void
478 {
479  std::vector<DetId> vid = m_caloGeom->getValidDetIds(); // Calo
480  for( std::vector<DetId>::const_iterator it = vid.begin(),
481  end = vid.end();
482  it != end; ++it )
483  {
484  const CaloCellGeometry::CornersVec& cor( m_caloGeom->getGeometry( *it )->getCorners());
485  unsigned int id = insert_id( it->rawId());
486  fillPoints( id, cor.begin(), cor.end());
487  }
488  for( const auto& hgc : m_hgcalGeoms ) {
489  const auto& hgcprod = *hgc;
490  const auto& vid = hgcprod.getValidDetIds();
491  for( const auto& id : vid ) {
492  uint32_t intid = insert_id( id.rawId() );
493  const auto& cor = hgcprod.getCorners( id );
494  fillPoints( intid, cor.begin(), cor.end() );
495  }
496  }
497 }
498 
499 unsigned int
501 {
502  ++m_current;
503  m_fwGeometry->idToName.push_back(FWRecoGeom::Info());
504  m_fwGeometry->idToName.back().id = rawid;
505 
506  return m_current;
507 }
508 
509 void
510 FWRecoGeometryESProducer::fillPoints( unsigned int id, std::vector<GlobalPoint>::const_iterator begin, std::vector<GlobalPoint>::const_iterator end )
511 {
512  unsigned int index( 0 );
513  for( std::vector<GlobalPoint>::const_iterator i = begin; i != end; ++i )
514  {
515  assert( index < 23 );
516  m_fwGeometry->idToName[id].points[index] = i->x();
517  m_fwGeometry->idToName[id].points[++index] = i->y();
518  m_fwGeometry->idToName[id].points[++index] = i->z();
519  ++index;
520  }
521 }
522 
523 
525 void
527 {
528  // Trapezoidal
529  const Bounds *b = &((det->surface ()).bounds ());
530  if( const TrapezoidalPlaneBounds *b2 = dynamic_cast<const TrapezoidalPlaneBounds *> (b))
531  {
532  std::array< const float, 4 > const & par = b2->parameters ();
533 
534  // These parameters are half-lengths, as in CMSIM/GEANT3
535  m_fwGeometry->idToName[id].shape[0] = 1;
536  m_fwGeometry->idToName[id].shape[1] = par [0]; // hBottomEdge
537  m_fwGeometry->idToName[id].shape[2] = par [1]; // hTopEdge
538  m_fwGeometry->idToName[id].shape[3] = par [2]; // thickness
539  m_fwGeometry->idToName[id].shape[4] = par [3]; // apothem
540  }
541  if( const RectangularPlaneBounds *b2 = dynamic_cast<const RectangularPlaneBounds *> (b))
542  {
543  // Rectangular
544  m_fwGeometry->idToName[id].shape[0] = 2;
545  m_fwGeometry->idToName[id].shape[1] = b2->width() * 0.5; // half width
546  m_fwGeometry->idToName[id].shape[2] = b2->length() * 0.5; // half length
547  m_fwGeometry->idToName[id].shape[3] = b2->thickness() * 0.5; // half thickness
548  }
549 
550  // Position of the DetUnit's center
551  GlobalPoint pos = det->surface().position();
552  m_fwGeometry->idToName[id].translation[0] = pos.x();
553  m_fwGeometry->idToName[id].translation[1] = pos.y();
554  m_fwGeometry->idToName[id].translation[2] = pos.z();
555 
556  // Add the coeff of the rotation matrix
557  // with a projection on the basis vectors
558  TkRotation<float> detRot = det->surface().rotation();
559  m_fwGeometry->idToName[id].matrix[0] = detRot.xx();
560  m_fwGeometry->idToName[id].matrix[1] = detRot.yx();
561  m_fwGeometry->idToName[id].matrix[2] = detRot.zx();
562  m_fwGeometry->idToName[id].matrix[3] = detRot.xy();
563  m_fwGeometry->idToName[id].matrix[4] = detRot.yy();
564  m_fwGeometry->idToName[id].matrix[5] = detRot.zy();
565  m_fwGeometry->idToName[id].matrix[6] = detRot.xz();
566  m_fwGeometry->idToName[id].matrix[7] = detRot.yz();
567  m_fwGeometry->idToName[id].matrix[8] = detRot.zz();
568 }
T xx() const
virtual char const * what() const
Definition: Exception.cc:141
void swap(ora::Record &rh, ora::Record &lh)
Definition: Record.h:70
virtual int nstrips() const =0
virtual const TrackerGeomDet * idToDetUnit(DetId) const
Return the pointer to the GeomDetUnit corresponding to a given DetId.
int i
Definition: DBlmapReader.cc:9
const std::vector< const CSCLayer * > & layers() const
Return all layers.
Definition: CSCChamber.h:57
const std::vector< const DTChamber * > & chambers() const
Return a vector of all Chamber.
Definition: DTGeometry.cc:85
const CSCWireTopology * wireTopology() const
JetCorrectorParameters::Record record
Definition: classes.h:7
FWRecoGeometryESProducer(const edm::ParameterSet &)
DTLayerId id() const
Return the DetId of this SL.
Definition: DTLayer.cc:46
const std::vector< const RPCRoll * > & rolls() const
Return a vector of all RPC rolls.
Definition: RPCGeometry.cc:67
assert(m_qm.get())
T y() const
Definition: PV3DBase.h:63
T yx() const
float cellWidth() const
Returns the cell width.
Definition: DTTopology.h:68
edm::ESHandle< GlobalTrackingGeometry > m_geomRecord
const ChamberContainer & chambers() const
Return a vector of all chambers.
Definition: CSCGeometry.cc:106
double wireSpacing() const
const StripTopology & specificTopology() const
Definition: RPCRoll.cc:107
const TrackerGeometry * m_trackerGeom
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:40
boost::shared_ptr< FWRecoGeometry > m_fwGeometry
void setWhatProduced(T *iThis, const es::Label &iLabel=es::Label())
Definition: ESProducer.h:115
int firstChannel() const
Returns the wire number of the first wire.
Definition: DTTopology.h:78
uint32_t rawId() const
get the raw id
Definition: DetId.h:43
T zx() const
void fillShapeAndPlacement(unsigned int id, const GeomDet *det)
T xy() const
int lastChannel() const
Returns the wire number of the last wire.
Definition: DTTopology.h:80
const DTTopology & specificTopology() const
Definition: DTLayer.cc:42
T zz() const
virtual float stripLength() const =0
float wireAngle() const
const DetContainer & detsTEC() const
T z() const
Definition: PV3DBase.h:64
const DetContainer & detsPXB() const
boost::shared_ptr< FWRecoGeometry > produce(const FWRecoGeometryRecord &)
T zy() const
const std::vector< const GEMEtaPartition * > & etaPartitions() const
Return a vector of all GEM eta partitions.
Definition: GEMGeometry.cc:63
DetId geographicalId() const
The label of this GeomDet.
Definition: GeomDet.h:77
#define end
Definition: vmac.h:37
float cellHeight() const
Returns the cell height.
Definition: DTTopology.h:70
T yy() const
#define ADD_PIXEL_TOPOLOGY(rawid, detUnit)
unsigned int insert_id(unsigned int id)
const DetContainer & detsTIB() const
const std::vector< ME0EtaPartition const * > & etaPartitions() const
Return a vector of all ME0 eta partitions.
Definition: ME0Geometry.cc:57
Definition: DetId.h:18
const CSCStripTopology * topology() const
void fillPoints(unsigned int id, std::vector< GlobalPoint >::const_iterator begin, std::vector< GlobalPoint >::const_iterator end)
std::vector< edm::ESHandle< HGCalGeometry > > m_hgcalGeoms
virtual float stripOffset(void) const
int channels() const
Returns the number of wires in the layer.
Definition: DTTopology.h:75
const std::vector< const DTLayer * > & layers() const
Return a vector of all SuperLayer.
Definition: DTGeometry.cc:95
#define ADD_SISTRIP_TOPOLOGY(rawid, detUnit)
double b
Definition: hdecay.h:120
edm::ESHandle< CaloGeometry > m_caloGeom
#define begin
Definition: vmac.h:30
T xz() const
virtual float pitch() const =0
const DetContainer & detsPXF() const
const DetContainer & detsTOB() const
const RotationType & rotation() const
Definition: Bounds.h:22
float cellLenght() const
Definition: DTTopology.h:73
T x() const
Definition: PV3DBase.h:62
const PositionType & position() const
T yz() const
const CSCLayerGeometry * geometry() const
Definition: CSCLayer.h:47
const DetContainer & detsTID() const
virtual const TrackerGeomDet * idToDet(DetId) const