Basic Methods for Computing Special Functions

This paper gives an overview of methods for the numerical evaluation of special functions, that is, the functions that arise in many problems from mathematical physics, engineering, probability theory, and other applied sciences. We consider in detail a selection of basic methods which are frequently used in the numerical evaluation of special functions: converging and asymptotic series, including Chebyshev expansions, linear recurrence relations, and numerical quadrature. Several other methods are available and some of these will be discussed in less detail. We give examples of recent software for special functions where these methods are used. We mention a list of new publications on computational aspects of special functions available on our website

Measurement of the azimuthal anisotropy of Y(1S) and Y(2S) mesons in PbPb collisions at √S^{S}NN = 5.02 TeV

The second-order Fourier coefficients (υ$_{2}$) characterizing the azimuthal distributions of Υ(1S) and Υ(2S) mesons produced in PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV are studied. The Υmesons are reconstructed in their dimuon decay channel, as measured by the CMS detector. The collected data set corresponds to an integrated luminosity of 1.7 nb$^{-1}$. The scalar product method is used to extract the υ$_{2}$ coefficients of the azimuthal distributions. Results are reported for the rapidity range |y| < 2.4, in the transverse momentum interval 0 < p$_{T}$ < 50 GeV/c, and in three centrality ranges of 10–30%, 30–50% and 50–90%. In contrast to the J/ψ mesons, the measured υ$_{2}$ values for the Υ mesons are found to be consistent with zero

Measurement of prompt D0^{0} and D‾\overline{D}0^{0} meson azimuthal anisotropy and search for strong electric fields in PbPb collisions at root SNN\sqrt{S_{NN}} = 5.02 TeV

The strong Coulomb field created in ultrarelativistic heavy ion collisions is expected to produce a rapiditydependent difference (Av2) in the second Fourier coefficient of the azimuthal distribution (elliptic flow, v2) between D0 (uc) and D0 (uc) mesons. Motivated by the search for evidence of this field, the CMS detector at the LHC is used to perform the first measurement of Av2. The rapidity-averaged value is found to be (Av2) = 0.001 ? 0.001 (stat)? 0.003 (syst) in PbPb collisions at ?sNN = 5.02 TeV. In addition, the influence of the collision geometry is explored by measuring the D0 and D0mesons v2 and triangular flow coefficient (v3) as functions of rapidity, transverse momentum (pT), and event centrality (a measure of the overlap of the two Pb nuclei). A clear centrality dependence of prompt D0 meson v2 values is observed, while the v3 is largely independent of centrality. These trends are consistent with expectations of flow driven by the initial-state geometry. ? 2021 The Author. Published by Elsevier B.V. This is an open access article under the CC BY licens

Measurement of the Y(1S) pair production cross section and search for resonances decaying to Y(1S)μ⁺μ⁻ in proton-proton collisions at √s = 13 TeV

The fiducial cross section for Y(1S)pair production in proton-proton collisions at a center-of-mass energy of 13TeVin the region where both Y(1S)mesons have an absolute rapidity below 2.0 is measured to be 79 ± 11 (stat) ±6 (syst) ±3 (B)pbassuming the mesons are produced unpolarized. The last uncertainty corresponds to the uncertainty in the Y(1S)meson dimuon branching fraction. The measurement is performed in the final state with four muons using proton-proton collision data collected in 2016 by the CMS experiment at the LHC, corresponding to an integrated luminosity of 35.9fb$^{-1}$. This process serves as a standard model reference in a search for narrow resonances decaying to Y(1S)μ$^{+}$μ$^{-}$ in the same final state. Such a resonance could indicate the existence of a tetraquark that is a bound state of two bquarks and two b̅ antiquarks. The tetraquark search is performed for masses in the vicinity of four times the bottom quark mass, between 17.5 and 19.5GeV, while a generic search for other resonances is performed for masses between 16.5 and 27GeV. No significant excess of events compatible with a narrow resonance is observed in the data. Limits on the production cross section times branching fraction to four muons via an intermediate Y(1S)resonance are set as a function of the resonance mass

Search for heavy resonances decaying to a top quark and a bottom quark in the lepton+jets final state in proton–proton collisions at 13 TeV

A search is presented for narrow heavy resonances decaying to a top quark and a bottom quark using data collected by the CMS experiment at √s = 13TeV in 2016. The data set analyzed corresponds to an integrated luminosity of 35.9 fb−1. Final states that include a single lepton (e, μ), multiple jets, and missing transverse momentum are analyzed. No evidence is found for the production of a W′ boson, and the production of right-handed W′ bosons is excluded at 95% confidence level for masses up to 3.6 TeV depending on the scenario considered. Exclusion limits for W′ bosons are also presented as a function of their coupling strength to left- and right-handed fermions. These limits on a W′ boson decaying via a top and a bottom quark are the most stringent published to date

Measurement of angular parameters from the decay B⁰ → K⁎0^{⁎0} μ⁺ μ⁻ in proton–proton collisions at √s 8 TeV

Angular distributions of the decay B⁰  → K$^{⁎0}$ μ⁺ μ⁻ are studied using a sample of proton–proton collisions at √s=8TeV collected with the CMS detector at the LHC, corresponding to an integrated luminosity of 20.5fb⁻¹ . An angular analysis is performed to determine the P₁ and P$^{\u27}$₅ parameters, where the P$^{\u27}$₅ parameter is of particular interest because of recent measurements that indicate a potential discrepancy with the standard model predictions. Based on a sample of 1397 signal events, the P₁ and P$^{\u27}$₅ parameters are determined as a function of the dimuon invariant mass squared. The measurements are in agreement with predictions based on the standard model