TauFinder.hh

Go to the documentation of this file.

#ifndef RIVET_TauFinder_HH
#define RIVET_TauFinder_HH

#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/UnstableFinalState.hh"

namespace Rivet {



  /// @brief Convenience finder of unstable taus
  /// @todo Inherit directly from ParticleFinder, not FinalState
  class TauFinder : public FinalState {
  public:

    enum DecayType { ANY=0, LEPTONIC=1, HADRONIC };
    
    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(DecayType decaytype, const Cut& cut=Cuts::open()) {
      setName("TauFinder");
      //_theParticles.clear();
      _dectype=decaytype;
      //cout << "Init tf with "<< _dectype << endl;
      addProjection(UnstableFinalState(cut), "UFS");
    }


    /// Clone on the heap.
    virtual const Projection* clone() const {
      return new TauFinder(*this);
    }


    const Particles& taus() const { return _theParticles; }


  protected:

    /// Apply the projection on the supplied event.
    void project(const Event& e);

    /// Compare with other projections.
    virtual int compare(const Projection& p) const;


  private:

    /// The decaytype enum
    DecayType _dectype;
  
  };


}


#endif