reco::helper::VirtualJetProducerHelper Namespace Reference

Functions
double	intersection (double r12)
double	intersection (double r12, double r23, double r13)
double	intersection (double r12, double r23, double r13, double a12, double a23, double a13)

Function Documentation

double reco::helper::VirtualJetProducerHelper::intersection ( double  r12 )

inline

Definition at line 14 of file VirtualJetProducerHelper.h.

References M_PI, SiStripPI::max, and mathSSE::sqrt().

Referenced by compareHistograms::create_dif(), BeautifulSoup.BeautifulStoneSoup::endData(), RectangularEtaPhiTrackingRegion::estimator(), intersection(), FKDTree< TYPE, numberOfDimensions >::search(), edm::serviceregistry::ServicesManager::ServicesManager(), and VirtualJetProducer::writeJets().

                                             {
        if (r12 == 0)
          return M_PI;
        if (r12 >= 2)
          return 0;
        return 2 * std::acos(0.5 * r12) - 0.5 * r12 * sqrt(std::max(0.0, 4 - r12 * r12));
      }

double reco::helper::VirtualJetProducerHelper::intersection ( double  r12, double  r23, double  r13  )

inline

Definition at line 23 of file VirtualJetProducerHelper.h.

References trackingPlots::common, intersection(), M_PI, SiStripPI::max, mathSSE::sqrt(), and groupFilesInBlocks::temp.

                                                                     {
        if (r12 >= 2 || r23 >= 2 || r13 >= 2)
          return 0;
        const double r12_2 = r12 * r12;
        const double r13_2 = r13 * r13;
        const double temp = (r12_2 + r13_2 - r23 * r23);
        const double T2 = std::max(0.0, 4 * r12_2 * r13_2 - temp * temp);
        const double common = 0.5 * (intersection(r12) + intersection(r13) + intersection(r23) - M_PI + 0.5 * sqrt(T2));
        return common;
      }

double reco::helper::VirtualJetProducerHelper::intersection ( double  r12, double  r23, double  r13, double  a12, double  a23, double  a13  )

inline

Definition at line 33 of file VirtualJetProducerHelper.h.

References trackingPlots::common, M_PI, SiStripPI::max, mathSSE::sqrt(), and groupFilesInBlocks::temp.

                                                                                                         {
        if (r12 >= 2 || r23 >= 2 || r13 >= 2)
          return 0;
        const double r12_2 = r12 * r12;
        const double r13_2 = r13 * r13;
        const double temp = (r12_2 + r13_2 - r23 * r23);
        const double T2 = std::max(0.0, 4 * r12_2 * r13_2 - temp * temp);
        const double common = 0.5 * (a12 + a13 + a23 - M_PI + 0.5 * sqrt(T2));
        return common;
      }