fluxscaleParseLog
Return to directory of Todd's CASA extensions
Extracts the spectral index fit parameters from the ASCII file produced by the casa task
fluxscale. See also
fluxscaleMean to operate on the dictionary, and
fluxscaleLogMean to extract mean flux densities from a series of executions.
Usage:
au.fluxscaleParseLog(logfile, field=None)
Inputs:
- field: use this parameter (a name string, e.g. 'J1924-2914') in order to pick a specific source, in the case that multiple sources are in the logfile. Only the first matching line containing the spidx string will be returned.
Returns:
- a list of the 4 items, the first 3 of which are needed by the standard='manual' option of the casa task setjy:
- flux density list: [I,0,0,0]
- reference frequency as a string (e.g. '44.5507GHz')
- spectral index as a float
- spectral index uncertainty
Examples:
CASA <8>: au.fluxscaleParseLog('flux.cal.txt')
Out[8]: ([0.846707, 0, 0, 0], '24.6475GHz', -0.432606)
CASA <9>: au.fluxscaleParseLog('flux.cal.txt',field='J1924')
Out[9]: ([12.3011, 0, 0, 0], '24.6475GHz', -0.378065)
Contents of log file:
# Flux density for J1832-1035 in SpW=0 (freq=2.36849e+10 Hz) is: 0.853721 +/- 0.0197959 (SNR = 43.1262, N = 25)
# Flux density for J1832-1035 in SpW=1 (freq=2.37169e+10 Hz) is: 0.86203 +/- 0.0198578 (SNR = 43.4101, N = 25)
# Flux density for J1832-1035 in SpW=2 (freq=2.38609e+10 Hz) is: 0.857586 +/- 0.0221797 (SNR = 38.6653, N = 25)
# Flux density for J1832-1035 in SpW=3 (freq=2.39569e+10 Hz) is: 0.859422 +/- 0.0213741 (SNR = 40.2085, N = 25)
# Flux density for J1832-1035 in SpW=4 (freq=2.40769e+10 Hz) is: 0.859113 +/- 0.0198957 (SNR = 43.1809, N = 25)
# Flux density for J1832-1035 in SpW=5 (freq=2.41329e+10 Hz) is: 0.857724 +/- 0.0226778 (SNR = 37.8222, N = 25)
# Flux density for J1832-1035 in SpW=6 (freq=2.45009e+10 Hz) is: 0.851167 +/- 0.0194383 (SNR = 43.7881, N = 25)
# Flux density for J1832-1035 in SpW=7 (freq=2.45249e+10 Hz) is: 0.849605 +/- 0.0211481 (SNR = 40.174, N = 25)
# Flux density for J1832-1035 in SpW=8 (freq=2.49198e+10 Hz) is: 0.842319 +/- 0.020341 (SNR = 41.4098, N = 25)
# Flux density for J1832-1035 in SpW=9 (freq=2.49518e+10 Hz) is: 0.840599 +/- 0.0215761 (SNR = 38.9598, N = 24)
# Flux density for J1832-1035 in SpW=10 (freq=2.49758e+10 Hz) is: 0.84174 +/- 0.0208597 (SNR = 40.3525, N = 24)
# Flux density for J1832-1035 in SpW=11 (freq=2.50478e+10 Hz) is: 0.838325 +/- 0.0221262 (SNR = 37.8883, N = 24)
# Flux density for J1832-1035 in SpW=12 (freq=2.52878e+10 Hz) is: 0.836628 +/- 0.0262306 (SNR = 31.8952, N = 24)
# Flux density for J1832-1035 in SpW=13 (freq=2.53198e+10 Hz) is: 0.838218 +/- 0.0229182 (SNR = 36.5743, N = 24)
# Flux density for J1832-1035 in SpW=14 (freq=2.56798e+10 Hz) is: 0.831277 +/- 0.0235539 (SNR = 35.2925, N = 24)
# Flux density for J1832-1035 in SpW=15 (freq=2.58718e+10 Hz) is: 0.828188 +/- 0.027608 (SNR = 29.9981, N = 23)
# Flux density for J1924-2914 in SpW=0 (freq=2.36849e+10 Hz) is: 12.55 +/- 0.0576067 (SNR = 217.857, N = 25)
# Flux density for J1924-2914 in SpW=1 (freq=2.37169e+10 Hz) is: 12.5135 +/- 0.0578109 (SNR = 216.456, N = 25)
# Flux density for J1924-2914 in SpW=2 (freq=2.38609e+10 Hz) is: 12.4194 +/- 0.0565607 (SNR = 219.577, N = 25)
# Flux density for J1924-2914 in SpW=3 (freq=2.39569e+10 Hz) is: 12.3557 +/- 0.0558014 (SNR = 221.422, N = 25)
# Flux density for J1924-2914 in SpW=4 (freq=2.40769e+10 Hz) is: 12.3553 +/- 0.0588083 (SNR = 210.094, N = 25)
# Flux density for J1924-2914 in SpW=5 (freq=2.41329e+10 Hz) is: 12.2701 +/- 0.0585931 (SNR = 209.411, N = 25)
# Flux density for J1924-2914 in SpW=6 (freq=2.45009e+10 Hz) is: 12.3153 +/- 0.0596394 (SNR = 206.496, N = 25)
# Flux density for J1924-2914 in SpW=7 (freq=2.45249e+10 Hz) is: 12.3949 +/- 0.0609402 (SNR = 203.394, N = 25)
# Flux density for J1924-2914 in SpW=8 (freq=2.49198e+10 Hz) is: 12.5108 +/- 0.198246 (SNR = 63.1075, N = 25)
# Flux density for J1924-2914 in SpW=9 (freq=2.49518e+10 Hz) is: 12.3846 +/- 0.0515656 (SNR = 240.171, N = 24)
# Flux density for J1924-2914 in SpW=10 (freq=2.49758e+10 Hz) is: 12.3967 +/- 0.0511078 (SNR = 242.56, N = 24)
# Flux density for J1924-2914 in SpW=11 (freq=2.50478e+10 Hz) is: 12.2587 +/- 0.052966 (SNR = 231.444, N = 24)
# Flux density for J1924-2914 in SpW=12 (freq=2.52878e+10 Hz) is: 12.1087 +/- 0.0526555 (SNR = 229.962, N = 24)
# Flux density for J1924-2914 in SpW=13 (freq=2.53198e+10 Hz) is: 12.064 +/- 0.0537709 (SNR = 224.358, N = 24)
# Flux density for J1924-2914 in SpW=14 (freq=2.56798e+10 Hz) is: 12.0506 +/- 0.0596768 (SNR = 201.931, N = 24)
# Flux density for J1924-2914 in SpW=15 (freq=2.58718e+10 Hz) is: 12.0492 +/- 0.0620363 (SNR = 194.229, N = 23)
# Fitted spectrum for J1832-1035 with fitorder=1: Flux density = 0.846707 +/- 0.000751838 (freq=24.6475 GHz) spidx=-0.432606 +/- 0.033
3691
# Fitted spectrum for J1924-2914 with fitorder=1: Flux density = 12.3011 +/- 0.0237904 (freq=24.6475 GHz) spidx=-0.378065 +/- 0.070374
9
--
ToddHunter - 2013-08-07