Pseudorapidity (eta) distributions of charged particles produced in proton-proton collisions at a centre-of-mass energy of 8 TeV are measured in the ranges abs(eta) < 2.2 and 5.3 < abs(eta) < 6.4 covered by the CMS and TOTEM detectors, respectively. The data correspond to an integrated luminosity of 45 inverse microbarns. Measurements are presented for three event categories. The most inclusive category is sensitive to 91-96% of the total inelastic proton-proton cross section. The other two categories are disjoint subsets of the inclusive sample that are either enhanced or depleted in single diffractive dissociation events. The data are compared to models used to describe high-energy hadronic interactions. None of the models considered provide a consistent description of the measured distributions.

Measurements are presented of the single-diffractive dijet cross section and the diffractive cross section as a function of the proton fractional momentum loss $\xi$ and the four-momentum transfer squared $t$. Both processes pp$\to$pX and pp$\to$Xp, ie with the proton scattering to either side of the interaction point, are measured, where X includes at least two jets; the results of the two processes are averaged. The analyses are based on data collected simultaneously with the CMS and TOTEM detectors at the LHC in proton-proton collisions at $\sqrt{s}$= 8 TeV during a dedicated run with $\beta^{\ast} =$ 90 m at low instantaneous luminosity and correspond to an integrated luminosity of 37.5 nb$^{-1}$. The single-diffractive dijet cross section $\sigma^\mathrm{pX}_{\mathrm{jj}}$, in the kinematic region $\xi \lt$ 0.1, 0.03 $\lt |$t$| \lt 1$ GeV$^2$, with at least two jets with transverse momentum $p_\mathrm{T} >$ 40 GeV, and pseudorapidity $|\eta| \lt$ 4.4, is 21.7$\pm$0.9 (stat)$^{+3.0}_{-3.3}$ (syst) $\pm$ 0.9 (lum) nb. The ratio of the single-diffractive to inclusive dijet yields, normalised per unit of $\xi$, is presented as a function of $x$, the longitudinal momentum fraction of the proton carried by the struck parton. The ratio in the kinematic region defined above, for $x$ values in the range $-$2.9 $\leq \log_{10} x \leq$$-$1.6, is $R = (\sigma^\mathrm{pX}_{\mathrm{jj}}/\Delta\xi)/\sigma_{\mathrm{jj}} =$ 0.025$\pm$0.001 (stat) $\pm$ 0.003 (syst), where $\sigma^\mathrm{pX}_{\mathrm{jj}}$ and $\sigma_{\mathrm{jj}}$ are the single-diffractive and inclusive dijet cross sections, respectively. The results are compared with predictions from models of diffractive and nondiffractive interactions. Monte Carlo predictions based on the HERA diffractive parton distribution functions agree well with the data when corrected for the effect of soft rescattering between the spectator partons.

A search for exclusive two-photon production via photon exchange in proton-proton collisions, pp $\to$ p$\gamma\gamma$p with intact protons, is presented. The data correspond to an integrated luminosity of 9.4 fb$^{-1}$ collected in 2016 using the CMS and TOTEM detectors at a center-of-mass energy of 13 TeV at the LHC. Events are selected with a diphoton invariant mass above 350 GeV and with both protons intact in the final state, to reduce backgrounds from strong interactions. The events of interest are those where the invariant mass and rapidity calculated from the momentum losses of the forward-moving protons matches the mass and rapidity of the central, two-photon system. No events are found that satisfy this condition. Interpreting this result in an effective dimension-8 extension of the standard model, the first limits are set on the two anomalous four-photon coupling parameters. If the other parameter is constrained to its standard model value, the limits at 95% CL are $\lvert\zeta_1\rvert$ $\lt$ 2.9 $\times$ 10$^{-13}$ GeV$^{-4}$ and $\lvert\zeta_2\rvert$ $\lt$ 6.0 $\times$ 10$^{-13}$ GeV$^{-4}$.

A search is performed for exclusive high-mass $\gamma\gamma$$\to$ WW and $\gamma\gamma$$\to$ ZZ production in proton-proton collisions using intact forward protons reconstructed in near-beam detectors, with both weak bosons decaying into boosted and merged jets. The analysis is based on a sample of proton-proton collisions collected by the CMS and TOTEM experiments at $\sqrt{s}$ = 13 TeV, corresponding to an integrated luminosity of 100 fb$^{-1}$. No excess above the standard model background prediction is observed, and upper limits are set on the pp $\to$ pWWp and pp $\to$ pZZp cross sections in a fiducial region defined by the diboson invariant mass $m$(VV) $\lt$ 1 TeV (with V = W, Z) and proton fractional momentum loss 0.04 $\lt$$\xi$$\lt$ 0.20. The results are interpreted as new limits on dimension-6 and dimension-8 anomalous quartic gauge couplings.

A search is presented for high-mass exclusive diphoton production via photon-photon fusion in proton-proton collisions at $\sqrt{s}$ = 13 TeV in events where both protons survive the interaction. The analysis utilizes data corresponding to an integrated luminosity of 103 fb$^{-1}$ collected in 2016-2018 with the central CMS detector and the CMS and TOTEM precision proton spectrometer (PPS). Events that have two photons with high transverse momenta ($p_\mathrm{T}^\gamma > $ 100 GeV), back-to-back in azimuth, and with a large diphoton invariant mass ($m_{\gamma\gamma} \gt$ 350 GeV) are selected. To remove the dominant inclusive diphoton backgrounds, the kinematic properties of the protons detected in PPS are required to match those of the central diphoton system. Only events having opposite-side forward protons detected with a fractional momentum loss between 0.035 and 0.15 (0.18) for the detectors on the negative (positive) side of CMS are considered. One exclusive diphoton candidate is observed for an expected background of 1.1 events. Limits at 95% confidence level are derived for the four-photon anomalous coupling parameters $\lvert\zeta_1\rvert \lt$ 0.073 TeV$^{-4}$ and $\lvert\zeta_2\rvert \lt$ 0.15 TeV$^{-4}$, using an effective field theory. Additionally, upper limits are placed on the production of axion-like particles with coupling strength to photons $f^{-1}$ that varies from 0.03 TeV$^{-1}$ to 1 TeV$^{-1}$ over the mass range from 500 to 2000 GeV.

A search for the central exclusive production of top quark-antiquark pairs ($\mathrm{t\bar{t}}$) is performed for the first time using proton-tagged events in proton-proton collisions at the LHC at a centre-of-mass energy of 13 TeV. The data correspond to an integrated luminosity of 29.4 fb$^{-1}$. The $\mathrm{t\bar{t}}$ decay products are reconstructed using the central CMS detector, while forward protons are measured in the CMS-TOTEM precision proton spectrometer. An observed (expected) upper bound on the production cross section of 0.59 (1.14) pb is set at 95% confidence level, for collisions of protons with fractional momentum losses between 2 and 20%.

The central exclusive production of charged-hadron pairs in pp collisions at a centre-of-mass energy of 13\TeV is examined, based on data collected in a special high-$\beta^*$ run of the LHC. The nonresonant continuum processes are studied with the invariant mass of the centrally produced two-pion system in the resonance-free region, $m_{\pi^+\pi^-}$$\lt$ 0.7 GeV or $m_{\pi^+\pi^-}$$\gt$ 1.8 GeV. Differential cross sections as functions of the azimuthal angle between the surviving protons, squared exchanged four-momenta, and $m_{\pi^+\pi^-}$ are measured in a wide region of scattered proton transverse momenta, between 0.2 and 0.8 GeV, and for pion rapidities $\lvert y\rvert$$\lt$ 2. A rich structure of interactions related to double-pomeron exchange is observed. A parabolic minimum in the distribution of the two-proton azimuthal angle is observed for the first time. It can be interpreted as an effect of additional pomeron exchanges between the protons from the interference between the bare and the rescattered amplitudes. After model tuning, various physical quantities are determined that are related to the pomeron cross section, proton-pomeron and meson-pomeron form factors, pomeron trajectory and intercept, and coefficients of diffractive eigenstates of the proton.

A generic search is presented for the associated production of a Z boson or a photon with an additional unspecified massive particle X, pp $\to$ pp + Z/$\gamma$ + X, in proton-tagged events from proton-proton collisions at $\sqrt{s}$ = 13 TeV, recorded in 2017 with the CMS detector and the CMS-TOTEM precision proton spectrometer. The missing mass spectrum is analysed in the 600-1600 GeV range and a fit is performed to search for possible deviations from the background expectation. No significant excess in data with respect to the background predictions has been observed. Model-independent upper limits on the visible production cross section of pp $\to$ pp + Z/$\gamma$ + X are set.

The SciBooNE Collaboration has performed a search for charged current coherent pion production from muon neutrinos scattering on carbon, $\nu_{\mu}$ $^{12}C \to \mu^{-12}C \pi^+$, with two distinct data samples. No evidence for coherent pion production is observed. We set 90% confidence level upper limits on the cross section ratio of charged current coherent pion production to the total charged current cross section at $ 0.67 \times 10^{-2}$ at mean neutrino energy 1.1 GeV and 1.36\times 10^{-2} at mean neutrino energy 2.2 GeV.

The SciBooNE Collaboration reports inclusive neutral current neutral pion production by a muon neutrino beam on a polystyrene target (C8H8). We obtain (7.7 \pm 0.5(stat.) \pm 0.5 (sys.)) x 10^(-2) as the ratio of the neutral current neutral pion production to total charged current cross section; the mean energy of neutrinos producing detected neutral pions is 1.1 GeV. The result agrees with the Rein-Sehgal model implemented in our neutrino interaction simulation program with nuclear effects. The spectrum shape of the neutral pion momentum and angle agree with the model. We also measure the ratio of the neutral current coherent pion production to total charged current cross section to be (0.7 \pm 0.4) x 10^(-2).

We measure triangular flow relative to the reaction plane at 3 GeV center-of-mass energy in Au+Au collisions at the BNL Relativistic Heavy Ion Collider. A significant $v_3$ signal for protons is observed, which increases for higher rapidity, higher transverse momentum, and more peripheral collisions. The triangular flow is essentially rapidity-odd with a slope at mid-rapidity, $dv_3/dy|_{(y=0)}$, opposite in sign compared to the slope for directed flow. No significant $v_3$ signal is observed for charged pions and kaons. Comparisons with models suggest that a mean field potential is required to describe these results, and that the triangular shape of the participant nucleons is the result of stopping and nuclear geometry.

We present data from Fermilab experiment E781 (SELEX) on the hadroproduction asymmetry for anti-Lambda_c compared to Lambda_c+ as a function of xF and pt2 distributions for Lambda_c+. These data were measured in the same apparatus using incident pi-, sigma- beams at 600 GeV/c and proton beam at 540 GeV/c. The asymmetry is studied as a function of xF. In the forward hemisphere with xF >= 0.2 both baryon beams exhibit very strong preference for producing charm baryons rather than charm antibaryons, while the pion beam asymmetry is much smaller. In this energy regime the results show that beam fragments play a major role in the kinematics of Lambda_c formation, as suggested by the leading quark picture.

Total cross sections for Sigma- and pi- on beryllium, carbon, polyethylene and copper as well as total cross sections for protons on beryllium and carbon have been measured in a broad momentum range around 600GeV/c. These measurements were performed with a transmission technique adapted to the SELEX hyperon-beam experiment at Fermilab. We report on results obtained for hadron-nucleus cross sections and on results for sigma_tot(Sigma- N) and sigma_tot(pi- N), which were deduced from nuclear cross sections.

The Sigma^- mean squared charge radius has been measured in the space-like Q^2 range 0.035-0.105 GeV^2/c^2 by elastic scattering of a Sigma^- beam off atomic electrons. The measurement was performed with the SELEX (E781) spectrometer using the Fermilab hyperon beam at a mean energy of 610 GeV/c. We obtain <r^2> = (0.61 +/- 0.12 (stat.) +/- 0.09 (syst.)) fm^2. The proton and pi^- charge radii were measured as well and are consistent with results of other experiments. Our result agrees with the recently measured strong interaction radius of the Sigma^-.

The production of $D_s^-$ relative to $D_s^+$ as a function of $x_F $ with 600 GeV/c $\Sigma^-$ beam is measured in the interval $0.15 < x_F < 0.7$ by the SELEX (E781) experiment at Fermilab. The integrated charge asymmetries with 600 GeV/c $\Sigma^-$ beam ($0.53\pm0.06$) and $\pi^-$ beam ($0.06\pm0.11$) are also compared. The results show the $\Sigma^-$ beam fragments play a role in the production of $D_s^-$, as suggested by the leading quark model.

We report on the first observation of open charm production in neutral current deep inelastic neutrino scattering as seen in the NuTeV detector at Fermilab. The production rate is shown to be consistent with a pure gluon-$% Z^{0}$ boson production model, and the observed level of charm production is used to determine the effective charm mass. As part of our analysis, we also obtain a new measurement for the proton-nucleon charm production cross section at $\sqrt{s}=38.8$ GeV.

We present measurements of the semi-inclusive cross sections for νμ- and ν¯μ-nucleon deep inelastic scattering interactions with two oppositely charged muons in the final state. These events dominantly arise from the production of a charm quark during the scattering process. The measurement was obtained from the analysis of 5102 νμ-induced and 1458 ν¯μ-induced events collected with the NuTeV detector exposed to a sign-selected beam at the Fermilab Tevatron. We also extract a cross-section measurement from a reanalysis of 5030 νμ-induced and 1060 ν¯μ-induced events collected from the exposure of the same detector to a quad-triplet beam by the Chicago Columbia Fermilab Rochester (CCFR) experiment. The results are combined to obtain the most statistically precise measurement of neutrino-induced dimuon production cross sections to date. These measurements should be of broad use to phenomenologists interested in the dynamics of charm production, the strangeness content of the nucleon, and the Cabibbo-Kobayashi-Maskawa matrix element Vcd.

We present evidence for the diffractive processes nu_mu Fe -> mu^- D_s^+ (D_s^*+) Fe and nubar_mu Fe -> mu^+ D_s^- (D_s^*-) Fe using the Fermilab SSQT neutrino beam and the Lab E neutrino detector. We observe the neutrino trident reactions nu_mu Fe -> nu_mu mu^- mu^+ Fe and nubar_mu Fe -> nubar_mu mu^+ mu^- Fe at rates consistent with Standard Model expectations. We see no evidence for neutral-current production of J/psi via either diffractive or deep inelastic scattering mechanisms.

The NuTeV experiment at Fermilab has used a sign-selected neutrino beam to perform a search for the lepton number violating process $\bar{\nu}_mu e^- \to \mu^- \bar{\nu}_e$, and to measure the cross-section of the Standard Model inverse muon decay process $\nu_{\mu} e^- \to \mu^- \nu_e$. NuTeV measures the inverse muon decay asymptotic cross-section $\sigma/E$ to be 13.8 $\pm$ 1.2 $\pm$ 1.4 x $10^{-42} cm^2$/GeV. The experiment also observes no evidence for lepton number violation and places one of the most restrictive limits on the LNV/IMD cross-section ratio at $\sigma (\bar{\nu}_{\mu} e^- \to \mu^- \bar{\nu}_e) /\sigma (\nu_{\mu}e^- \to \mu^- \nu_e$) $\le$ 1.7% at 90% C.L. for V-A couplings and $\le$ 0.6% for scalar couplings.

We present a new measurement of the difference between the nucleon strange and antistrange quark distributions from dimuon events recorded by the NuTeV experiment at Fermilab. This analysis is the first to use a complete next to leading order QCD d escription of charm production from neutrino scattering. Dimuon events in neutrino deep inelastic scattering allow direct and independent study of the strange and antistrange content of the nucleon. We find a positive strange asymmetry with a significance of 1.6sigma . We also report a new measurement of the charm mass.

The NuTeV experiment at Fermilab has obtained a unique high statistics sample of neutrino and anti-neutrino interactions using its high-energy sign-selected beam. We present a measurement of the differential cross section for charged-current neutrino and anti-neutrino scattering from iron. Structure functions, F_2(x,Q^2) and xF_3(x,Q^2), are determined by fitting the inelasticity, y, dependence of the cross sections. This measurement has significantly improved systematic precision as a consequence of more precise understanding of hadron and muon energy scales.

Measurements of the ratio of Drell-Yan yields from an 800 \rm{GeV/c} proton beam incident on liquid hydrogen and deuterium targets are reported. Approximately 360,000 Drell-Yan muon pairs remained after all cuts on the data. From these data, the ratio of anti-down ($\bar{d}$) to anti-up ($\bar{u}$) quark distributions in the proton sea is determined over a wide range in Bjorken-$x$. These results confirm previous measurements by E866 and extend them to lower $x$. From these data, $(\bar{d}-\bar{u})$ and $\int(\bar{d}-\bar{u})dx$ are evaluated for $0.015&lt;x&lt;0.35$. These results are compared with parameterizations of various parton distribution functions, models and experimental results from NA51, NMC, and HERMES.

We present a measurement of the polarization observed for bottomonium states produced in p-Cu collisions at sqrt(s)=38.8 GeV. The angular distribution of the decay dimuons of the Upsilon(1S) state show no polarization at small xF and pT but significant positive transverse production polarization for either pT > 1.8 GeV/c or for xF > 0.35. The Upsilon(2S+3S) unresolved states show a large transverse production polarization at all values of xF and pT measured. These observations are compared with an NRQCD calculation that predicts a transverse polarization in bottomonium production arising from quark-antiquark fusion and gluon-gluon fusion diagrams.

A precise measurement of the ratios of the Drell-Yan cross section per nucleon for an 800 GeV/c proton beam incident on Be, Fe and W targets is reported. The behavior of the Drell-Yan ratios at small target parton momentum fraction is well described by an existing fit to the shadowing observed in deep-inelastic scattering. The cross section ratios as a function of the incident parton momentum fraction set tight limits on the energy loss of quarks passing through a cold nucleus.

A precise measurement of the ratio of Drell-Yan yields from an 800 GeV/c proton beam incident on hydrogen and deuterium targets is reported. Over 140,000 Drell-Yan muon pairs with dimuon mass M_{mu+ mu-} >= 4.5 GeV/c^2 were recorded. From these data, the ratio of anti-down (dbar) to anti-up (ubar) quark distributions in the proton sea is determined over a wide range in Bjorken-x. A strong x dependence is observed in the ratio dbar/ubar, showing substantial enhancement of dbar with respect to ubar for x<0.2. This result is in fair agreement with recent parton distribution parameterizations of the sea. For x>0.2, the observed dbar/ubar ratio is much nearer unity than given by the parameterizations.

We report a high statistics measurement of Upsilon production with an 800 GeV/c proton beam on hydrogen and deuterium targets. The dominance of the gluon-gluon fusion process for Upsilon production at this energy implies that the cross section ratio, $\sigma (p + d \to \Upsilon) / 2\sigma (p + p\to \Upsilon)$, is sensitive to the gluon content in the neutron relative to that in the proton. Over the kinematic region 0 < x_F < 0.6, this ratio is found to be consistent with unity, in striking contrast to the behavior of the Drell-Yan cross section ratio $\sigma(p+d)_{DY}/2\sigma(p+p)_{DY}$. This result shows that the gluon distributions in the proton and neutron are very similar. The Upsilon production cross sections are also compared with the p+d and p+Cu cross sections from earlier measurements.

Measurements of the suppression of the yield per nucleon of J/Psi and Psi' production for 800 GeV/c protons incident on heavy relative to light nuclear targets have been made with very broad coverage in xF and pT. The observed suppression is smallest at xF values of 0.25 and below and increases at larger values of xF. It is also strongest at small pT. Substantial differences between the Psi' and J/Psi are observed for the first time in p-A collisions. The suppression for the Psi' is stronger than that for the J/Psi for xF near zero, but becomes comparable to that for the J/Psi for xF > 0.6.

Results on two-particle $\Delta\eta\Delta\phi$ correlations in inelastic p+p interactions at 20, 31, 40, 80, and 158~GeV/c are presented. The measurements were performed using the large acceptance NA61/SHINE hadron spectrometer at the CERN Super Proton Synchrotron. The data show structures which can be attributed mainly to effects of resonance decays, momentum conservation, and quantum statistics. The results are compared with the EPOS and UrQMD models.

Measurements of inclusive spectra and mean multiplicities of $\pi^\pm$, K$^\pm$, p and $\bar{\textrm{p}}$ produced in inelastic p+p interactions at incident projectile momenta of 20, 31, 40, 80 and 158 GeV/c ($\sqrt{s} = $ 6.3, 7.7, 8.8, 12.3 and 17.3 GeV, respectively) were performed at the CERN Super Proton Synchrotron using the large acceptance NA61/SHINE hadron spectrometer. Spectra are presented as function of rapidity and transverse momentum and are compared to predictions of current models. The measurements serve as the baseline in the NA61/SHINE study of the properties of the onset of deconfinement and search for the critical point of strongly interacting matter.

The measurement of $K^{*}(892)^0$ resonance production via its $K^{+}\pi^{-}$ decay mode in inelastic p+p collisions at beam momentum 158 GeV/$c$ ($\sqrt{s_{NN}}=17.3$ GeV) is presented. The data were recorded by the NA61/SHINE hadron spectrometer at the CERN Super Proton Synchrotron. The $\textit{template}$ method was used to extract the $K^{*}(892)^0$ signal and double differential transverse momentum and rapidity spectra were obtained. The full phase-space mean multiplicity of $K^{*}(892)^0$ mesons was found to be $(78.44 \pm 0.38 \mathrm{(stat)} \pm 6.0 \mathrm{(sys)) \cdot 10^{-3}}$. The NA61/SHINE results are compared with the E$_{POS}$1.99 and Hadron Resonance Gas models as well as with world data from p+p and nucleus-nucleus collisions.

Pentaquark states have been extensively investigated theoretically in the context of the constituent quark model. In this paper experimental searches in the $Xi^{-}\pi^{-}$, $Xi^{-}\pi^{+}$, $Xi^{+}\pi^{-}$ and $Xi^{+}\pi^{+}$ invariant mass spectra in proton-proton interactions at $\sqrt{s}$=17.3 GeV are presented. Previous possible evidence from the NA49 collaboration of the existence of a narrow $Xi^{-}\pi^{-}$ baryon resonance in p+p interactions is not confirmed with almost 10 times greater event statistics. The search was performed using the NA61/SHINE detector which reuses the main components of the NA49 apparatus. No signal was observed with either the selection cuts of NA49 or a newly optimised set.

The production cross section of 30.92 GeV/$c$ protons on carbon is measured by the NA61/SHINE spectrometer at the CERN SPS by means of beam attenuation in a copy (replica) of the 90-cm-long target of the T2K neutrino oscillation experiment. The employed method for direct production cross-section estimation minimizes model corrections for elastic and quasi-elastic interactions. The obtained production cross section is $\sigma_\mathrm{prod}~=~227.6~\pm~0.8\mathrm{(stat)}~_{-~3.2}^{+~1.9}\mathrm{(sys)}~{-~0.8}\mathrm{(mod)}$ mb. It is in agreement with previous NA61/SHINE results obtained with a thin carbon target, while providing improved precision with a total fractional uncertainty of less than 2$\%$. This direct measurement is performed to reduce the uncertainty on the T2K neutrino flux prediction associated with the re-weighting of the interaction rate of neutrino-yielding hadrons.

This paper presents several measurements of total production cross sections and total inelastic cross sections for the following reactions: $\pi^{+}$+C, $\pi^{+}$+Al, $K^{+}$+C, $K^{+}$+Al at 60 GeV/c, $\pi^{+}$+C and $\pi^{+}$+Al at 31 GeV/c . The measurements were made using the NA61/SHINE spectrometer at the CERN SPS. Comparisons with previous measurements are given and good agreement is seen. These interaction cross sections measurements are a key ingredient for neutrino flux prediction from the reinteractions of secondary hadrons in current and future accelerator-based long-baseline neutrino experiments.

The physics goal of the strong interaction program of the NA61/SHINE experiment at the CERN Super Proton Synchrotron (SPS) is to study the phase diagram of hadronic matter by a scan of particle production in collisions of nuclei with various sizes at a set of energies covering the SPS energy range. This paper presents differential inclusive spectra of transverse momentum, transverse mass and rapidity of $\pi^{-}$ mesons produced in $central$ ${}^{40}$Ar+${}^{45}$Sc collisions at beam momenta of 13$A$, 19$A$, 30$A$, 40$A$, 75$A$ and 150$A$ GeV/$c$. Energy and system size dependence of parameters of these distributions -- mean transverse mass, the inverse slope parameter of transverse mass spectra, width of the rapidity distribution and mean multiplicity -- are presented and discussed. Furthermore, the dependence of the ratio of the mean number of produced pions to the mean number of wounded nucleons on the collision energy was derived. The results are compared to predictions of several models.

The NA61/SHINE experiment at the CERN Super Proton Synchrotron (SPS) studies the onset of deconfinement in hadron matter by a scan of particle production in collisions of nuclei with various sizes at a set of energies covering the SPS energy range. This paper presents results on inclusive double-differential spectra, transverse momentum and rapidity distributions and mean multiplicities of $\pi^\pm$, $K^\pm$, $p$ and $\bar{p}$ produced in the 20$\%$ most $central$$^7$Be+$^9$Be collisions at beam momenta of 19$A$, 30$A$, 40$A$, 75$A$ and 150$A$ GeV/$c$. The energy dependence of the $K^\pm$/$\pi^\pm$ ratios as well as of inverse slope parameters of the $K^\pm$ transverse mass distributions are close to those found in inelastic $p$+$p$ reactions. The new results are compared to the world data on $p$+$p$ and Pb+Pb collisions as well as to predictions of the EPOS, UrQMD, AMPT, PHSD and SMASH models.

This paper presents measurements of production cross sections and inelastic cross sections for the following reactions: 60 GeV/$c$ protons with C, Be, Al targets and 120 GeV/$c$ protons with C and Be targets. The analysis was performed using the NA61/SHINE spectrometer at the CERN SPS. First measurements were obtained using protons at 120 GeV/$c$, while the results for protons at 60 GeV/$c$ were compared with previously published measurements. These interaction cross section measurements are critical inputs for neutrino flux prediction in current and future accelerator-based long-baseline neutrino experiments.

Precise knowledge of hadron production rates in the generation of neutrino beams is necessary for accelerator-based neutrino experiments to achieve their physics goals. NA61/SHINE, a large-acceptance hadron spectrometer, has recorded hadron+nucleus interactions relevant to ongoing and future long-baseline neutrino experiments at Fermi National Accelerator Laboratory. This paper presents three analyses of interactions of 60 GeV/$c$ $\pi^+$ with thin, fixed carbon and beryllium targets. Integrated production and inelastic cross sections were measured for both of these reactions. In an analysis of strange, neutral hadron production, differential production multiplicities of $K^0_{S}$, $\Lambda$ and anti-$\Lambda$ were measured. Lastly, in an analysis of charged hadron production, differential production multiplicities of $\pi^+$, $\pi^-$, $K^+$, $K^-$ and protons were measured. These measurements will enable long-baseline neutrino experiments to better constrain predictions of their neutrino flux in order to achieve better precision on their neutrino cross section and oscillation measurements.

Results on $\phi$ meson production in inelastic p+p collisions at CERN SPS energies are presented. They are derived from data collected by the NA61/SHINE fixed target experiment, by means of invariant mass spectra fits in the $\phi \to K^+K^-$ decay channel. They include the first ever measured double differential spectra of $\phi$ mesons as a function of rapidity $y$ and transverse momentum $p_T$ for proton beam momenta of 80 GeV/c and 158 GeV/c, as well as single differential spectra of $y$ or $p_T$ for beam momentum of 40 GeV/c. The corresponding total $\phi$ yields per inelastic p+p event are obtained. These results are compared with existing data on $\phi$ meson production in p+p collisions. The comparison shows consistency but superior accuracy of the present measurements. The emission of $\phi$ mesons in p+p reactions is confronted with that occurring in Pb+Pb collisions, and the experimental results are compared with model predictions. It appears that none of the considered models can properly describe all the experimental observables.

Measurements of the $\pi^{\pm}$, $K^{\pm}$, and proton double differential yields emitted from the surface of the 90-cm-long carbon target (T2K replica) were performed for the incoming 31 GeV/c protons with the NA61/SHINE spectrometer at the CERN SPS using data collected during 2010 run. The double differential $\pi^{\pm}$ yields were measured with increased precision compared to the previously published NA61/SHINE results, while the $K^{\pm}$ and proton yields were obtained for the first time. A strategy for dealing with the dependence of the results on the incoming proton beam profile is proposed. The purpose of these measurements is to reduce significantly the (anti)neutrino flux uncertainty in the T2K long-baseline neutrino experiment by constraining the production of (anti)neutrino ancestors coming from the T2K target.

The production of $\Xi(1321)^{-}$ and $\overline{\Xi}(1321)^{+}$ hyperons in inelastic p+p interactions is studied in a fixed target experiment at a beam momentum of 158 GeV/textitc. Double differential distributions in rapidity y and transverse momentum $p_{T}$ are obtained from a sample of 33M inelastic events. They allow to extrapolate the spectra to full phase space and to determine the mean multiplicity of both $\Xi^{-}$ and $\overline{\Xi}^{+}$. The rapidity and transverse momentum spectra are compared to transport model predictions. The $\Xi^{-}$ mean multiplicity in inelastic p+p interactions at 158~\GeVc is used to quantify the strangeness enhancement in A+A collisions at the same centre-of-mass energy per nucleon pair.

The NA61/SHINE collaboration studies at the CERN Super Proton Synchrotron (SPS) the onset of deconfinement in hadronic matter by the measurement of particle production in collisions of nuclei with various sizes at a set of energies covering the SPS energy range. This paper presents results on inclusive double-differential spectra and mean multiplicities of $\pi^{-}$ mesons produced in the 5\% most \textit{central} $^7$Be+$^9$Be collisions at beam momenta of 19$A$, 30$A$, 40$A$, 75$A$ and 150$A$ GeV/$c$ obtained by the so-called $h^-$ method which does not require any particle identification. The shape of the transverse mass spectra differs from the shapes measured in central Pb+Pb collisions and inelastic p+p interactions. The normalized width of the rapidity distribution decreases with increasing collision energy and is in between the results for inelastic nucleon-nucleon and central Pb+Pb collisions. The mean multiplicity of pions per wounded nucleon in \textit{central} $^7$Be+$^9$Be collisions is close to that in central Pb+Pb collisions up to 75$A$GeV/$c$. However, at the top SPS energy the result is close to the one for inelastic nucleon-nucleon interactions. The results are discussed in the context of predictions for the onset of deconfinement at the CERN SPS collision energies.

The NA61/SHINE experiment at the CERN SPS is performing a uniqe study of the phase diagram of strongly interacting matter by varying collision energy and nuclear mass number of colliding nuclei. In central Pb+Pb collisions the NA49 experiment found structures in the energy dependence of several observables in the CERN SPS energy range that had been predicted for the transition to a deconfined phase. New measurements of NA61/SHINE find intriguing similarities in p+p interactions for which no deconfinement transition is expected at SPS energies. Possible implications will be discussed.

The largest sample ever recorded of $\numub$ charged-current quasi-elastic (CCQE, $\numub + p \to \mup + n$) candidate events is used to produce the minimally model-dependent, flux-integrated double-differential cross section $\frac{d^{2}\sigma}{dT_\mu d\uz}$ for $\numub$ incident on mineral oil. This measurement exploits the unprecedented statistics of the MiniBooNE anti-neutrino mode sample and provides the most complete information of this process to date. Also given to facilitate historical comparisons are the flux-unfolded total cross section $\sigma(E_\nu)$ and single-differential cross section $\frac{d\sigma}{d\qsq}$ on both mineral oil and on carbon by subtracting the $\numub$ CCQE events on hydrogen. The observed cross section is somewhat higher than the predicted cross section from a model assuming independently-acting nucleons in carbon with canonical form factor values. The shape of the data are also discrepant with this model. These results have implications for intra-nuclear processes and can help constrain signal and background processes for future neutrino oscillation measurements.

The MiniBooNE experiment at Fermilab reports a total excess of $638.0 \pm 132.8$ electron-like events ($4.8 \sigma$) from a data sample corresponding to $18.75 \times 10^{20}$ protons-on-target in neutrino mode, which is a 46\% increase in the data sample with respect to previously published results, and $11.27 \times 10^{20}$ protons-on-target in antineutrino mode. The additional statistics allow several studies to address questions on the source of the excess. First, we provide two-dimensional plots in visible energy and cosine of the angle of the outgoing lepton, which can provide valuable input to models for the event excess. Second, we test whether the excess may arise from photons that enter the detector from external events or photons exiting the detector from $\pi^0$ decays in two model independent ways. Beam timing information shows that almost all of the excess is in time with neutrinos that interact in the detector. The radius distribution shows that the excess is distributed throughout the volume, while tighter cuts on the fiducal volume increase the significance of the excess. We conclude that models of the event excess based on entering and exiting photons are disfavored.

The MiniBooNE Collaboration reports first results of a search for $\nu_e$ appearance in a $\nu_\mu$ beam. With two largely independent analyses, we observe no significant excess of events above background for reconstructed neutrino energies above 475 MeV. The data are consistent with no oscillations within a two neutrino appearance-only oscillation model.

We present a measurement of the $\nu_e$-interaction rate in the MicroBooNE detector that addresses the observed MiniBooNE anomalous low-energy excess (LEE). The approach taken isolates neutrino interactions consistent with the kinematics of charged-current quasi-elastic (CCQE) events. The topology of such signal events has a final state with 1 electron, 1 proton, and 0 mesons ($1e1p$). Multiple novel techniques are employed to identify a $1e1p$ final state, including particle identification that use two methods of deep-learning-based image identification, and event isolation using a boosted decision-tree ensemble trained to recognize two-body scattering kinematics. This analysis selects 25 $\nu_e$-candidate events in the reconstructed neutrino energy range of 200--1200 MeV, while $29.0 \pm 1.9_\text{(sys)} \pm 5.4_\text{(stat)}$ are predicted when using $\nu_\mu$ CCQE interactions as a constraint. We use a simplified model to translate the MiniBooNE LEE observation into a prediction for a $\nu_e$ signal in MicroBooNE. A $\Delta \chi^2$ test statistic, based on the combined Neyman--Pearson $\chi^2$ formalism, is used to define frequentist confidence intervals for the LEE signal strength. Using this technique, in the case of no LEE signal, we expect this analysis to exclude a normalization factor of 0.75 (0.98) times the median MiniBooNE LEE signal strength at 90% ($2\sigma$) confidence level, while the MicroBooNE data yield an exclusion of 0.25 (0.38) times the median MiniBooNE LEE signal strength at 90% ($2\sigma$) confidence

We report results from a search for neutrino-induced neutral current (NC) resonant $\Delta$(1232) baryon production followed by $\Delta$ radiative decay, with a $\langle0.8\rangle$~GeV neutrino beam. Data corresponding to MicroBooNE's first three years of operations (6.80$\times$10$^{20}$ protons on target) are used to select single-photon events with one or zero protons and without charged leptons in the final state ($1\gamma1p$ and $1\gamma0p$, respectively). The background is constrained via an in-situ high-purity measurement of NC $\pi^0$ events, made possible via dedicated $2\gamma1p$ and $2\gamma0p$ selections. A total of 16 and 153 events are observed for the $1\gamma1p$ and $1\gamma0p$ selections, respectively, compared to a constrained background prediction of $20.5 \pm 3.65 \text{(sys.)} $ and $145.1 \pm 13.8 \text{(sys.)} $ events. The data lead to a bound on an anomalous enhancement of the normalization of NC $\Delta$ radiative decay of less than $2.3$ times the predicted nominal rate for this process at the 90% confidence level (CL). The measurement disfavors a candidate photon interpretation of the MiniBooNE low-energy excess as a factor of $3.18$ times the nominal NC $\Delta$ radiative decay rate at the 94.8% CL, in favor of the nominal prediction, and represents a greater than $50$-fold improvement over the world's best limit on single-photon production in NC interactions in the sub-GeV neutrino energy range

We report a measurement of the energy-dependent total charged-current cross section $\sigma\left(E_\nu\right)$ for inclusive muon neutrinos scattering on argon, as well as measurements of flux-averaged differential cross sections as a function of muon energy and hadronic energy transfer ($\nu$). Data corresponding to 5.3$\times$10$^{19}$ protons on target of exposure were collected using the MicroBooNE liquid argon time projection chamber located in the Fermilab Booster Neutrino Beam with a mean neutrino energy of approximately 0.8~GeV. The mapping between the true neutrino energy $E_\nu$ and reconstructed neutrino energy $E^{rec}_\nu$ and between the energy transfer $\nu$ and reconstructed hadronic energy $E^{rec}_{had}$ are validated by comparing the data and Monte Carlo (MC) predictions. In particular, the modeling of the missing hadronic energy and its associated uncertainties are verified by a new method that compares the $E^{rec}_{had}$ distributions between data and an MC prediction after constraining the reconstructed muon kinematic distributions, energy and polar angle, to those of data. The success of this validation gives confidence that the missing energy in the MicroBooNE detector is well-modeled and underpins first-time measurements of both the total cross section $\sigma\left(E_\nu\right)$ and the differential cross section $d\sigma/d\nu$ on argon.

We report a search for an anomalous excess of inclusive charged-current (CC) $\nu_e$ interactions using the Wire-Cell event reconstruction package in the MicroBooNE experiment, which is motivated by the previous observation of a low-energy excess (LEE) of electromagnetic events from the MiniBooNE experiment. With a single liquid argon time projection chamber detector, the measurements of $\nu_{\mu}$ CC interactions as well as $\pi^0$ interactions are used to constrain signal and background predictions of $\nu_e$ CC interactions. A data set collected from February 2016 to July 2018 corresponding to an exposure of 6.369 $\times$ 10$^{20}$ protons on target from the Booster Neutrino Beam at FNAL is analyzed. With $x$ representing an overall normalization factor and referred to as the LEE strength parameter, we select 56 fully contained $\nu_e$ CC candidates while expecting 69.6 $\pm$ 8.0 (stat.) $\pm$ 5.0 (sys.) and 103.8 $\pm$ 9.0 (stat.) $\pm$ 7.4 (sys.) candidates after constraints for the absence (eLEE$_{x=0}$) of the median signal strength derived from the MiniBooNE observation and the presence (eLEE$_{x=1}$) of that signal strength, respectively. Under a nested hypothesis test using both rate and shape information in all available channels, the best-fit $x$ is determined to be 0 (eLEE$_{x=0}$) with a 95.5% confidence level upper limit of $x$ at 0.502. Under a simple-vs-simple hypotheses test, the eLEE$_{x=1}$ hypothesis is rejected at 3.75$\sigma$, while the eLEE$_{x=0}$ hypothesis is shown to be consistent with the observation at 0.45$\sigma$. In the context of the eLEE model, the estimated 68.3% confidence interval of the $\nu_e$ hypothesis to explain the LEE observed in the MiniBooNE experiment is disfavored at a significance level of more than 2.6$\sigma$ (3.0$\sigma$) considering MiniBooNE's full (statistical) uncertainties.

The fixed-target MIPP experiment, Fermilab E907, was designed to measure the production of hadrons from the collisions of hadrons of momenta ranging from 5 to 120 GeV/c on a variety of nuclei. These data will generally improve the simulation of particle detectors and predictions of particle beam fluxes at accelerators. The spectrometer momentum resolution is between 3 and 4%, and particle identification is performed for particles ranging between 0.3 and 80 GeV/c using $dE/dx$, time-of-flight and Cherenkov radiation measurements. MIPP collected $1.42 \times10^6$ events of 120 GeV Main Injector protons striking a target used in the NuMI facility at Fermilab. The data have been analyzed and we present here charged pion yields per proton-on-target determined in bins of longitudinal and transverse momentum between 0.5 and 80 GeV/c, with combined statistical and systematic relative uncertainties between 5 and 10%.