97 ''' look at cluster and new clusters'''
101 print(
"number of clusters = " + str(clusterList.getEntries()) +
" (old) VS " + str(clusterListNew.getEntries()) +
" (new)")
103 for d
in clusterList:
104 trueList = d.getRelationsTo(
'SVDTrueHits')
108 truePos = trueList[0].getU()
110 truePos = trueList[0].getV()
112 self.size.Fill(d.getSize())
113 self.time.Fill(d.getClsTime())
114 self.ff.Fill(d.getFirstFrame())
115 self.charge.Fill(d.getCharge() / 1000)
116 self.SNR.Fill(d.getSNR())
117 self.position.Fill(d.getPosition())
118 self.positionSigma.Fill(d.getPositionSigma() * 1e4)
120 self.positionS1.Fill(d.getPosition())
121 self.positionS1Sigma.Fill(d.getPositionSigma() * 1e4)
122 self.positionPull1.Fill((d.getPosition() - truePos) / d.getPositionSigma())
124 self.positionS2.Fill(d.getPosition())
125 self.positionS2Sigma.Fill(d.getPositionSigma() * 1e4)
126 self.positionPull2.Fill((d.getPosition() - truePos) / d.getPositionSigma())
128 self.positionS3.Fill(d.getPosition())
129 self.positionS3Sigma.Fill(d.getPositionSigma() * 1e4)
130 self.positionPull3.Fill((d.getPosition() - truePos) / d.getPositionSigma())
132 self.test[d.getSensorID().getLayerNumber() - 3].Fill(d.getSensorID().getLadderNumber(),
133 d.getSensorID().getSensorNumber() + isU)
134 self.positionPull.Fill((d.getPosition() - truePos) / d.getPositionSigma())
136 for d
in clusterListNew:
139 trueListNew = d.getRelationsTo(
'SVDTrueHits')
142 truePos = trueListNew[0].getU()
144 truePos = trueListNew[0].getV()
146 self.sizeNew.Fill(d.getSize())
147 self.timeNew.Fill(d.getClsTime())
148 self.ffNew.Fill(d.getFirstFrame())
149 self.chargeNew.Fill(d.getCharge() / 1000)
150 self.SNRNew.Fill(d.getSNR())
151 self.positionNew.Fill(d.getPosition())
152 self.positionSigmaNew.Fill(d.getPositionSigma() * 1e4)
154 self.positionS1New.Fill(d.getPosition())
155 self.positionS1SigmaNew.Fill(d.getPositionSigma() * 1e4)
156 self.positionPull1New.Fill((d.getPosition() - truePos) / d.getPositionSigma())
158 self.positionS2New.Fill(d.getPosition())
159 self.positionS2SigmaNew.Fill(d.getPositionSigma() * 1e4)
160 self.positionPull2New.Fill((d.getPosition() - truePos) / d.getPositionSigma())
162 self.positionS3New.Fill(d.getPosition())
163 self.positionS3SigmaNew.Fill(d.getPositionSigma() * 1e4)
164 self.positionPull3New.Fill((d.getPosition() - truePos) / d.getPositionSigma())
166 self.testNew[d.getSensorID().getLayerNumber() - 3].Fill(d.getSensorID().getLadderNumber(),
167 d.getSensorID().getSensorNumber() + isU)
168 self.positionPullNew.Fill((d.getPosition() - truePos) / d.getPositionSigma())
A (simplified) python wrapper for StoreArray.