With ongoing advancements in nuclear theory and experimentation, together with a growing body of neutron star (NS) observations, a wealth of information on the equation of state (EOS) for matter at extreme densities has become accessible. Here, we utilize a hybrid EOS formulation that combines an empirical parametrization centered around the nuclear saturation density with a generic three-segment piecewise polytrope model at higher densities. We incorporate data derived from chiral effective field theory (LEFT), perturbative quantum chromodynamics (pQCD), and experiments such as PREX-II and CREX. Furthermore, we examine the influence of a total of 70 NS mass measurements up to April 2023, as well as simultaneous mass and radius measurements derived from the x-ray emission from surface hot spots on NSs. Additionally, we consider constraints on tidal properties inferred from the gravitational waves emitted by coalescing NS binaries. To integrate this extensive and varied array of constraints, we utilize a hierarchical Bayesian statistical framework to simultaneously deduce the EOS and the distribution of NS masses. We find that incorporating data from LEFT significantly tightens the constraints on the EOS of NSs near or below the nuclear saturation density. However, constraints derived from pQCD computations and nuclear experiments such as PREX-II and CREX have minimal impact. Taking into account all available data, we derive estimates for key parameters characterizing the EOS of dense nuclear matter. Specifically, we determine the slope (L) and curvature (Ksym) of the symmetry energy to be 54+10-10 and -158+73 90% credibility, respectively. Additionally, we infer the radius and tidal deformability of an NS with a mass of 1.4 solar masses (M circle dot) to be 12.34+0.43 maximum mass of a nonrotating NS to be 2.22+0.21 -0.19M circle dot with 90% credibility.