VariableMassCutOff.h

// -*- C++ -*-
#ifndef Herwig_VariableMassCutOff_H
#define Herwig_VariableMassCutOff_H
//
// This is the declaration of the VariableMassCutOff class.
//
 
#include "SudakovCutOff.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class VariableMassCutOff: public SudakovCutOff {
 
public:
 
  virtual const vector<Energy> & virtualMasses(const IdList & ids);
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
private:
 
  Energy kinematicCutOff(Energy scale, Energy mq) const {
        return max((scale -a_*mq)/b_,c_);
  }
 
 
  VariableMassCutOff & operator=(const VariableMassCutOff &) = delete;
 
private:
 
  double a_ = 0.3;
 
  double b_ = 2.3;
 
  Energy c_ = 0.3_GeV;
 
  Energy kinCutoffScale_ = 2.3_GeV;
 
 
 
};
 
}
 
#endif /* Herwig_VariableMassCutOff_H */