icecube / skyreader

An API for Results Produced by SkyDriver & the Skymap Scanner
MIT License
0 stars 0 forks source link

Improved area calculation, mock ts maps for circularized contours, and neutrino floor implementation #34

Closed G-Sommani closed 1 month ago

G-Sommani commented 2 months ago

This pull request would introduce three main changes/new implementations:

  1. Improved area calculation. The estimation with Gauss-Green's theorem takes now into account the spherical space. The estimation by counting pixels is therefore now dropped.
  2. Mock ts maps for circularized contours. A new parameter, circularized_ts_map, is now introduced. As default it is False, but if set to True it produces a circularized ts map with a bi-variate Gaussian shape and centered on the best-fit direction. The standard deviation of the bi-variate Gaussian is estimated from circular_err90 such as the contour with 90% coverage is at a radius equal to circular_err90.
  3. Neutrino floor implementation. Also in this case there is a new parameter, neutrino_floor, by default set as False. If set into True it checks the area of the 90% contour. If it is smaller than what expected in case of neutrino floor (as implemented offline, a bi-variate Gaussian with sigma=0.2 degrees), it sets circularized_ts_map to True and circular_err90 to the correct value for having sigma=0.2 degrees.