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