def make_ratio_integ(): """ This is an alternate implementation to that already included in make_ratiotem_cubesims """ npixflat = bmasksm_rs.include().sum(axis=0) sum321 = cube321.with_mask(bmasksm_rs).sum(axis=0) sum303 = cube303.with_mask(bmasksm_rs).sum(axis=0) integ303321b = sum321 / sum303 error303 = cube303.with_mask(~bmasksm_rs).std(axis=0) error321 = cube321.with_mask(~bmasksm_rs).std(axis=0) eintegratio303321 = ((integ303321b.value**2 * (error303.value**2 / (sum303.value / npixflat)**2 + error321.value**2 / (sum321.value / npixflat)**2))**0.5) ok = (integ303321b.value / eintegratio303321) > 3 integ303321ok = integ303321b.value integ303321ok[~ok] = np.nan hdu = cube321[0, :, :].hdu hdu.data = integ303321ok hdu.writeto(paths.hpath('H2CO_321220_to_303202_integ_smoothmask.fits'), clobber=True) hdu.data = eintegratio303321 hdu.writeto(paths.hpath('H2CO_321220_to_303202_einteg_smoothmask.fits'), clobber=True)
def make_ratio_max(): npixflat = bmasksm_rs.include().sum(axis=0) max321 = cube321.with_mask(bmasksm_rs).max(axis=0) max303 = cube303.with_mask(bmasksm_rs).max(axis=0) maxratio303321b = max321/max303 error303 = cube303.with_mask(~bmasksm_rs).std(axis=0) error321 = cube321.with_mask(~bmasksm_rs).std(axis=0) eintegratio303321 = ((maxratio303321b.value**2 * (error303.value**2/(max303.value)**2 + error321.value**2/(max321.value)**2))**0.5) ok = (maxratio303321b.value/eintegratio303321) > 3 maxratio303321ok = maxratio303321b.value maxratio303321ok[~ok] = np.nan hdu = cube321[0,:,:].hdu hdu.data = maxratio303321ok hdu.writeto(paths.hpath('H2CO_321220_to_303202_maxratio_smoothmask.fits'),clobber=True) hdu.data = eintegratio303321 hdu.writeto(paths.hpath('H2CO_321220_to_303202_emaxratio_smoothmask.fits'),clobber=True)
def make_ratio_integ(): """ This is an alternate implementation to that already included in make_ratiotem_cubesims """ npixflat = bmasksm_rs.include().sum(axis=0) sum321 = cube321.with_mask(bmasksm_rs).sum(axis=0) sum303 = cube303.with_mask(bmasksm_rs).sum(axis=0) integ303321b = sum321/sum303 error303 = cube303.with_mask(~bmasksm_rs).std(axis=0) error321 = cube321.with_mask(~bmasksm_rs).std(axis=0) eintegratio303321 = ((integ303321b.value**2 * (error303.value**2/(sum303.value/npixflat)**2 + error321.value**2/(sum321.value/npixflat)**2))**0.5) ok = (integ303321b.value/eintegratio303321) > 3 integ303321ok = integ303321b.value integ303321ok[~ok] = np.nan hdu = cube321[0,:,:].hdu hdu.data = integ303321ok hdu.writeto(paths.hpath('H2CO_321220_to_303202_integ_smoothmask.fits'),clobber=True) hdu.data = eintegratio303321 hdu.writeto(paths.hpath('H2CO_321220_to_303202_einteg_smoothmask.fits'),clobber=True)
def make_ratio_max(): npixflat = bmasksm_rs.include().sum(axis=0) max321 = cube321.with_mask(bmasksm_rs).max(axis=0) max303 = cube303.with_mask(bmasksm_rs).max(axis=0) maxratio303321b = max321 / max303 error303 = cube303.with_mask(~bmasksm_rs).std(axis=0) error321 = cube321.with_mask(~bmasksm_rs).std(axis=0) eintegratio303321 = ((maxratio303321b.value**2 * (error303.value**2 / (max303.value)**2 + error321.value**2 / (max321.value)**2))**0.5) ok = (maxratio303321b.value / eintegratio303321) > 3 maxratio303321ok = maxratio303321b.value maxratio303321ok[~ok] = np.nan hdu = cube321[0, :, :].hdu hdu.data = maxratio303321ok hdu.writeto(paths.hpath('H2CO_321220_to_303202_maxratio_smoothmask.fits'), clobber=True) hdu.data = eintegratio303321 hdu.writeto(paths.hpath('H2CO_321220_to_303202_emaxratio_smoothmask.fits'), clobber=True)