#include <GlobalDetRodRangeZPhi.h>

Public Types
typedef std::pair< float, float >	Range

Public Member Functions
	GlobalDetRodRangeZPhi (const Plane &RodPlane)

Range	phiRange () const

Range	zRange () const

Private Attributes
Range	thePhiRange

Range	theZRange

Detailed Description

Implementation class for PhiZMeasurementEstimator etc.

Definition at line 12 of file GlobalDetRodRangeZPhi.h.

Member Typedef Documentation

typedef std::pair<float, float> GlobalDetRodRangeZPhi::Range

Definition at line 14 of file GlobalDetRodRangeZPhi.h.

Constructor & Destructor Documentation

GlobalDetRodRangeZPhi::GlobalDetRodRangeZPhi ( const Plane & RodPlane )

Definition at line 8 of file GlobalDetRodRangeZPhi.cc.

References Surface::bounds(), PVValHelper::dx, PVValHelper::dy, PVValHelper::dz, Bounds::length(), SiStripPI::max, min(), PV3DBase< T, PVType, FrameType >::perp(), phi, Geom::phiLess(), phimax, phimin, thePhiRange, theZRange, Bounds::thickness(), Surface::toGlobal(), Bounds::width(), and z.

                                                                {
   float dx = plane.bounds().width() / 2.;
   float dy = plane.bounds().length() / 2.;
   float dz = plane.bounds().thickness() / 2.;
 
   // rods may be inverted (actually are in every other layer), so have to find out the
   // orientation of the local frame
   float deltaZ =
       (plane.toGlobal(LocalPoint(0, 0, -dz)).perp() < plane.toGlobal(LocalPoint(0, 0, dz)).perp()) ? -dz : dz;
 
   const std::array<Surface::GlobalPoint, 4> corners{{plane.toGlobal(LocalPoint(-dx, -dy, deltaZ)),
                                                      plane.toGlobal(LocalPoint(-dx, dy, deltaZ)),
                                                      plane.toGlobal(LocalPoint(dx, -dy, deltaZ)),
                                                      plane.toGlobal(LocalPoint(dx, dy, deltaZ))}};
 
   float phimin = corners[0].phi();
   float phimax = phimin;
 
   float zmin = corners[0].z();
   float zmax = zmin;
 
   for (const auto& corner : corners) {
     float phi = corner.phi();
     if (Geom::phiLess(phi, phimin))
       phimin = phi;
     if (Geom::phiLess(phimax, phi))
       phimax = phi;
 
     float z = corner.z();
     zmin = std::min(zmin, z);
     zmax = std::max(zmax, z);
   }
 
   theZRange.first = zmin;
   theZRange.second = zmax;
   thePhiRange.first = phimin;
   thePhiRange.second = phimax;
 }