evselect table=PNclean.fits withspectrumset=yes spectrumset=PNsource_spectrum.fits energycolumn=PI spectralbinsize=5 withspecranges=yes specchannelmin=0 specchannelmax=20479 expression='(FLAG==0) && (PATTERN<=4) && ((X,Y) IN circle(26162.487, 23845.058, 341.85327))'


evselect table=PNclean.fits withspectrumset=yes spectrumset=PNbackground_spectrum.fits  energycolumn=PI spectralbinsize=5 withspecranges=yes specchannelmin=0 specchannelmax=20479 expression='(FLAG==0) && (PATTERN<=4) && ((X,Y) IN circle(22930.487, 26917.058, 2244.8385))'


backscale spectrumset=PNsource_spectrum.fits badpixlocation=PNclean.fits
backscale spectrumset=PNbackground_spectrum.fits badpixlocation=PNclean.fits


rmfgen spectrumset=PNsource_spectrum.fits rmfset=PN.rmf

arfgen spectrumset=PNsource_spectrum.fits arfset=PN.arf withrmfset=yes rmfset=PN.rmf  badpixlocation=PNclean.fits detmaptype=psf

specgroup spectrumset=PNsource_spectrum.fits mincounts=20 oversample=3 rmfset=PN.rmf arfset=PN.arf backgndset=PNbackground_spectrum.fits groupedset=PN_spectrum_grp.fits




#evselect table=MOS1clean.fits withspectrumset=yes spectrumset=MOS1source_spectrum.fits \
#   energycolumn=PI spectralbinsize=5 withspecranges=yes specchannelmin=0 specchannelmax=11999 \
#   expression='#XMMEA_EM && (PATTERN<=12) && ((X,Y) IN circle(26716.298,27824.699,744.11197))'

#evselect table=MOS1clean.fits withspectrumset=yes spectrumset=MOS1background_spectrum.fits \
#   energycolumn=PI spectralbinsize=5 withspecranges=yes specchannelmin=0 specchannelmax=11999 \
#   expression='#XMMEA_EM && (PATTERN<=12) && ((X,Y) IN circle(24642.224,24009.884,3091.7389))'

#backscale spectrumset=MOS1source_spectrum.fits badpixlocation=MOS1clean.fits
#backscale spectrumset=MOS1background_spectrum.fits badpixlocation=MOS1clean.fits

#rmfgen spectrumset=MOS1source_spectrum.fits rmfset=m1_e14_im_pall_c.rmf

#arfgen spectrumset=MOS1source_spectrum.fits arfset=MOS1.arf withrmfset=yes #rmfset=m1_e14_im_pall_c.rmf \
#   badpixlocation=MOS1clean.fits detmaptype=psf

#specgroup spectrumset=MOS1source_spectrum.fits mincounts=20 oversample=3 rmfset=m1_e14_im_pall_c.rmf \
#  arfset=MOS1.arf backgndset=MOS1background_spectrum.fits groupedset=MOS1_spectrum_grp.fits





#evselect table=MOS2clean.fits withspectrumset=yes spectrumset=MOS2source_spectrum.fits \
#   energycolumn=PI spectralbinsize=5 withspecranges=yes specchannelmin=0 specchannelmax=11999 \
#   expression='#XMMEA_EM && (PATTERN<=12) && ((X,Y) IN circle(26716.298,27824.699,779.00359))'

#evselect table=MOS2clean.fits withspectrumset=yes spectrumset=MOS2background_spectrum.fits \
#   energycolumn=PI spectralbinsize=5 withspecranges=yes specchannelmin=0 specchannelmax=11999 \
#   expression='#XMMEA_EM && (PATTERN<=12) && ((X,Y) IN circle(24494.076,24135.81,3063.5446))'

#backscale spectrumset=MOS2source_spectrum.fits badpixlocation=MOS2clean.fits
#backscale spectrumset=MOS2background_spectrum.fits badpixlocation=MOS2clean.fits

#rmfgen spectrumset=MOS2source_spectrum.fits rmfset=m2_e14_im_pall_c.rmf

#arfgen spectrumset=MOS2source_spectrum.fits arfset=MOS2.arf withrmfset=yes #rmfset=m2_e14_im_pall_c.rmf \
#   badpixlocation=MOS2clean.fits detmaptype=psf

#specgroup spectrumset=MOS2source_spectrum.fits mincounts=20 oversample=3 rmfset=m2_e14_im_pall_c.rmf \
#  arfset=MOS2.arf backgndset=MOS2background_spectrum.fits groupedset=MOS2_spectrum_grp.fits