  We present a novel construction of the QCD equation of state (EoS) at finite
baryon density. Our work combines a recently proposed resummation scheme for
lattice QCD results with the universal critical behavior at the QCD critical
point. This allows us to obtain a family of equations of state in the range $0
\leq \mu_B \leq 700$ MeV and 25 MeV $\leq T \leq 800$ MeV, which match lattice
QCD results near $\mu_B=0$ while featuring a critical point in the 3D Ising
model universality class. The position of the critical point can be chosen
within the range accessible to beam-energy scan heavy-ion collision
experiments. The strength of the singularity and the shape of the critical
region are parameterized using a standard parameter set. We impose stability
and causality constraints and discuss the available ranges of critical point
parameter choices, finding that they extend beyond earlier parametric QCD EoS
proposals. We present thermodynamic observables, including baryon density,
pressure, entropy density, energy density, baryon susceptibility and speed of
sound, that cover a wide range in the QCD phase diagram relevant for
experimental exploration.
