MEGG2GG.h

// -*- C++ -*-
//
// MEGG2GG.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2019 Leif Lonnblad
//
// ThePEG is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef ThePEG_MEGG2GG_H
#define ThePEG_MEGG2GG_H
// This is the declaration of the MEGG2GG class.
 
#include "ThePEG/MatrixElement/ME2to2QCD.h"
 
namespace ThePEG {
 
class MEGG2GG: public ME2to2QCD {
 
public:
 
  virtual double me2() const;
 
  virtual void getDiagrams() const;
 
  virtual Selector<const ColourLines *>
  colourGeometries(tcDiagPtr diag) const;
 
  virtual Selector<DiagramIndex> diagrams(const DiagramVector & dv) const;
 
protected:
 
  double colA1() const { return sqr(1.0 + sHat()/tHat()) + 0.5; }
 
  double colB1() const { return sqr(1.0 + uHat()/sHat()) + 0.5; }
 
  double colC1() const { return sqr(1.0 + tHat()/uHat()) + 0.5; }
 
  double colA2() const { return 0.5 + sqr(1.0 + tHat()/sHat()); }
 
  double colB2() const { return 0.5 + sqr(1.0 + sHat()/uHat()); }
 
  double colC2() const { return 0.5 + sqr(1.0 + uHat()/tHat()); }
 
public:
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
private:
 
  static NoPIOClassDescription<MEGG2GG> initMEGG2GG;
 
  MEGG2GG & operator=(const MEGG2GG &) = delete;
 
};
 
}
 
 
namespace ThePEG {
 
template <>
struct BaseClassTrait<MEGG2GG,1>: public ClassTraitsType {
  typedef ME2to2QCD NthBase;
};
 
template <>
struct ClassTraits<MEGG2GG>: public ClassTraitsBase<MEGG2GG> {
  static string className() { return "ThePEG::MEGG2GG"; }
  static string library() { return "MEQCD.so"; }
};
 
}
 
#endif /* ThePEG_MEGG2GG_H */