<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Aptos;
panose-1:2 11 0 4 2 2 2 2 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
font-size:12.0pt;
font-family:"Aptos",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
span.EmailStyle19
{mso-style-type:personal-reply;
font-family:"Aptos",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;
mso-ligatures:none;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-US" link="blue" vlink="purple" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt">Hi all,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Dave Gaskell setup a new simc branch for us for generating the distributions and rates for the electron arm DIS electron.
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">I am passing this along to anyone who wanted to play with it. I will use this to start generating the comparison plots for our kinematics in the spectrometer quantities.
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">-Holly<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<div id="mail-editor-reference-message-container">
<div>
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal" style="margin-bottom:12.0pt"><b><span style="color:black">From:
</span></b><span style="color:black">Dave Gaskell <gaskelld@jlab.org><br>
<b>Date: </b>Monday, December 9, 2024 at 4:22</span><span style="font-family:"Arial",sans-serif;color:black"> </span><span style="color:black">PM<br>
<b>To: </b>Holly Szumila-Vance <hszumila@jlab.org><br>
<b>Subject: </b>electron singles rates<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt">Hi Holly -<br>
<br>
I made a branch of mc-single-arm that adds the cross section/weight to <br>
the output so one can use it similar to simc to get the rates/yield.<br>
<br>
<a href="https://github.com/gaskelld/mc-single-arm/tree/LAD">https://github.com/gaskelld/mc-single-arm/tree/LAD</a><br>
<br>
I've generated a cross section table using externals for the LAD <br>
kinematics, and interpolate the values in the table to get a cross <br>
section to put in the tree. New variables in the tree are:<br>
<br>
xsecv - vertex cross section in ub/sr/GeV<br>
xsecr - radiated cross section in ub/sr/GeV<br>
weight = xsecr*jacobian<br>
<br>
where the jacobian is the small transformation from cartesian <br>
(xprime,yprime) to spherical angles.<br>
<br>
To get the yield, just weight all the events by *weight*normfac", where <br>
"normfac" is printed out in the text output file. That will give you <br>
the number of counts for 1 mC of charge.<br>
<br>
For the HMS/SHMS kinematics you sent me, I get:<br>
<br>
P=4.4 GeV<br>
<br>
HMS:<br>
theta = 17 deg: 244 Hz<br>
theta=13.5 deg: 899 Hz<br>
<br>
SHMS<br>
theta=17 deg: 277 Hz<br>
theta=13.5 deg: 1050 Hz<br>
<br>
for 1 uA of beam current.<br>
<br>
I checked these values against my stand-alone program that ignores <br>
radiative corrections, and they are in same ballpark.<br>
<br>
Cheers,<br>
Dave<o:p></o:p></span></p>
</div>
</div>
</div>
</div>
</div>
</body>
</html>