Evgueni Tcherniaev

Geant4 is a software toolkit for the simulation of the passage of particles through matter. It is used by a large number of experiments and projects in a variety of application domains, including high energy physics, astrophysics and space science, medical physics and radiation protection. Its functionality and modeling capabilities continue to be extended, while its performance is enhanced. An overview of recent developments in diverse areas of the toolkit is presented. These include performance optimization for complex setups; improvements for the propagation in fields; new options for event biasing; and additions and improvements in geometry, physics processes and interactive capabilities

The TOTEM collaboration has measured the proton–proton total cross section at $$\sqrt{s}=13~\hbox {TeV}$$ with a luminosity-independent method. Using dedicated $$\beta ^{*}=90~\hbox {m}$$ beam optics, the Roman Pots were inserted very close to the beam. The inelastic scattering rate has been measured by the T1 and T2 telescopes during the same LHC fill. After applying the optical theorem the total proton–proton cross section is $$\sigma _\mathrm{tot}=(110.6~\pm ~3.4$$ ) mb, well in agreement with the extrapolation from lower energies. This method also allows one to derive the luminosity-independent elastic and inelastic cross sections: $$\sigma _\mathrm{el}=(31.0~\pm ~1.7)~\hbox {mb}$$ and $$\sigma _\mathrm{inel}=(79.5~\pm ~1.8)~\hbox {mb}$$ .

Abstract The TOTEM experiment at the LHC has performed the first measurement at $$\sqrt{s} = 13\,\mathrm{TeV}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:msqrt><mml:mi>s</mml:mi></mml:msqrt><mml:mo>=</mml:mo><mml:mn>13</mml:mn><mml:mspace /><mml:mi>TeV</mml:mi></mml:mrow></mml:math> of the $$\rho $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>ρ</mml:mi></mml:math> parameter, the real to imaginary ratio of the nuclear elastic scattering amplitude at $$t=0$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:math> , obtaining the following results: $$\rho = 0.09 \pm 0.01$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>ρ</mml:mi><mml:mo>=</mml:mo><mml:mn>0.09</mml:mn><mml:mo>±</mml:mo><mml:mn>0.01</mml:mn></mml:mrow></mml:math> and $$\rho = 0.10 \pm 0.01$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>ρ</mml:mi><mml:mo>=</mml:mo><mml:mn>0.10</mml:mn><mml:mo>±</mml:mo><mml:mn>0.01</mml:mn></mml:mrow></mml:math> , depending on different physics assumptions and mathematical modelling. The unprecedented precision of the $$\rho $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>ρ</mml:mi></mml:math> measurement, combined with the TOTEM total cross-section measurements in an energy range larger than $$10\,\mathrm{TeV}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mn>10</mml:mn><mml:mspace /><mml:mi>TeV</mml:mi></mml:mrow></mml:math> (from 2.76 to $$13\,\mathrm{TeV}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mn>13</mml:mn><mml:mspace /><mml:mi>TeV</mml:mi></mml:mrow></mml:math> ), has implied the exclusion of all the models classified and published by COMPETE. The $$\rho $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>ρ</mml:mi></mml:math> results obtained by TOTEM are compatible with the predictions, from other theoretical models both in the Regge-like framework and in the QCD framework, of a crossing-odd colourless 3-gluon compound state exchange in the t -channel of the proton–proton elastic scattering. On the contrary, if shown that the crossing-odd 3-gluon compound state t -channel exchange is not of importance for the description of elastic scattering, the $$\rho $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>ρ</mml:mi></mml:math> value determined by TOTEM would represent a first evidence of a slowing down of the total cross-section growth at higher energies. The very low-| t | reach allowed also to determine the absolute normalisation using the Coulomb amplitude for the first time at the LHC and obtain a new total proton–proton cross-section measurement $$\sigma _{\mathrm{tot}} = (110.3 \pm 3.5)\,\mathrm{mb}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:msub><mml:mi>σ</mml:mi><mml:mi>tot</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mn>110.3</mml:mn><mml:mo>±</mml:mo><mml:mn>3.5</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mspace /><mml:mi>mb</mml:mi></mml:mrow></mml:math> , completely independent from the previous TOTEM determination. Combining the two TOTEM results yields $$\sigma _{\mathrm{tot}} = (110.5 \pm 2.4)\,\mathrm{mb}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:msub><mml:mi>σ</mml:mi><mml:mi>tot</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mn>110.5</mml:mn><mml:mo>±</mml:mo><mml:mn>2.4</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mspace /><mml:mi>mb</mml:mi></mml:mrow></mml:math> .

Abstract The TOTEM collaboration has measured the elastic proton-proton differential cross section $$\mathrm{d}\sigma /\mathrm{d}t$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>d</mml:mi><mml:mi>σ</mml:mi><mml:mo>/</mml:mo><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:math> at $$\sqrt{s}=13$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:msqrt><mml:mi>s</mml:mi></mml:msqrt><mml:mo>=</mml:mo><mml:mn>13</mml:mn></mml:mrow></mml:math> TeV LHC energy using dedicated $$\beta ^{*}=90$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:msup><mml:mi>β</mml:mi><mml:mrow><mml:mrow /><mml:mo>∗</mml:mo></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:mn>90</mml:mn></mml:mrow></mml:math> m beam optics. The Roman Pot detectors were inserted to 10 $$\sigma $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>σ</mml:mi></mml:math> distance from the LHC beam, which allowed the measurement of the range [0.04 GeV $$^{2}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msup><mml:mrow /><mml:mn>2</mml:mn></mml:msup></mml:math> ; 4 GeV $$^{2}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msup><mml:mrow /><mml:mn>2</mml:mn></mml:msup></mml:math> $$]$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mo>]</mml:mo></mml:math> in four-momentum transfer squared | t |. The efficient data acquisition allowed to collect about 10 $$^{9}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msup><mml:mrow /><mml:mn>9</mml:mn></mml:msup></mml:math> elastic events to precisely measure the differential cross-section including the diffractive minimum (dip), the subsequent maximum (bump) and the large-| t | tail. The average nuclear slope has been found to be $$B=(20.40 \pm 0.002^{\mathrm{stat}} \pm 0.01^{\mathrm{syst}})~$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>B</mml:mi><mml:mo>=</mml:mo><mml:mo>(</mml:mo><mml:mn>20.40</mml:mn><mml:mo>±</mml:mo><mml:mn>0</mml:mn><mml:mo>.</mml:mo><mml:msup><mml:mn>002</mml:mn><mml:mi>stat</mml:mi></mml:msup><mml:mo>±</mml:mo><mml:mn>0</mml:mn><mml:mo>.</mml:mo><mml:msup><mml:mn>01</mml:mn><mml:mi>syst</mml:mi></mml:msup><mml:mo>)</mml:mo><mml:mspace /></mml:mrow></mml:math> GeV $$^{-2}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msup><mml:mrow /><mml:mrow><mml:mo>-</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:math> in the | t |-range 0.04–0.2 GeV $$^{2}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msup><mml:mrow /><mml:mn>2</mml:mn></mml:msup></mml:math> . The dip position is $$|t_{\mathrm{dip}}|=(0.47 \pm 0.004^{\mathrm{stat}} \pm 0.01^{\mathrm{syst}})~$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mrow><mml:mo>|</mml:mo></mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mi>dip</mml:mi></mml:msub><mml:mrow><mml:mo>|</mml:mo><mml:mo>=</mml:mo></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mn>0.47</mml:mn><mml:mo>±</mml:mo><mml:mn>0</mml:mn><mml:mo>.</mml:mo><mml:msup><mml:mn>004</mml:mn><mml:mi>stat</mml:mi></mml:msup><mml:mo>±</mml:mo><mml:mn>0</mml:mn><mml:mo>.</mml:mo><mml:msup><mml:mn>01</mml:mn><mml:mi>syst</mml:mi></mml:msup><mml:mo>)</mml:mo></mml:mrow><mml:mspace /></mml:mrow></mml:math> GeV $$^{2}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msup><mml:mrow /><mml:mn>2</mml:mn></mml:msup></mml:math> . The differential cross section ratio at the bump vs. at the dip $$R=1.77\pm 0.01^{\mathrm{stat}}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>R</mml:mi><mml:mo>=</mml:mo><mml:mn>1.77</mml:mn><mml:mo>±</mml:mo><mml:mn>0</mml:mn><mml:mo>.</mml:mo><mml:msup><mml:mn>01</mml:mn><mml:mi>stat</mml:mi></mml:msup></mml:mrow></mml:math> has been measured with high precision. The series of TOTEM elastic pp measurements show that the dip is a permanent feature of the pp differential cross-section at the TeV scale.

Abstract The proton–proton elastic differential cross section $${\mathrm{d}}\sigma /{\mathrm{d}}t$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>d</mml:mi><mml:mi>σ</mml:mi><mml:mo>/</mml:mo><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:math> has been measured by the TOTEM experiment at $$\sqrt{s}=2.76\hbox { TeV}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:msqrt><mml:mi>s</mml:mi></mml:msqrt><mml:mo>=</mml:mo><mml:mn>2.76</mml:mn><mml:mspace /><mml:mtext>TeV</mml:mtext></mml:mrow></mml:math> energy with $$\beta ^{*}=11\hbox { m}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:msup><mml:mi>β</mml:mi><mml:mrow><mml:mrow /><mml:mo>∗</mml:mo></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:mn>11</mml:mn><mml:mspace /><mml:mtext>m</mml:mtext></mml:mrow></mml:math> beam optics. The Roman Pots were inserted to 13 times the transverse beam size from the beam, which allowed to measure the differential cross-section of elastic scattering in a range of the squared four-momentum transfer (| t |) from 0.36 to $$0.74\hbox { GeV}^{2}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mn>0.74</mml:mn><mml:mspace /><mml:msup><mml:mtext>GeV</mml:mtext><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:math> . The differential cross-section can be described with an exponential in the | t |-range between 0.36 and $$0.54\hbox { GeV}^{2}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mn>0.54</mml:mn><mml:mspace /><mml:msup><mml:mtext>GeV</mml:mtext><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:math> , followed by a diffractive minimum (dip) at $$|t_{\mathrm{dip}}|=(0.61\pm 0.03)\hbox { GeV}^{2}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mrow><mml:mo>|</mml:mo></mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mi>dip</mml:mi></mml:msub><mml:mrow><mml:mo>|</mml:mo><mml:mo>=</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mn>0.61</mml:mn><mml:mo>±</mml:mo><mml:mn>0.03</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mspace /></mml:mrow><mml:msup><mml:mtext>GeV</mml:mtext><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:math> and a subsequent maximum (bump). The ratio of the $${\mathrm{d}}\sigma /{\mathrm{d}}t$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>d</mml:mi><mml:mi>σ</mml:mi><mml:mo>/</mml:mo><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:math> at the bump and at the dip is $$1.7\pm 0.2$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mn>1.7</mml:mn><mml:mo>±</mml:mo><mml:mn>0.2</mml:mn></mml:mrow></mml:math> . When compared to the proton–antiproton measurement of the D0 experiment at $$\sqrt{s} = 1.96\hbox { TeV}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:msqrt><mml:mi>s</mml:mi></mml:msqrt><mml:mo>=</mml:mo><mml:mn>1.96</mml:mn><mml:mspace /><mml:mtext>TeV</mml:mtext></mml:mrow></mml:math> , a significant difference can be observed. Under the condition that the effects due to the energy difference between TOTEM and D0 can be neglected, the result provides evidence for the exchange of a colourless C-odd three-gluon compound state in the t -channel of the proton–proton and proton–antiproton elastic scattering.

The double-differential production cross-section of positive pions, $d^2\sigma^{\pi^{+}}/dpd\Omega$, measured in the HARP experiment is presented. The incident particles are 8.9 GeV/c protons directed onto a beryllium target with a nominal thickness of 5% of a nuclear interaction length. The measured cross-section has a direct impact on the prediction of neutrino fluxes for the MiniBooNE and SciBooNE experiments at Fermilab. After cuts, 13 million protons on target produced about 96,000 reconstructed secondary tracks which were used in this analysis. Cross-section results are presented in the kinematic range 0.75 GeV/c < $p_{\pi}$ < 6.5 GeV/c and 30 mrad < $\theta_{\pi}$ < 210 mrad in the laboratory frame.

A precision measurement of the double-differential production cross-section, d2σπ+/dpdΩ, for pions of positive charge, performed in the HARP experiment is presented. The incident particles are protons of 12.9 GeV/c momentum impinging on an aluminium target of 5% nuclear interaction length. The measurement of this cross-section has a direct application to the calculation of the neutrino flux of the K2K experiment. After cuts, 210 000 secondary tracks reconstructed in the forward spectrometer were used in this analysis. The results are given for secondaries within a momentum range from 0.75 to 6.5 GeV/c, and within an angular range from 30 mrad to 210 mrad. The absolute normalization was performed using prescaled beam triggers counting protons on target. The overall scale of the cross-section is known to better than 6%, while the average point-to-point error is 8.2%.

HARP is a high-statistics, large solid angle experiment to measure hadron production using proton and pion beams with momenta between 1.5 and 15 GeV/c impinging on many different solid and liquid targets from low to high Z. The experiment, located in the T9 beam of the CERN PS, took data in 2001 and 2002. For the measurement of momenta of produced particles and for the identification of particle types, the experiment includes a large-angle spectrometer, based on a Time Projection Chamber and a system of Resistive Plate Chambers, and a forward spectrometer equipped with a set of large drift chambers, a threshold Cherenkov detector, a time-of-flight wall and an electromagnetic calorimeter. The large angle system uses a solenoidal magnet, while the forward spectrometer is based on a dipole magnet. Redundancy in particle identification has been sought, to enable the cross-calibration of efficiencies and to obtain a few percent overall accuracy in the cross-section measurements. Detector construction, operation and initial physics performances are reported. In addition, the full chain for data recording and analysis, from trigger to the software framework, is described.

Measurements of the double-differential π± production cross section in the momentum range 100⩽p⩽800 MeV/c and angle range 0.35⩽θ⩽2.15 rad in proton-beryllium, proton-carbon, proton-aluminium, proton-copper, proton-tin, proton-tantalum, and proton-lead collisions are presented. The data were taken with the large-acceptance HARP detector in the T9 beam line of the CERN PS. The pions were produced by proton beams in a momentum range from 3 to 12.9 GeV/c hitting a target with a thickness of 5% of a nuclear interaction length. Tracking and identification of the produced particles was performed by using a small-radius cylindrical Time Projection Chamber (TPC) placed inside a solenoidal magnet. Incident particles were identified by an elaborate system of beam detectors. Results are obtained for the double-differential cross sections d2σ/(dpdθ) at six incident proton beam momenta [3, 5, 8, and 8.9 GeV/c (Be only) and 12 and 12.9 GeV/c (Al only)]. They are based on a complete correction of static and dynamic distortions of tracks in the HARP TPC, which allows the complete statistics of the collected data set to be used. The results include and supersede our previously published results and are compatible with these. Results are compared with the GEANT4 and MARS Monte Carlo simulation.26 MoreReceived 27 December 2007DOI:https://doi.org/10.1103/PhysRevC.77.055207©2008 American Physical Society

A measurement of the double-differential cross-section for the production of charged pions in proton--tantalum collisions emitted at large angles from the incoming beam direction is presented. The data were taken in 2002 with the HARP detector in the T9 beam line of the CERN PS. The pions were produced by proton beams in a momentum range from 3 \GeVc to 12 \GeVc hitting a tantalum target with a thickness of 5% of a nuclear interaction length. The angular and momentum range covered by the experiment ($100 \MeVc \le p < 800 \MeVc$ and $0.35 \rad \le \theta <2.15 \rad$) is of particular importance for the design of a neutrino factory. The produced particles were detected using a small-radius cylindrical time projection chamber (TPC) placed in a solenoidal magnet. Track recognition, momentum determination and particle identification were all performed based on the measurements made with the TPC. An elaborate system of detectors in the beam line ensured the identification of the incident particles. Results are shown for the double-differential cross-sections ${{\mathrm{d}^2 \sigma}} / {{\mathrm{d}p\mathrm{d}\theta}}$ at four incident proton beam momenta (3 \GeVc, 5 \GeVc, 8 \GeVc and 12 \GeVc). In addition, the pion yields within the acceptance of typical neutrino factory designs are shown as a function of beam momentum. The measurement of these yields within a single experiment eliminates most systematic errors in the comparison between rates at different beam momenta and between positive and negative pion production.

The results of the measurements of the double-differential production cross-sections of pions in p-C and $\pi^\pm$-C interactions using the forward spectrometer of the HARP experiment are presented. The incident particles are 12 GeV/c protons and charged pions directed onto a carbon target with a thickness of 5% of a nuclear interaction length. For p-C interactions the analysis is performed using 100035 reconstructed secondary tracks, while the corresponding numbers of tracks for $\pi^-$-C and $\pi^+$-C analyses are 106534 and 10122 respectively. Cross-section results are presented in the kinematic range 0.5 GeV/c $\leq p_{\pi} <$ 8 GeV/c and 30 mrad $\leq \theta_{\pi} <$ 240 mrad in the laboratory frame. The measured cross-sections have a direct impact on the precise calculation of atmospheric neutrino fluxes and on the improved reliability of extensive air shower simulations by reducing the uncertainties of hadronic interaction models in the low energy range.

A measurement of the double-differential π± production cross-section in proton–carbon, proton–copper and proton–tin collisions in the range of pion momentum 100 MeV/c≤p<800 MeV/c and angle 0.35 rad≤θ<2.15 rad is presented. The data were taken with the HARP detector in the T9 beam line of the CERN PS. The pions were produced by proton beams in a momentum range from 3 GeV/c to 12 GeV/c hitting a target with a thickness of 5% of a nuclear interaction length. The tracking and identification of the produced particles was done using a small-radius cylindrical time projection chamber (TPC) placed in a solenoidal magnet. An elaborate system of detectors in the beam line ensured the identification of the incident particles. Results are shown for the double-differential cross-sections d2σ/dpdθ at four incident proton beam momenta (3 GeV/c, 5 GeV/c, 8 GeV/c and 12 GeV/c).

Measurements of the double-differential charged pion production cross-section in the range of momentum 0.5 GeV/c < p < 8.0 GeV/c and angle 0.025 rad < theta <0.25 rad in collisions of protons on beryllium, carbon, nitrogen, oxygen, aluminium, copper, tin, tantalum and lead are presented. The data were taken with the large acceptance HARP detector in the T9 beam line of the CERN PS. Incident particles were identified by an elaborate system of beam detectors. The data were taken with thin targets of 5% of a nuclear interaction length. The tracking and identification of the produced particles was performed using the forward system of the HARP experiment. Results are obtained for the double-differential cross section mainly at four incident proton beam momenta (3 GeV/c, 5 GeV/c, 8 GeV/c and 12 GeV/c). Measurements are compared with the GEANT4 and MARS Monte Carlo generators. A global parametrization is provided as an approximation of all the collected datasets which can serve as a tool for quick yields estimates.

Abstract The TOTEM collaboration at the CERN LHC has measured the differential cross-section of elastic proton–proton scattering at $$\sqrt{s} = 8\,\mathrm{TeV}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:msqrt> <mml:mi>s</mml:mi> </mml:msqrt> <mml:mo>=</mml:mo> <mml:mn>8</mml:mn> <mml:mspace /> <mml:mi>TeV</mml:mi> </mml:mrow> </mml:math> in the squared four-momentum transfer range $$0.2\,\mathrm{GeV^{2}}&lt; |t| &lt; 1.9\,\mathrm{GeV^{2}}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mn>0.2</mml:mn> <mml:mspace /> <mml:msup> <mml:mi>GeV</mml:mi> <mml:mn>2</mml:mn> </mml:msup> <mml:mo>&lt;</mml:mo> <mml:mrow> <mml:mo>|</mml:mo> <mml:mi>t</mml:mi> <mml:mo>|</mml:mo> </mml:mrow> <mml:mo>&lt;</mml:mo> <mml:mn>1.9</mml:mn> <mml:mspace /> <mml:msup> <mml:mi>GeV</mml:mi> <mml:mn>2</mml:mn> </mml:msup> </mml:mrow> </mml:math> . This interval includes the structure with a diffractive minimum (“dip”) and a secondary maximum (“bump”) that has also been observed at all other LHC energies, where measurements were made. A detailed characterisation of this structure for $$\sqrt{s} = 8\,\mathrm{TeV}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:msqrt> <mml:mi>s</mml:mi> </mml:msqrt> <mml:mo>=</mml:mo> <mml:mn>8</mml:mn> <mml:mspace /> <mml:mi>TeV</mml:mi> </mml:mrow> </mml:math> yields the positions, $$|t|_{\mathrm{dip}} = (0.521 \pm 0.007)\,\mathrm{GeV^2}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mo>|</mml:mo> <mml:mi>t</mml:mi> <mml:mo>|</mml:mo> </mml:mrow> <mml:mi>dip</mml:mi> </mml:msub> <mml:mo>=</mml:mo> <mml:mrow> <mml:mo>(</mml:mo> <mml:mn>0.521</mml:mn> <mml:mo>±</mml:mo> <mml:mn>0.007</mml:mn> <mml:mo>)</mml:mo> </mml:mrow> <mml:mspace /> <mml:msup> <mml:mi>GeV</mml:mi> <mml:mn>2</mml:mn> </mml:msup> </mml:mrow> </mml:math> and $$|t|_{\mathrm{bump}} = (0.695 \pm 0.026)\,\mathrm{GeV^2}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mo>|</mml:mo> <mml:mi>t</mml:mi> <mml:mo>|</mml:mo> </mml:mrow> <mml:mi>bump</mml:mi> </mml:msub> <mml:mo>=</mml:mo> <mml:mrow> <mml:mo>(</mml:mo> <mml:mn>0.695</mml:mn> <mml:mo>±</mml:mo> <mml:mn>0.026</mml:mn> <mml:mo>)</mml:mo> </mml:mrow> <mml:mspace /> <mml:msup> <mml:mi>GeV</mml:mi> <mml:mn>2</mml:mn> </mml:msup> </mml:mrow> </mml:math> , as well as the cross-section values, $$\left. {\mathrm{d}\sigma /\mathrm{d}t}\right| _{\mathrm{dip}} = (15.1 \pm 2.5)\,\mathrm{{\mu b/GeV^2}}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:msub> <mml:mfenced> <mml:mrow> <mml:mi>d</mml:mi> <mml:mi>σ</mml:mi> <mml:mo>/</mml:mo> <mml:mi>d</mml:mi> <mml:mi>t</mml:mi> </mml:mrow> </mml:mfenced> <mml:mi>dip</mml:mi> </mml:msub> <mml:mo>=</mml:mo> <mml:mrow> <mml:mo>(</mml:mo> <mml:mn>15.1</mml:mn> <mml:mo>±</mml:mo> <mml:mn>2.5</mml:mn> <mml:mo>)</mml:mo> </mml:mrow> <mml:mspace /> <mml:mrow> <mml:mi>μ</mml:mi> <mml:mi>b</mml:mi> <mml:mo>/</mml:mo> <mml:msup> <mml:mi>GeV</mml:mi> <mml:mn>2</mml:mn> </mml:msup> </mml:mrow> </mml:mrow> </mml:math> and $$\left. {\mathrm{d}\sigma /\mathrm{d}t}\right| _{\mathrm{bump}} = (29.7 \pm 1.8)\,\mathrm{{\mu b/GeV^2}}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:msub> <mml:mfenced> <mml:mrow> <mml:mi>d</mml:mi> <mml:mi>σ</mml:mi> <mml:mo>/</mml:mo> <mml:mi>d</mml:mi> <mml:mi>t</mml:mi> </mml:mrow> </mml:mfenced> <mml:mi>bump</mml:mi> </mml:msub> <mml:mo>=</mml:mo> <mml:mrow> <mml:mo>(</mml:mo> <mml:mn>29.7</mml:mn> <mml:mo>±</mml:mo> <mml:mn>1.8</mml:mn> <mml:mo>)</mml:mo> </mml:mrow> <mml:mspace /> <mml:mrow> <mml:mi>μ</mml:mi> <mml:mi>b</mml:mi> <mml:mo>/</mml:mo> <mml:msup> <mml:mi>GeV</mml:mi> <mml:mn>2</mml:mn> </mml:msup> </mml:mrow> </mml:mrow> </mml:math> , for the dip and the bump, respectively.

Measurements of the double-differential π± production cross-section in the range of momentum 100 MeV/c≤p< 800 MeV/c and angle 0.35 rad ≤θ< 2.15 rad in proton–beryllium, proton–aluminium and proton–lead collisions are presented. The data were taken with the HARP detector in the T9 beam line of the CERN PS. The pions were produced by proton beams in a momentum range from 3 GeV/c to 12.9 GeV/c hitting a target with a thickness of 5% of a nuclear interaction length. The tracking and identification of the produced particles was performed using a small-radius cylindrical time projection chamber (TPC) placed inside a solenoidal magnet. Incident particles were identified by an elaborate system of beam detectors. Results are obtained for the double-differential cross-sections d2σ/dpdθ at six incident proton beam momenta (3 GeV/c, 5 GeV/c, 8 GeV/c, 8.9 GeV/c (Be only), 12 GeV/c and 12.9 GeV/c (Al only)) and compared to previously available data.

The Geant4 Visualisation System is a multi-driver graphics system designed to serve the Geant4 Simulation Toolkit. It is aimed at the visualisation of Geant4 data, primarily detector descriptions and simulated particle trajectories and hits. It can handle a variety of graphical technologies simultaneously and interchangeably, allowing the user to choose the visual representation most appropriate to requirements. It conforms to the low-level Geant4 abstract graphical user interfaces and introduces new abstract classes from which the various drivers are derived and that can be straightforwardly extended, for example, by the addition of a new driver. It makes use of an extendable class library of models and filters for data representation and selection. The Geant4 Visualisation System supports a rich set of interactive commands based on the Geant4 command system. It is included in the Geant4 code distribution and maintained and documented like other components of Geant4.

Measurements of the double-differential π± production cross-section in the range of momentum 0.5 GeV/c⩽p⩽8.0 GeV/c and angle 0.025rad⩽θ⩽0.25rad in interactions of charged pions on beryllium, carbon, aluminium, copper, tin, tantalum and lead are presented. These data represent the first experimental campaign to systematically measure forward pion hadroproduction. The data were taken with the large acceptance HARP detector in the T9 beam line of the CERN PS. Incident particles, impinging on a 5% nuclear interaction length target, were identified by an elaborate system of beam detectors. The tracking and identification of the produced particles was performed using the forward spectrometer of the HARP detector. Results are obtained for the double-differential cross-sections d2σ/dpdΩ mainly at four incident pion beam momenta (3 GeV/c, 5 GeV/c, 8 GeV/c and 12 GeV/c). The measurements are compared with the GEANT4 and MARS Monte Carlo simulation.

Measurements of the double-differential pi+/- production cross-section in the range of momentum 100 MeV/c <= p <= 800 MeV/c and angle 0.35 rad <= theta <= 2.15 rad using pi+/- beams incident on beryllium, aluminium, carbon, copper, tin, tantalum and lead targets are presented. The data were taken with the large acceptance HARP detector in the T9 beam line of the CERN Proton Synchrotron. The secondary pions were produced by beams in a momentum range from 3 GeV/c to 12.9 GeV/c hitting a solid target with a thickness of 5% of a nuclear interaction length. The tracking and identification of the produced particles was performed using a small-radius cylindrical time projection chamber (TPC) placed inside a solenoidal magnet. Incident particles were identified by an elaborate system of beam detectors. Results are obtained for the double-differential cross-sections d2sigma/dpdtheta at six incident beam momenta. Data at 3 GeV/c, 5 GeV/c, 8 GeV/c, and 12 GeV/c are available for all targets while additional data at 8.9 GeV/c and 12.9 GeV/c were taken in positive particle beams on Be and Al targets, respectively. The measurements are compared with several generators of GEANT4 and the MARS Monte Carlo simulation.

The particle identification (PID) methods used for the calculation of secondary pion yields with the HARP forward spectrometer are presented. Information from time of flight and Cherenkov detectors is combined using likelihood techniques. The efficiencies and purities associated with the different PID selection criteria are obtained from the data. For the proton–aluminium interactions at 12.9 GeV/c incident momentum, the PID efficiencies for positive pions are 86% in the momentum range below 2 GeV/c, 92% between 2 and 3 GeV/c and 98% in the momentum range above 3 GeV/c. The purity of the selection is better than 92% for all momenta. Special emphasis has been put on understanding the main error sources. The final PID uncertainty on the pion yield is 3.3%.

Measurements of double-differential charged pion production cross-sections in interactions of 12 GeV/c protons on O_2 and N_2 thin targets are presented in the kinematic range 0.5 GeV/c < p_{\pi} < 8 GeV/c and 50 mrad < \theta_{\pi} < 250 mrad (in the laboratory frame) and are compared with p--C results. For p--N_2 (p--O_2) interactions the analysis is performed using 38576 (7522) reconstructed secondary pions. The analysis uses the beam instrumentation and the forward spectrometer of the HARP experiment at CERN PS. The measured cross-sections have a direct impact on the precise calculation of atmospheric neutrino fluxes and on the improved reliability of extensive air shower simulations by reducing the uncertainties of hadronic interaction models in the low energy range. In particular, the present results allow the common hypothesis that p--C data can be used to predict the p--N_2 and p--O_2 pion production cross-sections to be tested.

The HARP system of resistive plate chambers (RPCs) was designed to perform particle identification by the measurement of the difference in the time-of-flight of different particles. In previous papers an apparent discrepancy was shown between the response of the RPCs to minimum ionizing pions and heavily ionizing protons. Using the kinematics of elastic scattering off a hydrogen target a controlled beam of low momentum recoil protons was directed onto the RPC chambers. With this method the trajectory and momentum, and hence the time-of-flight of the protons can be precisely predicted without need for a measurement of momentum of the protons. It is demonstrated that the measurement of the time-of-arrival of particles by the thin gas-gap glass RPC system of the HARP experiment depends on the primary ionization deposited by the particle in the detector.

In the HARP experiment the large-angle spectrometer is using a cylindrical TPC as main tracking and particle identification detector. The momentum scale of reconstructed tracks in the TPC is the most important systematic error for the majority of kinematic bins used for the HARP measurements of the double-differential production cross-section of charged pions in proton interactions on nuclear targets at large angle. The HARP TPC operated with a number of hardware shortfalls and operational mistakes. Thus it was important to control and characterize its momentum calibration. While it was not possible to enter a direct particle beam into the sensitive volume of the TPC to calibrate the detector, a set of physical processes and detector properties were exploited to achieve a precise calibration of the apparatus. In the following we recall the main issues concerning the momentum measurement in the HARP TPC, and describe the cross-checks made to validate the momentum scale. As a conclusion, this analysis demonstrates that the measurement of momentum is correct within the published precision of 3%.

Measurements of the double-differential proton production cross-section ${d}^{2}\ensuremath{\sigma}/dpd\ensuremath{\Omega}$ in the range of momentum $0.5 \mathrm{GeV}/c\ensuremath{\leqslant}p&lt;8.0 \mathrm{GeV}/c$ and angle $0.05 \text{rad}\ensuremath{\leqslant}\ensuremath{\theta}&lt;0.25 \text{rad}$ in collisions of charged pions and protons on beryllium, carbon, aluminium, copper, tin, tantalum, and lead are presented. The data were taken with the large acceptance HARP detector in the T9 beam line of the CERN Proton Synchrotron. Incident particles were identified by an elaborate system of beam detectors and impinged on a target of $5%$ of a nuclear interaction length. The tracking and identification of the produced particles was performed using the forward spectrometer of the HARP experiment. Results are obtained for the double-differential cross-sections mainly at four incident beam momenta ($3,5,8$, and $12$ $\text{GeV}/c$). Measurements are compared with predictions of the geant4 and mars Monte Carlo generators.

The TOTEM collaboration has measured the proton-proton total cross section at $\sqrt{s}=13$ TeV with a luminosity-independent method. Using dedicated $β^{*}=90$ m beam optics, the Roman Pots were inserted very close to the beam. The inelastic scattering rate has been measured by the T1 and T2 telescopes during the same LHC fill. After applying the optical theorem the total proton-proton cross section is $σ_{\rm tot}=(110.6 \pm 3.4$) mb, well in agreement with the extrapolation from lower energies. This method also allows one to derive the luminosity-independent elastic and inelastic cross sections: $σ_{\rm el} = (31.0 \pm 1.7)$ mb and $σ_{\rm inel} = (79.5 \pm 1.8)$ mb.

The HARP Collaboration has presented measurements of the double-differential ${\ensuremath{\pi}}^{\ifmmode\pm\else\textpm\fi{}}$ production cross section in the range of momentum $100 \text{MeV}/c\ensuremath{\leqslant}p\ensuremath{\leqslant}800 \text{MeV}/c$ and angle $0.35 \text{rad}\ensuremath{\leqslant}\ensuremath{\theta}\ensuremath{\leqslant}2.15 \text{rad}$ with proton beams hitting thin nuclear targets. In many applications the extrapolation to long targets is necessary. In this article the analysis of data taken with long (one interaction length) solid cylindrical targets made of carbon, tantalum, and lead is presented. The data were taken with the large-acceptance HARP detector in the T9 beam line of the CERN proton synchrotron. The secondary pions were produced by beams of protons with momenta of 5, 8, and $12\text{GeV}/c$. The tracking and identification of the produced particles were performed using a small-radius cylindrical time projection chamber placed inside a solenoidal magnet. Incident protons were identified by an elaborate system of beam detectors. Results are obtained for the double-differential yields per target nucleon ${\mathrm{d}}^{2}\ensuremath{\sigma}/\mathrm{d}p\mathrm{d}\ensuremath{\theta}$. The measurements are compared with predictions of the MARS and GEANT4 Monte Carlo simulations.

Full detector simulation was among the largest CPU consumer in all CERN experiment software stacks for the first two runs of the Large Hadron Collider (LHC). In the early 2010's, the projections were that simulation demands would scale linearly with luminosity increase, compensated only partially by an increase of computing resources. The extension of fast simulation approaches to more use cases, covering a larger fraction of the simulation budget, is only part of the solution due to intrinsic precision limitations. The remainder corresponds to speeding-up the simulation software by several factors, which is out of reach using simple optimizations on the current code base. In this context, the GeantV R&amp;D project was launched, aiming to redesign the legacy particle transport codes in order to make them benefit from fine-grained parallelism features such as vectorization, but also from increased code and data locality. This paper presents extensively the results and achievements of this R&amp;D, as well as the conclusions and lessons learnt from the beta prototype.

Abstract Full detector simulation was among the largest CPU consumers in all CERN experiment software stacks for the first two runs of the Large Hadron Collider. In the early 2010s, it was projected that simulation demands would scale linearly with increasing luminosity, with only partial compensation from increasing computing resources. The extension of fast simulation approaches to cover more use cases that represent a larger fraction of the simulation budget is only part of the solution, because of intrinsic precision limitations. The remainder corresponds to speeding up the simulation software by several factors, which is not achievable by just applying simple optimizations to the current code base. In this context, the GeantV R&amp;D project was launched, aiming to redesign the legacy particle transport code in order to benefit from features of fine-grained parallelism, including vectorization and increased locality of both instruction and data. This paper provides an extensive presentation of the results and achievements of this R&amp;D project, as well as the conclusions and lessons learned from the beta version prototype.

Current Total Skin Electron Therapy (TSET) Stanford technique for cutaneous lymphoma, established in the 70's, involves a unique irradiation setup, i.e. patient's position and beam arrangement, for all patients with ensuing great variability in dose distribution and difficult dose optimization. A Geant4-based simulation has been developed to explore the possibility of personalizing the dose to each patient's anatomy. To achieve this optimization of the treatment method, this project enrolls different aspects of the clinical and computational techniques: starting with the knowledge of the experimental parameters involving TSET practice, passing through an innovative approach to model the patient's anatomy, a precise description of the electron beam and a validated configuration of the physics models handling the interactions of the electrons and of secondary particles. The Geant4-based simulation models the patient as a tessellated solid derived from the optical scan of her/his body, realistically reproduces the irradiation environment in detail and calculates the energy deposition corresponding to each facet of the patient's scanned surface. The resulting three-dimensional dose distribution constitutes the basis for the personalization of the medical treatement as appropriate to each patient's specific characteristics.

The SkinScan project develops computational tools for a challenging radiotherapy technique -total skin irradiation with electrons for the treatment of cutaneous lymphomas. It addresses a complex problem domain, which still rests on a radiation therapy methodology developed in the 80's of the past century, by means of an innovative computational approach that involves a multidisciplinary team of physicists, medical experts, computer scientists and bioengineers. In addition to its intrinsic scientific difficulties, the project has suffered hardships associated with the computational environment, the software development process, the heterogeneous background of the team and unforeseeable external constraints, which have threatened its success. We review the obstacles encountered in the computational development and we analyze in retrospect the effectiveness of the steps taken to address them. We highlight how technical and methodological aspects are intertwined and show how issues in software engineering hindered the achievement of scientific results. The lessons learned from this project represent valuable food for thought for other computational developments. Finally, we present an overview of SkinScan scientific results and identify the key factors which enabled them.

The TOTEM experiment at the LHC has performed the first measurement at $\sqrt{s} = 13$ TeV of the $ρ$ parameter, the real to imaginary ratio of the nuclear elastic scattering amplitude at $t=0$, obtaining the following results: $ρ= 0.09 \pm 0.01$ and $ρ= 0.10 \pm 0.01$, depending on different physics assumptions and mathematical modelling. The unprecedented precision of the $ρ$ measurement, combined with the TOTEM total cross-section measurements in an energy range larger than 10 TeV (from 2.76 to 13 TeV), has implied the exclusion of all the models classified and published by COMPETE. The $ρ$ results obtained by TOTEM are compatible with the predictions, from alternative theoretical models both in the Regge-like framework and in the QCD framework, of a colourless 3-gluon bound state exchange in the $t$-channel of the proton-proton elastic scattering. On the contrary, if shown that the 3-gluon bound state $t$-channel exchange is not of importance for the description of elastic scattering, the $ρ$ value determined by TOTEM would represent a first evidence of a slowing down of the total cross-section growth at higher energies. The very low-$|t|$ reach allowed also to determine the absolute normalisation using the Coulomb amplitude for the first time at the LHC and obtain a new total proton-proton cross-section measurement $σ_{tot} = 110.3 \pm 3.5$ mb, completely independent from the previous TOTEM determination. Combining the two TOTEM results yields $σ_{tot} = 110.5 \pm 2.4$ mb.