ME2to2QCD.h

// -*- C++ -*-
//
// ME2to2QCD.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2019 Leif Lonnblad
//
// ThePEG is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef ThePEG_ME2to2QCD_H
#define ThePEG_ME2to2QCD_H
// This is the declaration of the ME2to2QCD class.

#include "ThePEG/MatrixElement/ME2to2Base.h"

namespace ThePEG {

class ME2to2QCD: public ME2to2Base {

public:

  ME2to2QCD()
    : theMaxFlavour(5), theKfac(1.0), theKfacA(1.0), useInterference(true) {}

  virtual ~ME2to2QCD();

public:

  virtual unsigned int orderInAlphaS() const;

  virtual unsigned int orderInAlphaEW() const;

  double comfac() const;

  int maxFlavour() const { return theMaxFlavour; }

  double Kfac() const { return theKfac; }

  double KfacA() const { return theKfacA >= 0.0? theKfacA: theKfac; }

  bool interference() const { return useInterference; }

  bool isQuark(const ParticleData & p) const {
    return ( p.id() && abs(p.id()) <= maxFlavour() );
  }

  tcPDPtr quark(int i) const;

public:

  void persistentOutput(PersistentOStream & os) const;

  void persistentInput(PersistentIStream & is, int version);

  static void Init();

private:

  int theMaxFlavour;

  double theKfac;

  double theKfacA;

  bool useInterference;

private:

  static AbstractClassDescription<ME2to2QCD> initME2to2QCD;

  ME2to2QCD & operator=(const ME2to2QCD &) = delete;

};

}


namespace ThePEG {

template <>
struct BaseClassTrait<ME2to2QCD,1>: public ClassTraitsType {
  typedef ME2to2Base NthBase;
};

template <>
struct ClassTraits<ME2to2QCD>: public ClassTraitsBase<ME2to2QCD> {
  static string className() { return "ThePEG::ME2to2QCD"; }
};

}

#endif /* ThePEG_ME2to2QCD_H */