29 #include "fastjet/ClusterSequencePassiveArea.hh"
30 #include "fastjet/ClusterSequenceVoronoiArea.hh"
32 FASTJET_BEGIN_NAMESPACE
41 void ClusterSequencePassiveArea::_initialise_and_run_PA (
42 const JetDefinition & jet_def_in,
43 const GhostedAreaSpec & area_spec,
44 const bool & writeout_combinations) {
48 ClusterSequenceVoronoiArea csva(_jets,jet_def_in,VoronoiAreaSpec(1.0));
51 transfer_from_sequence(csva);
53 _resize_and_zero_AA();
54 for (
unsigned i = 0; i < _history.size(); i++) {
55 int ijetp = _history[i].jetp_index;
56 if (ijetp != Invalid) {
57 _average_area[i] = csva.area(_jets[ijetp]);
58 _average_area_4vector[i] = csva.area_4vector(_jets[ijetp]);
65 JetDefinition tmp_jet_def = jet_def_in;
66 tmp_jet_def.set_jet_finder(cambridge_for_passive_algorithm);
67 tmp_jet_def.set_extra_param(sqrt(area_spec.mean_ghost_kt()));
68 _initialise_and_run_AA(tmp_jet_def, area_spec, writeout_combinations);
69 _jet_def = jet_def_in;
73 _initialise_and_run_AA(jet_def_in, area_spec, writeout_combinations);
76 jet_def_in.plugin()->supports_ghosted_passive_areas()) {
80 double ghost_sep_scale_store = jet_def_in.plugin()->ghost_separation_scale();
81 jet_def_in.plugin()->set_ghost_separation_scale(sqrt(area_spec.mean_ghost_kt()));
82 _initialise_and_run_AA(jet_def_in, area_spec, writeout_combinations);
85 jet_def_in.plugin()->set_ghost_separation_scale(ghost_sep_scale_store);
89 _initialise_and_run_1GPA(jet_def_in, area_spec, writeout_combinations);
95 double ClusterSequencePassiveArea::empty_area (
const Selector & selector)
const {
98 return ClusterSequenceAreaBase::empty_area(selector);
100 return ClusterSequence1GhostPassiveArea::empty_area(selector);
105 FASTJET_END_NAMESPACE