Return to directory of CASA extensions in analysisUtils.py
The toddTools package is an extension of simple python tasks for radio/millimeter/submillimeter astronomy and star formation studies mean to run inside the monolithic CASA environment.
It can also be imported into plain python, but tasks that require CASA tasks or tools will not function in that environment. Some of the tasks are reliant upon functions in the analysisUtils package of ALMA.
While toddTools is not (yet) publicly-available, subsets of the code will be made available upon reasonable request to the author.
A broken link here simply means that the command exists and works, but I have not yet written the help page for that command. But there is always python inline help available.
A
accretionEnergy -- Computes gravitational potential energy of a piece of gas falling onto an object
accretionLuminosity -- Computes gravitational potential luminosity of a piece of gas falling toward an object
accretionRadius -- Computes radius required to achieve given accretion rate, central mass, and accretion luminosity
accretionRate -- Computes accretion rate required to achieve accretion luminosity given central mass and radius
addLineToScripts -- Adds a line to the bottom of a script, or list of scripts
addLinks -- Adds soft links to a directory based on existing subdirectory names
adjustDivergingLut -- Removes the middle N values of a diverging RGB color lookup table and stretches the remainder together (to make better outflow palettes)
angularSeparationOfImagePeaks -- Finds the peak in two images and computes the angular separation in arcsec, and fractional beams
antennaCombinations -- Computes the number of N-antenna subsets of the specified array size
antennaVelocityDifference -- Computes maximum differential velocity of Earth rotation for a given east-west baseline and latitude
aperturePhotometry -- Uses rg.wpolygon to approximate (up to) two circles of different radii, and subtracts the flux of one from the other (works for continuum and cubes)
appendParameterToScript -- Takes a script of N lines containing a call of the same task and appends an additional parameter to the end of each one.
argmedian -- Returns a list of indices of an array corresponding to the value(s) equal to or closest to the median
atacCombinations -- Computes the fraction of 4-antenna combinations possible in the ATAC.HILSE if each fully-populated AVCC can select 1 of its 4 attached antennas
atmosphericVariation -- Computes the percentage variation in transmission and Tsky across an image cube's spectral axis
cardinalShiftRegion -- Reads a polygon region from a CRTF file, shifts it in each of the 4 cardinal directions and writes 4 new CRTF files (up, down, left, right)
checkPipelineChanchunks -- Compiles list of chanchunk values used for mfs and cube imaging in a list of pipeline runs
checkPipelineChecksource -- Checks if a checksource is included in the first EB of each pipeline run in a list of pipeline runs in a directory
checkPipelineFieldnames -- Checks if any pipeline image cube names contain integer field names and returns a dictionary of them
checkPipelineIntegrationTimes -- Returns a dictionary keyed by intent with values being the integration time in the first vis, and a Boolean on whether they differ
computeMom0peakOverMad -- Computes the ratio of the imstat 'max' over 'medabsdevmed' of a specified image
computeRcluster -- Compute the radius and mean position of a sehttpdrt://science.nrao.edu/facilities/alma/science_sustainability/na-alma-develpment-fy2023-call-for-study-proposalst of points.
concatSMAdataset -- Uses offsetScanNumbers to concatenate several measurement sets exported from MIR/IDL using autofits and previously filled to ms using MIRSFITStoCASA
constellation -- Calls astropy's get_constellation for a specified RA/Dec
convertBenchmarkScript -- Converts a regression script of findContinuum to a new ALMA cycle
convertSainiCSVFilesToFEICsquintFormat -- Converts per-Band .csv files dumped from per-Band .xls files created by Christopher at CDL (for Saini) into a single multi-Band csv file compatible with the original NA FEIC csv file
convolveChannel -- convolves a sinc function with the Hann window, and normalizes the resulting spectrum to unity
convolveChannelAveraging -- convolves a sinc function with the Hann window, then convolves with 2,4,8, or 16 delta functions and normalizes the resulting spectrum to unity
correctionsForRandomGauss --Computes corrections for the median and MAD of a population as judged by a sample from one end of the percentile distribution (developed for findContinuum)
correlateRuns -- Compares execution times of two sets of pipeline run pairs
correlateVisstatDifferences -- Computes fractional difference of per-scan/per-spw visibility statistics between two pipeline runs
countLinesPythonCode -- Counts lines of code in *.py files in specified directory tree
crtfMaskSolidAngle -- Computes solid angle of a mask image (such as that built from a crtf file using the CASA task makemask)
findContinuumForPixel -- Examines the intensity histogram of a spectrum and returns the channels that are close to the continuum level as located by the rolloff point
findContinuumHistogram -- Takes a cube and list of pixels, builds a histogram of flux density to identify continuum level, then finds the best continuum channels (calls findContinuumForPixel or findContinuumForPixelSCL)
findFrequencyForTau -- compute frequency for a specific dust opacity and value of beta
findGreyBodyFluxRatio -- computed the required tau to produce the observed flux ratio between two frequencies
findIntersection -- Finds the intersection point for 2 lines, each defined by a pair of points
findOuterAnnulusForPBCube -- Finds the scaled up values of a size mitigated cube that correspond to 0.2-0.3
imcombine -- Combines non-frequency-overlapping subcubes with the same frequency grid and spatial dimension into a single cube.
imfitBrightnessTemperature -- Reads a list of files of format imfitparseCLB and finds the one with peak intensity, then uses its fitted flux density and fitted size
imfitCube -- Runs imfit on all, or a list of channels in a CASA cube, generating a plot with errorbars
imfitEstimateChanges -- Computes the percentage change in the parameters between two files of imfit estimates
imfitIterate -- Runs imfit successively with the new estimates produced from the previous iteration until the result stops changing
imfitResultsToLatex -- Takes a list of imfit results and produces a Latex table
immoments4panel -- Builds moment images 0,1,2,8 of a CASA cube, then plots them on a 4-panel page
importfitshdu -- imports multiple HDUs from a single FITS image into multiple casa images
importsma -- imports one or more SMA UVFITS files that have been processed by uti_to_topo in MIR IDL
imstatBackground -- run imstat on the specified box region, then on the 4 adjacent boxes to determine the lowest rms in all 5 boxes
imstatCube -- compute the rms and standard deviation of chunks of a CASA image cube, and plots the result
imstatparse -- converts a dictionary returned by imstat into a formatted string
imsubimageRadec -- runs imsubimage centered on a specified RA Dec
imviewOverlayEllipse -- Overlay an ellipse using CASA imview -- obsoleted by au.imviewField(ellipse)
inclinationError -- computes the uncertainty on the inclination of a circular disk based on the apparent major and minor axes of its elliptical projection
maserUpperLimit -- Computes lower limit to the maser increase factor when only upper limits are available from one epoch
maskChannel -- Mask a channel or range of channels in a CASA image cube
matchImfitEstimates -- Copies imfit estimates files for a prior cube that match in velocity to a new cube
meanEpoch -- Uses au.yearFraction to compute mean epoch of a list of dates
meanPosition -- Takes two sexagesimal directions and computes the mean of both axes to produce a new direction
measuredProperMotion -- Takes a measured shift in RA seconds and Dec arcseconds between 2 epochs, and computes proper motion and extrapolates for a third epoch
pipelineNantennas -- Returns number of antennas in first EB of a pipeline run
pipelinePhasecalSolint -- Checks for solint specified gaincal as a value in 's'econds, rather than 'int' or 'inf'
pipelinePlotmsMPIFraction -- Searches a pipeline casa log file, counts number of lines with "Executing plotms" and counts again for lines also containing "MPI" earlier on the line
plotLinesOfCodeVsTime -- Plots the lines of code vs. cvs version number and vs. time
plotMaxHoldSpectrum -- Plots the maximum of a uv spectrum output by plotms text file tool
plotmsSWARM -- overlays the spws from many measurement sets into one plotms spectrum
plotNoiseCorrectionFactor -- plots the noise correction factor applied by the pipeline as defined in Todd's 13/Jul/2016 comment in CAS-8534 for the Hann window used by the ALMA baseline correlator (used for Fig 15 in the ALMA pipeline heuristics paper in PASP 2023)
plotRegionSpectrum -- computes and plots the spectrum from the unmasked pixels in a CASA region file
plotRobustGrid -- Plots the results file from buildRobustWeightingGrid
version -- returns the CVS version number of toddTools.py
visibleSky -- compute the fraction of the sky ever visible from a specified latitude and horizon limit
vislistFromFluxcsv -- Returns list of measurement sets contained a flux.csv file W
widthOfMaskArcsec -- Finds the width of the smallest central square that circumscribes a mask
wsuspw -- For a specified spw width (in units of FS) and desired spectral resolution, returns a dictionary keyed by nchan with values being a sub-dictionary of information for the parent setup and each CAF multiple that resultin an integer number of output channels of equal sensitivity
wsuspwMatrix Produces a terminal display (space-delimited for readability) and a text file (comma-delimited for import into spreadsheets) of possible spw combinations for each element in a list of number of adjacent frequency slices. Z
zams -- Computes luminosity and radius for a given stellar mass using Hosokawa's fits to Schaller et al. 1992