Measurement of the residual energy of muons in the Gran Sasso underground Laboratories

The MACRO detector was located in the Hall B of the Gran Sasso underground Laboratories under an average rock overburden of 3700 hg/cm^2. A transition radiation detector composed of three identical modules, covering a total horizontal area of 36 m^2, was installed inside the empty upper part of the detector in order to measure the residual energy of muons. This paper presents the measurement of the residual energy of single and double muons crossing the apparatus. Our data show that double muons are more energetic than single ones. This measurement is performed over a standard rock depth range from 3000 to 6500 hg/cm^2.Comment: 28 pages, 9 figure

Muon Energy Estimate Through Multiple Scattering with the Macro Detector

Muon energy measurement represents an important issue for any experiment addressing neutrino induced upgoing muon studies. Since the neutrino oscillation probability depends on the neutrino energy, a measurement of the muon energy adds an important piece of information concerning the neutrino system. We show in this paper how the MACRO limited streamer tube system can be operated in drift mode by using the TDC's included in the QTPs, an electronics designed for magnetic monopole search. An improvement of the space resolution is obtained, through an analysis of the multiple scattering of muon tracks as they pass through our detector. This information can be used further to obtain an estimate of the energy of muons crossing the detector. Here we present the results of two dedicated tests, performed at CERN PS-T9 and SPS-X7 beam lines, to provide a full check of the electronics and to exploit the feasibility of such a multiple scattering analysis. We show that by using a neural network approach, we are able to reconstruct the muon energy for $E_\mu<$40 GeV. The test beam data provide an absolute energy calibration, which allows us to apply this method to MACRO data.Comment: 25 pages, 11 figures, Submitted to Nucl. Instr. & Meth.

Low energy atmospheric muon neutrinos in MACRO

We present the measurement of two event samples induced by atmospheric $\nu_\mu$ of average energy $\bar {E}_\nu \sim 4 GeV$. In the first sample, the neutrino interacts inside the MACRO detector producing an upward-going muon leaving the apparatus. The ratio of the number of observed to expected events is $0.57 \pm0.05_{stat} \pm0.06_{syst} \pm0.14_{theor}$ with an angular distribution similar to that expected from the Bartol atmospheric neutrino flux. The second is a mixed sample of internally produced downward-going muons and externally produced upward-going muons stopping inside the detector. These two subsamples are selected by topological criteria; the lack of timing information makes it impossible to distinguish stopping from downgoing muons. The ratio of the number of observed to expected events is $0.71 \pm 0.05_{stat} \pm0.07_{syst} \pm0.18_{theor}$ . Using the ratio of the two subsamples (for which most theoretical uncertainties cancel) we can test the pathlength dependence of the oscillation hypothesis. The probability of agreement with the no-oscillation hypothesis is 5% . The deviations of our observations from the expectations has a preferred interpretation in terms of $\nu_\mu$ oscillations with maximal mixing and $\Delta m^2 \sim 10^{-3} \div 10^{-2} eV^2$. These parameters are in agreement with our results from upward throughgoing muons, induced by $\nu_\mu$ of much higher energies.Comment: 7 pages, 6 figures. Submitted to Phys. Lett.

Measurement of the atmospheric neutrino-induced upgoing muon flux using MACRO

We present a measurement of the flux of neutrino-induced upgoing muons (~100 GeV) using the MACRO detector. The ratio of the number of observed to expected events integrated over all zenith angles is 0.74 +/- 0.036 (stat) +/- 0.046(systematic) +/- 0.13 (theoretical). The observed zenith distribution for -1.0 < cos(theta) < -0.1 does not fit well with the no oscillation expectation, giving a maximum probability for chi^2 of 0.1%. The acceptance of the detector has been extensively studied using downgoing muons, independent analyses and Monte-Carlo simulations. The other systematic uncertainties cannot be the source of the discrepancies between the data and expectations. We have investigated whether the observed number of events and the shape of the zenith distribution can be explained by a neutrino oscillation hypothesis. Fitting either the flux or zenith distribution independently yields mixing parameters of sin^2 (2theta)=1.0 and delta m^2 of a few times 10^-3 eV^2. However, the observed zenith distribution does not fit well with any expectations giving a maximum probability for chi^2 of 5% for the best oscillation hypothesis, and the combined probability for the shape and number of events is 17%. We conclude that these data favor a neutrino oscillation hypothesis, but with unexplained structure in the zenith distribution not easily explained by either the statistics or systematics of the experiment.Comment: 7 pages (two-column) with 4 figure

Search for diffuse neutrino flux from astrophysical sources with MACRO

Many galactic and extragalactic astrophysical sources are currently considered promising candidates as high energy neutrino emitters. Astrophysical neutrinos can be detected as upward-going muons produced in charged-current interactions with the medium surrounding the detector. The expected neutrino fluxes from various models start to dominate on the atmospheric neutrino background at neutrino energies above some tens of TeV. We present the results of a search for an excess of high energy upward-going muons among the sample of data collected by MACRO during ~5.8 years of effective running time. No significant evidence for this signal was found. As a consequence, an upper limit on the flux of upward-going muons from high-energy neutrinos was set at the level of 1.7 10^(-14) cm^(-2) s^(-1) sr^(-1). The corresponding upper limit for the diffuse neutrino flux was evaluated assuming a neutrino power law spectrum. Our result was compared with theoretical predictions and upper limits from other experiments.Comment: 19 pages, 8 figures, 2 table

The primary cosmic ray composition between 10**15 and 10**16 eV from Extensive Air Showers electromagnetic and TeV muon data

The cosmic ray primary composition in the energy range between 10**15 and 10**16 eV, i.e., around the "knee" of the primary spectrum, has been studied through the combined measurements of the EAS-TOP air shower array (2005 m a.s.l., 10**5 m**2 collecting area) and the MACRO underground detector (963 m a.s.l., 3100 m w.e. of minimum rock overburden, 920 m**2 effective area) at the National Gran Sasso Laboratories. The used observables are the air shower size (Ne) measured by EAS-TOP and the muon number (Nmu) recorded by MACRO. The two detectors are separated on average by 1200 m of rock, and located at a respective zenith angle of about 30 degrees. The energy threshold at the surface for muons reaching the MACRO depth is approximately 1.3 TeV. Such muons are produced in the early stages of the shower development and in a kinematic region quite different from the one relevant for the usual Nmu-Ne studies. The measurement leads to a primary composition becoming heavier at the knee of the primary spectrum, the knee itself resulting from the steepening of the spectrum of a primary light component (p, He). The result confirms the ones reported from the observation of the low energy muons at the surface (typically in the GeV energy range), showing that the conclusions do not depend on the production region kinematics. Thus, the hadronic interaction model used (CORSIKA/QGSJET) provides consistent composition results from data related to secondaries produced in a rapidity region exceeding the central one. Such an evolution of the composition in the knee region supports the "standard" galactic acceleration/propagation models that imply rigidity dependent breaks of the different components, and therefore breaks occurring at lower energies in the spectra of the light nuclei.Comment: Submitted to Astroparticle Physic

Spin asymmetry A_1^d and the spin-dependent structure function g_1^d of the deuteron at low values of x and Q^2

We present a precise measurement of the deuteron longitudinal spin asymmetry A_1^d and of the deuteron spin-dependent structure function g_1^d at Q^2 < 1 GeV^2 and 4*10^-5 < x < 2.5*10^-2 based on the data collected by the COMPASS experiment at CERN during the years 2002 and 2003. The statistical precision is tenfold better than that of the previous measurement in this region. The measured A_1^d and g_1^d are found to be consistent with zero in the whole range of x.Comment: 17 pages, 10 figure

Gluon polarization in the nucleon from quasi-real photoproduction of high-pT hadron pairs

We present a determination of the gluon polarization Delta G/G in the nucleon, based on the helicity asymmetry of quasi-real photoproduction events, Q^2<1(GeV/c)^2, with a pair of large transverse-momentum hadrons in the final state. The data were obtained by the COMPASS experiment at CERN using a 160 GeV polarized muon beam scattered on a polarized 6-LiD target. The helicity asymmetry for the selected events is = 0.002 +- 0.019(stat.) +- 0.003(syst.). From this value, we obtain in a leading-order QCD analysis Delta G/G=0.024 +- 0.089(stat.) +- 0.057(syst.) at x_g = 0.095 and mu^2 =~ 3 (GeV}/c)^2.Comment: 10 pages, 3 figure

Associations of body size with all-cause and cause-specific mortality in healthy older adults

In the general population, body mass index (BMI) and waist circumference are recognized risk factors for several chronic diseases and all-cause mortality. However, whether these associations are the same for older adults is less clear. The association of baseline BMI and waist circumference with all-cause and cause-specific mortality was investigated in 18,209 Australian and US participants (mean age: 75.1 ± 4.5 years) from the ASPirin in Reducing Events in the Elderly (ASPREE) study, followed up for a median of 6.9 years (IQR: 5.7, 8.0). There were substantially different relationships observed in men and women. In men, the lowest risk of all-cause and cardiovascular mortality was observed with a BMI in the range 25.0–29.9 kg/m2 [HR25-29.9 vs 21–24.9 kg/m2: 0.85; 95% CI, 0.73–1.00] while the highest risk was in those who were underweight [HRBMI <21 kg/m2 vs BMI 21–24.9 kg/m2: 1.82; 95% CI 1.30–2.55], leading to a clear U-shaped relationship. In women, all-cause mortality was highest in those with the lowest BMI leading to a J-shaped relationship (HRBMI <21 kg/m2 vs BMI 21–24.9 kg/m2: 1.64; 95% CI 1.26–2.14). Waist circumference showed a weaker relationship with all-cause mortality in both men and women. There was little evidence of a relationship between either index of body size and subsequent cancer mortality in men or women, while non-cardiovascular non-cancer mortality was higher in underweight participants. For older men, being overweight was found to be associated with a lower risk of all-cause mortality, while among both men and women, a BMI in the underweight category was associated with a higher risk. Waist circumference alone had little association with all-cause or cause-specific mortality risk. Trial registration ASPREE https://ClinicalTrials.gov number NCT01038583.Prudence R. Carr, Katherine L. Webb, Johannes T. Neumann, Le T. P. Thao, Lawrence J. Beilin, Michael E. Ernst, Bernadette Fitzgibbon, Danijela Gasevic, Mark R. Nelson, Anne B. Newman, Suzanne G. Orchard, Alice Owen, Christopher M. Reid, Nigel P. Stocks, Andrew M. Tonkin, Robyn L. Woods, John J. McNei

The Deuteron Spin-dependent Structure Function g1d and its First Moment

We present a measurement of the deuteron spin-dependent structure function g1d based on the data collected by the COMPASS experiment at CERN during the years 2002-2004. The data provide an accurate evaluation for Gamma_1^d, the first moment of g1d(x), and for the matrix element of the singlet axial current, a0. The results of QCD fits in the next to leading order (NLO) on all g1 deep inelastic scattering data are also presented. They provide two solutions with the gluon spin distribution function Delta G positive or negative, which describe the data equally well. In both cases, at Q^2 = 3 (GeV/c)^2 the first moment of Delta G is found to be of the order of 0.2 - 0.3 in absolute value.Comment: fits redone using MRST2004 instead of MRSV1998 for G(x), correlation matrix adde