 |
 |
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 ######################################################################## Use REPLY-ALL to reply to list To unsubscribe from the HPS-SOFTWARE list, click the following link: https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=HPS-SOFTWARE&A=1