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, hpstanc_transforms::max, and mathSSE::sqrt().

Referenced by BeautifulSoup.BeautifulStoneSoup::endData(), RectangularEtaPhiTrackingRegion::estimator(), intersection(), 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 22 of file VirtualJetProducerHelper.h.

References trackingPlots::common, intersection(), M_PI, hpstanc_transforms::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 32 of file VirtualJetProducerHelper.h.

References trackingPlots::common, M_PI, hpstanc_transforms::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;
}