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