<div dir="ltr"><div>Kei,</div><div><br></div><div>You wrote,</div><div>| If we are not losing any particles, I don't understand<br>
</div>| how this can happen.<br><div><br></div><div>If you look back again at my earlier message, I explain how it can go either up or down. There is another more obvious explanation than extra hadrons from interactions of extra beam photons, and that is that the random numbers that seed all events (except the first one in the file) are completely different if you change anything about how the simulation proceeds. Thus, you can gain or lose tracks because of interactions, decays, etc. which happen (or don't happen) differently when the random number changes.</div>
<div><br></div><div>-Richard Jones</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On Fri, Mar 14, 2014 at 12:49 PM, Kei Moriya <<a href="mailto:kmoriya@indiana.edu">kmoriya@indiana.edu</a>> wrote:<br><blockquote class="gmail_quote">
<br>
Hi Richard,<br>
<br>
Within the program I check that the momentum difference between<br>
all particles I select and the original photon and proton<br>
is less than 1 MeV, and the difference in energy is less than<br>
0.1 MeV. I would assume, since as you say EM background<br>
does not produce many hadrons, that we wouldn't get any pions<br>
that satisfy these tight conditions. The other thing is that<br>
if you look at the numbers I got, the no EM setting<br>
has *more* events than with the background turned on.<br>
If we are not losing any particles, I don't understand<br>
how this can happen.<font color="#888888"><br>
<br>
Kei</font><div class="HOEnZb"><div class="h5"><br>
<br>
On 3/14/14 12:38 PM, Richard Jones wrote:<br>
<blockquote class="gmail_quote">
Kei,<br>
<br>
You wrote,<br>
| If people could give me comments on why I don't get<br>
| exactly the same number of events when running with<br>
| different EM backgrounds, I would be grateful.<br>
<br>
The short answer is that hdgeant appends new particles to the end of the<br>
event Monte Carlo record (Reaction tag) depending on the secondaries that<br>
are produced during tracking. You cannot lose particles from the Monte<br>
Carlo record, but you can gain more. This can either increase or decrease<br>
your event count in a given channel, either by moving events into the<br>
channel by adding particles that it was missing, allowing it to satisfy the<br>
particle counts for the channel, or by moving events out the channel by<br>
pushing the particle counts outside the ranges required for the channel.<br>
EM background is very poor in hadrons, so this explains why your different<br>
channel counts differ by only a couple dozen when you change the bg<br>
conditions.<br>
<br>
-Richard Jones<br>
<br>
<br>
On Fri, Mar 14, 2014 at 12:24 PM, Kei Moriya <<a href="mailto:kmoriya@indiana.edu">kmoriya@indiana.edu</a>> wrote:<br>
<br>
<blockquote class="gmail_quote">
<br>
I've updated the code to select specific reactions.<br>
The files are in<br>
trunk/sim-recon/src/programs/<u></u>Utilities/hddm_select_events<br>
<br>
You can run on non-REST files with<br>
hddm_select_events -s 5 -i [input hddm] -o [output hddm]<br>
To run over REST files use the -r option.<br>
The option -s [selectType] specifies which events you<br>
want to select. 5 is for p pi+ pi-, 6 is for p pi+ pi- pi0.<br>
<br>
If people could give me comments on why I don't get<br>
exactly the same number of events when running with<br>
different EM backgrounds, I would be grateful.<br>
<br>
Kei<br>
______________________________<u></u>_________________<br>
Halld-offline mailing list<br>
<a href="mailto:Halld-offline@jlab.org">Halld-offline@jlab.org</a><br>
<a href="https://mailman.jlab.org/mailman/listinfo/halld-offline">https://mailman.jlab.org/<u></u>mailman/listinfo/halld-offline</a><br>
<br>
</blockquote>
<br>
</blockquote>
</div></div></blockquote></div><br></div>