Browse Source

Merge pull request 'Implemented plotting of raw autocorrelations' (#8) from raw_quality_control into dev_spectral_line

Reviewed-on: #8
Ilse Bemmel, van 5 months ago
parent
commit
c031a32c69
  1. 76
      EVN_spectral_line.ipynb

76
EVN_spectral_line.ipynb

@ -74,7 +74,7 @@
"outputs": [],
"source": [
"# Point the mypath variable to where your data live\n",
"mypath='/home/jupyter/work/Spectral_line/Datacheck/EVN_CASA_pipeline/'\n",
"mypath='/home/jupyter/work/Spectral_line/RawQualityControl/EVN_CASA_pipeline/'\n",
"sys.path.append(mypath)"
]
},
@ -159,7 +159,7 @@
"outputs": [],
"source": [
"#Pair of target and associated phase calibrator source\n",
"target1 = 'J1849+3024'\n",
"target1 = 'J2202+4216'\n",
"phasecal1 = '1848+283'\n",
"\n",
"#Calibrated visibility data for the pair above\n",
@ -373,13 +373,6 @@
" flagmanager(vis=vis_line,mode='restore',versionname='precal_flags')"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "markdown",
"metadata": {},
@ -418,6 +411,71 @@
" print(\"A file with listobs information is already present\")"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Quality control of raw auto-correlations"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Make plots of the uncalibrated auto-correlations, phase and amplitude versus time and frequency, for the target source. Use the line pass for this. Plot for all stations."
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"# Implmented this with a single correlator pass, so the variable \"vis\" includes the full spectral resolution\n",
"# Need to find an elegant way to handle this, since for multiple correlator passes, the variable \"vis_line\"\n",
"# contains the full spectral resolution\n",
"\n",
"# For now only plot LL correlation, CASA has no option 'HALF' like AIPS to display the average of L/R\n",
"# Need to have bigger points in the plots, for some reason this makes the plotting even slower...\n",
"\n",
"#Amp vs frequency\n",
"plotms(vis=vis,xaxis='frequency',yaxis='amp',ydatacolumn='data',field=target1,\n",
" correlation='ll',coloraxis='antenna1',antenna='*&&&')"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"# Amp vs time\n",
"plotms(vis=vis,xaxis='time',yaxis='amp',ydatacolumn='data',field=target1,\n",
" correlation='ll',coloraxis='antenna1',antenna='*&&&')"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"# Phase vs frequency\n",
"plotms(vis=vis,xaxis='frequency',yaxis='phase',ydatacolumn='data',field=target1,\n",
" correlation='ll',coloraxis='antenna1',antenna='*&&&')"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"# Phase vs time\n",
"plotms(vis=vis,xaxis='time',yaxis='phase',ydatacolumn='data',field=target1,\n",
" correlation='ll',coloraxis='antenna1',antenna='*&&&')"
]
},
{
"cell_type": "markdown",
"metadata": {},

Loading…
Cancel
Save