JetSplitter.hh

Go to the documentation of this file.

 
#pragma once
 
#include <vector>
 
#include "fastjet/JetDefinition.hh"
#include "fastjet/PseudoJet.hh"
 
namespace HEJ {
  struct RNG;
 
  class JetSplitter {
  public:
    struct SplitResult {
      std::vector<fastjet::PseudoJet> constituents;
      double weight;
    };
 
 
    JetSplitter(fastjet::JetDefinition jet_def, double min_pt);
 
 
    SplitResult split(
        fastjet::PseudoJet const & j2split, int ncons, RNG & ran
    ) const;
 
    static constexpr double R_FACTOR = 5./3.;
 
  private:
    SplitResult Split2(fastjet::PseudoJet const & j2split, RNG & ran) const;
 
    double min_jet_pt_;
    fastjet::JetDefinition jet_def_;
  };
} // namespace HEJ