39 return StatusCode::SUCCESS;
61 ATH_MSG_DEBUG (
" BEGIN EVENT ========================================== " );
69 size_t nEtaClusWidthCnt = 0, nPhiClusWidthCnt = 0;
72 int clusCount[33][9], sigclusCount[33][9];
73 for(
size_t kl = 0;
kl < 33;
kl++ )
75 for(
size_t km = 0;
km < 9;
km++ )
77 clusCount[
kl][
km] = 0;
78 sigclusCount[
kl][
km] = 0;
104 int ns = sectorNo < 0 ? sectorNo*(-1) : sectorNo+16;
105 int nl = (measuresPhi ? wireLayer : wireLayer+4);
110 ATH_MSG_DEBUG(
" ns = " <<
ns <<
"\tm_nl = " << nl <<
"\tm_sec = " << sectorNo <<
"\t m_lay= " << wireLayer <<
"\tmPhi = " << measuresPhi);
118 int xfac = measuresPhi ? -1 : 1;
123 std::vector<float> xVals;
128 if (measuresPhi) nPhiClusWidthCnt++;
129 else nEtaClusWidthCnt++;
130 fill(
"CscRdoMonitor", raw_clus_width, secLayer, measuresPhi_mon, measuresEta);
133 for (
size_t n = 0;
n < raw_clus_width;
n++ )
139 fill(
"CscRdoMonitor", stripId, secLayer);
142 xVals.push_back(stripId);
158 diff = n_max - n_min;
161 if(
diff > diff_max) diff_max =
diff;
163 ATH_MSG_DEBUG (
n <<
" Max = " << n_max <<
" Min = " << n_min <<
" Diff = " <<
diff );
167 ATH_MSG_DEBUG (
" End loop over strips======================" );
172 fill(
"CscRdoMonitor", diff_max_mon);
176 std::vector<int> tmp_xVals_signal;
177 std::vector<int> tmp_xVals_noise;
178 std::vector<float> tmp_secLayer_signal;
179 std::vector<float> tmp_secLayer_noise;
183 sigclusCount[
ns][nl]++;
184 for(
size_t nf = 0;
nf < xVals.size();
nf++)
186 tmp_xVals_signal.push_back(xVals[
nf]);
187 tmp_secLayer_signal.push_back(secLayer);
192 for(
size_t nf = 0;
nf < xVals.size();
nf++)
194 tmp_xVals_noise.push_back(xVals[
nf]);
195 tmp_secLayer_noise.push_back(secLayer);
215 fill(
"CscRdoMonitor", isSignal, isNoise, secLayer, tmp_secLayer_sig_mon, tmp_secLayer_noise_mon, tmp_xVals_signal_mon, tmp_xVals_noise_mon, raw_clus_width, sideC_signal, sideA_signal, signal_measuresEta, signal_measuresEta_sideC, signal_measuresEta_sideA, signal_measuresPhi, signal_measuresPhi_sideC, signal_measuresPhi_sideA, noise_measuresEta, noise_measuresPhi);
220 ATH_MSG_DEBUG (
" End loop over clusters======================" );
223 fill(
"CscRdoMonitor", nPhiClusWidthCnt_mon, nEtaClusWidthCnt_mon);
225 int numeta = 0, numphi = 0, numetasignal = 0, numphisignal = 0;
226 for(
int kl = 1;
kl < 33;
kl++ )
229 for(
int km = 1;
km < 9;
km++ )
231 int lay = (
km > 4 &&
km < 9) ?
km-4 :
km;
232 bool mphi = (
km > 0 &&
km < 5) ?
true :
false;
233 std::string wlay = mphi ?
"Phi-Layer " :
"Eta-Layer: ";
236 int tmp_counts_phi = 0, tmp_counts_eta = 0;
240 float secLayer =
kl-16 + 0.2 * (lay - 1) + 0.1;
249 else tmp_counts_phi =
count;
259 else tmp_counts_eta =
count;
273 fill(
"CscRdoMonitor", count_mon, secLayer_count, mphi_mon, mphi_scount, scount_mon, tmp_counts_phi_mon, mphi_inv, mphi_inv_scount, tmp_counts_eta_mon);
285 fill(
"CscRdoMonitor", numphi_mon, numeta_mon, numphisignal_mon, numetasignal_mon, diffnumphi, diffnumeta);
288 ATH_MSG_DEBUG (
" END EVENT ============================================ " );
290 ATH_MSG_DEBUG (
"CSCRdoMon::fillHistograms reports success" );
292 return StatusCode::SUCCESS;