<div dir="ltr"><div><br></div><div>OK, </div><div><br></div><div>I think 1/2 hour is not sufficient, so how about 2:00pm.   (I&#39;ll skip my student office hours.)</div><div><br></div><div>JP?</div><div>Dustin? </div><div>

Oscar?</div><div><br></div><div>Ellie and Patricia have already said they are free this time, and anytime today is unfortunately difficult for Narbe.</div><div><br></div><div>thanks</div><div class="gmail_extra">

<br></div><div class="gmail_extra">-Karl<br clear="all"><div><br></div>
<br><br><div class="gmail_quote">On Wed, May 1, 2013 at 10:40 AM, Oscar  Rondon-Aramayo &lt;<a href="mailto:or@cms.mail.virginia.edu">or@cms.mail.virginia.edu</a>&gt; wrote:<br><blockquote class="gmail_quote">

Hi Karl,<br>
<br>
We met with Dustin  last evening and after going in detail over the formulas<br>
for the ratio Npol/Nun, we found that the unpolarized sigma_N, sigma_D and<br>
sigma_He (see my last email) can be collected in one group, which cancels<br>
with the denominator (all unpolarized), leaving a term<br>
sigma_D*Azz*Pzz/denominator, which I realized can be written as f*Azz*Pzz, f<br>
= dilution factor.<br>
<br>
With the dilution factor, the formulas in Dustin&#39;s third row of equalities<br>
in his Observables2 report, which are valid only for pure D (the HERMES<br>
case), can also be used for ND3 targets.<br>
<br>
In summary, we can take the ratio of the pol to unpol counts, which takes us<br>
to Azz, at the price of the dilution factor and its error, plus the need to<br>
use some form of F1 to get b1 from Azz, or the difference, which takes us<br>
directly to b1.<br>
<br>
In both cases the systematic errors, other than the charge and detector<br>
efficiency are normalizations, and since the error on Pzz is expected to<br>
dominate, it really is a matter of taste, once we have the numbers on hand.<br>
We&#39;ll surely try both.<br>
<br>
For the statistical errors, f enters in the Azz time estimate, but Q*A*l*pf<br>
enter in the difference (Pzz is in both). I need to do some numbers yet<br>
(everyone should try) to compare the two approaches.<br>
<br>
Finally, today we have the SANE meeting at 3:30, so I can join b1 from 1:00<br>
to 3:30.<br>
<br>
Cheers,<br>
<br>
Oscar<br>
<div class=""><div class="h5"><br>
On Wed, 1 May 2013 09:18:48 -0400<br>
  Karl Slifer &lt;<a href="mailto:karl.slifer@unh.edu">karl.slifer@unh.edu</a>&gt; wrote:<br>
&gt; Hi all,<br>
&gt;<br>
&gt; The methodology is the central question and I think we have to resolve any<br>
&gt; lingering doubts today.  I highly encourage that everyone really read<br>
&gt; Oscar&#39;s note (Eq 19 and 20) and his last email before we discuss today.<br>
&gt;<br>
&gt; I would really not like to delay till tomorrow if possible since time is<br>
&gt;so<br>
&gt; tight. I hope we can get a majority to participate at 3pm.  Please let me<br>
&gt; know if you can&#39;t.<br>
&gt;<br>
&gt; -Karl<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; ---<br>
&gt; Karl J. Slifer<br>
&gt; Assistant Professor<br>
&gt; University of New Hampshire<br>
&gt; Telephone : <a href="tel:603-722-0695">603-722-0695</a><br>
&gt;<br>
&gt;<br>
&gt; On Tue, Apr 30, 2013 at 5:59 PM, O. A. Rondon &lt;<a href="mailto:or@virginia.edu">or@virginia.edu</a>&gt; wrote:<br>
&gt;<br>
&gt;&gt; Hi Dustin,<br>
&gt;&gt;<br>
&gt;&gt; Dustin Keller wrote:<br>
&gt;&gt; &gt; You can only benefit from the systematic reduction if you us Azz as<br>
&gt;&gt; &gt; discussed yesterday.  But at this point I am not partial.<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; dustin<br>
&gt;&gt; &gt;<br>
&gt;&gt;<br>
&gt;&gt; In the experiment, we only have counts. What we need to show to the PAC<br>
&gt;&gt; is how we go from the counts Npol and Nu, to Azz or b1. A measured<br>
&gt;&gt; quantity needs to be on one side and physics on the other. Lets say we<br>
&gt;&gt; start with your ratio Npol/Nu - 1 = Pzz*Azz, which only requires Pzz &gt;0.<br>
&gt;&gt;<br>
&gt;&gt; We need to prove that the lhs reproduces the rhs. We have, in general,<br>
&gt;&gt; N = Q*e*A*l*sigma. But since N are counts from everything in the target,<br>
&gt;&gt; it is not a simple matter of canceling quantities that stay the same<br>
&gt;&gt; when the polarization changes:<br>
&gt;&gt;<br>
&gt;&gt; Npol = Qpol*epol*Apol*lpol*sigma_pol<br>
&gt;&gt;      = Qpol*epol*Apol*lpol*[(sigma_N+3*sigma_Dpol)*pf + sigma_He*(1-pf)]<br>
&gt;&gt;<br>
&gt;&gt; Nu   = Qu*eu*Au*lu*[(sigma_N+3*sigma_D)*pf + sigma_He*(1-pf)]<br>
&gt;&gt;<br>
&gt;&gt; sigma_N and sigma_He are the same, always unpol. And<br>
&gt;&gt; sigma_Dpol = sigma_D(1+Pzz*Azz).<br>
&gt;&gt;<br>
&gt;&gt; Then, since Apol = Au = A, and lpol = lu = l,<br>
&gt;&gt;<br>
&gt;&gt; Npol/Nu =<br>
&gt;&gt; (Qpol/Qu)*(epol/eu)*[(sigma_N+3sigma_D(1+ Azz*Pzz))*pf+..)]/[(sigma_N+..<br>
&gt;&gt;<br>
&gt;&gt; where I just put ..., because I don&#39;t see how it can be simplified to<br>
&gt;&gt; just leave Azz*Pzz + 1, to equal the rhs.<br>
&gt;&gt;<br>
&gt;&gt; On the other hand, if instead of taking the ratio Npol/Nu first, we take<br>
&gt;&gt; the difference first, it&#39;s indeed possible to isolate the required<br>
&gt;&gt; Pzz*b1 on on side, like I do in my draft, eq. (19) or (20). And in<br>
&gt;&gt; fact, we don&#39;t even need to bother with Azz, because we get b1 without<br>
&gt;&gt; having to multiply Azz by F1, introducing one more systematic error.<br>
&gt;&gt;<br>
&gt;&gt; So, in summary, once one substitutes all the ingredients for your sigmas<br>
&gt;&gt; we get, or ought to get, eq.(19) or (20) back.<br>
&gt;&gt;<br>
&gt;&gt; In both of those equations, the systematics for Pzz, A, and l(pf) are<br>
&gt;&gt; normalization factors, just like we want them to be, for control of<br>
&gt;&gt; systematics, but the charge and the detector efficiency are not common<br>
&gt;&gt; factors, they depend on the period when the data are taken, either pol.<br>
&gt;&gt; or unpol.<br>
&gt;&gt;<br>
&gt;&gt; My point is that for the proposal, we must spell this all out, to give<br>
&gt;&gt; explicit sources of errors, and to calculate times or statistical errors<br>
&gt;&gt; correctly. For example, the statistical error must be sqrt(Npol + N_U) ~<br>
&gt;&gt; sqrt(2N), because it is just the error of a difference, etc.<br>
&gt;&gt;<br>
&gt;&gt; We need to have a consensus on how the method is going to be described<br>
&gt;&gt; in the proposal, which needs to be done in the most precise way to avoid<br>
&gt;&gt; any confusion due to ambiguities.<br>
&gt;&gt;<br>
&gt;&gt; Cheers,<br>
&gt;&gt;<br>
&gt;&gt; Oscar<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; &gt; On Tue, 30 Apr 2013, O. A. Rondon wrote:<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt; Hi,<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; Since I couldn&#39;t stay until the end of the meeting, and I don&#39;t think<br>
&gt;&gt; &gt;&gt; there will be minutes of it, I would like to share some ideas for the<br>
&gt;&gt; &gt;&gt; proposal draft.<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; Basically, what we need is an equation with the measured quantity on<br>
&gt;&gt;one<br>
&gt;&gt; &gt;&gt; side and b1 or Azz on the other. Based on what I think the consensus<br>
&gt;&gt; &gt;&gt; was, to measure polarized minus unpolarized counts on a single cup with<br>
&gt;&gt; &gt;&gt; the target field aligned along the beam, I&#39;ve updated the draft of my<br>
&gt;&gt; &gt;&gt; method, see subsection 0.2, which discusses this. Eq. (19) or eq. (20)<br>
&gt;&gt; &gt;&gt; meet the conditions stated above. This is the approach I would<br>
&gt;&gt;subscribe<br>
&gt;&gt; &gt;&gt; to, unless there is another version that is shown to also represent the<br>
&gt;&gt; &gt;&gt; procedure, which should be circulated as soon as possible. The draft<br>
&gt;&gt; &gt;&gt; is here<br>
&gt;&gt; &gt;&gt; <a href="http://twist.phys.virginia.edu/~or/b1/b1_method-v2.pdf">http://twist.phys.virginia.edu/~or/b1/b1_method-v2.pdf</a><br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; Cheers,<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; Oscar<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; _______________________________________________<br>
&gt;&gt; &gt;&gt; b1_ana mailing list<br>
&gt;&gt; &gt;&gt; <a href="mailto:b1_ana@jlab.org">b1_ana@jlab.org</a><br>
&gt;&gt; &gt;&gt; <a href="https://mailman.jlab.org/mailman/listinfo/b1_ana">https://mailman.jlab.org/mailman/listinfo/b1_ana</a><br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; _______________________________________________<br>
&gt;&gt; b1_ana mailing list<br>
&gt;&gt; <a href="mailto:b1_ana@jlab.org">b1_ana@jlab.org</a><br>
&gt;&gt; <a href="https://mailman.jlab.org/mailman/listinfo/b1_ana">https://mailman.jlab.org/mailman/listinfo/b1_ana</a><br>
&gt;&gt;<br>
<br>
_______________________________________________<br>
b1_ana mailing list<br>
<a href="mailto:b1_ana@jlab.org">b1_ana@jlab.org</a><br>
<a href="https://mailman.jlab.org/mailman/listinfo/b1_ana">https://mailman.jlab.org/mailman/listinfo/b1_ana</a><br>
</div></div></blockquote></div><br></div></div>