MENeutralCurrentDIS.h

// -*- C++ -*-
#ifndef HERWIG_MENeutralCurrentDIS_H
#define HERWIG_MENeutralCurrentDIS_H
//
// This is the declaration of the MENeutralCurrentDIS class.
//
 
#include "DISBase.h"
#include "ThePEG/Helicity/Vertex/AbstractFFVVertex.fh"
#include "Herwig/MatrixElement/ProductionMatrixElement.h"
#include "ThePEG/Helicity/WaveFunction/SpinorWaveFunction.h"
#include "ThePEG/Helicity/WaveFunction/SpinorBarWaveFunction.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class MENeutralCurrentDIS: public DISBase {
 
public:
 
  MENeutralCurrentDIS();
 
  virtual unsigned int orderInAlphaS() const;
 
  virtual unsigned int orderInAlphaEW() const;
 
  virtual double me2() 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:
 
  double helicityME(vector<SpinorWaveFunction>    & f1 ,
                    vector<SpinorWaveFunction>    & f2,
                    vector<SpinorBarWaveFunction> & a1 ,
                    vector<SpinorBarWaveFunction> & a2,
                    bool lorder, bool qorder,
                    bool me) const;
 
 
  inline unsigned int gammaZOption() const {return _gammaZ;}
 
  virtual double A(tcPDPtr lin, tcPDPtr lout, tcPDPtr qin, tcPDPtr qout,
                   Energy2 scale) const;
 
protected:
 
  virtual void doinit();
 
protected:
 
  virtual IBPtr clone() const {return new_ptr(*this);}
 
  virtual IBPtr fullclone() const {return new_ptr(*this);}
 
private:
 
  MENeutralCurrentDIS & operator=(const MENeutralCurrentDIS &) = delete;
 
private:
 
  AbstractFFVVertexPtr _theFFZVertex;
 
  AbstractFFVVertexPtr _theFFPVertex;
 
  tcPDPtr _z0;
 
  tcPDPtr _gamma;
 
  int _minflavour;
 
  int _maxflavour;
 
  unsigned int _gammaZ;
 
  ProductionMatrixElement _me;
 
  double _sinW;
 
  double _cosW;
 
  Energy2 _mz2;
 
};
 
}
 
#endif /* HERWIG_MENeutralCurrentDIS_H */