From decb76d661569de6b5056d4d3ed9dd1a2f453980 Mon Sep 17 00:00:00 2001 From: Roman Krivonos Date: Fri, 10 Oct 2025 15:51:45 +0300 Subject: [PATCH] spectra --- scripts/02_esas_image.py | 4 +- scripts/03_esas_cheese.py | 4 +- ...04_esass_spectra.py => 04_esas_spectra.py} | 27 +++++++-- scripts/05_esas_back.py | 58 ------------------- ...up_spectra.py => 05_esas_group_spectra.py} | 4 +- 5 files changed, 28 insertions(+), 69 deletions(-) rename scripts/{04_esass_spectra.py => 04_esas_spectra.py} (60%) delete mode 100755 scripts/05_esas_back.py rename scripts/{06_esas_group_spectra.py => 05_esas_group_spectra.py} (91%) diff --git a/scripts/02_esas_image.py b/scripts/02_esas_image.py index 8993992..04a1e93 100755 --- a/scripts/02_esas_image.py +++ b/scripts/02_esas_image.py @@ -29,7 +29,7 @@ print("Arches root path: {}".format(root_path)) archive_dir=root_path+'/data/archive' events_dir=root_path+'/data/processed' events_oot_dir=root_path+'/data/processed-oot' -products_dir=root_path+'/products' +products_dir=root_path+'/products/esas' ds9reg_dir=root_path+'/data/ds9reg' fast_track=True @@ -48,7 +48,7 @@ for obsid in files: if(obsid in skip): continue - work_dir = init_work_dir(obsid) + work_dir = init_work_dir(obsid, products_dir=products_dir) print(f"\n*** jump to {work_dir} ***\n") os.chdir(work_dir) diff --git a/scripts/03_esas_cheese.py b/scripts/03_esas_cheese.py index 067ce35..5b1eafd 100755 --- a/scripts/03_esas_cheese.py +++ b/scripts/03_esas_cheese.py @@ -27,7 +27,7 @@ print("Arches root path: {}".format(root_path)) archive_dir=root_path+'/data/archive' events_dir=root_path+'/data/processed' events_oot_dir=root_path+'/data/processed-oot' -products_dir=root_path+'/products' +products_dir=root_path+'/products/esas' ds9reg_dir=root_path+'/data/ds9reg' create_folder(products_dir) @@ -44,7 +44,7 @@ for obsid in files: if(obsid in skip): continue - work_dir = init_work_dir(obsid) + work_dir = init_work_dir(obsid, products_dir=products_dir) print(f"\n*** jump to {work_dir} ***\n") os.chdir(work_dir) diff --git a/scripts/04_esass_spectra.py b/scripts/04_esas_spectra.py similarity index 60% rename from scripts/04_esass_spectra.py rename to scripts/04_esas_spectra.py index 0f0030e..cebbc02 100755 --- a/scripts/04_esass_spectra.py +++ b/scripts/04_esas_spectra.py @@ -27,7 +27,7 @@ print("Arches root path: {}".format(root_path)) archive_dir=root_path+'/data/archive' events_dir=root_path+'/data/processed' events_oot_dir=root_path+'/data/processed-oot' -products_dir=root_path+'/products' +products_dir=root_path+'/products/esas' ds9reg_dir=root_path+'/data/ds9reg' create_folder(products_dir) @@ -38,19 +38,36 @@ t.run() files = glob.glob(archive_dir+'/*') + for obsid in files: obsid = os.path.basename(obsid) if(obsid in skip): continue - work_dir = init_work_dir(obsid) + work_dir = init_work_dir(obsid, products_dir=products_dir) print(f"\n*** jump to {work_dir} ***\n") os.chdir(work_dir) + + elow='350' + ehigh='1100' + + key='mos1S001' + run_mosspectra(key,'T T F T T F T',elow=elow,ehigh=ehigh) + run_mosback( key,'T T F T T F T') + run_bkgimsky( key,elow=elow,ehigh=ehigh) - run_mosspectra('mos1S001','T T F T T F T') - run_mosspectra('mos2S002','T T T T T T T') - run_pnspectra('pnS003','T T T T') + key='mos2S002' + run_mosspectra(key,'T T T T T T T',elow=elow,ehigh=ehigh) + run_mosback( key,'T T T T T T T') + run_bkgimsky( key,elow=elow,ehigh=ehigh) + + key='pnS003' + run_pnspectra(key,'T T T T',elow=elow,ehigh=ehigh) + run_pnback( key,'T T T T') + run_bkgimsky( key,elow=elow,ehigh=ehigh) + + #sys.exit() diff --git a/scripts/05_esas_back.py b/scripts/05_esas_back.py deleted file mode 100755 index c804350..0000000 --- a/scripts/05_esas_back.py +++ /dev/null @@ -1,58 +0,0 @@ -#!/usr/bin/env python - -from pysas.wrapper import Wrapper as w -import os, sys -from os.path import dirname -import inspect -import glob - -import os.path -from os import path -import subprocess -import numpy as np -import matplotlib.pyplot as plt -from astropy.io import fits -from astropy.table import Table -from matplotlib.colors import LogNorm - -import pyds9 - -import arches -from arches.utils import * -from arches.config import * - -root_path=dirname(dirname(dirname(inspect.getfile(arches)))) -print("Arches root path: {}".format(root_path)) - -archive_dir=root_path+'/data/archive' -events_dir=root_path+'/data/processed' -events_oot_dir=root_path+'/data/processed-oot' -products_dir=root_path+'/products' -ds9reg_dir=root_path+'/data/ds9reg' - -create_folder(products_dir) - -inargs = ['--version'] -t = w('sasver', inargs) -t.run() - -files = glob.glob(archive_dir+'/*') - -for obsid in files: - obsid = os.path.basename(obsid) - - if(obsid in skip): - continue - - work_dir = init_work_dir(obsid) - print(f"\n*** jump to {work_dir} ***\n") - os.chdir(work_dir) - - run_mosback('mos1S001','T T F T T F T') - run_mosback('mos2S002','T T T T T T T') - run_pnback('pnS003','T T T T') - - - - - diff --git a/scripts/06_esas_group_spectra.py b/scripts/05_esas_group_spectra.py similarity index 91% rename from scripts/06_esas_group_spectra.py rename to scripts/05_esas_group_spectra.py index 0e69dbd..6683a5e 100755 --- a/scripts/06_esas_group_spectra.py +++ b/scripts/05_esas_group_spectra.py @@ -27,7 +27,7 @@ print("Arches root path: {}".format(root_path)) archive_dir=root_path+'/data/archive' events_dir=root_path+'/data/processed' events_oot_dir=root_path+'/data/processed-oot' -products_dir=root_path+'/products' +products_dir=root_path+'/products/esas' ds9reg_dir=root_path+'/data/ds9reg' create_folder(products_dir) @@ -44,7 +44,7 @@ for obsid in files: if(obsid in skip): continue - work_dir = init_work_dir(obsid) + work_dir = init_work_dir(obsid, products_dir=products_dir) print(f"\n*** jump to {work_dir} ***\n") os.chdir(work_dir)