InvisibleFinalState.hh

// -*- C++ -*-
#ifndef RIVET_InvisibleFinalState_HH
#define RIVET_InvisibleFinalState_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/PromptFinalState.hh"
#include "Rivet/Tools/Logging.hh"
 
namespace Rivet {
 
 
  enum class OnlyDirect { YES, NO };
  using OnlyPrompt = OnlyDirect;

  class InvisibleFinalState : public FinalState {
  public:


    InvisibleFinalState(OnlyDirect requirepromptness = OnlyDirect::NO,
                        TauDecaysAs taudecays = TauDecaysAs::NONPROMPT,
                        MuDecaysAs mudecays = MuDecaysAs::NONPROMPT)
        : _requirePromptness(requirepromptness == OnlyDirect::YES),
          _taudecays(taudecays == TauDecaysAs::PROMPT),
          _mudecays(mudecays == MuDecaysAs::PROMPT) {
      setName("InvisibleFinalState");
      declare(FinalState(), "FS");
    }

    RIVET_DEFAULT_PROJ_CLONE(InvisibleFinalState);


    using Projection::operator=;
 

    void requirePromptness(bool acc = true,
                           TauDecaysAs taudecays = TauDecaysAs::NONPROMPT,
                           MuDecaysAs mudecays = MuDecaysAs::NONPROMPT) {
      _requirePromptness = acc;
      _taudecays = (taudecays == TauDecaysAs::PROMPT);
      _mudecays = (mudecays == MuDecaysAs::PROMPT);
    }

    void project(const Event& e);

    CmpState compare(const Projection& p) const;
 
 
  protected:
 
    bool _requirePromptness, _taudecays, _mudecays;
  };
 
 
}
 
#endif