Array combination for CS(2-1) (12m + 7m + TP)

[pyhsiehATalma]

This is general issue for tracking the 12m+7m+TP combination for the CS(2-1) line data.

• 12m: spw 33
• 7m: spw 24
• tp: spw 25

Update: 2024 Jan 31 12m+7m+TP 41 fields combination.

Proposal ID | Obs ID | 7m MOUS ID | 12m MOUS ID | issue | Note -- | -- | -- | -- | -- | -- cmz01 | a | X1590_X30ac | X1590_X30aa | #242 |   cmz02 | b | X15b4_X37 | X15a0_Xa6 | #150 | imcontsub was applied to CS2-1 cube (-300~300 km/s), the continuum subtration looks fine. cmz03 | c | X15a0_Xae | X15a0_Xac | #256 |   cmz04 | d | X15a0_Xb4 | X15a0_Xb2 | #247 |   cmz05 | e | X15b4_X39 | X15a0_Xb8 | #147 |   cmz06 | f | X15a0_Xc0 | X15a0_Xbe | #266 | imcontsub was applied to CS2-1 cube (-300~300 km/s), the continuum subtration looks fine. cmz07 | g | X15a0_Xc6 | X15a0_Xc4 | #134 |   cmz08 | h | X15a0_Xcc | X15a0_Xca | #142 |   cmz09 | i | X15b4_Xc5 | X15a0_Xd0 | #138 |   cmz10 | j | X15a0_Xd8 | X15b4_X3d | #117 |   cmz11 | k | X15a0_Xde | X15a0_Xdc | #152 |   cmz12 | l | X15a0_Xe4 | X15a0_Xe2 | #137 |   cmz13 | m | X15b4_X43 | X15b4_X41 | #118 | Sgr B2, imcontsub was applied to 12m cmz14 | n | X15a0_Xf0 | X15a0_Xee | #231 | size mitigated, need re-run pipeline cmz15 | o | X15a0_Xf6 | X15a0_Xf4 | 272 |   cmz16 | p | X15a0_Xfc | X15a0_Xfa | #149 |   cmz17 | q | X15a0_X102 | X15a0_X100 | #109 |   cmz18 | r | X15a0_X108 | X15a0_X106 | #263 | imcontsub was applied to CS2-1 cube (-300~300 km/s), the continuum subtration looks fine. cmz19 | s | X15a0_X10e | X15a0_X10c | #264 | divergence, i am doing re-imaing with full cube. cmz20 | t | X15a0_X114 | X15a0_X112 | #113 |   cmz21 | u | X15a0_X11a | X15a0_X118 | #241 |   cmz22 | v | X15a0_X120 | X15a0_X11e | #235 |   cmz23 | w | X15a0_X126 | X15a0_X124 | #124 |   cmz24 | x | X15a0_X12c | X15a0_X12a | #125 |   cmz25 | y | X15a0_X132 | X15a0_X130 | #40 |   cmz26 | z | X15a9_X12d3 | X15a0_X136 | #131 |   cmz27 | aa | X15a0_X13e | X15a0_X13c | #251 | imcontsub was applied to full spw 33, the continuum subtration looks fine. cmz28 | ab | X15a0_X144 | X15a0_X142 | #129 |   cmz29 | ac | X15a0_X14a | X15a0_X148 | #232 |   cmz30 | ad | X15a0_X150 | X15a0_X14e | #28 |   cmz31 | ae | X15a0_X156 | X15a0_X154 | #132 |   cmz32 | af | X15a0_X15c | X15a0_X15a | #255 |   cmz33 | ag | X15a0_X162 | X15a0_X160 | #136 |   cmz34 | ah | X15a0_X168 | X15a0_X166 | #238 |   cmz35 | ai | X15a0_X16e | X15a0_X16c | #130 |   cmz36 | aj | X15b4_X47 | X15a0_X172 | #123 |   cmz37 | ak | X15a9_X12d5 | X15b4_X4b | #120 |   cmz38 | al | X15a0_X180 | X15a0_X17e | #126 |   cmz39 | am | X15a0_X186 | X15a0_X184 | #259 | divergence, i am doing re-imaing with CS(2-1) cube (-300~300 km/s) cmz40 | an | X15a0_X18c | X15a0_X18a | #133 |   cmz41 | ao | X15a0_X192 | X15a0_X190 | #41 | 12m divergence at channel 645, low level and at edge, mayby ok cmz42 | ap | X15a0_X198 | X15a0_X196 | #43 |   cmz43 | aq | X15a0_X19e | X15a0_X19c | #140 | Sgr A*, divergence, imcontsub was applied to 12m, i am also doing re-imaing with CS(2-1) cube (-300~300 km/s) cmz44 | ar | X15b4_X4f | X15a0_X1a2 | #114 |   cmz45 | as | X15a0_X1aa | X15a0_X1a8 | #145 |   Jan 18 7m+TP 45m fields combination. Files are under /upload/Mosaic-CS2-1/ACA_TP_Combination**, - Combined cubes of individual fields: /Individual Cubes - MOM8 of individual fields: /Individual MOM8 - MOM0 of individual fields: /Individual MOM0 - Mosaic combined cube and moments: /Mosaic Note that the individual combined cubes were created using the full cube, while the mosaic combined cube was generated specifically for the CS(2-1) line with a velocity range of -500 to +500 km/s. Proposal ID | Obs ID | Name | vel. offset | 7m MOUS ID | issue | Divergence | Using Product | Residual Cont. | TP MOUS ID | issue -- | -- | -- | -- | -- | -- | -- | -- | -- | -- | -- cmz01 | a | - | 40 | X1590_X30ac | #214 | no |   |   | X1590_X30ae | #215 cmz02 | b | Sgr B2 | 40 | X15b4_X37 | #162 | no |   | yes | X15a0_Xaa | #252 cmz03 | c | - | 0 | X15a0_Xae | #14 | yes | yes |   | X15a0_Xb0 | #175 cmz04 | d | - | 30 | X15a0_Xb4 | #228 | ch=597-602, 613 | yes |   | X15a0_Xb6 | #227 cmz05 | e | - | -30 | X15b4_X39 | #168 | no |   |   | X15a0_Xbc | #285 cmz06 | f | - | 0 | X15a0_Xc0 | #254 | no |   | yes | X15a0_Xc2 | #315 cmz07 | g | - | 30 | X15a0_Xc6 | #161 | no |   |   | X15a0_Xc8 | #197 cmz08 | h | - | 0 | X15a0_Xcc | #151 | no |   |   | X15a0_Xce | #177 cmz09 | i | - | 30 | X15b4_Xc5 | #191 | no |   | yes | X15a0_Xd4 | #187 cmz10 | j | - | 30 | X15a0_Xd8 | #181 | no |   |   | X15b4_X3f | #324 cmz11 | k | - | 30 | X15a0_Xde | #182 | no |   |   | X15a0_Xe0 | #193 cmz12 | l | - | 0 | X15a0_Xe4 | #153 | no |   |   | X15a0_Xe6 | #224 cmz13 | m | - | 30 | X15b4_X43 | #110 | no |   | yes | X15b4_X45 | #170 cmz14 | n | - | 30 | X15a0_Xf0 | #225 | no |   |   | X15a0_Xf2 | #226 cmz15 | o | - | 0 | X15a0_Xf6 | #257 | no |   | yes | X15a0_Xf8 | #278 cmz16 | p | - | 0 | X15a0_Xfc | #148 | no |   | yes | X15a0_Xfe | #183 cmz17 | q | - | 0 | X15a0_X102 | #171 | no |   |   | X15a0_X104 | #284 cmz18 | r | - | 30 | X15a0_X108 | #240 | no |   | yes | X15a0_X10a | #246 cmz19 | s | - | 0 | X15a0_X10e | #243 | no |   |   | X15a0_X110 | #262 cmz20 | t | - | 0 | X15a0_X114 | #8 | ch=701-703 |   |   | X15a0_X116 | #172 cmz21 | u | - | 30 | X15a0_X11a | #216 | no |   |   | X15a0_X11c | #217 cmz22 | v | - | 30 | X15a0_X120 | #218 | no |   | yes | X15a0_X122 | #219 cmz23 | w | - | 30 | X15a0_X126 | #144 | no |   |   | X15a0_X128 | #173 cmz24 | x | - | 0 | X15a0_X12c | #135 | no |   |   | X15a0_X12e | #164 cmz25 | y | - | 0 | X15a0_X132 | #4 | no |   |   | X15a0_X134 | #27 cmz26 | z | - | 40 | X15a9_X12d3 | #115 | no |   |   | X15a0_X13a | #236 cmz27 | aa | - | 0 | X15a0_X13e | #237 | no |   | yes | X15a0_X140 | #229 cmz28 | ab | - | 0 | X15a0_X144 | #141 | no |   |   | X15a0_X146 | #167 cmz29 | ac | - | 0 | X15a0_X14a | #220 | no |   |   | X15a0_X14c | #221 cmz30 | ad | - | 0 | X15a0_X150 | #5 | no |   |   | X15a0_X152 | #25 cmz31 | ae | - | 0 | X15a0_X156 | #186 | ch=652,653,655 | yes |   | X15a0_X158 | #180 cmz32 | af | - | 0 | X15a0_X15c | #244 | no |   | yes | X15a0_X15e | #271 cmz33 | ag | - | 0 | X15a0_X162 | #174 | no |   |   | X15a0_X164 | #190 cmz34 | ah | - | 30 | X15a0_X168 | #230 | ch=596, 615, 617 | yes |   | X15a0_X16a | #233 cmz35 | ai | - | 0 | X15a0_X16e | #143 | no |   |   | X15a0_X170 | #169 cmz36 | aj | - | -30 | X15b4_X47 | #119 | no |   |   | X15b4_X49 | #239 cmz37 | ak | - | 30 | X15a9_X12d5 | #116 | no |   |   | X15b4_X4d | #121 cmz38 | al | - | 0 | X15a0_X180 | #139 | no |   |   | X15a0_X182 | #166 cmz39 | am | - | 30 | X15a0_X186 | #265 | no |   | yes | X15a0_X188 | #245 cmz40 | an | - | 30 | X15a0_X18c | #222 | no |   |   | X15a0_X18e | #184 cmz41 | ao | Brick | 0 | X15a0_X192 | #6 | ch=631 | need re-imaging |   | X15a0_X194 | #29 cmz42 | ap | - | 0 | X15a0_X198 | #158 | no |   |   | X3577_X6c5 | #325 cmz43 | aq | Sgr A* | 0 | X15a0_X19e | #33 | no |   | yes | X15a0_X1a0 | #223 cmz44 | ar | - | 40 | X15b4_X4f | #159 | no |   |   | X15a0_X1a6 | #317 cmz45 | as | - | 0 | X15a0_X1aa | #165 | yes, ch=616 | yes |   | X15a0_X1ac | #322 **2023 July 18/December 28:** 12m+7m+tp 37 fields mom8 and 7m+tp 44 fields mom8 maps. note that the 0, 1, 2... is not the field number. 37 Fields Proposal ID | Obs ID |   | 7m MOUS ID | 12m MOUS ID |   | TP MOUS ID -- | -- | -- | -- | -- | -- | -- cmz01 | a | 0 | X1590_X30ac | X1590_X30aa |   | X1590_X30ae cmz02 | b | 1 | X15b4_X37 | X15a0_Xa6 |   | X15a0_Xaa cmz04 | d | 2 | X15a0_Xb4 | X15a0_Xb2 |   | X15a0_Xb6 cmz05 | e | 3 | X15b4_X39 | X15a0_Xb8 |   | X15a0_Xbc cmz06 | f | 4 | X15a0_Xc0 | X15a0_Xbe |   | X15a0_Xc2 cmz07 | g | 5 | X15a0_Xc6 | X15a0_Xc4 |   | X15a0_Xc8 cmz08 | h | 6 | X15a0_Xcc | X15a0_Xca |   | X15a0_Xce cmz09 | i | 7 | X15b4_Xc5 | X15a0_Xd0 |   | X15a0_Xd4 cmz10 | j | 8 | X15a0_Xd8 | X15b4_X3d |   | X15b4_X3f cmz11 | k | 9 | X15a0_Xde | X15a0_Xdc |   | X15a0_Xe0 cmz12 | l | 10 | X15a0_Xe4 | X15a0_Xe2 |   | X15a0_Xe6 cmz14 | n | 11 | X15a0_Xf0 | X15a0_Xee |   | X15a0_Xf2 cmz15 | o | 12 | X15a0_Xf6 | X15a0_Xf4 |   | X15a0_Xf8 cmz16 | p | 13 | X15a0_Xfc | X15a0_Xfa |   | X15a0_Xfe cmz17 | q | 14 | X15a0_X102 | X15a0_X100 |   | X15a0_X104 cmz21 | u | 15 | X15a0_X11a | X15a0_X118 |   | X15a0_X11c cmz22 | v | 16 | X15a0_X120 | X15a0_X11e |   | X15a0_X122 cmz23 | w | 17 | X15a0_X126 | X15a0_X124 |   | X15a0_X128 cmz24 | x | 18 | X15a0_X12c | X15a0_X12a |   | X15a0_X12e cmz25 | y | 19 | X15a0_X132 | X15a0_X130 |   | X15a0_X134 cmz26 | z | 20 | X15a9_X12d3 | X15a0_X136 |   | X15a0_X13a cmz27 | aa | 21 | X15a0_X13e | X15a0_X13c |   | X15a0_X140 cmz28 | ab | 22 | X15a0_X144 | X15a0_X142 |   | X15a0_X146 cmz29 | ac | 23 | X15a0_X14a | X15a0_X148 |   | X15a0_X14c cmz30 | ad | 24 | X15a0_X150 | X15a0_X14e |   | X15a0_X152 cmz31 | ae | 25 | X15a0_X156 | X15a0_X154 |   | X15a0_X158 cmz33 | ag | 26 | X15a0_X162 | X15a0_X160 |   | X15a0_X164 cmz34 | ah | 27 | X15a0_X168 | X15a0_X166 |   | X15a0_X16a cmz35 | ai | 28 | X15a0_X16e | X15a0_X16c |   | X15a0_X170 cmz36 | aj | 29 | X15b4_X47 | X15a0_X172 |   | X15b4_X49 cmz37 | ak | 30 | X15a9_X12d5 | X15b4_X4b |   | X15b4_X4d cmz38 | al | 31 | X15a0_X180 | X15a0_X17e |   | X15a0_X182 cmz41 | ao | 32 | X15a0_X192 | X15a0_X190 |   | X15a0_X194 cmz42 | ap | 33 | X15a0_X198 | X15a0_X196 |   | X3577_X6c5 cmz43 | aq | 34 | X15a0_X19e | X15a0_X19c |   | X15a0_X1a0 cmz44 | ar | 35 | X15b4_X4f | X15a0_X1a2 |   | X15a0_X1a6 cmz45 | as | 36 | X15a0_X1aa | X15a0_X1a8 |   | X15a0_X1ac

[pyhsiehATalma]

The testing script is to feather 12m data with existing 7m+TP cubes (made in another script). The next step is to combine these two scripts. Basically they are doing the same thing. Note that primary beam corrected files (*.pbcor) should be used. The maps I show below used non-primary corrected one (just an error).

intensity peak maps (moment 8) of these 25 fields are shown below.

[pyhsiehATalma]

Mosaic intensity peak map of 25 fields is shown below. Upper panel shows 12m+7m+TP, lower panel shows 7m+TP.

Zoomed-in maps are shown below. Left panel shows 12m+7m+TP, right panel shows 7m+TP.

[pyhsiehATalma]

just a note: recording hardisk memory required and time for calculation.

[pyhsiehATalma]

Trouble shooting 1) If you see this error when running casa task feather, it means one of image has multiple beam information in the header. Sometimes it happens, and the solution is to use casa task imsmooth to restore the image to a common beam. imsmooth(imagename='input.image', kernel='c',outfile='out.image')

2023-07-20 04:44:12 SEVERE feather::imager::feather() Caught exception: (casacore/images/Images/ImageBeamSet.cc : 115) Failed AlwaysAssert chan >=0 && chan < Int(nchan()) && stokes >= 0 && stokes < Int(nstokes()) 2023-07-20 04:44:12 SEVERE feather::::casa Task feather raised an exception of class RuntimeError with the following message: 2023-07-20 04:44:12 SEVERE feather::imager::feather() Caught exception: (casacore/images/Images/ImageBeamSet.cc : 115) Failed AlwaysAssert chan >=0 && chan < Int(nchan()) && stokes >= 0 && stokes < Int(nstokes()) 2023-07-20 04:44:12 SEVERE feather::::casa Exception Reported: Error in feather: 2023-07-20 04:44:12 SEVERE feather::imager::feather() Caught exception: (casacore/images/Images/ImageBeamSet.cc : 115) Failed AlwaysAssert chan >=0 && chan < Int(nchan()) && stokes >= 0 && stokes < Int(nstokes()) 962.793-2023-07-20 05:25:02 SEVERE feather::imager::feather() Caught exception: (casacore/images/Images/ImageBeamSet.cc : 115) Failed AlwaysAssert chan >=0 && chan < Int(nchan()) && stokes >= 0 && stokes < Int(nstokes())

2) axes orders of two images don't match (I will copy the error messages later). Sometimes the axes orders of input images are different due to different casa version, so feather also fails. The solution is to use task imtrans to change the orders.

for example, imtrans(imagename='input.image',order='0132',out='output.image')

There are more errors, TBC.

[pyhsiehATalma]

12m+7m+tp 37 fields mom8 and 7m+tp 44 fields mom8 maps.

[pyhsiehATalma]

2024 Jan 17 45 Fields ACA+TP combined maps. Upper: CS(2-1) MOM0 Lower: CS(2-1) MOM8 Note: The divergence in most of the fields has been resolved; for those that could not be resolved, product images (with no divergence) were used for combination. Therefore, this 45-field map does not have divergence, except for channel 637 in field ao (X15a0_X192), which shows a slight divergence.

Residual of Continuum Emission

13 Fields show continuum residual. Below I show the maps integrated from -500km/s - -400 km/s (line-free region in the CS(2-1) ACA+TP mosaic cube). Field aq, am, af, aa, o, p, v, x, m, r, i, f, b. Lower: residual of continuum emission Upper: residual of continuum emission (contours) overlaid on the CS(2-1) mom0 map. The continuum residual is weaker than CS(2-1) map.

Update of 12m+ACA+TP combined maps. TBC...

[pyhsiehATalma]

2024 Jan 24 12m+7m+TP combination, problems of the 12m data.

Missing in Mosaic | Obs ID | 12m MOUS ID | Notes | issue | Assigned To |   -- | -- | -- | -- | -- | -- | -- no included in mosaic | m | X15b4_X41 | yes and contsub issue | #118 | Nazar | Divergence fixed, constub issue remains no included in mosaic | s | X15a0_X10c | divergence, Pei-Ying is re-imaging. | #264 | Qizhou / Pei Ying (in progress) |     | aa | X15a0_X13c | divergence was fixed, but continuum is not subtracted? | #251 | Katha | imcontsub over full spw no included in mosaic | aq | X15a0_X19c | divergence at channel 638 | #140 | Dan (in progress) |   no included in mosaic | am | X15a0_X184 | there are two cubes | #259 | Dan (done) |     |   |   | Sgr_A_star_sci33.cube.I.manual.pbcor.fits |   |   | only 1917 channels, no divergence   |   |   | uid___A001_X15a0_X184.Sgr_A_star_sci33.cube.I.manual.image.pbcor |   |   | divergence around ch=626   | n | X15a0_Xee | bandwidth=1917 channel, need re-run pipeline | #231 |   |     | ad | X15a0_X14e | 5" beam | #28 |   | 2" is doable for narrow bandwidth   | l | X15a0_Xe2 | multibeam | #137 |   |     | ao | X15a0_X190 | divergence at channel 645 | #41 |   |     | f | X15a0_Xbe | contsub issue? | #266 |   | imcontsub over narrow bandwidth   | r | X15a0_X106 | contsub issue? | #263 |   | imcontsub over narrow bandwidth   | b | X15a0_Xa6 | absorption? mask may be needed for negative emission | #150 |   |  

[pyhsiehATalma]

moment 8 maps of 41 EB 12m+7m+TP combined maps.

[pyhsiehATalma]

I started working on mosaic of 45 fields, but note that final beam size will be 3". This is because several fields - X15a, X184, Xac have larger beam sizes (2.758"x1.964", 2.837"x1.92", 2.526"x1.738").

Also, I need to subtract the residual of the continuum on 7m and 12m if there is any, but I will leave this task for next week. Mosaic first!

I am also re-imaging on,

X15a0_X10c: full band X15a0_X19c: narrow bandwidth X15a0_X184: narrow bandwidth

[pyhsiehATalma]

Full CMZ 12m+7m+TP map. Upper: moment 0. Lower: moment 8. bandwidth of mosaic CS(2-1) cube: -250 km/s - 250 km/s beam: 3"

Just my personal idea: those black dots seem natural, likely due to foreground contamination, especialy in band 3.

Remaining Issues Proposal ID	Obs ID	7m MOUS ID	12m MOUS ID	issue	Iusse	For the purpose of combination
cmz02	b	X15b4_X37	X15a0_Xa6	#150	residual of continuum emission	I applied imcontsub to the CS2-1 cube (-300~300 km/s), the continuum subtration looks fine.
cmz06	f	X15a0_Xc0	X15a0_Xbe	#266	residual of continuum emission	I applied imcontsub to the CS2-1 cube, the continuum subtration looks fine.
cmz13	m	X15b4_X43	X15b4_X41	#118	residual of continuum emission	Sgr B2, I applied imcontsub to the CS2-1 cube (-300~300 km/s), slightly improved.
cmz14	n	X15a0_Xf0	X15a0_Xee	#231	size mitigation	need re-run pipeline
cmz18	r	X15a0_X108	X15a0_X106	#263	residual of continuum emission	I applied imcontsub to the CS2-1 cube (-300~300 km/s), the continuum subtration looks fine.
cmz19	s	X15a0_X10e	X15a0_X10c	#264	divergence	I am doing re-imaing with full cube, parallel tclean is ongoing
cmz27	aa	X15a0_X13e	X15a0_X13c	#251	continuum was not subtracted in the cube	I applied imcontsub to the full spw 33, the continuum subtration looks fine.
cmz39	am	X15a0_X186	X15a0_X184	#259	divergence	I am doing re-imaing with the CS(2-1) cube (-300~300 km/s)
cmz41	ao	X15a0_X192	X15a0_X190	#41	low-level divergence	no action item
cmz43	aq	X15a0_X19e	X15a0_X19c	#140	divergence and residual of continuum	Sgr A*, I applied imcontsub to 12m, i am also doing re-imaing with CS(2-1) cube (-300~300 km/s)