Galactic Plane Survey simulation's bug with the new ctools-0.10.0
|Assigned To:||Knödlseder Jürgen||% Done:|
Hi, I’m Tan. I had a bug when executing the Galactic Plane Survey script you gave me (as follows) with the new ctools-0.10.0.
It seems that the GPS is not suitable for the new version of ctools-0.10.0, because of the changes in ctmodel (I guess). Please see the attachment for details. So can you fix the GPS script ?
Thank you very much.
#3 Updated by Knödlseder Jürgen over 6 years ago
- File make_gps.py added
You are right. The
get_obs() method has been removed by accident. I did now put this method back, and everything seems to work on my side. See http://cta.irap.omp.eu/ctools/download.html to learn how to get the latest version from gitlab.
Please note that in the meanwhile, the
GSkymap class was replaced by
GSkyMap, hence the script needs to be modified. I attached the actual version of the script that works
#5 Updated by Dang Viet Tan over 6 years ago
- File gps_model_ics_1-10TeV.fits added
- File gps_model_pi0_1-10TeV.fits added
- File gps_model_1-10TeV.fits added
- File gps_model_bkg_1-10TeV.fits added
- File gps_model_cygnus_1-10TeV.fits added
- File gps_model_pwne_1-10TeV.fits added
- File gps_model_snrs_1-10TeV.fits added
- File gps_model_tevcat_1-10TeV.fits added
Thank you for fixing it. I’ve gotten the latest version from gitlab and already installed it successfully.
I also checked your fixed file - make_gps.py but there’s only 1 line in it, and it didn’t work. So I fixed my file by replacing
GSkyMap (only 2 places in the script). The script has worked, but the output files seem to be wrong. These files of
gps_model_*.fits have just only 1 pixel (see attachment). Can you check this bug?
#9 Updated by Knödlseder Jürgen over 6 years ago
- File make_gps.py added
There was another regression in addition to the removal of the
get_obs() method (see #1609), sorry for that. It should now be fixed (just get the latest code from gitlab using the
I also updated the script that works on my side.
Can you check that everything is also fine on your side?
#10 Updated by Dang Viet Tan over 6 years ago
- File ics.jpg added
- File pi0.jpg added
- File map_ics.fits added
- File map_pi0.fits added
- File model_ics.xml added
- File model_pi0.xml added
I’ve already checked it on my side and the results seems good. Thank you very much.
By the way, I still have problem with simulating ics and pi0 map by ctobssim. It didn’t work even with the new version of ctools. When I check your map (in the directory “models” in gps.tar.gz), I found that the maps inconsistent with Galactic coordinates (see attachment). Maybe that’s why I cannot execute ctobssim by using your XML and FITS files (attachment). Can you check this problem for me ?
#12 Updated by Knödlseder Jürgen over 6 years ago
Dang Viet Tan wrote:
I’m sorry that I was wrong about the coordinates mentioned above. The maps are fine.
But I still cannot execute GPS by ctobssim (it was aborted immediately). Can you check it one more time and show me the output file of ics and pi0? Thank you.
Could you provide me with the command sequence that you try to execute so that I can reproduce it here and follow this up?
#13 Updated by Dang Viet Tan over 6 years ago
- File ctobssim.py added
Yes, you can see the executed file as attachment (ctobssim.py). For easy, I just specify only one point on Galactic plane.
I’m trying to simulate some celestial objects on Galactic plane now. I have a presentation at the Astronomical Society in this March, so it’s very helpful if you can support me fast. Thank you a lot.
#14 Updated by Knödlseder Jürgen over 6 years ago
I looked into this and there was indeed still a bug in a GammaLib class (
GModelSpatialDiffuseMap) for maps that were not normalized, i.e. the ones we use for the KSP simulation.
I fixed the bug and you can check out the latest code from gitlab.
I hope this solves all you problems.
Thanks for your feedback, it’s really useful in finding problems.