30 # define ADD_PIXEL_TOPOLOGY( rawid, detUnit ) \
31 const PixelGeomDetUnit* det = dynamic_cast<const PixelGeomDetUnit*>( detUnit ); \
34 const PixelTopology* topo = &det->specificTopology(); \
35 m_fwGeometry->idToName[rawid].topology[0] = topo->nrows(); \
36 m_fwGeometry->idToName[rawid].topology[1] = topo->ncolumns(); \
39 # define ADD_SISTRIP_TOPOLOGY( rawid, detUnit ) \
40 const StripGeomDetUnit* det = dynamic_cast<const StripGeomDetUnit*>( detUnit ); \
43 const StripTopology* topo = dynamic_cast<const StripTopology*>( &det->specificTopology() ); \
44 m_fwGeometry->idToName[rawid].topology[0] = 0; \
45 m_fwGeometry->idToName[rawid].topology[1] = topo->nstrips(); \
46 m_fwGeometry->idToName[rawid].topology[2] = topo->stripLength(); \
47 if( const RadialStripTopology* rtop = dynamic_cast<const RadialStripTopology*>( &(det->specificType().specificTopology()) ) ) \
49 m_fwGeometry->idToName[rawid].topology[0] = 1; \
50 m_fwGeometry->idToName[rawid].topology[3] = rtop->yAxisOrientation(); \
51 m_fwGeometry->idToName[rawid].topology[4] = rtop->originToIntersection(); \
52 m_fwGeometry->idToName[rawid].topology[5] = rtop->phiOfOneEdge(); \
53 m_fwGeometry->idToName[rawid].topology[6] = rtop->angularWidth(); \
55 else if( dynamic_cast<const RectangularStripTopology*>( &(det->specificType().specificTopology()) ) ) \
57 m_fwGeometry->idToName[rawid].topology[0] = 2; \
58 m_fwGeometry->idToName[rawid].topology[3] = topo->pitch(); \
60 else if( dynamic_cast<const TrapezoidalStripTopology*>( &(det->specificType().specificTopology()) ) ) \
62 m_fwGeometry->idToName[rawid].topology[0] = 3; \
63 m_fwGeometry->idToName[rawid].topology[3] = topo->pitch(); \
76 boost::shared_ptr<FWRecoGeometry>
115 for(
auto it = cscGeometry->
chambers().begin(),
129 for( std::vector< const CSCLayer* >::const_iterator lit = chamber->
layers().begin(),
130 lend = chamber->
layers().end();
138 unsigned int current =
insert_id( rawid );
167 for(
auto it = dtGeometry->
chambers().begin(),
182 for(
auto it = dtGeometry->
layers().begin(),
190 unsigned int rawid = layer->
id().
rawId();
220 for(
auto it = rpcGeom->
rolls().begin(),
242 m_fwGeometry->extraDet.Add(
new TNamed(
"RE4",
"RPC endcap station 4"));
262 unsigned int rawid = roll->geographicalId().rawId();
280 m_fwGeometry->extraDet.Add(
new TNamed(
"GEM",
"GEM muon detector"));
284 m_fwGeometry->extraDet.Add(
new TNamed(
"GE2",
"GEM endcap station 2"));
291 edm::LogInfo(
"FWRecoGeometry") <<
"failed to produce GEM geometry " << exception.
what() << std::endl;
310 unsigned int rawid = roll->geographicalId().rawId();
327 m_fwGeometry->extraDet.Add(
new TNamed(
"ME0",
"ME0 muon detector"));
331 edm::LogInfo(
"FWRecoGeometry") <<
"failed to produce ME0 geometry " << exception.
what() << std::endl;
347 unsigned int rawid = detid.
rawId();
368 unsigned int rawid = detid.
rawId();
389 unsigned int rawid = detid.
rawId();
410 unsigned int rawid = detid.
rawId();
431 unsigned int rawid = detid.
rawId();
452 unsigned int rawid = detid.
rawId();
464 std::vector<DetId> vid =
m_caloGeom->getValidDetIds();
465 for( std::vector<DetId>::const_iterator it = vid.begin(),
470 unsigned int id =
insert_id( it->rawId());
488 unsigned int index( 0 );
489 for( std::vector<GlobalPoint>::const_iterator
i = begin;
i !=
end; ++
i )
508 std::array< const float, 4 >
const & par = b2->parameters ();
521 m_fwGeometry->idToName[id].shape[1] = b2->width() * 0.5;
522 m_fwGeometry->idToName[id].shape[2] = b2->length() * 0.5;
523 m_fwGeometry->idToName[id].shape[3] = b2->thickness() * 0.5;
virtual char const * what() const
void swap(ora::Record &rh, ora::Record &lh)
virtual int nstrips() const =0
virtual const TrackerGeomDet * idToDetUnit(DetId) const
Return the pointer to the GeomDetUnit corresponding to a given DetId.
const std::vector< const CSCLayer * > & layers() const
Return all layers.
const std::vector< const DTChamber * > & chambers() const
Return a vector of all Chamber.
void addPixelForwardGeometry(void)
const CSCWireTopology * wireTopology() const
void addTIBGeometry(void)
void addTECGeometry(void)
float yCentreOfStripPlane() const
virtual ~FWRecoGeometryESProducer(void)
FWRecoGeometryESProducer(const edm::ParameterSet &)
DTLayerId id() const
Return the DetId of this SL.
const std::vector< const RPCRoll * > & rolls() const
Return a vector of all RPC rolls.
void addRPCGeometry(void)
float cellWidth() const
Returns the cell width.
edm::ESHandle< GlobalTrackingGeometry > m_geomRecord
const ChamberContainer & chambers() const
Return a vector of all chambers.
double wireSpacing() const
const StripTopology & specificTopology() const
const TrackerGeometry * m_trackerGeom
const Plane & surface() const
The nominal surface of the GeomDet.
boost::shared_ptr< FWRecoGeometry > m_fwGeometry
void setWhatProduced(T *iThis, const es::Label &iLabel=es::Label())
int firstChannel() const
Returns the wire number of the first wire.
uint32_t rawId() const
get the raw id
void fillShapeAndPlacement(unsigned int id, const GeomDet *det)
const DepRecordT & getRecord() const
void addPixelBarrelGeometry(void)
int lastChannel() const
Returns the wire number of the last wire.
const DTTopology & specificTopology() const
virtual float stripLength() const =0
const DetContainer & detsTEC() const
void addCaloGeometry(void)
void addCSCGeometry(void)
float yAxisOrientation() const
const DetContainer & detsPXB() const
void addGEMGeometry(void)
boost::shared_ptr< FWRecoGeometry > produce(const FWRecoGeometryRecord &)
const std::vector< const GEMEtaPartition * > & etaPartitions() const
Return a vector of all GEM eta partitions.
DetId geographicalId() const
The label of this GeomDet.
float centreToIntersection() const
float phiOfOneEdge() const
float cellHeight() const
Returns the cell height.
float angularWidth() const
#define ADD_PIXEL_TOPOLOGY(rawid, detUnit)
unsigned int insert_id(unsigned int id)
const DetContainer & detsTIB() const
void addTIDGeometry(void)
const std::vector< ME0EtaPartition const * > & etaPartitions() const
Return a vector of all ME0 eta partitions.
const CSCStripTopology * topology() const
void fillPoints(unsigned int id, std::vector< GlobalPoint >::const_iterator begin, std::vector< GlobalPoint >::const_iterator end)
virtual float stripOffset(void) const
int channels() const
Returns the number of wires in the layer.
const std::vector< const DTLayer * > & layers() const
Return a vector of all SuperLayer.
void addME0Geometry(void)
#define ADD_SISTRIP_TOPOLOGY(rawid, detUnit)
edm::ESHandle< CaloGeometry > m_caloGeom
virtual float pitch() const =0
const DetContainer & detsPXF() const
const DetContainer & detsTOB() const
const RotationType & rotation() const
const PositionType & position() const
const CSCLayerGeometry * geometry() const
void addTOBGeometry(void)
const DetContainer & detsTID() const
virtual const TrackerGeomDet * idToDet(DetId) const