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;
 
  enum {
    kEEhalf = 7324,
    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