EcalEndcapNumberingScheme.cc

Go to the documentation of this file.

// File: EcalEndcapNumberingScheme.cc
// Description: Numbering scheme for endcap electromagnetic calorimeter
#include "Geometry/EcalCommonData/interface/EcalEndcapNumberingScheme.h"
#include "DataFormats/EcalDetId/interface/EEDetId.h"
 
#include <iostream>
#include <iomanip>
 
EcalEndcapNumberingScheme::EcalEndcapNumberingScheme() : EcalNumberingScheme() {
  edm::LogInfo("EcalGeom") << "Creating EcalEndcapNumberingScheme";
}
 
EcalEndcapNumberingScheme::~EcalEndcapNumberingScheme() {
  edm::LogInfo("EcalGeom") << "Deleting EcalEndcapNumberingScheme";
}
uint32_t EcalEndcapNumberingScheme::getUnitID(const EcalBaseNumber& baseNumber) const {
  const uint32_t nLevels(baseNumber.getLevels());
 
  //LogDebug("EcalGeom") << "ECalEndcapNumberingScheme geometry levels = " << nLevels;
 
  if (7 > nLevels) {
    edm::LogWarning("EcalGeom") << "ECalEndcapNumberingScheme::getUnitID(): "
                                << "Not enough levels found in EcalBaseNumber ( " << nLevels << ") Returning 0";
    return 0;
  }
 
  if (nLevels <= 10) {
    // Static geometry
 
    if (baseNumber.getLevels() < 1) {
      edm::LogWarning("EcalGeom") << "EalEndcaplNumberingScheme::getUnitID: No "
                                  << "level found in EcalBaseNumber Returning 0";
      return 0;
    }
 
    int PVid = baseNumber.getCopyNumber(0);
    int MVid = 1;
    if (baseNumber.getLevels() > 1)
      MVid = baseNumber.getCopyNumber(1);
    else
      edm::LogWarning("EcalGeom") << "ECalEndcapNumberingScheme::getUnitID: Null"
                                  << " pointer to alveole ! Use default id=1";
 
    int zside = baseNumber.getCopyNumber("EREG");
    zside = 2 * (1 - zside) + 1;
    int module_number = MVid;
    int crystal_number = PVid;
 
    uint32_t intindex = EEDetId(module_number, crystal_number, zside, EEDetId::SCCRYSTALMODE).rawId();
 
    //LogDebug("EcalGeom") << "EcalEndcapNumberingScheme: zside = "  << zside
    //                     << " super crystal = " << module_number << " crystal = "
    //                     << crystal_number << " packed index = 0x" << std::hex
    //                     << intindex << std::dec;
 
    return intindex;
  } else {
    // algorithmic geometry
 
    const int ic(baseNumber.getCopyNumber(0) % 100);  // crystal #, 0-44
    const int icx(ic / 10);
    const int icy(ic % 5);
    const int is(baseNumber.getCopyNumber(2) % 100);  // supercrystal #, 0-99
    const int isx(is / 10);
    const int isy(is % 10);
 
    const int iq(3 - 2 * baseNumber.getCopyNumber(3));  // quadrant #, -1, +1
    const int id(3 - 2 * baseNumber.getCopyNumber(5));  // dee      #, -1, +1
 
    const int iz(3 - 2 * baseNumber.getCopyNumber(7));  // z: -1, +1
 
    const int ix(50 + id * iz * (5 * isx + icx + 1) - (id * iz - 1) / 2);  // x: 1-100
    const int iy(50 + id * iq * (5 * isy + icy + 1) - (id * iq - 1) / 2);  // y: 1-100
 
    const uint32_t idet(DetId(DetId::Ecal, EEDetId::Subdet) | (((0 < iz ? 0x4000 : 0)) + (ix << 7) + iy));
 
    //*************************** ERROR CHECKING **********************************
 
    if (0 > icx || 4 < icx || 0 > icy || 4 < icy) {
      edm::LogWarning("EdalGeom") << "ECalEndcapNumberingScheme::getUnitID(): "
                                  << "****************** Bad crystal number = " << ic
                                  << ", Volume Name = " << baseNumber.getLevelName(0);
      return 0;
    }
 
    if (0 > isx || 9 < isx || 0 > isy || 9 < isy) {
      edm::LogWarning("EdalGeom") << "ECalEndcapNumberingScheme::getUnitID(): "
                                  << "****************** Bad supercrystal number = " << is
                                  << ", Volume Name = " << baseNumber.getLevelName(3);
      return 0;
    }
 
    if (1 != iq && -1 != iq) {
      edm::LogWarning("EdalGeom") << "ECalEndcapNumberingScheme::getUnitID(): "
                                  << "****************** Bad quadrant number = " << iq
                                  << ", Volume Name = " << baseNumber.getLevelName(4);
      return 0;
    }
 
    if (1 != id && -1 != id) {
      edm::LogWarning("EdalGeom") << "ECalEndcapNumberingScheme::getUnitID(): "
                                  << "****************** Bad dee number = " << id
                                  << ", Volume Name = " << baseNumber.getLevelName(6);
      return 0;
    }
 
    if (-1 != iz && 1 != iz) {
      edm::LogWarning("EdalGeom") << "ECalEndcapNumberingScheme::getUnitID(): "
                                  << "****************** Bad z-end number = " << iz
                                  << ", Volume Name = " << baseNumber.getLevelName(8);
      return 0;
    }
 
    if (!EEDetId::validDetId(ix, iy, iz)) {
      edm::LogWarning("EdalGeom") << "ECalEndcapNumberingScheme::getUnitID(): "
                                  << "****************** Bad DetId numbers = " << ix << ", " << iy << ", " << iz;
      return 0;
    }
 
    //*************************************************************************************
    /*
      edm::LogWarning("EdalGeom") << "ECalEndcapNumberingScheme::getUnitID(): "
                  <<std::dec<< ix << ", " << iy << ", " <<iq << ", " << id << ", " << iz << ", " << std::hex << idet              ;
 
      edm::LogWarning("EdalGeom") << "ECalEndcapNumberingScheme::EEDetId: "
                  << EEDetId(idet)              ;
*/
    return idet;
  }
}