DTGeometryBuilderFromCondDB Class Reference

#include <DTGeometryBuilderFromCondDB.h>

Public Member Functions
void	build (const std::shared_ptr< DTGeometry > &theGeometry, const RecoIdealGeometry &rig)
	DTGeometryBuilderFromCondDB ()
virtual	~DTGeometryBuilderFromCondDB ()

Private Types
typedef ReferenceCountingPointer< Plane >	RCPPlane

Private Member Functions
DTChamber *	buildChamber (const DetId &id, const RecoIdealGeometry &rig, size_t idt) const
DTLayer *	buildLayer (DTSuperLayer *sl, const DetId &id, const RecoIdealGeometry &rig, size_t idt) const
DTSuperLayer *	buildSuperLayer (DTChamber *chamber, const DetId &id, const RecoIdealGeometry &rig, size_t idt) const
RCPPlane	plane (const std::vector< double >::const_iterator tranStart, const std::vector< double >::const_iterator rotStart, Bounds *bounds) const

Detailed Description

Description:

detailed description

Author

: Stefano Lacaprara - INFN LNL stefa.nosp@m.no.l.nosp@m.acapr.nosp@m.ara@.nosp@m.pd.in.nosp@m.fn.i.nosp@m.t

Modification:

Definition at line 40 of file DTGeometryBuilderFromCondDB.h.

Member Typedef Documentation

RCPPlane

typedef ReferenceCountingPointer<Plane> DTGeometryBuilderFromCondDB::RCPPlane

private

Definition at line 58 of file DTGeometryBuilderFromCondDB.h.

Constructor & Destructor Documentation

DTGeometryBuilderFromCondDB()

DTGeometryBuilderFromCondDB::DTGeometryBuilderFromCondDB

(

)

Definition at line 36 of file DTGeometryBuilderFromCondDB.cc.

{}

~DTGeometryBuilderFromCondDB()

DTGeometryBuilderFromCondDB::~DTGeometryBuilderFromCondDB ( )

virtual

Definition at line 39 of file DTGeometryBuilderFromCondDB.cc.

{}

Member Function Documentation

build()

void DTGeometryBuilderFromCondDB::build	(	const std::shared_ptr< DTGeometry > &	theGeometry,
		const RecoIdealGeometry &	rig
	)

Definition at line 42 of file DTGeometryBuilderFromCondDB.cc.

References relativeConstraints::chamber, RecoIdealGeometry::detIds(), and RecoIdealGeometry::shapeStart().

Referenced by DTGeometryESModule::setupDBGeometry().

                                                                                                                  {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("DTGeometry") << "DTGeometryBuilderFromCondDB ";
#endif
  const std::vector<DetId>& detids(rig.detIds());
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("DTGeometry") << "size " << detids.size();
#endif

  size_t idt = 0;
  DTChamber* chamber(nullptr);
  DTSuperLayer* sl(nullptr);
  while (idt < detids.size()) {
    //copy(par.begin(), par.end(), ostream_iterator<double>(std::cout," "));
    if (int(*(rig.shapeStart(idt))) == 0) {  // a Chamber
      // add the provious chamber which by now has been updated with SL and
      // layers
      if (chamber)
        theGeometry->add(chamber);
      // go for the actual one
      DTChamberId chid(detids[idt]);
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("DTGeometry") << "CH: " << chid;
#endif
      chamber = buildChamber(chid, rig, idt);
    } else if (int(*(rig.shapeStart(idt))) == 1) {  // a SL
      DTSuperLayerId slid(detids[idt]);
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("DTGeometry") << "  SL: " << slid;
#endif
      sl = buildSuperLayer(chamber, slid, rig, idt);
      theGeometry->add(sl);
    } else if (int(*(rig.shapeStart(idt))) == 2) {  // a Layer
      DTLayerId lid(detids[idt]);
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("DTGeometry") << "    LAY: " << lid;
#endif
      DTLayer* lay = buildLayer(sl, lid, rig, idt);
      theGeometry->add(lay);
    } else {
      edm::LogVerbatim("DTGeometry") << "What is this?";
    }
    ++idt;
  }
  if (chamber)
    theGeometry->add(chamber);  // add the last chamber
}

buildChamber()

DTChamber * DTGeometryBuilderFromCondDB::buildChamber ( const DetId &  id, const RecoIdealGeometry &  rig, size_t  idt  ) const

private

SL the definition of length, width, thickness depends on the local reference frame of the Det

Definition at line 98 of file DTGeometryBuilderFromCondDB.cc.

References relativeConstraints::chamber, hcalRecHitTable_cff::detId, dtGeometryBuilder::getRecPlaneBounds(), RecoIdealGeometry::rotStart(), RecoIdealGeometry::shapeStart(), and RecoIdealGeometry::tranStart().

                                                                                                                    {
  DTChamberId detId(id);

  // width is along local X
  // length is along local Y
  // length z      dimension - constant 125.55 cm
  // thickness is along local Z
  // radial thickness - almost constant about 18 cm
  RCPPlane surf(
      plane(rig.tranStart(idt), rig.rotStart(idt), dtGeometryBuilder::getRecPlaneBounds(++rig.shapeStart(idt))));

  DTChamber* chamber = new DTChamber(detId, surf);

  return chamber;
}

buildLayer()

DTLayer * DTGeometryBuilderFromCondDB::buildLayer ( DTSuperLayer *  sl, const DetId &  id, const RecoIdealGeometry &  rig, size_t  idt  ) const

private

Definition at line 139 of file DTGeometryBuilderFromCondDB.cc.

References DTSuperLayer::add(), cms::cuda::assert(), angle_units::operators::convertMmToCm(), dtGeometryBuilder::getRecPlaneBounds(), DTSuperLayer::id(), RecoIdealGeometry::rotStart(), RecoIdealGeometry::shapeStart(), and RecoIdealGeometry::tranStart().

                                                                   {
  DTLayerId layId(id);

  // Layer specific parameter (size)
  // r-phi  dimension - different in different chambers
  // z      dimension - constant 126.8 cm
  // radial thickness - almost constant about 20 cm

  auto shapeStartPtr = rig.shapeStart(idt);
  RCPPlane surf(
      plane(rig.tranStart(idt), rig.rotStart(idt), dtGeometryBuilder::getRecPlaneBounds((shapeStartPtr + 1))));

  // Loop on wires
  int firstWire = static_cast<int>(*(shapeStartPtr + 4));       //par[4]);
  int WCounter = static_cast<int>(*(shapeStartPtr + 5));        //par[5]);
  double sensibleLength = convertMmToCm(*(shapeStartPtr + 6));  //par[6] in cm;
  DTTopology topology(firstWire, WCounter, sensibleLength);

  DTLayerType layerType;

  DTLayer* layer = new DTLayer(layId, surf, topology, layerType, sl);
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("DTGeometry") << "adding layer " << layer->id() << " to sl " << sl->id();
#endif
  assert(sl);
  sl->add(layer);
  return layer;
}

buildSuperLayer()

DTSuperLayer * DTGeometryBuilderFromCondDB::buildSuperLayer ( DTChamber *  chamber, const DetId &  id, const RecoIdealGeometry &  rig, size_t  idt  ) const

private

Definition at line 115 of file DTGeometryBuilderFromCondDB.cc.

References cms::cuda::assert(), relativeConstraints::chamber, dtGeometryBuilder::getRecPlaneBounds(), RecoIdealGeometry::rotStart(), RecoIdealGeometry::shapeStart(), and RecoIdealGeometry::tranStart().

                                                                             {
  DTSuperLayerId slId(id);

  // r-phi  dimension - different in different chambers
  // z      dimension - constant 126.8 cm
  // radial thickness - almost constant about 5 cm

  // Ok this is the slayer position...
  RCPPlane surf(
      plane(rig.tranStart(idt), rig.rotStart(idt), dtGeometryBuilder::getRecPlaneBounds(++rig.shapeStart(idt))));

  DTSuperLayer* slayer = new DTSuperLayer(slId, surf, chamber);

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("DTGeometry") << "adding slayer " << slayer->id() << " to chamber " << chamber->id();
#endif
  assert(chamber);
  chamber->add(slayer);
  return slayer;
}

plane()

DTGeometryBuilderFromCondDB::RCPPlane DTGeometryBuilderFromCondDB::plane ( const std::vector< double >::const_iterator  tranStart, const std::vector< double >::const_iterator  rotStart, Bounds *  bounds  ) const

private

Definition at line 171 of file DTGeometryBuilderFromCondDB.cc.

                                                                                               {
  // extract the position
  const Surface::PositionType posResult(*(tranStart), *(tranStart + 1), *(tranStart + 2));
  // now the rotation
  Surface::RotationType rotResult(*(rotStart + 0),
                                  *(rotStart + 1),
                                  *(rotStart + 2),
                                  *(rotStart + 3),
                                  *(rotStart + 4),
                                  *(rotStart + 5),
                                  *(rotStart + 6),
                                  *(rotStart + 7),
                                  *(rotStart + 8));

  return RCPPlane(new Plane(posResult, rotResult, bounds));
}