1
0
forked from erosita/uds
2025-04-21 12:44:25 +03:00
..
2023-03-09 21:10:39 +03:00
upd
2024-11-29 17:31:43 +03:00
2024-12-02 11:14:30 +03:00
upd
2024-11-29 17:31:43 +03:00
2025-04-21 12:44:25 +03:00
2025-04-21 12:44:25 +03:00
2025-04-21 12:44:25 +03:00
2025-04-21 12:44:25 +03:00
2024-03-18 14:46:38 +03:00
log
2025-02-03 16:39:33 +03:00
2024-05-15 16:21:15 +03:00
2024-03-18 14:46:38 +03:00
2024-03-18 14:46:38 +03:00
2023-04-25 19:16:05 +03:00

Подготовка рабочего окружения:

source <MY PATH>/venv/bin/activate.csh
source <MY PATH>/eSASS4EDR/bin/esass-init.csh

01_init_events.py

Создает начальные списки событий и помещает их в uds/data/processed Оригинальные файлы со списками событий задаются в файлах uds/data/evtlists/*.txt

02_merge_events.py

Создает объедененный список событий и помещает его в uds/products. Этот список событий нужен, в основном для извлечения спектров с помощью srctool.

Попутно этот скрипт унифицирует оригинальные списки событий для последующей обработки. А именно, корректируются слова OBS_MODE=POINING/SURVEY в зависимости от типа наблюдения и производится центрирование на одни и те же координаты с помощью команды radec2xy.

Для запуска адаптивного сглаживания do_adapt = True требуется запустить окружение ciao, так как нужна команда dmimgadapt

conda activate ciao-4.15
source <MY PATH>/eSASS4EDR/bin/esass-init.csh

03_init_obs.py

  1. Подготавливает списки событий в разных энергетических диапазонах.
  2. Запускает erbox в три этапа, чтобы получить рабочий список источников для ermldet.
  3. Запускает ermldet
  4. Делает кросс-корреляцию с каталогом Gaia-unWISE do_cross_match=True, которая создает три файла: .cross -- все пересечения, и .ref / .src -- входные каталоги для последующей команды wcs_match
  5. Делает матрицу преобразования координат и корректирует списки событий. Для запуска командwcs_match/wcs_update требуется запустить окружение ciao (см. выше)

04_mosaics.py

Создает сборные изображения (мозайки) в разных энергетических диапазонах.

05_scrtool.py

Запускает scrtool для самого широкого канала 0.2-10 кэВ, чтобы спектры имели самое полное покрытие по энергиям. Список источников берется из 0.3-2.3 кэВ.

Вычисляет ECF для всех диапазонов.

Делает принудительную фотометрию в выбранных каналах (параметрforced=True). Внимание! ermldet из eSASS4EDR не делает ассимитричные ошибки на потоки. Мы запускаем более последнюю версию ermldet (v1.56/2.18 esass_200412 Jul 2 12:04:46 2022). Для этого используется параметр local_run=True, который высвечивает какую команду надо запустить на другой машине и ждет ввода.