Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 12:53:33

0001 #ifndef GeometryVector_PhiInterval_H
0002 #define GeometryVector_PhiInterval_H
0003 #include "DataFormats/GeometryVector/interface/Basic3DVector.h"
0004 #include "DataFormats/Math/interface/normalizedPhi.h"
0005 
0006 class PhiInterval {
0007 public:
0008   PhiInterval(float phi1, float phi2) {
0009     phi2 = proxim(phi2, phi1);
0010     constexpr float c1 = 2. * M_PI;
0011     if (phi2 < phi1)
0012       phi2 += c1;
0013     auto dphi = 0.5f * (phi2 - phi1);
0014     auto phi = phi1 + dphi;
0015     x = std::cos(phi);
0016     y = std::sin(phi);
0017     dcos = std::cos(dphi);
0018   }
0019 
0020   PhiInterval(float ix, float iy, float dphi) {
0021     auto norm = 1.f / std::sqrt(ix * ix + iy * iy);
0022     x = ix * norm;
0023     y = iy * norm;
0024     dcos = std::cos(dphi);
0025   }
0026 
0027   template <typename T>
0028   bool inside(Basic3DVector<T> const& v) const {
0029     return inside(v.x(), v.y());
0030   }
0031 
0032   bool inside(float ix, float iy) const { return ix * x + iy * y > dcos * std::sqrt(ix * ix + iy * iy); }
0033 
0034 private:
0035   float x, y;
0036   float dcos;
0037 };
0038 
0039 #endif