<div dir="ltr">Hello Igal et. al.,<div><br></div><div>Please create an issue on github describing this problem, and I will submit a fix in a side branch. Once you have verified that this works, it will be merged into the master branch.</div><div><br></div><div>-Richard Jones</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jan 26, 2021 at 2:56 PM Igal Jaegle <<a href="mailto:ijaegle@jlab.org">ijaegle@jlab.org</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">



<div dir="ltr">
<p></p>
<div style="background-color:rgb(255,235,156);width:100%;border-style:none;border-color:rgb(250,235,204);border-width:1pt;padding:10pt;font-size:11pt;line-height:12pt;font-family:Calibri;color:rgb(0,0,0);text-align:left">
<span style="color:rgb(156,101,0)"></span>*Message sent from a system outside of UConn.*</div>
<br>
<p></p>
<div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Below an e-mail of Drew noting a bug in the tagger counter.</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
tks ig.<br>
</div>
<div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div id="gmail-m_-7682848818724345140appendonsend"></div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_-7682848818724345140divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> Drew Smith <<a href="mailto:andrew.p.smith@duke.edu" target="_blank">andrew.p.smith@duke.edu</a>><br>
<b>Sent:</b> Tuesday, January 26, 2021 2:34 PM<br>
<b>To:</b> Alexander Somov <<a href="mailto:somov@jlab.org" target="_blank">somov@jlab.org</a>>; Igal Jaegle <<a href="mailto:ijaegle@jlab.org" target="_blank">ijaegle@jlab.org</a>>; Tyler J Hague <<a href="mailto:tjhague@ncat.edu" target="_blank">tjhague@ncat.edu</a>><br>
<b>Subject:</b> [EXTERNAL] No TAGH Hits in Simulation</font>
<div> </div>
</div>
<div dir="ltr">
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Hi Guys,</div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
I believe I found the reason why there are no tagger hodoscope hits produced in my simulations for the region below the microscope minimum energy. As you recall, we decided to set the scaled_energy_range values for hodoscope counters 128-178 (the counters whose
 energies would otherwise span the microscope region) to 0. However, this introduces a problem in the 'addTaggerPhoton' function in hdgeant4 (defined in GlueXDetectorTAG.cc) which is responsible for generating tagger hits for beam photons. Specifically, if
 you look at lines 138-149 of GlueXDetectorTAG.cc it becomes clear that because this way of setting the scaled_energy_range table for the hodoscope will not work.</div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
For example, if there is a beam photon at 6 GeV, these lines of code will give a starting guess for the corresponding hodoscope counter of (line 139): int i = 274 * (6 GeV - 11.6061 GeV) / (0 GeV - 11.6061 GeV) = 132 (this particular example uses the endpoint
 energy for run 61321). Next, the code compares whether the beam photon is less than the minimum energy for this tagger counter. If so, then it increments the counter number by 1 and re-checks this until it finds a counter where this is no longer true. However,
 for counter 132 the minimum energy is 0 GeV because the scaled_energy_range is set to 0. So, it does not increase the counter number until the actual counter is found, and in the end no tagger hits are constructed.</div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Arial,Helvetica,sans-serif;font-size:12pt"><br>
</span></div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Arial,Helvetica,sans-serif;font-size:12pt">The question then is how can we correct this? If the scaled_energy_range table is reverted to its previous state, then this will also be a problem (albeit a lesser problem).
 For example, with the old scaled_energy_table, if I generate events with gen_primex_compton for tagger hodoscope counter 200, then the photons will not be reconstructed in this counter but would instead be reconstructed in one or more of the counters between
 128 and 178. The photon beam energy would be approximately right, but we would lose quite a lot of events because of gaps between adjacent counters. </span></div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Arial,Helvetica,sans-serif;font-size:12pt"><br>
</span></div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Arial,Helvetica,sans-serif;font-size:12pt">Alternatively, I can make a change to the hdgeant4 code where the hodoscope hits are registered to avoid this specific problem and allow us to keep counters 128-178 set to
 0. However, I'm sure other GlueX people would not want me making changes to their simulation code, so I would probably only want to make this change to my own local installation of hdgeant4. </span></div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Arial,Helvetica,sans-serif;font-size:12pt"><br>
</span></div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Does anyone have suggestions on how to best handle this problem? </div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Igal,</div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
In your simulations for the eta photoproduction, did you see tagger hodoscope hits registered for beam energies between ~4-7 GeV? I think below 4 GeV there should be no problem, and above the microscope energies there should also be no problem. </div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Arial,Helvetica,sans-serif;font-size:12pt"><br>
</span></div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Arial,Helvetica,sans-serif;font-size:12pt"><br>
</span></div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Arial,Helvetica,sans-serif;font-size:12pt">Best,</span></div>
<div style="font-family:Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Arial,Helvetica,sans-serif;font-size:12pt">-- Drew </span></div>
</div>
</div>
</div>
</div>

</blockquote></div>