TauFinder.hh

#ifndef RIVET_TauFinder_HH
#define RIVET_TauFinder_HH
 
#include "Rivet/DressedLepton.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/UnstableParticles.hh"
 
namespace Rivet {
 

  enum class TauDecay { ANY = 0, ALL = 0, LEPTONIC, HADRONIC };
 

  class TauFinder : public FinalState {
  public:
 
    static bool isHadronic(const Particle& tau) {
      assert(tau.abspid() == PID::TAU);
      return any(tau.stableDescendants(), isHadron);
    }
 
    static bool isLeptonic(const Particle& tau) {
      return !isHadronic(tau);
    }
 
 
    TauFinder(TauDecay decaymode = TauDecay::ANY,
              LeptonOrigin origin = LeptonOrigin::ANY,
              const Cut& cut = Cuts::open()) {
      setName("TauFinder");
      _decmode = decaymode;
      _origin = origin;
      declare(UnstableParticles(cut), "UFS");
    }
 
 
    TauFinder(TauDecay decaymode, const Cut& cut, LeptonOrigin origin = LeptonOrigin::ANY)
        : TauFinder(decaymode, origin, cut) { }
 
 
    TauFinder(const Cut& cut, TauDecay decaymode = TauDecay::ANY, LeptonOrigin origin = LeptonOrigin::ANY)
        : TauFinder(decaymode, origin, cut) { }
 

    RIVET_DEFAULT_PROJ_CLONE(TauFinder);

    using Projection::operator=;
 
 
    const Particles& taus() const {
      return _theParticles;
    }
 
 
  protected:

    void project(const Event& e);

    virtual CmpState compare(const Projection& p) const;
 
 
  protected:

    TauDecay _decmode;

    LeptonOrigin _origin;
  };
 
 
  // Alias
  using Taus = TauFinder;
 
 
}
 
 
#endif