#------------------------------------------------ # IC2574, Raster, Night 1, 1st February 2008 run #------------------------------------------------ ------------------------------------------------------------------------------- INITIAL REDUCTION USING NEAREST SPREAD: makecube in="a20080211_000*_flagspike.sdf" autogrid=true crota=40 pixsize=7.2761 sparse=false spread='nearest' specbounds=\"-370.0 440.0\" out=ic2574_prelim20080211_wholenight reset # Projection parameters used: # CRPIX1 = 46.2020408213139 # CRPIX2 = 29.8653061240911 # CRVAL1 = 157.099583333333 ( RA = 10:28:23.9 ) # CRVAL2 = 68.4122222222222 ( Dec = 68:24:44 ) # CDELT1 = -0.00202113893296984 ( -7.3 arcsec ) # CDELT2 = 0.00202113893296984 ( 7.3 arcsec ) # CROTA2 = 40 # # Output cube pixel bounds: ( -44:46, -28:30, -957:956 ) # Output cube WCS bounds: # Right ascension: 10:27:13.1 -> 10:29:35.0 # Declination: 68:18:28 -> 68:31:03 # Optical velocity (USB): -370.0955 -> 440.3164 km/s mfittrend ic2574_prelim20080211_wholenight ranges=\"-370 -70 130 440\" order=3 axis=3 out=ic2574_prelim20080211_wholenight_base3 sub ic2574_prelim20080211_wholenight ic2574_prelim20080211_wholenight_base3 ic2574_prelim20080211_wholenight_sub3 compave ic2574_prelim20080211_wholenight_sub3 ic2574_prelim20080211_wholenight_sub3_20kms compress="[1,1,47]" stats 'ic2574_prelim20080211_wholenight_sub3_20kms(~35,~35,:)' | grep Standard # Standard deviation : 0.0159319 stats 'ic2574_prelim20080211_wholenight_sub3_20kms(~35,~35,-370.0:-70.0)' | grep Standard # Standard deviation : 0.0157069 stats 'ic2574_prelim20080211_wholenight_sub3_20kms(~35,~35,130.0:440.0)' | grep Standard # Standard deviation : 0.0141058 # We are well under our target rms after only three MSB repeats, so the remaining MSBs are not required. No CO detection seen yet. ------------------------------------------------------------------------------- NOW SINCSINC: makecube in="a*" autogrid=true crota=40 pixsize=7.2761 sparse=false spread='sincsinc' params='[1,2]' specbounds=\"-370.0 440.0\" out=ic2574_sincsinc reset # Projection parameters used: # CRPIX1 = 46.4489795863628 # CRPIX2 = 30.067346945405 # CRVAL1 = 157.097916666667 ( RA = 10:28:23.5 ) # CRVAL2 = 68.4122222222222 ( Dec = 68:24:44 ) # CDELT1 = -0.00202113893296984 ( -7.3 arcsec ) # CDELT2 = 0.00202113893296984 ( 7.3 arcsec ) # CROTA2 = 40 # # Output cube pixel bounds: ( -44:46, -28:30, -957:956 ) # Output cube WCS bounds: # Right ascension: 10:27:13.1 -> 10:29:35.0 # Declination: 68:18:28 -> 68:31:03 # Optical velocity (USB): -370.0955 -> 440.3164 km/s # # # WARNING: Cannot create output Tsys values. chpix in=ic2574_sincsinc out=ic2574_sincsinc_spikebad section=\":,:,-444:-431\" newval=bad ndfcopy in=ic2574_sincsinc_spikebad'(-36:29,-11:14,)' out=ic2574_sincsinc_trim mfittrend ic2574_sincsinc_trim ranges=\"-370 -70 130 440\" order=3 axis=3 out=ic2574_sincsinc_trim_base3 # Fitting NDF Axis 3, using ranges in current co-ordinate Frame (pixels): # -369.6732 : -70.40656 km/s (-957 : -249) # 129.8608 : 439.8919 km/s (224 : 955) sub ic2574_sincsinc_trim ic2574_sincsinc_trim_base3 ic2574_sincsinc_trim_sub3 compave ic2574_sincsinc_trim_sub3 ic2574_sincsinc_trim_sub3_20kms compress="[1,1,47]" stats 'ic2574_sincsinc_trim_sub3_20kms(~35,~35,:)' | grep Standard # Standard deviation : 0.0145549 stats 'ic2574_sincsinc_trim_sub3_20kms(~35,~35,-370.0:-70.0)' | grep Standard # Standard deviation : 0.0142807 stats 'ic2574_sincsinc_trim_sub3_20kms(~35,~35,130.0:440.0)' | grep Standard # Standard deviation : 0.0128754 ------------------------------------------------------------------------------- BASELINE MASK PRODUCTION: # First produce a smoothed cube to identify the line and make masks with mfittrend block in=ic2574_sincsinc_trim out=ic2574_sincsinc_trim_vsmooth 'axes=[1,3] box=[1,25]' block in=ic2574_sincsinc_trim_vsmooth out=ic2574_sincsinc_trim_smooth 'axes=[1,2] box=[3,3]' mfittrend in=ic2574_sincsinc_trim_smooth axis=3 order=3 out=ic2574_sincsinc_trim_smooth_base3 auto mask=ic2574_sincsinc_trim_smooth_basemsk3 method=single variance=false subtract=false 'clip=[2,3]' mfittrend in=ic2574_sincsinc_trim_smooth axis=3 order=3 out=ic2574_sincsinc_trim_smooth_base3g auto mask=ic2574_sincsinc_trim_smooth_basemsk3g method=global variance=false subtract=false 'clip=[2,2.5,3]' add ic2574_sincsinc_trim_smooth_basemsk3 ic2574_sincsinc_trim_smooth_basemsk3g ic2574_trim_basemask0 chpix in=ic2574_trim_basemask0 out=ic2574_trim_basemask0_end1 section=\":,:,-957:-249\" newval=0 chpix in=ic2574_trim_basemask0_end1 out=ic2574_trim_basemask0_end2 section=\":,:,224:955\" newval=0 cadd ic2574_trim_basemask0_end2 1 ic2574_trim_basemask1 mult ic2574_sincsinc_trim ic2574_trim_basemask1 ic2574_sincsinc_trim_masked mfittrend ic2574_sincsinc_trim_masked ranges=\"-370 440\" order=3 axis=3 out=ic2574_sincsinc_trim_masked_base3 # Fitting NDF Axis 3, using ranges in current co-ordinate Frame (pixels): # -369.6732 : 439.8919 km/s (-957 : 955) sub ic2574_sincsinc_trim ic2574_sincsinc_trim_masked_base3 ic2574_sincsinc_trim_submask3 compave ic2574_sincsinc_trim_submask3 ic2574_sincsinc_trim_submask3_20kms compress="[1,1,47]" stats 'ic2574_sincsinc_trim_submask3_20kms(~35,~35,:)' | grep Standard # Standard deviation : 0.013902 stats 'ic2574_sincsinc_trim_submask3_20kms(~35,~35,-370.0:-70.0)' | grep Standard # Standard deviation : 0.0146078 stats 'ic2574_sincsinc_trim_submask3_20kms(~35,~35,130.0:440.0)' | grep Standard # Standard deviation : 0.0130904 ------------------------------------------------------------------------------- CREATE MOMENT MAPS: mkdir Moment_sincsinc cp ic2574_sincsinc_trim_submask3.sdf Moment_sincsinc/ic2574_sincsinc_trim_submask3_rembsl.sdf cp ic2574_trim_basemask1.sdf Moment_sincsinc/ic2574_sincsinc_trim_submask3_msk.sdf cp ../NGC3031/Moment_sincsinc/REDUCE_SCIENCE_GRADIENT Moment_sincsinc/ cd Moment_sincsinc *Edit RSG Script: emacs REDUCE_SCIENCE_GRADIENT & set object = "ic2574_sincsinc_trim_submask3" set vrange = "low=-70 high=130" # Quick jumps to repeat parts of the reduction #goto STEP1a #goto STEP1b #goto STEP2 goto STEP3 foreach rmslev (2.0 2.5 3.0 3.5 4.0 4.5 5.0) **Remove trim function: # Multiply baselined cube with clump mask ${kap}/mult ${sdf} ${sdf}_findclump${rmslab} ${sdf}_clmsk${rmslab} # Trim the noisy edges # ${kap}/ndfcopy temp'(-26:26,-36:40,)' ${sdf}_clmsk${rmslab} \rm -f ${sdf}_mask_temp1.sdf ${sdf}_mask_temp2.sdf >& /dev/null ./REDUCE_SCIENCE_GRADIENT # RMS MASKED SMOOTHED CUBE(-24:17,-6:9,-899:898): 0.00710530998471045 Rotate intvel maps: rotate ic2574_sincsinc_trim_submask3_intvel20 ic2574_sincsinc_trim_submask3_intvel20_rot angle=! rotate ic2574_sincsinc_trim_submask3_intvel25 ic2574_sincsinc_trim_submask3_intvel25_rot angle=! rotate ic2574_sincsinc_trim_submask3_intvel30 ic2574_sincsinc_trim_submask3_intvel30_rot angle=! rotate ic2574_sincsinc_trim_submask3_intvel35 ic2574_sincsinc_trim_submask3_intvel35_rot angle=! rotate ic2574_sincsinc_trim_submask3_intvel40 ic2574_sincsinc_trim_submask3_intvel40_rot angle=! rotate ic2574_sincsinc_trim_submask3_intvel45 ic2574_sincsinc_trim_submask3_intvel45_rot angle=! rotate ic2574_sincsinc_trim_submask3_intvel50 ic2574_sincsinc_trim_submask3_intvel50_rot angle=! (contoured totint maps and velocity maps then made using gaia) ------------------------------------------------------------------------------- MAKE CUBES DISPLAYING S/N RATIO: collapse in=ic2574_sincsinc_trim_submask3.sdf low=-370 high=-70 axis=vopt out=sigma_low estimator=sigma reset collapse in=ic2574_sincsinc_trim_submask3.sdf low=130 high=440 axis=vopt out=sigma_high estimator=sigma reset block ic2574_sincsinc_trim_submask3_20kms.sdf ic2574_sincsinc_trim_submask3_20kms_bin5.sdf axes="[2,3]" box="[1,5]" estimator=mean block ic2574_sincsinc_trim_submask3_20kms.sdf ic2574_sincsinc_trim_submask3_20kms_bin9.sdf axes="[2,3]" box="[1,9]" estimator=mean manic in=sigma_low out=sigma_low_3d axes="[1,2,0]" lbound=-19 ubound=40 reset wcscopy ndf=sigma_low_3d like=ic2574_sincsinc_trim_submask3_20kms_bin5.sdf confirm=false reset cdiv sigma_low_3d 6.85565 sigma_low_20kms div ic2574_sincsinc_trim_submask3_20kms_bin5.sdf sigma_low_20kms.sdf ic2574_sincsinc_trim_submask3_20kms_bin5_SN.sdf div ic2574_sincsinc_trim_submask3_20kms_bin9.sdf sigma_low_20kms.sdf ic2574_sincsinc_trim_submask3_20kms_bin9_SN.sdf ------------------------------------------------------------------------------- MAKE TOTINT20 NOISE and SN MAPS: mkdir noise_maps cd noise_maps cp ../sigma_high.sdf . cp ../sigma_low.sdf . cp ../Moment_sincsinc/ic2574_sincsinc_trim_submask3_rembsl_findclump20.sdf . cp ../Moment_sincsinc/ic2574_sincsinc_trim_submask3_totint20.sdf . add sigma_high sigma_low out=sigma_sum cdiv sigma_sum.sdf 2 sigma_avg collapse in= ic2574_sincsinc_trim_submask3_rembsl_findclump20.sdf out=test0 axis=vopt estimator=sum wlim=0.0 (must make mask to multiply by collapsed image so values are real, not integer, otherwise sqrt function below will round to nearest integer) cmult sigma_avg 0 test1 cadd test1 1.0 test2 mult test2 test0 test3 (now have same collapsed image with real, not integer values) maths exp="sqrt(IA)" out=test4 IA=test3 cmult test4 0.42 test5 mult sigma_avg test5 ic2574_totint20_noise div ic2574_sincsinc_trim_submask3_totint20.sdf ic2574_totint20_noise ic2574_totint20_SN (contoured SN maps then made using gaia) ------------------------------------------------------------------------------- TRIM OUT NOISE FROM TOTINT20 MAP: *Try trimming regions with SN < 2.0: thresh ic2574_totint20_SN.sdf thrlo=2 newlo=bad thrhi=500 newhi=bad out=ic2574_totint20_SN_thresh2 # There were 398 elements changed in the DATA array below the threshold. # There were 0 elements changed in the DATA array above the threshold. cmult ic2574_totint20_SN_thresh2.sdf 0 ic2574_zeroall_2 cadd ic2574_zeroall_2 1 ic2574_msk_2 mult ic2574_msk_2.sdf ic2574_sincsinc_trim_submask3_totint20.sdf ic2574_totint20_noisetrim2 *Try trimming regions with SN < 1.5: thresh ic2574_totint20_SN.sdf thrlo=1.5 newlo=bad thrhi=500 newhi=bad out=ic2574_totint20_SN_thresh15 # There were 357 elements changed in the DATA array below the threshold. # There were 0 elements changed in the DATA array above the threshold. cmult ic2574_totint20_SN_thresh15.sdf 0 ic2574_zeroall_15 cadd ic2574_zeroall_15.sdf 1 ic2574_msk_15 mult ic2574_msk_15.sdf ic2574_sincsinc_trim_submask3_totint20.sdf ic2574_totint20_noisetrim_15 ------------------------------------------------------------------------------- MAKE TOTINT25 NOISE and SN MAPS: cp ../Moment_sincsinc/ic2574_sincsinc_trim_submask3_rembsl_findclump25.sdf . cp ../Moment_sincsinc/ic2574_sincsinc_trim_submask3_totint25.sdf . collapse in= ic2574_sincsinc_trim_submask3_rembsl_findclump25.sdf out=test6 axis=vopt estimator=sum wlim=0.0 mult test2 test6 test7 maths exp="sqrt(IA)" out=test8 IA=test7 cmult test8 0.42 test9 mult sigma_avg test9 ic2574_totint25_noise div ic2574_sincsinc_trim_submask3_totint25.sdf ic2574_totint25_noise ic2574_totint25_SN (contoured SN maps then made using gaia)