<div dir="ltr">Kei,<div><br></div><div>You wrote, </div><div>|  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>
</div><div><br></div><div>The short answer is that hdgeant appends new particles to the end of the event Monte Carlo record (Reaction tag) depending on the secondaries that are produced during tracking.  You cannot lose particles from the Monte Carlo record, but you can gain more.  This can either increase or decrease your event count in a given channel, either by moving events into the channel by adding particles that it was missing, allowing it to satisfy the particle counts for the channel, or by moving events out the channel by pushing the particle counts outside the ranges required for the channel.  EM background is very poor in hadrons, so this explains why your different channel counts differ by only a couple dozen when you change the bg conditions.</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:24 PM, Kei Moriya <<a href="mailto:kmoriya@indiana.edu">kmoriya@indiana.edu</a>> wrote:<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/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>
_______________________________________________<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/mailman/listinfo/halld-offline</a><br>
</blockquote></div><br></div>