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;
void add_ghosts(std::vector< PseudoJet > &) const
push a set of ghost 4-momenta onto the back of the vector of PseudoJets
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.
class that is intended to hold a full definition of the jet clusterer