EcalEtaPhiRegion.cc

Go to the documentation of this file.

#include "RecoEcal/EgammaCoreTools/interface/EcalEtaPhiRegion.h"

EcalEtaPhiRegion::EcalEtaPhiRegion(double etaLow, double etaHigh, double phiLow, double phiHigh):
  etaLow_(etaLow),etaHigh_(etaHigh),phiLow_(phiLow),phiHigh_(phiHigh)
{
  // put phi in range -pi to pi
  if(phiLow_ > Geom::pi()) phiLow_ -= Geom::twoPi();
  if(phiLow_ < -Geom::pi()) phiLow_ += Geom::twoPi();
  if(phiHigh_ > Geom::pi()) phiHigh_ -= Geom::twoPi();
  if(phiHigh_ < -Geom::pi()) phiHigh_ += Geom::twoPi();
}

bool EcalEtaPhiRegion::inRegion(const GlobalPoint& position) const {

  double phi = position.phi();
  double phihightemp = phiHigh_;

  if(phihightemp<phiLow_) phihightemp += Geom::twoPi();

  // put phi in range -pi to pi
  if(phi > Geom::pi()) phi -= Geom::twoPi();
  if(phi < -Geom::pi()) phi += Geom::twoPi();
  if(phi<phiLow_) phi += Geom::twoPi();


  return (position.eta() > etaLow_ && position.eta() < etaHigh_ &&
      phi > phiLow_ && phi < phihightemp);

}