Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Gtobssim in fermitools-2.0.0 cannot handle small (float64) values in source maps #3

Open
Alex-Broughton opened this issue Dec 2, 2020 · 1 comment
Assignees

Comments

@Alex-Broughton
Copy link

Alex-Broughton commented Dec 2, 2020

With fermitools-2.0.0/fermitools-data-0.18, and using python3, I attempted to run gtobssim with Mapcube sources. The largest of the values of the data in the source files where of order 1e-7). In the file they are stored as float64 numbers, but I believe gtobssim treats them as floats if I am not mistaken. The result is that it treats these values as zero, which throws the error below.

`$ gtobssim emin=1.e3 emax=5.e5 startdate="2008-01-01 00:00:00" infile="srcmdl_gtobssim.xml" srclist="srclist.txt" scfile=none simtime=604800 use_ac=no irfs="P8R3_SOURCE_V2" evroot=root seed=293060
Simulation start time (seconds wrt MET 0)[INDEF] INDEF

Caught St13runtime_error at the top level: MapCube: There are negative or zero-valued pixels in the FITS image read from/data/users/abrought/ML_Project/Galactic_Center_Machine_Learning/Gamma_Ray_Sims/Process_GALPROP_Output/zeroed/zeroed_masked_reprojected_pion_decay_annuli_0_5.fits`

@dagorym
Copy link
Member

dagorym commented Dec 8, 2020

@Alex-Broughton could you attach the input files you are using so we can try to reproduce this?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants