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 */