DipoleMIOperator.h

// -*- C++ -*-
//
// DipoleMIOperator.h is a part of Herwig - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-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_DipoleMIOperator_H
#define HERWIG_DipoleMIOperator_H
//
// This is the declaration of the DipoleMIOperator class.
//
 
#include "Herwig/MatrixElement/Matchbox/InsertionOperators/MatchboxInsertionOperator.h"
#include "Herwig/MatrixElement/Matchbox/Base/MatchboxMEBase.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class DipoleMIOperator: public MatchboxInsertionOperator {
 
public:
 
  DipoleMIOperator();
 
  virtual ~DipoleMIOperator();
 
public:
 
  virtual void setXComb(tStdXCombPtr xc);
 
  virtual bool apply(const cPDVector&) const;
 
  bool apply(tcPDPtr) const;
 
  vector<int> NLightJetVec() const;
 
  vector<int> NHeavyJetVec() const;
 
  vector<int> NLightBornVec() const;
 
  vector<int> NHeavyBornVec() const;
 
  vector<int> NLightProtonVec() const;
 
  virtual double me2() const;
 
  virtual double oneLoopDoublePole() const;
 
  virtual double oneLoopSinglePole() const;
 
public:
  
  template <class T>
  inline T rootOfKallen (T a, T b, T c) const {
    return sqrt( a*a + b*b + c*c - 2.*( a*b+a*c+b*c ) ); }
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
private:
 
  double CA;
 
  double CF;
 
  double gammaQuark;
 
  double gammaGluon;
  
  double betaZero;
 
  double KQuark;
 
  double KGluon;
  
private:
 
  // double Vj(const ParticleData&, const ParticleData&, Energy2, double, bool=false) const;
  double Vj(const ParticleData&, const ParticleData&, Energy2, double, double, int, bool=false) const;
 
  double VsDoublePole(const ParticleData&, const ParticleData&) const;
 
  double VsSinglePole(const ParticleData&, const ParticleData&, Energy2) const;
 
  double GammaQuark(const ParticleData&) const;
  
  double GammaGluon() const;
 
  double GammaQuarkSinglePole(const ParticleData&) const;
  
  double GammaGluonSinglePole() const;
 
private:
 
  DipoleMIOperator & operator=(const DipoleMIOperator &) = delete;
 
};
 
}
 
#endif /* HERWIG_DipoleMIOperator_H */