OneOneOneMassiveSplitFn.h

// -*- C++ -*-
//
// OneOneOneSplitFn.h is a part of Herwig - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2019 The Herwig Collaboration
//
// Herwig is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef HERWIG_OneOneOneMassiveSplitFn_H
#define HERWIG_OneOneOneMassiveSplitFn_H
//
// This is the declaration of the OneOneOneMassiveSplitFn class.
//
 
#include "Herwig/Shower/QTilde/SplittingFunctions/SplittingFunction.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class OneOneOneMassiveSplitFn: 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);
 
public:
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const {return new_ptr(*this);}
 
  virtual IBPtr fullclone() const {return new_ptr(*this);}
 
private:
 
  OneOneOneMassiveSplitFn & operator=(const OneOneOneMassiveSplitFn &) = delete;
 
};
 
}
 
#endif /* HERWIG_OneOneOneMassiveSplitFn_H */