MEGammaP2Jets.h

// -*- C++ -*-
#ifndef HERWIG_MEGammaP2Jets_H
#define HERWIG_MEGammaP2Jets_H
//
// This is the declaration of the MEGammaP2Jets class.
//
 
#include "Herwig/MatrixElement/HwMEBase.h"
#include "ThePEG/Helicity/Vertex/AbstractFFVVertex.h"
#include "Herwig/MatrixElement/ProductionMatrixElement.h"
#include "ThePEG/Helicity/WaveFunction/SpinorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/VectorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/SpinorBarWaveFunction.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class MEGammaP2Jets: public HwMEBase {
 
public:
 
  MEGammaP2Jets();
 
  virtual unsigned int orderInAlphaS() const;
 
  virtual unsigned int orderInAlphaEW() const;
 
  virtual double me2() const;
 
  virtual Energy2 scale() const;
 
  virtual void getDiagrams() const;
 
  virtual Selector<DiagramIndex> diagrams(const DiagramVector & dv) const;
 
  virtual Selector<const ColourLines *>
  colourGeometries(tcDiagPtr diag) const;
 
protected:
 
  double gammagluonME(vector<VectorWaveFunction> & gmin,
                      vector<VectorWaveFunction> & glin,
                      vector<SpinorBarWaveFunction> & fout, 
                      vector<SpinorWaveFunction> & aout,
                      bool calc) const;
 
  double gammaquarkME(vector<VectorWaveFunction> & gmin,
                      vector<SpinorWaveFunction> & fin,
                      vector<VectorWaveFunction> & gout,
                      vector<SpinorBarWaveFunction> & fout,
                      bool calc) const;
 
  double gammaantiquarkME(vector<VectorWaveFunction> & gmin,
                          vector<SpinorBarWaveFunction> & fin,
                          vector<VectorWaveFunction> & gout,
                          vector<SpinorWaveFunction> & fout,
                          bool calc) 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;
 
protected:
 
  virtual void doinit();
 
private:
 
  MEGammaP2Jets & operator=(const MEGammaP2Jets &) = delete;
 
private:
 
  AbstractFFVVertexPtr _gluonvertex;
 
  AbstractFFVVertexPtr _photonvertex;
 
  unsigned int _process;
 
  int _minflavour;
 
  int _maxflavour;
  
  ProductionMatrixElement _me;
 
};
 
}
 
#endif /* HERWIG_MEGammaP2Jets_H */