HalfHalfZeroEWSplitFn.h

// -*- C++ -*-
#ifndef Herwig_HalfHalfZeroEWSplitFn_H
#define Herwig_HalfHalfZeroEWSplitFn_H
//
// This is the declaration of the HalfHalfZeroEWSplitFn class.
//
 
#include "SplittingFunction.h"
#include "Herwig/Models/StandardModel/StandardModel.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class HalfHalfZeroEWSplitFn: 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 & gH, const IdList & ids) const;
 
  void getCouplings(double & gH, const IdList & ids, const Energy2 t) 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:
 
  HalfHalfZeroEWSplitFn & operator=(const HalfHalfZeroEWSplitFn &) = delete;
 
private:
 
  double ghqq_;
 
 
  tcHwSMPtr _theSM;
 
};
 
}
 
#endif /* Herwig_HalfHalfZeroEWSplitFn_H */