EEDetId.h

Go to the documentation of this file.

#ifndef ECALDETID_EEDETID_H
#define ECALDETID_EEDETID_H

#include <iosfwd>
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/EcalDetId/interface/EcalScDetId.h"
#include "DataFormats/EcalDetId/interface/EcalSubdetector.h"

class EEDetId : public DetId {
public:
  enum {
    Subdet = EcalEndcap
  };

  EEDetId() {}

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

  // fast
  EEDetId(int crystal_ix, int crystal_iy, int iz) : DetId(Ecal, EcalEndcap) {
    id_ |= (crystal_iy & 0x7f) | ((crystal_ix & 0x7f) << 7) | ((iz > 0) ? (0x4000) : (0));
  }
  // slow
  EEDetId(int i, int j, int iz, int mode);

  EEDetId(const DetId& id) : DetId(id) {}

  EEDetId& operator=(const DetId& id) {
    id_ = id.rawId();
    return *this;
  }

  static EcalSubdetector subdet() { return EcalEndcap; }

  int zside() const { return (id_ & 0x4000) ? (1) : (-1); }

  int ix() const { return (id_ >> 7) & 0x7F; }

  int iy() const { return id_ & 0x7F; }

  EcalScDetId sc() const {
    const int scEdge = 5;
    return EcalScDetId(1 + (ix() - 1) / scEdge, 1 + (iy() - 1) / scEdge, zside());
  }

  int isc() const;

  int ic() const;

  int iquadrant() const;

  bool positiveZ() const { return id_ & 0x4000; }

  int iPhiOuterRing() const;  // 1-360 else==0 if not on outer ring!

  static EEDetId idOuterRing(int iPhi, int zEnd);

  int hashedIndex() const {
    const uint32_t jx(ix());
    const uint32_t jd(2 * (iy() - 1) + (jx - 1) / 50);
    return ((positiveZ() ? kEEhalf : 0) + kdi[jd] + jx - kxf[jd]);
  }

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

  EEDetId offsetBy(int nrStepsX, int nrStepsY) const;

  EEDetId switchZSide() const;

  static DetId offsetBy(const DetId startId, int nrStepsX, int nrStepsY);
  static DetId switchZSide(const DetId startId);

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

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

  static bool isNextToBoundary(EEDetId id);

  static bool isNextToDBoundary(EEDetId id);

  static bool isNextToRingBoundary(EEDetId id);

  static EEDetId unhashIndex(int hi);

  static bool validHashIndex(int i) { return (i < kSizeForDenseIndexing); }

  static bool validDetId(int crystal_ix, int crystal_iy, int iz) {
    return crystal_ix >= IX_MIN && crystal_ix <= IX_MAX && crystal_iy >= IY_MIN && crystal_iy <= IY_MAX &&
           std::abs(iz) == 1 && (fastValidDetId(crystal_ix, crystal_iy) || slowValidDetId(crystal_ix, crystal_iy));
  }
  static bool slowValidDetId(int crystal_ix, int crystal_iy);

  static bool fastValidDetId(int crystal_ix, int crystal_iy) {
    float x = crystal_ix;
    float y = crystal_iy;
    float r = (x - 50.5f) * (x - 50.5f) + (y - 50.5f) * (y - 50.5f);
    return r > 12.f * 12.f && r < 48.f * 48.f;
  }

  static int distanceX(const EEDetId& a, const EEDetId& b);

  static int distanceY(const EEDetId& a, const EEDetId& b);

  static int isc(int iscCol,   // output is 1-316
                 int iscRow);  //

  static const int IX_MIN = 1;

  static const int IY_MIN = 1;

  static const int IX_MAX = 100;

  static const int IY_MAX = 100;

  static const int ISC_MIN = 1;

  static const int ICR_MIN = 1;

  static const int ISC_MAX = 316;

  static const int ICR_MAX = 25;

  static constexpr int kEEhalf = 7324;
  static constexpr int kSizeForDenseIndexing = 2 * kEEhalf;

  static const int XYMODE = 0;
  static const int SCCRYSTALMODE = 1;
private:
  bool isOuterRing() const;

  static bool isOuterRingXY(int ax, int ay);

  //Functions from B. Kennedy to retrieve ix and iy from SC and Crystal number

  static const int nCols = 10;
  static const int nCrys = 5; /* Number of crystals per row in SC */
  static const int QuadColLimits[nCols + 1];
  static const int iYoffset[nCols + 1];

  static const unsigned short kxf[2 * IY_MAX];
  static const unsigned short kdi[2 * IY_MAX];

  int ix(int iSC, int iCrys) const;
  int iy(int iSC, int iCrys) const;
  int ixQuadrantOne() const;
  int iyQuadrantOne() const;
};

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

#endif