DressedLepton.hh

// -*- C++ -*-
#ifndef RIVET_DressedLepton_HH
#define RIVET_DressedLepton_HH
 
#include "Rivet/Particle.hh"
#include "Rivet/Config/RivetCommon.hh"
 
namespace Rivet {
 
 
  enum class LeptonOrigin { PROMPT=1, NODECAY=1, ALL };
 
  enum class LeptonReco { ALL=0, ALL_BARE=0,
              ALL_DRESSED=1,
              DIRECT_BARE=2, PROMPT_BARE=2,
              DIRECT_DRESSED=3, PROMPT_DRESSED=3 };
 
  enum class DressingType { DR=0, CONE=0, CLUSTER=1, AKT=1 };
 
  enum class PhotonOrigin { NONE=0, PROMPT=1, NODECAY=1, ALL };
 
 
  class DressedLepton : public Particle {
  public:
 
    DressedLepton(const Particle& dlepton);
 
    DressedLepton(const Particle& lepton, const Particles& photons, bool momsum=true);
 
    void addPhoton(const Particle& p, bool momsum=true);
 
    const Particle& bareLepton() const;
 
    const Particles photons() const { return slice(constituents(), 1); }
 
  };
 
 
  using DressedLeptons = vector<DressedLepton>;
 
  template <typename T, typename FN = T(const ParticleBase&)>
  inline T sum(const DressedLeptons& c, FN&& fn, const T& start=T()) {
    auto f = std::function(std::forward<FN>(fn));
    T rtn = start;
    for (const auto& x : c) rtn += fn(x);
    return rtn;
  }
 
 
 
}
 
#endif