HGCalNumberingScheme Class Reference

#include <HGCalNumberingScheme.h>

Public Member Functions
uint32_t	getUnitID (int layer, int module, int cell, int iz, const G4ThreeVector &pos, double &wt)
	assigns the det id to a hit More...
	HGCalNumberingScheme (const HGCalDDDConstants &hgc, const DetId::Detector &det, const std::string &name, const std::string &fileName)
	HGCalNumberingScheme ()=delete
	~HGCalNumberingScheme ()

Private Member Functions
void	checkPosition (uint32_t index, const G4ThreeVector &pos, bool matchOnly, bool debug) const

Private Attributes
const DetId::Detector	det_
std::vector< int >	dumpCassette_
std::vector< int >	dumpDets_
int	firstLayer_
const HGCalDDDConstants &	hgcons_
std::vector< int >	indices_
const HGCalGeometryMode::GeometryMode	mode_
const std::string	name_

Detailed Description

Definition at line 19 of file HGCalNumberingScheme.h.

Constructor & Destructor Documentation

HGCalNumberingScheme() [1/2]

HGCalNumberingScheme::HGCalNumberingScheme	(	const HGCalDDDConstants &	hgc,
		const DetId::Detector &	det,
		const std::string &	name,
		const std::string &	fileName
	)

Definition at line 21 of file HGCalNumberingScheme.cc.

References edmScanValgrind::buffer, HGCalCassette::cassetteIndex(), det_, dumpCassette_, dumpDets_, MillePedeFileConverter_cfg::fileName, firstLayer_, edm::FileInPath::fullPath(), HGCalDDDConstants::getLayerOffset(), HGCalGeometryMode::Hexagon8, HGCalGeometryMode::Hexagon8Cassette, HGCalGeometryMode::Hexagon8File, HGCalGeometryMode::Hexagon8Full, HGCalGeometryMode::Hexagon8Module, hgcons_, indices_, hcalRecHitTable_cff::iphi, mps_monitormerge::items, nano_mu_digi_cff::layer, mode_, name_, relativeConstraints::ring, CaloSimUtils::splitString(), AlCaHLTBitMon_QueryRunRegistry::string, HGCalTileIndex::tileIndex(), HGCalDDDConstants::tileTrapezoid(), HGCalGeometryMode::Trapezoid, HGCalGeometryMode::TrapezoidCassette, HGCalGeometryMode::TrapezoidFile, HGCalGeometryMode::TrapezoidModule, HGCalDDDConstants::waferHexagon8File(), HGCalWaferIndex::waferIndex(), HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), and ecaldqm::zside().

    : hgcons_(hgc), mode_(hgc.geomMode()), det_(det), name_(name) {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCSim") << "Creating HGCalNumberingScheme for " << name_ << " Det " << det_ << " Mode " << mode_
                             << ":" << HGCalGeometryMode::Hexagon8Full << ":" << HGCalGeometryMode::Hexagon8 << ":"
                             << HGCalGeometryMode::Hexagon8File << ":" << HGCalGeometryMode::Hexagon8Module << ":"
                             << ":" << HGCalGeometryMode::Hexagon8Cassette << ":" << HGCalGeometryMode::Trapezoid << ":"
                             << HGCalGeometryMode::TrapezoidFile << ":" << HGCalGeometryMode::TrapezoidModule << ":"
                             << HGCalGeometryMode::TrapezoidCassette;
#endif
  firstLayer_ = hgcons_.getLayerOffset();
  if (!fileName.empty()) {
    edm::FileInPath filetmp1("SimG4CMS/Calo/data/" + fileName);
    std::string filetmp2 = filetmp1.fullPath();
    std::ifstream fInput(filetmp2.c_str());
    if (!fInput.good()) {
      edm::LogVerbatim("HGCalSim") << "Cannot open file " << filetmp2;
    } else {
      char buffer[80];
      while (fInput.getline(buffer, 80)) {
        std::vector<std::string> items = CaloSimUtils::splitString(std::string(buffer));
        if (items.size() == 3) {
          if (hgcons_.waferHexagon8File()) {
            int layer = std::atoi(items[0].c_str());
            int waferU = std::atoi(items[1].c_str());
            int waferV = std::atoi(items[2].c_str());
            indices_.emplace_back(HGCalWaferIndex::waferIndex(layer, waferU, waferV, false));
          } else if (hgcons_.tileTrapezoid()) {
            int layer = std::atoi(items[0].c_str());
            int ring = std::atoi(items[1].c_str());
            int iphi = std::atoi(items[2].c_str());
            indices_.emplace_back(HGCalTileIndex::tileIndex(layer, ring, iphi));
          }
        } else if (items.size() == 1) {
          int dumpdet = std::atoi(items[0].c_str());
          dumpDets_.emplace_back(dumpdet);
        } else if (items.size() == 4) {
          int idet = std::atoi(items[0].c_str());
          int layer = std::atoi(items[1].c_str());
          int zside = std::atoi(items[2].c_str());
          int cassette = std::atoi(items[3].c_str());
          dumpCassette_.emplace_back(HGCalCassette::cassetteIndex(idet, layer, zside, cassette));
        }
      }
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HGCalSim") << "Reads in " << indices_.size() << ":" << dumpDets_.size() << ":"
                                   << dumpCassette_.size() << " component information from " << filetmp2
                                   << " Layer Offset " << firstLayer_;
#endif
      fInput.close();
    }
  }
}

HGCalNumberingScheme() [2/2]

HGCalNumberingScheme::HGCalNumberingScheme ( )

delete

~HGCalNumberingScheme()

HGCalNumberingScheme::~HGCalNumberingScheme

(

)

Definition at line 78 of file HGCalNumberingScheme.cc.

                                            {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCSim") << "Deleting HGCalNumberingScheme";
#endif
}

Member Function Documentation

checkPosition()

void HGCalNumberingScheme::checkPosition ( uint32_t  index, const G4ThreeVector &  pos, bool  matchOnly, bool  debug  ) const

private

Definition at line 198 of file HGCalNumberingScheme.cc.

References funct::abs(), debug, det_, HGC3DClusterGenMatchSelector_cfi::dR, PVValHelper::dx, PVValHelper::dy, DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, hgcons_, hcalRecHitTable_cff::iphi, HGCalDDDConstants::locateCell(), HGCalDDDConstants::locateCellTrap(), match(), convertSQLiteXML::ok, diffTwoXMLs::r1, diffTwoXMLs::r2, HGCalDDDConstants::rangeR(), HGCalDDDConstants::rangeZ(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, HGCalDDDConstants::waferFromPosition(), HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), HGCalDDDConstants::waferZ(), geometryCSVtoXML::xx, geometryCSVtoXML::xy, testProducerWithPsetDescEmpty_cfi::z2, TkClusParameters_cff::zrange, and ecaldqm::zside().

Referenced by getUnitID().

                                                                                                                   {
  std::pair<float, float> xy;
  bool ok(false);
  double z1(0), tolR(14.0), tolZ(1.0);
  int lay(-1);
  if (index == 0) {
  } else if (DetId(index).det() == DetId::HGCalHSi) {
    HGCSiliconDetId id = HGCSiliconDetId(index);
    lay = id.layer();
    xy = hgcons_.locateCell(
        id.zside(), lay, id.waferU(), id.waferV(), id.cellU(), id.cellV(), false, true, false, false);
    z1 = hgcons_.waferZ(lay, false);
    ok = true;
    tolR = 14.0;
    tolZ = 1.0;
  } else if (DetId(index).det() == DetId::HGCalHSc) {
    HGCScintillatorDetId id = HGCScintillatorDetId(index);
    lay = id.layer();
    xy = hgcons_.locateCellTrap(id.zside(), lay, id.ietaAbs(), id.iphi(), false, false);
    z1 = hgcons_.waferZ(lay, false);
    ok = true;
    tolR = 50.0;
    tolZ = 5.0;
  }
  if (ok) {
    double r1 = std::sqrt(xy.first * xy.first + xy.second * xy.second);
    double r2 = pos.perp();
    double z2 = std::abs(pos.z());
    std::pair<double, double> zrange = hgcons_.rangeZ(false);
    std::pair<double, double> rrange = hgcons_.rangeR(z2, false);
    bool match = (std::abs(r1 - r2) < tolR) && (std::abs(z1 - z2) < tolZ);
    bool inok = ((r2 >= rrange.first) && (r2 <= rrange.second) && (z2 >= zrange.first) && (z2 <= zrange.second));
    bool outok = ((r1 >= rrange.first) && (r1 <= rrange.second) && (z1 >= zrange.first) && (z1 <= zrange.second));
    std::string ck = (((r1 < rrange.first - tolR) || (r1 > rrange.second + tolR) || (z1 < zrange.first - tolZ) ||
                       (z1 > zrange.second + tolZ))
                          ? "***** ERROR *****"
                          : "");
    if (matchOnly && match)
      ck = "";
    if (!(match && inok && outok) || debug) {
      edm::LogVerbatim("HGCSim") << "HGCalNumberingScheme::Detector " << det_ << " Layer " << lay << " R " << r2 << ":"
                                 << r1 << ":" << rrange.first << ":" << rrange.second << " Z " << z2 << ":" << z1 << ":"
                                 << zrange.first << ":" << zrange.second << " Match " << match << ":" << inok << ":"
                                 << outok << " " << ck;
      edm::LogVerbatim("HGCSim") << "Original " << pos.x() << ":" << pos.y() << " return " << xy.first << ":"
                                 << xy.second;
      if ((DetId(index).det() == DetId::HGCalEE) || (DetId(index).det() == DetId::HGCalHSi)) {
        int zside = (pos.z() > 0) ? 1 : -1;
        double wt(0), xx(zside * pos.x());
        int waferU, waferV, cellU, cellV, waferType;
        hgcons_.waferFromPosition(xx, pos.y(), zside, lay, waferU, waferV, cellU, cellV, waferType, wt, false, true);
        xy = hgcons_.locateCell(zside, lay, waferU, waferV, cellU, cellV, false, true, false, true);
        double dx = (xx - xy.first);
        double dy = (pos.y() - xy.second);
        double dR = std::sqrt(dx * dx + dy * dy);
        ck = (dR > tolR) ? " ***** ERROR *****" : "";
        edm::LogVerbatim("HGCSim") << "HGCalNumberingScheme " << HGCSiliconDetId(index) << " original position " << xx
                                   << ":" << pos.y() << " derived " << xy.first << ":" << xy.second << " Difference "
                                   << dR << ck;
      }
    }
  }
}

getUnitID()

uint32_t HGCalNumberingScheme::getUnitID	(	int	layer,
		int	module,
		int	cell,
		int	iz,
		const G4ThreeVector &	pos,
		double &	wt
	)

assigns the det id to a hit

Definition at line 84 of file HGCalNumberingScheme.cc.

References HGCalDDDConstants::assignCellTrap(), HGCalCassette::cassetteIndex(), checkPosition(), debug, TauDecayModes::dec, det_, hcalRecHitTable_cff::detId, dumpCassette_, dumpDets_, spr::find(), firstLayer_, HGCalDDDConstants::getParameter(), HGCalTypes::getUnpackedCellU(), HGCalTypes::getUnpackedCellV(), HGCalTypes::getUnpackedType(), HGCalTypes::getUnpackedU(), HGCalTypes::getUnpackedV(), HGCalGeometryMode::Hexagon8, HGCalGeometryMode::Hexagon8Cassette, HGCalGeometryMode::Hexagon8Module, HGCSiliconDetId::HGCalCoarseThick, hgcons_, indices_, nano_mu_digi_cff::layer, mode_, name_, HGCalParameters::waferInfo::part, DetId::rawId(), HGCalTileIndex::tileIndex(), HGCalDDDConstants::tileSiPM(), HGCalDDDConstants::tileTrapezoid(), HGCalDDDConstants::tileType(), HGCalDDDConstants::waferFromPosition(), HGCalTypes::WaferFull, HGCalDDDConstants::waferHexagon8(), HGCalDDDConstants::waferHexagon8File(), HGCalWaferIndex::waferIndex(), HGCalDDDConstants::waferInfo(), HGCalParameters::waferInfoMap_, HGCalDDDConstants::waferType(), HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), geometryCSVtoXML::xx, and ecaldqm::zside().

                                                                                                                      {
  // module is the copy number of the wafer as placed in the layer
  uint32_t index(0);
  wt = 1.0;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCSim") << "HGCalNumberingScheme:: input Layer " << layer << " Module " << module << " Cell "
                             << cell << " iz " << iz << " Position " << pos;
#endif
  if (hgcons_.waferHexagon8()) {
    int cellU(0), cellV(0), waferType(-1), waferU(0), waferV(0);
    bool debug(false);
    if (cell >= 0) {
      waferType = HGCalTypes::getUnpackedType(module);
      waferU = HGCalTypes::getUnpackedU(module);
      waferV = HGCalTypes::getUnpackedV(module);
      cellU = HGCalTypes::getUnpackedCellU(cell);
      cellV = HGCalTypes::getUnpackedCellV(cell);
    } else if (mode_ != HGCalGeometryMode::Hexagon8) {
      int zside = (pos.z() > 0) ? 1 : -1;
      double xx = zside * pos.x();
      int wU = HGCalTypes::getUnpackedU(module);
      int wV = HGCalTypes::getUnpackedV(module);
      if (!indices_.empty()) {
        int indx = HGCalWaferIndex::waferIndex(firstLayer_ + layer, wU, wV, false);
        if (std::find(indices_.begin(), indices_.end(), indx) != indices_.end())
          debug = true;
      }
      if (!dumpDets_.empty()) {
        if ((std::find(dumpDets_.begin(), dumpDets_.end(), det_) != dumpDets_.end()) &&
            (hgcons_.waferInfo(layer, wU, wV).part != HGCalTypes::WaferFull))
          debug = true;
      }
      if (!dumpCassette_.empty()) {
        int indx = HGCalWaferIndex::waferIndex(firstLayer_ + layer, wU, wV, false);
        auto ktr = hgcons_.getParameter()->waferInfoMap_.find(indx);
        if (ktr != hgcons_.getParameter()->waferInfoMap_.end()) {
          if (std::find(dumpCassette_.begin(),
                        dumpCassette_.end(),
                        HGCalCassette::cassetteIndex(det_, firstLayer_ + layer, zside, (ktr->second).cassette)) !=
              dumpCassette_.end())
            debug = true;
        }
      }
      hgcons_.waferFromPosition(xx, pos.y(), zside, layer, waferU, waferV, cellU, cellV, waferType, wt, false, debug);
    }
    if (waferType >= 0) {
      if (hgcons_.waferHexagon8File()) {
        int type = hgcons_.waferType(layer, waferU, waferV, true);
        if (type != waferType) {
#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("HGCSim") << "HGCalNumberingScheme:: " << name_ << " Layer|u|v|Index|module|cell " << layer
                                     << ":" << waferU << ":" << waferV << ":"
                                     << HGCalWaferIndex::waferIndex(layer, waferU, waferV, false) << ":" << module
                                     << ":" << cell << " has a type mismatch " << waferType << ":" << type;
#endif
          if (type != HGCSiliconDetId::HGCalCoarseThick)
            waferType = type;
        }
      }
      index = HGCSiliconDetId(det_, iz, waferType, layer, waferU, waferV, cellU, cellV).rawId();
      if (debug) {
        int zside = (pos.z() > 0) ? 1 : -1;
        double xx = zside * pos.x();
        edm::LogVerbatim("HGCSim") << "OK Input " << det_ << ":" << zside << ":" << layer << ":" << xx << ":" << pos.y()
                                   << ":" << module << ":" << cell << " WaferType " << waferType << " Wafer " << waferU
                                   << ":" << waferV << " Cell " << cellU << ":" << cellV << " wt " << wt << " Mode "
                                   << mode_;
      }
#ifdef EDM_ML_DEBUG
    } else {
      edm::LogVerbatim("HGCSim") << "Bad WaferType " << waferType << " for Layer:u:v " << layer << ":" << waferU << ":"
                                 << waferV;
#endif
    }
  } else if (hgcons_.tileTrapezoid()) {
    std::array<int, 3> id = hgcons_.assignCellTrap(pos.x(), pos.y(), pos.z(), layer, false);
    if (id[2] >= 0) {
      std::pair<int, int> typm = hgcons_.tileType(layer, id[0], 0);
      HGCScintillatorDetId detId(id[2], layer, iz * id[0], id[1], false, 0);
      if (typm.first >= 0) {
        detId.setType(typm.first);
        detId.setSiPM(typm.second);
      }
      index = detId.rawId();
      bool debug(false);
      if (!indices_.empty()) {
        int indx = HGCalTileIndex::tileIndex(layer, id[0], id[1]);
        if (std::find(indices_.begin(), indices_.end(), indx) != indices_.end())
          debug = true;
      }
      if (debug)
        edm::LogVerbatim("HGCSim") << "Radius/Phi " << id[0] << ":" << id[1] << " Type " << id[2] << ":" << typm.first
                                   << " SiPM " << typm.second << ":" << hgcons_.tileSiPM(typm.second) << " Layer "
                                   << layer << " z " << iz << " " << detId << " wt " << wt << " position " << pos
                                   << " R " << pos.perp();
#ifdef EDM_ML_DEBUG
    } else {
      edm::LogVerbatim("HGCSim") << "Radius/Phi " << id[0] << ":" << id[1] << " Type " << id[2] << " Layer|iz " << layer
                                 << ":" << iz << " ERROR";
#endif
    }
  }
#ifdef EDM_ML_DEBUG
  bool matchOnly = ((mode_ == HGCalGeometryMode::Hexagon8Module) || (mode_ == HGCalGeometryMode::Hexagon8Cassette));
  bool debug = hgcons_.waferHexagon8File();
  if (debug)
    edm::LogVerbatim("HGCSim") << "HGCalNumberingScheme::i/p " << det_ << ":" << layer << ":" << module << ":" << cell
                               << ":" << iz << ":" << pos.x() << ":" << pos.y() << ":" << pos.z() << " ID " << std::hex
                               << index << std::dec << " wt " << wt;
  checkPosition(index, pos, matchOnly, debug);
#endif
  return index;
}

Member Data Documentation

det_

const DetId::Detector HGCalNumberingScheme::det_

private

Definition at line 38 of file HGCalNumberingScheme.h.

Referenced by checkPosition(), getUnitID(), and HGCalNumberingScheme().

dumpCassette_

std::vector<int> HGCalNumberingScheme::dumpCassette_

private

Definition at line 43 of file HGCalNumberingScheme.h.

Referenced by getUnitID(), and HGCalNumberingScheme().

dumpDets_

std::vector<int> HGCalNumberingScheme::dumpDets_

private

Definition at line 42 of file HGCalNumberingScheme.h.

Referenced by getUnitID(), and HGCalNumberingScheme().

firstLayer_

int HGCalNumberingScheme::firstLayer_

private

Definition at line 40 of file HGCalNumberingScheme.h.

Referenced by getUnitID(), and HGCalNumberingScheme().

hgcons_

const HGCalDDDConstants& HGCalNumberingScheme::hgcons_

private

Definition at line 36 of file HGCalNumberingScheme.h.

Referenced by checkPosition(), getUnitID(), and HGCalNumberingScheme().

indices_

std::vector<int> HGCalNumberingScheme::indices_

private

Definition at line 41 of file HGCalNumberingScheme.h.

Referenced by getUnitID(), and HGCalNumberingScheme().

mode_

const HGCalGeometryMode::GeometryMode HGCalNumberingScheme::mode_

private

Definition at line 37 of file HGCalNumberingScheme.h.

Referenced by getUnitID(), and HGCalNumberingScheme().

name_

const std::string HGCalNumberingScheme::name_

private

Definition at line 39 of file HGCalNumberingScheme.h.

Referenced by getUnitID(), and HGCalNumberingScheme().