  We explore the QCD equation of state at finite baryon density through an
expansion in baryon number fugacity, making use of the recent lattice data on
the four leading Fourier coefficients of the expansion. A state-of-the-art
description of the lattice data at both imaginary $\mu_B$ and $\mu_B = 0$ is
provided by the cluster expansion model (CEM). The CEM pressure function has
three temperature dependent input parameters, which we fix by parameterizing
the available lattice QCD data. This results in a crossover model of the QCD
equation of state at finite baryon density, which can be used in fluid
dynamical simulations of heavy-ion collisions.
