ReweightEW.h

// -*- C++ -*-
#ifndef Herwig_ReweightEW_H
#define Herwig_ReweightEW_H
//
// This is the declaration of the ReweightEW class.
//
 
#include "ThePEG/MatrixElement/ReweightBase.h"
#include <array>
 
namespace Herwig {
 
using namespace ThePEG;
 
class ReweightEW: public ReweightBase {
 
public:
 
  ReweightEW();
 
  virtual ~ReweightEW();
 
public:
 
  virtual double weight() const;
 
  double lastS() const {return thelasts;} 
  double lastT() const {return thelastt;} 
  double lastK() const {return thelastk;} 
 
  void setSTK(double s, double t, double K); 
 
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
 
private:
 
  mutable double thelasts;
 
  mutable double thelastt;
 
  mutable double thelastk;
 
  std::array<std::array<double,6>,40001> tab;
 
  string filename;
 
public:
  double EWKFac(unsigned int f, double s, double t) const;
 
private:
  void inittable();
 
private:
 
  virtual void doinit();
 
  virtual void doinitrun();
 
private:
 
  ReweightEW & operator=(const ReweightEW &) = delete;
 
};
 
}
 
#endif /* Herwig_ReweightEW_H */