Unverified Commit d5fb6abd authored by Akke Viitanen's avatar Akke Viitanen
Browse files

add Li datapoint, re-normalize baseline

parent 862450fa

src/scripts/plots/plot_number_counts_gband.py

+15 −3
Original line number Diff line number Diff line
@@ -16,7 +16,8 @@ NCOLS = 1 # number of cols for subplot 
# SAVEFIG = "fig/number_counts_lsst_gband.pdf"  # output filename

# SAVEFIG = "fig/number_counts_lsst_gband_20250707.pdf"  # output filename

# SAVEFIG = "fig/number_counts_lsst_gband_20250811.pdf"  # output filename

SAVEFIG = "fig/number_counts_lsst_gband_20250930.pdf"  # output filename

# SAVEFIG = "fig/number_counts_lsst_gband_20250930.pdf"  # output filename

SAVEFIG = "fig/number_counts_lsst_gband_20251017.pdf"  # output filename

###############################################################################
@@ -129,6 +130,15 @@ class SinglePlotter: 
            lsst_fov = 9.6  # in deg2

            ax.plot(x, np.log10(y / lsst_fov), marker="o", fillstyle="none", color="C0", markersize=10)



            # Plot another example with using 18.55 as the FoV from the COSMOS dither pattern

            # (FoV estimated with topcat on 20251017 from baseline v4.0)

            ax.plot(x, np.log10(y / 18.55), marker="^", fillstyle="none", color="C0", markersize=10)



            # Plot Li+2025 (submitted) values, see their Table 1 for the g-band

            x = 26.7

            y = 349.0

            ax.plot(x, np.log10(y), marker="s", fillstyle="none", color="C0", markersize=10)



            # Plot literature values

            if True:

                if b == "g":
@@ -192,7 +202,8 @@ class SinglePlotter: 


            # ax.text(0.10, 0.90, f"lsst-${b}$", transform=ax.transAxes, ha="left", va="top")



            ax.set_xlim(10, 30)

            # ax.set_xlim(10, 30)

            ax.set_xlim(9, 30)

            # ax.set_xticks(np.arange(10, 40))



            ax.set_ylim(-1.5, None)
@@ -205,7 +216,8 @@ class SinglePlotter: 
        # plot_one(self.axes[1, 1], "z")

        # plot_one(self.axes[1, 2], "y")



        self.axes.legend(loc="lower right", fontsize="small")

        # self.axes.legend(loc="lower right", fontsize="small")

        self.axes.legend(loc="upper left", fontsize="small")

        self.axes.set_xlabel(r"$g$ [ABmag]")

        self.axes.set_ylabel(r"$\log N(<g)$ [1/deg$^2$]")

        self.fig.savefig(self.savefig, bbox_inches="tight")