DTGeometry.cc

Go to the documentation of this file.

#include <Geometry/DTGeometry/interface/DTGeometry.h>
#include <Geometry/CommonDetUnit/interface/GeomDet.h>

#include <algorithm>
#include <iostream>

DTGeometry::DTGeometry() {}

DTGeometry::~DTGeometry() { deallocate(); }

void DTGeometry::deallocate() {
  // delete all the chambers (which will delete the SL which will delete the
  // layers)
  for (auto& theChamber : theChambers)
    delete theChamber;
}

void DTGeometry::clear() {
  deallocate();

  theChambers.clear();
  theSuperLayers.clear();
  theLayers.clear();
  theMap.clear();
  theDetUnits.clear();
  theDets.clear();
  theDetTypes.clear();
  theDetUnitIds.clear();
  theDetIds.clear();
}

const DTGeometry::DetTypeContainer& DTGeometry::detTypes() const {
  // FIXME - fill it at runtime
  return theDetTypes;
}

void DTGeometry::add(DTChamber* ch) {
  theDets.emplace_back(ch);
  theChambers.emplace_back(ch);
  theMap.insert(DTDetMap::value_type(ch->geographicalId(), ch));
}

void DTGeometry::add(DTSuperLayer* sl) {
  theDets.emplace_back(sl);
  theSuperLayers.emplace_back(sl);
  theMap.insert(DTDetMap::value_type(sl->geographicalId(), sl));
}

void DTGeometry::add(DTLayer* l) {
  theDetUnits.emplace_back(l);
  theDets.emplace_back(l);
  theLayers.emplace_back(l);
  theMap.insert(DTDetMap::value_type(l->geographicalId(), l));
}

const DTGeometry::DetContainer& DTGeometry::detUnits() const { return theDetUnits; }

const DTGeometry::DetContainer& DTGeometry::dets() const { return theDets; }

const DTGeometry::DetIdContainer& DTGeometry::detUnitIds() const {
  // FIXME - fill it at runtime
  return theDetUnitIds;
}

const DTGeometry::DetIdContainer& DTGeometry::detIds() const {
  // FIXME - fill it at runtime
  return theDetIds;
}

const GeomDet* DTGeometry::idToDetUnit(DetId id) const { return dynamic_cast<const GeomDet*>(idToDet(id)); }

const GeomDet* DTGeometry::idToDet(DetId id) const {
  // Strip away wire#, if any!
  DTLayerId lId(id.rawId());
  DTDetMap::const_iterator i = theMap.find(lId);
  return (i != theMap.end()) ? i->second : nullptr;
}

const std::vector<const DTChamber*>& DTGeometry::chambers() const { return theChambers; }

const std::vector<const DTSuperLayer*>& DTGeometry::superLayers() const { return theSuperLayers; }

const std::vector<const DTLayer*>& DTGeometry::layers() const { return theLayers; }

const DTChamber* DTGeometry::chamber(const DTChamberId& id) const { return (const DTChamber*)(idToDet(id)); }

const DTSuperLayer* DTGeometry::superLayer(const DTSuperLayerId& id) const {
  return (const DTSuperLayer*)(idToDet(id));
}

const DTLayer* DTGeometry::layer(const DTLayerId& id) const { return (const DTLayer*)(idToDet(id)); }