RealQuadraticEquation.hpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file RealQuadraticEquation.hpp

// This file is part of the Acts project.
//
// Copyright (C) 2017-2018 CERN for the benefit of the Acts project
//
// This Source Code Form is subject to the terms of the Mozilla Public
// License, v. 2.0. If a copy of the MPL was not distributed with this
// file, You can obtain one at http://mozilla.org/MPL/2.0/.

// RealQuadraticEquation.h, Acts project

#pragma once
#include <cmath>
#include <utility>

namespace Acts {
namespace detail {


struct RealQuadraticEquation {
  double first = 0;
  double second = 0;
  int solutions = 0;

  RealQuadraticEquation(double alpha, double beta, double gamma) {
    double discriminant = beta * beta - 4 * alpha * gamma;
    if (discriminant >= 0) {
      solutions = (discriminant == 0) ? 1 : 2;
      double q = -0.5 * (beta + (beta > 0 ? std::sqrt(discriminant)
                                          : -std::sqrt(discriminant)));
      first = q / alpha;
      second = gamma / q;
    }
  }
};

}  // namespace detail
}  // namespace Acts