SISCone 3.0.5
|
class for the computation of jet areas. More...
#include <area.h>
Public Member Functions | |
Carea () | |
default ctor | |
~Carea () | |
default dtor | |
int | compute_areas (std::vector< Cmomentum > &_particles, double _radius, double _f, int _n_pass_max=0, Esplit_merge_scale _split_merge_scale=SM_pttilde, bool _hard_only=false) |
compute the jet areas from a given particle set. | |
int | compute_active_areas (std::vector< Cmomentum > &_particles, double _radius, double _f, int _n_pass_max=0, Esplit_merge_scale _split_merge_scale=SM_pttilde, bool _hard_only=false) |
compute the jet active areas from a given particle set. | |
int | compute_passive_areas (std::vector< Cmomentum > &_particles, double _radius, double _f, int _n_pass_max=0, Esplit_merge_scale _split_merge_scale=SM_pttilde) |
compute the jet passive areas from a given particle set. | |
![]() | |
Csiscone () | |
default ctor | |
~Csiscone () | |
default dtor | |
int | compute_jets (std::vector< Cmomentum > &_particles, double _radius, double _f, int _n_pass_max=0, double _ptmin=0.0, Esplit_merge_scale _split_merge_scale=SM_pttilde) |
compute the jets from a given particle set. | |
int | compute_jets_progressive_removal (std::vector< Cmomentum > &_particles, double _radius, int _n_pass_max=0, double _ptmin=0.0, Esplit_merge_scale _ordering_scale=SM_pttilde) |
compute the jets from a given particle set. | |
int | recompute_jets (double _f, double _ptmin=0.0, Esplit_merge_scale _split_merge_scale=SM_pttilde) |
recompute the jets with a different overlap parameter. | |
![]() | |
Cstable_cones () | |
default ctor | |
Cstable_cones (std::vector< Cmomentum > &_particle_list) | |
ctor with initialisation (sse init for details) | |
~Cstable_cones () | |
default dtor | |
void | init (std::vector< Cmomentum > &_particle_list) |
initialisation | |
int | get_stable_cones (double _radius) |
compute stable cones. | |
![]() | |
Cvicinity () | |
default constructor | |
Cvicinity (std::vector< Cmomentum > &_particle_list) | |
constructor with initialisation (see set_particle_list) | |
~Cvicinity () | |
default destructor | |
void | set_particle_list (std::vector< Cmomentum > &_particle_list) |
set the particle_list | |
void | build (Cmomentum *_parent, double _VR) |
build the vicinity list from the list of points. | |
![]() | |
Csplit_merge () | |
default ctor | |
~Csplit_merge () | |
default dtor | |
int | init (std::vector< Cmomentum > &_particles, std::vector< Cmomentum > *protocones, double R2, double ptmin=0.0) |
initialisation function | |
int | init_particles (std::vector< Cmomentum > &_particles) |
initialisation function for particle list | |
int | init_pleft () |
build initial list of left particles | |
int | set_pt_weighted_splitting (bool _use_pt_weighted_splitting) |
use a pt-dependent boundary for splitting When called with true, the criterium for splitting two protojets will be to compare D1^2/kt1^2 vs. | |
int | partial_clear () |
partial clearance | |
int | full_clear () |
full clearance | |
void | set_user_scale (const Cuser_scale_base *user_scale_in) |
associate a user-defined scale to order the stable cones | |
const Cuser_scale_base * | user_scale () const |
return the user-defined scale (NULL if none) | |
int | merge_collinear_and_remove_soft () |
build the list 'p_uncol_hard' from p_remain by clustering collinear particles and removing particles softer than stable_cone_soft_pt2_cutoff note that thins in only used for stable-cone detection so the parent_index field is unnecessary | |
int | add_protocones (std::vector< Cmomentum > *protocones, double R2, double ptmin=0.0) |
add a list of protocones | |
int | add_hardest_protocone_to_jets (std::vector< Cmomentum > *protocones, double R2, double ptmin=0.0) |
remove the hardest protocone and declare it a jet | |
int | perform (double overlap_tshold, double ptmin=0.0) |
really do the splitting and merging At the end, the vector jets is filled with the jets found. | |
int | save_contents (FILE *flux) |
save final jets | |
int | show () |
show jets/candidates status | |
Public Attributes | |
int | grid_size |
size of the grid we add soft particles on (N_soft=(grid_size^2)) | |
double | grid_eta_max |
maximal value of eta we add soft particles on | |
double | grid_shift |
fractional (random) displacement of the points om the grid | |
double | pt_soft |
pt of the soft particles added | |
double | pt_shift |
amplitude of the pt random shift | |
double | pt_soft_min |
pt_min used in SM to compute passive areas | |
std::vector< Cjet_area > | jet_areas |
jets with their areas | |
![]() | |
std::vector< std::vector< Cmomentum > > | protocones_list |
list of protocones found pass-by-pass (not filled by compute_jets_progressive_removal) | |
![]() | |
std::vector< Cmomentum > | protocones |
list of stable cones | |
hash_cones * | hc |
list of candidates | |
int | nb_tot |
total number of tested cones | |
![]() | |
Cmomentum * | parent |
parent vector | |
double | VR |
radius of the vicinity | |
double | VR2 |
squared radius of the vicinity | |
double | R |
normal radius | |
double | R2 |
squared normal radius | |
double | inv_R_EPS_COCIRC |
R / EPSILON_COCIRCULAR. | |
double | inv_R_2EPS_COCIRC |
R / (2*EPSILON_COCIRCULAR) | |
int | n_part |
number of particles | |
std::vector< Cmomentum > | plist |
the list of particles | |
std::vector< Cvicinity_inclusion > | pincluded |
the inclusion state of particles | |
Cvicinity_elm * | ve_list |
list of vicinity elements built from particle list (size=2*n) | |
std::vector< Cvicinity_elm * > | vicinity |
list of points in parent's vicinity | |
unsigned int | vicinity_size |
number of elements in vicinity | |
![]() | |
int | n |
number of particles | |
std::vector< Cmomentum > | particles |
list of particles | |
std::vector< double > | pt |
list of particles' pt | |
int | n_left |
numer of particles that does not belong to any jet | |
std::vector< Cmomentum > | p_remain |
list of particles remaining to deal with | |
std::vector< Cmomentum > | p_uncol_hard |
list of particles remaining with collinear clustering | |
int | n_pass |
index of the run | |
double | most_ambiguous_split |
minimal difference in squared distance between a particle and two overlapping protojets when doing a split (useful when testing approx. | |
std::vector< Cjet > | jets |
list of jets | |
int * | indices |
maximal size array for indices works | |
int | idx_size |
number of elements in indices1 | |
bool | merge_identical_protocones |
The following flag indicates that identical protocones are to be merged automatically each time around the split-merge loop and before anything else happens. | |
Csplit_merge_ptcomparison | ptcomparison |
member used for detailed comparisons of pt's | |
double | SM_var2_hardest_cut_off |
stop split–merge or progressive-removal when the squared SM_var of the hardest protojet is below this cut-off. | |
double | stable_cone_soft_pt2_cutoff |
pt cutoff for the particles to put in p_uncol_hard this is meant to allow removing soft particles in the stable-cone search. | |
Additional Inherited Members | |
![]() | |
static void | set_banner_stream (std::ostream *ostr) |
A call to this function modifies the stream used to print banners (by default cout). | |
static std::ostream * | banner_stream () |
returns a pointer to the stream to be used to print banners (cout by default) | |
![]() | |
static bool | init_done =false |
check random generator initialisation | |
![]() | |
void | append_to_vicinity (Cmomentum *v) |
append a particle to the 'vicinity' list after having tested it and computed the angular-ordering quantities | |
![]() | |
double | R |
cone radius | |
double | R2 |
cone radius SQUARED | |
![]() | |
double | pcx |
parent centre (eta) | |
double | pcy |
parent centre (phi) | |
class for the computation of jet areas.
This is the class user should use whenever you want to compute the jet area (passive and active). . It uses the SISCone algorithm to perform the jet analysis.
siscone::Carea::Carea | ( | ) |
default ctor
Definition at line 88 of file area.cpp.
References grid_eta_max, grid_shift, grid_size, pt_shift, pt_soft, and pt_soft_min.
int siscone::Carea::compute_active_areas | ( | std::vector< Cmomentum > & | _particles, |
double | _radius, | ||
double | _f, | ||
int | _n_pass_max = 0 , |
||
Esplit_merge_scale | _split_merge_scale = SM_pttilde , |
||
bool | _hard_only = false |
||
) |
compute the jet active areas from a given particle set.
The parameters of this method are the ones which control the jet clustering alghorithn. Note that the pt_min is not allowed here soince the jet-area determination involves soft particles/jets and thus is used internally. See compute_areas for paramters definition.
Definition at line 291 of file area.cpp.
References siscone::Csiscone::compute_jets(), grid_eta_max, grid_shift, grid_size, jet_areas, siscone::Csplit_merge::jets, siscone::Csplit_merge::n, pt_shift, pt_soft, pt_soft_min, siscone::Csplit_merge::SM_var2_hardest_cut_off, and twopi.
int siscone::Carea::compute_areas | ( | std::vector< Cmomentum > & | _particles, |
double | _radius, | ||
double | _f, | ||
int | _n_pass_max = 0 , |
||
Esplit_merge_scale | _split_merge_scale = SM_pttilde , |
||
bool | _hard_only = false |
||
) |
compute the jet areas from a given particle set.
The parameters of this method are the ones which control the jet clustering alghorithn. Note that the pt_min is not allowed here soince the jet-area determination involves soft particles/jets and thus is used internally.
_particles | list of particles |
_radius | cone radius |
_f | shared energy threshold for splitting&merging |
_n_pass_max | maximum number of passes (0=full search, the default) |
_split_merge_scale | the scale choice for the split-merge procedure NOTE: SM_pt leads to IR unsafety for some events with momentum conservation. SM_Et is IR safe but not boost invariant and not implemented(!) SM_mt is IR safe for hadronic events, but not for decays of two back-to-back particles of identical mass SM_pttilde is always IR safe, and also boost invariant (default) |
_hard_only | when this is set on, only hard jets are computed and not the purely ghosted jets (default: false) |
Definition at line 125 of file area.cpp.
References siscone::Csiscone::compute_jets(), grid_eta_max, grid_shift, grid_size, jet_areas, siscone::Csplit_merge::jets, siscone::Csplit_merge::n, pt_shift, pt_soft, pt_soft_min, siscone::Csiscone::recompute_jets(), siscone::Csplit_merge::SM_var2_hardest_cut_off, and twopi.
int siscone::Carea::compute_passive_areas | ( | std::vector< Cmomentum > & | _particles, |
double | _radius, | ||
double | _f, | ||
int | _n_pass_max = 0 , |
||
Esplit_merge_scale | _split_merge_scale = SM_pttilde |
||
) |
compute the jet passive areas from a given particle set.
The parameters of this method are the ones which control the jet clustering alghorithn. Note that the pt_min is not allowed here soince the jet-area determination involves soft particles/jets and thus is used internally. See compute_areas for paramters definition.
Definition at line 217 of file area.cpp.
References siscone::Csiscone::compute_jets(), grid_eta_max, grid_shift, grid_size, jet_areas, siscone::Csplit_merge::jets, siscone::Csplit_merge::n, pt_shift, pt_soft, pt_soft_min, siscone::Csplit_merge::stable_cone_soft_pt2_cutoff, and twopi.
double siscone::Carea::grid_eta_max |
maximal value of eta we add soft particles on
Definition at line 127 of file area.h.
Referenced by Carea(), compute_active_areas(), compute_areas(), and compute_passive_areas().
double siscone::Carea::grid_shift |
fractional (random) displacement of the points om the grid
Definition at line 128 of file area.h.
Referenced by Carea(), compute_active_areas(), compute_areas(), and compute_passive_areas().
int siscone::Carea::grid_size |
size of the grid we add soft particles on (N_soft=(grid_size^2))
Definition at line 126 of file area.h.
Referenced by Carea(), compute_active_areas(), compute_areas(), and compute_passive_areas().
std::vector<Cjet_area> siscone::Carea::jet_areas |
jets with their areas
Definition at line 135 of file area.h.
Referenced by compute_active_areas(), compute_areas(), and compute_passive_areas().
double siscone::Carea::pt_shift |
amplitude of the pt random shift
Definition at line 131 of file area.h.
Referenced by Carea(), compute_active_areas(), compute_areas(), and compute_passive_areas().
double siscone::Carea::pt_soft |
pt of the soft particles added
Definition at line 130 of file area.h.
Referenced by Carea(), compute_active_areas(), compute_areas(), and compute_passive_areas().
double siscone::Carea::pt_soft_min |
pt_min used in SM to compute passive areas
Definition at line 132 of file area.h.
Referenced by Carea(), compute_active_areas(), compute_areas(), and compute_passive_areas().