ReflectedRegionsFinder#

class gammapy.makers.ReflectedRegionsFinder(angle_increment='0.1 rad', min_distance='0 rad', min_distance_input='0.1 rad', max_region_number=10000, binsz='0.01 deg')[source]#

Bases: gammapy.makers.background.reflected.RegionsFinder

Find reflected regions.

This class is responsible for placing a reflected region for a given input region and pointing position. It converts to pixel coordinates internally assuming a tangent projection at center position.

If the center lies inside the input region, no reflected regions can be found.

If you want to make a background estimate for an IACT observation using the reflected regions method, see also ReflectedRegionsBackgroundMaker.

Parameters
angle_incrementAngle, optional

Rotation angle applied when a region falls in an excluded region.

min_distanceAngle, optional

Minimum rotation angle between two consecutive reflected regions.

min_distance_inputAngle, optional

Minimum rotation angle between the input region and the first reflected region.

max_region_numberint, optional

Maximum number of regions to use.

binszAngle

Bin size of the reference map used for region finding.

Examples

>>> from astropy.coordinates import SkyCoord, Angle
>>> from regions import CircleSkyRegion
>>> from gammapy.makers import ReflectedRegionsFinder
>>> pointing = SkyCoord(83.2, 22.7, unit='deg', frame='icrs')
>>> target_position = SkyCoord(80.2, 23.5, unit='deg', frame='icrs')
>>> theta = Angle(0.4, 'deg')
>>> on_region = CircleSkyRegion(target_position, theta)
>>> finder = ReflectedRegionsFinder(min_distance_input='1 rad')
>>> regions, wcs = finder.run(region=on_region, center=pointing)
>>> print(regions[0]) 
Region: CircleSkyRegion
center: <SkyCoord (ICRS): (ra, dec) in deg
    (83.19879005, 25.57300957)>
radius: 1438.3... arcsec

Create a new RegionFinder.

Methods Summary

run(region, center[, exclusion_mask])

Find reflected regions.

Methods Documentation

run(region, center, exclusion_mask=None)[source]#

Find reflected regions.

Parameters
regionSkyRegion

Region to rotate.

centerSkyCoord

Rotation point.

exclusion_maskWcsNDMap, optional

Exclusion mask. Regions intersecting with this mask will not be included in the returned regions. Default is None.

Returns
regionslist of SkyRegion

Reflected regions.

wcsWCS

WCS for the determined regions.