29 #include "fastjet/ClusterSequence1GhostPassiveArea.hh"
31 FASTJET_BEGIN_NAMESPACE
38 void ClusterSequence1GhostPassiveArea::_initialise_and_run_1GPA (
41 const bool & writeout_combinations) {
43 bool continue_running;
44 _initialise_AA(jet_def_in, area_spec, writeout_combinations, continue_running);
45 if (continue_running) {
47 _postprocess_AA(area_spec);
54 void ClusterSequence1GhostPassiveArea::_run_1GPA (
const GhostedAreaSpec & area_spec) {
56 vector<PseudoJet> input_jets(_jets);
59 vector<int> unique_tree;
62 valarray<double> lcl_average_area2(0.0, _average_area.size());
63 valarray<double> last_average_area(0.0, _average_area.size());
66 for (
int irepeat = 0; irepeat < area_spec.repeat(); irepeat++) {
69 vector<PseudoJet> all_ghosts;
73 for (
unsigned ig = 0; ig < all_ghosts.size(); ig++) {
74 vector<PseudoJet> some_ghosts;
75 some_ghosts.push_back(all_ghosts[ig]);
77 some_ghosts, area_spec.actual_ghost_area());
79 if (irepeat == 0 && ig == 0) {
82 _transfer_ghost_free_history(clust_seq);
84 unique_tree = unique_history_order();
88 _transfer_areas(unique_tree, clust_seq);
91 lcl_average_area2 += (_average_area - last_average_area)*
92 (_average_area - last_average_area);
93 last_average_area = _average_area;
95 _average_area2 = lcl_average_area2;
Like ClusterSequence with computation of the active jet area with the addition of explicit ghosts...
Parameters to configure the computation of jet areas using ghosts.
void add_ghosts(std::vector< PseudoJet > &) const
push a set of ghost 4-momenta onto the back of the vector of PseudoJets
class that is intended to hold a full definition of the jet clusterer