Demo_science/ASTE_RCW38-CI
[
Front page
] [
New
|
Page list
|
Search
|
Recent changes
|
Help
]
Start:
[[ASTE Wiki]]
* Demo [CI](3P1-3P0) map of a star forming region RCW38 †...
: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.j...
** 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...
** CASA data reduction [#ta8cd472]
- Measurement Set
|CENTER: ||CENTER: ||c
|CENTER: Scan &br; direction |CENTER: File |CENTER: Numbe...
|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 Obs...
|~| R38otfy.RCW_38.20161126060513.4.ms | 16 | <= Start ob...
|~| R38otfY.RCW_38.20161126075453.4.ms | 61 ||
|~| R38otfY.RCW_38.20161126101442.4.ms | 61 ||
-- Amplitude calibration with the chopper wheel method is...
-- The unit of amplitude in all measurement sets are Ta*.
- Spectral window
| 0, 1 (A1, A2) | [CI] (same signal (polarization) input ...
| 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', ...
blfunc='poly', order=1, ...
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 data...
// datacolumn = 'float_data' # name of data column t...
// antenna = '' # select data by antenna name...
// field = '' # select data by field IDs an...
// spw = '0:350~470;560~650,1:350~470;560~650' ...
// timerange = '' # select data by time range, ...
// scan = '' # select data by scan numbers...
// pol = '' # select data by polarization...
// intent = '' # select data by observationa...
// maskmode = 'list' # mode of setting additional ...
// blmode = 'fit' # baselining mode ["fit" or "...
// dosubtract = True # subtract baseline from inpu...
// blformat = 'text' # format(s) of file(s) in whi...
// bloutput = '' # name(s) of file(s) in which...
// blfunc = 'poly' # baseline model function
// order = 1 # order of baseline model fun...
// clipthresh = 3.0 # clipping threshold for iter...
// clipniter = 0 # maximum iteration number fo...
// showprogress = False # (NOT SUPPORTED YET) show pr...
// outfile = 'XXX.bl.ms' # name of output file
// overwrite = False # overwrite the output file i...
+ 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...
caltype='amp',
parameter=[sqrt(0.45)]) # Main beam efficiency of 4...
# See the ASTE home page (h...
// Task : gencal (Specify Calibration Values of Various T...
// Parameters (examples)
// vis = 'XXX.bl.ms' # Name of input vi...
// caltable = 'XXX.bl.mb.tbl’ # The new/existing ...
// caltype = 'amp' # The calibration ...
// spw = '' # Calibration spw(...
// antenna = '' # Calibration ante...
// pol = '' # Calibration pola...
// parameter = [sqrt(0.45)] # The calibration...
++ Applying the scaling table
$ default(applycal)
$ applycal(
vis=‘R38otfX.RCW_38.20161126043008.4.A1A2.bl.ms',
gaintable=['R38otfX.RCW_38.20161126043008.4.A1A2.bl....
calwt=[False])
// Task : applycal (Apply calibrations solution(s) to data)
// Parameters (examples)
// vis = '○○○.bl.ms' # Name of input visi...
// field = '' # Select field using field id...
// spw = '' # Select spectral window/chan...
// intent = '' # Select observing intent
// selectdata = True # Other data selection parame...
// timerange = '' # Select data based on time r...
// uvrange = '' # Select data within uvrange ...
// antenna = '' # Select data based on antenn...
// scan = '' # Scan number range
// observation = '' # Select by observation ID(s)
// msselect = '' # Optional complex data selec...
// docallib = False # Use callib or traditional c...
// gaintable = [] # Gain calibration table(s) t...
// gainfield = [] # Select a subset of calibrat...
// interp = [] # Interp type in time[,freq],...
// spwmap = [] # Spectral windows combinatio...
// calwt = [False] # Calibrate data weights per ...
// parang = False # Apply parallactic angle cor...
// applymode = '' # Calibration mode: ""="calfl...
// flagbackup = True # Automatically back up the s...
+ 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=...
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 ...
//Parameters (examples)
//nfiles = ['○○○.bl.ms','×××.bl.ms',......
//outfile = '○○○.int' #...
//overwrite = True # o...
//field = '' ...
//spw = '' ...
//antenna = '' ...
//scan = '' ...
//intent = '*ON_SOURCE*' # select d...
//mode = 'velocity' # ...
//nchan = 167 # num...
//start = '-49.7km/s' # Velo...
//width = '0.6km/s' # Chan...
//veltype = 'radio' # ve...
//outframe = '' ...
//gridfunction = 'GJINC' # g...
//imsize = [60,80] ...
//cell = ['6arcsec','6arcsec'] ...
//phasecenter = 'J2000 08:59:05.78 -047.31.00.020' ...
//projection = 'SIN' # m...
//ephemsrcname = '' # epheme...
//pointingcolumn = 'direction' # pointing...
//restfreq = '' #...
//stokes = '' #...
//minweight = 0.1 # Min...
//brightnessunit = '' # O...
//clipminmax = False # Clip ...
+ 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']...
//mode = 'fft_mask' ...
//direction = [0,90] # scan...
//maskwidth = [5.0,3.0] # mask widt...
//tmax = 0.0 # maximu...
//tmin = 0.0 # minim...
//outfile = '○○○.int' # name of...
//overwrite = True # overwrite...
+ 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 in...
//fitsimage = '○○○.fits' # Name ...
//velocity = True # Use ...
//optical = False # Use ...
//bitpix = -32 # B...
//minpix = 0 # Mi...
//maxpix = -1 # Max...
//overwrite = True # Overwr...
//dropstokes = True # Drop the...
//stokeslast = True # Put Sto...
//history = True # Writ...
//dropdeg = False # Drop al...
** 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 sho...
** Acknowledgement and references in publications [#l3e2b...
*** Acknowledgement [#xd3de2b1]
- Please include the following statement in the acknowled...
#br
The ASTE demo science data collection was partially suppo...
*** References [#ra523fff]
- Please include the following paper for the reference to...
#br
Izumi, N. et al. 2021, PASJ, 73, 174
End:
[[ASTE Wiki]]
* Demo [CI](3P1-3P0) map of a star forming region RCW38 †...
: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.j...
** 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...
** CASA data reduction [#ta8cd472]
- Measurement Set
|CENTER: ||CENTER: ||c
|CENTER: Scan &br; direction |CENTER: File |CENTER: Numbe...
|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 Obs...
|~| R38otfy.RCW_38.20161126060513.4.ms | 16 | <= Start ob...
|~| R38otfY.RCW_38.20161126075453.4.ms | 61 ||
|~| R38otfY.RCW_38.20161126101442.4.ms | 61 ||
-- Amplitude calibration with the chopper wheel method is...
-- The unit of amplitude in all measurement sets are Ta*.
- Spectral window
| 0, 1 (A1, A2) | [CI] (same signal (polarization) input ...
| 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', ...
blfunc='poly', order=1, ...
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 data...
// datacolumn = 'float_data' # name of data column t...
// antenna = '' # select data by antenna name...
// field = '' # select data by field IDs an...
// spw = '0:350~470;560~650,1:350~470;560~650' ...
// timerange = '' # select data by time range, ...
// scan = '' # select data by scan numbers...
// pol = '' # select data by polarization...
// intent = '' # select data by observationa...
// maskmode = 'list' # mode of setting additional ...
// blmode = 'fit' # baselining mode ["fit" or "...
// dosubtract = True # subtract baseline from inpu...
// blformat = 'text' # format(s) of file(s) in whi...
// bloutput = '' # name(s) of file(s) in which...
// blfunc = 'poly' # baseline model function
// order = 1 # order of baseline model fun...
// clipthresh = 3.0 # clipping threshold for iter...
// clipniter = 0 # maximum iteration number fo...
// showprogress = False # (NOT SUPPORTED YET) show pr...
// outfile = 'XXX.bl.ms' # name of output file
// overwrite = False # overwrite the output file i...
+ 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...
caltype='amp',
parameter=[sqrt(0.45)]) # Main beam efficiency of 4...
# See the ASTE home page (h...
// Task : gencal (Specify Calibration Values of Various T...
// Parameters (examples)
// vis = 'XXX.bl.ms' # Name of input vi...
// caltable = 'XXX.bl.mb.tbl’ # The new/existing ...
// caltype = 'amp' # The calibration ...
// spw = '' # Calibration spw(...
// antenna = '' # Calibration ante...
// pol = '' # Calibration pola...
// parameter = [sqrt(0.45)] # The calibration...
++ Applying the scaling table
$ default(applycal)
$ applycal(
vis=‘R38otfX.RCW_38.20161126043008.4.A1A2.bl.ms',
gaintable=['R38otfX.RCW_38.20161126043008.4.A1A2.bl....
calwt=[False])
// Task : applycal (Apply calibrations solution(s) to data)
// Parameters (examples)
// vis = '○○○.bl.ms' # Name of input visi...
// field = '' # Select field using field id...
// spw = '' # Select spectral window/chan...
// intent = '' # Select observing intent
// selectdata = True # Other data selection parame...
// timerange = '' # Select data based on time r...
// uvrange = '' # Select data within uvrange ...
// antenna = '' # Select data based on antenn...
// scan = '' # Scan number range
// observation = '' # Select by observation ID(s)
// msselect = '' # Optional complex data selec...
// docallib = False # Use callib or traditional c...
// gaintable = [] # Gain calibration table(s) t...
// gainfield = [] # Select a subset of calibrat...
// interp = [] # Interp type in time[,freq],...
// spwmap = [] # Spectral windows combinatio...
// calwt = [False] # Calibrate data weights per ...
// parang = False # Apply parallactic angle cor...
// applymode = '' # Calibration mode: ""="calfl...
// flagbackup = True # Automatically back up the s...
+ 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=...
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 ...
//Parameters (examples)
//nfiles = ['○○○.bl.ms','×××.bl.ms',......
//outfile = '○○○.int' #...
//overwrite = True # o...
//field = '' ...
//spw = '' ...
//antenna = '' ...
//scan = '' ...
//intent = '*ON_SOURCE*' # select d...
//mode = 'velocity' # ...
//nchan = 167 # num...
//start = '-49.7km/s' # Velo...
//width = '0.6km/s' # Chan...
//veltype = 'radio' # ve...
//outframe = '' ...
//gridfunction = 'GJINC' # g...
//imsize = [60,80] ...
//cell = ['6arcsec','6arcsec'] ...
//phasecenter = 'J2000 08:59:05.78 -047.31.00.020' ...
//projection = 'SIN' # m...
//ephemsrcname = '' # epheme...
//pointingcolumn = 'direction' # pointing...
//restfreq = '' #...
//stokes = '' #...
//minweight = 0.1 # Min...
//brightnessunit = '' # O...
//clipminmax = False # Clip ...
+ 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']...
//mode = 'fft_mask' ...
//direction = [0,90] # scan...
//maskwidth = [5.0,3.0] # mask widt...
//tmax = 0.0 # maximu...
//tmin = 0.0 # minim...
//outfile = '○○○.int' # name of...
//overwrite = True # overwrite...
+ 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 in...
//fitsimage = '○○○.fits' # Name ...
//velocity = True # Use ...
//optical = False # Use ...
//bitpix = -32 # B...
//minpix = 0 # Mi...
//maxpix = -1 # Max...
//overwrite = True # Overwr...
//dropstokes = True # Drop the...
//stokeslast = True # Put Sto...
//history = True # Writ...
//dropdeg = False # Drop al...
** 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 sho...
** Acknowledgement and references in publications [#l3e2b...
*** Acknowledgement [#xd3de2b1]
- Please include the following statement in the acknowled...
#br
The ASTE demo science data collection was partially suppo...
*** References [#ra523fff]
- Please include the following paper for the reference to...
#br
Izumi, N. et al. 2021, PASJ, 73, 174
Page:
![[PukiWiki]](http://www2.nao.ac.jp/~aste-www/aste_wiki/obs.jpg "[PukiWiki]")
