HFNoseNumberingScheme.cc

Go to the documentation of this file.

// File: HFNoseNumberingScheme.cc
// Description: Numbering scheme for HFNose detector
#include "SimG4CMS/Calo/interface/HFNoseNumberingScheme.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <iostream>

//#define EDM_ML_DEBUG

HFNoseNumberingScheme::HFNoseNumberingScheme(const HGCalDDDConstants& hgc) :
  hgcons_(hgc), mode_(hgc.geomMode()) {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCSim") << "Creating HFNoseNumberingScheme";
#endif
}

uint32_t HFNoseNumberingScheme::getUnitID(int layer, int module, int cell,
                      int iz, const G4ThreeVector &pos,
                      double& wt) {
  // module is the copy number of the wafer as placed in the layer
  uint32_t index(0);
  wt = 1.0;
  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   = HFNoseDetId(iz,waferType,layer,waferU,waferV,cellU,cellV).rawId();
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HFNSim") << "OK WaferType " << waferType << " Wafer "
                   << waferU << ":" << waferV << " Cell "
                   << cellU << ":" << cellV;
  } else {
    edm::LogVerbatim("HFNSim") << "Bad WaferType " << waferType;
#endif
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HFNSim") << "HFNoseNumberingScheme::i/p " << layer << ":"
                 << module << ":" << cell << ":" << iz << ":" 
                 << pos.x() << ":" << pos.y() << ":" << pos.z() 
                 << " ID " << std::hex << index << std::dec
                 << " wt " << wt;
#endif
  return index;
}