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

second

Browse Source
This commit is contained in:
Roman Krivonos 2024-04-13 15:05:44 +03:00
parent ce39dcc9c6
commit c1c726e728
4 changed files with 31 additions and 108 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
README.md
View File

@ -23,13 +23,13 @@ git clone git@heagit.cosmos.ru:integral/ridge.git
``` ```
cd ridge cd ridge
pip install --editable uds/ pip install --editable ridge/
``` ```
Обратите внимание на параметр **--editable**, он позволяет вам редактировать исходный код данного пакета и сразу его выполнять. Если вы не планируете модифицировать локальную копию кода, можете убрать этот параметр. Обратите внимание на параметр **--editable**, он позволяет вам редактировать исходный код данного пакета и сразу его выполнять. Если вы не планируете модифицировать локальную копию кода, можете убрать этот параметр.
После работы можете удалить проект: После работы можете удалить проект:
``` pip uninstall uds ``` ``` pip uninstall ridge ```
## Работа с данными ## Работа с данными

132
ridge/ridge/config.py
View File

@ -1,109 +1,33 @@
from pathlib import Path """ Определение областей "Галактика" и "Внегалактика" """
bmax=30.0
lmax=90.0
work_dir = Path('/path/to/some/logical/parent/dir') """ Определение допустимых значений для фазы орбиты ИНТЕГРАЛа """
phmin=0.2
phmax=0.8
""" Минимальный диапазон наблюдаемых фаз при калибровке модели фона на внегалактике """
phase_diff=0.25
""" Рабочий диапазон орбит """
revmin=64
revmax=2800
nmin=3
nmax=10
""" Максимальное смещение от Краба, где калибруется поток """
crab_sep_max = 4.5
""" Минимально допустимое число наблюдений Краба в одной орбите """
crab_nmax=10
""" Минимальное число наблюдений в одном пикселе карты """
nscw_min=10
""" """
Координаты сентрального кадра, к которому будут Пиксель карты со значением ошибки выше этого персентиля будет отброшен
приводиться изображения всех списков событий SEM означает standatd error on mean (~RMS/sqrt(n))
""" """
ra_cen=34.5342131 sem_cut=98
de_cen=-4.7956710
"""
Словарь камер со списком наблюдений каждой камеры.
Номера камер должны быть отсортированы
"""
keylist_tm={'1':['tm1_obs_1',],
'5':['tm5_obs_1','tm5_obs_2',],
'6':['tm6_obs_1','tm6_obs_2_badpix','tm6_obs_3_badpix',
'tm6_park_1','tm6_park_2','tm6_park_3','tm6_park_4',
'tm6_scan_1','tm6_scan_2','tm6_scan_3','tm6_scan_4'],
'7':['tm7_obs_1','tm7_obs_2',]}
"""
Примерные центры изображений каждого наблюдения.
Требуется для астрокоррекции. Будет расчитываться матрица сдвигов и поворотов,
так вот, повороты будут проводиться вокруг данных координат.
"""
wcslist={'tm1_obs_1':[34.7279760,-5.0680267],
'tm5_obs_1':[34.7351248,-4.4407314],
'tm5_obs_2':[34.8748997,-4.4871658],
'tm7_obs_1':[35.0015120,-4.7602124],
'tm7_obs_2':[34.9810029,-4.5915912],
'tm6_obs_1':[34.4227062,-4.7207170],
'tm6_obs_2_badpix':[34.7272339,-4.4294153],
'tm6_obs_3_badpix':[34.8750291,-4.4708468],
'tm6_park_1':[35.0544951,-4.0613619],
'tm6_park_2':[35.0558675,-4.0683084],
'tm6_park_3':[35.0565263,-4.0583538],
'tm6_park_4':[35.0602986,-4.0622220],
'tm6_scan_1':[34.5405596,-4.8088748],
'tm6_scan_2':[34.5405596,-4.8088748],
'tm6_scan_3':[34.5405596,-4.8088748],
'tm6_scan_4':[34.5405596,-4.8088748]}
""" like in the paper (Table 1) """
obslist={'tm1_obs_1':'N01',
'tm5_obs_1':'N02',
'tm5_obs_2':'N03',
'tm7_obs_1':'N15',
'tm7_obs_2':'N16',
'tm6_obs_1':'N12',
'tm6_obs_2_badpix':'N13',
'tm6_obs_3_badpix':'N14',
'tm6_park_1':'N04',
'tm6_park_2':'N06',
'tm6_park_3':'N08',
'tm6_park_4':'N10',
'tm6_scan_1':'N05',
'tm6_scan_2':'N07',
'tm6_scan_3':'N09',
'tm6_scan_4':'N11'}
""" Это просто индекс диапазона для выходных файлов. """
eband=[ "0", "1", "2", "3", "4", "5", "6", "7", "8"]
""" Диапазоны энергий. """
emin_ev=[ 300, 300, 600, 2300, 200, 300, 5000, 500, 1000]
emax_ev=[2300, 600, 2300, 5000, 10000,8000, 8000, 1000, 2000]
emin_kev=[0.3, 0.3, 0.6, 2.3, 0.2, 0.3, 5.0, 0.5, 1.0]
emax_kev=[2.3, 0.6, 2.3, 5.0, 10.0, 8.0, 8.0, 1.0, 2.0]
#ecf = [ 1.0, 1.0, 1.0, 1.0, 1.0, 1.0]
# something is wrong here
#ecf = [9.7817E+11, 3.2982E+12, 1.3903E+12, 2.3322E+12, 5.2022E+11, 5.8453E+11, 3.8468E+12]
"""
*** en0 ecf 9.7817E+11 +/- 2.4606E+10 2.52% N=17
*** en1 ecf 3.2982E+12 +/- 8.2963E+10 2.52% N=17
*** en2 ecf 1.3903E+12 +/- 3.5036E+10 2.52% N=17
*** en3 ecf 2.3322E+12 +/- 5.8717E+10 2.52% N=17
*** en4 ecf 5.2022E+11 +/- 1.3110E+10 2.52% N=17
*** en5 ecf 5.8453E+11 +/- 1.4743E+10 2.52% N=17
"""
# finally use Pavel's file ../data/ECF/ecf_tbabspow_g2nh0.02.pkl
"""
for e in [(0.3, 2.3), (0.3, 0.6), (0.6, 2.3), (2.3, 5.0), (5.0, 8.0)]:
print(f'{ecf[e]:g}')
"""
ecf = [1.0911e+12, # (0.3, 2.3)
1.07252e+12, # (0.3, 0.6)
1.08963e+12, # (0.6, 2.3)
1.14674e+11, # (2.3, 5.0)
1.0,
1.0,
2.77581e+10, # (5.0, 8.0)
1354632916123.6191, # (0.5, 1.0) 4XMM-DR12 EP2 band
1014764099304.4968] # (1.0, 2.0) 4XMM-DR12 EP3 band
outfile_post='.fits'
"""
Pavel Medvedev:
0.3-2.3: 9.135819435325375e-13
0.3-0.6: 9.160477830652834e-13
0.6-2.3: 9.201664167869427e-13
2.3-5.0: 8.721504826794627e-12
"""

2
scripts/01_init_events.py
View File

@ -9,7 +9,7 @@ import inspect
from os.path import dirname from os.path import dirname
import uds import ridge
from uds.config import * from uds.config import *
from uds.utils import * from uds.utils import *

1
scripts/README.md
View File

@ -1,7 +1,6 @@
Подготовка рабочего окружения: Подготовка рабочего окружения:
``` ```
source <MY PATH>/venv/bin/activate.csh source <MY PATH>/venv/bin/activate.csh
source <MY PATH>/eSASS4EDR/bin/esass-init.csh
``` ```