aboutsummaryrefslogtreecommitdiffstats
path: root/mass_profile
diff options
context:
space:
mode:
Diffstat (limited to 'mass_profile')
-rwxr-xr-xmass_profile/cooling_time2.sh580
1 files changed, 0 insertions, 580 deletions
diff --git a/mass_profile/cooling_time2.sh b/mass_profile/cooling_time2.sh
deleted file mode 100755
index 775e175..0000000
--- a/mass_profile/cooling_time2.sh
+++ /dev/null
@@ -1,580 +0,0 @@
-#!/bin/sh
-#
-unalias -a
-export LC_COLLATE=C
-##
-## based on `ciao_r500avgt'
-## for calculating the `cooling time'
-## within (0.-0.048 r500) region
-##
-## Junhua Gu
-## August 22, 2012
-##
-## Weitian LI
-## 2013-04-28
-##
-##
-## Change logs:
-## 2016-05-28:
-## * Remove 'COSCALC', just use 'cosmo_calc'
-##
-
-## about, used in `usage' {{{
-VERSION="v2.0"
-UPDATE="2016-05-28"
-## about }}}
-
-## error code {{{
-ERR_USG=1
-ERR_DIR=11
-ERR_EVT=12
-ERR_BKG=13
-ERR_REG=14
-ERR_ASOL=21
-ERR_BPIX=22
-ERR_PBK=23
-ERR_MSK=24
-ERR_BKGTY=31
-ERR_SPEC=32
-ERR_ARF=51
-ERR_RMF=52
-ERR_UNI=61
-## error code }}}
-
-## usage, help {{{
-case "$1" in
- -[hH]*|--[hH]*)
- printf "usage:\n"
- printf " `basename $0` evt=<evt2_clean> r500=<r500_kpc> regin=<input_reg> regout=<output_reg> bkgd=<blank_evt | lbkg_reg | bkg_spec> nh=<nH> z=<redshift> arf=<warf_file> rmf=<wrmf_file> [ grpcmd=<grppha_cmd> log=<log_file> ]\n"
- printf "\nversion:\n"
- printf "${VERSION}, ${UPDATE}\n"
- exit ${ERR_USG}
- ;;
-esac
-## usage, help }}}
-
-## default parameters {{{
-# default `event file' which used to match `blanksky' files
-#DFT_EVT="_NOT_EXIST_"
-DFT_EVT="`ls evt2*_clean.fits 2> /dev/null`"
-# default `bkgd', use `bkgcorr_blanksky*' corrected bkg spectrum
-DFT_BKGD="`ls bkgcorr_blanksky_*.pi 2> /dev/null`"
-# default basedir
-DFT_BASEDIR="../.."
-# default info_json pattern
-DFT_JSON_PAT="*_INFO.json"
-# default `radial region file'
-#DFT_REG_IN="_NOT_EXIST_"
-DFT_REG_IN="rspec.reg"
-# default output region file (0.1-0.5 r500 region)
-DFT_REG_OUT="cooling.reg"
-# default ARF/RMF, the one of the outmost region
-DFT_ARF="`ls r1_*.warf 2> /dev/null`"
-DFT_RMF="`ls r1_*.wrmf 2> /dev/null`"
-
-# default `group command' for `grppha'
-#DFT_GRP_CMD="group 1 128 2 129 256 4 257 512 8 513 1024 16"
-DFT_GRP_CMD="group min 20"
-# default `log file'
-DFT_LOGFILE="cooling_`date '+%Y%m%d'`.log"
-
-INNER=0.0
-OUTER=0.048
-CT_RES="cooling_results.txt"
-# default output xspec scripts
-XSPEC_SCRIPT="xspec_cooling.xcm"
-# deproj xspec script, generated by `deproj_spectra'
-# from which get `nh' and `redshift'
-XSPEC_DEPROJ="xspec_deproj.xcm"
-## default parameters }}}
-
-## functions {{{
-# process commandline arguments
-# cmdline arg format: `KEY=VALUE'
-getopt_keyval() {
- until [ -z "$1" ]
- do
- key=${1%%=*} # extract key
- val=${1#*=} # extract value
- keyval="${key}=\"${val}\""
- echo "## getopt: eval '${keyval}'"
- eval ${keyval}
- shift # shift, process next one
- done
-}
-
-## background renormalization (BACKSCAL) {{{
-# renorm background according to particle background
-# energy range: 9.5-12.0 keV (channel: 651-822)
-CH_LOW=651
-CH_HI=822
-pb_flux() {
- punlearn dmstat
- COUNTS=`dmstat "$1[channel=${CH_LOW}:${CH_HI}][cols COUNTS]" | grep -i 'sum:' | awk '{ print $2 }'`
- punlearn dmkeypar
- EXPTIME=`dmkeypar $1 EXPOSURE echo=yes`
- BACK=`dmkeypar $1 BACKSCAL echo=yes`
- # fix `scientific notation' bug for `bc'
- EXPTIME_B=`echo ${EXPTIME} | sed 's/[eE]/\*10\^/' | sed 's/+//'`
- BACK_B=`echo "( ${BACK} )" | sed 's/[eE]/\*10\^/' | sed 's/+//'`
- PB_FLUX=`echo "scale = 16; ${COUNTS} / ${EXPTIME_B} / ${BACK_B}" | bc -l`
- echo ${PB_FLUX}
-}
-
-bkg_renorm() {
- # $1: src spectrum, $2: back spectrum
- PBFLUX_SRC=`pb_flux $1`
- PBFLUX_BKG=`pb_flux $2`
- BACK_OLD=`dmkeypar $2 BACKSCAL echo=yes`
- BACK_OLD_B=`echo "( ${BACK_OLD} )" | sed 's/[eE]/\*10\^/' | sed 's/+//'`
- BACK_NEW=`echo "scale = 16; ${BACK_OLD_B} * ${PBFLUX_BKG} / ${PBFLUX_SRC}" | bc -l`
- printf "\`$2': BACKSCAL:\n"
- printf " ${BACK_OLD} --> ${BACK_NEW}\n"
- punlearn dmhedit
- dmhedit infile=$2 filelist=none operation=add \
- key=BACKSCAL value=${BACK_NEW} comment="old value: ${BACK_OLD}"
-}
-## bkg renorm }}}
-## functions end }}}
-
-## parameters {{{
-# process cmdline args using `getopt_keyval'
-getopt_keyval "$@"
-
-## check log parameters {{{
-if [ ! -z "${log}" ]; then
- LOGFILE="${log}"
-else
- LOGFILE=${DFT_LOGFILE}
-fi
-printf "## use logfile: \`${LOGFILE}'\n"
-[ -e "${LOGFILE}" ] && mv -fv ${LOGFILE} ${LOGFILE}_bak
-TOLOG="tee -a ${LOGFILE}"
-echo "process script: `basename $0`" >> ${LOGFILE}
-echo "process date: `date`" >> ${LOGFILE}
-## log }}}
-
-# check given parameters
-# process `nh' and `redshift' {{{
-if [ ! -r "${XSPEC_DEPROJ}" ]; then
- printf "## the script \`${XSPEC_DEPROJ}' generated by \`deproj_spectra' NOT found\n" | ${TOLOG}
- printf "## please supply the value of \`nh' and \`redshift' manully\n" | ${TOLOG}
- read -p "> value of nh: " N_H
- read -p "> value of redshift: " REDSHIFT
-else
- # get `nh' and `redshift' from xspec script
- LN=`grep -n 'projct\*wabs\*apec' ${XSPEC_DEPROJ} | tail -n 1 | cut -d':' -f1`
- # calc the line number of which contains `nh'
- LN_NH=`expr ${LN} + 4`
- NH_XCM=`head -n ${LN_NH} ${XSPEC_DEPROJ} | tail -n 1 | awk '{ print $1 }'`
- # calc the line number of `redshift'
- LN_RS=`expr ${LN} + 7`
- RS_XCM=`head -n ${LN_RS} ${XSPEC_DEPROJ} | tail -n 1 | awk '{ print $1 }'`
- printf "## get value of nh: \`${NH_XCM}' (from \`${XSPEC_DEPROJ}')\n" | ${TOLOG}
- printf "## get value of redshift: \`${RS_XCM}' (from \`${XSPEC_DEPROJ}')\n" | ${TOLOG}
-
- ## if `nh' and `redshift' supplied in cmdline, then use them
- if [ ! -z "${nh}" ]; then
- N_H=${nh}
- else
- N_H=${NH_XCM}
- fi
- # redshift
- if [ ! -z "${z}" ]; then
- REDSHIFT=${z}
- else
- REDSHIFT=${RS_XCM}
- fi
-fi
-printf "## use nH: ${N_H}\n" | ${TOLOG}
-printf "## use redshift: ${REDSHIFT}\n" | ${TOLOG}
-# nh & redshift }}}
-
-# check basedir & json file
-if [ -d "${basedir}" ]; then
- BASEDIR=${basedir}
-else
- BASEDIR=${DFT_BASEDIR}
-fi
-if [ ! -z "${json}" ] && [ -r "${BASEDIR}/${json}" ]; then
- JSON_FILE="${BASEDIR}/${json}"
-elif [ `ls ${BASEDIR}/${DFT_JSON_PAT} 2> /dev/null | wc -l` -eq 1 ]; then
- JSON_FILE=`ls ${BASEDIR}/${DFT_JSON_PAT}`
-else
- read -p "> JSON_file: " JSON_FILE
- if [ ! -r "${JSON_FILE}" ]; then
- printf "ERROR: cannot access given \`${JSON_FILE}'\n"
- exit ${ERR_JSON}
- fi
-fi
-printf "## use json_file: \`${JSON_FILE}'\n" | ${TOLOG}
-
-# process `r500' {{{
-R500_RAW=`grep '"R500.*kpc' ${JSON_FILE} | sed 's/.*"R500.*":\ //' | sed 's/\ *,$//'`
-if [ ! -z "${r500}" ]; then
- R500_RAW=${r500}
-fi
-if [ -z "${R500_RAW}" ]; then
- printf "## input R500 followed with unit, e.g.: 800kpc, 400pix\n"
- read -p "> value of \`R500' (in pixel/kpc): " R500_RAW
-fi
-R500_VAL=`echo "${R500_RAW}" | tr -d 'a-zA-Z, '`
-R500_UNI=`echo "${R500_RAW}" | tr -d '0-9, '`
-printf "## get \`R500': ${R500_VAL} in unit \`${R500_UNI}'\n" | ${TOLOG}
-
-# if in kpc, convert to pix
-case "${R500_UNI}" in
- [pP]*)
- printf "## units in \`pixel', conversion not needed\n" | ${TOLOG}
- R500_PIX_B=`echo ${R500_VAL} | sed 's/[eE]/\*10\^/' | sed 's/+//'`
- ;;
- *)
- printf "## units in \`kpc', convert to \`Chandra pixel'\n" | ${TOLOG}
- KPC_PER_PIX=`cosmo_calc ${REDSHIFT} | grep 'kpc.*pix' | tr -d 'a-zA-Z_#=(),:/ '`
- # convert scientific notation for `bc'
- KPC_PER_PIX_B=`echo ${KPC_PER_PIX} | sed 's/[eE]/\*10\^/' | sed 's/+//'`
- printf "## calculated \`kpc/pixel': ${KPC_PER_PIX_B}\n"
- R500_VAL_B=`echo ${R500_VAL} | sed 's/[eE]/\*10\^/' | sed 's/+//'`
- R500_PIX_B=`echo "scale = 4; ${R500_VAL_B} / ( ${KPC_PER_PIX_B} )" | bc -l`
- ;;
-esac
-# calc (inner-outer R500)
-R_IN=`echo "scale = 4; ${INNER} * ${R500_PIX_B}" | bc -l`
-R_OUT=`echo "scale = 4; ${OUTER} * ${R500_PIX_B}" | bc -l`
-printf "## R500 in units pixel: ${R500_PIX_B}\n" | ${TOLOG}
-printf "## (${INNER}-${OUTER} R500) range in pixel: ${R_IN} - ${R_OUT}\n" | ${TOLOG}
-# r500 }}}
-
-# check evt file
-if [ -r "${evt}" ]; then
- EVT=${evt}
-elif [ -r "${DFT_EVT}" ]; then
- EVT=${DFT_EVT}
-else
- read -p "> clean evt2 file: " EVT
- if [ ! -r "${EVT}" ]; then
- printf "ERROR: cannot access given \`${EVT}' evt file\n"
- exit ${ERR_EVT}
- fi
-fi
-printf "## use evt file: \`${EVT}'\n" | ${TOLOG}
-
-# input and output region files {{{
-if [ -r "${regin}" ]; then
- REG_IN="${regin}"
-elif [ -r "${DFT_REG_IN}" ]; then
- REG_IN=${DFT_REG_IN}
-else
- read -p "> previous used radial spec regfile: " REG_IN
- if [ ! -r "${REG_IN}" ]; then
- printf "ERROR: cannot access given \`${REG_IN}' region file\n"
- exit ${ERR_REG}
- fi
-fi
-printf "## use previous regfile: \`${REG_IN}'\n" | ${TOLOG}
-if [ ! -z "${regout}" ]; then
- REG_OUT="${regout}"
-else
- REG_OUT=${DFT_REG_OUT}
-fi
-[ -e "${REG_OUT}" ] && mv -fv ${REG_OUT} ${REG_OUT}_bak
-printf "## set output regfile (0.1-0.5 r500 region): \`${REG_OUT}'\n" | ${TOLOG}
-
-# get center position from `regin'
-# only consider `pie' or `annulus'-shaped region
-TMP_REG=`grep -iE '(pie|annulus)' ${REG_IN} | head -n 1`
-XC=`echo ${TMP_REG} | tr -d ' ' | awk -F'[(),]' '{ print $2 }'`
-YC=`echo ${TMP_REG} | tr -d ' ' | awk -F'[(),]' '{ print $3 }'`
-printf "## get center coord: (${XC},${YC})\n" | ${TOLOG}
-# region files }}}
-
-# check given bkgd, determine background {{{
-if [ -r "${bkgd}" ]; then
- BKGD=${bkgd}
-elif [ -r "${DFT_BKGD}" ]; then
- BKGD=${DFT_BKGD}
-else
- read -p "> background (blanksky_evt | lbkg_reg | bkg_spec): " BKGD
- if [ ! -r "${BKGD}" ]; then
- printf "ERROR: cannot access given \`${BKGD}'\n"
- exit ${ERR_BKG}
- fi
-fi
-printf "## use bkgd: \`${BKGD}'\n" | ${TOLOG}
-# determine bkg type: blanksky, lbkg_reg, bkg_spec ?
-# according to file type first: text / FITS
-# if FITS, then get values of `HDUCLAS1' and `OBJECT'
-if file -bL ${BKGD} | grep -qi 'text'; then
- printf "## given \`${BKGD}' is a \`text file'\n"
- printf "## use it as local bkg region file\n"
- printf "## use *LOCAL BKG SPEC*\n" | ${TOLOG}
- # just set flags, extract spectrum later
- USE_LBKG_REG=YES
- USE_BLANKSKY=NO
- USE_BKG_SPEC=NO
-elif file -bL ${BKGD} | grep -qi 'FITS'; then
- printf "## given \`${BKGD}' is a \`FITS file'\n"
- # get FITS header keyword
- HDUCLAS1=`dmkeypar ${BKGD} HDUCLAS1 echo=yes`
- if [ "${HDUCLAS1}" = "EVENTS" ]; then
- # event file
- printf "## given file is \`event'\n"
- # check if `blanksky' or `stowed bkg'
- BKG_OBJ=`dmkeypar ${BKGD} OBJECT echo=yes`
- if [ "${BKG_OBJ}" = "BACKGROUND DATASET" ] || [ "${BKG_OBJ}" = "ACIS STOWED" ]; then
- # valid bkg evt file
- printf "## given FITS file is a valid bkgrnd file\n"
- printf "## use *BLANKSKY*\n" | ${TOLOG}
- USE_BLANKSKY=YES
- USE_LBKG_REG=NO
- USE_BKG_SPEC=NO
- # specify `BLANKSKY'
- BLANKSKY=${BKGD}
- else
- # invalid bkg evt file
- printf "ERROR: invalid bkg evt file given\n"
- exit ${ERR_BKGTY}
- fi
- elif [ "${HDUCLAS1}" = "SPECTRUM" ]; then
- # spectrum file
- printf "## given file is \`spectrum'\n"
- printf "## use *BKG SPECTRUM*\n" | ${TOLOG}
- USE_BKG_SPEC=YES
- USE_BLANKSKY=NO
- USE_LBKG_REG=NO
- # specify `BKG_SPEC'
- BKG_SPEC=${BKGD}
- else
- # other type
- printf "ERROR: other type FITS given\n"
- exit ${ERR_BKGTY}
- fi
-else
- printf "ERROR: given \`${BKGD}' type UNKNOWN\n"
- exit ${ERR_BKGTY}
-fi
-# bkgd }}}
-
-# check `arf' and `rmf' {{{
-if [ -r "${arf}" ]; then
- ARF=${arf}
-elif [ -r "${DFT_ARF}" ]; then
- ARF=${DFT_ARF}
-else
- read -p "> provide the ARF to use: " ARF
- if [ ! -r "${ARF}" ]; then
- printf "ERROR: cannot access given \`${ARF}'\n"
- exit ${ERR_ARF}
- fi
-fi
-printf "## use ARF: \`${ARF}'\n" | ${TOLOG}
-# rmf
-if [ -r "${rmf}" ]; then
- RMF=${rmf}
-elif [ -r "${DFT_RMF}" ]; then
- RMF=${DFT_RMF}
-else
- read -p "> provide the RMF to use: " RMF
- if [ ! -r "${RMF}" ]; then
- printf "ERROR: cannot access given \`${RMF}'\n"
- exit ${ERR_RMF}
- fi
-fi
-printf "## use RMF: \`${RMF}'\n" | ${TOLOG}
-# arf & rmf }}}
-
-# check given `grpcmd'
-if [ ! -z "${grpcmd}" ]; then
- GRP_CMD="${grpcmd}"
-else
- GRP_CMD="${DFT_GRP_CMD}"
-fi
-printf "## use grppha cmd: \`${GRP_CMD}'\n" | ${TOLOG}
-## parameters }}}
-
-##################################################
-#### main
-## D_A
-D_A_CM=`cosmo_calc ${REDSHIFT} | grep 'Angular_diameter_distance' | awk -F'[=,]' '{ print $2 }' | tr -d '[cm] '
-printf "D_A_CM(${REDSHIFT})= ${D_A_CM}\n"
-
-## region related {{{
-## generate the needed region file
-printf "generate the output region file ...\n"
-cat > ${REG_OUT} << _EOF_
-# Region file format: CIAO version 1.0
-pie(${XC},${YC},${R_IN},${R_OUT},0,360)
-_EOF_
-
-## open the evt file to verify or modify
-printf "## check the generated pie region ...\n"
-printf "## if modified, save with the same name \`${REG_OUT}' (overwrite)\n"
-ds9 ${EVT} -regions ${REG_OUT} -cmap sls -bin factor 4
-
-## check the (modified) region (pie region end angle)
-printf "check the above region (for pie region end angle) ...\n"
-INVALID=`grep -i 'pie' ${REG_OUT} | awk -F'[,()]' '$7 > 360'`
-if [ "x${INVALID}" != "x" ]; then
- printf "*** WARNING: there are pie regions' END_ANGLE > 360\n" | ${TOLOG}
- printf "*** will to fix ...\n"
- mv -fv ${REG_OUT} ${REG_OUT}_tmp
- # using `awk' to fix
- awk -F'[,()]' '{
- if ($7 > 360) {
- printf "%s(%.2f,%.2f,%.2f,%.2f,%.2f,%.2f)\n", $1,$2,$3,$4,$5,$6,($7-360)
- }
- else {
- print $0
- }
- }' ${REG_OUT}_tmp > ${REG_OUT}
- rm -f ${REG_OUT}_tmp
-fi
-## region related }}}
-
-## generate spectrum {{{
-# check counts
-punlearn dmlist
-CNT_RC=`dmlist infile="${EVT}[sky=region(${REG_OUT})][energy=700:7000]" opt=block | grep 'EVENTS' | awk '{ print $8 }'`
-if [ ${CNT_RC} -lt 500 ]; then
- F_WC=true
- WC="LOW_COUNTS"
- printf "*** WARNING: counts_in_0.048R500=${CNT_RC} < 500 ***\n"
-fi
-
-# object
-printf "extract object spectrum \`${AVGT_SPEC}' ...\n"
-AVGT_SPEC="${REG_OUT%.reg}.pi"
-AVGT_SPEC_GRP="${AVGT_SPEC%.pi}_grp.pi"
-punlearn dmextract
-dmextract infile="${EVT}[sky=region(${REG_OUT})][bin PI]" \
- outfile="${AVGT_SPEC}" wmap="[bin det=8]" clobber=yes
-# group spectrum
-printf "group object spectrum ...\n"
-grppha infile="${AVGT_SPEC}" outfile="${AVGT_SPEC_GRP}" \
- comm="${GRP_CMD} & exit" clobber=yes > /dev/null
-
-# background
-printf "generate the background spectrum ...\n"
-AVGT_BKG="${AVGT_SPEC%.pi}_bkg.pi"
-if [ "${USE_BLANKSKY}" = "YES" ]; then
- # use blanksky as background file
- printf "extract spectrum from blanksky ...\n"
- punlearn dmextract
- dmextract infile="${BLANKSKY}[sky=region(${REG_OUT})][bin PI]" \
- outfile=${AVGT_BKG} wmap="[bin det=8]" clobber=yes
-elif [ "${USE_LBKG_REG}" = "YES" ]; then
- printf "extract local background ...\n"
- punlearn dmextract
- dmextract infile="${EVT}[sky=region(${BKGD})][bin PI]" \
- outfile=${AVGT_BKG} wmap="[bin det=8]" clobber=yes
-elif [ "${USE_BKG_SPEC}" = "YES" ]; then
- printf "copy specified background spectrum ...\n"
- cp -fv ${BKG_SPEC} ${AVGT_BKG}
-fi
-
-printf "renormalize the background ...\n"
-bkg_renorm ${AVGT_SPEC} ${AVGT_BKG}
-## spectrum }}}
-
-## generate XSPEC script {{{
-printf "generate a XSPEC script ...\n"
-[ -e "${XSPEC_SCRIPT}" ] && mv -fv ${XSPEC_SCRIPT} ${XSPEC_SCRIPT}_bak
-cat > ${XSPEC_SCRIPT} << _EOF_
-## XSPEC script
-## spectrum analysis to get the average temperatue with (0-0.048 R500)
-##
-## generated by: \``basename $0`'
-## date: \``date`'
-##
-
-# xspec settings
-statistic chi
-abund grsa
-query yes
-
-# data
-data ${AVGT_SPEC_GRP}
-response ${RMF}
-arf ${ARF}
-backgrnd ${AVGT_BKG}
-
-# fitting range
-ignore bad
-ignore 0.0-0.7,7.0-**
-
-# plot related
-setplot energy
-
-method leven 10 0.01
-xsect bcmc
-cosmo 70 0 0.73
-xset delta 0.01
-systematic 0
-
-# model
-model wabs*apec
- ${N_H} -0.001 0 0 100000 1e+06
- 1.0 0.01 0.008 0.008 64 64
- 0.4 0.001 0 0 5 5
- ${REDSHIFT} -0.01 -0.999 -0.999 10 10
- 0.0 0.01 0 0 1e+24 1e+24
-
-## xspec script end
-
-proc calc_cooling_time {} {
- set rout ${R_OUT}
- set d_a_cm ${D_A_CM}
- fit 1000
- tclout param 4
- set z [ lindex \$xspec_tclout 0 ]
- tclout param 2
- set T [ lindex \$xspec_tclout 0 ]
- tclout param 5
- set norm [ lindex \$xspec_tclout 0 ]
- newpar 1 0
- dummyrsp .001 100
- flux .001 100
-
- tclout flux
- set flux [ lindex \$xspec_tclout 0 ]
- puts "flux(0.01-100kev): \$flux"
- set rout_cm [ expr \$rout*.492/3600/180*3.14159*\$d_a_cm ]
- set V [ expr 4./3.*3.14159*\$rout_cm*\$rout_cm*\$rout_cm ]
- set nenh [ expr \$norm*1E14*4*3.14159*\$d_a_cm*\$d_a_cm*(1+\$z)*(1+\$z)*(1+\$z)*(1+\$z)/\$V ]
- set d_l_cm [ expr \$d_a_cm*(1+\$z)*(1+\$z) ]
- set ne_np_ratio 1.2
- set ne [ expr sqrt(\$nenh*\$ne_np_ratio) ]
- set lx [ expr \$flux*4*3.14159*\$d_l_cm*\$d_l_cm ]
- set kb 1.602E-9
- set ct [ expr 3./2.*(\$ne+\$ne/\$ne_np_ratio)*\$kb*\$T*\$V/\$lx ]
- set ct_gyr [ expr \$ct/(3600*24*365.25*1E9) ]
- puts "Cooling_time= \$ct_gyr Gyr"
-}
-
-fit 1000
-calc_cooling_time
-
-tclexit
-_EOF_
-## xspec script }}}
-
-## invoke xspec to calc
-if [ "x${F_WC}" = "xtrue" ]; then
- printf "\n*** WC: LOW_COUNTS ***\n"
- printf "*** WARNING: counts_in_0.048R500=${CNT_RC} < 500 ***\n"
-else
- [ -e "${CT_RES}" ] && mv -f ${CT_RES} ${CT_RES}_bak
- printf "invoking XSPEC to calculate cooling time ...\n"
- xspec - ${XSPEC_SCRIPT} | tee ${CT_RES}
-
- OBS_ID=`grep '"Obs.*ID' ${JSON_FILE} | awk -F':' '{ print $2 }' | tr -d ' ,'`
- OBJ_NAME=`grep '"Source\ Name' ${JSON_FILE} | awk -F':' '{ print $2 }' | sed -e 's/\ *"//' -e 's/"\ *,$//'`
- CT=`grep -i '^Cooling_time' ${CT_RES} | awk '{ print $2 }'`
-
- printf "\n" | tee -a ${CT_RES}
- printf "# OBS_ID,OBJ_NAME,CT_gyr\n" | tee -a ${CT_RES}
- printf "# $OBS_ID,$OBJ_NAME,$CT\n" | tee -a ${CT_RES}
-fi
-
-exit 0
-