next up previous contents
Next: Electron telescope measurements Up: Amplitude ratio method Previous: Amplitude ratio method

Electron scintillator measurements

Shown in Fig. 6.5 is an example of the electron time spectrum and fit with two exponential functions with a constant background. Plotted with error bars are the first electrons detected by electron pair scintillators for a target of 1000 T$\cdot l$ pure hydrogen, while the solid curve is a fit in the time interval of [0.02, 6] $\mu $s.


  
Figure 6.5: An example of electron time spectrum and fit with two exponential functions with a constant background. Plotted with error bars are the first electrons detected by electron pair scintillators, while the solid curve is a fit in the time interval of [0.02, 6] $\mu $s.

The resulting fit amplitudes (A1/Q1, A2, and bkgd), normalized to GMU, and the lifetimes ( and ) are given in Table 6.6, for the fit from each detector (Ege, En1, and En2) as well as that from the first hit in the three detectors (1st). Also shown is the reduced muon stopping fraction in hydrogen ( ) derived from A2/(A1/Q1 + A2) with the Huff factor for Au of Q1=0.85. The data of runs 1650 and 1654, taken with the same conditions, were analyzed both separately and together.


 
Table 6.6: The result of a fit to the electron time spectrum with two exponential functions with a constant background, and the muon stopping fraction to hydrogen ( ) derived from the amplitude ratio, A2/(A1/Q1 + A2) where Q1=0.85 is the Huff factor for Au. The fit amplitudes (A1/Q1, A2, bkgd) are normalized to GMU. Fits to the spectrum from each detector (Ege, En1, En2) as well as that for the first hit in the three detectors (1st) are listed.
Run Det. A1/Q1 A2 bkgd $\chi ^{2}$/dof(cl)  
  (10-2) (10-2) (10-4) (ns) ($\mu $s) (%)    

1650

Ege 5.06(13) 1.760(9) 2.4(4) 76.9(20) 2.030(23) 25.81(50) 0.99 (71%)  
  En1 4.86(13) 1.734(9) 2.3(4) 76.0(21) 2.027(24) 26.28(52) 1.02 (18%)  
  En2 4.86(12) 1.563(9) 2.9(4) 79.0(21) 2.041(26) 24.34(48) 1.05 (3.8%)  
  1st 14.71(22) 5.061(17) 6.9(7) 77.1(12) 2.029(15) 25.59(30) 1.09 (.03%)  

1654

Ege 4.72(12) 1.770(10) 2.5(4) 80.7(22) 2.003(23) 27.28(51) 1.01 (32%)  
  En1 4.71(12) 1.722(9) 2.1(4) 77.8(21) 2.043(24) 26.75(51) 0.99 (69%)  
  En2 4.95(13) 1.584(9) 3.0(4) 74.9(20) 2.012(24) 24.26(48) 1.04 (7.0%)  
  1st 14.30(21) 5.081(16) 6.8(7) 77.6(12) 2.017(13) 26.22(29) 1.02 (26%)  

1650

Ege 4.88(9) 1.765(7) 2.4(3) 78.8(15) 2.016(16) 26.56(36) 0.99 (62%)  
+ En1 4.77(9) 1.728(7) 2.2(3) 76.9(15) 2.035(17) 26.52(37) 1.04 (8.6%)  
1654 En2 4.90(9) 1.574(7) 2.9(3) 76.9(15) 2.025(19) 24.30(36) 1.13 (.01%)  
  1st 14.50(16) 5.071(12) 6.8(5) 77.4(9) 2.023(10) 25.91(21) 1.09 (.04%)  
 

Fits are mostly of satisfactory quality, but in some cases we observed rather low confidence levels. This could be in part due to beam related background which is correlated to the 23 MHz (43 ns) cyclotron RF cycle. This periodic background, when fitted with a straight line, would increase the total $\chi ^{2}$, but should not affect the fit parameters for the signals, since it is the average over the relevant time scale that affects the fit results[*].

Our fitted short lifetime is reasonably consistent with 74.31.5 ns for gold given by Suzuki, Measday and Roalsvig [215], but the long lifetime is somewhat shorter than 2.195 $\mu $s for hydrogen [215].

In the single hit mode[*], if the detection efficiency is high, there could be a distortion in the recorded accidental background, since the detector is more likely active at earlier time. But the individual detector fits and the fit to the first hit among them (``1st'') are consistent with each other, indicating that the effect of accidental background in distorting the fit is negligible (``1st'' can be considered as one single-hit detector with an efficiency three times as large as a single detector, hence if an efficiency dependent effect was important in the fit, it would show up as the difference between the 1st and other detectors). In fact, fits with explicitly non-constant backgrounds were tried, but the stopping fraction was found rather insensitive to the background model.

The dependence of the derived stopping fraction was tested by changing the fit region from [0.02;6] $\mu $s to [0.02; 9.5] $\mu $s. While the H2 component lifetime was slightly increased ($\sim 2$%), the variation in stopping fraction was negligible.

We have also fitted the later time t>1 $\mu $s, at which the gold signal is negligible, to a single exponential and a background. In these fits, in which 100 ns bin size was used (averaging out the RF structure), a lifetime of 2.11(1) $\mu $s was obtained with a confidence level of 8%, a lifetime closer to, yet still significantly smaller than, the literature value for muons in hydrogen. An exponential background was also tried, giving $\mu $s with a background lifetime of 0.0829 $\mu $s with a confidence level of 11%[*]. The amplitude for the hydrogen component thus derived is somewhat smaller than the ones in Table 6.6, but by no more than 3%.

Alternatively, fits were tried with the lifetime fixed to 2.195 $\mu $s and with a constant background (varying the amplitude) or an exponential background (varying the lifetime and the amplitude). Neither background gave a satisfactory fit, suggesting that the deviation of the lifetime is unlikely due to a trivial error in modeling the background, thus pointing to the existence of one or more muon loss mechanisms. This could include a small amount of non-hydrogenic contamination in the target, which would not affect the extraction of the stopping fraction, as long as it is well approximated with an exponential function. Emission of muonic protium from the solid layer, recently observed by our collaboration for the first time [216], could be partly responsible for the discrepancy also. According to Wozniak [217], several percent of muonic protium atoms may be emitted back to the gold foil, whereby the muon transfers and is captured by the gold nucleus.


 
Table 6.7: Sensitivity of the reduced stopping fraction to the shift of time zero ( ) for the first electron time spectrum for the sum of runs 1650 and 1654. is the relative shift with respect to the value. The variation in is smaller than the last digit shown.
(ns) (%) cl (%)
-4 27.07 (23) 4.4% 1.09 0.053
-2 26.45 (22) 2.1% 1.09 0.045
-0 25.91 (21) 0 % 1.09 0.043
+2 25.38 (21) -2.0% 1.09 0.041
+4 24.85 (20) -4.1% 1.09 0.039
 

The sensitivity of to the shift of the time zero definition ( ) was also investigated. An example for the first electron spectrum is given in Table 6.7. The fit to the individual detectors as well as to individual runs gives similar trends. The error 4 ns is a rather conservative estimate of the shift, and it is more likely to be less than 2 ns. The variation in the amplitude for the hydrogen component in this time scale is, nevertheless, completely negligible, hence the change in comes from that in the gold component. In any case, this alone cannot explain the 20% difference between method 1 and method 2, as will be discussed later.

While for Ege and En1 are consistent with each other, there is a systematic difference between En2 and the others, with En2 being relatively lower by nearly 10%. This discrepancy will be addressed later in relation to the acceptance difference for the upstream and downstream targets.

We note that agreement between Ege and En1, which were symmetrically located across the target (beam-left and beam-right, respectively), exclude the possible effects of polarized muon precession about the vertical axis.


next up previous contents
Next: Electron telescope measurements Up: Amplitude ratio method Previous: Amplitude ratio method