From c1c726e72881a1d941772bc6f5f246350213d6a6 Mon Sep 17 00:00:00 2001 From: Roman Krivonos Date: Sat, 13 Apr 2024 15:05:44 +0300 Subject: [PATCH] second --- README.md | 4 +- ridge/ridge/config.py | 132 ++++++++------------------------------ scripts/01_init_events.py | 2 +- scripts/README.md | 1 - 4 files changed, 31 insertions(+), 108 deletions(-) diff --git a/README.md b/README.md index d5bc14e..a90d86f 100644 --- a/README.md +++ b/README.md @@ -23,13 +23,13 @@ git clone git@heagit.cosmos.ru:integral/ridge.git ``` cd ridge -pip install --editable uds/ +pip install --editable ridge/ ``` Обратите внимание на параметр **--editable**, он позволяет вам редактировать исходный код данного пакета и сразу его выполнять. Если вы не планируете модифицировать локальную копию кода, можете убрать этот параметр. После работы можете удалить проект: -``` pip uninstall uds ``` +``` pip uninstall ridge ``` ## Работа с данными diff --git a/ridge/ridge/config.py b/ridge/ridge/config.py index ff0a2e0..5036fd8 100644 --- a/ridge/ridge/config.py +++ b/ridge/ridge/config.py @@ -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 -de_cen=-4.7956710 +sem_cut=98 - -""" -Словарь камер со списком наблюдений каждой камеры. -Номера камер должны быть отсортированы -""" -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 -""" diff --git a/scripts/01_init_events.py b/scripts/01_init_events.py index c8d9fc3..9229ee2 100755 --- a/scripts/01_init_events.py +++ b/scripts/01_init_events.py @@ -9,7 +9,7 @@ import inspect from os.path import dirname -import uds +import ridge from uds.config import * from uds.utils import * diff --git a/scripts/README.md b/scripts/README.md index 20c756f..0968533 100644 --- a/scripts/README.md +++ b/scripts/README.md @@ -1,7 +1,6 @@ Подготовка рабочего окружения: ``` source /venv/bin/activate.csh -source /eSASS4EDR/bin/esass-init.csh ```