MantisBT - Hall D Offline
|View Issue Details|
|ID||Project||Category||View Status||Date Submitted||Last Update|
|0000103||Hall D Offline||General||public||2011-06-29 11:10||2011-07-18 09:36|
|Summary||0000103: Reconstructed number of particles non-deterministic|
|Description||From e-mail by Will Levine on 6/28:|
I noticed recently that I've been getting different results when running reconstruction multiple times on the same file. It seems like that for a small number of events, the number of particles reconstructed varies. I think this problem will occur with any input, but I've uploaded one particular hddm file where I've seen the problem (http://www-meg.phys.cmu.edu/wilevine/hdgeant_smeared.hddm [^] ). I think the simplest way to see the problem is to run the reconstruction with the phys_tree plugin twice (i.e.
hd_root -PPLUGINS=phys_tree hdgeant_smeared.hddm -o output.1.root
hd_root -PPLUGINS=phys_tree hdgeant_smeared.hddm -o output.2.root
) and compare the output. In particular, differences should be seen in the Npip, Nphoton, etc. histograms.
For the record I'm using svn revision 7974 and the problem seems to occur in both SIMD and non-SIMD builds.
If anyone can confirm or diagnose this issue, that would be helpful.
|Steps To Reproduce|
|Additional Information||Beni reproduced the issue and sent this email:|
Hi Will, David,
I duplicated what you see. I look for each event how many photons are reconstructed. the result is
shown in the two attached files. The diff_sidebyside.dat file shows for each event the number of
reconstructed photons on the left for the first reconstruction on the right for the second reconstruction.
the diff.dat file contains only those event where a difference is found. This should help for further
|Tags||No tags attached.|
|2011-06-29 11:10||davidl||New Issue|
|2011-06-29 14:44||pmatt||Status||new => assigned|
|2011-06-29 14:44||pmatt||Assigned To||=> pmatt|
|2011-07-18 09:36||davidl||Note Added: 0000146|
|2011-07-18 09:36||davidl||Status||assigned => resolved|
|2011-07-18 09:36||davidl||Resolution||open => fixed|