This little script is used to compute the mean|median energy of events based on how they have been spatially grouped via the input map file.
So the basic idea is that a user creates a map (2D image) whose pixels
values identify spatial groups. This could be the output of say
dmnautilus
or some other binning algorithm (eg contour_bin
).
The script identifies the events in each group, computes the specified statistic for the specified column, and then outputs an image with the stat value replacing the mapID value, eg a median_energy map.
statmap acis_evt.fits[ccd_id=7,energy=500:2000] cbin.map mean_energy.map col=energy stat=mean
This script requires CIAO be installed.