MTDGeomBuilderFromGeometricTimingDet Class Reference

#include <MTDGeomBuilderFromGeometricTimingDet.h>

Public Member Functions
MTDGeometry *	build (const GeometricTimingDet *gd, const PMTDParameters &ptp, const MTDTopology *tTopo)

Private Member Functions
void	buildGeomDet (MTDGeometry *)
void	buildPixel (std::vector< const GeometricTimingDet *> const &, MTDGeometry *, GeomDetType::SubDetector det, const PMTDParameters &ptp)
PlaneBuilderFromGeometricTimingDet::ResultType	buildPlaneWithMaterial (const GeometricTimingDet *gd, double scaleFactor=1.) const

Private Attributes
std::map< std::string, const MTDGeomDetType * >	theMTDDetTypeMap
const MTDTopology *	theTopo

Detailed Description

Definition at line 16 of file MTDGeomBuilderFromGeometricTimingDet.h.

Member Function Documentation

build()

MTDGeometry * MTDGeomBuilderFromGeometricTimingDet::build	(	const GeometricTimingDet *	gd,
		const PMTDParameters &	ptp,
		const MTDTopology *	tTopo
	)

Definition at line 36 of file MTDGeomBuilderFromGeometricTimingDet.cc.

References cms::cuda::assert(), buildGeomDet(), buildPixel(), AlCaHLTBitMon_QueryRunRegistry::comp, GeometricTimingDet::components(), TauDecayModes::dec, GeometricTimingDet::deepComponents(), mps_fire::i, LogDebug, LogTrace, MTDDetId::mtdSubDetector(), theMTDDetTypeMap, theTopo, GeomDetEnumerators::TimingBarrel, GeomDetEnumerators::TimingEndcap, PbPb_ZMuSkimMuonDPG_cff::tracker, and susybsm::HSCParticleType::unknown.

Referenced by MTDDigiGeometryESModule::produce().

                                                                                   {
  theMTDDetTypeMap.clear();

  MTDGeometry* tracker = new MTDGeometry(gd);
  std::vector<const GeometricTimingDet*> comp;
  gd->deepComponents(comp);

  if (tTopo)
    theTopo = tTopo;

  //define a vector which associate to the detid subdetector index -1 (from 0 to 5) the GeometridDet enumerator to be able to know which type of subdetector it is

  std::vector<GeometricTimingDet::GTDEnumType> gdsubdetmap(
      2, GeometricTimingDet::unknown);  // hardcoded "2" should not be a surprise...
  GeometricTimingDet::ConstGeometricTimingDetContainer subdetgd = gd->components();

  LogDebug("SubDetectorGeometricTimingDetType") << "MTD GeometricTimingDet enumerator values of the subdetectors";
  for (unsigned int i = 0; i < subdetgd.size(); ++i) {
    MTDDetId mtdid(subdetgd[i]->geographicalId());
    assert(mtdid.mtdSubDetector() > 0 && mtdid.mtdSubDetector() < 3);
    gdsubdetmap[mtdid.mtdSubDetector() - 1] = subdetgd[i]->type();
    LogTrace("SubDetectorGeometricTimingDetType") << "MTD subdet " << i << " type " << subdetgd[i]->type() << " detid "
                                                  << std::hex << subdetgd[i]->geographicalId().rawId() << std::dec
                                                  << " subdetid " << subdetgd[i]->geographicalId().subdetId();
  }

  std::vector<const GeometricTimingDet*> dets[2];
  std::vector<const GeometricTimingDet*>& btl = dets[0];
  btl.reserve(comp.size());
  std::vector<const GeometricTimingDet*>& etl = dets[1];
  etl.reserve(comp.size());

  for (auto& i : comp) {
    MTDDetId mtdid(i->geographicalId());
    dets[mtdid.mtdSubDetector() - 1].emplace_back(i);
  }

  //loop on all the six elements of dets and firstly check if they are from pixel-like detector and call buildPixel, then loop again and check if they are strip and call buildSilicon. "unknown" can be filled either way but the vector of GeometricTimingDet must be empty !!
  // this order is VERY IMPORTANT!!!!! For the moment I (AndreaV) understand that some pieces of code rely on pixel-like being before strip-like

  // now building the Pixel-like subdetectors
  for (unsigned int i = 0; i < 2; ++i) {
    if (gdsubdetmap[i] == GeometricTimingDet::BTL)
      buildPixel(dets[i], tracker, GeomDetEnumerators::SubDetector::TimingBarrel, ptp);
    if (gdsubdetmap[i] == GeometricTimingDet::ETL)
      buildPixel(dets[i], tracker, GeomDetEnumerators::SubDetector::TimingEndcap, ptp);
  }

  buildGeomDet(tracker);  //"GeomDet"

  verifyDUinTG(*tracker);

  return tracker;
}

buildGeomDet()

void MTDGeomBuilderFromGeometricTimingDet::buildGeomDet ( MTDGeometry *  tracker )

private

Definition at line 174 of file MTDGeomBuilderFromGeometricTimingDet.cc.

References mps_fire::i, and PbPb_ZMuSkimMuonDPG_cff::tracker.

Referenced by build().

                                                                            {
  auto const& gdu = tracker->detUnits();
  auto const& gduId = tracker->detUnitIds();

  for (u_int32_t i = 0; i < gdu.size(); i++) {
    tracker->addDet(gdu[i]);
    tracker->addDetId(gduId[i]);
    string gduTypeName = gdu[i]->type().name();
  }
}

buildPixel()

void MTDGeomBuilderFromGeometricTimingDet::buildPixel ( std::vector< const GeometricTimingDet *> const &  gdv, MTDGeometry *  tracker, GeomDetType::SubDetector  det, const PMTDParameters &  ptp  )

private

Definition at line 93 of file MTDGeomBuilderFromGeometricTimingDet.cc.

References MTDPixelTopologyBuilder::build(), buildPlaneWithMaterial(), Exception, mps_fire::i, LogDebug, AlCaHLTBitMon_QueryRunRegistry::string, GeomDetEnumerators::subDetGeom, submitPVValidationJobs::t, groupFilesInBlocks::temp, theMTDDetTypeMap, GeomDetEnumerators::TimingBarrel, GeomDetEnumerators::TimingEndcap, PbPb_ZMuSkimMuonDPG_cff::tracker, and PMTDParameters::vitems_.

Referenced by build().

{
  LogDebug("BuildingGeomDetUnits") << " MTD Pixel type. Size of vector: " << gdv.size()
                                   << " GeomDetType subdetector: " << det
                                   << " logical subdetector: " << GeomDetEnumerators::subDetGeom[det]
                                   << " big pix per ROC x: " << 0 << " y: " << 0 << " is upgrade: " << true;

  // this is a hack while we put things into the DDD
  int ROCrows(0), ROCcols(0), ROCSx(0), ROCSy(0);
  int GAPxInterpad(0), GAPxBorder(0), GAPyInterpad(0), GAPyBorder(0);
  switch (det) {
    case GeomDetType::SubDetector::TimingBarrel:
      GAPxInterpad = ptp.vitems_[0].vpars_[4];  // Value given in microns
      GAPxBorder = ptp.vitems_[0].vpars_[5];    // Value given in microns
      GAPyInterpad = ptp.vitems_[0].vpars_[6];  // Value given in microns
      GAPyBorder = ptp.vitems_[0].vpars_[7];    // Value given in microns
      ROCrows = ptp.vitems_[0].vpars_[8];
      ROCcols = ptp.vitems_[0].vpars_[9];
      ROCSx = ptp.vitems_[0].vpars_[10];
      ROCSy = ptp.vitems_[0].vpars_[11];
      break;
    case GeomDetType::SubDetector::TimingEndcap:
      GAPxInterpad = ptp.vitems_[1].vpars_[4];
      GAPxBorder = ptp.vitems_[1].vpars_[5];
      GAPyInterpad = ptp.vitems_[1].vpars_[6];
      GAPyBorder = ptp.vitems_[1].vpars_[7];
      ROCrows = ptp.vitems_[1].vpars_[8];
      ROCcols = ptp.vitems_[1].vpars_[9];
      ROCSx = ptp.vitems_[1].vpars_[10];
      ROCSy = ptp.vitems_[1].vpars_[11];
      break;
      break;
    default:
      throw cms::Exception("UnknownDet") << "MTDGeomBuilderFromGeometricTimingDet got a weird detector: " << det;
  }

  switch (det) {
    case GeomDetEnumerators::TimingBarrel:
      tracker->setOffsetDU(0);
      break;
    case GeomDetEnumerators::TimingEndcap:
      tracker->setOffsetDU(1);
      break;
    default:
      throw cms::Exception("MTDGeomBuilderFromGeometricTimingDet") << det << " is not a timing detector!";
  }

  for (auto i : gdv) {
    std::string const& detName = i->name();
    if (theMTDDetTypeMap.find(detName) == theMTDDetTypeMap.end()) {
      std::unique_ptr<const Bounds> bounds(i->bounds());

      PixelTopology* t = MTDPixelTopologyBuilder().build(
          &*bounds, ROCrows, ROCcols, ROCSx, ROCSy, GAPxInterpad, GAPxBorder, GAPyInterpad, GAPyBorder);

      theMTDDetTypeMap[detName] = new MTDGeomDetType(t, detName, det);
      tracker->addType(theMTDDetTypeMap[detName]);
    }

    PlaneBuilderFromGeometricTimingDet::ResultType plane = buildPlaneWithMaterial(i);
    GeomDetUnit* temp = new MTDGeomDetUnit(&(*plane), theMTDDetTypeMap[detName], i->geographicalID());

    tracker->addDetUnit(temp);
    tracker->addDetUnitId(i->geographicalID());
  }
  switch (det) {
    case GeomDetEnumerators::TimingBarrel:
      tracker->setEndsetDU(0);
      break;
    case GeomDetEnumerators::TimingEndcap:
      tracker->setEndsetDU(1);
      break;
    default:
      throw cms::Exception("MTDGeomBuilderFromGeometricTimingDet") << det << " is not a timing detector!";
  }
}

buildPlaneWithMaterial()

PlaneBuilderFromGeometricTimingDet::ResultType MTDGeomBuilderFromGeometricTimingDet::buildPlaneWithMaterial ( const GeometricTimingDet *  gd, double  scaleFactor = 1.  ) const

private

Definition at line 185 of file MTDGeomBuilderFromGeometricTimingDet.cc.

References PlaneBuilderFromGeometricTimingDet::plane(), GeometricTimingDet::radLength(), isotrackNtupler::scale, and GeometricTimingDet::xi().

Referenced by buildPixel().

                                                      {
  PlaneBuilderFromGeometricTimingDet planeBuilder;
  PlaneBuilderFromGeometricTimingDet::ResultType plane = planeBuilder.plane(gd);
  //
  // set medium properties (if defined)
  //
  plane->setMediumProperties(MediumProperties(gd->radLength() * scale, gd->xi() * scale));

  return plane;
}

Member Data Documentation

theMTDDetTypeMap

std::map<std::string, const MTDGeomDetType*> MTDGeomBuilderFromGeometricTimingDet::theMTDDetTypeMap

private

Definition at line 30 of file MTDGeomBuilderFromGeometricTimingDet.h.

Referenced by build(), and buildPixel().

theTopo

const MTDTopology* MTDGeomBuilderFromGeometricTimingDet::theTopo

private

Definition at line 31 of file MTDGeomBuilderFromGeometricTimingDet.h.

Referenced by build().