MEPP2GammaJet.h

// -*- C++ -*-
//
// MEPP2GammaJet.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_MEPP2GammaJet_H
#define HERWIG_MEPP2GammaJet_H
//
// This is the declaration of the MEPP2GammaJet 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;
using namespace ThePEG::Helicity;
 
class MEPP2GammaJet: public HwMEBase {
 
public:
 
  MEPP2GammaJet();
 
  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;
 
  virtual void constructVertex(tSubProPtr);
 
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const { return new_ptr(*this); }
 
  virtual IBPtr fullclone() const { return new_ptr(*this); }
 
protected:
 
  virtual void doinit();
 
  virtual void rebind(const TranslationMap & trans)
   ;
 
  virtual IVector getReferences();
 
private:
 
  double qqbarME(vector<SpinorWaveFunction> & fin, vector<SpinorBarWaveFunction> & ain,
                 vector<VectorWaveFunction> & gout, vector<VectorWaveFunction>    & pout,
                 bool me) const;
 
  double qgME(vector<SpinorWaveFunction> & fin,vector<VectorWaveFunction>     & gin,
              vector<VectorWaveFunction> & pout,vector<SpinorBarWaveFunction> & fout,
              bool me) const;
 
  double qbargME(vector<SpinorBarWaveFunction> & ain, vector<VectorWaveFunction> & gin,
                 vector<VectorWaveFunction> & pout, vector<SpinorWaveFunction> & aout,
                 bool me) const;
  
private:
 
  MEPP2GammaJet & operator=(const MEPP2GammaJet &) = delete;
 
private:
 
  AbstractFFVVertexPtr _gluonvertex;
 
  AbstractFFVVertexPtr _photonvertex;
 
  int _maxflavour;
 
  unsigned int _processopt;
  
  ProductionMatrixElement _me;
 
  double scalePreFactor_;
};
 
}
 
#endif /* HERWIG_MEPP2GammaJet_H */