ProjectingSampler.h

// -*- C++ -*-
//
// ProjectingSampler.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_ProjectingSampler_H
#define Herwig_ProjectingSampler_H
//
// This is the declaration of the ProjectingSampler class.
//

#include "BinSampler.h"
#include "BinnedStatistics.h"

namespace Herwig {

using namespace ThePEG;

class ProjectingSampler: public Herwig::BinSampler {

public:

  ProjectingSampler();

  virtual ~ProjectingSampler();

public:

  virtual void generate(bool noMaxInfo = false);

  virtual void initialize(bool progress);

  void adapt();

public:

  virtual void select(double weight);

  virtual void accept();

  virtual void reject();

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:

  ProjectingSampler & operator=(const ProjectingSampler &);

  bool theFirstIteration;

  unsigned long theNIterations;

  double theEnhancementFactor;

  unsigned int theNBins;

  double theEpsilon;

  unsigned long theLastNPoints;

  vector<BinnedStatistics> theProjections;

  double theLastValue;

  double theWeightThreshold;

};

}

#endif /* Herwig_ProjectingSampler_H */