Hello hps-software,

I have recently noticed something odd with reconstruction for Monte Carlo, and wanted to see if any recon experts could comment on it.

I ran readout over a pure tri-trig MC sample (with ~500 empty events spacing to prevent overlap) using the 2016 pair 1 steering file.

/steering-files/src/main/resources/org/hps/steering/readout/PhysicsRun2016TrigPairs1.lcsim

The only change that I made was to disable noise on the SVT and calorimeter.

Doing this produces results as I would expect, and I see that every event registers at least two GTP clusters during readout, as it must for a pair trigger.

However, when I run recon, I find that many events have either no clusters or only one cluster after recon, and that the recon energies are often a lot lower than the readout energies. Consider for example a few of the following event logs:

Event 1276
Raw Hits:
Hit at  126989 and time    0 ns and position < -190.1,    30.0,  1472.1>.
ADC:   120   120   120   120   120   120   120   120   120   120   120   129   149   165   175   178   175   170   163   156   150   144   139   135   131   128   126   125   123   123   122   121   121   121   121   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120
Hit at  127245 and time    0 ns and position < -174.8,    30.0,  1472.3>.
ADC:   120   120   120   120   120   120   120   120   120   120   120   123   131   137   141   142   141   139   137   134   131   129   127   126   124   123   122   122   121   121   121   121   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120
Hit at 4067597 and time    0 ns and position <  290.5,   -45.1,  1471.9>.
ADC:   120   120   120   120   120   120   120   120   120   120   120   123   130   136   140   141   140   138   136   133   131   129   127   125   124   123   122   122   121   121   121   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120
Hit at 4133133 and time    0 ns and position <  290.5,   -30.1,  1471.9>.
ADC:   120   120   120   120   120   120   120   120   120   120   120   127   144   159   167   169   167   163   157   151   146   140   136   132   130   127   125   124   123   122   122   121   121   121   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120
Converted Hits:
Hit at  126989 with energy  0.124 GeV ( 0.124 GeV) at time  -0.78 ns and position < -190.1,    30.0,  1472.1>.
Truth Hits: 1
Hit at  127245 with energy  0.078 GeV ( 0.078 GeV) at time   2.37 ns and position < -174.8,    30.0,  1472.3>.
Truth Hits: 1
Hit at 4067597 with energy  0.025 GeV ( 0.025 GeV) at time  28.39 ns and position <  290.5,   -45.1,  1471.9>.
Truth Hits: 7
Hit at 4133133 with energy  0.071 GeV ( 0.071 GeV) at time  40.91 ns and position <  290.5,   -30.1,  1471.9>.
Truth Hits: 7
GTP Clusters:
Cluster at  126989 with energy  0.703 GeV and 2 hits at time 1352007.75 ns and position < -185.9,    30.0,  1472.2>.
Cluster at 4133133 with energy  0.617 GeV and 2 hits at time 1352007.25 ns and position <  290.5,   -34.5,  1471.9>.
Clusters:
None


Event 1787
Raw Hits:
Hit at   66829 and time    0 ns and position <  109.5,    30.0,  1473.4>.
ADC:   120   120   120   120   120   120   120   120   120   120   120   131   155   175   186   190   187   180   172   164   156   149   143   137   133   130   128   126   124   123   122   122   121   121   121   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120
Hit at  132365 and time    0 ns and position <  109.5,    45.0,  1473.5>.
ADC:   120   120   120   120   120   120   120   120   120   120   120   125   137   147   153   154   153   150   146   142   138   134   131   129   127   125   124   123   122   122   121   121   121   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120
Hit at 4061453 and time    0 ns and position <  -54.3,   -60.1,  1473.5>.
ADC:   120   120   120   120   120   120   120   120   120   120   120   124   134   141   146   147   146   143   140   137   134   131   129   127   125   124   123   122   122   121   121   121   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120
Hit at 4126989 and time    0 ns and position <  -54.3,   -45.1,  1473.4>.
ADC:   120   120   120   120   120   120   120   120   120   120   120   129   150   167   177   180   177   172   165   158   151   145   139   135   131   129   127   125   124   123   122   121   121   121   121   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120
Converted Hits:
Hit at   66829 with energy  0.081 GeV ( 0.081 GeV) at time  38.48 ns and position <  109.5,    30.0,  1473.4>.
Truth Hits: 2
Hit at  132365 with energy  0.052 GeV ( 0.052 GeV) at time  38.02 ns and position <  109.5,    45.0,  1473.5>.
Truth Hits: 2
Hit at 4061453 with energy  0.037 GeV ( 0.037 GeV) at time  37.93 ns and position <  -54.3,   -60.1,  1473.5>.
Truth Hits: 1
Hit at 4126989 with energy  0.170 GeV ( 0.170 GeV) at time  -0.77 ns and position <  -54.3,   -45.1,  1473.4>.
Truth Hits: 1
GTP Clusters:
Cluster at 4126989 with energy  0.768 GeV and 2 hits at time 1863007.88 ns and position <  -54.3,   -49.7,  1473.5>.
Cluster at   66829 with energy  0.921 GeV and 2 hits at time 1863008.00 ns and position <  109.5,    34.9,  1473.4>.
Clusters:
Cluster at   66829 with energy  0.133 GeV and 2 hits at time  38.48 ns and position <  103.5,    35.5,  1393.6>.



Event 2477
Raw Hits:
Hit at  328973 and time    0 ns and position <  109.5,    90.0,  1473.6>.
ADC:   120   120   120   120   120   120   120   120   120   120   120   129   150   168   177   180   178   172   165   158   151   145   140   135   132   129   127   125   124   123   122   121   121   121   121   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120
Hit at  329229 and time    0 ns and position <  124.4,    90.0,  1473.6>.
ADC:   120   120   120   120   120   120   120   120   120   120   120   124   133   140   144   146   145   142   139   136   133   131   128   126   125   124   123   122   122   121   121   121   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120
Hit at 3996685 and time    0 ns and position <   -9.6,   -75.1,  1473.6>.
ADC:   120   120   120   120   120   120   120   120   120   120   120   128   146   160   168   171   169   164   158   152   146   141   136   133   130   127   126   124   123   122   122   121   121   121   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120
Hit at 4062221 and time    0 ns and position <   -9.6,   -60.1,  1473.5>.
ADC:   120   120   120   120   120   120   120   120   120   120   120   125   137   147   153   155   153   150   146   142   138   134   131   129   127   125   124   123   122   122   121   121   121   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120   120
Converted Hits:
Hit at  328973 with energy  0.073 GeV ( 0.073 GeV) at time  40.86 ns and position <  109.5,    90.0,  1473.6>.
Truth Hits: 5
Hit at  329229 with energy  0.112 GeV ( 0.112 GeV) at time  -0.34 ns and position <  124.4,    90.0,  1473.6>.
Truth Hits: 3
Hit at 3996685 with energy  0.178 GeV ( 0.178 GeV) at time  -0.18 ns and position <   -9.6,   -75.1,  1473.6>.
Truth Hits: 3
Hit at 4062221 with energy  0.044 GeV ( 0.044 GeV) at time  41.09 ns and position <   -9.6,   -60.1,  1473.5>.
Truth Hits: 2
GTP Clusters:
Cluster at  328973 with energy  0.760 GeV and 2 hits at time 2553008.00 ns and position <  114.0,    90.0,  1473.6>.
Cluster at 3996685 with energy  0.755 GeV and 2 hits at time 2553007.25 ns and position <   -9.6,   -69.0,  1473.5>.
Clusters:
None

This is a bit confusing to me - I would expect that we should always see at least two clusters. I suspect that this is related to whatever is causing the drop in energy. This seems to occur on the hit level, and presumably some of the clusters and/or hits are coming up below whatever the cut-off threshold is for inclusion and thus being dropped. I'm not completely sure, though. Can any recon experts comment on what is going on here?

Thanks,

Kyle