HcalZDCDetId.cc

Go to the documentation of this file.

#include "DataFormats/HcalDetId/interface/HcalZDCDetId.h"
#include "FWCore/Utilities/interface/Exception.h"

const int HcalZDCDetId::kZDCChannelMask;
const int HcalZDCDetId::kZDCSectionMask;
const int HcalZDCDetId::kZDCSectionOffset;
const int HcalZDCDetId::kZDCZsideMask;
const int HcalZDCDetId::kZDCRPDMask;
const int HcalZDCDetId::SubdetectorId;

HcalZDCDetId::HcalZDCDetId() : DetId() {
}


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

HcalZDCDetId::HcalZDCDetId(Section section, bool true_for_positive_eta, int channel) : DetId(DetId::Calo,SubdetectorId) {
  if (section == RPD) {
    id_ |= (Unknown&kZDCSectionMask)<<kZDCSectionOffset;
    id_ |= kZDCRPDMask;
    id_ |= ((channel-1)&kZDCChannelMask);
  } else {
    id_ |= (section&kZDCSectionMask)<<kZDCSectionOffset;
    id_ |= (channel&kZDCChannelMask);
  }
  if (true_for_positive_eta) id_ |= kZDCZsideMask;
}

HcalZDCDetId::HcalZDCDetId(const DetId& gen) {
  if (!gen.null() && (gen.det()!=Calo || gen.subdetId()!=SubdetectorId)) {
    throw cms::Exception("Invalid DetId") 
      << "Cannot initialize ZDCDetId from " << std::hex << gen.rawId() << std::dec; 
  }
  id_ = gen.rawId();
}

HcalZDCDetId& HcalZDCDetId::operator=(const DetId& gen) {
  if (!gen.null() && (gen.det()!=Calo || gen.subdetId()!=SubdetectorId)) {
    throw cms::Exception("Invalid DetId")
      << "Cannot assign ZDCDetId from " << std::hex << gen.rawId() << std::dec; 
  }
  id_=gen.rawId();
  return *this;
}

HcalZDCDetId::Section HcalZDCDetId::section() const {
  if (id_ & kZDCRPDMask) 
    return RPD;
  else
    return (Section)((id_>>kZDCSectionOffset)&kZDCSectionMask);
}

int HcalZDCDetId::depth() const {
  const int se ( section() ) ;
  if      (se == EM)  return 1;
  else if (se == HAD) return (channel()+2);
  else if (se == RPD) return 2;
  else                return channel();
}

int HcalZDCDetId::channel() const {
  const int se ( section() ) ;
  if (se == RPD)  return (1+(id_&kZDCChannelMask));
  else            return (id_&kZDCChannelMask);
}

uint32_t HcalZDCDetId::denseIndex() const {
  const int se ( section() ) ;
  uint32_t  di = (channel() - 1 + 
          ( se == RPD ? 2*kDepRun1 + ( zside()<0 ? 0 : kDepRPD ) :
            ( ( zside()<0 ? 0 : kDepRun1 ) + 
              ( se == HAD  ? kDepEM :
            ( se == LUM ? kDepEM + kDepHAD : 0 ) ) ) ) );
  return di;
}

HcalZDCDetId HcalZDCDetId::detIdFromDenseIndex( uint32_t di ) {
  if ( validDenseIndex( di ) ) {
    bool     lz (false);
    uint32_t dp(0);
    Section  se(Unknown);
    if (di >= 2*kDepRun1) {
      lz          = (di >= (kDepRun1+kDepTot));
      se          = RPD;
      dp          = 1 + ((di-2*kDepRun1) % kDepRPD);
    } else {
      lz          = (di >= kDepRun1);
      uint32_t in = (di%kDepRun1) ;
      se          = (in<kDepEM ? EM :
             (in<kDepEM+kDepHAD ? HAD : LUM) );
      dp          = (EM == se ? in+1 :
             (HAD == se ? in-kDepEM+1 : in-kDepEM-kDepHAD+1) );
    }
    return HcalZDCDetId( se, lz, dp ) ;
  } else{
    return HcalZDCDetId() ;
  }
}

bool HcalZDCDetId::validDetId( Section se , int  dp ) {
  bool flag = (dp >= 1 && (
               ( ( se == EM      ) &&
                 ( dp <= kDepEM  )    ) ||
               ( ( se == HAD     ) &&
                 ( dp <= kDepHAD )    ) ||
               ( ( se == LUM     ) &&
                 ( dp <= kDepLUM )    ) ||
               ( ( se == RPD     ) &&
                 ( dp <= kDepRPD )    ) 
               )
           );
  return flag;
}

std::ostream& operator<<(std::ostream& s,const HcalZDCDetId& id) {
  s << "(ZDC" << ((id.zside()==1)?("+"):("-"));
  switch (id.section()) {
  case(HcalZDCDetId::EM)  : s << " EM ";  break;
  case(HcalZDCDetId::HAD) : s << " HAD "; break;
  case(HcalZDCDetId::LUM) : s << " LUM "; break;
  case(HcalZDCDetId::RPD) : s << " RPD "; break;
  default : s << " UNKNOWN ";
  } 
  return s << id.channel() << "," << id.depth() << ')';
}