"Figure 3 (right)" of "Measurement of the production cross section for Z + b jets in proton-proton collisions at $\sqrt{s}$ = 13 TeV"

The measurements of the cross section of the \PZ boson, decaying to dielectrons or dimuons, in association with at least one bottom quark jet are performed with proton-proton collision data at $\sqrt{s}=$ 13\TeV. The data sample corresponds to an integrated luminosity of 137.1\fbinv, collected by the CMS experiment at the LHC. The integrated cross sections of \PZ~$+\geq$~1~\cPqb~jet\xspace and \PZ~$+\geq$~2~\cPqb~jets\xspace are reported for the electron, muon, and combined channels with full Run 2 data. The measured integrated cross sections are 6.52 $\pm$ 0.04 (stat) $\pm$ 0.40 (syst) $\pm$ 0.14 (theo) pb for \PZ~$+\geq$~1~\cPqb~jet\xspace and 0.65 $\pm$ 0.03 (stat) $\pm$ 0.07 (syst) $\pm$ 0.02 (theo) pb for \PZ~$+\geq$~2~\cPqb~jets\xspace. The differential cross section distributions are measured as a function of various kinematic observables that are useful for precision tests of the perturbative quantum chromodynamics predictions. The ratios of integrated and differential cross sections of \PZ~$+\geq$~2~\cPqb~jets\xspace and \PZ~$+\geq$~1~\cPqb~jet\xspace processes are also determined. The value of the integrated cross section ratio measured in the combined channel is 0.100 $\pm$ 0.005 (stat) $\pm$ 0.007 (syst) $\pm$ 0.003 (theo). All of the measurements are compared with predictions from Monte Carlo simulations.