ESDetId.h

Go to the documentation of this file.

#ifndef ECALDETID_ESDETID_H
#define ECALDETID_ESDETID_H

#include <iosfwd>
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
#include "FWCore/Utilities/interface/Exception.h"

class ESDetId : public DetId {
public:
  enum { Subdet = EcalPreshower };
  ESDetId() {}
  ESDetId(uint32_t rawid) : DetId(rawid) {}

  ESDetId(int strip, int ixs, int iys, int plane, int iz, bool doverify = false) : DetId(Ecal, EcalPreshower) {
    id_ |= (strip & 0x3F) | ((ixs & 0x3F) << 6) | ((iys & 0x3F) << 12) | (((plane - 1) & 0x1) << 18) |
           ((iz > 0) ? (1 << 19) : (0));
    if (doverify)
      verify(strip, ixs, iys, plane, iz);
  }

  ESDetId(const DetId& id);
  ESDetId& operator=(const DetId& id);

  EcalSubdetector subdet() const { return EcalSubdetector(subdetId()); }
  int zside() const { return (id_ & 0x80000) ? (1) : (-1); }
  int plane() const { return ((id_ >> 18) & 0x1) + 1; }
  int six() const { return (id_ >> 6) & 0x3F; }
  int siy() const { return (id_ >> 12) & 0x3F; }
  int strip() const { return (id_ & 0x3F); }
  int hashedIndex() const;

  uint32_t denseIndex() const { return hashedIndex(); }

  static bool validDenseIndex(uint32_t din) { return validHashIndex(din); }

  static ESDetId detIdFromDenseIndex(uint32_t din) { return unhashIndex(din); }

  static ESDetId unhashIndex(int hi);
  static bool validHashIndex(int hi) { return (hi < kSizeForDenseIndexing); }
  static bool validDetId(int istrip, int ixs, int iys, int iplane, int iz);
  static void verify(int istrip, int ixs, int iys, int iplane, int iz);

  static const int IX_MIN = 1;
  static const int IY_MIN = 1;
  static const int IX_MAX = 40;
  static const int IY_MAX = 40;
  static const int ISTRIP_MIN = 1;
  static const int ISTRIP_MAX = 32;
  static const int PLANE_MIN = 1;
  static const int PLANE_MAX = 2;
  static const int IZ_NUM = 2;

private:
  enum {
    kXYMAX = 1072,
    kXYMIN = 1,
    kXMAX = 40,
    kYMAX = 40,
    kXMIN = 1,
    kYMIN = 1,
    // now normalize to A-D notation for ease of use
    kNa = IZ_NUM,
    kNb = PLANE_MAX - PLANE_MIN + 1,
    kNc = kXYMAX - kXYMIN + 1,
    kNd = ISTRIP_MAX - ISTRIP_MIN + 1,
    kLd = kNd,
    kLc = kLd * kNc,
    kLb = kLc * kNb,
    kLa = kLb * kNa
  };

  static const unsigned short hxy1[kXMAX][kYMAX];
  static const unsigned short hx1[kXYMAX];
  static const unsigned short hy1[kXYMAX];

  static const unsigned short hxy2[kXMAX][kYMAX];
  static const unsigned short hx2[kXYMAX];
  static const unsigned short hy2[kXYMAX];

public:
  static constexpr int kSizeForDenseIndexing = kLa;
};

std::ostream& operator<<(std::ostream&, const ESDetId& id);

inline ESDetId::ESDetId(const DetId& gen) : DetId(gen) {
#ifdef EDM_ML_DEBUG
  if (!gen.null() && (gen.det() != Ecal || gen.subdetId() != EcalPreshower)) {
    throw cms::Exception("InvalidDetId");
  }
#endif
}

inline ESDetId& ESDetId::operator=(const DetId& gen) {
#ifdef EDM_ML_DEBUG
  if (!gen.null() && (gen.det() != Ecal || gen.subdetId() != EcalPreshower)) {
    throw cms::Exception("InvalidDetId");
  }
#endif
  id_ = gen.rawId();
  return *this;
}

#endif