Hello, I wanted to add a little more information for anyone willing to look into this. The original tritrig file I used can be found at: /work/hallb/hps/ > mccaky/tritrigv2_5mrad_z-0pt5mm_HPS-PhysicsRun2016-v5-3-fieldmap_globalAlign_9999-spaced.slcio The readout and recon files can be found at: /work/hallb/hps/mccaky/tritrigReadout_Pair1.slcio > /work/hallb/hps/mccaky/tritrigRecon_Pair1.slcio The relevant collections in these files are: - Raw Hits >> Collection Name: "EcalReadoutHits"; Type: RawTrackerHit - Converted Hits >> Collection Name: "EcalCalHits"; Type: CalorimeterHit - GTP Clusters >> Collection Name: "EcalClustersGTP"; Type: Cluster - Clusters >> Collection Name: "EcalClusters "; Type: Cluster These files were created with > > /steering-files/src/main/resources/org/hps/steering/readout/PhysicsRun2016TrigPairs1.lcsim for readout (noise was disabled for the SVT and calorimeter) and /steering-files/src/main/resources/org/hps/steering/recon/PhysicsRun2016FullReconMC.lcsim with no changes for recon. Thanks, Kyle On Tue, Mar 27, 2018 at 2:28 PM, Kyle McCarty <[log in to unmask]> wrote: > 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