mygaus_G0p9.xpos, mygaus_G0p9.ypos = skycoord_to_pixel(source_center_G0p9, on.wcs)
mygaus_G0p9.xpos.val += 0.5
mygaus_G0p9.ypos.val += 0.5
#Arc_source
arc_source=SkyCoord(0.130, -0.139, unit='deg', frame="galactic")
mygaus_arcsource = Gauss2D("Arc Source")
mygaus_arcsource.xpos, mygaus_arcsource.ypos = skycoord_to_pixel(arc_source, on.wcs)
# Modele pour le bkg base sur la carte de fond
load_table_model("bkg", "off_maps" + str(E1) + "_" + str(E2) + "_TeV.fits")
bkg.ampl = 1
load_table_model("exposure", "exp_maps" + str(E1) + "_" + str(E2) + "_TeV.fits")
exposure.ampl = 1
# Modele de CS
CS_map = SkyMap.read("CStot.fits")
cs_reproj = CS_map.reproject(on)
cs_reproj.data[np.where(np.isnan(cs_reproj.data))] = 0
cs_reproj.data[np.where(cs_reproj.data < 30)] = 0
cs_reproj.write("cs_map_reproj.fits", clobber=True)
load_table_model("CS", "cs_map_reproj.fits")
#set_full_model(bkg + psf_SgrA(mygaus_SgrA)+psf_G0p9(mygaus_G0p9))
#set_full_model(bkg + psf_SgrA(exposure*mygaus_SgrA)+psf_G0p9(exposure*mygaus_G0p9))
#set_full_model(bkg + psf_SgrA(exposure*(mygaus_SgrA+large_gaus * CS))+psf_G0p9(exposure*mygaus_G0p9))
set_full_model(bkg + psf_SgrA(mygaus_SgrA + large_gaus * CS) + psf_G0p9(mygaus_G0p9))
#set_full_model(bkg + psf_SgrA(mygaus_SgrA + large_gaus * CS + central_gauss) + psf_G0p9(mygaus_G0p9))
#set_full_model(bkg + psf_SgrA(mygaus_SgrA + large_gaus * CS + mygaus_arcsource) + psf_G0p9(mygaus_G0p9))
mygaus_SgrA.fwhm = 1
freeze(mygaus_SgrA.fwhm)
mygaus_G0p9.fwhm = 1
