JetSplitter.hh

Go to the documentation of this file.

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