<div dir="ltr">Hello Ken and Peter,<div><br></div><div>With the loss of the beam this week, I have finally had a moment to go back and tackle my backlog of unanswered questions regarding the beamline. Back in early January, you wrote to me with a specific question about how to translate between our two angle conventions. Now, with the help of your animated gif, I am able to provide what I think is a definitive answer. First some definitions.</div><div><ul><li>I call the Glasgow secondary angle theta_ortho because it is called the orthogonal angle in the animated gif labels. Its counterpart in the Glasgow convention I will call theta_edge, which is the Glasgow angle that sets the location of the primary (2,-2,0) coherent edge. So theta_ortho, theta_edge, and phi are the three angles that set the orientation of the diamond in the Glasgow code.</li><li>The UConn code uses a theta_small for the angle that sets the location of the coherent edge, and the complementary tilt angle I call theta_large. Of course one may set theta_large < theta_small, they are just labels that indicate their intended use in the code. So theta_large, theta_small, and phi are the three angles that set the orientation of the diamond in the Glasgow code.</li></ul><div>After a detailed study, I now conclude that our two definitions are identical, that is theta_ortho == theta_large and theta_edge == theta_small in terms of their definitions. To show that this actually works, I have done a special exercise with my code, and created an animated gif that I can compare directly with yours. Recall that the UConn code is a generator, not an analytical formula, so it takes a lot of Monte Carlo to get sufficient statistics for a detailed comparison. Because of this, I created a toy integrator that takes the kernel from my generator and creates an approximate curve to describe the distribution that it would generate in the large N limit. There are some approximations that I used in this.</div></div><div><ol><li>Only lattice vectors with h,k,l in the range [-10, 10] are included. Your code generates peaks outside this range, perhaps as large as [-20, 20], but I thought that [-10,10] would be sufficient to test the equivalence of our angle definitions.</li><li>All of the edges are smeared with the same smearing function as the (2,2,0) vector would be. This simplifies the code that generates the peaks, but underestimates the broadening of the peaks at higher order. Those should be smeared by an additional factor of the order multiplier: eq. (4,0,0) edge is under-smeared by factor sqrt(2), (6,6,0) is under-smeared by factor sqrt(18), etc. But the edges should all be in the right place, just somewhat sharper than they would be in real life.</li><li>We are probably using somewhat different values for some physical constants, because our edges are not at exactly the same place for the same angle settings. But the patterns are the same, and the evolution of the peaks is unmistakably identical. I give my exact values for the physical constants below, so you can compare with what you have in your code.</li></ol><div>Here is a link to my animated gif, so you can have a look for yourself and see if you agree.</div></div><div><a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__drive.google.com_file_d_1F-2DhT7rM8ikzJVNMDmODWsCXKH391mnZF_view-3Fusp-3Dsharing&d=DwMFaQ&c=lz9TcOasaINaaC3U7FbMev2lsutwpI4--09aP8Lu18s&r=lneuPoR_Y_qxqVHZY6VsYU8W43i3GcK7nXR2epWlm0Y&m=U5TV4QD4zS6w7-iHN4pT2AmN9IqSKGEtBOY0a9P7JHA&s=wBPkKRyDPtY4mC8H426xeRxfU8SMQuicZarjIg5fLko&e=">https://drive.google.com/file/d/1F-hT7rM8ikzJVNMDmODWsCXKH391mnZF/view?usp=sharing</a><br></div><div><br></div><div>My physical constants are: (energy,momentum in GeV, length in m)</div><div><ul><li>mElectron = 0.510998910e-3;</li><li>hbarc = 0.1973269718e-15;</li><li>diamond lattice_constant = 3.5668e-10;</li><li>diamond thickness = 5e-5;</li><li>collimator diameter = 0.005;</li><li>radiator-collimator distance = 76;</li><li>endpoint energy = 12;</li><li>theta_small = 3.5e-3;</li></ul><div>From now on, I propose that we decide on a common naming scheme for these angles. I am happy to call them theta_ortho and theta_edge, if that works for you. With this settled, we should be ready to start production of some simulations. Please send me your order.</div><div><br></div><div>-Richard Jones</div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Mar 5, 2018 at 9:30 AM, Ken Livingston <span dir="ltr"><<a href="mailto:Kenneth.Livingston@glasgow.ac.uk" target="_blank">Kenneth.Livingston@glasgow.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div text="#000000" bgcolor="#FFFFFF">
<p>Hi Richard,</p>
<p>I'm on strike today (literally, not metaphorically), so won't
attend. Did you manage to have a look at the reason for the
discrepancy between your crystal angles and mine yet?<br>
</p>
Regards,<br>
Ken<div><div class="h5"><br>
<br>
<br>
<div class="m_-1035765733938230348moz-cite-prefix">On 05/03/18 10:29, Richard Jones wrote:<br>
</div>
</div></div><blockquote type="cite"><div><div class="h5">
<div dir="ltr">Dear colleagues,
<div><br>
</div>
<div>Please remember our photon beam meeting this morning
starting at 11:30. I have created an agenda page at the usual
place on the wiki. I hope that I will be able to arrive at
JLab in time for the start of the meeting. If not, please
start without me, I will come as soon as I get in.</div>
<div><br>
</div>
<div><a href="https://halldweb.jlab.org/wiki/index.php/BLTWG_Meeting_03/05/2018" target="_blank">https://halldweb.jlab.org/<wbr>wiki/index.php/BLTWG_Meeting_<wbr>03/05/2018</a><br>
</div>
<div><br>
</div>
<div>-Richard Jones</div>
</div>
<br>
<fieldset class="m_-1035765733938230348mimeAttachmentHeader"></fieldset>
<br>
</div></div><pre>______________________________<wbr>_________________
Halld-tagger mailing list
<a class="m_-1035765733938230348moz-txt-link-abbreviated" href="mailto:Halld-tagger@jlab.org" target="_blank">Halld-tagger@jlab.org</a>
<a class="m_-1035765733938230348moz-txt-link-freetext" href="https://mailman.jlab.org/mailman/listinfo/halld-tagger" target="_blank">https://mailman.jlab.org/<wbr>mailman/listinfo/halld-tagger</a></pre>
</blockquote>
<br>
</div>
</blockquote></div><br></div>