EcalBarrelTopology.cc

Go to the documentation of this file.

#include "Geometry/CaloTopology/interface/EcalBarrelTopology.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"


EBDetId EcalBarrelTopology::incrementIeta(const EBDetId& id) const {
  if (!(*theGeom_).getSubdetectorGeometry(DetId::Ecal,EcalBarrel)->present(id))
    return EBDetId(0);
      
  EBDetId nextPoint;
  if (id.ieta()==-1) 
    {
      if (EBDetId::validDetId(1,id.iphi()))
    nextPoint=EBDetId (1,id.iphi());
      else
    return EBDetId(0);
    }
  else
    {
      if (EBDetId::validDetId(id.ieta()+1,id.iphi()))
    nextPoint=EBDetId(id.ieta()+1,id.iphi());
      else
    return EBDetId(0);
    }
  if ((*theGeom_).getSubdetectorGeometry(DetId::Ecal,EcalBarrel)->present(nextPoint))
    return nextPoint;
  else
    return EBDetId(0);
}

EBDetId EcalBarrelTopology::decrementIeta(const EBDetId& id) const {
  
  if (!(*theGeom_).getSubdetectorGeometry(DetId::Ecal,EcalBarrel)->present(id))
    return EBDetId(0);
  
  EBDetId nextPoint;
  if (id.ieta()==1)
    { 
      if (EBDetId::validDetId(-1,id.iphi()))
    nextPoint=EBDetId(-1,id.iphi());
      else
    return EBDetId(0);
    }
  else
    { 
      if (EBDetId::validDetId(id.ieta()-1,id.iphi()))
    nextPoint=EBDetId(id.ieta()-1,id.iphi());
      else
    return EBDetId(0);
    }
  
  if ((*theGeom_).getSubdetectorGeometry(DetId::Ecal,EcalBarrel)->present(nextPoint))
    return nextPoint;
  else
    return EBDetId(0);
} 


EBDetId EcalBarrelTopology::incrementIphi(const EBDetId& id) const {
  if (!(*theGeom_).getSubdetectorGeometry(DetId::Ecal,EcalBarrel)->present(id))
    return EBDetId(0);
  
  EBDetId nextPoint;
  
  if (id.iphi()==EBDetId::MAX_IPHI) 
    {
      if (EBDetId::validDetId(id.ieta(),EBDetId::MIN_IPHI))
    nextPoint=EBDetId(id.ieta(),EBDetId::MIN_IPHI);
      else
    return EBDetId(0);
    }
  else
    {
      if (EBDetId::validDetId(id.ieta(),id.iphi()+1)) 
    nextPoint=EBDetId(id.ieta(),id.iphi()+1);
      else
    return EBDetId(0);
    }
  
  if ((*theGeom_).getSubdetectorGeometry(DetId::Ecal,EcalBarrel)->present(nextPoint))
    return nextPoint;
  else
    return EBDetId(0);
} 


EBDetId EcalBarrelTopology::decrementIphi(const EBDetId& id) const {
  if (!(*theGeom_).getSubdetectorGeometry(DetId::Ecal,EcalBarrel)->present(id))
    return EBDetId(0);
  
  EBDetId nextPoint;
  
  if (id.iphi()==EBDetId::MIN_IPHI)
    { 
      if (EBDetId::validDetId(id.ieta(),EBDetId::MAX_IPHI))
    nextPoint=EBDetId(id.ieta(),EBDetId::MAX_IPHI);
      else
    return EBDetId(0);
    }
  else
    {
      if (EBDetId::validDetId(id.ieta(),id.iphi()-1))
    nextPoint=EBDetId(id.ieta(),id.iphi()-1);
      else
        return EBDetId(0);
    }  
  
  if ((*theGeom_).getSubdetectorGeometry(DetId::Ecal,EcalBarrel)->present(nextPoint))
    return nextPoint;
  else
    return EBDetId(0);
}