orca.tasks.run_image
====================

.. py:module:: orca.tasks.run_image

.. autoapi-nested-parse::

   Imaging pipeline with model subtraction and peeling.

   Implements a complete imaging workflow including:
   - Calibration application (delay + bandpass)
   - Antenna flagging
   - TTCal peeling
   - WSClean imaging (dirty + phase-centered)
   - Source blanking and uvsub

   This version includes split after uvsub for cleaner output.



Attributes
----------

.. autoapisummary::

   orca.tasks.run_image.DELAY_TABLE
   orca.tasks.run_image.BANDPASS_TBL
   orca.tasks.run_image.BAD_CORR_NUMS
   orca.tasks.run_image.JULIA_BIN
   orca.tasks.run_image.TTCAL_SCRIPT
   orca.tasks.run_image.SOURCES_JSON
   orca.tasks.run_image.WSCLEAN_DIRTY
   orca.tasks.run_image.WSCLEAN_PHASEC
   orca.tasks.run_image.BLANK_SIZE_DEG
   orca.tasks.run_image.LOG
   orca.tasks.run_image.ap


Functions
---------

.. autoapisummary::

   orca.tasks.run_image.copy_ms
   orca.tasks.run_image.find_freq
   orca.tasks.run_image.apply_calibration
   orca.tasks.run_image.run_flagging
   orca.tasks.run_image.run_ttcal_zest
   orca.tasks.run_image.wsclean_run
   orca.tasks.run_image.plot_dirty
   orca.tasks.run_image.blank_source
   orca.tasks.run_image.subtract_and_split
   orca.tasks.run_image.rephase_center
   orca.tasks.run_image.process_ms


Module Contents
---------------

.. py:data:: DELAY_TABLE
   :value: '/lustre/pipeline/calibration/delay/2024-05-24/20240524_delay.delay'


.. py:data:: BANDPASS_TBL
   :value: '/lustre/pipeline/calibration/bandpass/{freq}/2024-05-24/11/bandpass_concat.{freq}_11.bandpass'


.. py:data:: BAD_CORR_NUMS
   :value: (3, 12, 14, 17, 28, 31, 33, 34, 41, 44, 51, 56, 57, 79, 80, 87, 92, 117, 124, 126, 127, 137,...


.. py:data:: JULIA_BIN
   :value: '/opt/devel/pipeline/envs/julia060/bin/julia'


.. py:data:: TTCAL_SCRIPT
   :value: '/opt/devel/pipeline/envs/julia060/bin/ttcal.jl'


.. py:data:: SOURCES_JSON

.. py:data:: WSCLEAN_DIRTY
   :value: ['-pol', 'IV', '-size', '4096', '4096', '-scale', '0.03125', '-niter', '20000', '-weight',...


.. py:data:: WSCLEAN_PHASEC
   :value: ['-pol', 'IV', '-size', '48', '48', '-scale', '0.0833333', '-niter', '1000', '-beam-size',...


.. py:data:: BLANK_SIZE_DEG
   :value: 0.5


.. py:data:: LOG

.. py:function:: copy_ms(src: pathlib.Path, dst_dir: pathlib.Path) -> pathlib.Path

.. py:function:: find_freq(ms_path: pathlib.Path) -> str

   Return the first path component like '73MHz'.


.. py:function:: apply_calibration(ms: pathlib.Path, freq: str) -> None

.. py:function:: run_flagging(ms: pathlib.Path) -> None

.. py:function:: run_ttcal_zest(ms: pathlib.Path) -> None

.. py:function:: wsclean_run(ms: pathlib.Path, out_dir: pathlib.Path, prefix: str, args: list[str], tag: str)

.. py:function:: plot_dirty(dirty_fits: pathlib.Path, png_out: pathlib.Path)

.. py:function:: blank_source(model_fits: pathlib.Path, out_casa: pathlib.Path)

.. py:function:: subtract_and_split(ms_src: pathlib.Path, casa_model: pathlib.Path, workdir: pathlib.Path) -> pathlib.Path

   Returns path to *_UVSUB_CORR.ms* (split, corrected column).


.. py:function:: rephase_center(ms: pathlib.Path, ra_deg: float, dec_deg: float)

.. py:function:: process_ms(ms_path: str | pathlib.Path, workdir: str | pathlib.Path) -> pathlib.Path

.. py:data:: ap