triangle.h

//
// QCDLoop 2016
//
// Authors: Stefano Carrazza: stefano.carrazza@cern.ch
//          Keith Ellis: keith.ellis@durham.ac.uk
//          Giulia Zanderighi: giulia.zanderighi@cern.ch

#pragma once

#include "topology.h"

namespace ql
{
  template<typename TOutput = complex, typename TMass = double, typename TScale = double>
  class Triangle: public Topology<TOutput,TMass,TScale>
  {
  public:
    Triangle();  
    ~Triangle(); 

    void integral(vector<TOutput> &res, TScale const& mu2, vector<TMass> const& m, vector<TScale> const& p);

    void T0(vector<TOutput> &res, TMass const (&xpi)[6], int const& massive) const;

    void T1(vector<TOutput> &res, TScale const& mu2, TScale const& p) const;

    void T2(vector<TOutput> &res, TScale const& mu2, TScale const& p2, TScale const& p3) const;

    void T3(vector<TOutput> &res, TScale const& mu2, TMass const& m, TScale const& p1, TScale const& p2) const;

    void T4(vector<TOutput> &res, TScale const& mu2, TMass const& m, TScale const& p2) const;

    void T5(vector<TOutput> &res, TScale const& mu2, TMass const& m) const;

    void T6(vector<TOutput> &res, TScale const& mu2, TMass const& m2, TMass const& m3, TScale const& p2) const;

  private:
    void TriSort(TScale (&psq)[3], TMass (&msq)[3]) const;

    void TriSort2(TMass const (&xpi)[6], TMass (&ypi)[6]) const;

    void SnglSort(TScale (&psq)[3]) const;

    void TIN0(TOutput &res, TMass const (&xpi)[6]) const;

    void TIN1(TOutput &res, TMass const (&xpi)[6], TMass const (&sxpi)[6], int const& massive) const;

    void TIN2(TOutput &res, TMass const (&xpi)[6], TMass const (&sxpi)[6], int const& massive) const;

    void TIN3(TOutput &res, TMass const (&xpi)[6], TMass const (&sxpi)[6], int const& massive) const;

    void TINDNS(TOutput &res, TMass const (&xpi)[6]) const;

    void TINDNS1(TOutput &res, TMass const (&xpi)[6]) const;

    void TINDNS2(TOutput &res, TMass const (&xpi)[6]) const;

    TOutput Kallen(TOutput const& p1, TOutput const& p2, TOutput const& p3) const;

    TOutput Kallen2(TOutput const& p1, TOutput const& p2, TOutput const& p3) const;

  };
}