Single Spin Asymmetry Analysis » History » Version 2

« Previous - Version 2/5 (diff) - Next » - Current version
Stephen Wood, 11/14/2018 01:53 PM
How to modify replay script

Single Spin Asymmetry Analysis


During the fall 2018 running of the Kaon-LT experiment, polarized beam was delivered to Hall B. According to Dave Gaskell, at the energy we ran at, the spin precession in the Hall C arc was such that we got nearly as much polarization as Hall B, so we had a polarization of about 85%. This means that, for hadron production out of the plane defined by the beam and scattered electron, we should see an asymmetry in the yield between the two beam helicity states. This asymmetry is known variously as the single spin asymmetry, Sigma_LT', or the fifth structure function.


Hall B paper exclusive pion production single spin asymmetries, Park et. al. 2008

How to add helicity to ROOT files.

The helicity of the beam potentially reversed every 1/30 of a second according to a pseudorandom sequence. The helicity signals delivered to the hall are delayed by 8*(1/3) second. This delayed helicity signal is recorded in several ADCs which are read each event. In order to extract the helicity and undo the delayed reporting of helicity, a recent version of hcana which includes the class THcHelicity is needed. In order to add the true beam helicity to each event in the root tree, make sure that the following are implemented in the replay script and analysis def file.

After the trigger apparatus is defined, add the following lines to the analysis script:

// Add helicity detector to grigger apparatus
THcHelicity* helicity = new THcHelicity("helicity","Helicity Detector");

In the def file (coin_production.def) add these lines:

block T.helicity.hel
block T.helicity.helrep
block T.helicity.mps

T.helicity.hel will be +/-1 indicating the actual helicity of the beam, event-by-event. If it is zero, then either the event was during the MPS (settling time between helicity reversals.) or the analyzer was unable to decode the helicity properly. helrep is the reported T.helicity.helrep which is delayed. T.helicity.mps is 1 if the event was during the MPS.


  • Get the history of the status of the Insertable Half Wave Plate (EPICS signal IGL1I00DI24_24M). The IHWP is a passive optical element on the injector laser table that reverses the helicity of the electron beam, but does not reverse the reported helicity. Experimental asymmetries will be reversed when the IHWP is in, so the analysis must reverse the asymmetries for runs when this IHWP was in. It is possible that the state of the IHWP was changed in the middle of one of runs. Any such runs need to be identified and analyzed specially.
  • Get the history of Hall B beam polarization measurements. Translate these polarizations into Hall C polarizations.
  • FIgure out how to bin and display results. Should asymmetries be presented as a function of t or Theta_hq (cm or lab). How finely can we bin in t or theta?
  • For pion asymmetries. Determine accidental background. This asymmetry from accidental backgrounds should be zero, but this background will reduce or dilute the observed asymmetry.
  • For Kaon asymmetries it is more complicated. As well as a dilution from accidentals, the radiative tail from exclusive pions misidentified as Kaons contributes to the observed Kaon asymmetries. The asymmetry from these pions will not be zero. In order to correct for this to get Kaon asymmetries, the accepted yield of the tail and its asymmetry needs to be determined.
  • Determine if it is possible to get any meaningful asymmetries from the baryons beyond the Lambda.

a_pion.png - Beam asymmetry of pions vs. phi (10.6 KB) Stephen Wood, 11/14/2018 04:11 PM

a_kaon.png - Beam asymmetry of "kaons" vs. phi (9.13 KB) Stephen Wood, 11/14/2018 04:12 PM