bubble.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 Bubble: public Topology<TOutput,TMass,TScale>
  {
  public:
    Bubble();  
    ~Bubble(); 

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

    void BB0(vector<TOutput> &res, TScale const& mu2, TMass const& m0, TMass const& m1, TScale const& s) const;

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

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

    void BB3(vector<TOutput> &res, TScale const& mu2, TMass const& s) const;

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

    void BB5(vector<TOutput> &res, TScale const& mu2, TMass const& m0, TMass const& m1) const;
    };
}