def calc_coldens_products(nhi, av, dgr, nhi_error=0.0, av_error=0.0, dgr_error=0.0, ): from myimage_analysis import calculate_nhi, calculate_noise_cube, \ calculate_sd, calculate_nh2, calculate_nh2_error # Calculate N(H2) and error # --------------------------------------------------------------------------- nh2 = calculate_nh2(nhi_image=nhi, av_image=av, dgr=dgr) # nh2 = (av / dgr - nhi) / 2 # nh2_error = ((nh2 * ((av_error / av)**2 + (dgr_error / dgr)**2)**0.5)**2 \ # - nhi_error**2.0)**0.5 / 2 comp_av_error = av_error / dgr comp_dgr_error = - av / dgr**2 * dgr_error comp_nhi_error = nhi_error nh2_error = 0.5 * (comp_av_error**2 + \ comp_nhi_error**2)**0.5 nh = nhi + nh2 nh_error = (nh2_error**2 + nhi_error**2)**0.5 # Convert to column density to surface density # --------------------------------------------------------------------------- hi_sd = calculate_sd(nhi, sd_factor=1/1.25) hi_sd_error = calculate_sd(nhi_error, sd_factor=1/1.25) h2_sd = calculate_sd(nh2, sd_factor=1/0.625) h2_sd_error = calculate_sd(nh2_error, sd_factor=1/0.625) h_sd = hi_sd + h2_sd h_sd_error = (hi_sd_error**2 + h2_sd_error**2)**0.5 # Write ratio between H2 and HI # --------------------------------------------------------------------------- rh2 = h2_sd / hi_sd comp_hi_error = - h2_sd / hi_sd**2 * hi_sd_error comp_h2_error = h2_sd_error / hi_sd rh2_error = (comp_hi_error**2 + comp_h2_error**2)**0.5 return ((nhi, nh2, nh, hi_sd, h2_sd, h_sd, rh2), (nhi_error, nh2_error, nh_error, hi_sd_error, h2_sd_error, h_sd_error, rh2_error))
def calc_coldens_products( nhi, av, dgr, nhi_error=0.0, av_error=0.0, dgr_error=0.0, ): from myimage_analysis import calculate_nhi, calculate_noise_cube, \ calculate_sd, calculate_nh2, calculate_nh2_error # Calculate N(H2) and error # --------------------------------------------------------------------------- nh2 = calculate_nh2(nhi_image=nhi, av_image=av, dgr=dgr) # nh2 = (av / dgr - nhi) / 2 # nh2_error = ((nh2 * ((av_error / av)**2 + (dgr_error / dgr)**2)**0.5)**2 \ # - nhi_error**2.0)**0.5 / 2 comp_av_error = av_error / dgr comp_dgr_error = -av / dgr**2 * dgr_error comp_nhi_error = nhi_error nh2_error = 0.5 * (comp_av_error**2 + \ comp_nhi_error**2)**0.5 nh = nhi + nh2 nh_error = (nh2_error**2 + nhi_error**2)**0.5 # Convert to column density to surface density # --------------------------------------------------------------------------- hi_sd = calculate_sd(nhi, sd_factor=1 / 1.25) hi_sd_error = calculate_sd(nhi_error, sd_factor=1 / 1.25) h2_sd = calculate_sd(nh2, sd_factor=1 / 0.625) h2_sd_error = calculate_sd(nh2_error, sd_factor=1 / 0.625) h_sd = hi_sd + h2_sd h_sd_error = (hi_sd_error**2 + h2_sd_error**2)**0.5 # Write ratio between H2 and HI # --------------------------------------------------------------------------- rh2 = h2_sd / hi_sd comp_hi_error = -h2_sd / hi_sd**2 * hi_sd_error comp_h2_error = h2_sd_error / hi_sd rh2_error = (comp_hi_error**2 + comp_h2_error**2)**0.5 return ((nhi, nh2, nh, hi_sd, h2_sd, h_sd, rh2), (nhi_error, nh2_error, nh_error, hi_sd_error, h2_sd_error, h_sd_error, rh2_error))