Ejemplo n.º 1
0
    def test_plot_lc(self):
        n = 100
        start, end = 0.1, 200.
        time = np.exp(np.linspace(np.log(start), np.log(end), n))

        popov = Popov('test', R=400., M=15., Mni=0.04, E=1.)
        popov.plot_Lbol(time)
        plt.show()
Ejemplo n.º 2
0
    def test_plot_lc(self):
        n = 100
        start, end = 0.1, 200.
        time = np.exp(np.linspace(np.log(start), np.log(end), n))

        popov = Popov('test', R=400., M=15., Mni=0.04, E=1.)
        popov.plot_Lbol(time)
        plt.show()
Ejemplo n.º 3
0
    def test_popov_SN1999em_emcee(self):
        n = 100
        start, end = 0.1, 200.
        jd_shift = 20.
        dm = -29.38  # D = 7.5e6 pc
        # dm = -30.4  # D = 12.e6 pc
        time = np.exp(np.linspace(np.log(start), np.log(end), n))
        popov = Popov('test', R=450., M=15., Mni=0.04, E=0.7)
        lc_m = popov.LCBol(time)

        # ax = popov.plot_Lbol(time)
        # sn1999em.plot_ubv(ax, path=sn1999em.sn_path, jd_shift=jd_shift, mshift=dm)
        # plt.show()
        curves_o = sn1999em.read_curves()
        lc_o = curves_o.get('V')
        lc_o.mshift = dm
        print("Run: find tshift with bayesian: obs band %s with %s ..." %
              (lc_o.Band.Name, popov))
        fitter = FitMCMC()
        fitter.is_debug = True
        res = fitter.fit_lc(lc_o, lc_m)
        tshift, tsigma = res.tshift, res.tsigma
        print("Result: tshift= %s tsigma= %s ..." % (tshift, tsigma))

        ax = popov.plot_Lbol(time)
        lcp.lc_plot(lc_m, ax)
        lcp.lc_plot(lc_o, ax, is_line=False)
        # sn1999em.plot_ubv(ax, path=sn1999em.sn_path, jd_shift=-tshift, mshift=dm)
        plt.show()
Ejemplo n.º 4
0
    def test_fit_time_popov_SN1999em(self):
        jd_shift = 20.
        dm = -29.38  # D = 7.5e6 pc
        # dm = -30.4  # D = 12.e6 pc
        curves = sn1999em.read_curves()
        lc = curves.get('V')
        lc.mshift = dm

        time = lc.Time - lc.Tmin
        # time = np.exp(np.linspace(np.log(start), np.log(end), n))
        popov = Popov('test', R=450., M=15., Mni=0.04, E=0.7)
        mags = popov.MagBol(time)

        # fit
        tshift = myfit(mags, lc)

        # plot
        ax = popov.plot_Lbol(time)
        x = lc.Time + tshift
        # x = lc.Time + jd_shift + res
        y = lc.Mag
        ax.plot(x,
                y,
                label='%s SN 1999em' % lc.Band.Name,
                ls=":",
                color='red',
                markersize=8,
                marker="o")
        plt.show()
Ejemplo n.º 5
0
    def test_popov_SN87A(self):
        n = 100
        start, end = 0.1, 1200.
        jd_shift = -2446850.  # moment of explosion SN 1987A, Hamuy 1988, doi:10.1086/114613
        dm = -18.6  # D = 7.5e6 pc
        time = np.exp(np.linspace(np.log(start), np.log(end), n))

        popov = Popov('test', R=50., M=15., Mni=0.07, E=1.)
        ax = popov.plot_Lbol(time)
        sn87a.plot_ubv(ax, path=sn87a.sn_path, jd_shift=jd_shift, mshift=dm)
        plt.show()
Ejemplo n.º 6
0
 def test_popov_rednova(self):
     n = 100
     start, end = 0.1, 90.
     jd_shift = 10.
     dm = 24.38  # D = 7.5e5 pc
     # dm = -30.4  # D = 12.e6 pc
     time = np.exp(np.linspace(np.log(start), np.log(end), n))
     popov = Popov('test', R=35., M=6., Mni=0., E=0.008)
     ax = popov.plot_Lbol(time)
     rednova.plot_ubv(ax, path=rednova.sn_path, jd_shift=jd_shift, mshift=-dm)
     plt.show()
Ejemplo n.º 7
0
 def test_popov_SN1999em(self):
     n = 100
     start, end = 0.1, 200.
     jd_shift = 20.
     dm = -29.38  # D = 7.5e6 pc
     # dm = -30.4  # D = 12.e6 pc
     time = np.exp(np.linspace(np.log(start), np.log(end), n))
     popov = Popov('test', R=450., M=15., Mni=0.04, E=0.7)
     ax = popov.plot_Lbol(time)
     sn1999em.plot_ubv(ax, path=sn1999em.sn_path, jd_shift=jd_shift, mshift=dm)
     plt.show()
Ejemplo n.º 8
0
    def test_popov_SN87A(self):
        n = 100
        start, end = 0.1, 1200.
        jd_shift = -2446850.  # moment of explosion SN 1987A, Hamuy 1988, doi:10.1086/114613
        dm = -18.6  # D = 7.5e6 pc
        time = np.exp(np.linspace(np.log(start), np.log(end), n))

        popov = Popov('test', R=50., M=15., Mni=0.07, E=1.)
        ax = popov.plot_Lbol(time)
        sn87a.plot_ubv(ax, path=sn87a.sn_path, jd_shift=jd_shift, mshift=dm)
        plt.show()
Ejemplo n.º 9
0
 def test_popov_rednova(self):
     n = 100
     start, end = 0.1, 90.
     jd_shift = -2500000 + 42963.
     dm = -24.38  # D = 7.5e5 pc
     # dm = -30.4  # D = 12.e6 pc
     time = np.exp(np.linspace(np.log(start), np.log(end), n))
     popov = Popov('test', R=35., M=1., Mni=0., E=0.002)
     ax = popov.plot_Lbol(time)
     rednova.plot_ubv(ax,
                      path=rednova.sn_path,
                      jd_shift=jd_shift,
                      mshift=dm)
     plt.show()
Ejemplo n.º 10
0
 def test_popov_SN1999em(self):
     n = 100
     start, end = 0.1, 200.
     jd_shift = 20.
     dm = -29.38  # D = 7.5e6 pc
     # dm = -30.4  # D = 12.e6 pc
     time = np.exp(np.linspace(np.log(start), np.log(end), n))
     popov = Popov('test', R=450., M=15., Mni=0.04, E=0.7)
     ax = popov.plot_Lbol(time)
     sn1999em.plot_ubv(ax,
                       path=sn1999em.sn_path,
                       jd_shift=jd_shift,
                       mshift=dm)
     plt.show()
Ejemplo n.º 11
0
    def test_fit_time_popov_SN1999em(self):
        jd_shift = 20.
        dm = -29.38  # D = 7.5e6 pc
        # dm = -30.4  # D = 12.e6 pc
        curves = sn1999em.read_curves()
        lc = curves.get('V')
        lc.mshift = dm

        time = lc.Time - lc.tmin
        # time = np.exp(np.linspace(np.log(start), np.log(end), n))
        popov = Popov('test', R=450., M=15., Mni=0.04, E=0.7)
        mags = popov.MagBol(time)

        # fit
        tshift = myfit(mags, lc)

        # plot
        ax = popov.plot_Lbol(time)
        x = lc.Time + tshift
        # x = lc.Time + jd_shift + res
        y = lc.Mag
        ax.plot(x, y, label='%s SN 1999em' % lc.Band.Name,
                ls=".", color='red', markersize=8, marker="o")
        plt.show()