ছায়াপথ স্তবক/চাও

উইকিবই থেকে
;;;;; Observation ID 11714

; Processing the evt1 and evt2 file

cd home/md/chandra/11714
chandra_repro  . ./new
cd ./new
dmcopy "acisf11714_000N001_repro_evt2.fits[energy=500:2500]" img11714.fits

; Processing the background file

cd ..
tar -xvf acis_bkgrnd_4.3.1.tar.gz
mkdir bkgrnd
cd data/chandra/acis/bkgrnd
ls
cp acis0iD2005-09-01bkgrnd_ctiN0001.fits acis2iD2005-09-01bkgrnd_ctiN0001.fits acis6iD2005-09-01bkgrnd_ctiN0001.fits acis1iD2005-09-01bkgrnd_ctiN0001A.fits acis3iD2005-09-01bkgrnd_ctiN0001.fits /home/md/chandra/bkgrnd
cd /home/md/chandra/bkgrnd
dmmerge "acis0iD2005-09-01bkgrnd_ctiN0001.fits, acis2iD2005-09-01bkgrnd_ctiN0001.fits, acis6iD2005-09-01bkgrnd_ctiN0001.fits, acis1iD2005-09-01bkgrnd_ctiN0001A.fits, acis3iD2005-09-01bkgrnd_ctiN0001.fits" bg_merge.fits
cp bg_merge.fits ../11714/new
cp bg_merge.fits ../11716/new
cp bg_merge.fits ../12037/new

; Subtracting cleaned background from cleaned image

cd ../11714/new
dmkeypar acisf11714_000N001_repro_evt2.fits GAINFILE echo+
dmkeypar bg_merge.fits GAINFILE echo+
dmlist bg_merge.fits header |grep PNT
dmmakepar acisf11714_000N001_repro_evt2.fits event_header.par
more event_header.par
grep _pnt event_header.par > event_pnt.par
more event_pnt.par
cat event_pnt.par
chmod +w bg_merge.fits
dmreadpar event_pnt.par "bg_merge.fits[events]" clobber+
dmlist bg_merge.fits header |grep PNT
punlearn reproject_events
pset reproject_events infile=bg_merge.fits
pset reproject_events outfile=bg11714.fits
pset reproject_events match=acisf11714_000N001_repro_evt2.fits
pset reproject_events aspect=pcadf376598961N001_asol1.fits
pset reproject_events random=0
plist reproject_events
reproject_events clobber+
dmcopy "bg11714.fits[energy=500:2500]" bg11714_cleaned.fits

dmkeypar "acisf11714_000N001_repro_evt2.fits" EXPOSURE echo+
dmkeypar "bg_merge.fits" EXPOSURE echo+
dmimgcalc infile=bg11714_cleaned.fits infile2=none out=bg_texp_11714.fits op=add weight=0.059350279
dmimgcalc img11714.fits bg_texp_11714.fits subimg11714.fits sub clobber+

; Exposure maps

dmcopy "subimg11714.fits[x=2700:5700,y=2700:5600]" cut_img11714.fits clobber+
dmkeypar acisf11714_000N001_repro_evt2.fits DETNAM echo+
merge_all evtfile="acisf11714_000N001_repro_evt2.fits" asol="pcadf376598961N001_asol1.fits" chip="0,1,2,3,6" refcoord="acisf11714_000N001_repro_evt2.fits" xygrid="2700:5700:#3000,2700:5600:#2900" energy=1 merged="merged_output_11714.fits" expmap="emap11714.fits" expcorr="exposure_corrected_11714.fits" clobber+
dmkeypar "acisf11714_000N001_repro_evt2.fits" EXPOSURE echo+

;;;;; Observation ID 11716

; Processing the evt1 and evt2 file

cd home/md/chandra/11716
chandra_repro  . ./new
cd ./new
dmcopy "acisf11716_000N001_repro_evt2.fits[energy=500:2500][status=0]" img11716.fits

; Subtracting cleaned background from cleaned image

cd ../11716/new
dmkeypar acisf11716_000N001_repro_evt2.fits GAINFILE echo+
dmkeypar bg_merge.fits GAINFILE echo+
dmlist bg_merge.fits header |grep PNT
dmmakepar acisf11716_000N001_repro_evt2.fits event_header.par
more event_header.par
grep _pnt event_header.par > event_pnt.par
more event_pnt.par
cat event_pnt.par
chmod +w bg_merge.fits
dmreadpar event_pnt.par "bg_merge.fits[events]" clobber+
dmlist bg_merge.fits header |grep PNT
punlearn reproject_events
reproject_events infile="bg_merge-fits" outfile="bg11716.fits" match="acisf11716_000N001_repro_evt2.fits" aspect="pcadf371565034N001_asol1.fits" random="0" clobber+
dmcopy "bg11716.fits[energy=500:2500]" bg11716_cleaned.fits

dmkeypar "acisf11716_000N001_repro_evt2.fits" EXPOSURE echo+
dmkeypar "bg_merge.fits" EXPOSURE echo+
dmimgcalc infile=bg11716_cleaned.fits infile2=none out=bg_texp_11716.fits op=add weight=0.025575389
dmimgcalc img11716.fits bg_texp_11716.fits subimg11716.fits sub clobber+

; Exposure maps

dmcopy "subimg11716.fits[x=2700:5700,y=2700:5600]" cut_img11716.fits clobber+
dmkeypar acisf11714_000N001_repro_evt2.fits DETNAM echo+
merge_all evtfile="acisf11716_000N001_repro_evt2.fits" asol="pcadf371565034N001_asol1.fits" chip="0,1,2,3,6" refcoord="acisf11716_000N001_repro_evt2.fits" xygrid="2700:5700:#3000,2700:5600:#2900" energy=1 merged="merged_output_11716.fits" expmap="emap11716.fits" expcorr="exposure_corrected_11716.fits" clobber+
dmkeypar "acisf11716_000N001_repro_evt2.fits" EXPOSURE echo+

;;;;; Observation ID 12037

; Processing the evt1 and evt2 file

cd home/md/chandra/12037
chandra_repro  . ./new
cd ./new
dmcopy "acisf12037_000N001_repro_evt2.fits[energy=500:2500][status=0]" img12037.fits

; Subtracting cleaned background from cleaned image

cd ../12037/new
dmkeypar acisf12037_000N001_repro_evt2.fits GAINFILE echo+
dmkeypar bg_merge.fits GAINFILE echo+
dmlist bg_merge.fits header |grep PNT
dmmakepar acisf12037_000N001_repro_evt2.fits event_header.par
more event_header.par
grep _pnt event_header.par > event_pnt.par
more event_pnt.par
cat event_pnt.par
chmod +w bg_merge.fits
dmreadpar event_pnt.par "bg_merge.fits[events]" clobber+
punlearn reproject_events
reproject_events infile="bg_merge.fits" outfile="bg12037.fits" match="acisf12037_000N001_repro_evt2.fits" aspect="pcadf376424703N001_asol1.fits" random="0" clobber+
dmcopy "bg12037.fits[energy=500:2500]" bg12037_cleaned.fits

dmkeypar "acisf12037_000N001_repro_evt2.fits" EXPOSURE echo+
dmkeypar "bg_merge.fits" EXPOSURE echo+
dmimgcalc infile=bg12037_cleaned.fits infile2=none out=bg_texp_12037.fits op=add weight=0.054603099
dmimgcalc img12037.fits bg_texp_12037.fits subimg12037.fits sub clobber+

; Exposure maps

dmcopy "subimg12037.fits[x=2700:5700,y=2700:5600]" cut_img12037.fits clobber+
dmkeypar acisf12037_000N001_repro_evt2.fits DETNAM echo+
merge_all evtfile="acisf12037_000N001_repro_evt2.fits" asol="pcadf376424703N001_asol1.fits" chip="0,1,2,3,6" refcoord="acisf12037_000N001_repro_evt2.fits" xygrid="2700:5700:#3000,2700:5600:#2900" energy=1 merged="merged_ourput_12037.fits" expmap="emap12037.fits" expcorr="exposure_corrected_12037.fits" clobber+
dmkeypar "acisf12037_000N001_repro_evt2.fits" EXPOSURE echo+

; Mosaic

cd /home/md/chandra
mkdir mosaic
stk_build infile="/home/md/chandra/11714/new/subimg11714.fits,/home/md/chandra/11716/new/subimg11716.fits,/home/md/chandra/12037/new/subimg12037.fits" outfile=bkgrnd_subtracted_images.lis clobber+
reproject_image_grid infile=@bkgrnd_subtracted_images.lis outfile=mosaic_subimg.fits xsize=3000 ysize=2900 xcenter="49.92" ycenter="41.56" pixelsize=0.5\" theta="0" method=sum clobber+
stk_build infile="/home/md/chandra/11714/new/emap11714.fits,/home/md/chandra/11716/new/emap11716.fits,/home/md/chandra/12037/new/emap12037.fits" outfile=exposure_maps.lis clobber+
reproject_image infile=@exposure_maps.lis matchfile=mosaic_subimg.fits outfile=mosaic_emap.fits method=average clobber+
dmimgcalc mosaic_subimg.fits mosaic_emap.fits mosaic.fits div clobber+
91992.932984718+39641.852314833+84634.804039914=216269.589339465
dmimgcalc infile=mosaic.fits op="imgout=img1*216269.589339465" out="abell426.fits" mode=h