EcalPreshowerTopology.cc

Go to the documentation of this file.

#include "Geometry/CaloTopology/interface/EcalPreshowerTopology.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include <stdexcept>

ESDetId EcalPreshowerTopology::incrementIy(const ESDetId& id) const {
  ESDetId nextPoint(0);
  if (id == nextPoint)
    return nextPoint;

  //Strips orientend along x direction for plane 2
  if (id.plane() == 2) {
    if (id.strip() < 32) {
      //Incrementing just strip number
      if (ESDetId::validDetId(id.strip() + 1, id.six(), id.siy(), id.plane(), id.zside()))
        nextPoint = ESDetId(id.strip() + 1, id.six(), id.siy(), id.plane(), id.zside());
    } else {
      //Changing wafer
      if (ESDetId::validDetId(1, id.six(), id.siy() + 1, id.plane(), id.zside()))
        nextPoint = ESDetId(1, id.six(), id.siy() + 1, id.plane(), id.zside());
    }
  }
  //Strips orientend along y direction for plane 1
  else if (id.plane() == 1) {
    //Changing wafer
    if (ESDetId::validDetId(id.strip(), id.six(), id.siy() + 1, id.plane(), id.zside()))
      nextPoint = ESDetId(id.strip(), id.six(), id.siy() + 1, id.plane(), id.zside());
  }

  return nextPoint;
}

ESDetId EcalPreshowerTopology::decrementIy(const ESDetId& id) const {
  ESDetId nextPoint(0);
  if (id == nextPoint)
    return nextPoint;

  //Strips orientend along x direction for plane 2
  if (id.plane() == 2) {
    if (id.strip() > 1) {
      //Decrementing just strip number
      if (ESDetId::validDetId(id.strip() - 1, id.six(), id.siy(), id.plane(), id.zside()))
        nextPoint = ESDetId(id.strip() - 1, id.six(), id.siy(), id.plane(), id.zside());
    } else {
      //Changing wafer
      if (ESDetId::validDetId(32, id.six(), id.siy() - 1, id.plane(), id.zside()))
        nextPoint = ESDetId(32, id.six(), id.siy() - 1, id.plane(), id.zside());
    }
  }
  //Strips orientend along y direction for plane 1
  else if (id.plane() == 1) {
    //Changing wafer
    if (ESDetId::validDetId(id.strip(), id.six(), id.siy() - 1, id.plane(), id.zside()))
      nextPoint = ESDetId(id.strip(), id.six(), id.siy() - 1, id.plane(), id.zside());
  }
  return nextPoint;
}

ESDetId EcalPreshowerTopology::incrementIx(const ESDetId& id) const {
  ESDetId nextPoint(0);
  if (id == nextPoint)
    return nextPoint;

  //Strips orientend along x direction for plane 2
  if (id.plane() == 2) {
    //Changing wafer
    if (ESDetId::validDetId(id.strip(), id.six() + 1, id.siy(), id.plane(), id.zside()))
      nextPoint = ESDetId(id.strip(), id.six() + 1, id.siy(), id.plane(), id.zside());
  }
  //Strips orientend along y direction for plane 1
  else if (id.plane() == 1) {
    if (id.strip() < 32) {
      //Incrementing just strip number
      if (ESDetId::validDetId(id.strip() + 1, id.six(), id.siy(), id.plane(), id.zside()))
        nextPoint = ESDetId(id.strip() + 1, id.six(), id.siy(), id.plane(), id.zside());
    } else {
      //Changing wafer
      if (ESDetId::validDetId(1, id.six() + 1, id.siy(), id.plane(), id.zside()))
        nextPoint = ESDetId(1, id.six() + 1, id.siy(), id.plane(), id.zside());
    }
  }

  return nextPoint;
}

ESDetId EcalPreshowerTopology::decrementIx(const ESDetId& id) const {
  ESDetId nextPoint(0);
  if (id == nextPoint)
    return nextPoint;

  //Strips orientend along x direction for plane 2
  if (id.plane() == 2) {
    //Changing wafer
    if (ESDetId::validDetId(id.strip(), id.six() - 1, id.siy(), id.plane(), id.zside()))
      nextPoint = ESDetId(id.strip(), id.six() - 1, id.siy(), id.plane(), id.zside());
  }
  //Strips orientend along y direction for plane 1
  else if (id.plane() == 1) {
    if (id.strip() > 1) {
      //Decrementing just strip number
      if (ESDetId::validDetId(id.strip() - 1, id.six(), id.siy(), id.plane(), id.zside()))
        nextPoint = ESDetId(id.strip() - 1, id.six(), id.siy(), id.plane(), id.zside());
    } else {
      //Changing wafer
      if (ESDetId::validDetId(32, id.six() - 1, id.siy(), id.plane(), id.zside()))
        nextPoint = ESDetId(32, id.six() - 1, id.siy(), id.plane(), id.zside());
    }
  }
  return nextPoint;
}

ESDetId EcalPreshowerTopology::incrementIz(const ESDetId& id) const {
  ESDetId nextPoint(0);
  if (id == nextPoint)
    return nextPoint;

  if (ESDetId::validDetId(id.strip(), id.six(), id.siy(), id.plane() + 1, id.zside()))
    nextPoint = ESDetId(id.strip(), id.six(), id.siy(), id.plane() + 1, id.zside());

  return nextPoint;
}

ESDetId EcalPreshowerTopology::decrementIz(const ESDetId& id) const {
  ESDetId nextPoint(0);
  if (id == nextPoint)
    return nextPoint;

  if (ESDetId::validDetId(id.strip(), id.six(), id.siy(), id.plane() - 1, id.zside()))
    nextPoint = ESDetId(id.strip(), id.six(), id.siy(), id.plane() - 1, id.zside());

  return nextPoint;
}