Today's Progress 29. Apr. 2010

Final evaluation of per-stopped-K- formation rate of 4&LambdaHe -> d + d decay

Here,we try to finalize the per-stopped-K- formation rate of 4&LambdaHe -> d + d decay. The preliminary results are given in 2007 Dec 1st report, but now we consider the exact E549/E570 geometry, and the energy loss correction efficiency as well.

Raw dd coincidence spectra

Since deuteron selection procedure has been updated ( 2007 Dec 19th report) after the latest report of dd coincidence spectra ( 2007 Dec 1st report), we present here the final dd-coincidence spectra. After that, we try to give an exact normalization.
dd momenta of both arms.
Correlation of dd momenta
Opening angle of dd 3-momenta
dd invariant mass (red-shaded:cos(dd) < -0.990). Shaded 57 events are identified as dd decay events.
Deuteron 3-momentum VS dd opening angle
Total 3-momentum of detected dd pair
Invariant mass VS total 3-momentum of dd
Invariant mass VS opening angle

Acceptance study for the hypernuclear cascade processes

Now, we are interested in to estimate exactly the acceptance of a cascade of hypernuclear formation,

(4He-K-)atomic ->4&LambdaHe+&pi- (255.4 MeV/c) ,

and its non-mesonic two-body weak decay,

4&LambdaHe -> d+d (1 : 571.8 MeV/c),

4&LambdaHe -> t+p (2 : 508.4 MeV/c),

4&LambdaHe -> 3He+n (3 : 507.4 MeV/c).

The particle and nulcear mass table for the calculation of decay momentum is given below. It should be noted that the nuclear masses defined in GEANT3 is sometimes different from the exact values taken from here .

Particle/Nuclear Mass table

Particle/Nucleus mass (MeV/c2)
4He 3727.40841
3He 2808.41374
t 2808.94327
d 1875.62764
p 938.27203
n 939.56536
&Lambda 1115.683
K- 493.677
&pi- 139.57018
4&LambdaHe 3921.707 (B.E. measured from 3He + &Lambda at rest is 2.39 MeV, taken from M. Julic et al., Nucl. Phys B52 (1973) 1.)

The normalization of dd mode

Since the dd decay mode was observed, the exact normalization is performed by

, where

  • P- : Free-decay rate of atomic K- by meta-stable states, 3.5 +- 0.5 %.
  • Br(K- -> μ&nu) : Branching ratio to Kμ2 decay, 63.55 %.
  • &epsilon-&mu : Physical acceptance of &mu-, 7.55 %.
  • &epsilon-delay : Efficiency by delayed timing gate 1.2 nsec, 89.1%.
  • &Sigmacycles C- : Total number of observed dd decay events on the dd invariant mass spectrum: 57
  • N-Kμ2 : cycle-by-cycle value of detected &mu- number under the condition, Tdelay > 1.2 nsec. 111561/111704/44124 for E549/E570-1/E570-2.
  • &epsilon-dd : cycle-by-cycle value of physical acceptance of dd, to be finalized here.
  • In order to estimate physical acceptance &epsilon-dd defined as,

    Ndetected = Nproduced*&epsilon-dd,

    we perform a Monte-Carlo simulation of the hyper-nuclear formation and its dd decay, and conditions of the simulation are as shown below. The physical acceptance includes
  • Geometrical acceptance by PA-PDC-PB-NT (double),
  • In-flight reaction (double),
  • Analysis efficiency of deuteron selection by 1/&beta VS total energy at PB + NT (double),
  • Analysis efficiency of energy loss correction (double),
  • PA-PB timing gate, 45 nsec by E549/E570 electronics (double),
  • and
  • Efficiency of the opening-angle selection, cos(dd)<-0.990.
  • Since we count the dd pair on the invariant mass spectrum, the uncertainty of the 4&LambdaHe lifetime only very slightly affects the opening-angle selection.

    The upper limit evaluation of tp mode

    On the other hand, tp decay mode was not observed. Then, the upper limit with a certain confidence lecel is evaluated by

    , where
  • C-U.L. : Upper limit of Gaussian peak area adopting corresponding confidence level on the t-coincident proton momentum spectrum at 508 MeV/c (for the value, see a separate report),
  • &epsilon-tp : cycle-by-cycle value of physical acceptance of tp, to be finalized here.
  • In order to estimate physical acceptance &epsilon-tp defined as,

    Ndetected = Nproduced*&epsilon-tp,

    we perform a Monte-Carlo simulation of hyper-nuclear formation and its tp decay, and conditions of the simulation are as shown below. The physical acceptance includes
  • Geometrical acceptance by PA-PDC-PB-NT for p,
  • In-flight reaction of p,
  • Analysis efficiency of proton selection by 1/&beta VS total energy at PB + NT,
  • Analysis efficiency of energy loss correction for p,
  • PA-PB timing gate, 45 nsec for p, by E549/E570 electronics,
  • and
  • Fraction of hypernuclear decay at rest,
  • Efficiency of the opening-angle selection, cos(tp)<-0.990.
  • It should be noted that the uncertainty of the life time now largely affects the acceptance, since it involves a non-negligible change of fraction of at-rest decay.

    Conditions of the Monte-Carlo.

    software GEANT3.21(partice tracking tuned for slow particles)
    Formation reaction (K--4He)atomic -> 4&LambdaHe + &pi- (254.6 MeV/c with 4&LambdaHe mass shown below), 100%
    4&LambdaHe mass M3He(2809.23)+M&Lambda(1115.684)-2.39 = 3922.5 (MeV/c2)
    4&LambdaHe decay modes/bracnhing ratio on the simulation d+d:t+p:3He+n=40%(573.2MeV/c):40%(509.1MeV/c):20%(507.4MeV/c)
    4&LambdaHe Life time 257 psec : figure / 250 +- 18 psec : acceptance
    (K--4He)atomic life time 6.6*10-17 sec (Γ=10 eV)
    Generated number of (K--4He)atomic / 4&LambdaHe 2*107
    Target center (-0.3,0.,1.3):E549 (-0.5,0.5,1.2):E570
    x/y generation point distribution 4.0 cm &sigma Gaussian centered at (x,y)=(-0.3,0.)
    z generation point distribution uniform
    Multiple scattering on(Moliere)
    Energy loss straggling on(Gauss/Landau/Vavilov are internally selected adequately)
    In-flight nuclear reaction of p/d on(GHEISHA, unique selection for in-flight reaction loss for d)
    Coincidense time gate for PA-PB 45 nsec
    Birk's coefficients0.013/(MeV/cm), 9.6E-06/(MeV/cm)2 for PA/PB/NT
    Time resolution of PA/PB 60/90 psec (~108.2 psec in total), fixed at MIP-measured value.
    Proton/Deuteron selection Simulated with 1/&beta VS total Light output (with infinite energy resolution).
    PDC hit Only 1 particle is accepted per chamber. Multi track events (p+&pi-,etc) are eliminated.
    Inefficiency of energy loss correction Newly considered. As the reaction vertex, exact position is adopted. The direction vector at the PDC position is adopted as the direction vector to simulate properly the experimental effect.

    Properties of the process

    Momentum distribution of 4&LambdaHe at the decay. ~86.4% decay at rest.
    TOF of 4&LambdaHe until stopping if they decay at rest.
    Range spectrum of 4&LambdaHe. Its tipical range is ~0.5 mm in 4He.
    Decay timing spectrum. The decay timing is defined as the time difference between the formation of K--4He atom and the decay of 4&LambdaHe regardless the status (in-flight or at rest) at the decay. The lifetime is found to be 257 psec by the fitting of an exponential function.
    Range spectra for 3H and 3He in super-fluid 4He. Tritons can fire the PA, but 3He cannot.

    Results

    As already shown, d+d decay channel was observed by back-to-back coincidence of dd pair. As expected from the range spectrum, p+t channel would be barely seen with the triton stops on PA leaving a charged particle track on PDC, in coincidence with a proton with 508 MeV/c momentum detected on the opposite arm. Therefore, we try to evaluate the physical acceptance of ....

    1. dd events - back-to-back dd coincidence on PDC/PA/PB,

    2. pt events - p on PDC/PA/PB and t on PDC/PA on the opossite arm.

    In these conditions, simuated decay vertex and cos(dd/pt) distributions with the geometrical bias are exhibited below. The red spectra are from stopped 4&LambdaHe decay, while the blue are from in-flight events.
    Simulated decay vertex and cos(dd) distribution for dd events. Event selection procedure of deuteron is applied on both arm as done for data.
    Simulated decay vertex and cos(pt) distribution for pt events. proton slection procedure is applied as done for data, and t on PA is required on the opposite arm.
    Simulated deuteron momentum distribution by 4&LambdaHe->dd events on the left arm. The top is constructed from the event set if a deuteron is detected and identified with successfull energy loss correction on the left arm. In the bottom, we further require detection and identification of another d in the counter-PDC-PA-PB R arm, which directly corresponds to the measured momentum spectrum presented above. The blue shows the whole events, while the red shows the events from stopped 4&LambdaHe decay. The Gaussian fitting to the bottom-red gives ~2.0 MeV/c as the standard deviation. Linear-scale and log-scale figures are presented simultaneously.
    Simulated proton momentum distribution by 4&LambdaHe->tp events. The top is constructed from the event set if a proton is detected and identified with successfull energy loss correction. In the bottom, we further require detection of t in the counter-PDC-PA arm. The blue shows the whole events, while the red shows the events from stopped 4&LambdaHe decay. The Gaussian fitting to the bottom-red gives ~4.8 MeV/c as the standard deviation.

    The step-by-step evaluation of &epsilondd and &epsilontp are tabulatted below:

    Generated and Accepted Event Numbers for dd mode

    Starting from the column 4, we add additional selection as described on the top row. The numbers rounded by brackets are from decay at rest. &epsilondd is just defined by the bold number in the 8th colmun divided by the number in the 3rd colmun.
    Cycle Life time (psec) generated event No. of dd decayPDC-PA-PB hit of d PDC-PA-PB hit of another d time gate (45 nsec) + d selection on analysis (double) Successful energy loss correction (double) cos(dd)<-0.99 &epsilondd (%)
    E5492327.63051E+06650304(550444)223544(205111)203434(187051)199988(184042)187403(183630)2.456
    E5492507.63429E+06649406(555910)225301(208034)205168(189755)201741(186720)189958(186321)2.488
    E5492577.63669E+06650645(559749)226593(209722)206411(191409)202952(188355)191319(187908)2.505
    E5492687.63644E+06649973(562708)226281(210275)206050(191827)202674(188834)191725(188394)2.511
    E5702327.58286E+06636601(539080)209310(192650)188505(173821)184781(170521)173353(170015)2.286
    E5702507.57735E+06636998(545510)210741(195079)189808(176021)186048(172703)175394(172151)2.315
    E5702577.57948E+06637705(548961)211302(196179)190181(176847)186521(173576)176151(173030)2.324
    E5702687.58212E+06637878(552292)212166(197620)190785(177932)187117(174631)177065(174052)2.335

    Generated and Accepted Event Numbers for tp mode

    Starting from the column 4, we add additional selection as described on the top row. The numbers rounded by brackets are from decay at rest. &epsilontp is then just defined by the bold number in the 8th colmun divided by the number in the 3rd colmun. The statistics of tp coincidence events are not sufficient here, so we prepare 25 times higher statistics data (Br: 40%->99% / event number: 2×107->2×108) later.
    Cycle Life time (psec) generated event No of tp decay.PDC-PA-PB hit of p PDC-PA hit of t time gate (45 nsec) + selection on analysis (p only) Successful energy loss correction (p only) cos(tp)<-0.99 &epsilontp(%)
    E5492327.63546E+06535167(454178)30259(23249)26244(20201)26036(20045)20373(19431)2.545E-01
    E5492507.63493E+06533210(458144)29737(23419)25772(20382)25588(20233)20461(19621)2.570E-01
    E5492577.63293E+06533087(459672)29730(23434)25817(20380)25620(20228)20504(19649)2.574E-01
    E5492687.63342E+06533432(462909)29822(23819)25821(20692)25631(20538)20746(19948)2.613E-01
    E5702327.58012E+06530726(450900)21837(16140)18944(14024)18767(13899)14136(13427)1.771E-01
    E5702507.57949E+06529807(454815)22146(16667)19355(14609)19175(14470)14699(13979)1.844E-01
    E5702577.58236E+06530861(458013)21583(16400)18644(14199)18475(14074)14210(13600)1.794E-01
    E5702687.57706E+06529216(459550)21810(16708)19044(14650)18858(14510)14699(14014)1.850E-01

    Finalized per-stopped-K- branching ratio of dd mode

    Substituting all numerical values, we finally obtain the per-stopped-K- branching fraction of the hyper-nuclear cascade reaction,

    (4He-K-)atomic-> &pi-+4&LambdaHe,

    4&LambdaHe-> d+d,

    as

    Br((4He-K-)atomic-> &pi-+4&LambdaHe ; 4&LambdaHe -> d+d)=1.34±0.19(syst)±0.18(stat)×10-5 (per stopped K-).

    As the source of the systematic error, we considered only the error attributed to the free-decay rate of K- in meta-stable orbits. Relative error attributed to the quantitiy is more than one order larger than other possible sources, and hence the simple estimation is obviously justified.