71 int useAllHistos = 0, subdet =5;
78 MonitorElement* ZS_HE1=0, *ZS_seqHE1=0, *ZS_HE2=0, *ZS_seqHE2=0, *ZS_HE3=0, *ZS_seqHE3=0;
79 MonitorElement* map_depth1 =0, *map_depth2 =0, *map_depth3 =0, *map_depth4 =0;
82 MonitorElement* emap_depth1 =0, *emap_depth2 =0, *emap_depth3 =0, *emap_depth4 =0;
84 MonitorElement* occupancy_seqHE1 =0, *occupancy_seqHE2 =0, *occupancy_seqHE3 =0;
88 MonitorElement* emean_seqHE1 =0, *emean_seqHE2 =0, *emean_seqHE3 =0;
100 MonitorElement* occupancy_map_HE1 =0, *occupancy_map_HE2 =0, *occupancy_map_HE3 =0;
107 MonitorElement* RMS_vs_ieta_HE1 =0, *RMS_vs_ieta_HE2 =0, *RMS_vs_ieta_HE3 =0;
110 MonitorElement* occupancy_vs_ieta_HB1 =0, *occupancy_vs_ieta_HB2 =0;
111 MonitorElement* occupancy_vs_ieta_HE1 =0, *occupancy_vs_ieta_HE2 =0, *occupancy_vs_ieta_HE3 =0;
112 MonitorElement* occupancy_vs_ieta_HF1 =0, *occupancy_vs_ieta_HF2 =0;
115 MonitorElement* RecHit_StatusWord_HB =0, *RecHit_StatusWord_HE=0, *RecHit_StatusWord_HO =0, *RecHit_StatusWord_HF =0, *RecHit_StatusWord_HF67 =0;
116 MonitorElement* RecHit_Aux_StatusWord_HB =0, *RecHit_Aux_StatusWord_HE=0, *RecHit_Aux_StatusWord_HO =0, *RecHit_Aux_StatusWord_HF =0;
118 for(
unsigned int ih=0; ih<hcalMEs.size(); ih++){
119 if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_map_depth1") ==0 ){
120 useAllHistos =1; subdet =6;
122 if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_vs_ieta_HB1") ==0 ){
126 if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_simple1D_HO") ==0 ){ ZS_HO = hcalMEs[ih]; }
127 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_sequential1D_HO") ==0 ){ ZS_seqHO = hcalMEs[ih]; }
128 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_simple1D_HB1") ==0 ){ ZS_HB1 = hcalMEs[ih]; }
129 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_sequential1D_HB1") ==0 ){ ZS_seqHB1 = hcalMEs[ih]; }
130 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_simple1D_HB2") ==0 ){ ZS_HB2 = hcalMEs[ih]; }
131 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_sequential1D_HB2") ==0 ){ ZS_seqHB2 = hcalMEs[ih]; }
132 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_simple1D_HF1") ==0 ){ ZS_HF1 = hcalMEs[ih]; }
133 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_sequential1D_HF1") ==0 ){ ZS_seqHF1 = hcalMEs[ih]; }
134 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_simple1D_HF2") ==0 ){ ZS_HF2 = hcalMEs[ih]; }
135 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_sequential1D_HF2") ==0 ){ ZS_seqHF2 = hcalMEs[ih]; }
136 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_simple1D_HE1") ==0 ){ ZS_HE1 = hcalMEs[ih]; }
137 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_sequential1D_HE1") ==0 ){ ZS_seqHE1 = hcalMEs[ih]; }
138 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_simple1D_HE2") ==0 ){ ZS_HE2 = hcalMEs[ih]; }
139 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_sequential1D_HE2") ==0 ){ ZS_seqHE2 = hcalMEs[ih]; }
140 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_simple1D_HE3") ==0 ){ ZS_HE3 = hcalMEs[ih]; }
141 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_sequential1D_HE3") ==0 ){ ZS_seqHE3 = hcalMEs[ih]; }
142 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_map_depth1") ==0 ){ map_depth1= hcalMEs[ih]; }
143 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_map_depth2") ==0 ){ map_depth2= hcalMEs[ih]; }
144 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_map_depth3") ==0 ){ map_depth3= hcalMEs[ih]; }
145 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_map_depth4") ==0 ){ map_depth4= hcalMEs[ih]; }
147 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"N_HF") ==0 ){ Nhf= hcalMEs[ih]; }
148 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emap_depth1") ==0 ){ emap_depth1= hcalMEs[ih]; }
149 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emap_depth2") ==0 ){ emap_depth2= hcalMEs[ih]; }
150 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emap_depth3") ==0 ){ emap_depth3= hcalMEs[ih]; }
151 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emap_depth4") ==0 ){ emap_depth4= hcalMEs[ih]; }
153 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occ_sequential1D_HB1") ==0 ){ occupancy_seqHB1= hcalMEs[ih]; }
154 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occ_sequential1D_HB2") ==0 ){ occupancy_seqHB2= hcalMEs[ih]; }
155 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occ_sequential1D_HE1") ==0 ){ occupancy_seqHE1= hcalMEs[ih]; }
156 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occ_sequential1D_HE2") ==0 ){ occupancy_seqHE2= hcalMEs[ih]; }
157 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occ_sequential1D_HE3") ==0 ){ occupancy_seqHE3= hcalMEs[ih]; }
158 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occ_sequential1D_HF1") ==0 ){ occupancy_seqHF1= hcalMEs[ih]; }
159 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occ_sequential1D_HF2") ==0 ){ occupancy_seqHF2= hcalMEs[ih]; }
160 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occ_sequential1D_HO") ==0 ){ occupancy_seqHO= hcalMEs[ih]; }
161 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_seq_HB1") ==0 ){ emean_seqHB1= hcalMEs[ih]; }
162 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_seq_HB2") ==0 ){ emean_seqHB2= hcalMEs[ih]; }
163 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_seq_HE1") ==0 ){ emean_seqHE1= hcalMEs[ih]; }
164 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_seq_HE2") ==0 ){ emean_seqHE2= hcalMEs[ih]; }
165 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_seq_HE3") ==0 ){ emean_seqHE3= hcalMEs[ih]; }
166 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_seq_HF1") ==0 ){ emean_seqHF1= hcalMEs[ih]; }
167 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_seq_HF2") ==0 ){ emean_seqHF2= hcalMEs[ih]; }
168 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_seq_HO") ==0 ){ emean_seqHO= hcalMEs[ih]; }
169 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_seq_HB1") ==0 ){ RMS_seq_HB1= hcalMEs[ih]; }
170 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_seq_HB2") ==0 ){ RMS_seq_HB2= hcalMEs[ih]; }
171 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_seq_HE1") ==0 ){ RMS_seq_HE1= hcalMEs[ih]; }
172 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_seq_HE2") ==0 ){ RMS_seq_HE2= hcalMEs[ih]; }
173 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_seq_HE3") ==0 ){ RMS_seq_HE3= hcalMEs[ih]; }
174 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_seq_HF1") ==0 ){ RMS_seq_HF1= hcalMEs[ih]; }
175 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_seq_HF2") ==0 ){ RMS_seq_HF2= hcalMEs[ih]; }
176 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_seq_HO") ==0 ){ RMS_seq_HO= hcalMEs[ih]; }
177 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_map_HB1") ==0 ){ occupancy_map_HB1= hcalMEs[ih]; }
178 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_map_HB2") ==0 ){ occupancy_map_HB2= hcalMEs[ih]; }
179 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_map_HE1") ==0 ){ occupancy_map_HE1= hcalMEs[ih]; }
180 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_map_HE2") ==0 ){ occupancy_map_HE2= hcalMEs[ih]; }
181 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_map_HE3") ==0 ){ occupancy_map_HE3= hcalMEs[ih]; }
182 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_map_HF1") ==0 ){ occupancy_map_HF1= hcalMEs[ih]; }
183 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_map_HF2") ==0 ){ occupancy_map_HF2= hcalMEs[ih]; }
184 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_map_HO") ==0 ){ occupancy_map_HO= hcalMEs[ih]; }
185 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_vs_ieta_HB1") ==0 ){ emean_vs_ieta_HB1= hcalMEs[ih]; }
186 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_vs_ieta_HB2") ==0 ){ emean_vs_ieta_HB2= hcalMEs[ih]; }
187 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_vs_ieta_HE1") ==0 ){ emean_vs_ieta_HE1= hcalMEs[ih]; }
188 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_vs_ieta_HE2") ==0 ){ emean_vs_ieta_HE2= hcalMEs[ih]; }
193 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_vs_ieta_HB1") ==0 ){ RMS_vs_ieta_HB1= hcalMEs[ih]; }
194 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_vs_ieta_HB2") ==0 ){ RMS_vs_ieta_HB2= hcalMEs[ih]; }
195 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_vs_ieta_HE1") ==0 ){ RMS_vs_ieta_HE1= hcalMEs[ih]; }
196 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_vs_ieta_HE2") ==0 ){ RMS_vs_ieta_HE2= hcalMEs[ih]; }
197 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_vs_ieta_HE3") ==0 ){ RMS_vs_ieta_HE3= hcalMEs[ih]; }
198 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_vs_ieta_HF1") ==0 ){ RMS_vs_ieta_HF1= hcalMEs[ih]; }
199 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_vs_ieta_HF2") ==0 ){ RMS_vs_ieta_HF2= hcalMEs[ih]; }
200 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_vs_ieta_HO") ==0 ){ RMS_vs_ieta_HO= hcalMEs[ih]; }
201 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_vs_ieta_HB1") ==0 ){ occupancy_vs_ieta_HB1= hcalMEs[ih]; }
202 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_vs_ieta_HB2") ==0 ){ occupancy_vs_ieta_HB2= hcalMEs[ih]; }
203 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_vs_ieta_HE1") ==0 ){ occupancy_vs_ieta_HE1= hcalMEs[ih]; }
204 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_vs_ieta_HE2") ==0 ){ occupancy_vs_ieta_HE2= hcalMEs[ih]; }
205 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_vs_ieta_HE3") ==0 ){ occupancy_vs_ieta_HE3= hcalMEs[ih]; }
206 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_vs_ieta_HF1") ==0 ){ occupancy_vs_ieta_HF1= hcalMEs[ih]; }
207 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_vs_ieta_HF2") ==0 ){ occupancy_vs_ieta_HF2= hcalMEs[ih]; }
208 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_vs_ieta_HO") ==0 ){ occupancy_vs_ieta_HO= hcalMEs[ih]; }
209 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"HcalRecHitTask_RecHit_StatusWord_HB") ==0 ){ RecHit_StatusWord_HB= hcalMEs[ih]; }
210 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"HcalRecHitTask_RecHit_StatusWord_HE") ==0 ){ RecHit_StatusWord_HE= hcalMEs[ih]; }
211 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"HcalRecHitTask_RecHit_StatusWord_HO") ==0 ){ RecHit_StatusWord_HO= hcalMEs[ih]; }
212 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"HcalRecHitTask_RecHit_StatusWord_HF") ==0 ){ RecHit_StatusWord_HF= hcalMEs[ih]; }
213 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"HcalRecHitTask_RecHit_StatusWord_HF67") ==0 ){ RecHit_StatusWord_HF67= hcalMEs[ih]; }
214 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"HcalRecHitTask_RecHit_Aux_StatusWord_HB") ==0 ){ RecHit_Aux_StatusWord_HB= hcalMEs[ih]; }
215 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"HcalRecHitTask_RecHit_Aux_StatusWord_HE") ==0 ){ RecHit_Aux_StatusWord_HE= hcalMEs[ih]; }
216 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"HcalRecHitTask_RecHit_Aux_StatusWord_HO") ==0 ){ RecHit_Aux_StatusWord_HO= hcalMEs[ih]; }
217 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"HcalRecHitTask_RecHit_Aux_StatusWord_HF") ==0 ){ RecHit_Aux_StatusWord_HF= hcalMEs[ih]; }
219 if( useAllHistos !=0 && useAllHistos !=1 )
return 0;
222 if (subdet==6 && useAllHistos) {
225 double emap_min[82][72][4][4];
238 for (
unsigned int i1 = 0; i1 < 82; i1++) {
239 for (
unsigned int i2 = 0; i2 < 72; i2++) {
241 int index = (i1-41) * 72 + i2;
243 double e = emap_min[i1][i2][0][0];
246 ZS_seqHB1->Fill(
double(index),e);
248 e = emap_min[i1][i2][1][0];
251 ZS_seqHB2->Fill(
double(index),e);
254 e = emap_min[i1][i2][0][1];
257 ZS_seqHE1->Fill(
double(index),e);
259 e = emap_min[i1][i2][1][1];
262 ZS_seqHE2->Fill(
double(index),e);
264 e = emap_min[i1][i2][2][1];
267 ZS_seqHE3->Fill(
double(index),e);
270 e = emap_min[i1][i2][3][2];
273 ZS_seqHO->Fill(
double(index),e);
276 e = emap_min[i1][i2][0][3];
279 ZS_seqHF1->Fill(
double(index),e);
282 e = emap_min[i1][i2][1][3];
285 ZS_seqHF2->Fill(
double(index),e);
288 for (
unsigned int i3 = 0; i3 < 4; i3++) {
289 double emin = 100000.;
290 for (
unsigned int i4 = 0; i4 < 4; i4++) {
297 if ( emin > emap_min [i1][i2][i3][i4])
298 emin = emap_min [i1][i2][i3][i4];
302 if( i3 == 0 && emin < 10000.) {
303 map_depth1->
Fill(
double(ieta),
double(i2),emin);
309 if( i3 == 1 && emin < 10000.)
310 map_depth2->Fill(
double(ieta),
double(i2),emin);
311 if( i3 == 2 && emin < 10000.)
312 map_depth3->Fill(
double(ieta),
double(i2),emin);
313 if( i3 == 3 && emin < 10000.)
314 map_depth4->Fill(
double(ieta),
double(i2),emin);
329 float fev = float (nevtot);
332 float sumphi_hb1, sumphi_hb2, sumphi_he1, sumphi_he2, sumphi_he3,
333 sumphi_ho, sumphi_hf1, sumphi_hf2;
343 for (
int i = 1;
i <= nx1;
i++) {
344 for (
int j = 1;
j <= ny1;
j++) {
347 cnorm = emap_depth2->getBinContent(
i,
j) / fev;
348 emap_depth2->setBinContent(
i,
j,cnorm);
349 cnorm = emap_depth3->getBinContent(
i,
j) / fev;
350 emap_depth3->setBinContent(
i,
j,cnorm);
351 cnorm = emap_depth4->getBinContent(
i,
j) / fev;
352 emap_depth4->setBinContent(
i,
j,cnorm);
357 for (
int i = 1;
i <= nx;
i++) {
367 for (
int j = 1;
j <= ny;
j++) {
369 int index = (
i-42) * ny +
j-1;
375 cnorm = occupancy_map_HB2->getBinContent(
i,
j) / fev;
376 occupancy_map_HB2->setBinContent(
i,
j,cnorm);
381 cnorm = occupancy_map_HE2->getBinContent(
i,
j) / fev;
382 occupancy_map_HE2->setBinContent(
i,
j,cnorm);
384 cnorm = occupancy_map_HE3->getBinContent(
i,
j) / fev;
385 occupancy_map_HE3->setBinContent(
i,
j,cnorm);
393 cnorm = occupancy_map_HF2->getBinContent(
i,
j) / fev;
394 occupancy_map_HF2->setBinContent(
i,
j,cnorm);
397 sumphi_hb2 += occupancy_map_HB2->getBinContent(
i,
j);
399 sumphi_he2 += occupancy_map_HE2->getBinContent(
i,
j);
400 sumphi_he3 += occupancy_map_HE3->getBinContent(
i,
j);
403 sumphi_hf2 += occupancy_map_HF2->getBinContent(
i,
j);
407 occupancy_seqHB1->
Fill(
double(index),cnorm);
408 occupancy_seqHB2->Fill(
double(index),cnorm);
409 occupancy_seqHE1->
Fill(
double(index),cnorm);
410 occupancy_seqHE2->Fill(
double(index),cnorm);
411 occupancy_seqHE3->Fill(
double(index),cnorm);
412 occupancy_seqHO->
Fill(
double(index),cnorm);
413 occupancy_seqHF1->
Fill(
double(index),cnorm);
414 occupancy_seqHF2->Fill(
double(index),cnorm);
419 if(ieta >=0 ) ieta +=1;
421 if(ieta >= -20 && ieta <= 20 )
424 if(ieta >= 40 || ieta <= -40 ) {phi_factor = 18.;}
428 if(ieta >= 0) ieta -= 1;
437 cnorm = sumphi_hb1 / phi_factor;
438 occupancy_vs_ieta_HB1->
Fill(
float(ieta), cnorm);
439 cnorm = sumphi_hb2 / phi_factor;
440 occupancy_vs_ieta_HB2->Fill(
float(ieta), cnorm);
441 cnorm = sumphi_he1 / phi_factor;
442 occupancy_vs_ieta_HE1->
Fill(
float(ieta), cnorm);
443 cnorm = sumphi_he2 / phi_factor;
444 occupancy_vs_ieta_HE2->Fill(
float(ieta), cnorm);
445 cnorm = sumphi_he3 / phi_factor;
446 occupancy_vs_ieta_HE3->Fill(
float(ieta), cnorm);
447 cnorm = sumphi_ho / phi_factor;
448 occupancy_vs_ieta_HO->
Fill(
float(ieta), cnorm);
449 cnorm = sumphi_hf1 / phi_factor;
450 occupancy_vs_ieta_HF1->
Fill(
float(ieta), cnorm);
451 cnorm = sumphi_hf2 / phi_factor;
452 occupancy_vs_ieta_HF2->Fill(
float(ieta), cnorm);
457 RMS_vs_ieta_HB1->
Fill(ieta,cnorm);
458 cnorm = emean_vs_ieta_HB2->getBinError(
i);
459 RMS_vs_ieta_HB2->Fill(ieta,cnorm);
461 RMS_vs_ieta_HE1->
Fill(ieta,cnorm);
462 cnorm = emean_vs_ieta_HE2->getBinError(
i);
463 RMS_vs_ieta_HE2->Fill(ieta,cnorm);
465 RMS_vs_ieta_HE3->Fill(ieta,cnorm);
467 RMS_vs_ieta_HO->
Fill(ieta,cnorm);
469 RMS_vs_ieta_HF1->
Fill(ieta,cnorm);
471 RMS_vs_ieta_HF2->Fill(ieta,cnorm);
479 for(
int ibin = 1; ibin <= nx; ibin++ ){
482 cnorm = emean_seqHB2->getBinError(ibin);
483 RMS_seq_HB2->setBinContent(ibin, cnorm);
488 for(
int ibin = 1; ibin <= nx; ibin++ ){
491 cnorm = emean_seqHE2->getBinError(ibin);
492 RMS_seq_HE2->setBinContent(ibin, cnorm);
493 cnorm = emean_seqHE3->getBinError(ibin);
494 RMS_seq_HE3->setBinContent(ibin, cnorm);
497 for(
int ibin = 1; ibin <= nx; ibin++ ){
500 cnorm = emean_seqHF2->getBinError(ibin);
501 RMS_seq_HF2->setBinContent(ibin, cnorm);
506 for (
int ibin = 1; ibin <= nx; ibin++) {
507 cnorm = RecHit_StatusWord_HB->
getBinContent(ibin) / (fev * 2592.);
510 cnorm = RecHit_StatusWord_HE->getBinContent(ibin) / (fev * 2592.);
511 RecHit_StatusWord_HE->setBinContent(ibin,cnorm);
513 cnorm = RecHit_StatusWord_HO->getBinContent(ibin) / (fev * 2160.);
514 RecHit_StatusWord_HO->setBinContent(ibin,cnorm);
516 cnorm = RecHit_StatusWord_HF->getBinContent(ibin) / (fev * 1728.);
517 RecHit_StatusWord_HF->setBinContent(ibin,cnorm);
519 cnorm = RecHit_Aux_StatusWord_HB->
getBinContent(ibin) / (fev * 2592.);
522 cnorm = RecHit_Aux_StatusWord_HE->getBinContent(ibin) / (fev * 2592.);
523 RecHit_Aux_StatusWord_HE->setBinContent(ibin,cnorm);
525 cnorm = RecHit_Aux_StatusWord_HO->getBinContent(ibin) / (fev * 2160.);
526 RecHit_Aux_StatusWord_HO->setBinContent(ibin,cnorm);
528 cnorm = RecHit_Aux_StatusWord_HF->getBinContent(ibin) / (fev * 1728.);
529 RecHit_Aux_StatusWord_HF->setBinContent(ibin,cnorm);
533 nx = RecHit_StatusWord_HF67->getNbinsX();
534 for (
int ibin = 1; ibin <= nx; ibin++) {
535 cnorm = RecHit_StatusWord_HF67->getBinContent(ibin) / (fev * 1728.);
536 RecHit_StatusWord_HF67->setBinContent(ibin,cnorm);
void setBinContent(int binx, double content)
set content of bin (1-D)
double getEntries(void) const
get # of entries
int getNbinsY(void) const
get # of bins in Y-axis
double getBinError(int binx) const
get uncertainty on content of bin (1-D) - See TH1::GetBinError for details
double getBinContent(int binx) const
get content of bin (1-D)
int getNbinsX(void) const
get # of bins in X-axis