FDR Analysis
Data and Tools used
Data is at
$HAWCROOT/maps/maps-20150604-aligned (this path will be called $PATH). README file states:
* These maps are made with 1486 days of HAWC-250 data, r2.01.01, with bins 1 to 9 optimized on the Crab by Paco. Bin 0 is not optimized.
* The map has been aligned using the Crab location per calibration epoch based on the work by the FearTheKitten team.
* GHW are tuned to the SNR at the Crab location.
* The PSF used for smoothing corresponds to the one obtained using expectation.py with the simulation file $HAWCROOT/sim/reco/aerie_svn_25030/aerie_HAWC-250-qe0.56_10kHz/sweets/2.63_3.45e-11_1000000_20_1.0/succeeded/2.63_3.45e-11_1000000_20_1.0_combined_rec.root with an extra smearing of 0.3 degree.
The codes used are at
http://private.hawc-observatory.org/svn/hawc/sandbox/tyapici/FDR and can be checked out via
svn. There are three main scripts written by Tolga Yapici; other scripts/codes exist in AERIE.
- getHistogramData.py: Script reads the fits files and gets the histogram for galactic plane (define in ???) and off-galactic plane. User needs to define two arguments.
- --inputs INPUT_FITS_FILE
- --output FILENAME
- processHistogramData.py: Script read the file generated by getHistogramData.pyand calculates sigma values for a set of q values. User needs to define two arguments. The output filename will be automatically generated using the input filename.
- --input INPUT_HISTO_FILENAME
- --conservative [0/1]
- drawHistogram.py:
All Skymap
The maps were generated with the commands below:
plotMollweide.py $PATH/SignifMaps/signif_hawc_111_250.fits.gz --milagro -m 0 [-c G]

Following Hao's definition of
Inner Galactic Plane , the map is separated into two regions:
|
Longitutde |
Latitude |
Galactic Plane [GP] |
15-50° |
± 4° |
Off Galactic Plane [OGP] |
elsewhere |
|
Analysis
Finding the cuts
The detailed map of the GP was generated with the command below:
plotMercator.py $PATH/SignifMaps/signif_hawc_111_250.fits.gz --origin 32.5 0 35 8 -c G --milagro --dpar 1 --dmer 5 -m 0
The same function can be used with
--xyrange
option.
The result is shown in the following figure.
The significance values were then extracted from the maps and put into a histogram.

On the plots, three of the significance cuts value corresponding to different q-values are plotted. Similar histograms are done for the Off Galactic Plane (shown below).

The significance cuts values for other q-values are tabulated below:
|
conservative c(m) |
non-conservative c(m) |
q |
significance |
|
OGP |
GP |
OGP |
GP |
0.5 |
4.14 |
2.10 |
2.80 |
0.45 |
0.1 |
4.48 |
2.80 |
3.80 |
1.75 |
0.05 |
4.65 |
3.04 |
3.97 |
2.10 |
0.01 |
4.98 |
3.51 |
4.48 |
2.75 |
0.001 |
5.49 |
4.16 |
4.98 |
3.51 |
0.0001 |
5.99 |
4.69 |
5.49 |
4.10 |
1e-05 |
6.32 |
5.15 |
5.82 |
4.69 |
Hot-spots with cuts
--
TolgaYapici - 11 Jun 2015