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 { NONE = 0, DIRECT = 1, PROMPT = 1, NODECAY = 1, DECAY = 2, ALL = 3, ANY = 3 };

  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, DIRECT = 1, PROMPT = 1, NODECAY = 1, DECAY = 2, ALL = 3, ANY = 3 };
 

  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