nbodykit.algorithms.fibercollisions¶
Classes
|
Run an angular FOF algorithm to determine fiber collision groups from an input catalog, and then assign fibers such that the maximum amount of object receive a fiber. |
-
class
nbodykit.algorithms.fibercollisions.
FiberCollisions
(ra, dec, collision_radius=0.017222222222222226, seed=None, degrees=True, comm=None)[source]¶ Run an angular FOF algorithm to determine fiber collision groups from an input catalog, and then assign fibers such that the maximum amount of object receive a fiber.
This amounts to determining the following population of objects:
- population 1:
the maximal “clean” sample of objects in which each object is not angularly collided with any other object in this subsample
- population 2:
the potentially-collided objects; these objects are those that are fiber collided + those that have been “resolved” due to multiple coverage in tile overlap regions
Results are computed when the object is inititalized. See the documenation of
run()
for the attributes storing the results.- Parameters
ra (array_like) – the right ascension coordinate column
dec (array_like) – the declination coordinate column
collision_radius (float, optional) – the size of the angular collision radius (in degrees); default is 62 arcseconds
seed (int, optional) – the random seed to use when determining which objects get fibers
degrees (bool, optional) – set to True if the units of
ra
anddec
are degrees
References
Methods
run
(self)Run the fiber assignment algorithm.
-
run
(self)[source]¶ Run the fiber assignment algorithm. This attaches the following attribute to the object:
Note
The
labels
attribute has a 1-to-1 correspondence with the rows in the input source.-
labels
¶ a CatalogSource that has the following columns:
- Label :
the group labels for each object in the input CatalogSource; label == 0 objects are not in a group
- Collided :
a flag array specifying which objects are collided, i.e., do not receive a fiber
- NeighborID :
for those objects that are collided, this gives the (global) index of the nearest neighbor on the sky (0-indexed) in the input catalog
source
, else it is set to -1
- Type
ArrayCatalog
; size:size
-