lst_math Namespace Reference

Classes
class	Helix
class	Hit

Functions
float	ATan2 (float y, float x)
Hit	getCenterFromThreePoints (Hit &hitA, Hit &hitB, Hit &hitC)
float	Phi_mpi_pi (float x)
float	ptEstimateFromRadius (float radius)

Function Documentation

ATan2()

float lst_math::ATan2 ( float  y, float  x  )

inline

Definition at line 23 of file lst_math.h.

References Matriplex::atan2(), M_PI, and x.

Referenced by CosmicGenFilterHelix::monitorEnd(), CosmicGenFilterHelix::monitorStart(), lst_math::Hit::setDerivedQuantities(), and metsig::significanceAlgo::significance().

                                       {
    if (x != 0)
      return atan2(y, x);
    if (y == 0)
      return 0;
    if (y > 0)
      return M_PI / 2;
    else
      return -M_PI / 2;
  };

getCenterFromThreePoints()

Hit lst_math::getCenterFromThreePoints ( Hit &  hitA, Hit &  hitB, Hit &  hitC  )

inline

Definition at line 154 of file lst_math.h.

References gather_cfg::cout, or, x, lst_math::Hit::x(), and lst_math::Hit::y().

Referenced by setGnnNtupleBranches().

                                                                       {
    //       C
    //
    //
    //
    //    B           d <-- find this point that makes the arc that goes throw a b c
    //
    //
    //     A

    // Steps:
    // 1. Calculate mid-points of lines AB and BC
    // 2. Find slopes of line AB and BC
    // 3. construct a perpendicular line between AB and BC
    // 4. set the two equations equal to each other and solve to find intersection

    float xA = hitA.x();
    float yA = hitA.y();
    float xB = hitB.x();
    float yB = hitB.y();
    float xC = hitC.x();
    float yC = hitC.y();

    float x_mid_AB = (xA + xB) / 2.;
    float y_mid_AB = (yA + yB) / 2.;

    //float x_mid_BC = (xB + xC) / 2.;
    //float y_mid_BC = (yB + yC) / 2.;

    bool slope_AB_inf = (xB - xA) == 0;
    bool slope_BC_inf = (xC - xB) == 0;

    bool slope_AB_zero = (yB - yA) == 0;
    bool slope_BC_zero = (yC - yB) == 0;

    float slope_AB = slope_AB_inf ? 0 : (yB - yA) / (xB - xA);
    float slope_BC = slope_BC_inf ? 0 : (yC - yB) / (xC - xB);

    float slope_perp_AB = slope_AB_inf or slope_AB_zero ? 0. : -1. / (slope_AB);
    //float slope_perp_BC = slope_BC_inf or slope_BC_zero ? 0. : -1. / (slope_BC);

    if ((slope_AB - slope_BC) == 0) {
      std::cout << " slope_AB_zero: " << slope_AB_zero << std::endl;
      std::cout << " slope_BC_zero: " << slope_BC_zero << std::endl;
      std::cout << " slope_AB_inf: " << slope_AB_inf << std::endl;
      std::cout << " slope_BC_inf: " << slope_BC_inf << std::endl;
      std::cout << " slope_AB: " << slope_AB << std::endl;
      std::cout << " slope_BC: " << slope_BC << std::endl;
      std::cout << "MathUtil::getCenterFromThreePoints() function the three points are in straight line!" << std::endl;
      return Hit();
    }

    float x =
        (slope_AB * slope_BC * (yA - yC) + slope_BC * (xA + xB) - slope_AB * (xB + xC)) / (2. * (slope_BC - slope_AB));
    float y = slope_perp_AB * (x - x_mid_AB) + y_mid_AB;

    return Hit(x, y, 0);
  };

Phi_mpi_pi()

float lst_math::Phi_mpi_pi ( float  x )

inline

Definition at line 8 of file lst_math.h.

References gather_cfg::cout, CommonMethods::isnan(), M_PI, and x.

Referenced by TTbarSpinCorrHepMCAnalyzer::analyze(), L1MetPfProducer::CalcMetHLS(), L1TStage2CaloLayer2Offline::fillEnergySums(), EcalClusterCrackCorrection::getValue(), lst_math::Helix::Helix(), reco::MustacheKernel::inDynamicDPhiWindow(), reco::MustacheKernel::inMustache(), egammaTools::localEcalClusterCoordsEB(), egammaTools::localEcalClusterCoordsEE(), EcalBasicClusterLocalContCorrection::operator()(), and lst_math::Hit::setDerivedQuantities().

                                   {
    if (std::isnan(x)) {
      std::cout << "MathUtil::Phi_mpi_pi() function called with NaN" << std::endl;
      return x;
    }

    while (x >= M_PI)
      x -= 2. * M_PI;

    while (x < -M_PI)
      x += 2. * M_PI;

    return x;
  };

ptEstimateFromRadius()

float lst_math::ptEstimateFromRadius ( float  radius )

inline

Definition at line 34 of file lst_math.h.

References CosmicsPD_Skims::radius.

Referenced by setGnnNtupleBranches().

{ return 2.99792458e-3 * 3.8 * radius; };