We calculate several diagonal and non-diagonal fluctuations of conserved charges in a system of 2+1+1 quark flavors with physical masses, on a lattice with size $48^3times12$. Higher order fluctuations at $mu_B=0$ are obtained as derivatives of the lower order ones, simulated at imaginary chemical potential. From these correlations and fluctuations we construct ratios of net-baryon number cumulants as functions of temperature and chemical potential, which satisfy the experimental conditions of strangeness neutrality and proton/baryon ratio. Our results qualitatively explain the behavior of the measured cumulant ratios by the STAR collaboration.