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_increment
Angle
, optional Rotation angle applied when a region falls in an excluded region.
- min_distance
Angle
, optional Minimal distance between two consecutive reflected regions
- min_distance_input
Angle
, optional Minimal distance from input region
- max_region_numberint, optional
Maximum number of regions to use
- binsz
Angle
Bin size of the reference map used for region finding.
- angle_increment
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