HalfHalfOneEWSplitFn.h

// -*- C++ -*-
#ifndef Herwig_HalfHalfOneEWSplitFn_H
#define Herwig_HalfHalfOneEWSplitFn_H
//
// This is the declaration of the HalfHalfOneEWSplitFn class.
//
 
#include "SplittingFunction.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class HalfHalfOneEWSplitFn: 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 & gL, double & gR, 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:
 
  HalfHalfOneEWSplitFn & operator=(const HalfHalfOneEWSplitFn &) = delete;
 
private:
 
  map<long,pair<double,double> > gZ_;
 
  double gWL_;
 
};
 
}
 
#endif /* Herwig_HalfHalfOneEWSplitFn_H */