ExSampler2.h

// -*- C++ -*-
//
// ExSampler.h is a part of Herwig++ - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2012 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_ExSampler_H
#define Herwig_ExSampler_H
//
// This is the declaration of the ExSampler class.
//

#include "BinSampler.h"
#include "Herwig++/Exsample2/exsample/generator.h"

namespace Herwig {

using namespace ThePEG;

  class ExSampler: public Herwig::BinSampler {

public:

  ExSampler();

  virtual ~ExSampler();

public:

  virtual bool isUnweighting() const { return true; }

  virtual bool compensating() const { return generator_.compensating(); }

  virtual void generate(bool noMaxInfo = false);

  virtual void initialize(bool progress);

  virtual void finalize(bool verbose);

public:

  virtual void reject() {
    GeneralStatistics::reject();
    generator_.reject();
  }

public:

  double evaluate(const vector<double>& p) const {
    double ret;
    try {
      ret = eventHandler()->dSigDR(p) / nanobarn;
    } catch (Veto&) {
      ret = 0.0;
    } catch (...) {
      throw;
    }
    return ret;
  }

  pair<vector<double>,vector<double> > support() const {
    vector<double> lower(dimension(),0.);
    vector<double> upper(dimension(),1.);
    return make_pair(lower,upper);
  }

  void start_presampling() { }

  void stop_presampling() { }

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:

  unsigned long presampling_points_;

  unsigned long freeze_grid_;

  double efficiency_threshold_;

  double gain_threshold_;

  exsample::generator<ExSampler,UseRandom> generator_;

private:

  ExSampler & operator=(const ExSampler &);

};

}

#endif /* Herwig_ExSampler_H */