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 dRmax,
                   const Cut& lcuts=Cuts::OPEN,
                   const Cut& llcuts=Cuts::OPEN,
                   DressingType dressing=DressingType::CONE);
 
 
    DileptonFinder(const FinalState& allfs,
           const Cut& lcuts,
           double masstarget,
                   PdgId pid,
           double minmass, double maxmass,
           double dRdress)
      : DileptonFinder(allfs, masstarget, dRdress,
                       lcuts && Cuts::abspid == pid,
                       Cuts::mass >= minmass && Cuts::mass < maxmass)
    {   }
 
 
    DileptonFinder(const FinalState& allfs,
           PdgId pid,
           double masstarget,
           double dRdress,
           const Cut& lcuts=Cuts::OPEN,
           double minmass=-DBL_MAX, double maxmass=DBL_MAX)
      : DileptonFinder(allfs, masstarget, dRdress,
                       lcuts && Cuts::abspid == pid,
                       Cuts::mass >= minmass && Cuts::mass < maxmass)
    {   }
 
 
    DileptonFinder(PdgId pid,
           double masstarget,
           double dRdress,
           const Cut& lcuts=Cuts::OPEN,
           double minmass=-DBL_MAX, double maxmass=DBL_MAX,
           LeptonOrigin whichleptons=LeptonOrigin::PROMPT,
           PhotonOrigin whichphotons=PhotonOrigin::NODECAY,
                   TauDecaysAs tauDecays=TauDecaysAs::PROMPT,
                   MuDecaysAs muDecays=MuDecaysAs::PROMPT,
                   DressingType dressing=DressingType::CONE)
      : DileptonFinder(masstarget, dRdress,
                       lcuts && Cuts::abspid == pid,
                       Cuts::mass >= minmass && Cuts::mass < maxmass,
                       whichleptons, whichphotons, tauDecays, muDecays, dressing)
    {   }
 
 
    DileptonFinder(double masstarget,
           double dRdress,
           const Cut& lcuts=Cuts::OPEN,
           double minmass=-DBL_MAX, double maxmass=DBL_MAX,
           LeptonOrigin whichleptons=LeptonOrigin::PROMPT,
           PhotonOrigin whichphotons=PhotonOrigin::NODECAY,
                   TauDecaysAs tauDecays=TauDecaysAs::PROMPT,
                   MuDecaysAs muDecays=MuDecaysAs::PROMPT,
                   DressingType dressing=DressingType::CONE)
      : DileptonFinder(masstarget, dRdress, lcuts,
                       Cuts::mass >= minmass && Cuts::mass < maxmass,
                       whichleptons, whichphotons, tauDecays, muDecays, dressing)
    {   }
 
 
    DileptonFinder(const FinalState& allfs,
                   double masstarget,
           double dRdress,
           const Cut& lcuts=Cuts::OPEN,
           double minmass=-DBL_MAX, double maxmass=DBL_MAX)
      : DileptonFinder(allfs, masstarget, dRdress, lcuts,
                       Cuts::mass >= minmass && Cuts::mass < maxmass)
    {   }
 
 
    RIVET_DEFAULT_PROJ_CLONE(DileptonFinder);
 
 
 
    const Particles& bosons() const { return particles(); }
    const Particle& boson() const { return bosons().front(); }
 
 
    const Particles& constituents() 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