ConstLossyFinalState.hh

// -*- C++ -*-
#ifndef RIVET_ConstLossyFinalState_HH
#define RIVET_ConstLossyFinalState_HH
 
#include "Rivet/Event.hh"
#include "Rivet/Particle.hh"
#include "Rivet/Projection.hh"
#include "Rivet/Config/RivetCommon.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/LossyFinalState.hh"
#include "Rivet/Tools/Logging.hh"
#include "Rivet/Tools/Random.hh"
 
namespace Rivet {
 

  class ConstRandomFilter {
  public:
 
    ConstRandomFilter(double lossFraction)
        : _lossFraction(lossFraction) {
      assert(_lossFraction >= 0);
    }
 
    // If operator() returns true, particle is deleted ("lost")
    bool operator()(const Particle&) {
      return rand01() < _lossFraction;
    }
 
    CmpState compare(const ConstRandomFilter& other) const {
      return cmp(_lossFraction, other._lossFraction);
    }
 
 
  protected:
 
    double _lossFraction;
  };
 

  class ConstLossyFinalState : public LossyFinalState<ConstRandomFilter> {
  public:


    ConstLossyFinalState(const FinalState& fsp, double lossfraction)
        : LossyFinalState<ConstRandomFilter>(fsp, ConstRandomFilter(lossfraction)) {
      setName("ConstLossyFinalState");
    }

    ConstLossyFinalState(double lossfraction, const Cut& c = Cuts::open())
        : LossyFinalState<ConstRandomFilter>(ConstRandomFilter(lossfraction), c) {
      setName("ConstLossyFinalState");
    }

    RIVET_DEFAULT_PROJ_CLONE(ConstLossyFinalState);

    using Projection::operator=;

  };
 
 
}
 
#endif