ColourMatrixElementCorrection.h

// -*- C++ -*-
//
// ColourMatrixElementCorrection.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 2 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef Herwig_ColourMatrixElementCorrection_H
#define Herwig_ColourMatrixElementCorrection_H
//
// This is the declaration of the ColourMatrixElementCorrection class.
//
 
#include "Herwig/Shower/Dipole/Base/DipoleSplittingReweight.h"
 
#include <tuple>
 
namespace Herwig {
 
using namespace ThePEG;
 
class ColourMatrixElementCorrection: public DipoleSplittingReweight {
 
public:
 
  ColourMatrixElementCorrection();
 
public:
 
  double cmec(const DipoleSplittingInfo&) const;
 
  virtual double evaluate(const DipoleSplittingInfo& s) const {
    return cmec(s);
  }
 
  virtual double hint(const DipoleSplittingInfo& s) const {
    if ( hintOnly(s) )
      return cmec(s);
    return abs(cmec(s))*lambda;
  }
 
  virtual bool hintOnly(const DipoleSplittingInfo& s) const {
    return cmec(s) > 0.;
  }
 
  virtual void reweightFactor(const double c) {
    assert(c > 0.0);
    lambda = c;
  }
 
  virtual void negativeScaling(const double c) {
    assert(c >= 0.0);
    negCMECScaling = c;
  }
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
// If needed, insert declarations of virtual function defined in the
// InterfacedBase class here (using ThePEG-interfaced-decl in Emacs).
 
private:
 
  double lambda;
 
  double negCMECScaling;
 
  ColourMatrixElementCorrection & operator=(const ColourMatrixElementCorrection &) = delete;
 
};
 
}
 
#endif /* Herwig_ColourMatrixElementCorrection_H */