EcalTrigTowerDetId.h

Go to the documentation of this file.

#ifndef ECALDETID_ECALTRIGTOWERDETID_H
#define ECALDETID_ECALTRIGTOWERDETID_H

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

class EcalTrigTowerDetId : public DetId {
public:
  EcalTrigTowerDetId();
  EcalTrigTowerDetId(uint32_t rawid);
  EcalTrigTowerDetId(int zside, EcalSubdetector subdet, int i, int j, int mode = SUBDETIJMODE);

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

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

  EcalSubdetector subDet() const { return (id_ & 0x4000) ? EcalBarrel : EcalEndcap; }

  int ietaAbs() const {
    /*       if ( subDet() == EcalBarrel) */
    return (id_ >> 7) & 0x7f;
    /*       else */
    /*  throw(std::runtime_error("EcalTrigTowerDetId: ietaAbs not applicable for this subDetector.")); */
  }

  int ieta() const {
    /*       if ( subDet() == EcalBarrel) */
    return zside() * ietaAbs();
    /*       else */
    /*  throw(std::runtime_error("EcalTrigTowerDetId: ieta not applicable for this subDetector.")); */
  }

  int iphi() const {
    /*       if ( subDet() == EcalBarrel) */
    return id_ & 0x7F;
    /*       else */
    /*  throw(std::runtime_error("EcalTrigTowerDetId: iphi not applicable for this subDetector.")); */
  }

  int iquadrant() const;

  int ix() const {
    if (subDet() == EcalEndcap)
      return (id_ >> 7) & 0x7f;
    else
      throw(std::runtime_error("EcalTrigTowerDetId: ix not applicable for this subDetector."));
  }

  int iy() const {
    if (subDet() == EcalEndcap)
      return id_ & 0x7F;
    else
      throw(std::runtime_error("EcalTrigTowerDetId: ix not applicable for this subDetector."));
  }

  int hashedIndex() const;

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

  static bool validDenseIndex(uint32_t din) { return (din < kSizeForDenseIndexing); }

  static EcalTrigTowerDetId detIdFromDenseIndex(uint32_t di);

  static bool validDetId(int iz, EcalSubdetector sd, int i, int j);

  int iDCC() const;

  int iTT() const;

  static const int MIN_I = 1;
  static const int MIN_J = 1;
  static const int MAX_I = 127;
  static const int MAX_J = 127;

  static const int kEBTowersInPhi = 4;              // per SM (in the Barrel)
  static const int kEBTowersPerSM = 68;             // per SM (in the Barrel)
  static const int kEBTowersInEta = 17;             // per SM (in the Barrel)
                                                    //  static const int kEETowersInEta = 11; // Endcap
  static const int kEETowersInPhiPerQuadrant = 18;  // per Quadrant (in the Endcap)

  // function modes for (int, int) constructor
  static const int SUBDETIJMODE = 0;
  static const int SUBDETDCCTTMODE = 1;

  static constexpr int kEETowersInPhiPerEndcap = 4 * kEETowersInPhiPerQuadrant, kEEOuterEta = 18, kEEInnerEta = 28,
                       kEETowersInEta = (kEEInnerEta - kEEOuterEta + 1), kEBHalfTowers = kEBTowersPerSM * 18,
                       kEBTotalTowers = kEBHalfTowers * 2,
                       kEETowersPerEndcap = kEETowersInEta * kEETowersInPhiPerEndcap - 72,
                       kEETotalTowers = kEETowersPerEndcap * 2, kSizeForDenseIndexing = kEBTotalTowers + kEETotalTowers;
};

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

#endif