
<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">In this particular event, are there any particles in DMCThrown that are shown as decaying from this pion? &nbsp;Or are they truncated from the DMCThrown output? &nbsp;Is it possible / would it be useful to re-enable their output for debugging purposes? &nbsp;Do/would they have a time/location in DMCThrown that would correlate with the extra hits? &nbsp;<div><br></div><div>&nbsp;- Paul</div><div><br></div><div><div>On May 29, 2014, at 6:47 PM, Richard Jones wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div dir="ltr">Hello Will and all,<div><br></div><div>Once a primary particle has entered the BCAL and had its TruthShower recorded, it is marked with a spot of red dye (colloquially). Any secondaries that are generated from interactions by that particle inherit its red spot, even if they escape from that module and enter another one. &nbsp;However, these late interactions typically are generated by delayed decays of muons from pi decays, aren't they? &nbsp;Typically that is what makes late hits in a calorimeter, if my memory serves me correctly.</div>

<div><br></div><div>-Richard J.</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, May 29, 2014 at 6:32 PM, Will Levine &lt;<a href="mailto:wilevine@andrew.cmu.edu">wilevine@andrew.cmu.edu</a>&gt; wrote:<br>

<blockquote class="gmail_quote">Hi offline,<br>

Following up on the issue of "late" energy deposition in hadronic<br>

showers in the BCAL.<br>

<br>

There are two different ways that BCAL truth information from hdgeant<br>

is stored: these can be accessed in JANA from the DBCALTruthShower and<br>

the DBCALTruthCell objects. The code that stores this information is<br>

in hitBarrelEMcal() in<br>

sim-recon/src/programs/Simulation/HDGeant/hitBCal.cc<br>

<br>

As I understand the code, DBCALTruthShower records every particle that<br>

geant tracks thru the BCAL (with energy greater than 1 MeV). One<br>

DBCALTruthShower for each unique particle in the BCAL. This includes<br>

all secondary shower particles. The particle type, initial energy, and<br>

initial position in the BCAL is recorded for each particle.<br>

<br>

DBCALTruthCell just records the energy deposited in a given BCAL cell<br>

without recording which particle was responsible for the energy.<br>

(Again this is "truth" information, before any detector resolution or<br>

readout effects are incorporated).<br>

<br>

These two sources of information are not always tremendously helpful.<br>

For example in one single-pion event that has some characteristic<br>

"late hits", we see the following DBCALTruthShower's:<br>

<br>

DBCALTruthShower:<br>

&nbsp;ptype: track: primary: &nbsp; &nbsp;phi: &nbsp; &nbsp; &nbsp;r: &nbsp; &nbsp; z: &nbsp; &nbsp; t: &nbsp; &nbsp; p: &nbsp; &nbsp; E:<br>

-------------------------------------------------------------------<br>

&nbsp; &nbsp; &nbsp;8 &nbsp; &nbsp; &nbsp;1 &nbsp; &nbsp; &nbsp; &nbsp;1 &nbsp;-1.863 &nbsp;65.042 &nbsp;185.3 &nbsp;4.780 &nbsp;0.571 &nbsp;0.588<br>

<br>

Only the primary pion is recorded. There is no record of any secondary<br>

tracks. But when we look at the TruthCell's in the same event:<br>

<br>

DBCALTruthCell [truncated]:<br>

&nbsp;module: layer: sector: &nbsp; &nbsp; E: &nbsp; &nbsp; &nbsp; t: zLocal:<br>

-----------------------------------------------<br>

&nbsp; &nbsp; &nbsp;34 &nbsp; &nbsp; &nbsp;1 &nbsp; &nbsp; &nbsp; 1 &nbsp;0.006 &nbsp; &nbsp;71.31 &nbsp; &nbsp;-6.0<br>

&nbsp; &nbsp; &nbsp;34 &nbsp; &nbsp; &nbsp;1 &nbsp; &nbsp; &nbsp; 3 &nbsp;0.027 &nbsp; &nbsp; 4.87 &nbsp; -24.4<br>

&nbsp; &nbsp; &nbsp;34 &nbsp; &nbsp; &nbsp;1 &nbsp; &nbsp; &nbsp; 4 &nbsp;0.002 &nbsp; &nbsp; 4.79 &nbsp; -26.5<br>

&nbsp; &nbsp; &nbsp;34 &nbsp; &nbsp; &nbsp;2 &nbsp; &nbsp; &nbsp; 1 &nbsp;0.006 &nbsp; &nbsp;33.59 &nbsp; -16.9<br>

&nbsp; &nbsp; &nbsp;34 &nbsp; &nbsp; &nbsp;2 &nbsp; &nbsp; &nbsp; 2 &nbsp;0.018 &nbsp; &nbsp; 5.07 &nbsp; -19.6<br>

&nbsp; &nbsp; &nbsp;34 &nbsp; &nbsp; &nbsp;2 &nbsp; &nbsp; &nbsp; 3 &nbsp;0.013 &nbsp; &nbsp; 4.98 &nbsp; -21.7<br>

&nbsp; &nbsp; &nbsp;34 &nbsp; &nbsp; &nbsp;3 &nbsp; &nbsp; &nbsp; 1 &nbsp;0.020 &nbsp; &nbsp;28.91 &nbsp; -15.4<br>

&nbsp; &nbsp; &nbsp;34 &nbsp; &nbsp; &nbsp;3 &nbsp; &nbsp; &nbsp; 1 &nbsp;0.005 &nbsp;4301.28 &nbsp; -19.1<br>

&nbsp; &nbsp; &nbsp;34 &nbsp; &nbsp; &nbsp;3 &nbsp; &nbsp; &nbsp; 2 &nbsp;0.023 &nbsp; &nbsp; 5.19 &nbsp; -16.6<br>

&nbsp; &nbsp; &nbsp;34 &nbsp; &nbsp; &nbsp;4 &nbsp; &nbsp; &nbsp; 1 &nbsp;0.029 &nbsp; &nbsp; 6.45 &nbsp; -13.1<br>

&nbsp; &nbsp; &nbsp;34 &nbsp; &nbsp; &nbsp;4 &nbsp; &nbsp; &nbsp; 1 &nbsp;0.019 &nbsp; &nbsp;73.07 &nbsp; -13.0<br>

&nbsp; &nbsp; &nbsp;34 &nbsp; &nbsp; &nbsp;4 &nbsp; &nbsp; &nbsp; 1 &nbsp;0.004 &nbsp;2219.77 &nbsp; -13.3<br>

&nbsp; &nbsp; &nbsp;34 &nbsp; &nbsp; &nbsp;4 &nbsp; &nbsp; &nbsp; 1 &nbsp;0.039 &nbsp;4301.12 &nbsp; -15.0<br>

&nbsp; &nbsp; &nbsp;34 &nbsp; &nbsp; &nbsp;5 &nbsp; &nbsp; &nbsp; 1 &nbsp;0.017 &nbsp; &nbsp; 8.24 &nbsp; &nbsp;-9.6<br>

&nbsp; &nbsp; &nbsp;34 &nbsp; &nbsp; &nbsp;5 &nbsp; &nbsp; &nbsp; 1 &nbsp;0.001 &nbsp;4301.19 &nbsp; -17.1<br>

&nbsp; &nbsp; &nbsp;34 &nbsp; &nbsp; &nbsp;7 &nbsp; &nbsp; &nbsp; 1 &nbsp;0.011 &nbsp; &nbsp;10.44 &nbsp; -21.4<br>

<br>

The time is in the second-to-last column. The "on-time" hits have t ~<br>

5 ns. Where are all the late hits coming from? It is not clear. There<br>

must be some tracks that are not recorded in the TruthShower info.<br>

<br>

Even in events with a more extensive record of secondary tracks it's<br>

not clear what is responsible for a late hit:<br>

<br>

DBCALTruthShower:<br>

&nbsp;ptype: track: primary: &nbsp; &nbsp;phi: &nbsp; &nbsp; &nbsp;r: &nbsp; &nbsp; z: &nbsp; &nbsp; t: &nbsp; &nbsp; p: &nbsp; &nbsp; E:<br>

-------------------------------------------------------------------<br>

&nbsp; &nbsp; &nbsp;8 &nbsp; &nbsp; &nbsp;1 &nbsp; &nbsp; &nbsp; &nbsp;1 &nbsp;-2.072 &nbsp;65.042 &nbsp;186.4 &nbsp;4.812 &nbsp;0.572 &nbsp;0.589<br>

&nbsp; &nbsp; 13 &nbsp; &nbsp; 22 &nbsp; &nbsp; &nbsp; &nbsp;0 &nbsp;-2.083 &nbsp;65.969 &nbsp;188.5 &nbsp;5.185 &nbsp;0.284 &nbsp;0.982<br>

&nbsp; &nbsp; 14 &nbsp; &nbsp; 21 &nbsp; &nbsp; &nbsp; &nbsp;0 &nbsp;-2.083 &nbsp;65.969 &nbsp;188.5 &nbsp;5.185 &nbsp;0.040 &nbsp;0.939<br>

&nbsp; &nbsp; 13 &nbsp; &nbsp; 20 &nbsp; &nbsp; &nbsp; &nbsp;0 &nbsp;-2.083 &nbsp;65.969 &nbsp;188.5 &nbsp;5.185 &nbsp;0.173 &nbsp;0.955<br>

&nbsp; &nbsp; &nbsp;3 &nbsp; &nbsp; 19 &nbsp; &nbsp; &nbsp; &nbsp;0 &nbsp;-2.081 &nbsp;65.814 &nbsp;188.1 &nbsp;4.879 &nbsp;0.005 &nbsp;0.005<br>

&nbsp; &nbsp; &nbsp;3 &nbsp; &nbsp; 19 &nbsp; &nbsp; &nbsp; &nbsp;0 &nbsp;-2.083 &nbsp;65.851 &nbsp;188.2 &nbsp;4.884 &nbsp;0.002 &nbsp;0.002<br>

&nbsp; &nbsp; &nbsp;1 &nbsp; &nbsp; 19 &nbsp; &nbsp; &nbsp; &nbsp;0 &nbsp;-2.083 &nbsp;65.852 &nbsp;188.2 &nbsp;4.884 &nbsp;0.001 &nbsp;0.001<br>

<br>

Any ideas?<br>

<br>

Will<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>

_______________________________________________<br>Halld-offline mailing list<br><a href="mailto:Halld-offline@jlab.org">Halld-offline@jlab.org</a><br>https://mailman.jlab.org/mailman/listinfo/halld-offline</blockquote></div><br></body></html>