integral/ridge
1
0
Fork 0
generated from erosita/uds
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Spectra

Browse Source
This commit is contained in:
Roman Krivonos
2024-06-28 12:35:49 +03:00
parent ef604e08c8
commit 04843f27fc
72 changed files with 755 additions and 42 deletions
Download Patch File Download Diff File Expand all files Collapse all files

21
scripts/01_bgdmodel.sh
View File

@@ -26,3 +26,24 @@
./01_bgdmodel.py A10
./01_bgdmodel.py A11
./01_bgdmodel.py B01
./01_bgdmodel.py B02
./01_bgdmodel.py B03
./01_bgdmodel.py B04
./01_bgdmodel.py B05
./01_bgdmodel.py B06
./01_bgdmodel.py B07
./01_bgdmodel.py B08
./01_bgdmodel.py B09
./01_bgdmodel.py B10
./01_bgdmodel.py B11
./01_bgdmodel.py B12
./01_bgdmodel.py B13
./01_bgdmodel.py B14
./01_bgdmodel.py B15
./01_bgdmodel.py B16
./01_bgdmodel.py B17
./01_bgdmodel.py B18
./01_bgdmodel.py B19
./01_bgdmodel.py B20
./01_bgdmodel.py B21

9
scripts/01_crabmodel.py
View File

@@ -42,16 +42,17 @@ sigma=3
# for these bands, slope is taken from stacked profile
fixed_slope = ['E04','E05','E06','E07','E08','E09','E10','E11','E12','E15',
'A02','A03','A04','A05','A06','A07','A08','A09','A10','A11',]
'A02','A03','A04','A05','A06','A07','A08','A09','A10','A11',
'B02','B03','B04','B05','B06','B07','B08','B09','B10','B11','B12','B13','B14','B15','B16','B17','B18','B19','B20','B21']
# for these bands, slope is free for each revolution
free_slope = ['E01', 'E02', 'E03', 'A01']
free_slope = ['E01', 'E02', 'E03', 'A01','B01']
# for these bands, slope is fixed at constant (or if positive, which is not allowed)
const_slope = ['E10','E11','E12','A10','A11']
const_slope = ['E10','E11','E12','A10','A11','B18','B19','B20','B21']
# for stacked profile, skip orbits>800 for energy channels <30 keV
skip800 = ['E02','E03','A01']
skip800 = ['E02','E03','A01','B01']
# some static revs/scws to be removed
ignore_orbits=[352,834,912,1019,1021,1028,2275,2405,2493]

22
scripts/01_crabmodel.sh
View File

@@ -25,3 +25,25 @@
./01_crabmodel.py A09
./01_crabmodel.py A10
./01_crabmodel.py A11
./01_crabmodel.py B01
./01_crabmodel.py B02
./01_crabmodel.py B03
./01_crabmodel.py B04
./01_crabmodel.py B05
./01_crabmodel.py B06
./01_crabmodel.py B07
./01_crabmodel.py B08
./01_crabmodel.py B09
./01_crabmodel.py B10
./01_crabmodel.py B11
./01_crabmodel.py B12
./01_crabmodel.py B13
./01_crabmodel.py B14
./01_crabmodel.py B15
./01_crabmodel.py B16
./01_crabmodel.py B17
./01_crabmodel.py B18
./01_crabmodel.py B19
./01_crabmodel.py B20
./01_crabmodel.py B21

21
scripts/02_grxe_resid.sh
View File

@@ -26,6 +26,27 @@
./02_grxe_resid.py A10 ALL
./02_grxe_resid.py A11 ALL
./02_grxe_resid.py B01 ALL
./02_grxe_resid.py B02 ALL
./02_grxe_resid.py B03 ALL
./02_grxe_resid.py B04 ALL
./02_grxe_resid.py B05 ALL
./02_grxe_resid.py B06 ALL
./02_grxe_resid.py B07 ALL
./02_grxe_resid.py B08 ALL
./02_grxe_resid.py B09 ALL
./02_grxe_resid.py B10 ALL
./02_grxe_resid.py B11 ALL
./02_grxe_resid.py B12 ALL
./02_grxe_resid.py B13 ALL
./02_grxe_resid.py B14 ALL
./02_grxe_resid.py B15 ALL
./02_grxe_resid.py B16 ALL
./02_grxe_resid.py B17 ALL
./02_grxe_resid.py B18 ALL
./02_grxe_resid.py B19 ALL
./02_grxe_resid.py B20 ALL
./02_grxe_resid.py B21 ALL
#./02_grxe_resid.py E01 GAL

4
scripts/03_grxe_map.py
View File

@@ -41,8 +41,8 @@ from ridge.utils import *
from ridge.config import *
enkey = sys.argv[1]
key = sys.argv[2]
#key="ALL"
#key = sys.argv[2]
key="ALL"
fn='detcnts.{}.{}.resid.fits'.format(enkey,key)
d = fits.getdata(proddir+fn)

66
scripts/03_grxe_map.sh
View File

@@ -1,36 +1,36 @@
./03_grxe_map.py E02 ALL
./03_grxe_map.py E03 ALL
./03_grxe_map.py E04 ALL
./03_grxe_map.py E05 ALL
./03_grxe_map.py E06 ALL
./03_grxe_map.py E07 ALL
./03_grxe_map.py E08 ALL
./03_grxe_map.py E09 ALL
./03_grxe_map.py E10 ALL
./03_grxe_map.py E11 ALL
./03_grxe_map.py E12 ALL
./03_grxe_map.py E02
./03_grxe_map.py E03
./03_grxe_map.py E04
./03_grxe_map.py E05
./03_grxe_map.py E06
./03_grxe_map.py E07
./03_grxe_map.py E08
./03_grxe_map.py E09
./03_grxe_map.py E10
./03_grxe_map.py E11
./03_grxe_map.py E12
#./03_grxe_map.py E02 GAL
#./03_grxe_map.py E03 GAL
#./03_grxe_map.py E04 GAL
#./03_grxe_map.py E05 GAL
#./03_grxe_map.py E06 GAL
#./03_grxe_map.py E07 GAL
#./03_grxe_map.py E08 GAL
#./03_grxe_map.py E09 GAL
#./03_grxe_map.py E10 GAL
#./03_grxe_map.py E11 GAL
#./03_grxe_map.py E12 GAL
#./03_grxe_map.py E02 BKG
#./03_grxe_map.py E03 BKG
#./03_grxe_map.py E04 BKG
#./03_grxe_map.py E05 BKG
#./03_grxe_map.py E06 BKG
#./03_grxe_map.py E07 BKG
#./03_grxe_map.py E08 BKG
#./03_grxe_map.py E09 BKG
#./03_grxe_map.py E10 BKG
#./03_grxe_map.py E11 BKG
#./03_grxe_map.py E12 BKG
./03_grxe_map.py B01
./03_grxe_map.py B02
./03_grxe_map.py B03
./03_grxe_map.py B04
./03_grxe_map.py B05
./03_grxe_map.py B06
./03_grxe_map.py B07
./03_grxe_map.py B08
./03_grxe_map.py B09
./03_grxe_map.py B10
./03_grxe_map.py B11
./03_grxe_map.py B12
./03_grxe_map.py B13
./03_grxe_map.py B14
./03_grxe_map.py B15
./03_grxe_map.py B16
./03_grxe_map.py B17
./03_grxe_map.py B18
./03_grxe_map.py B19
./03_grxe_map.py B20
./03_grxe_map.py B21

28
scripts/03_grxe_spec.py
View File

@@ -40,6 +40,7 @@ n_bins = 80
sigma=3
plotme=False
"""
ebands0={#'E02':[0.0,0.0],
'E03':[0.0,0.0],
'E04':[0.0,0.0],
@@ -51,6 +52,31 @@ ebands0={#'E02':[0.0,0.0],
'E10':[0.0,0.0],
'E11':[0.0,0.0],
'E12':[0.0,0.0],}
"""
ebands0={'B01':[0.0,0.0],
'B02':[0.0,0.0],
'B03':[0.0,0.0],
'B04':[0.0,0.0],
'B05':[0.0,0.0],
'B06':[0.0,0.0],
'B07':[0.0,0.0],
'B08':[0.0,0.0],
'B09':[0.0,0.0],
'B10':[0.0,0.0],
'B11':[0.0,0.0],
'B12':[0.0,0.0],
'B13':[0.0,0.0],
'B14':[0.0,0.0],
'B15':[0.0,0.0],
'B16':[0.0,0.0],
'B17':[0.0,0.0],
'B18':[0.0,0.0],
'B19':[0.0,0.0],
'B20':[0.0,0.0],
'B21':[0.0,0.0],
}
#skey='Geminga'
if len(sys.argv) > 1:
@@ -130,7 +156,7 @@ for skey in skeys:
subprocess.run(["perl", "do_pha_igr_v3_mCrab.pl", fspec])
try:
for remfile in ["cols","cols1","cols2","header"]:
for remfile in ["cols","cols1","cols2","header",]:
os.remove(remfile)
except OSError:
pass

2
scripts/do_pha_igr_v3_mCrab.pl
View File

@@ -187,3 +187,5 @@ EOHEAD
close(HEADER);
system "fcreate headfile=header cdfile=cols datafile=$inp.dat outfile=$out.pha clobber=yes";
system "rm $inp.dat"