3c111
ra =04 18 21.3
dec=38 01 36
z=0.0485
Nh,Gal=2.91e21


==========================================================
REPROCESSING Chandra data using the new tool chandra_repro
==========================================================
# ---------------------------
# Check DATAMODE and LIVETIME
# ---------------------------
dmkeypar 14990/primary/*evt2* DATAMODE echo+
dmkeypar 14990/primary/*evt2* LIVETIME echo+

--> T92101s, VFAINT MODE

# ------------------
# Reprocess the data
# ------------------
punlearn chandra_repro
chandra_repro indir=14990 outdir=14990_new check_vf_pha+ verbose=3

# -----------------------
# Check relevant keywords
# -----------------------
dmkeypar 14990_new/acisf14990_repro_evt2.fits EXPOSURE echo+
dmkeypar 14990_new/acisf14990_repro_evt2.fits DATAMODE echo+
dmkeypar 14990_new/acisf14990_repro_evt2.fits DETNAM   echo+

--> VFAINT ACIS-56789 92101s

# ----------------------------------------------------
# Select good grades and data in the 0.5-7keV interval
# ----------------------------------------------------
punlearn dmcopy
dmcopy "14990_new/acisf14990_repro_evt2.fits[EVENTS][grade=0,2,3,4,6,status=0,energy=500:7000]" 14990_new/obs14990_057keV_repro_evt2.fits

# ----------------------------------------------------
# Extract images in the soft (0.5-2keV), hard (2-7keV) 
# and full (0.5-7keV) bands. Bin=2: 1 pixel~1 arcsec
# ----------------------------------------------------
punlearn dmcopy
dmcopy "14990_new/obs14990_057keV_repro_evt2.fits[ccd_id=7][bin X=2,Y=2]" 3c111_057keV_ima_bin2.fits
dmcopy "14990_new/obs14990_057keV_repro_evt2.fits[ccd_id=7][bin X=2,Y=2][energy=500:2000]" 3c111_052keV_ima_bin2.fits
dmcopy "14990_new/obs14990_057keV_repro_evt2.fits[ccd_id=7][bin X=2,Y=2][energy=2000:7000]" 3c111_27keV_ima_bin2.fits




===================
SPECTRAL EXTRACTION
===================
# ------------------------
# Visualize the event file
# ------------------------
ds9 obs14990_057keV_repro_evt2.fits.gz &

# -----------------------------------------------------
# Select the region files (CIAO format, physical scale) 
# for the source and the background
# -----------------------------------------------------
xemacs jet.reg &
# Region file format: CIAO version 1.0
rotbox(4006.1,4170.1,23.35,8.75,330)

xemacs back.reg &
# Region file format: CIAO version 1.0
circle(3955,4163,23.6)

#============================================================================
# Run SPECEXTRACT to create ACIS spectra and responses matrices
# correct=yes -> energy-dependent point-source aperture correction in the ARF
# weight=no   -> produces a point-like ARF
#============================================================================
punlearn specextract
pset specextract infile="obs14990_057keV_repro_evt2.fits.gz[sky=region(jet.reg)]"
pset specextract bkgfile="obs14990_057keV_repro_evt2.fits.gz[sky=region(back.reg)]"
pset specextract outroot=3c111_jet
pset specextract asp=pcadf474179932N001_asol1.fits.gz
pset specextract pbkfile=acisf474180206N001_pbk0.fits.gz
pset specextract mskfile=acisf14990_000N001_msk1.fits.gz
pset specextract badpixfile=acisf14990_repro_bpix1.fits.gz
pset specextract weight=yes
pset specextract correct=no
pset specextract verbose=2
pset specextract grouptype=NONE
pset specextract binspec=NONE
specextract

#-----------------------------------
# Grouping the data to 20 counts/bin
#-----------------------------------
grppha 3c111_jet.pi 3c111_jet_c20.pi comm="group min 20 & exit"


#--------------
# Running xspec
#--------------
xspec
data 3c111_jet_c20.pi
cpd /xw
setp ene
query yes
abund angr
ign **-0.3 7.-**
ign bad
pl lda

================================== The End ==================================