METRP2to2.h

// -*- C++ -*-
//
// METRP2to2.h is a part of Herwig - A multi-purpose Monte Carlo event generator
// Copyright (C) 2009-2019 The Herwig Collaboration
//
// Herwig is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef HERWIG_METRP2to2_H
#define HERWIG_METRP2to2_H
//
// This is the declaration of the METRP2to2 class.
//
// The implementation of this process is based upon hep-ph/0112161 by G.F. Giudice, R. Rattazzi, J.D. Wells.
 
#include "Herwig/MatrixElement/HwMEBase.h"
#include "ThePEG/Repository/UseRandom.h"
#include "Herwig/Utilities/Interpolator.h"
 
 
namespace Herwig {
using namespace ThePEG;
 
class METRP2to2: public HwMEBase {
 
public:
 
  METRP2to2();
 
  virtual unsigned int orderInAlphaS() const { return 0; }
 
  virtual unsigned int orderInAlphaEW() const { return 0; }
 
  virtual Energy2 scale() const;
 
  virtual void getDiagrams() const;
 
  virtual Selector<DiagramIndex> diagrams(const DiagramVector & dv) const;
 
  virtual Selector<const ColourLines *>
  colourGeometries(tcDiagPtr diag) const;
 
  virtual double me2() const;
 
 
 
protected:
 
  virtual void doinit();
 
  virtual void rebind(const TranslationMap & trans);
 
  virtual IVector getReferences();
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
protected:
 
  InvEnergy bccalc(Energy2 s) const;
 
  double A_ny(Energy2 s, Energy2 t) const;
 
  double fpoint(double x) const;
 
  double fnyasympt(double y) const;
 
  void setup_interpolator();
  
private:
 
  METRP2to2 & operator=(const METRP2to2 &) = delete;
 
private:
 
  Interpolator<double, double>::Ptr _interpol;
  
  unsigned int _maxflavour;
 
  unsigned int _ndim;
 
  Energy _planckmass;
 
 
  unsigned int _process;
};
 
}
 
#endif /* HERWIG_METRP2to2_H */