ALMA Renormalization Correction 6.6.1: Difference between revisions

From CASA Guides
Jump to navigationJump to search
Ekeller (talk | contribs)
No edit summary
Ekeller (talk | contribs)
No edit summary
 
(8 intermediate revisions by the same user not shown)
Line 2: Line 2:
ALMA observations of sources with strong line emission are affected by amplitude calibration issues caused by the ALMA normalization strategy, as described in the following Knowledgebase Articles: [https://help.almascience.org/kb/articles/what-are-the-amplitude-calibration-issues-caused-by-alma-s-normalization-strategy] and [https://help.almascience.org/kb/articles/what-errors-could-originate-from-the-correlator-spectral-normalization-and-tsys-calibration].
ALMA observations of sources with strong line emission are affected by amplitude calibration issues caused by the ALMA normalization strategy, as described in the following Knowledgebase Articles: [https://help.almascience.org/kb/articles/what-are-the-amplitude-calibration-issues-caused-by-alma-s-normalization-strategy] and [https://help.almascience.org/kb/articles/what-errors-could-originate-from-the-correlator-spectral-normalization-and-tsys-calibration].


If a correction factor of more than 2% is calculated for a dataset, it will be corrected during data reduction for all newly observed datasets. If you have a dataset where this correction has not been applied (this should be indicated in the Pipeline weblog or Manual Data Reduction Script), the following example script will calculate and apply the corrections to the dataset:
If a correction factor of more than 2% is calculated for a dataset, it will be corrected during data reduction for all newly observed datasets. If you have a dataset where this correction has not been applied (this should be indicated in the Pipeline weblog or Manual Data Reduction Script), the following example script will calculate and apply the corrections to the dataset. Modify the vislist parameter with a list of the ASDM(s) to check.
 
<source lang="python">
<source lang="python">
import os
__rethrow_casa_exceptions = True
__rethrow_casa_exceptions = True
vislist = ['uid___A002_Xe5aacf_X244e.ms']
for ii in range(1,len(vislist)+1):
    sessionIDs.append(f'session_{ii}')
h_init()
h_init()
try:
try:
     hifa_importdata(vis=['uid___A002_Xe3da01_X18fa.ms'], session=['session_1'], ocorr_mode='ca')
     hifa_importdata(vis=vislist, session=sessionIDs, ocorr_mode='ca')
     hifa_renorm(createcaltable=True, atm_auto_exclude=True)  
     hifa_renorm(createcaltable=True, atm_auto_exclude=True)  
     applycal(vis='uid___A002_Xe3da01_X18fa.ms',
     hif_mstransform()
      caltable = 'uid___A002_Xe3da01_X18fa.ms.hifa_renorm.s19_0.ampcal.tbl',
    for vis in vislist:
      intent='OBSERVE_TARGET#ON_SOURCE',
        targetsVis = vis.split('.')[0] + '_targets.ms'
      interp = 'nearest,linear'
        applycal(vis=targetsVis,
      calwt = True)
            gaintable = f'{vis}.hifa_renorm.s2_0.ampcal.tbl',
    hif_mstransform(pipelinemode='automatic')
            intent='OBSERVE_TARGET#ON_SOURCE',
            interp = 'nearest,linear',
            calwt = True)
finally:
finally:
     h_save()
     h_save()
</source>
</source>
This script should be run in the calibrated/working/ directory after you've restored the calibrated data. The resulting uid*.ms file will be corrected if the renormalization factor is >2%. The calculated renormalization factors along with plots can be viewed in the weblog produced by the script. For more information, please refer to the [https://almascience.nrao.edu/documents-and-tools/cycle10/alma_pipeline_users_guide_2023 Pipeline User's Guide]. Please note that the hifa_renorm task is only available in CASA version 6.2.1-7 and Pipeline version 2021.2.0.128 and above. The parameter atm_auto_exclude is only available in CASA version 6.4.1.12 with Pipeline version 2022.2.0.64 and above. The application of renorm was moved from within the hifa_renorm task to hif_applycal in CASA verion 6.6.1 and Pipeline version 2024.1.0.8. For older versions of CASA and the pipeline, please see [[ALMA_Renormalization_Correction]].
 
This script should be run in the calibrated/working/ directory after you've restored the calibrated data. The resulting uid*_targets.ms file will be corrected if the renormalization factor is >2%. Make sure to use the CORRECTED column. The calculated renormalization factors along with plots can be viewed in the weblog produced by the script. You are encouraged to reach out to [https://help.almascience.org the Helpdesk] for assistance in reviewing renormalization plots.
 
For more information, please refer to the [https://almascience.nrao.edu/documents-and-tools/cycle10/alma_pipeline_users_guide_2023 Pipeline User's Guide]. Please note that the hifa_renorm task is only available in CASA version 6.2.1-7 and Pipeline version 2021.2.0.128 and above. The parameter atm_auto_exclude is only available in CASA version 6.4.1.12 with Pipeline version 2022.2.0.64 and above. The application of renorm was moved from within the hifa_renorm task to hif_applycal in CASA verion 6.6.1 and Pipeline version 2024.1.0.8. For older versions of CASA and the pipeline, please see [[ALMA Renormalization Correction 6.5.4]].

Latest revision as of 14:23, 3 December 2024

Applying renormalization corrections to data

ALMA observations of sources with strong line emission are affected by amplitude calibration issues caused by the ALMA normalization strategy, as described in the following Knowledgebase Articles: [1] and [2].

If a correction factor of more than 2% is calculated for a dataset, it will be corrected during data reduction for all newly observed datasets. If you have a dataset where this correction has not been applied (this should be indicated in the Pipeline weblog or Manual Data Reduction Script), the following example script will calculate and apply the corrections to the dataset. Modify the vislist parameter with a list of the ASDM(s) to check.

import os
__rethrow_casa_exceptions = True
vislist = ['uid___A002_Xe5aacf_X244e.ms']
for ii in range(1,len(vislist)+1):
    sessionIDs.append(f'session_{ii}')
h_init()
try:
    hifa_importdata(vis=vislist, session=sessionIDs, ocorr_mode='ca')
    hifa_renorm(createcaltable=True, atm_auto_exclude=True) 
    hif_mstransform()
    for vis in vislist:
        targetsVis = vis.split('.')[0] + '_targets.ms'
        applycal(vis=targetsVis,
            gaintable = f'{vis}.hifa_renorm.s2_0.ampcal.tbl',
            intent='OBSERVE_TARGET#ON_SOURCE',
            interp = 'nearest,linear',
            calwt = True)
finally:
    h_save()

This script should be run in the calibrated/working/ directory after you've restored the calibrated data. The resulting uid*_targets.ms file will be corrected if the renormalization factor is >2%. Make sure to use the CORRECTED column. The calculated renormalization factors along with plots can be viewed in the weblog produced by the script. You are encouraged to reach out to the Helpdesk for assistance in reviewing renormalization plots.

For more information, please refer to the Pipeline User's Guide. Please note that the hifa_renorm task is only available in CASA version 6.2.1-7 and Pipeline version 2021.2.0.128 and above. The parameter atm_auto_exclude is only available in CASA version 6.4.1.12 with Pipeline version 2022.2.0.64 and above. The application of renorm was moved from within the hifa_renorm task to hif_applycal in CASA verion 6.6.1 and Pipeline version 2024.1.0.8. For older versions of CASA and the pipeline, please see ALMA Renormalization Correction 6.5.4.