/data/refman/pasoursint/CMSSW_4_1_8_patch12/src/TrackingTools/GeomPropagators/src/RealQuadEquation.h

Go to the documentation of this file.

#ifndef RealQuadEquation_H
#define RealQuadEquation_H

#include <utility>
#include <cmath>

struct RealQuadEquation {

  bool hasSolution;
  double first;
  double second;

  RealQuadEquation( double A, double B, double C) {
    double D = B*B - 4*A*C;
    if (D<0) hasSolution = false;
    else {
      hasSolution = true;
      double q = -0.5*(B + (B>0 ? sqrt(D) : -sqrt(D)));
      first = q/A;
      second = C/q;
    }
  }

};

#endif