CMWHalfHalfOneSplitFn.h

  // -*- C++ -*-
  //
  // CMWHalfHalfOneSplitFn.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_CMWHalfHalfOneSplitFn_H
#define HERWIG_CMWHalfHalfOneSplitFn_H
  //
  // This is the declaration of the CMWHalfHalfOneSplitFn class.
  //
 
#include "HalfHalfOneSplitFn.h"
#include "Herwig/Shower/ShowerAlpha.h"
 
 
namespace Herwig {
  
  using namespace ThePEG;
  
  class CMWHalfHalfOneSplitFn: public HalfHalfOneSplitFn {
    
  public:
    
    virtual double P(const double z, const Energy2 t, const IdList & ids,
                     const bool mass, const RhoDMatrix & rho) const;
    
    virtual double ratioP(const double z, const Energy2 t, const IdList & ids,
                          const bool mass, const RhoDMatrix & rho) const;
    
    double Kg(Energy2 )const{
        //TODO: Should be t-dependent
      int Nf=5;//alpha_->Nf(t)
      return (3.*(67./18.-1./6.*sqr(Constants::pi))-5./9.*Nf);
    }
    
    
    
  public:
    
    void persistentOutput(PersistentOStream & os) const;
    
    void persistentInput(PersistentIStream & is, int version);
    
    static void Init();
    
  protected:
    
    virtual void doinit(){
      HalfHalfOneSplitFn::doinit();
    };
    
    
  private:
    
      // Pointer to the alpha_s object in use.
    ShowerAlphaPtr alpha_;
      // Provide information if the kernel is used for initial state.
      // as the pt definition contains an additional factor of z.
    bool isIS_=false;
    
  protected:
    
    virtual IBPtr clone() const {return new_ptr(*this);}
    
    virtual IBPtr fullclone() const {return new_ptr(*this);}
    
  private:
    
    CMWHalfHalfOneSplitFn & operator=(const CMWHalfHalfOneSplitFn &) = delete;
    
    
  };
  
}
 
#endif /* HERWIG_CMWHalfHalfOneSplitFn_H */