55 if (not
m_dsps.isValid())
return;
61 if (
m_dsps.getEntries() == 8736) {
69 for (
auto& aECLDsp :
m_dsps) {
71 int nadc = aECLDsp.getNADCPoints();
72 int cellID = aECLDsp.getCellId();
74 int crysID = cellID - 1;
75 std::vector<int> dspAv = aECLDsp.getDspA();
78 double dspMean = (std::accumulate(dspAv.begin(), dspAv.begin() + nadc, 0.0)) / nadc;
80 for (
int v = 0; v < nadc; v++) {
81 wpsum += pow(dspAv[v] - dspMean, 2);
83 double dspRMS =
sqrt(wpsum / nadc);
87 double dspBkgRate = ((pow(dspSigma, 2)) - (pow(
m_waveformNoise[crysID], 2))) / (2.53 * 1e-12);
93 rates.averageDspBkgRate[segmentNumber] += dspBkgRate;
98 rates.validDspRate =
true;
122 for (
int cid = 1; cid < 8737; cid++) {
125 const double phi = position.
Phi();
126 const double z = position.
Z();
129 if (phi > 0.7853 && phi < 2.356) {
132 }
else if (phi >= 2.356 || phi <= -2.356) {
135 }
else if (phi > -2.356 && phi < -0.7853) {
142 }
else if (cid < 7777) {
144 if (phi > 0.7853 && phi < 2.356) {
147 }
else if (phi >= 2.356 || phi <= -2.356) {
150 }
else if (phi > -2.356 && phi < -0.7853) {
158 if (phi > 0.7853 && phi < 2.356) {
161 }
else if (phi >= 2.356 || phi <= -2.356) {
164 }
else if (phi > -2.356 && phi < -0.7853) {
173 if (phi > 0.7853 && phi < 2.356) {
176 }
else if (phi >= 2.356 || phi <= -2.356) {
179 }
else if (phi > -2.356 && phi < -0.7853) {