DiskSectorBounds.cc

Go to the documentation of this file.

#include "DiskSectorBounds.h"

using namespace std;

bool DiskSectorBounds::inside( const Local3DPoint& p) const {
  
  // transform to system with local frame situated at disk center
  // and rotated  x/y axis
  Local3DPoint tmp( p.y()+theOffset, -p.x(), p.z());
  
  return 
    ( (tmp.z() >= theZmin) & (tmp.z() <= theZmax)  &
      (tmp.perp2() >= theRmin*theRmin) & (tmp.perp2() <= theRmax*theRmax) 
      ) &&
    (std::abs(tmp.barePhi()) <= thePhiExtH);
}

bool DiskSectorBounds::inside( const Local3DPoint& p, const LocalError& err, float scale) const {

  if ( (p.z() < theZmin) | (p.z() > theZmax)) return false;

  Local2DPoint tmp( p.x(), p.y()+ theOffset);
  auto perp2 = tmp.mag2();
  auto perp = std::sqrt(perp2);

  // this is not really correct, should consider also the angle of the error ellipse
   if (perp2 == 0) return scale*scale*err.yy()   > theRmin*theRmin;

   LocalError tmpErr( err.xx(), err.xy(), err.yy());
   LocalError rotatedErr = tmpErr.rotate(tmp.x(), tmp.y());
   // x direction in this system is now r, phi is y
 
   float deltaR = scale*std::sqrt(rotatedErr.xx());
   float deltaPhi = std::atan( scale*std::sqrt(rotatedErr.yy()/perp2));

   float tmpPhi = std::acos( tmp.y() / perp);

   // the previous code (tmpPhi <= thePhiExt + deltaPhi) was wrong
   return ( (perp >= std::max(theRmin-deltaR, 0.f)) & (perp <= theRmax+deltaR) )
     && (tmpPhi <= thePhiExtH + deltaPhi) ;  

}