SISCone 3.0.5
|
Go to the source code of this file.
Macros | |
#define | PT_TSHOLD 1000.0 |
program name More... | |
#define | EPSILON_COLLINEAR 1e-8 |
The following parameter controls collinear safety. More... | |
#define | EPSILON_COCIRCULAR 1e-12 |
The following parameter controls cocircular situations. More... | |
#define | EPSILON_SPLITMERGE 1e-12 |
The following define enables you to allow for identical protocones to be merged automatically after each split-merge step before anything else happens. More... | |
Variables | |
const double | twopi = 6.283185307179586476925286766559005768394 |
definition of 2*M_PI which is useful a bit everyhere! More... | |
#define EPSILON_COCIRCULAR 1e-12 |
#define EPSILON_COLLINEAR 1e-8 |
The following parameter controls collinear safety.
For the set of particles used in the search of stable cones, we gather particles if their distance in eta and phi is smaller than EPSILON_COLLINEAR.
NB: for things to behave sensibly one requires 1e-15 << EPSILON_COCIRCULAR << EPSILON_COLLINEAR << 1
among the scales that appear in practice (e.g. in deciding to use special strategies), we have EPSILON_COLLINEAR, EPSILON_COCIRCULAR, sqrt(EPSILON_COCIRCULAR) and EPSILON_COLLINEAR / EPSILON_COCIRCULAR.
#define EPSILON_SPLITMERGE 1e-12 |
The following define enables you to allow for identical protocones to be merged automatically after each split-merge step before anything else happens.
Whether this happens depends on teh value of the merge_identical_protocones flag in Csplit_merge.
It we allow such a merging and define allow MERGE_IDENTICAL_PROTOCONES_DEFAULT_TRUE then the 'merge_identical_protocones' flag in Csplit_merge to be set to 'true'. It may be manually reset to false in which case the merging of identical protocones (protojets) will be turned off.
Note that this merging identical protocones makes the algorithm infrared-unsafe, so it should be disabled except for testing purposes. if EPSILON_SPLITMERGE is defined then, during the split-merge step, when two jets are found with PTs that are identical to within a relative difference of EPSILON_SPLITMERGE they are compared element-by-element to see where the differences are, and one then uses pt1^2-pt2^2 = (pt1-pt2).(pt1+pt2) as an estimator of which is harder. NB: in unfortunate cases, this can take the split merge step up to N n * ln N time, though on normal events there don't seem to have been any major problems yet.
#define PT_TSHOLD 1000.0 |
program name
perform final stability check using the quadtree With the following define enabled, the final check for stability is performed using the quadtree rather than the explicit list of particles (see Cstable_cone::proceed_with_stability()) threshold for recomoutation of the cone (see Cstable_cones::update_cone()) When traversing cone candidates along the angular ordering, the momentum of the protojet candidate is computed incrementally from the particles that enter and leave the cone. When the cumulative change in "|px|+|py|" exceeds the cone "|px|+|py|" we explicitely recompute the cone contents
const double twopi = 6.283185307179586476925286766559005768394 |
definition of 2*M_PI which is useful a bit everyhere!
Definition at line 114 of file defines.h.
Referenced by siscone::Csplit_merge::add_hardest_protocone_to_jets(), siscone::Csplit_merge::add_protocones(), siscone::Cvicinity::append_to_vicinity(), siscone::Carea::compute_active_areas(), siscone::Carea::compute_areas(), siscone::Carea::compute_passive_areas(), siscone_spherical::CSphtheta_phi_range::CSphtheta_phi_range(), siscone_spherical::CSphsplit_merge::merge_collinear_and_remove_soft(), and siscone::Csplit_merge::merge_collinear_and_remove_soft().