source ~xraysw/cshrc  
         heainit
------------------------------------------------------------------
Below are notes from work on a different dataset, for illustration
------------------------------------------------------------------
cd /disk/xray3/gbr/jup_2014_XMM/0728960601/Jupiter

v120sas
printenv SAS_DIR
limit datasize unlimited

In /disk/xray3/gbr/jup_2014_XMM
cifbuild withobservationdate=YES observationdate="2014-04-15T20:55:00"
mv ccf.cif ccf_150414.cif

setenv SAS_CCF /disk/xray3/gbr/jup_2014_XMM/ccf_150414.cif
printenv SAS_CCF
/disk/xray3/gbr/jup_2014_XMM/ccf_150414.cif
ls /disk/xray3/gbr/jup_2014_XMM/ccf_150414.cif
/disk/xray3/gbr/jup_2014_XMM/ccf_150414.cif

setenv SAS_ODF /disk/xray3/gbr/jup_2014_XMM/0728960601/ODF
printenv SAS_ODF
/disk/xray3/gbr/jup_2014_XMM/0728960601/ODF
ls /disk/xray3/gbr/jup_2014_XMM/0728960601/ODF
2628_0728960601_SCX00000ATS_JUP.FIT  2628_0728960601_SCX00000ATS_SKY.FIT

xmmselect table=P0728960601PNS003PIEVLI0000.FIT &
Constructed 600x600 image with  (FLAG==0)&&(PATTERN<=4)
------------------------------------------
EPIC background filtering
http://xmm.esac.esa.int/external/xmm_user_support/documentation/sas_usg/USG/epicbkgfiltering.html

PN: < 0.4 ct/s for 
(PATTERN==0) && (PI in [10000:12000]) && #XMMEA_EP, on full FOV
(from http://xmm2.esac.esa.int/docs/documents/CAL-TN-0018.pdf)

If there are background flares:
Built background rate plot: Back_rates.ds/ps
tabgtigen table=Back_rates.ds gtiset=BKG_GTI.fits \
? expression='RATE < 0.4'
and BKG_GTI.fits got created

evselect table=P0728960801PNS003PIEVLI0000.FIT expression='gti(BKG_GTI.fits,TIME)' withfilteredset=y filteredset=GTI_flares_filtered_out.fits
------------------------------------------

Constructed 600x600 image with  (FLAG==0)&&(PATTERN<=4)&&(PI in [200:2000])
Set a circular region around N aurora, saved (Coordinate system WCS and 
Sexagesimal) in 
N_circle.reg    
circle(06:55:40.073,+23:07:08.13,9.3187")    9.31" appears to be the radius

Extracted counts histograms in 100 s bins for N and S spots and background: 

Displayed with fv (Grace does not work)
Saved ascii file (from Table data) in Fixed width columns

 counts samples in N_circle_200_2000_counts_100s_histo_t1.txt
 counts samples in back_circle_200_2000_counts_100s_histo_t1.txt

N starts at    5.139850654934E+08 and ends at 5.140261654934E+08 (412 bins)
back starts at 5.139908946936E+08 and ends at 5.140258946936E+08 (351 bins)
so curtailed N 5.139908654934E+08             5.140258654934E+08 (351 bins)

Ue IDL routine to check for periodicity

To check for IPT: Made contour images.
Box region of 5 pix side; b/g counts: 30, 24, 18, 40 (sum=112+/-10.6, 
average = 28 +/- 3)
IPT closest blob: 69 (-28 = 41 --> 13 sigma)
IPT farthest blob: 58 (-28 = 30 --> 10 sigma)
-----------------------------------------------------------------------