HGCalDetId.cc

Go to the documentation of this file.

#include "DataFormats/ForwardDetId/interface/HGCalDetId.h"
#include <iostream>

//#define DebugLog

HGCalDetId::HGCalDetId() : DetId() {
}

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

HGCalDetId::HGCalDetId(ForwardSubdetector subdet, int zp, int lay, int wafertype, int wafer, int cell) : DetId(Forward,subdet) {  

  if (wafertype < 0) wafertype = 0;
  if (cell>kHGCalCellMask || cell<0 || wafer>kHGCalWaferMask || wafer<0 || wafertype>kHGCalWaferTypeMask || lay>kHGCalLayerMask || lay<0) {
#ifdef DebugLog
    std::cout << "[HGCalDetId] request for new id for"
          << " layer=" << lay << ":" << kHGCalLayerMask
          << " @ zp=" << zp 
          << " wafer=" << wafer << ":" << kHGCalWaferMask
          << " waferType=" << wafertype << ":" << kHGCalWaferTypeMask
          << " cell=" << cell << ":" << kHGCalCellMask
          << " for subdet=" << subdet 
          << " has one or more fields out of bounds and will be reset" 
          << std::endl;
#endif
    zp = lay = wafertype = wafer = cell = 0;
  }
  id_ |= ((cell   & kHGCalCellMask)         << kHGCalCellOffset);
  id_ |= ((wafer  & kHGCalWaferMask)        << kHGCalWaferOffset);
  id_ |= ((wafertype & kHGCalWaferTypeMask) << kHGCalWaferTypeOffset);
  id_ |= ((lay    & kHGCalLayerMask)        << kHGCalLayerOffset);
  if (zp>0) id_ |= ((zp & kHGCalZsideMask)  << kHGCalZsideOffset);
}

HGCalDetId::HGCalDetId(const DetId& gen) {
  id_ = gen.rawId();
}

HGCalDetId& HGCalDetId::operator=(const DetId& gen) {
  id_ = gen.rawId();
  return (*this);
}

bool HGCalDetId::isValid(ForwardSubdetector subdet, int zp, int lay, int wafertype, int wafer, int cell) {

  bool ok = ((subdet == HGCEE || subdet == HGCHEF || subdet == HGCHEB) &&
         (cell >= 0 && cell <= kHGCalCellMask) && 
         (wafer >= 0 && wafer <= kHGCalWaferMask) &&
         (wafertype <= kHGCalWaferTypeMask) && 
         (lay >= 0 && lay <= kHGCalLayerMask) &&
         (zp == -1 || zp == 1));
#ifdef DebugLog
  if (!ok) 
    std::cout << "HGCalDetId: subdet " << subdet << ":" 
          << (subdet == HGCEE || subdet == HGCHEF || subdet == HGCHEB) 
          << " Cell " << cell << ":" << (cell >= 0 && cell <= kHGCalCellMask) 
          << " Wafer " << wafer << ":" << (wafer >= 0 && wafer <= kHGCalWaferMask) 
          << " WaferType " << wafertype << ":" << (wafertype <= kHGCalWaferTypeMask)
          << " Layer " << lay << ":" << (lay >= 0 && lay <= kHGCalLayerMask) 
          << " zp " << zp << ":" << (zp == -1 || zp == 1) << std::endl;
#endif
  return ok;
}



std::ostream& operator<<(std::ostream& s,const HGCalDetId& id) {
  return s << "isHGCal=" << id.isHGCal() << " subdet= " << id.subdetId()
       << " zpos=" << id.zside() << " layer=" << id.layer()  
       << " wafer type=" << id.waferType() << " wafer=" << id.wafer() 
       << " cell=" << id.cell();
}