DipoleMPKOperator.h

// -*- C++ -*-
//
// DipoleMPKOperator.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_DipoleMPKOperator_H
#define HERWIG_DipoleMPKOperator_H
//
// This is the declaration of the DipoleMPKOperator class.
//
 
#include "Herwig/MatrixElement/Matchbox/InsertionOperators/MatchboxInsertionOperator.h"
#include "Herwig/MatrixElement/Matchbox/Base/MatchboxMEBase.h"
#include "ThePEG/PDF/PDF.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class DipoleMPKOperator: public MatchboxInsertionOperator {
 
public:
 
  DipoleMPKOperator();
 
  virtual ~DipoleMPKOperator();
 
public:
 
  virtual void setXComb(tStdXCombPtr xc);
 
  virtual int nDimAdditional() const { return 1; }
 
  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 CrossSection dSigHatDR() const { 
    return sqr(hbarc) * me2() * lastBorn()->lastXComb().jacobian() / (2.*lastSHat());
  }
 
public:
 
  bool apply(tcPDPtr) const;
 
  bool applyNotMassless(tcPDPtr) const;
 
  virtual bool apply(const cPDVector&) const;
 
  double gammaSoft() const;
 
  double softLogByz(tcPDPtr p) const;
 
  double softLog(tcPDPtr p) const;
 
  double KBarqq() const;
 
  double KTildeqq() const;
 
  double Pqq() const;
 
  double KBarqg() const;
 
  double KTildeqg() const;
 
  double Pqg() const;
 
  double KBargq() const;
 
  double KTildegq() const;
 
  double Pgq() const;
 
  double KBargg() const;
 
  double KTildegg() const;
 
  double Pgg() const;
 
 
  double Ja_gQplus(double muQ2) const;
 
  double gammaSoft2(double muQ2) const;
 
  double Kscriptqq_q(Energy2 sja, Energy2 mj2) const;
 
  double Kscriptqg_q(Energy2 sja, Energy2 mj2) const;
 
  double Kscriptgq_q() const;
 
  double Kscriptgg_q(Energy2 sja, Energy2 mj2) const;
 
  double Kscriptbarqq_q(Energy2 Qja2, Energy2 mj2) const;
 
  double Kscriptbarqg_q(Energy2 Qja2, Energy2 mj2) const;
 
  double Kscriptbargq_q() const;
 
  double Kscriptbargg_q(Energy2 Qja2, Energy2 mj2) const;
 
 
  double JaNS_QQ(double muQ2) const;
 
  double Ja_QQzplus(double muQ2, int F, double zplus) const;
 
  // double Kscriptqq_g(Energy2 sja) const;
  double Kscriptqq_g(Energy2 sja, double lambda) const;
 
  double Kscriptqg_g() const;
 
  double Kscriptgq_g() const;
 
  // double Kscriptgg_g(Energy2 sja) const;
  double Kscriptgg_g(Energy2 sja, double lambda) const;
 
  // double Kscriptbarqq_g(Energy2 Qja2) const;
  double Kscriptbarqq_g(Energy2 Qja2, double lambda) const;
 
  double Kscriptbarqg_g() const;
 
  double Kscriptbargq_g() const;
 
  // double Kscriptbargg_g(Energy2 Qja2) const;
  double Kscriptbargg_g(Energy2 Qja2, double lambda) const;
 
 
  double sumParton(int id) const;
 
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 KQuark;
 
  double KGluon;
 
  mutable Energy2 scale;
 
  mutable tcPDFPtr pdf;
 
  mutable tcPDPtr particle;
 
  mutable double x;
 
  mutable double z;
 
  mutable map<pair<tcPDFPtr,tcPDPtr>,vector<double> > pdfCache;
 
  mutable tcPDPtr parton;
 
  double PDFx(tcPDPtr) const;
 
  double PDFxByz(tcPDPtr) const;
 
 
  double PDFxByzplus(tcPDPtr,int,double) const;
 
 
private:
 
  DipoleMPKOperator & operator=(const DipoleMPKOperator &) = delete;
 
};
 
}
 
#endif /* HERWIG_DipoleMPKOperator_H */