DileptonFinder.hh

// -*- C++ -*-
#ifndef RIVET_DileptonFinder_HH
#define RIVET_DileptonFinder_HH
 
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/LeptonFinder.hh"
#include "Rivet/Projections/VetoedFinalState.hh"
 
namespace Rivet {
 

  class DileptonFinder : public FinalState {
  public:


    DileptonFinder(double masstarget,
                   double dRdress,
                   const Cut& lcuts = Cuts::OPEN,
                   const Cut& llcuts = Cuts::OPEN,
                   LeptonOrigin leptonOrigin = LeptonOrigin::PROMPT,
                   PhotonOrigin photonOrigin = PhotonOrigin::NODECAY,
                   TauDecaysAs tauDecays = TauDecaysAs::PROMPT,
                   MuDecaysAs muDecays = MuDecaysAs::PROMPT,
                   DressingType dressing = DressingType::CONE);
 

    DileptonFinder(const FinalState& allfs,
                   double masstarget,
                   double dRdress,
                   const Cut& lcuts = Cuts::OPEN,
                   const Cut& llcuts = Cuts::OPEN,
                   DressingType dressing = DressingType::CONE);
 

    DileptonFinder(const FinalState& leptonfs,
                   const FinalState& photonfs,
                   double masstarget,
                   double dRdress,
                   const Cut& lcuts = Cuts::OPEN,
                   const Cut& llcuts = Cuts::OPEN,
                   DressingType dressing = DressingType::CONE);
 

    RIVET_DEFAULT_PROJ_CLONE(DileptonFinder);

 

    const Particles& bosons() const {
      return particles();
    }

    const Particle& boson() const {
      return bosons().front();
    }
 

    const Particles& constituents() const;

    const Particles& leptons() const {
      return constituents();
    }
 

    const VetoedFinalState& remainingFinalState() const;
 
 
  protected:

    void project(const Event& e);

    CmpState compare(const Projection& p) const;
 
 
  public:

    void clear() {
      _theParticles.clear();
    }

    using Projection::operator=;
 
 
  protected:

    double _masstarget;
  };
 

  using ZFinder = DileptonFinder;
  using LLFinder = DileptonFinder;
 
 
}
 
#endif