MTDTopology Class Reference

#include <MTDTopology.h>

Classes
struct	ETLfaceLayout

Public Types
using	ETLValues = std::vector< ETLfaceLayout >

Public Member Functions
int	getMTDTopologyMode () const
size_t	hshiftETL (const uint32_t detid, const int horizontalShift) const
	MTDTopology (const int &topologyMode, const ETLValues &etl)
size_t	vshiftETL (const uint32_t detid, const int verticalShift, size_t &closest) const

Static Public Member Functions
static bool	orderETLSector (const GeomDet *&gd1, const GeomDet *&gd2)

Private Attributes
const ETLValues	etlVals_
const int	mtdTopologyMode_

Static Private Attributes
static constexpr size_t	failIndex_

Detailed Description

Definition at line 13 of file MTDTopology.h.

Member Typedef Documentation

ETLValues

using MTDTopology::ETLValues = std::vector<ETLfaceLayout>

Definition at line 23 of file MTDTopology.h.

Constructor & Destructor Documentation

MTDTopology()

MTDTopology::MTDTopology	(	const int &	topologyMode,
		const ETLValues &	etl
	)

Definition at line 5 of file MTDTopology.cc.

    : mtdTopologyMode_(topologyMode), etlVals_(etl) {}

Member Function Documentation

getMTDTopologyMode()

int MTDTopology::getMTDTopologyMode ( ) const

inline

Definition at line 27 of file MTDTopology.h.

References mtdTopologyMode_.

Referenced by ETLDetLayerGeometryBuilder::buildLayers(), BTLDetLayerGeometryBuilder::buildLayers(), MTDThresholdClusterizer::clusterize(), MTDThresholdClusterizer::copy_to_buffer(), MtdTruthAccumulator::finalizeEvent(), BTLDeviceSim::getHitsResponse(), MTDTimeCalib::getTimeCalib(), mtd::MTDGeomUtil::globalPosition(), hshiftETL(), mtd::MTDGeomUtil::pixelInModule(), mtd::MTDGeomUtil::position(), mtd::MTDGeomUtil::sensorModuleId(), and vshiftETL().

{ return mtdTopologyMode_; }

hshiftETL()

size_t MTDTopology::hshiftETL	(	const uint32_t	detid,
		const int	horizontalShift
	)		const

Definition at line 23 of file MTDTopology.cc.

References MTDTopologyMode::btlv3etlv8, ALCARECOPPSCalTrackBasedSel_cff::detid, ETLDetId::discSide(), etlVals_, failIndex_, getMTDTopologyMode(), ETLDetId::modType(), callgraph::module, ETLDetId::module(), and ETLDetId::sensor().

Referenced by MTDDetSector::hshift().

                                                                                   {
  ETLDetId start_mod(detid);

  if (horizontalShift == 0) {
    edm::LogWarning("MTDTopology") << "asking of a null horizotalShift in ETL";
    return failIndex_;
  }
  int hsh = horizontalShift > 0 ? 1 : -1;

  int sensor = start_mod.sensor();
  int module = start_mod.module();
  uint32_t modtyp = start_mod.modType();
  uint32_t discside = start_mod.discSide();
  int geomDetIndex;

  // distinguish numbering in prev8 / v8 geometries
  if (getMTDTopologyMode() == static_cast<int>(MTDTopologyMode::Mode::btlv3etlv8)) {
    geomDetIndex = 2 * (module - 1) + sensor;
  } else {
    geomDetIndex = module;
  }

  // ilayout number coincides at present with disc face, use this

  size_t iHome = (modtyp == etlVals_[discside].idDetType1_) ? 0 : 1;
  size_t iLeft = (etlVals_[discside].idDetType1_ == 1) ? 0 : 1;

  // for left type modules the position according to the default order is module - 1, for the rigth type modules the total number of left modules must be added

  size_t nmodOffset = (modtyp == 1) ? 0 : etlVals_[discside].start_copy_[iLeft].back() - 1;

  for (size_t iloop = 0; iloop < etlVals_[discside].start_copy_[iHome].size() - 1; iloop++) {
    if (geomDetIndex >= etlVals_[discside].start_copy_[iHome][iloop] &&
        geomDetIndex < etlVals_[discside].start_copy_[iHome][iloop + 1]) {
      if (geomDetIndex + hsh >= etlVals_[discside].start_copy_[iHome][iloop] &&
          geomDetIndex + hsh < etlVals_[discside].start_copy_[iHome][iloop + 1]) {
        return geomDetIndex + hsh - 1 + nmodOffset;
      }
      break;
    }
  }

  return failIndex_;
}

orderETLSector()

bool MTDTopology::orderETLSector ( const GeomDet *&  gd1, const GeomDet *&  gd2  )

static

Definition at line 8 of file MTDTopology.cc.

References GeomDet::geographicalId(), MTDDetId::mtdRR(), and DetId::rawId().

Referenced by ETLDetLayerGeometryBuilder::buildLayer().

                                                                         {
  ETLDetId det1(gd1->geographicalId().rawId());
  ETLDetId det2(gd2->geographicalId().rawId());

  if (det1.mtdRR() != det2.mtdRR()) {
    return det1.mtdRR() < det2.mtdRR();
  } else if (det1.modType() != det2.modType()) {
    return det1.modType() < det2.modType();
  } else if (det1.module() != det2.module()) {
    return det1.module() < det2.module();
  } else {
    return det1.sensor() < det2.sensor();
  }
}

vshiftETL()

size_t MTDTopology::vshiftETL	(	const uint32_t	detid,
		const int	verticalShift,
		size_t &	closest
	)		const

Definition at line 68 of file MTDTopology.cc.

References MTDTopologyMode::btlv3etlv8, pv::closest(), ALCARECOPPSCalTrackBasedSel_cff::detid, ETLDetId::discSide(), etlVals_, failIndex_, getMTDTopologyMode(), ETLDetId::modType(), callgraph::module, ETLDetId::module(), ETLDetId::sensor(), and findQualityFiles::size.

Referenced by MTDDetSector::vshift().

                                                                                                  {
  closest = failIndex_;

  ETLDetId start_mod(detid);

  if (verticalShift == 0) {
    edm::LogWarning("MTDTopology") << "asking of a null verticalShift in ETL";
    return failIndex_;
  }
  int vsh = verticalShift > 0 ? 1 : -1;

  int sensor = start_mod.sensor();
  int module = start_mod.module();
  uint32_t modtyp = start_mod.modType();
  uint32_t discside = start_mod.discSide();
  int geomDetIndex;

  // distinguish numbering in prev8 / v8 geometries
  if (getMTDTopologyMode() == static_cast<int>(MTDTopologyMode::Mode::btlv3etlv8)) {
    geomDetIndex = 2 * (module - 1) + sensor;
  } else {
    geomDetIndex = module;
  }

  // ilayout number coincides at present with disc face, use this

  size_t iHome = (modtyp == etlVals_[discside].idDetType1_) ? 0 : 1;
  size_t iOther = (iHome == 0) ? 1 : 0;
  size_t iLeft = (etlVals_[discside].idDetType1_ == 1) ? 0 : 1;

  // for right type modules the offset of the total number of left modules needs to be added,
  // what matters here is the other type, i.e. if the starting module is left the vertical shift moves towards a right type, and viceversa

  size_t nmodOffset = (modtyp == 1) ? etlVals_[discside].start_copy_[iLeft].back() - 1 : 0;

  size_t iBin(etlVals_[discside].start_copy_[iHome].size());  // never allowed
  for (size_t iloop = 0; iloop < etlVals_[discside].start_copy_[iHome].size() - 1; iloop++) {
    if (geomDetIndex >= etlVals_[discside].start_copy_[iHome][iloop] &&
        geomDetIndex < etlVals_[discside].start_copy_[iHome][iloop + 1]) {
      iBin = iloop;
      break;
    }
  }

  if (iBin == etlVals_[discside].start_copy_[iHome].size()) {
    edm::LogWarning("MTDTopology") << "Module number not compatible with layout, abort";
    return failIndex_;
  }

  // define the interval of interest for the other type according to the vertical shift sign

  int iBinOther(iBin);
  if (iHome == 0 && vsh < 0) {
    iBinOther = iBin - 1;
  }
  if (iHome == 1 && vsh > 0) {
    iBinOther = iBin + 1;
  }
  if (iBinOther < 0 || iBinOther >= static_cast<int>(etlVals_[discside].start_copy_[iOther].size()) - 1) {
    return failIndex_;
  }

  // determine the position of the other type corresponding to the same column of the home type

  int vpos = etlVals_[discside].offset_[iHome][iBin] + geomDetIndex - etlVals_[discside].start_copy_[iHome][iBin] + 1;
  if (vpos <= etlVals_[discside].offset_[iOther][iBinOther]) {
    closest = etlVals_[discside].start_copy_[iOther][iBinOther];
  } else if (vpos > etlVals_[discside].offset_[iOther][iBinOther] +
                        etlVals_[discside].start_copy_[iOther][iBinOther + 1] -
                        etlVals_[discside].start_copy_[iOther][iBinOther] ||
             (vpos == etlVals_[discside].offset_[iOther][iBinOther] +
                          etlVals_[discside].start_copy_[iOther][iBinOther + 1] -
                          etlVals_[discside].start_copy_[iOther][iBinOther] &&
              iBinOther + 1 == static_cast<int>(etlVals_[discside].start_copy_[iOther].size()))) {
    closest = etlVals_[discside].start_copy_[iOther][iBinOther + 1] - 1;
  }
  if (closest < failIndex_) {
    closest = closest + nmodOffset - 1;
    return failIndex_;
  } else {
    // number of module shifted by 1 wrt the position in the array (i.e. module 1 has index 0)
    return etlVals_[discside].start_copy_[iOther][iBinOther] + vpos - 1 -
           etlVals_[discside].offset_[iOther][iBinOther] + nmodOffset - 1;
  }
}

Member Data Documentation

etlVals_

const ETLValues MTDTopology::etlVals_

private

Definition at line 42 of file MTDTopology.h.

Referenced by hshiftETL(), and vshiftETL().

failIndex_

constexpr size_t MTDTopology::failIndex_

staticprivate

Initial value:

=
      std::numeric_limits<unsigned int>::max()

Definition at line 44 of file MTDTopology.h.

Referenced by hshiftETL(), and vshiftETL().

mtdTopologyMode_

const int MTDTopology::mtdTopologyMode_

private

Definition at line 40 of file MTDTopology.h.

Referenced by getMTDTopologyMode().