EcalScDetId.cc

Go to the documentation of this file.

#include "DataFormats/EcalDetId/interface/EcalScDetId.h"
#include "FWCore/Utilities/interface/Exception.h"

#include <ostream>
#include <cassert>
#include <mutex>

short EcalScDetId::xyz2HashedIndex[EcalScDetId::IX_MAX][EcalScDetId::IY_MAX][EcalScDetId::nEndcaps];

EcalScDetId EcalScDetId::hashedIndex2DetId[kSizeForDenseIndexing];


EcalScDetId::EcalScDetId() : DetId() {
}

EcalScDetId::EcalScDetId(uint32_t rawid) : DetId(rawid) {
}

EcalScDetId::EcalScDetId(int ix, int iy, int iz) : DetId(Ecal,EcalEndcap) 
{
  if(!validDetId(ix,iy,iz))
    {
      throw cms::Exception("InvalidDetId") << "EcalScDetId:  Cannot create object.  Indexes out of bounds \n" 
                                           << "x = " << ix << " y = " << iy << " z = " << iz;
    }
  const int scBit = 1<<15; //bit set to 1 to distinguish from crystal id (EEDetId)
  //                         and for a reasonale behaviour of DetId ccomparison operators.
  id_|=(iy&0x7f)|((ix&0x7f)<<7)|((iz>0)?(1<<14):(0))|scBit;
}
  
EcalScDetId::EcalScDetId(const DetId& gen) {
  if (!gen.null() && (gen.det()!=Ecal || gen.subdetId()!=EcalEndcap)) {
    throw cms::Exception("InvalidDetId"); 
  }
  id_=gen.rawId();
}
  
EcalScDetId& EcalScDetId::operator=(const DetId& gen) {
  if (!gen.null() && ( gen.det()!=Ecal || gen.subdetId()!=EcalEndcap )) {
    throw cms::Exception("InvalidDetId"); 
  }
  id_=gen.rawId();
  return *this;
}
  
int EcalScDetId::iquadrant() const {
  const int xMiddle = IX_MAX/2; //y = 0 between xMiddle and xMiddle+1
  const int yMiddle = IY_MAX/2; //x = 0 between yMiddle and yMiddle+1
  if (iy()>yMiddle){// y>0
    if(ix()>xMiddle)   //             A y        
      return 1;        //             |          
    else           //      Q2     |    Q1        
      return 2;        //             |          
  } else{// y<0        //   ----------o---------> x   
    if(ix()>xMiddle)   //             |           
      return 4;        //      Q3     |    Q4       
    else           //             |               
      return 3;
  }
  //Should never be reached
  return -1;
}  

bool EcalScDetId::validDetId(int iX, int iY, int iZ) {
  static const char endcapMap[401] = {
    "       XXXXXX       "
    "    XXXXXXXXXXXX    "
    "   XXXXXXXXXXXXXX   "
    "  XXXXXXXXXXXXXXXX  "
    " XXXXXXXXXXXXXXXXXX "
    " XXXXXXXXXXXXXXXXXX "             //    Z
    " XXXXXXXXXXXXXXXXXX "             //     x-----> X
    "XXXXXXXXXXXXXXXXXXXX"             //     |
    "XXXXXXXXX  XXXXXXXXX"             //     |
    "XXXXXXXX    XXXXXXXX"//_          //     |
    "XXXXXXXX    XXXXXXXX"             //     V Y
    "XXXXXXXXX  XXXXXXXXX"
    "XXXXXXXXXXXXXXXXXXXX"
    " XXXXXXXXXXXXXXXXXX "
    " XXXXXXXXXXXXXXXXXX "
    " XXXXXXXXXXXXXXXXXX "
    "  XXXXXXXXXXXXXXXX  "
    "   XXXXXXXXXXXXXX   "
    "    XXXXXXXXXXXX    "
    "       XXXXXX       "};

  return abs(iZ)==1 && endcapMap[iX-1+(iY-1)*20]!=' ';
}

std::ostream& operator<<(std::ostream& s,const EcalScDetId& id) {
  return s << "(EE iz " << ((id.zside()>0)?("+ "):("- ")) << " ix " << id.ix() << " , iy " << id.iy() << ')';
}

//NOTE: When looping is possible in constexpr, this should be changed to one
static std::once_flag initializedFlag;
void EcalScDetId::checkHashedIndexMap(){
  std::call_once(initializedFlag, []() 
  {
    int hashedIndex = -1;
    for(int iZ = -1; iZ <= +1; iZ+=2){
      for(int iY = IY_MIN; iY <= IY_MAX; ++iY){
    for(int iX = IX_MIN; iX <= IX_MAX; ++iX){
      if(validDetId(iX,iY,iZ)){
        xyz2HashedIndex[iX-IX_MIN][iY-IY_MIN][iZ>0?1:0] = ++hashedIndex;
        assert((unsigned)hashedIndex < sizeof(hashedIndex2DetId)/sizeof(hashedIndex2DetId[0]));
        hashedIndex2DetId[hashedIndex] = EcalScDetId(iX, iY, iZ);
      }
    }
      }
    }
  });
}