![[PukiWiki]](http://www2.nao.ac.jp/~aste-www/aste_wiki/obs.jpg "[PukiWiki]")
#author("2021-04-05T21:22:51+09:00","default:astewww","astewww") #author("2021-08-04T10:44:24+09:00","default:astewww","astewww") [[ASTE Wiki]] * Demo [CI](3P1-3P0) map of a star forming region RCW38 † [#v961704a] :Last update| -- 2021-03-08 copyed from the ASTE Web site on Google. :Table of contents| #contents ---- ** Data archive [#rfc2e34c] - [[CASA data (image FITS, CASA MS2)>http://aste.nao.ac.jp/DemoData/CASA.tar]] (~1.5 GB) ** Data information [#i092b09f] | Target | RCW38 | | Observation date | 2016/11/26 | | Receiver | Band8 | | Spectrometer | MAC | | Mode | OTF mapping | | Frequency | 492.160656 GHz* | | Bandwidth | 512 MHz (Wide Band mode) | | Channel spacing | 0.5 MHz | | Velocity resolution | 0.3 km/s | | Beam size (FWHM) | 17 arcsec | | Mapping size | 6 arcmin x 8 arcmin | | Total mapping time | ~ 6 hrs | | Tsys | 800 - 1400 K | | PWV | ~ 0.5 mm | | Opacity (220 GHz) | 0.02 - 0.07 | - Literature value for frequency of [CI] is 492.16065 GHz. <= The difference (6 kHz) corresponds to ~ 3 m/s. ** CASA data reduction [#ta8cd472] - Measurement Set |CENTER: ||CENTER: ||c |CENTER: Scan &br; direction |CENTER: File |CENTER: Number &br; of row |CENTER: Note |h |X| R38otfX.RCW_38.20161126043008.4.ms | 81 || |~| R38otfX.RCW_38.20161126065817.4.ms | 81 || |~| R38otfX.RCW_38.20161126091618.4.ms | 81 || |Y| R38otfY.RCW_38.20161126052449.4.ms | 49 | <= Stop Observation at 49 row | |~| R38otfy.RCW_38.20161126060513.4.ms | 16 | <= Start observation from 46 row | |~| R38otfY.RCW_38.20161126075453.4.ms | 61 || |~| R38otfY.RCW_38.20161126101442.4.ms | 61 || -- Amplitude calibration with the chopper wheel method is already finished. -- The unit of amplitude in all measurement sets are Ta*. - Spectral window | 0, 1 (A1, A2) | [CI] (same signal (polarization) input into 0 and 1 (A1 and A2)) | | 2, 3 (A3, A4) | no data | - Software -- CASA version 5.1.1 *** Data reduction [#l2612c7d] + Baseline subtraction $ default(sdbaseline) $ sdbaseline( infile='R38otfX.RCW_38.20161126043008.4.ms', datacolumn='float_data', spw='0:350~470;560~650,1:350~470;560~650', # 350-470 ch and 560-650 ch correspond to -46.44~-9.93 km/s and 17.45~44.84 km/s, respectively. blfunc='poly', order=1, # fit function to subtract baselines outfile='R38otfX.RCW_38.20161126043008.4.A1A2.bl.ms', overwrite=True) // Task : sdbaseline (Fit/subtract a spectral baseline) // Parameters (examples) // infile = 'XXX.ms' # name of input SD dataset // datacolumn = 'float_data' # name of data column to be used ["data", "float_data", or "corrected"] // antenna = '' # select data by antenna name or ID, e.g. "PM03" // field = '' # select data by field IDs and names, e.g. "3C2*" (""=all) // spw = '0:350~470;560~650,1:350~470;560~650' # select data by IF IDs (spectral windows), e.g. "3,5,7" (""=all) // timerange = '' # select data by time range, e.g. "09:14:0~09:54:0" (""=all) (see examples in help) // scan = '' # select data by scan numbers, e.g. "21~23" (""=all) // pol = '' # select data by polarization IDs, e.g. "XX,YY" (""=all) // intent = '' # select data by observational intent, e.g. "*ON_SOURCE*" (""=all) // maskmode = 'list' # mode of setting additional channel masks // blmode = 'fit' # baselining mode ["fit" or "apply"] // dosubtract = True # subtract baseline from input data [True, False] // blformat = 'text' # format(s) of file(s) in which best-fit parameters are written // bloutput = '' # name(s) of file(s) in which best-fit parameters are written // blfunc = 'poly' # baseline model function // order = 1 # order of baseline model function // clipthresh = 3.0 # clipping threshold for iterative fitting // clipniter = 0 # maximum iteration number for iterative fitting // showprogress = False # (NOT SUPPORTED YET) show progress status for large data [True, False] (NOT SUPPORTED YET) // outfile = 'XXX.bl.ms' # name of output file // overwrite = False # overwrite the output file if already exists [True, False] + Scaling to convert Ta* to Tmb ++ Making table for scaling $ default(gencal) $ gencal( vis='R38otfX.RCW_38.20161126043008.4.A1A2.bl.ms', caltable=‘R38otfX.RCW_38.20161126043008.4.A1A2.bl.mb.tbl', caltype='amp', parameter=[sqrt(0.45)]) # Main beam efficiency of 45 % is used. # See the ASTE home page (https://alma.mtk.nao.ac.jp/aste/instruments/telescope.html). // Task : gencal (Specify Calibration Values of Various Types) // Parameters (examples) // vis = 'XXX.bl.ms' # Name of input visibility file // caltable = 'XXX.bl.mb.tbl’ # The new/existing calibration table // caltype = 'amp' # The calibration type: 'amp','ph','sbd','mbd','antpos','antposvla','tsys','evlagain', 'opac','gc','gceff','eff','tecim' // spw = '' # Calibration spw(s) selection // antenna = '' # Calibration antenna(s) selection // pol = '' # Calibration polarizations(s) selection // parameter = [sqrt(0.45)] # The calibration values ++ Applying the scaling table $ default(applycal) $ applycal( vis=‘R38otfX.RCW_38.20161126043008.4.A1A2.bl.ms', gaintable=['R38otfX.RCW_38.20161126043008.4.A1A2.bl.mb.tbl'], calwt=[False]) // Task : applycal (Apply calibrations solution(s) to data) // Parameters (examples) // vis = '○○○.bl.ms' # Name of input visibility file // field = '' # Select field using field id(s) or field name(s) // spw = '' # Select spectral window/channels // intent = '' # Select observing intent // selectdata = True # Other data selection parameters // timerange = '' # Select data based on time range // uvrange = '' # Select data within uvrange (default units meters) // antenna = '' # Select data based on antenna/baseline // scan = '' # Scan number range // observation = '' # Select by observation ID(s) // msselect = '' # Optional complex data selection (ignore for now) // docallib = False # Use callib or traditional cal apply parameters // gaintable = [] # Gain calibration table(s) to apply on the fly // gainfield = [] # Select a subset of calibrators from gaintable(s) // interp = [] # Interp type in time[,freq], per gaintable. default==linear,linear // spwmap = [] # Spectral windows combinations to form for gaintables(s) // calwt = [False] # Calibrate data weights per gaintable. // parang = False # Apply parallactic angle correction // applymode = '' # Calibration mode: ""="calflag","calflagstrict","trial","flagonly","flagonlystrict", or "calonly" // flagbackup = True # Automatically back up the state of flags before the run? + Imaging $ default(sdimaging) $ sdimaging( infiles='R38otfX.RCW_38.20161126043008.4.A1A2.bl.ms', 'R38otfX.RCW_38.20161126065817.4.A1A2.bl.ms', 'R38otfX.RCW_38.20161126091618.4.A1A2.bl.ms'], outfile='R38otfX.RCW_38.20161126.4.A1A2.bl.int', intent='*ON_SOURCE*', gridfunction='GJINC', cell=['6arcsec','6arcsec'], mode='velocity', nchan=167, start='-49.7km/s',width='0.6km/s', overwrite=True, imsize=[60,80], phasecenter='J2000 08:59:05.78 -047.31.00.020') -- We made two imaging data (Xscan and Yscan) //Task : sdimaging (Imaging for total power and spectral data) //Parameters (examples) //nfiles = ['○○○.bl.ms','×××.bl.ms',...] # a list of names of input SD Measurementsets (only MS is allowed for this task) //outfile = '○○○.int' # name of output image //overwrite = True # overwrite the output file if already exists [True, False] //field = '' # select data by field IDs and names, e.g. "3C2*" (""=all) //spw = '' # select data by IF IDs (spectral windows), e.g. "3,5,7" (""=all) //antenna = '' # select data by antenna names or IDs, e.g, "PM03" ("" = all antennas) //scan = '' # select data by scan numbers, e.g. "21~23" (""=all) //intent = '*ON_SOURCE*' # select data by observational intent, e.g. "*ON_SOURCE*" (""=all) //mode = 'velocity' # spectral gridding type ["channel", "frequency", "velocity"] //nchan = 167 # number of channels (planes) in output image (-1=all) //start = '-49.7km/s' # Velocity of first channel: e.g “0.0km/s" (""=first channel in first SpW of MS) //width = '0.6km/s' # Channel width e.g "-1.0km/s" (""=width of first channel in first SpW of MS) //veltype = 'radio' # velocity definition ["radio", "optical", "true" or "relativistic"] //outframe = '' # velocity frame of output image ["lsrk", "lsrd", "bary", "geo", "topo", "galacto", "lgroup", "cab"] (""=current frame or LSRK for multiple-MS inputs) //gridfunction = 'GJINC' # gridding function for imaging ["BOX", "SF", "PB", "GAUSS" or "GJINC"] (see description in help) //imsize = [60,80] # x and y image size in pixels, e.g., [64,64]. Single value: same for both spatial axes ([] = number of pixels to cover whole pointings in MSes) //cell = ['6arcsec','6arcsec'] # x and y cell size, (e.g., ["8arcsec","8arcsec"], default unit arcmin. ("" = 1/3 of FWHM of primary beam) //phasecenter = 'J2000 08:59:05.78 -047.31.00.020' # image center direction: position or field index, e.g., "J2000 17:30:15.0 -25.30.00.0". ("" = the center of pointing directions in MSes) //projection = 'SIN' # map projection type //ephemsrcname = '' # ephemeris source name, e.g. "MARS" //pointingcolumn = 'direction' # pointing data column to use ["direction", "target", “pointing_offset", "source_offset" or "encoder"] //restfreq = '' # rest frequency to assign to image, e.g., "114.5GHz" //stokes = '' # stokes parameters or polarization types to image, e.g. "I", "XX" //minweight = 0.1 # Minimum weight ratio to use //brightnessunit = '' # Overwrite the brightness unit in image ('' = respect the unit in MS) ['K' or 'Jy/beam'] //clipminmax = False # Clip minimum and maximum value from each pixel + Removing scanning effect (Basket-weave) $ default(sdfixscan) $ sdfixscan( infiles=['R38otfX.RCW_38.20161126.4.A1A2.bl.int', 'R38otfY.RCW_38.20161126.4.A1A2.bl.int'], direction=[0,90], maskwidth=[5.0,3.0], outfile='R38otfXY.RCW_38.20161126.4.A1A2.bl.int', overwrite=True) //Task : sdfixscan (Task for single-dish image processing) //Parameters (examples) //infiles = ['Xscan.int','Yscan.int'] # list of name of input SD images (FITS or CASA image) //mode = 'fft_mask' # image processing mode ["fft_mask", "model"] //direction = [0,90] # scan direction (p.a.) counterclockwise from the horizontal axis in unit of degree //maskwidth = [5.0,3.0] # mask width for Basket-Weaving (on percentage) //tmax = 0.0 # maximum threshold value for processing //tmin = 0.0 # minimum threshold value for processing //outfile = '○○○.int' # name of output file //overwrite = True # overwrite the output file if already exists [True, False] + Making Fits file $ default(exportfits) $ exportfits( imagename='R38otfXY.RCW_38.20161126.4.A1A2.bl.int', fitsimage='RCW_38.20161126.4.CI.casa.fits', velocity=True,dropstokes=True) //Task : exportfits (Convert a CASA image to a FITS file) //Parameters (examples) //imagename = '○○○.int' # Name of input CASA image //fitsimage = '○○○.fits' # Name of output image FITS file //velocity = True # Use velocity (rather than frequency) as spectral axis //optical = False # Use the optical (rather than radio) velocity convention //bitpix = -32 # Bits per pixel //minpix = 0 # Minimum pixel value (if minpix > maxpix, value is automatically determined) //maxpix = -1 # Maximum pixel value (if minpix > maxpix, value is automatically determined) //overwrite = True # Overwrite pre-existing imagename //dropstokes = True # Drop the Stokes axis? //stokeslast = True # Put Stokes axis last in header? //history = True # Write history to the FITS image? //dropdeg = False # Drop all degenerate axes (e.g. Stokes and/or Frequency)? ** Sample images [#ta9c1cdf] - Integrated intensity map &br; &ref(CIdemo-image-casa.png,,50%); -- Smoothing: gaussian, radius = 3 pixel (18 arcsec) -- The black filled circles at the lower right corner show the beam size (FWHM) of the Band8 (at 492 GHz) ** Acknowledgement and references in publications [#l3e2b53c] *** Acknowledgement [#xd3de2b1] - Please include the following statement in the acknowledgement, when you publish your studies using the ASTE demo science data: #br The ASTE demo science data collection was partially supported by JSPS KAKENHI Grant Nos. 15H02074 and 18H05441. *** References [#ra523fff] - Please include the following paper for the reference to the ASTE demo science data. #br Izumi, N. et al. 2001, PASJ, 73, 174 Izumi, N. et al. 2021, PASJ, 73, 174
