Hi Tim,
Look in 'lrn_event.cc', line 1003 for the deuterium case.
H2 gas only has vector states: 1,2,3,4, and we only injected the
pure states 1,3 in our experiment (with no mixing at low B field).
See my thesis p. 71 for a description of the H and D states. So the
target + state charge is in inhib_charge_p(1,1) and - charge in
inhib_charge_p(3,3). I can't remember why unpolarized got stuck in
(1,2), maybe for legacy reasons?
For D2, three mixtures of states were injected: states 1 and 6,
both with m_I=1, states 3 and 4, both with m_I=-1, and states 2 and
5, with m_I=0, but p_zz=-2. These are shown in the DST {inhib,scalr}
_charge_{p,m} matrices below. Note there is a separate matrix for
each beam helicity. In contrast, the init.C matrix 'charge[0]' is a
2x5 matrix with the row for electron helicity and the column for
target state.
So your table is correct for H2 if you set P_zz=0 everywhere and
P_z=0 is actually unpolarized (buffer gas system). You should add an
extra row with the DST matrix elements below for D2.
--hope this helps, Chris
lrn_event.cc:995
if (targ_gas == 1) { // H2
scpp = scalr_charge_p(1,1); bgpp = inhib_charge_p(1,1);
scpm = scalr_charge_p(3,3); bgpm = inhib_charge_p(3,3);
scp0 = scalr_charge_p(1,2); bgp0 = inhib_charge_p(1,2);
scmp = scalr_charge_m(1,1); bgmp = inhib_charge_m(1,1);
scmm = scalr_charge_m(3,3); bgmm = inhib_charge_m(3,3);
scm0 = scalr_charge_m(1,2); bgm0 = inhib_charge_m(1,2);
}
else if (targ_gas == 2) { // D2
scpp = scalr_charge_p(1,6); bgpp = inhib_charge_p(1,6);
scpm = scalr_charge_p(3,4); bgpm = inhib_charge_p(3,4);
scp0 = scalr_charge_p(2,5); bgp0 = inhib_charge_p(2,5);
scmp = scalr_charge_m(1,6); bgmp = inhib_charge_m(1,6);
scmm = scalr_charge_m(3,4); bgmm = inhib_charge_m(3,4);
scm0 = scalr_charge_m(2,5); bgm0 = inhib_charge_m(2,5);
scp0p = scalr_charge_p(3,6); bgp0p = inhib_charge_p(3,6);
scm0p = scalr_charge_m(3,6); bgm0p = inhib_charge_m(3,6);
}
http://blast.lns.mit.edu/PUBLIC_RESULTS/THESES/Crawford_thesis_revA.pdf
_______________________________________
TA-53/MPF-1/D111 P-23 MS H803
LANL, Los Alamos, NM 87545
505-665-9804(o) 665-4121(f) 662-0639(h)
_______________________________________
On Feb 9, 2006, at 10:23:36, Timothy Paul Smith wrote:
> Hello People,
>
> I am trying to understand charge matrix in the DST files.
>
>
> From Chi:
> Chi>> Aaron once has a short documentation about the design I
> remember.
>
> I can not find this - any ideas anyone?
>
> From Chris:
> Chris>> and m_ij is the charge on combined states 'i,j',
> Chris>> for the 6 deuterium states.
>
> Is this the |1> = F=3/2, m+f=+3/2, m_1=+1, m_s=+1/2
> |2> ...
>
> states?
>
>
>
>
>
> I spent some time with Vataliy yesterday and he gave me the
> following translation for the LR files.
>
> (Vitaliy's thesis)
> (eq 5.3) (eq 5.4)
>
> n(h, P_z, P_zz) lr matrix(#,#)
> n( 1, 1, 1) matrix(0,0)
> n( 1,-1, 1) matrix(0,1)
> n( 1, 0,-2) matrix(0,2)
> n(-1, 1, 1) matrix(1,0)
> n(-1,-1, 1) matrix(1,1)
> n(-1, 0,-2) matrix(1,2)
>
>
> I am trying to understand these interms of the matrix
> fScaler->inhib_charge_p(#i,#j)
>
> here the _p -> h=+1 and _n -> h=-1
>
> inhib_charge_p(j,k) = recon.fEpics->inhib_charge[0][j][k];
> inhib_charge[(int)fScaler->fScalerHel][s1][s2] += inhbq;
>
> it looks to me that s1 and s2 are spin bits from the ABS?
>
>
> In the code "lrn_event" I see
>
> bgpp = fill->inhib_charge_p(1,1);
> bgpm = fill->inhib_charge_p(3,3);
> bgp0 = fill->inhib_charge_p(1,2);
> bgmp = fill->inhib_charge_m(1,1);
> bgmm = fill->inhib_charge_m(3,3);
> bgm0 = fill->inhib_charge_m(1,2);
>
>
> from init.C
>
> ch_tmp(0,0) += bgpp; // +/+ +/+/+
> ch_tmp(0,1) += bgpm; // +/- +/-/+
> ch_tmp(1,0) += bgmp; // -/+ -/+/+
> ch_tmp(1,1) += bgmm; // -/- -/-/+
> ch_tmp(0,2) += bgp0; // +/u +/0/-
> ch_tmp(1,2) += bgm0; // -/u -/0/-
>
> so can I do this?:
>
> inhib_charge_p(0,0) = total charge with beam h=+1
>
> n(h, P_z, P_zz) lr matrix(#,#) DST
> n( 1, 1, 1) matrix(0,0) inhib_charge_p(1,1) n( 1,-1,
> 1) matrix(0,1) inhib_charge_p(3,3)
> n( 1, 0,-2) matrix(0,2) inhib_charge_p(1,2)
> n(-1, 1, 1) matrix(1,0) inhib_charge_m(1,1)
> n(-1,-1, 1) matrix(1,1) inhib_charge_m(3,3)
> n(-1, 0,-2) matrix(1,2) inhib_charge_m(1,2)
>
> --
>
> _____________________________________________________________
> Timothy Paul Smith Assistant Research Professor
> Dartmouth College, Physics (603)646-9346
> 6127 Wilder, Hanover, NH 03755 timothy.p.smith@dartmouth.edu
> http://he3.dartmouth.edu/TimSmith
>
This archive was generated by hypermail 2.1.2 : Mon Feb 24 2014 - 14:07:33 EST