EcalTrigTowerDetId.cc

Go to the documentation of this file.

#include "DataFormats/EcalDetId/interface/EcalTrigTowerDetId.h"
#include "FWCore/Utilities/interface/Exception.h"
#include <cassert>

EcalTrigTowerDetId::EcalTrigTowerDetId() {}

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

EcalTrigTowerDetId::EcalTrigTowerDetId(int zside, EcalSubdetector subDet, int i, int j, int mode)
    : DetId(Ecal, EcalTriggerTower) {
  int tower_i = 0;
  int tower_j = 0;

  if (mode == SUBDETIJMODE) {
    tower_i = i;
    tower_j = j;
  } else if (mode == SUBDETDCCTTMODE) {
    throw cms::Exception("InvalidDetId")
        << "EcalTriggerTowerDetId:  Cannot create object. SUBDETDCCTTMODE not yet implemented.";
  } else
    throw cms::Exception("InvalidDetId") << "EcalTriggerTowerDetId:  Cannot create object.  Unknown mode for (int, "
                                            "EcalSubdetector, int, int) constructor.";

  if (tower_i > MAX_I || tower_i < MIN_I || tower_j > MAX_J || tower_j < MIN_J)
    throw cms::Exception("InvalidDetId") << "EcalTriggerTowerDetId:  Cannot create object.  Indexes out of bounds.";

  id_ |= ((zside > 0) ? (0x8000) : (0x0)) | ((subDet == EcalBarrel) ? (0x4000) : (0x0)) | (tower_i << 7) |
         (tower_j & 0x7F);
}

EcalTrigTowerDetId::EcalTrigTowerDetId(const DetId& gen) {
  if (!gen.null() && (gen.det() != Ecal || gen.subdetId() != EcalTriggerTower)) {
    throw cms::Exception("InvalidDetId");
  }
  id_ = gen.rawId();
}

EcalTrigTowerDetId& EcalTrigTowerDetId::operator=(const DetId& gen) {
  if (!gen.null() && (gen.det() != Ecal || gen.subdetId() != EcalTriggerTower)) {
    throw cms::Exception("InvalidDetId");
  }
  id_ = gen.rawId();
  return *this;
}

//New SM numbering scheme. Avoids discontinuity in phi crossing \eta=0
int EcalTrigTowerDetId::iDCC() const {
  if (subDet() == EcalBarrel) {
    //Correction since iphi is uniformized with HB convention
    int iphi_simple = iphi() + 2;
    if (iphi_simple > 72)
      iphi_simple = iphi_simple % 72;
    int id = (iphi_simple - 1) / kEBTowersInPhi + 1;
    if (zside() < 0)
      id += 18;
    return id;
  } else
    throw cms::Exception("MethodNotImplemented") << "EcalTriggerTowerDetId: iDCC not yet implemented";
}

int EcalTrigTowerDetId::iTT() const {
  if (subDet() == EcalBarrel) {
    int ie = ietaAbs() - 1;
    int ip;
    int iphi_simple = iphi() + 2;
    if (iphi_simple > 72)
      iphi_simple = iphi_simple % 72;
    if (zside() < 0) {
      ip = ((iphi_simple - 1) % kEBTowersInPhi) + 1;
    } else {
      ip = kEBTowersInPhi - ((iphi_simple - 1) % kEBTowersInPhi);
    }

    return (ie * kEBTowersInPhi) + ip;
  } else
    throw cms::Exception("MethodNotImplemented") << "EcalTriggerTowerDetId: iTT not yet implemented";
}

int EcalTrigTowerDetId::iquadrant() const {
  if (subDet() == EcalEndcap)
    return int((iphi() - 1) / kEETowersInPhiPerQuadrant) + 1;
  else
    throw cms::Exception("MethodNotApplicable") << "EcalTriggerTowerDetId: iquadrant not applicable";
}

bool EcalTrigTowerDetId::validDetId(int iz, EcalSubdetector sd, int i, int j) {
  return (
      1 == abs(iz) && 0 < i && 0 < j && kEETowersInPhiPerEndcap >= j &&
      ((EcalBarrel == sd && kEBTowersInEta >= i) || (EcalEndcap == sd && kEEOuterEta <= i && kEEInnerEta >= i &&
                                                     (27 > i || ((0 > iz && 0 == j % 2) || (0 < iz && 1 == j % 2))))));
}

int EcalTrigTowerDetId::hashedIndex() const {
  const unsigned int iea(ietaAbs());
  const unsigned int iph(iphi());
  return (subDet() == EcalBarrel ? (iDCC() - 1) * kEBTowersPerSM + iTT() - 1
                                 : kEBTotalTowers + ((zside() + 1) / 2) * kEETowersPerEndcap +
                                       ((iea < 27 ? iea : 27) - kEEOuterEta) * kEETowersInPhiPerEndcap +
                                       (iea < 27 ? iph :  // for iphi=27,28 only half TT present, odd for EE-, even EE+
                                            (iea - 27) * kEETowersInPhiPerEndcap / 2 + (iph + 1) / 2) -
                                       1);
}

EcalTrigTowerDetId EcalTrigTowerDetId::detIdFromDenseIndex(uint32_t di) {
  const EcalSubdetector sd(di < kEBTotalTowers ? EcalBarrel : EcalEndcap);
  const int iz(di < kEBTotalTowers ? (di < kEBHalfTowers ? 1 : -1)
                                   : (di - kEBTotalTowers < kEETowersPerEndcap ? -1 : 1));
  int i;
  int j;
  if (di < kEBTotalTowers)  // barrel
  {
    const unsigned int itt(di % kEBTowersPerSM);
    const unsigned int idc(di / kEBTowersPerSM);
    j = (idc % 18) * kEBTowersInPhi + ((1 + iz) / 2) * kEBTowersInPhi - iz * (itt % kEBTowersInPhi) + 1 - (1 + iz) / 2 -
        2;
    if (j < 1)
      j += 72;
    i = 1 + itt / kEBTowersInPhi;
  } else {
    const int eonly((di - kEBTotalTowers) % kEETowersPerEndcap);
    i = kEEOuterEta + eonly / kEETowersInPhiPerEndcap;
    j = 1 + eonly % kEETowersInPhiPerEndcap;
    if (27 == i)  // last two rings have half of normal phi elementes
    {
      if (j > kEETowersInPhiPerEndcap / 2) {
        ++i;
        j -= kEETowersInPhiPerEndcap / 2;
      }
      j = 2 * j;
      if (0 < iz)
        --j;
    }
  }
  assert(validDetId(iz, sd, i, j));
  return EcalTrigTowerDetId(iz, sd, i, j);
}

#include <ostream>
std::ostream& operator<<(std::ostream& s, const EcalTrigTowerDetId& id) {
  return s << "(EcalTT subDet " << ((id.subDet() == EcalBarrel) ? ("Barrel") : ("Endcap")) << " iz "
           << ((id.zside() > 0) ? ("+ ") : ("- ")) << " ieta " << id.ietaAbs() << " iphi " << id.iphi() << ')';
}