The production of Upsilon(1S), Upsilon(2S) and Upsilon(3S) mesons decaying into the dimuon final state is studied with the LHCb detector using a data sample corresponding to an integrated luminosity of 3.3 pb^{-1} collected in proton-proton collisions at a centre-of-mass energy of sqrt{s}=2.76 TeV. The differential production cross-sections times dimuon branching fractions are measured as functions of the Upsilon transverse momentum and rapidity, over the ranges p_T<15 GeV/c and 2.0<y<4.5. The total cross-sections in this kinematic region, assuming unpolarised production, are measured to be sigma(pp -> Upsilon(1S) X) x B(Upsilon(1S) -> mu+mu-) = 1.111 +/- 0.043 +/- 0.044 nb, sigma(pp -> Upsilon(2S) X) x B(Upsilon(2S) -> mu+mu-) = 0.264 +/- 0.023 +/- 0.011 nb, sigma(pp -> Upsilon(3S) X) x B(Upsilon(3S) -> mu+mu-) = 0.159 +/- 0.020 +/- 0.007 nb, where the first uncertainty is statistical and the second systematic.

The production of Upsilon(1S), Upsilon(2S) and Upsilon(3S) mesons in proton-proton collisions at the centre-of-mass energy of sqrt(s)=7 TeV is studied with the LHCb detector. The analysis is based on a data sample of 25 pb-1 collected at the Large Hadron Collider. The Upsilon mesons are reconstructed in the decay mode Upsilon -&gt; mu+ mu- and the signal yields are extracted from a fit to the mu+ mu- invariant mass distributions. The differential production cross-sections times dimuon branching fractions are measured as a function of the Upsilon transverse momentum pT and rapidity y, over the range pT &lt; 15 GeV/c and 2.0 &lt; y &lt; 4.5. The cross-sections times branching fractions, integrated over these kinematic ranges, are measured to be sigma(pp -&gt; Upsilon(1S) X) x B(Upsilon(1S)-&gt;mu+ mu-) = 2.29 {\pm} 0.01 {\pm} 0.10 -0.37 +0.19 nb, sigma(pp -&gt; Upsilon(2S) X) x B(Upsilon(2S)-&gt;mu+ mu-) = 0.562 {\pm} 0.007 {\pm} 0.023 -0.092 +0.048 nb, sigma(pp -&gt; Upsilon(3S) X) x B(Upsilon(3S)-&gt;mu+ mu-) = 0.283 {\pm} 0.005 {\pm} 0.012 -0.048 +0.025 nb, where the first uncertainty is statistical, the second systematic and the third is due to the unknown polarisation of the three Upsilon states.

The production cross-sections of $\Upsilon(1S)$, $\Upsilon(2S)$ and $\Upsilon(3S)$ mesons in proton-proton collisions at $\sqrt{s}$= 13 TeV are measured with a data sample corresponding to an integrated luminosity of $277 \pm 11$ $\rm pb^{-1}$ recorded by the LHCb experiment in 2015. The $\Upsilon$ mesons are reconstructed in the decay mode $\Upsilon\to\mu^{+}\mu^{-}$. The differential production cross-sections times the dimuon branching fractions are measured as a function of the $\Upsilon$ transverse momentum, $p_{\rm T}$, and rapidity, $y$, over the range $0 < p_{\rm T}< 30$ GeV/c and $2.0 < y < 4.5$. The ratios of the cross-sections with respect to the LHCb measurement at $\sqrt{s}$= 8 TeV are also determined. The measurements are compared with theoretical predictions based on NRQCD.

Production of $B_c^+$ mesons in proton-proton collisions at a center-of-mass energy of 8 TeV is studied with data corresponding to an integrated luminosity of $2.0~{\rm fb}^{-1}$ recorded by the LHCb experiment. The ratio of production cross-sections times branching fractions between the $B_c^+\to J/\psi \pi^+$ and $B^+\to J/\psi K^+$ decays is measured as a function of transverse momentum and rapidity in the regions $0 < p_{\rm T} < 20~{\rm GeV}/c$ and $2.0 < y < 4.5$. The ratio in this kinematic range is measured to be $(0.683\pm0.018\pm0.009)\%$, where the first uncertainty is statistical and the second systematic.

The production of $J/\psi$ mesons in proton-proton collisions at $\sqrt{s}$= 7 tev is studied with the LHCb detector at the LHC. The differential cross-section for prompt $J/\psi$ production is measured as a function of the $J/\psi$ transverse momentum p_{\rm T}\ and rapidity $y$ in the fiducial region $p_{\rm T}\in[0;14]$ GeV/c and $y\in[2.0;4.5]$. The differential cross-section and fraction of \j psi\ from $b$-hadron decays are also measured in the same p_{\rm T}\ and $y$ range s. The analysis is based on a data sample corresponding to an integrated luminos ity of 5.2 pb$^{-1}$. The measured cross-sections integrated over the fiducial region are $10.52\pm 0. 04\pm 1.40^{+1.64}_{-2.20}\mu b$ for $\mathrm{prompt}~J/\psi$ production and $1.14 \pm 0. 01\pm 0.16 \mu b$ for $J/\psi$ from $b$-hadron decays, where the first uncertainty is statistical and the second systematic. The $\mathrm{prompt}~J/\psi$ production cross-section is obtained assuming no $J/\psi$ polarisation and the third error indicates the acceptance uncertainty due to this assumption.

The production of $J/\psi$ mesons with rapidity $1.5<y<4.0$ or $-5.0<y<-2.5$ and transverse momentum $p_\mathrm{T}<14 \mathrm{GeV}/c$ is studied with the LHCb detector in proton-lead collisions at a nucleon-nucleon centre-of-mass energy $\sqrt{s_{NN}}=5 \mathrm{TeV}$. The analysis is based on a data sample corresponding to an integrated luminosity of about $1.6 \mathrm{nb}^{-1}$. For the first time the nuclear modification factor and forward-backward production ratio are determined separately for prompt $J/\psi$ mesons and $J/\psi$ from $b$-hadron decays. Clear suppression of prompt $J/\psi$ production with respect to proton-proton collisions at large rapidity is observed, while the production of $J/\psi$ from $b$-hadron decays is less suppressed. These results show good agreement with available theoretical predictions. The measurement shows that cold nuclear matter effects are important for interpretations of the related quark-gluon plasma signatures in heavy-ion collisions.

A measurement of the $Z(\rightarrow\mu^+\mu^-)$+jet production cross-section in $pp$ collisions at a centre-of-mass energy $\sqrt{s} = 7$ TeV is presented. The analysis is based on an integrated luminosity of $1.0\,\text{fb}^{-1}$ recorded by the LHCb experiment. Results are shown with two jet transverse momentum thresholds, 10 and 20 GeV, for both the overall cross-section within the fiducial volume, and for six differential cross-section measurements. The fiducial volume requires that both the jet and the muons from the Z boson decay are produced in the forward direction ($2.0<\eta<4.5$). The results show good agreement with theoretical predictions at the second-order expansion in the coupling of the strong interaction.

The production of $\Upsilon$ mesons in $pp$ collisions at $\sqrt=7$ and $8\,\mathrm{TeV}$ is studied with the LHCb detector using data samples corresponding to an integrated luminosity of $1\,\mathrm{fb}^{-1}$ and $2\,\mathrm{fb}^{-1}$ respectively. The production cross-sections and ratios of cross-sections are measured as functions of the meson transverse momentum $p_T$ and rapidity $y$, for $p_T<30\,\mathrm{GeV}/c$} and $2.0<y<4.5$.

The polarization of the $\Upsilon(1S)$, $\Upsilon(2S)$ and $\Upsilon(3S) $mesons, produced in $pp$ collisions at centre-of-mass energies $\sqrt{s}$=7 and 8TeV, is measured using data samples collected by the LHCb experiment, corresponding to integrated luminosities of 1 and 2fb$^{-1}$, respectively. The measurements are performed in three polarization frames, using $\Upsilon\to\mu^+\mu^-$ decays in the kinematic region of the transverse momentum $p_{T}(\Upsilon)<30GeV/c$, and rapidity $2.2<y(\Upsilon)<4.5$. No large polarization is observed.

The production of $\Upsilon (nS)$ mesons ($n=1,2,3$) in $p$Pb and Pb$p$ collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{NN}}=8.16$ TeV is measured by the LHCb experiment, using a data sample corresponding to an integrated luminosity of 31.8 nb$^{-1}$. The $\Upsilon (nS)$ mesons are reconstructed through their decays into two opposite-sign muons. The measurements comprise the differential production cross-sections of the $\Upsilon (1S)$ and $\Upsilon (2S)$ states, their forward-to-backward ratios and nuclear modification factors, performed as a function of the transverse momentum $p_{\mathrm{T}}$ and rapidity in the nucleon-nucleon centre-of-mass frame $y^*$ of the $\Upsilon (nS)$ states, in the kinematic range $p_{\rm{T}}<25$ GeV/$c$ and $1.5<y^*<4.0$ ($-5.0<y^*<-2.5$) for $p$Pb (Pb$p$) collisions. In addition, production cross-sections for $\Upsilon (3S)$ are measured integrated over phase space and the production ratios between all three $\Upsilon (nS)$ states are determined. The measurements are compared to theoretical predictions and suppressions for quarkonium in $p$Pb collisions are observed.

Exclusive $J/\psi$ and $\psi(2S)$ vector meson production has been observed in the dimuon channel using the LHCb detector. The cross-section times branching fractions to two muons with pseudorapidities between 2.0 and 4.5 are measured to be \begin{equation*} \sigma_{pp\rightarrow J/\psi (\rightarrow \mu^{+} \mu^{-})} (2.0 <\eta_{\mu^{\pm}}< 4.5) = 307 \pm 21 \pm 36~\text{pb}, \end{equation*} \begin{equation*} \sigma_{pp\rightarrow \psi(2S) (\rightarrow \mu^{+} \mu^{-})} (2.0 <\eta_{\mu^{\pm}}< 4.5) = 7.8 \pm 1.3 \pm 1.0~\text{pb}, \end{equation*} where the first uncertainties are statistical and the second are systematic. The measurements are found to be in good agreement with results from previous experiments and theoretical predictions. The $J/\psi$ photoproduction cross-section has been measured as a function of the photon-proton centre-of-mass energy. The results are consistent with measurements obtained at HERA and confirm a similar power law behaviour for the photoproduction cross-section.

The production of J/psi and Upsilon mesons in pp collisions at sqrt(s) = 8 TeV is studied with the LHCb detector. The J/psi and Upsilon mesons are reconstructed in the mu+mu- decay mode and the signal yields are determined with a fit to the mu+mu- invariant mass distributions. The analysis is performed in the rapidity range 2.0<y<4.5 and transverse momentum range 0<p_T<14(15) GeV/c of the J/psi(Upsilon) mesons. The J/psi and Upsilon production cross-sections and the fraction of J/psi mesons from b-hadron decays are measured as a function of the meson p_T and y.

The cross-section for inclusive phi meson production in pp collisions at a centre-of-mass energy of sqrt(s) = 7 TeV has been measured with the LHCb detector at the Large Hadron Collider. The differential cross-section is measured as a function of the phi transverse momentum p_T and rapidity y in the region 0.6 < p_T < 5.0 GeV/c and 2.44 < y < 4.06. The cross-section for inclusive phi production in this kinematic range is sigma(pp -> phi X) = 1758 pm 19(stat) ^{+43}_{-14}(syst) pm 182(scale) microbarn, where the first systematic uncertainty depends on the p_T and y region and the second is related to the overall scale. Predictions based on the Pythia 6.4 generator underestimate the cross-section.

Production cross-sections of prompt charm mesons are measured using data from $pp$ collisions at the LHC at a centre-of-mass energy of $5\,$TeV. The data sample corresponds to an integrated luminosity of $8.60\pm0.33\,$pb$^{-1}$ collected by the LHCb experiment. The production cross-sections of $D^0$, $D^+$, $D_s^+$, and $D^{*+}$ mesons are measured in bins of charm meson transverse momentum, $p_{\text{T}}$, and rapidity, $y$. They cover the rapidity range $2.0<y<4.5$ and transverse momentum ranges $0 < p_{\text{T}} < 10\, \text{GeV}/c$ for $D^0$ and $D^+$ and $1 < p_{\text{T}} < 10\, \text{GeV}/c$ for $D_s^+$ and $D^{*+}$ mesons. The inclusive cross-sections for the four mesons, including charge-conjugate states, within the range of $1 < p_{\text{T}} < 8\, \text{GeV}/c$ are determined to be \sigma(pp\rightarrow D^0 X) = 1004 \pm 3 \pm 54\,\mu\text{b} \sigma(pp\rightarrow D^+ X) = 402 \pm 2 \pm 30\,\mu\text{b} \sigma(pp\rightarrow D_s^+ X) = 170 \pm 4 \pm 16\,\mu\text{b} \sigma(pp\rightarrow D^{*+} X)= 421 \pm 5 \pm 36\,\mu\text{b} where the uncertainties are statistical and systematic, respectively.

The production cross-sections of $J/\psi$ mesons in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=5$ TeV are measured using a data sample corresponding to an integrated luminosity of $9.13\pm0.18~\text{pb}^{-1}$, collected by the LHCb experiment. The cross-sections are measured differentially as a function of transverse momentum, $p_{\text{T}}$, and rapidity, $y$, and separately for $J/\psi$ mesons produced promptly and from beauty hadron decays (nonprompt). With the assumption of unpolarised $J/\psi$ mesons, the production cross-sections integrated over the kinematic range $0<p_{\text{T}}<20~\text{GeV}/c$ and $2.0<y<4.5$ are $8.154\pm0.010\pm0.283~\mu\text{b}$ for prompt $J/\psi$ mesons and $0.820\pm0.003\pm0.034~\mu\text{b}$ for nonprompt $J/\psi$ mesons, where the first uncertainties are statistical and the second systematic. These cross-sections are compared with those at $\sqrt{s}=8$ TeV and $13$ TeV, and are used to update the measurement of the nuclear modification factor in proton-lead collisions for $J/\psi$ mesons at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\text{NN}}}=5$ TeV. The results are compared with theoretical predictions.

The production cross-sections of B mesons are measured in pp collisions at a centre-of-mass energy of 7 TeV, using data collected with the LHCb detector corresponding to a integrated luminosity of 0.36 fb-1. The B+, B0 and Bs0 mesons are reconstructed in the exclusive decays B+ -> J/psi K+, B0 -> J/psi K*0 and Bs0 -> J/psi phi, with J/psi -> mu+ mu-, K*0 -> K+ pi- and phi -> K+ K-. The differential cross-sections are measured as functions of B meson transverse momentum pT and rapidity y, in the range 0 < pT < 40 GeV/c and 2.0 < y < 4.5. The integrated cross-sections in the same pT and y ranges, including charge-conjugate states, are measured to be sigma(pp -> B+ + X) = 38.9 +- 0.3 (stat.) +- 2.5 (syst.) +- 1.3 (norm.) mub, sigma(pp -> B0 + X) = 38.1 +- 0.6 (stat.) +- 3.7 (syst.) +- 4.7 (norm.) mub, sigma(pp -> Bs0 + X) = 10.5 +- 0.2 (stat.) +- 0.8 (syst.) +- 1.0 (norm.) mub, where the third uncertainty arises from the pre-existing branching fraction measurements.

The production of $J/\psi$ mesons in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=13$ TeV is studied with the \lhcb detector. Cross-section measurements are performed as a function of the transverse momentum $p_\mathrm{T}$ and the rapidity $y$ of the $J/\psi$ meson in the region $p_\mathrm{T}<14\mathrm{GeV}/c$ and $2.0<y<4.5$, for both prompt $J/\psi$ mesons and $J/\psi$ mesons from $b$-hadron decays. The production cross-sections integrated over the kinematic coverage are $15.03\pm 0.03\pm 0.94$ $\mu$b for prompt $J/\psi$ and $2.25\pm 0.01\pm 0.14$ $\mu$b for $J/\psi$ from $b$-hadron decays, assuming zero polarization of the $J/\psi$ meson. The first uncertainties are statistical and the second systematic. The cross-section reported for $J/\psi$ mesons from $b$-hadron decays is used to extrapolate to a total $b\bar{b}$ cross-section. The ratios of the cross-sections with respect to $\sqrt{s}=8$ TeV are also determined.

Production cross-sections of prompt charm mesons are measured with the first data from $pp$ collisions at the LHC at a centre-of-mass energy of $13\,\mathrm{TeV}$. The data sample corresponds to an integrated luminosity of $4.98 \pm 0.19\,\mathrm{pb}^{-1}$ collected by the LHCb experiment. The production cross-sections of $D^{0}$, $D^{+}$, $D_{s}^{+}$, and $D^{*+}$ mesons are measured in bins of charm meson transverse momentum, $p_{\mathrm{T}}$, and rapidity, $y$, and cover the range $0 < p_{\mathrm{T}} < 15\,\mathrm{GeV}/c$ and $2.0 < y < 4.5$. The inclusive cross-sections for the four mesons, including charge conjugation, within the range of $1 < p_{\mathrm{T}} < 8\,\mathrm{GeV}/c$ are found to be \begin{equation} \sigma(pp \to D^{0} X) = 2072 \pm 2 \pm 124\,\mu\mathrm{b}\\ \sigma(pp \to D^{+} X) = 834 \pm 2 \pm \phantom{1}78\,\mu\mathrm{b}\\ \sigma(pp \to D_{s}^{+} X) = 353 \pm 9 \pm \phantom{1}76\,\mu\mathrm{b}\\ \sigma(pp \to D^{*+} X) = 784 \pm 4 \pm \phantom{1}87\,\mu\mathrm{b} \end{equation} where the uncertainties are due to statistical and systematic uncertainties, respectively.

The inclusive $D_s^{\pm}$ production asymmetry is measured in $pp$ collisions collected by the LHCb experiment at centre-of-mass energies of $\sqrt{s} =7$ and 8 TeV. Promptly produced $D_s^{\pm}$ mesons are used, which decay as $D_s^{\pm}\to\phi\pi^{\pm}$, with $\phi\to K^+K^-$. The measurement is performed in bins of transverse momentum, $p_{\rm T}$, and rapidity, $y$, covering the range $2.5

A precision measurement of the $Z$ boson production cross-section at $\sqrt{s} = 13$ TeV in the forward region is presented, using $pp$ collision data collected by the LHCb detector, corresponding to an integrated luminosity of 5.1 fb$^{-1}$. The production cross-section is measured using $Z\rightarrow\mu^+\mu^-$ events within the fiducial region defined as pseudorapidity $2.0<\eta<4.5$ and transverse momentum $p_{T}>20$ GeV/$c$ for both muons and dimuon invariant mass $60<M_{\mu\mu}<120$ GeV/$c^2$. The integrated cross-section is determined to be $\sigma (Z \rightarrow \mu^+ \mu^-)$ = 196.4 $\pm$ 0.2 $\pm$ 1.6 $\pm$ 3.9~pb, where the first uncertainty is statistical, the second is systematic, and the third is due to the luminosity determination. The measured results are in agreement with theoretical predictions within uncertainties.

A measurement of the cross-section for pp$ \rightarrow$Z$ \rightarrow$e$^+$e$^-$ is presented using data at $\sqrt{s}=7$ TeV corresponding to an integrated luminosity of 0.94 fb$^{-1}$. The process is measured within the kinematic acceptance $p_{\mathrm{T}}&gt;20$GeV/$c$ and $2&lt;\eta&lt;4.5$ for the daughter electrons and dielectron invariant mass in the range 60-120 GeV/$c^2$. The cross-section is determined to be $$\sigma(pp \rightarrow Z \rightarrow e^+ e^- )=76.0\pm0.8\pm2.0\pm2.6{\rm pb}$$ where the first uncertainty is statistical, the second is systematic and the third is the uncertainty in the luminosity. The measurement is performed as a function of Z rapidity and as a function of an angular variable which is closely related to the Z transverse momentum. The results are compared with previous LHCb measurements and with theoretical predictions from QCD.

A measurement of the production cross-section for $Z$ bosons that decay to muons is presented. The data were recorded by the LHCb detector during $pp$ collisions at a centre-of-mass energy of 7 TeV, and correspond to an integrated luminosity of 1.0 fb$^{-1}$. The cross-section is measured for muons in the pseudorapidity range $2.0 < \eta < 4.5$ with transverse momenta $p_{T} > 20$ GeV/c. The dimuon mass is restricted to $60 < M_{\mu^{+}\mu^{-}} < 120$ GeV/c$^{2}$. The measured cross-section is $$\sigma_{Z\rightarrow\mu^{+}\mu^{-}} = (76.0 \pm 0.3 \pm 0.5 \pm 1.0 \pm 1.3) \, \text{pb}$$ where the uncertainties are due to the sample size, systematic effects, the beam energy and the luminosity. This result is in good agreement with theoretical predictions at next-to-next-to-leading order in perturbative quantum chromodynamics. The cross-section is also measured differentially as a function of kinematic variables of the $Z$ boson. Ratios of the production cross-sections of electroweak bosons are presented using updated LHCb measurements of $W$ boson production. A precise test of the Standard Model is provided by the measurement of the ratio $$\frac{\sigma_{W^{+}\rightarrow\mu^{+}\nu_{\mu}} + \sigma_{W^{-}\rightarrow\mu^{-}\bar{\nu}_{\mu}}}{\sigma_{Z\rightarrow\mu^{+}\mu^{-}}} = 20.63\pm0.09\pm0.12\pm0.05,$$ where the uncertainty due to luminosity cancels.

A measurement of the production cross-section of Z bosons in pp collisions at $\sqrt{s} = 13$ TeV is presented using dimuon and dielectron final states in LHCb data. The cross-section is measured for leptons with pseudorapidities in the range $2.0 < \eta < 4.5$, transverse momenta $p_\text{T} > 20$ GeV and dilepton invariant mass in the range $60<m(\ell\ell)<120$ GeV. The integrated cross-section from averaging the two final states is \begin{equation*}\sigma_{\text{Z}}^{\ell\ell} = 194.3 \pm 0.9 \pm 3.3 \pm 7.6\text{ pb,}\end{equation*} where the first uncertainty is statistical, the second is due to systematic effects, and the third is due to the luminosity determination. In addition, differential cross-sections are measured as functions of the Z boson rapidity, transverse momentum and the angular variable $\phi^*_\eta$.

A study is presented of central exclusive production of $\Upsilon(nS)$ states, where the $\Upsilon(nS)$ resonances decay to the $\mu^+\mu^-$ final state, using $pp$ collision data recorded by the LHCb experiment. The cross-section is measured in the rapidity range $2<y(\Upsilon)<4.5$ where the muons are reconstructed in the pseudorapidity range $2<\eta(\mu^\pm)<4.5$. The data sample corresponds to an integrated luminosity of 2.9 fb$^{-1}$ and was collected at centre-of-mass energies of $7$ TeV and $8$ TeV. The measured $\Upsilon(1S)$ and $\Upsilon(2S)$ production cross-sections are \begin{eqnarray} \sigma(pp \to p\Upsilon(1S)p) &=& 9.0 \pm 2.1 \pm 1.7\textrm{ pb and}\nonumber\\ \sigma(pp \to p\Upsilon(2S)p) &=& 1.3 \pm 0.8 \pm 0.3\textrm{ pb},\nonumber \end{eqnarray} where the first uncertainties are statistical and the second are systematic. The $\Upsilon(1S)$ cross-section is also measured as a function of rapidity and is found to be in good agreement with Standard Model predictions. An upper limit is set at 3.4 pb at the 95\% confidence level for the exclusive $\Upsilon(3S)$ production cross-section, including possible contamination from $\chi_b(3P)\to\Upsilon(3S)\gamma$ decays.

The differential cross-section as a function of rapidity has been measured for the exclusive production of $J/\psi$ and $\psi(2S)$ mesons in proton-proton collisions at $\sqrt{s}=7$ TeV, using data collected by the LHCb experiment, corresponding to an integrated luminosity of 930 pb$^{-1}$. The cross-sections times branching fractions to two muons having pseudorapidities between 2.0 and 4.5 are measured to be $$\begin{array}{rl} \sigma_{pp\rightarrow J/\psi\rightarrow{\mu^+}{\mu^-}}(2.0<\eta_{\mu^\pm }<4.5)=&291\pm 7\pm19 {\rm \ pb},\\ \sigma_{pp\rightarrow\psi(2S)\rightarrow{\mu^+}{\mu^-}}(2.0<\eta_{\mu^\pm}<4.5)=&6.5\pm 0.9\pm 0.4 {\rm \ pb},\end{array}$$ where the first uncertainty is statistical and the second is systematic. The measurements agree with next-to-leading order QCD predictions as well as with models that include saturation effects.