cosmicrays / hermes

HERMES is a publicly available computational framework for the line of sight integration over galactic radiative processes which creates sky maps in the HEALPix-compatibile format.
GNU General Public License v3.0
23 stars 9 forks source link

Issue in reproducing Figure 6 #43

Closed brijeshkanodia closed 1 year ago

brijeshkanodia commented 1 year ago

I was trying to reproduce Pion HI decay at 10 GeV, Figure 6 of your paper. Somehow I couldn't generate the same scale for the colour scheme. Overall the plot looks the same; the scaling is different for the colours, even after taking care of units and taking the log scale. I am attaching the scales I am getting for the Pion HI decay at 10 GeV Can you share the code you used to generate the plot or suggest something I might not follow?

carmeloevoli commented 1 year ago

Hi, how do you set the min and max of the log scale?

carmeloevoli commented 1 year ago

In fact the plot looks quite ok, just the colorbar is wrong. Could you print the min and max of your map?

brijeshkanodia commented 1 year ago

vmin and vmax are 45.66 and 52.69 which seems to be okay with the colourbar I am getting.

brijeshkanodia commented 1 year ago

Hi, how do you set the min and max of the log scale?

i just converted the base_units to log scale and changed the skymap computed in hermes_map_Pi0_HI to log scale

carmeloevoli commented 1 year ago

you may find my routines to obtain the plots as in the paper here: https://github.com/carmeloevoli/hermes-cmaps/tree/main/plots