SonnenscheinEngine.h

Go to the documentation of this file.

/*
   Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
 */
 
#ifndef SONNENSCEINENGINE_H_
#define SONNENSCEINENGINE_H_
 
#include <vector>
#include <string>
#include <utility>
 
#include "TLorentzVector.h"
 
namespace top {
  class SonnenscheinEngine {
  public:
    SonnenscheinEngine();
 
    ~SonnenscheinEngine();
 
    std::vector<std::pair<TLorentzVector, TLorentzVector> > solve(const TLorentzVector& me_t_lp,
                                                                  const TLorentzVector& me_t_b, double mass_t,
                                                                  double mass_wp, const TLorentzVector& me_tbar_lm,
                                                                  const TLorentzVector& me_tbar_bbar, double mass_tbar,
                                                                  double mass_wm, double me_mex, double me_mey) const;
  private:
    std::vector<double> solveQuartic(double h0, double h1, double h2, double h3, double h4) const;
 
    std::vector<double> solveCubic(double s1, double s2, double s3) const;
 
    std::vector<double> solveQuadratic(double a, double b, double c) const;
 
    double sign(double a) const;
 
    bool m_debug;
  };
}
 
#endif