evselect table=3885_0863941801_EPN_S003_ImagingEvts.ds withrateset=Y rateset=ratePN.fits  maketimecolumn=Y timebinsize=100 makeratecolumn=Y expression='#XMMEA_EP && (PI>10000&&PI<12000) && (PATTERN==0)'

tabgtigen table=ratePN.fits expression='RATE<=0.4' gtiset=PNgti.fits

evselect table=3885_0863941801_EPN_S003_ImagingEvts.ds withfilteredset=Y filteredset=PNclean.fits destruct=Y keepfilteroutput=T expression='#XMMEA_EP && gti(PNgti.fits,TIME) && (PI>150)'

evselect table=PNclean.fits imagebinning=binSize imageset=PNimage.fits withimageset=yes    xcolumn=X ycolumn=Y ximagebinsize=80 yimagebinsize=80

#evselect table=3726_0843440101_EMOS1_S001_ImagingEvts.ds withrateset=Y rateset=rateMOS1.fits \
#   maketimecolumn=Y timebinsize=100 makeratecolumn=Y \
#   expression='#XMMEA_EM && (PI>10000) && (PATTERN==0) '

#tabgtigen table=rateMOS1.fits expression='RATE<=0.35' gtiset=MOS1gti.fits

#evselect table=3726_0843440101_EMOS1_S001_ImagingEvts.ds withfilteredset=Y #filteredset=MOS1clean.fits \
#   destruct=Y keepfilteroutput=T \
#   expression='#XMMEA_EM && gti(MOS1gti.fits,TIME) && (PI>150)'

#evselect table=MOS1clean.fits imagebinning=binSize imageset=MOS1image.fits withimageset=yes    #xcolumn=X ycolumn=Y ximagebinsize=80 yimagebinsize=80

#evselect table=3726_0843440101_EMOS2_S002_ImagingEvts.ds withrateset=Y rateset=rateMOS2.fits \
#   maketimecolumn=Y timebinsize=100 makeratecolumn=Y \
#   expression='#XMMEA_EM && (PI>10000) && (PATTERN==0) '

#tabgtigen table=rateMOS2.fits expression='RATE<=0.35' gtiset=MOS2gti.fits

#evselect table=3726_0843440101_EMOS2_S002_ImagingEvts.ds withfilteredset=Y #filteredset=MOS2clean.fits \
#   destruct=Y keepfilteroutput=T \
#   expression='#XMMEA_EM && gti(MOS2gti.fits,TIME) && (PI>150)'

#evselect table=MOS2clean.fits imagebinning=binSize imageset=MOS2image.fits withimageset=yes    #xcolumn=X ycolumn=Y ximagebinsize=80 yimagebinsize=80