#include <DiskSectorBounds.h>

Inheritance diagram for DiskSectorBounds:

Public Member Functions
virtual Bounds *	clone () const

	DiskSectorBounds (float rmin, float rmax, float zmin, float zmax, float phiExt)

float	innerRadius () const

virtual bool	inside (const Local3DPoint &p) const
	Determine if the point is inside the bounds. More...

virtual bool	inside (const Local3DPoint &p, const LocalError &err, float scale) const
	Determine if a point is inside the bounds, taking error into account. More...

virtual float	length () const

float	outerRadius () const

float	phiHalfExtension () const

virtual float	thickness () const

virtual float	width () const

Public Member Functions inherited from Bounds
void	computeSpan (Surface const &plane)

virtual bool	inside (const Local2DPoint &p) const

virtual bool	inside (const Local2DPoint &p, float tollerance) const
	Determine if a point is inside the bounds, taking a tollerance into account. More...

virtual bool	inside (const Local2DPoint &p, const LocalError &err, float scale=1.f) const
	Determine if a 2D point is inside the bounds, taking error into account. More...

std::pair< float, float > const &	phiSpan () const

std::pair< float, float > const &	rSpan () const

virtual float	widthAtHalfLength () const
	Width at half length. Useful for e.g. pitch definition. More...

std::pair< float, float > const &	zSpan () const

virtual	~Bounds ()

Private Attributes
float	theOffset

float	thePhiExtH

float	theRmax

float	theRmin

float	theZmax

float	theZmin

Detailed Description

Definition at line 13 of file DiskSectorBounds.h.

Constructor & Destructor Documentation

DiskSectorBounds::DiskSectorBounds	(	float	rmin,
		float	rmax,
		float	zmin,
		float	zmax,
		float	phiExt
	)

inline

Definition at line 16 of file DiskSectorBounds.h.

References std::swap(), theOffset, thePhiExtH, theRmax, theRmin, theZmax, and theZmin.

Referenced by clone().

                                                                                    : 
      theRmin(rmin), theRmax(rmax), theZmin(zmin), theZmax(zmax), thePhiExtH(0.5f*phiExt) {
      assert(thePhiExtH>0);
      if ( theRmin > theRmax) std::swap( theRmin, theRmax);
      if ( theZmin > theZmax) std::swap( theZmin, theZmax);
      theOffset = theRmin +  0.5f*(theRmax-theRmin);
    }

Member Function Documentation

virtual Bounds* DiskSectorBounds::clone ( void ) const

inlinevirtual

Implements Bounds.

Definition at line 34 of file DiskSectorBounds.h.

References DiskSectorBounds().

                                  { 
      return new DiskSectorBounds(*this);
    }

float DiskSectorBounds::innerRadius ( void ) const

inline

Definition at line 38 of file DiskSectorBounds.h.

References theRmin.

Referenced by BoundDiskSector::innerRadius().

38 {return theRmin;}

DiskSectorBounds::theRmin

float theRmin

Definition: DiskSectorBounds.h:43

bool DiskSectorBounds::inside ( const Local3DPoint & ) const

virtual

Determine if the point is inside the bounds.

Implements Bounds.

Definition at line 5 of file DiskSectorBounds.cc.

References funct::abs(), tmp, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by thickness().

                                                           {
   
   // 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 &	,
		const LocalError &	,
		float	scale
	)		const

virtual

Determine if a point is inside the bounds, taking error into account.

Implements Bounds.

Definition at line 18 of file DiskSectorBounds.cc.

References MuonCkfTrajectoryBuilder_cfi::deltaPhi, deltaR(), f, hpstanc_transforms::max, perp(), perp2(), LocalError::rotate(), mathSSE::sqrt(), tmp, PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), LocalError::xy(), PV3DBase< T, PVType, FrameType >::y(), LocalError::yy(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                               {
 
   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) ;  
 
 }