def test_width_0():
geom, image, clean_mask = create_sample_image("30d")
hillas = hillas_parameters(geom[clean_mask], image[clean_mask])
hillas.width = 0 * u.m
conc = concentration_parameters(geom, image, hillas)
assert conc.core == 0
def test_concentration():
geom, image, clean_mask = create_sample_image("30d")
hillas = hillas_parameters(geom[clean_mask], image[clean_mask])
conc = concentration_parameters(geom, image, hillas)
assert 0.1 <= conc.cog <= 0.3
assert 0.04 <= conc.pixel <= 0.2
assert 0.3 <= conc.core <= 0.6
def test_no_pixels_near_cog():
geom, image, clean_mask = create_sample_image("30d")
hillas = hillas_parameters(geom[clean_mask], image[clean_mask])
# remove pixels close to cog from the cleaning pixels
clean_mask &= ((geom.pix_x - hillas.x)**2 +
(geom.pix_y - hillas.y)**2) >= (2 * geom.pixel_width**2)
conc = concentration_parameters(geom[clean_mask], image[clean_mask],
hillas)
assert conc.cog == 0