OneOneOneQEDSplitFn.h

// -*- C++ -*-
#ifndef Herwig_OneOneOneQEDSplitFn_H
#define Herwig_OneOneOneQEDSplitFn_H
//
// This is the declaration of the OneOneOneQEDSplitFn class.
//
 
#include "SplittingFunction.h"
#include "Herwig/Models/StandardModel/StandardModel.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class OneOneOneQEDSplitFn: public SplittingFunction {
 
public:
 
  virtual bool accept(const IdList & ids) const;
 
  virtual double P(const double z, const Energy2 t, const IdList & ids,
                   const bool mass, const RhoDMatrix & rho) const;
 
  virtual double overestimateP(const double z, const IdList & ids) const;
 
  virtual double ratioP(const double z, const Energy2 t, const IdList & ids,
                        const bool mass, const RhoDMatrix & rho) const;
 
  virtual double integOverP(const double z, const IdList & ids,
                            unsigned int PDFfactor=0) const;
 
  virtual double invIntegOverP(const double r, const IdList & ids,
                               unsigned int PDFfactor=0) const;
 
  virtual vector<pair<int,Complex> >
  generatePhiForward(const double z, const Energy2 t, const IdList & ids,
              const RhoDMatrix &);
 
  virtual vector<pair<int,Complex> >
  generatePhiBackward(const double z, const Energy2 t, const IdList & ids,
                      const RhoDMatrix &);
 
  virtual DecayMEPtr matrixElement(const double z, const Energy2 t,
                                   const IdList & ids, const double phi, bool timeLike);
 
protected:
 
  void getCouplings(double & gvvv, const IdList & ids) 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:
 
  OneOneOneQEDSplitFn & operator=(const OneOneOneQEDSplitFn &) = delete;
 
private:
 
  double gWWG_;
 
  tcHwSMPtr _theSM;
};
 
}
 
#endif /* Herwig_OneOneOneQEDSplitFn_H */