tkDetUtil Namespace Reference

Functions
float	calculatePhiWindow (const MeasurementEstimator::Local2DVector &imaxDistance, const TrajectoryStateOnSurface &ts, const Plane &plane)
float	computeWindowSize (const GeomDet *det, const TrajectoryStateOnSurface &tsos, const MeasurementEstimator &est)
bool	overlapInPhi (float phi, const GeomDet &det, float phiWindow)
bool	overlapInPhi (GlobalPoint crossPoint, const GeomDet &det, float phiWindow)

Function Documentation

float tkDetUtil::calculatePhiWindow	(	const MeasurementEstimator::Local2DVector &	imaxDistance,
		const TrajectoryStateOnSurface &	ts,
		const Plane &	plane
	)

Definition at line 20 of file TkDetUtil.cc.

References funct::abs(), PV3DBase< T, PVType, FrameType >::barePhi(), constexpr, f, mps_fire::i, likely, TrajectoryStateOnSurface::localPosition(), M_PI, SiStripPI::max, particleFlowClusterHGC_cfi::maxDistance, min(), Plane::normalVector(), Geom::phiLess(), phimax, phimin, trackingPOGFilters_cfi::phiWindow, Geom::pi(), mathSSE::sqrt(), Surface::toGlobal(), RecoTauCleanerPlugins::tolerance, GloballyPositioned< T >::toLocal(), PV2DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::x(), globals_cff::x1, PV2DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by computeWindowSize(), and overlapInPhi().

  {
    MeasurementEstimator::Local2DVector maxDistance(std::abs(imaxDistance.x()),std::abs(imaxDistance.y()));

    constexpr float tolerance=1.e-6;
    LocalPoint start = ts.localPosition();
    //     std::cout << "plane z " << plane.normalVector() << std::endl;
    float dphi=0;
    if likely(std::abs(1.f-std::abs(plane.normalVector().z()))<tolerance) {
      auto ori = plane.toLocal(GlobalPoint(0.,0.,0.));
      auto xc = std::abs(start.x() - ori.x());
      auto yc = std::abs(start.y() - ori.y());
      
      if (yc<maxDistance.y() && xc<maxDistance.x()) return M_PI;

      auto hori = yc>maxDistance.y(); // quadrant 1 (&2), otherwiase quadrant 1&4
      auto y0 = hori ? yc + std::copysign(maxDistance.y(), xc - maxDistance.x()) : xc - maxDistance.x();
      auto x0 = hori ? xc - maxDistance.x() : -yc - maxDistance.y();
      auto y1 = hori ? yc - maxDistance.y() :  xc - maxDistance.x();
      auto x1 = hori ? xc + maxDistance.x() : -yc + maxDistance.y();

      auto sp = (x0*x1+y0*y1)/std::sqrt((x0*x0+y0*y0)*(x1*x1+y1*y1));
      sp = std::min(std::max(sp,-1.f),1.f);
      dphi = std::acos(sp);      

      return dphi;
    }
    

    // generic algo
    float corners[]  =  { plane.toGlobal(LocalPoint( start.x()+maxDistance.x(), start.y()+maxDistance.y() )).barePhi(),
              plane.toGlobal(LocalPoint( start.x()-maxDistance.x(), start.y()+maxDistance.y() )).barePhi(),
              plane.toGlobal(LocalPoint( start.x()-maxDistance.x(), start.y()-maxDistance.y() )).barePhi(),
              plane.toGlobal(LocalPoint( start.x()+maxDistance.x(), start.y()-maxDistance.y() )).barePhi() 
    };
    
    float phimin = corners[0];
    float phimax = phimin;
    for ( int i = 1; i<4; i++) {
      float cPhi = corners[i];
      if ( Geom::phiLess(cPhi, phimin) ) { phimin = cPhi; }
      if ( Geom::phiLess( phimax, cPhi) ) { phimax = cPhi; }
    }
    float phiWindow = phimax - phimin;
    if ( phiWindow < 0.) { phiWindow +=  2.*Geom::pi();}
    // std::cout << "phiWindow " << phiWindow << ' ' << dphi << ' ' << dphi-phiWindow  << std::endl;
    return phiWindow;
  }

float tkDetUtil::computeWindowSize	(	const GeomDet *	det,
		const TrajectoryStateOnSurface &	tsos,
		const MeasurementEstimator &	est
	)

Definition at line 9 of file TkDetUtil.cc.

References calculatePhiWindow(), particleFlowClusterHGC_cfi::maxDistance, MeasurementEstimator::maximalLocalDisplacement(), and GeomDet::surface().

Referenced by TECLayer::groupedCompatibleDetsV(), TIDRing::groupedCompatibleDetsV(), CompositeTECWedge::groupedCompatibleDetsV(), Phase2EndcapRing::groupedCompatibleDetsV(), and overlapInPhi().

  {
    const Plane& startPlane = det->surface();  
    auto maxDistance =  est.maximalLocalDisplacement( tsos, startPlane);
    return std::copysign(calculatePhiWindow(maxDistance, tsos, startPlane),maxDistance.x()); 
  }

bool tkDetUtil::overlapInPhi ( float  phi, const GeomDet &  det, float  phiWindow  )

inline

Definition at line 19 of file TkDetUtil.h.

References Surface::phiSpan(), rangesIntersect(), and GeomDet::surface().

Referenced by overlapInPhi(), CompositeTECWedge::searchNeighbors(), TIDRing::searchNeighbors(), and Phase2EndcapRing::searchNeighbors().

                                                                      {
    std::pair<float,float> phiRange(phi-phiWindow, phi+phiWindow);
    return rangesIntersect( phiRange, det.surface().phiSpan(), PhiLess());
  }

bool tkDetUtil::overlapInPhi ( GlobalPoint  crossPoint, const GeomDet &  det, float  phiWindow  )

inline

Definition at line 26 of file TkDetUtil.h.

References PV3DBase< T, PVType, FrameType >::barePhi(), calculatePhiWindow(), computeWindowSize(), particleFlowClusterHGC_cfi::maxDistance, overlapInPhi(), and trackingPOGFilters_cfi::phiWindow.

                                                                                  {
    return overlapInPhi(crossPoint.barePhi(), det,phiWindow);
  }