#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)

	~HGCalNumberingScheme ()

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

	HGCalNumberingScheme ()=delete

Private Attributes
DetId::Detector	det_

const HGCalDDDConstants &	hgcons_

const HGCalGeometryMode::GeometryMode	mode_

std::string	name_

Detailed Description

Definition at line 16 of file HGCalNumberingScheme.h.

Constructor & Destructor Documentation

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

Definition at line 11 of file HGCalNumberingScheme.cc.

References det_, and name_.

     : hgcons_(hgc), mode_(hgc.geomMode()), det_(det), name_(name) {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCSim") << "Creating HGCalNumberingScheme for " << name_ << " Det " << det_;
 #endif
 }

HGCalNumberingScheme::~HGCalNumberingScheme ( )

Definition at line 20 of file HGCalNumberingScheme.cc.

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

HGCalNumberingScheme::HGCalNumberingScheme ( )

privatedelete

Member Function Documentation

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

private

Definition at line 90 of file HGCalNumberingScheme.cc.

References funct::abs(), det_, pfClusterIsolation_cfi::drMax, dzMax, DetId::HGCalHSc, DetId::HGCalHSi, hgcons_, HGCalDDDConstants::locateCell(), HGCalDDDConstants::locateCellTrap(), match(), convertSQLiteXML::ok, diffTwoXMLs::r1, diffTwoXMLs::r2, HGCalDDDConstants::rangeR(), HGCalDDDConstants::rangeZ(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, HGCalDDDConstants::waferFromPosition(), HGCalDDDConstants::waferZ(), geometryCSVtoXML::xx, and geometryCSVtoXML::xy.

Referenced by getUnitID().

                                                                                        {
   std::pair<float, float> xy;
   bool ok(false);
   double z1(0), drMax(10.0), dzMax(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(lay, id.waferU(), id.waferV(), id.cellU(), id.cellV(), false, true);
     z1 = hgcons_.waferZ(lay, false);
     ok = true;
     drMax = 10.0;
     dzMax = 1.0;
   } else if (DetId(index).det() == DetId::HGCalHSc) {
     HGCScintillatorDetId id = HGCScintillatorDetId(index);
     lay = id.layer();
     xy = hgcons_.locateCellTrap(lay, id.ietaAbs(), id.iphi(), false);
     z1 = hgcons_.waferZ(lay, false);
     ok = true;
     drMax = 50.0;
     dzMax = 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) < drMax) && (std::abs(z1 - z2) < dzMax);
     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 - 10.0) || (r1 > rrange.second +10.0) ||
                (z1 < zrange.first - 5.0) || (z1 > zrange.second + 5.0)) ?
               "***** ERROR *****" : "");
     if (!(match && inok && outok)) {
       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::HGCalHSi) {
         double wt = 0, xx = ((pos.z() > 0) ? pos.x() : -pos.x());
         int waferU, waferV, cellU, cellV, waferType;
         hgcons_.waferFromPosition(xx, pos.y(), lay, waferU, waferV, cellU, cellV, waferType, wt, true);
         xy = hgcons_.locateCell(lay, waferU, waferV, cellU, cellV, false, true, true);
         edm::LogVerbatim("HGCSim") 
       << "HGCalNumberingScheme " << HGCSiliconDetId(index) << " position " 
       << xy.first << ":" << xy.second;
       }
     }
   }
 }

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 26 of file HGCalNumberingScheme.cc.

References HGCalDDDConstants::assignCellTrap(), checkPosition(), TauDecayModes::dec, det_, HGCalGeometryMode::Hexagon8, HGCalGeometryMode::Hexagon8Full, hgcons_, mode_, DetId::rawId(), HGCalGeometryMode::Trapezoid, HGCalDDDConstants::waferFromPosition(), and geometryCSVtoXML::xx.

                                                                                                                       {
   // 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 << " Mode "
     << mode_ << ":" << HGCalGeometryMode::Hexagon8Full << ":"
     << HGCalGeometryMode::Hexagon8 << ":" << HGCalGeometryMode::Trapezoid;
 #endif
   if ((mode_ == HGCalGeometryMode::Hexagon8Full) || (mode_ == HGCalGeometryMode::Hexagon8)) {
     int cellU(0), cellV(0), waferType(-1), waferU(0), waferV(0);
     if (cell >= 0) {
       waferType = module / 1000000;
       waferU = module % 100;
       if ((module / 10000) % 10 > 0)
         waferU = -waferU;
       waferV = (module / 100) % 100;
       if ((module / 100000) % 10 > 0)
         waferV = -waferV;
       cellU = cell % 100;
       cellV = (cell / 100) % 100;
     } else if (mode_ == HGCalGeometryMode::Hexagon8Full) {
       double xx = (pos.z() > 0) ? pos.x() : -pos.x();
       hgcons_.waferFromPosition(xx, pos.y(), layer, waferU, waferV, cellU, cellV, waferType, wt);
     }
     if (waferType >= 0) {
       index = HGCSiliconDetId(det_, iz, waferType, layer, waferU, waferV, cellU, cellV).rawId();
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HGCSim")
     << "OK WaferType " << waferType << " Wafer " << waferU << ":" << waferV 
     << " Cell " << cellU << ":" << cellV;
     } else {
       edm::LogVerbatim("HGCSim") << "Bad WaferType " << waferType;
 #endif
     }
   } else if (mode_ == HGCalGeometryMode::Trapezoid) {
     std::array<int, 3> id = hgcons_.assignCellTrap(pos.x(), pos.y(), pos.z(), layer, false);
     if (id[2] >= 0) {
       index = HGCScintillatorDetId(id[2], layer, iz * id[0], id[1]).rawId();
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HGCSim") 
     << "Radius/Phi " << id[0] << ":" << id[1] << " Type " << id[2] 
     << " Layer|iz " << layer << ":" << iz << " " 
     << HGCScintillatorDetId(index);
     } else {
       edm::LogVerbatim("HGCSim") 
     << "Radius/Phi " << id[0] << ":" << id[1] << " Type " << id[2] 
     << " Layer|iz " << layer << ":" << iz << " ERROR";
 #endif
     }
   }
 #ifdef EDM_ML_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);
   }
 #endif
   return index;
 }