ciaover /bin/rm -rf 13858 download_chandra_obsid 13858 gunzip 13858/primary/*gz gunzip 13858/secondary/*gz download_obsid_caldb 13858 ./CALDB export CALDB=`pwd`/CALDB export CALDBCONFIG=`pwd`/CALDB/software/tools/caldb.config export CALDBALIS=none acroread 13858/axaff13858N001_VV001_vv2.pdf & echo /bin/ls -1 13858/*/*evt*.fits ds9 13858/primary/acisf13858N001_evt2.fits -scale log & xpaset -p ds9 crosshair 4103 4115 physical xpaset -p ds9 saveimage png ds9_01.png display < ds9_01.png xpaget ds9 crosshair wcs fk5 sexagesimal xpaset -p ds9 mode region xpaset -p ds9 bin factor 4 xpaset -p ds9 saveimage png ds9_exercise04_a.png display < ds9_exercise04_a.png xpaset -p ds9 bin factor 8 xpaset -p ds9 saveimage png ds9_exercise04_b.png display < ds9_exercise04_b.png xpaset -p ds9 bin factor 1 xpaset -p ds9 zoom 0.125 xpaset -p ds9 saveimage png ds9_exercise04_c.png display < ds9_exercise04_c.png xpaset -p ds9 zoom to 1 xpaset -p ds9 bin filter 'energy=500:7000' xpaset -p ds9 saveimage png ds9_exercise05_a.png display < ds9_exercise05_a.png xpaset -p ds9 bin filter 'energy=2000:10000' xpaset -p ds9 saveimage png ds9_exercise05_extra1.png display < ds9_exercise05_extra1.png dmlist 13858/primary/acisf13858N001_evt2.fits header | grep TST xpaset -p ds9 bin filter 'time=456520000:456530000' xpaset -p ds9 saveimage png ds9_exercise05_extra2.png display < ds9_exercise05_extra2.png xpaset -p ds9 quit prism 13858/primary/acisf13858N001_evt2.fits & import -window prism prism_exercise7.png display < prism_exercise7.png xpaset -p prism quit dmlist 13858/primary/acisf13858N001_evt2.fits header,clean | \ egrep '^._TARG|^.*_NOM|^.*_PNT|^SIM_|^DTCOR|^ASCDSVER|^CALDBVER' dmkeypar 13858/primary/acisf13858N001_evt2.fits TIMEDEL echo+ dmmakepar 13858/primary/acisf13858N001_evt2.fits[events] \ dmmakepar_exercise08.par clob+ pdump dmmakepar_exercise08.par | \ egrep -i '^._TARG|^.*_NOM|^.*_PNT|^SIM_|^DTCOR|^ASCDSVER|^CALDBVER' plist dmmakepar_exercise08.par | \ egrep -i '^._TARG|^.*_NOM|^.*_PNT|SIM_|DTCOR|ASCDSVER|CALDBVER' pget dmmakepar_exercise08.par \ sim_x sim_y sim_z ra_pnt dec_pnt roll_pnt ra_nom dec_nom roll_nom \ dtcor caldbver ascdsver dmhistory 13858/primary/acisf13858N001_evt2.fits acis_process_events chandra_repro 13858 out= clob+ dmdiff 13858/primary/acisf13858N001_evt2.fits 13858/repro/acisf13858_repro_evt2.fits data- || echo dmkeypar 13858/primary/acisf13858N001_evt2.fits CALDBVER echo+ dmkeypar 13858/repro/acisf13858_repro_evt2.fits CALDBVER echo+ check_ciao_caldb dmlist 13858/primary/acisf13858N001_evt2.fits counts dmlist 13858/repro/acisf13858_repro_evt2.fits counts fluximage 13858/repro/acisf13858_repro_evt2.fits out=acisf13858 bin=1 bands=broad clob+ mode=h mkpsfmap acisf13858_broad_thresh.img acisf13858_broad.psfmap \ energy=1.4967 ecf=0.393 mode=h clob+ pset wavdetect \ infile=acisf13858_broad_thresh.img \ expfile=acisf13858_broad_thresh.expmap \ psffile=acisf13858_broad.psfmap \ scales='1.4 2 4 8 12 16 24 32' \ outfile=acisf13858_wav.src \ scellfile=acisf13858_wav.cell \ imagefile=acisf13858_wav.recon \ defnbkg=acisf13858_wav.nbkg \ clob+ wavdetect mode=h ds9 acisf13858_broad.psfmap -scale linear -cmap standard -zoom to fit \ -saveimage png ds9_exercise10_a.png -quit display < ds9_exercise10_a.png ds9 acisf13858_broad_thresh.img -block 2 -scale log -region acisf13858_wav.src \ -scale limits 0 20 -smooth -saveimage png ds9_exercise10_b.png -quit display < ds9_exercise10_b.png wavdetect scales='1 4 16 64' mode=h clobber=yes outfile=fewscales.src ds9 acisf13858_broad_thresh.img -block 2 -scale log -region fewscales.src \ -scale limits 0 20 -smooth -saveimage png ds9_exercise10_ec1.png -quit display < ds9_exercise10_ec1.png wavdetect sigthresh=1e-7 mode=h clobber=yes outfile=lowersigthresh.src ds9 acisf13858_broad_thresh.img -block 2 -scale log -region lowersigthresh.src \ -scale limits 0 20 -smooth -saveimage png ds9_exercise10_ec2.png -quit display < ds9_exercise10_ec2.png mkpsfmap acisf13858_broad_thresh.img ecf90.psfmap \ energy=1.4967 ecf=0.90 mode=h clob+ wavdetect mode=h clobber=yes outfile=ecf90.src psffile=ecf90.psfmap ds9 acisf13858_broad_thresh.img -block 2 -scale log -region lowersigthresh.src \ -scale limits 0 20 -smooth -saveimage png ds9_exercise10_ec3.png -quit display < ds9_exercise10_ec3.png mkpsfmap acisf13858_broad_thresh.img 5keV.psfmap \ energy=5.0 ecf=0.90 mode=h clob+ wavdetect mode=h clobber=yes outfile=5keV.src psffile=5keV.psfmap ds9 acisf13858_broad_thresh.img -block 2 -scale log -region 5keV.src \ -scale limits 0 20 -smooth -saveimage png ds9_exercise10_ec4.png -quit display < ds9_exercise10_ec4.png celldetect acisf13858_broad_thresh.img \ expstk=acisf13858_broad_thresh.expmap \ psffile=ecf90.psfmap \ out=acisf13858_cell.src \ clob+ mode=h ds9 acisf13858_broad_thresh.img -block 2 -scale log -region acisf13858_cell.src \ -scale limits 0 20 -smooth -saveimage png ds9_exercise10_ec5.png -quit display < ds9_exercise10_ec5.png vtpdetect acisf13858_broad_thresh.img \ exp=acisf13858_broad_thresh.expmap \ out=acisf13858_vtp.src \ clob+ mode=h ds9 acisf13858_broad_thresh.img -block 2 -scale log -region acisf13858_vtp.src[src_region] \ -scale limits 0 20 -smooth -saveimage png ds9_exercise10_ec6.png -quit display < ds9_exercise10_ec6.png ds9 acisf13858_broad_thresh.img -scale log -pan to 4096.5 4096.5 physical & echo "physical;circle(4104.5,4115.5,10)" | xpaset ds9 region xpaset -p ds9 saveimage png ds9_exercise11_a.png xpaget ds9 regions -format ds9 -system wcs -skyformat sexagesimal > ds9_src.reg xpaset -p ds9 regions delete all display < ds9_exercise11_a.png cat ds9_src.reg echo "physical;annulus(4105.5,4113.5,45,125) # background" | xpaset ds9 region xpaset -p ds9 saveimage png ds9_exercise11_b.png xpaget ds9 regions -format ciao -system physical > ds9_bkg.reg xpaset -p ds9 regions delete all display < ds9_exercise11_b.png cat ds9_bkg.reg xpaset -p ds9 quit dmcopy "acisf13858_broad_thresh.img[sky=region(ds9_src.reg)]" dmcopy_e12_src.fits clob+ dmlist dmcopy_e12_src.fits blocks,cols dmstat dmcopy_e12_src.fits cen+ sig- med+ dmcopy 13858/repro/acisf13858_repro_evt2.fits"[sky=region(ds9_src.reg)][bin sky=1]" \ dmcopy_e12_evtsrc.fits clob+ dmlist dmcopy_e12_evtsrc.fits blocks,cols dmstat dmcopy_e12_evtsrc.fits cen+ sig- med+ $ASCDS_INSTALL/bin/ds9 acisf13858_broad_thresh.img -scale log -pan to 4096.5 4096.5 physical & xpaset -p ds9 regions ds9_bkg.reg xpaset -p ds9 regions select all xpaset -p ds9 regions background xpaset -p ds9 regions ds9_src.reg xpaset -p ds9 saveimage png exercise_13_ds9a.png display < exercise_13_ds9a.png xpaset -p ds9 regions select all xpaset -p ds9 analysis task "{Net Counts}" import -window 'Net Counts' exercise_13_a.png display < exercise_13_a.png xpaset -p ds9 regions delete all echo "physical; circle(4104.4434,4115.8478,17.145457)" | xpaset ds9 regions echo "physical; annulus(4104.5,4115.5,26.305893,125) # background" | xpaset ds9 regions xpaset -p ds9 saveimage png exercise_13_ds9b.png display < exercise_13_ds9b.png xpaset -p ds9 regions select all xpaset -p ds9 analysis task "{Net Counts}" import -window 'Net Counts' exercise_13_b.png display < exercise_13_b.png xpaset -p ds9 quit punlearn dmcoords dmcoords 13858/repro/acisf13858_repro_evt2.fits op=sky x=4104.4755245 y=4115.8174825 celfmt=hms verb=0 pget dmcoords ra dec punlearn srcflux srcflux 13858/repro/acisf13858_repro_evt2.fits "09:14:49.088, +08:53:21.16" exercise14_step2 clob+ mode=h punlearn srcflux srcflux 13858/repro/acisf13858_repro_evt2.fits "09:14:49.088, +08:53:21.16" exercise14_step3 \ psfmethod=arfcorr \ clob+ mode=h verbose=0 cat exercise14_step3_summary.txt punlearn srcflux srcflux 13858/repro/acisf13858_repro_evt2.fits "09:14:49.088, +08:53:21.16" exercise14_step4 \ psfmethod=arfcorr \ model="xsbbody.black_body" paramvals="black_body.kT=1" \ clob+ mode=h verbose=0 cat exercise14_step4_summary.txt punlearn srcflux srcflux 13858/repro/acisf13858_repro_evt2.fits "09:14:49.088, +08:53:21.16" exercise14_step5 \ psfmethod=arfcorr \ model="xsbbody.black_body" paramvals="black_body.kT=1" \ band="hard" \ clob+ mode=h verbose=0 cat exercise14_step5_summary.txt punlearn srcflux srcflux 13858/repro/acisf13858_repro_evt2.fits "09:14:49.088, +08:53:21.16" exercise14_step_ec1 \ psfmethod=quick \ clob+ mode=h verbose=0 cat exercise14_step_ec1_summary.txt source $ASCDS_INSTALL/marx-5.4.0/setup_marx.sh punlearn srcflux srcflux 13858/repro/acisf13858_repro_evt2.fits "09:14:49.088, +08:53:21.16" exercise14_step_ec2 \ psfmethod=marx \ clob+ mode=h verbose=0 cat exercise14_step_ec2_summary.txt punlearn srcflux srcflux 13858/repro/acisf13858_repro_evt2.fits "09:14:49.088, +08:53:21.16" exercise14_step_ec4 \ psfmethod=quick band="2.0:10.0:3.0" \ clob+ mode=h verbose=0 cat exercise14_step_ec4_summary.txt source $ASCDS_INSTALL/marx-5.4.0/setup_marx.sh dmkeypar exercise14_step2_broad.flux net_photflux_aper echo+ curl \ -F email=${USER}@head.cfa.harvard.edu \ -F coords=cel \ -F ra=09:14:49.088 \ -F dec=+08:53:21.16 \ -F asol=obi \ -F obsid=13858 \ -F obinum=0\ -F niter=10\ -F randseed=32767 \ -F energy=mono \ -F mono=2.3 \ -F flux=0.000221 \ http://saotrace.cfa.harvard.edu/cgi-bin/runwrapper curl -o rays.tar.gz ftp://saotrace.cfa.harvard.edu/pub/ChaRT/HRMA_ra138.70453_dec8.88921_en2.3_flux0.000221_dithered.tar.gz tar xvfz rays.tar.gz simulate_psf 13858/repro/acisf13858_repro_evt2.fits \ chart_sim ra=138.70453 dec=8.88921 \ simulator=file rayfile="HRMA*rays.fits" \ minsize=128 mode=h ds9 acisf13858_broad_thresh.img -scale log -pan to 4096.5 4096.5 physical \ -zoom to 4 \ chart_sim.psf -pan to 4096.5 4096.5 physical \ -view colorbar no \ -saveimage png ds9_exercise15_a.png -quit display < ds9_exercise15_a.png simulate_psf 13858/repro/acisf13858_repro_evt2.fits \ marx_sim ra=138.70453 dec=8.88921 \ simulator=marx \ minsize=128 mode=h \ mono=2.3 flux=0.000221 numiter=10 ds9 acisf13858_broad_thresh.img -scale log -pan to 4096.5 4096.5 physical \ -zoom to 4 \ marx_sim.psf -pan to 4096.5 4096.5 physical \ -view colorbar no \ -saveimage png ds9_exercise15_b.png -quit display < ds9_exercise15_b.png simulate_psf 13858/repro/acisf13858_repro_evt2.fits \ marxA_sim ra=138.70453 dec=8.88921 \ simulator=marx \ minsize=128 mode=h \ mono=2.3 flux=0.000221 numiter=10 verb=0 readout+ simulate_psf 13858/repro/acisf13858_repro_evt2.fits \ marxB_sim ra=138.70453 dec=8.88921 \ simulator=marx \ minsize=128 mode=h \ mono=2.3 flux=0.000221 numiter=10 verb=0 pileup+ ext- simulate_psf 13858/repro/acisf13858_repro_evt2.fits \ marxC_sim ra=138.70453 dec=8.88921 \ simulator=marx \ minsize=128 mode=h \ mono=2.3 flux=0.000221 numiter=10 verb=0 blur=0.2 ds9 marx_sim.psf -scale log -pan to 4105.5 4113.5 physical \ -zoom to 4 \ marxA_sim.psf -pan to 4105.5 4113.5 physical \ marxB_sim.psf -pan to 4105.5 4113.5 physical \ marxC_sim.psf -pan to 4105.5 4113.5 physical \ -view colorbar no \ -saveimage png ds9_exercise15_ec1.png -quit display < ds9_exercise15_ec1.png arestore "acisf13858_broad_thresh.img[sky=bounds(region(ds9_bkg.reg))]" \ chart_sim.psf arestore.out num=150 clob+ ds9 acisf13858_broad_thresh.img -scale log -pan to 4105.5 4113.5 physical \ -zoom to 4 \ chart_sim.psf -pan to 4105.5 4113.5 physical \ arestore.out -pan to 4105.5 4113.5 physical \ -view colorbar no \ -tile mode column \ -saveimage png ds9_exercise15_ec2.png -quit display < ds9_exercise15_ec2.png /bin/rm -rf 13736 download_chandra_obsid 13736 evt1,bpix,flt,mtl,msk,dtf,bias,pbk,asol,stat download_obsid_caldb 13736 ./CALDB chandra_repro 13736 out= fluximage 13736 out=tricolor band=csc bin=8 clob+ ds9 -rgb -view colorbar no \ -red tricolor_soft_flux.img -scale log -smooth yes -smooth radius 1 \ -green tricolor_medium_flux.img -scale log -smooth yes -smooth radius 1 \ -blue tricolor_hard_flux.img -scale log -smooth yes -smooth radius 1 \ -saveimage png ds9_exercise_16.png -quit display < ds9_exercise_16.png dmimg2jpg \ infile=tricolor_soft_flux.img \ greenfile=tricolor_medium_flux.img \ bluefile=tricolor_hard_flux.img \ outfile=dmimg2jpg_exercise13.jpg \ clob+ display < dmimg2jpg_exercise13.jpg # Should convolve image and exposure map separately! aconvolve tricolor_soft_flux.img tricolor_soft_flux_gaus.img "lib:gaus(2,5,1,1,1)" meth=slide clob+ aconvolve tricolor_medium_flux.img tricolor_medium_flux_gaus.img "lib:gaus(2,5,1,1,1)" meth=slide clob+ aconvolve tricolor_hard_flux.img tricolor_hard_flux_gaus.img "lib:gaus(2,5,1,1,1)" meth=slide clob+ dmimg2jpg \ infile=tricolor_soft_flux_gaus.img \ greenfile=tricolor_medium_flux_gaus.img \ bluefile=tricolor_hard_flux_gaus.img \ outfile=dmimg2jpg_exercise13_aconvolve.jpg \ clob+ display < dmimg2jpg_exercise13_aconvolve.jpg # Csmooth really, really wants integer counts, so we smooth the counts image not the flux'ed image. csmooth tricolor_soft_thresh.img none tricolor_soft_flux_csm.img sigmin=3 sclmax=20 mode=h clob+ csmooth tricolor_medium_thresh.img none tricolor_medium_flux_csm.img sigmin=3 sclmax=20 mode=h clob+ csmooth tricolor_hard_thresh.img none tricolor_hard_flux_csm.img sigmin=3 sclmax=20 mode=h clob+ dmimg2jpg \ infile=tricolor_soft_flux_csm.img \ greenfile=tricolor_medium_flux_csm.img \ bluefile=tricolor_hard_flux_csm.img \ outfile=dmimg2jpg_exercise13_csm.jpg \ clob+ display < dmimg2jpg_exercise13_csm.jpg dmimgadapt tricolor_soft_thresh.img tricolor_soft_flux_cone.img cone min=1 max=20 num=100 radscale=linear counts=25 clob+ mode=h dmimgadapt tricolor_medium_thresh.img tricolor_medium_flux_cone.img min=1 max=20 num=100 radscale=linear counts=25 clob+ mode=h dmimgadapt tricolor_hard_thresh.img tricolor_hard_flux_cone.img min=1 max=20 num=100 radscale=linear counts=25 clob+ mode=h dmimg2jpg \ infile=tricolor_soft_flux_cone.img \ greenfile=tricolor_medium_flux_cone.img \ bluefile=tricolor_hard_flux_cone.img \ outfile=dmimg2jpg_exercise13_cone.jpg \ clob+ display < dmimg2jpg_exercise13_cone.jpg /bin/rm -rf 635 download_chandra_obsid 635 evt1,bpix,flt,mtl,msk,dtf,bias,pbk,asol,stat download_obsid_caldb 635 ./CALDB chandra_repro 635 out= dmstat "635/repro/acisf00635_repro_evt2.fits[cols time]" cen- sig- med- python -c 'a=72039166;b=72141142;dt3=(b-a)/3.0;print(a+0.0*dt3,a+1.0*dt3,a+2.0*dt3,a+3.0*dt3)' dmcopy "635/repro/acisf00635_repro_evt2.fits[time=72039166.0:72073158.0]" 635/repro/begin_evt.fits clob+ dmcopy "635/repro/acisf00635_repro_evt2.fits[time=72073158.0:72107150.0]" 635/repro/middle_evt.fits clob+ dmcopy "635/repro/acisf00635_repro_evt2.fits[time=72107150.0:72141142.0]" 635/repro/end_evt.fits clob+ fluximage 635/repro/begin_evt.fits rhooph_begin bin=6 band=broad clob+ verb=0 fluximage 635/repro/middle_evt.fits rhooph_middle bin=6 band=broad clob+ verb=0 fluximage 635/repro/end_evt.fits rhooph_end bin=6 band=broad clob+ verb=0 dmimg2jpg \ infile=rhooph_begin_broad_thresh.img \ greenfile=rhooph_middle_broad_thresh.img \ bluefile=rhooph_end_broad_thresh.img \ outfile=rhooph_time_slice.jpg \ clob+ display < rhooph_time_slice.jpg e=635/repro/acisf00635_repro_evt2.fits f=635/repro/acisf00635_repro_fov1.fits dmcopy "${e}[sky=region($f),energy=500:7000][bin sky=6,time=::#3]" 635_time_cube.fits clob+ ds9 -frame delete -scale log -rgbcube 635_time_cube.fits -saveimage ds9_ex16_timecube.png -quit display < ds9_ex16_timecube.png dmextract 13858/repro/acisf13858_repro_evt2.fits"[sky=circle(4104.5,4115.5,10),energy=500:7000][bin time=::100]" \ out=dme_100s.lc op=ltc1 \ bkg=13858/repro/acisf13858_repro_evt2.fits"[sky=region(ds9_bkg.reg),energy=500:7000]" \ mode=h clob+ cat << EOM > chips_ex17a.py from pychips import * from pychips.hlui import * load_preferences() make_figure("dme_100s.lc[cols dt,net_rate]") print_window("out_bin100.png","export.clobber=True") quit() EOM python chips_ex17a.py display < out_bin100.png for tt in 200 500 1000 1 do dmextract 13858/repro/acisf13858_repro_evt2.fits"[sky=circle(4104.5,4115.5,10),energy=500:7000][bin time=::${tt}]" \ out=dme_${tt}s.lc op=ltc1 \ bkg=13858/repro/acisf13858_repro_evt2.fits"[sky=region(ds9_bkg.reg),energy=500:7000]" \ mode=h clob+ done cat << EOM > chips_ex17b.py from pychips import * from pychips.hlui import * load_preferences() strip_chart(4) set_current_plot("plot1") add_curve("dme_200s.lc[cols dt,net_rate]") set_current_plot("plot2") add_curve("dme_500s.lc[cols dt,net_rate]") set_current_plot("plot3") add_curve("dme_1000s.lc[cols dt,net_rate]") set_current_plot("plot4") add_curve("dme_1s.lc[cols dt,net_rate]") print_window("out_bin_strip.png","export.clobber=True") quit() EOM python chips_ex17b.py display < out_bin_strip.png dmlist dme_100s.lc cols dmlist 13858/repro/acisf13858_repro_evt2.fits blocks,subspace glvary 13858/repro/acisf13858_repro_evt2.fits"[sky=circle(4104.5,4115.5,10),energy=500:7000]" \ out=glvary_src.prob lc=glvary_src.lc clob+ mode=h cat << EOM > chips_ex18a.py from pychips import * from pychips.hlui import * load_preferences() make_figure("glvary_src.lc[cols time,count_rate]") print_window("glvary_src.png","export.clobber=True") quit() EOM python chips_ex18a.py display < glvary_src.png glvary 13858/repro/acisf13858_repro_evt2.fits"[sky=region(ds9_bkg.reg),energy=500:7000]" \ out=glvary_bkg.prob lc=glvary_bkg.lc clob+ mode=h cat << EOM > chips_ex18b.py from pychips import * from pychips.hlui import * load_preferences() make_figure("glvary_bkg.lc[cols time,count_rate]") print_window("glvary_bkg.png","export.clobber=True") quit() EOM python chips_ex18b.py display < glvary_bkg.png punlearn ardlib acis_set_ardlib 13858/repro/acisf13858_repro_bpix1.fits abs- dither_region \ 13858/repro/pcadf456520092N001_asol1.fits \ "circle(4104.5,4115.5,10)" \ dr_src.out \ wcs= 13858/repro/acisf13858_repro_evt2.fits \ mask= 13858/repro/acisf13858_000N001_msk1.fits \ mode=h clob+ dmlist dr_src.out cols cat << EOM > chips_ex18_drout.py from pychips import * from pychips.hlui import * load_preferences() make_figure("dr_src.out[cols time,fracarea]") print_window("dr_src.png","export.clobber=True") quit() EOM python chips_ex18_drout.py display < dr_src.png dither_region \ 13858/repro/pcadf456520092N001_asol1.fits \ "region(ds9_bkg.reg)" \ dr_bkg.out \ wcs= 13858/repro/acisf13858_repro_evt2.fits \ mask= 13858/repro/acisf13858_000N001_msk1.fits \ mode=h clob+ maxpix=10000 cat << EOM > chips_ex18_droutb.py from pychips import * from pychips.hlui import * load_preferences() make_figure("dr_bkg.out[cols time,fracarea]") print_window("dr_bkg.png","export.clobber=True") quit() EOM python chips_ex18_droutb.py display < dr_bkg.png punlearn dmcoords dmcoords 13858/repro/acisf13858_repro_evt2.fits op=sky x=4104.5 y=4115.5 verb=0 pget dmcoords x y ra dec chip_id chipx chipy python -c 'from region import *;wav=CXCRegion("acisf13858_wav.src");out=field()-wav;out.write("exclude.reg",fits=True, clobber=True)' dmcopy 13858/repro/acisf13858_repro_evt2.fits"[sky=region(exclude.reg)]" 13858_evt.holes clob+ ds9 13858_evt.holes -scale log -bin factor 8 -saveimage png ds9_ex19.png -quit display < ds9_ex19.png dmextract 13858_evt.holes'[ccd_id=7][bin time=::259.28]' bkg_lc op=ltc1 clob+ cat << EOM > chips_ex19.py from pychips import * from pychips.hlui import * load_preferences() make_figure("bkg_lc[cols dt,count_rate]") print_window("chips_ex19.png","export.clobber=True") quit() EOM python chips_ex19.py display < chips_ex19.png deflare bkg_lc bkg_lc_flare method=clean save="bkg_lc_flare" convert bkg_lc_flare.ps -trim bkg_lc_flare.png display < bkg_lc_flare.png dmcopy 13858/repro/acisf13858_repro_evt2.fits"[@bkg_lc_flare]" 13858_deflare.evt clob+ dmlist 13858/repro/acisf13858_repro_evt2.fits counts dmlist 13858_deflare.evt counts dmextract 13858_evt.holes'[ccd_id=7][bin time=::50]' bkg_lc_50s op=ltc1 clob+ deflare bkg_lc_50s bkg_lc_flare_50s method=clean save="bkg_lc_flare_50s" plot- verb=0 convert bkg_lc_flare_50s.ps -trim bkg_lc_flare50s.png display < bkg_lc_flare50s.png deflare bkg_lc_50s bkg_lc_sigmaclip_50s method=sigma save="bkg_lc_sigmaclip_50s" plot- verb=0 < /dev/null convert bkg_lc_sigmaclip_50s.ps -trim bkg_lc_sigmaclip50s.png display < bkg_lc_sigmaclip50s.png dmextract 13858_evt.holes'[ccd_id=6][bin time=::259]' bkg_6lc op=ltc1 clob+ deflare bkg_6lc bkg_6lc_flare method=clean save="bkg_6lc_flare" plot- verb=0 convert bkg_6lc_flare.ps -trim bkg_6lc_flare.png display < bkg_6lc_flare.png punlearn specextract specextract 13858/repro/acisf13858_repro_evt2.fits"[sky=circle(4104.5,4115.5,10)]" \ out=specextract_01 \ bkgfile=13858/repro/acisf13858_repro_evt2.fits"[sky=region(ds9_bkg.reg)]" \ mode=h clob+ punlearn specextract specextract 13858/repro/acisf13858_repro_evt2.fits"[sky=circle(4104.5,4115.5,10)]" \ out=specextract_02 \ bkgfile=13858/repro/acisf13858_repro_evt2.fits"[sky=region(ds9_bkg.reg)]" \ weight=no \ mode=h clob+ punlearn specextract specextract 13858/repro/acisf13858_repro_evt2.fits"[sky=circle(4104.5,4115.5,10)]" \ out=specextract_03 \ bkgfile=13858/repro/acisf13858_repro_evt2.fits"[sky=region(ds9_bkg.reg)]" \ correctpsf=yes \ mode=h clob+ /bin/ls specextra*pi dmdiff specextract_01.pi specextract_02.pi || echo dmdiff specextract_01.pi specextract_03.pi || echo dmdiff specextract_01.arf"[#row=1:10]" specextract_02.arf"[#row=1:10]" comment- || echo dmstat specextract_01.arf"[cols specresp]" | egrep 'min|max' dmstat specextract_02.arf"[cols specresp]" | egrep 'min|max' dmstat specextract_03.arf"[cols specresp]" | egrep 'min|max' dmlist specextract_01.pi,specextract_01_bkg.pi header,clean | egrep 'BACKSCAL|EXPOSURE' rmfimg specextract_01.rmf specextract_01.rmf.img clob+ ds9 specextract_01.rmf.img -scale log -zoom to fit -view colorbar no -cmap heat \ -saveimage png ds9_exercise20.png -quit display < ds9_exercise20.png