79 int useAllHistos = 0, subdet =5;
86 MonitorElement* ZS_HE1=0, *ZS_seqHE1=0, *ZS_HE2=0, *ZS_seqHE2=0, *ZS_HE3=0, *ZS_seqHE3=0;
87 MonitorElement* map_depth1 =0, *map_depth2 =0, *map_depth3 =0, *map_depth4 =0;
90 MonitorElement* emap_depth1 =0, *emap_depth2 =0, *emap_depth3 =0, *emap_depth4 =0;
92 MonitorElement* occupancy_seqHE1 =0, *occupancy_seqHE2 =0, *occupancy_seqHE3 =0;
96 MonitorElement* emean_seqHE1 =0, *emean_seqHE2 =0, *emean_seqHE3 =0;
108 MonitorElement* occupancy_map_HE1 =0, *occupancy_map_HE2 =0, *occupancy_map_HE3 =0;
115 MonitorElement* RMS_vs_ieta_HE1 =0, *RMS_vs_ieta_HE2 =0, *RMS_vs_ieta_HE3 =0;
118 MonitorElement* occupancy_vs_ieta_HB1 =0, *occupancy_vs_ieta_HB2 =0;
119 MonitorElement* occupancy_vs_ieta_HE1 =0, *occupancy_vs_ieta_HE2 =0, *occupancy_vs_ieta_HE3 =0;
120 MonitorElement* occupancy_vs_ieta_HF1 =0, *occupancy_vs_ieta_HF2 =0;
123 MonitorElement* RecHit_StatusWord_HB =0, *RecHit_StatusWord_HE=0, *RecHit_StatusWord_HO =0, *RecHit_StatusWord_HF =0, *RecHit_StatusWord_HF67 =0;
124 MonitorElement* RecHit_Aux_StatusWord_HB =0, *RecHit_Aux_StatusWord_HE=0, *RecHit_Aux_StatusWord_HO =0, *RecHit_Aux_StatusWord_HF =0;
126 for(
unsigned int ih=0; ih<hcalMEs.size(); ih++){
127 if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_map_depth1") ==0 ){
128 useAllHistos =1; subdet =6;
130 if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_vs_ieta_HB1") ==0 ){
134 if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_simple1D_HO") ==0 ){ ZS_HO = hcalMEs[ih]; }
135 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_sequential1D_HO") ==0 ){ ZS_seqHO = hcalMEs[ih]; }
136 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_simple1D_HB1") ==0 ){ ZS_HB1 = hcalMEs[ih]; }
137 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_sequential1D_HB1") ==0 ){ ZS_seqHB1 = hcalMEs[ih]; }
138 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_simple1D_HB2") ==0 ){ ZS_HB2 = hcalMEs[ih]; }
139 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_sequential1D_HB2") ==0 ){ ZS_seqHB2 = hcalMEs[ih]; }
140 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_simple1D_HF1") ==0 ){ ZS_HF1 = hcalMEs[ih]; }
141 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_sequential1D_HF1") ==0 ){ ZS_seqHF1 = hcalMEs[ih]; }
142 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_simple1D_HF2") ==0 ){ ZS_HF2 = hcalMEs[ih]; }
143 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_sequential1D_HF2") ==0 ){ ZS_seqHF2 = hcalMEs[ih]; }
144 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_simple1D_HE1") ==0 ){ ZS_HE1 = hcalMEs[ih]; }
145 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_sequential1D_HE1") ==0 ){ ZS_seqHE1 = hcalMEs[ih]; }
146 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_simple1D_HE2") ==0 ){ ZS_HE2 = hcalMEs[ih]; }
147 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_sequential1D_HE2") ==0 ){ ZS_seqHE2 = hcalMEs[ih]; }
148 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_simple1D_HE3") ==0 ){ ZS_HE3 = hcalMEs[ih]; }
149 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_sequential1D_HE3") ==0 ){ ZS_seqHE3 = hcalMEs[ih]; }
150 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_map_depth1") ==0 ){ map_depth1= hcalMEs[ih]; }
151 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_map_depth2") ==0 ){ map_depth2= hcalMEs[ih]; }
152 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_map_depth3") ==0 ){ map_depth3= hcalMEs[ih]; }
153 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"ZSmin_map_depth4") ==0 ){ map_depth4= hcalMEs[ih]; }
155 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"N_HF") ==0 ){ Nhf= hcalMEs[ih]; }
156 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emap_depth1") ==0 ){ emap_depth1= hcalMEs[ih]; }
157 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emap_depth2") ==0 ){ emap_depth2= hcalMEs[ih]; }
158 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emap_depth3") ==0 ){ emap_depth3= hcalMEs[ih]; }
159 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emap_depth4") ==0 ){ emap_depth4= hcalMEs[ih]; }
161 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occ_sequential1D_HB1") ==0 ){ occupancy_seqHB1= hcalMEs[ih]; }
162 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occ_sequential1D_HB2") ==0 ){ occupancy_seqHB2= hcalMEs[ih]; }
163 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occ_sequential1D_HE1") ==0 ){ occupancy_seqHE1= hcalMEs[ih]; }
164 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occ_sequential1D_HE2") ==0 ){ occupancy_seqHE2= hcalMEs[ih]; }
165 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occ_sequential1D_HE3") ==0 ){ occupancy_seqHE3= hcalMEs[ih]; }
166 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occ_sequential1D_HF1") ==0 ){ occupancy_seqHF1= hcalMEs[ih]; }
167 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occ_sequential1D_HF2") ==0 ){ occupancy_seqHF2= hcalMEs[ih]; }
168 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occ_sequential1D_HO") ==0 ){ occupancy_seqHO= hcalMEs[ih]; }
169 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_seq_HB1") ==0 ){ emean_seqHB1= hcalMEs[ih]; }
170 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_seq_HB2") ==0 ){ emean_seqHB2= hcalMEs[ih]; }
171 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_seq_HE1") ==0 ){ emean_seqHE1= hcalMEs[ih]; }
172 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_seq_HE2") ==0 ){ emean_seqHE2= hcalMEs[ih]; }
173 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_seq_HE3") ==0 ){ emean_seqHE3= hcalMEs[ih]; }
174 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_seq_HF1") ==0 ){ emean_seqHF1= hcalMEs[ih]; }
175 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_seq_HF2") ==0 ){ emean_seqHF2= hcalMEs[ih]; }
176 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_seq_HO") ==0 ){ emean_seqHO= hcalMEs[ih]; }
177 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_seq_HB1") ==0 ){ RMS_seq_HB1= hcalMEs[ih]; }
178 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_seq_HB2") ==0 ){ RMS_seq_HB2= hcalMEs[ih]; }
179 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_seq_HE1") ==0 ){ RMS_seq_HE1= hcalMEs[ih]; }
180 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_seq_HE2") ==0 ){ RMS_seq_HE2= hcalMEs[ih]; }
181 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_seq_HE3") ==0 ){ RMS_seq_HE3= hcalMEs[ih]; }
182 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_seq_HF1") ==0 ){ RMS_seq_HF1= hcalMEs[ih]; }
183 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_seq_HF2") ==0 ){ RMS_seq_HF2= hcalMEs[ih]; }
184 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_seq_HO") ==0 ){ RMS_seq_HO= hcalMEs[ih]; }
185 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_map_HB1") ==0 ){ occupancy_map_HB1= hcalMEs[ih]; }
186 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_map_HB2") ==0 ){ occupancy_map_HB2= hcalMEs[ih]; }
187 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_map_HE1") ==0 ){ occupancy_map_HE1= hcalMEs[ih]; }
188 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_map_HE2") ==0 ){ occupancy_map_HE2= hcalMEs[ih]; }
189 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_map_HE3") ==0 ){ occupancy_map_HE3= hcalMEs[ih]; }
190 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_map_HF1") ==0 ){ occupancy_map_HF1= hcalMEs[ih]; }
191 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_map_HF2") ==0 ){ occupancy_map_HF2= hcalMEs[ih]; }
192 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_map_HO") ==0 ){ occupancy_map_HO= hcalMEs[ih]; }
193 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_vs_ieta_HB1") ==0 ){ emean_vs_ieta_HB1= hcalMEs[ih]; }
194 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_vs_ieta_HB2") ==0 ){ emean_vs_ieta_HB2= hcalMEs[ih]; }
195 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_vs_ieta_HE1") ==0 ){ emean_vs_ieta_HE1= hcalMEs[ih]; }
196 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"emean_vs_ieta_HE2") ==0 ){ emean_vs_ieta_HE2= hcalMEs[ih]; }
201 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_vs_ieta_HB1") ==0 ){ RMS_vs_ieta_HB1= hcalMEs[ih]; }
202 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_vs_ieta_HB2") ==0 ){ RMS_vs_ieta_HB2= hcalMEs[ih]; }
203 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_vs_ieta_HE1") ==0 ){ RMS_vs_ieta_HE1= hcalMEs[ih]; }
204 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_vs_ieta_HE2") ==0 ){ RMS_vs_ieta_HE2= hcalMEs[ih]; }
205 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_vs_ieta_HE3") ==0 ){ RMS_vs_ieta_HE3= hcalMEs[ih]; }
206 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_vs_ieta_HF1") ==0 ){ RMS_vs_ieta_HF1= hcalMEs[ih]; }
207 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_vs_ieta_HF2") ==0 ){ RMS_vs_ieta_HF2= hcalMEs[ih]; }
208 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"RMS_vs_ieta_HO") ==0 ){ RMS_vs_ieta_HO= hcalMEs[ih]; }
209 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_vs_ieta_HB1") ==0 ){ occupancy_vs_ieta_HB1= hcalMEs[ih]; }
210 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_vs_ieta_HB2") ==0 ){ occupancy_vs_ieta_HB2= hcalMEs[ih]; }
211 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_vs_ieta_HE1") ==0 ){ occupancy_vs_ieta_HE1= hcalMEs[ih]; }
212 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_vs_ieta_HE2") ==0 ){ occupancy_vs_ieta_HE2= hcalMEs[ih]; }
213 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_vs_ieta_HE3") ==0 ){ occupancy_vs_ieta_HE3= hcalMEs[ih]; }
214 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_vs_ieta_HF1") ==0 ){ occupancy_vs_ieta_HF1= hcalMEs[ih]; }
215 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_vs_ieta_HF2") ==0 ){ occupancy_vs_ieta_HF2= hcalMEs[ih]; }
216 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"occupancy_vs_ieta_HO") ==0 ){ occupancy_vs_ieta_HO= hcalMEs[ih]; }
217 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"HcalRecHitTask_RecHit_StatusWord_HB") ==0 ){ RecHit_StatusWord_HB= hcalMEs[ih]; }
218 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"HcalRecHitTask_RecHit_StatusWord_HE") ==0 ){ RecHit_StatusWord_HE= hcalMEs[ih]; }
219 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"HcalRecHitTask_RecHit_StatusWord_HO") ==0 ){ RecHit_StatusWord_HO= hcalMEs[ih]; }
220 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"HcalRecHitTask_RecHit_StatusWord_HF") ==0 ){ RecHit_StatusWord_HF= hcalMEs[ih]; }
221 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"HcalRecHitTask_RecHit_StatusWord_HF67") ==0 ){ RecHit_StatusWord_HF67= hcalMEs[ih]; }
222 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"HcalRecHitTask_RecHit_Aux_StatusWord_HB") ==0 ){ RecHit_Aux_StatusWord_HB= hcalMEs[ih]; }
223 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"HcalRecHitTask_RecHit_Aux_StatusWord_HE") ==0 ){ RecHit_Aux_StatusWord_HE= hcalMEs[ih]; }
224 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"HcalRecHitTask_RecHit_Aux_StatusWord_HO") ==0 ){ RecHit_Aux_StatusWord_HO= hcalMEs[ih]; }
225 else if( strcmp(hcalMEs[ih]->getName().c_str(),
"HcalRecHitTask_RecHit_Aux_StatusWord_HF") ==0 ){ RecHit_Aux_StatusWord_HF= hcalMEs[ih]; }
227 if( useAllHistos !=0 && useAllHistos !=1 )
return 0;
230 if (subdet==6 && useAllHistos) {
233 double emap_min[82][72][4][4];
246 for (
unsigned int i1 = 0; i1 < 82; i1++) {
247 for (
unsigned int i2 = 0; i2 < 72; i2++) {
249 int index = (i1-41) * 72 + i2;
251 double e = emap_min[i1][i2][0][0];
254 ZS_seqHB1->Fill(
double(index),e);
256 e = emap_min[i1][i2][1][0];
259 ZS_seqHB2->Fill(
double(index),e);
262 e = emap_min[i1][i2][0][1];
265 ZS_seqHE1->Fill(
double(index),e);
267 e = emap_min[i1][i2][1][1];
270 ZS_seqHE2->Fill(
double(index),e);
272 e = emap_min[i1][i2][2][1];
275 ZS_seqHE3->Fill(
double(index),e);
278 e = emap_min[i1][i2][3][2];
281 ZS_seqHO->Fill(
double(index),e);
284 e = emap_min[i1][i2][0][3];
287 ZS_seqHF1->Fill(
double(index),e);
290 e = emap_min[i1][i2][1][3];
293 ZS_seqHF2->Fill(
double(index),e);
296 for (
unsigned int i3 = 0; i3 < 4; i3++) {
297 double emin = 100000.;
298 for (
unsigned int i4 = 0; i4 < 4; i4++) {
305 if ( emin > emap_min [i1][i2][i3][i4])
306 emin = emap_min [i1][i2][i3][i4];
310 if( i3 == 0 && emin < 10000.) {
311 map_depth1->
Fill(
double(ieta),
double(i2),emin);
317 if( i3 == 1 && emin < 10000.)
318 map_depth2->Fill(
double(ieta),
double(i2),emin);
319 if( i3 == 2 && emin < 10000.)
320 map_depth3->Fill(
double(ieta),
double(i2),emin);
321 if( i3 == 3 && emin < 10000.)
322 map_depth4->Fill(
double(ieta),
double(i2),emin);
337 float fev = float (nevtot);
340 float sumphi_hb1, sumphi_hb2, sumphi_he1, sumphi_he2, sumphi_he3,
341 sumphi_ho, sumphi_hf1, sumphi_hf2;
351 for (
int i = 1;
i <= nx1;
i++) {
352 for (
int j = 1;
j <= ny1;
j++) {
355 cnorm = emap_depth2->getBinContent(
i,
j) / fev;
356 emap_depth2->setBinContent(
i,
j,cnorm);
357 cnorm = emap_depth3->getBinContent(
i,
j) / fev;
358 emap_depth3->setBinContent(
i,
j,cnorm);
359 cnorm = emap_depth4->getBinContent(
i,
j) / fev;
360 emap_depth4->setBinContent(
i,
j,cnorm);
365 for (
int i = 1;
i <= nx;
i++) {
375 for (
int j = 1;
j <= ny;
j++) {
377 int index = (
i-42) * ny +
j-1;
383 cnorm = occupancy_map_HB2->getBinContent(
i,
j) / fev;
384 occupancy_map_HB2->setBinContent(
i,
j,cnorm);
389 cnorm = occupancy_map_HE2->getBinContent(
i,
j) / fev;
390 occupancy_map_HE2->setBinContent(
i,
j,cnorm);
392 cnorm = occupancy_map_HE3->getBinContent(
i,
j) / fev;
393 occupancy_map_HE3->setBinContent(
i,
j,cnorm);
401 cnorm = occupancy_map_HF2->getBinContent(
i,
j) / fev;
402 occupancy_map_HF2->setBinContent(
i,
j,cnorm);
405 sumphi_hb2 += occupancy_map_HB2->getBinContent(
i,
j);
407 sumphi_he2 += occupancy_map_HE2->getBinContent(
i,
j);
408 sumphi_he3 += occupancy_map_HE3->getBinContent(
i,
j);
411 sumphi_hf2 += occupancy_map_HF2->getBinContent(
i,
j);
415 occupancy_seqHB1->
Fill(
double(index),cnorm);
416 occupancy_seqHB2->Fill(
double(index),cnorm);
417 occupancy_seqHE1->
Fill(
double(index),cnorm);
418 occupancy_seqHE2->Fill(
double(index),cnorm);
419 occupancy_seqHE3->Fill(
double(index),cnorm);
420 occupancy_seqHO->
Fill(
double(index),cnorm);
421 occupancy_seqHF1->
Fill(
double(index),cnorm);
422 occupancy_seqHF2->Fill(
double(index),cnorm);
427 if(ieta >=0 ) ieta +=1;
429 if(ieta >= -20 && ieta <= 20 )
432 if(ieta >= 40 || ieta <= -40 ) {phi_factor = 18.;}
436 if(ieta >= 0) ieta -= 1;
445 cnorm = sumphi_hb1 / phi_factor;
446 occupancy_vs_ieta_HB1->
Fill(
float(ieta), cnorm);
447 cnorm = sumphi_hb2 / phi_factor;
448 occupancy_vs_ieta_HB2->Fill(
float(ieta), cnorm);
449 cnorm = sumphi_he1 / phi_factor;
450 occupancy_vs_ieta_HE1->
Fill(
float(ieta), cnorm);
451 cnorm = sumphi_he2 / phi_factor;
452 occupancy_vs_ieta_HE2->Fill(
float(ieta), cnorm);
453 cnorm = sumphi_he3 / phi_factor;
454 occupancy_vs_ieta_HE3->Fill(
float(ieta), cnorm);
455 cnorm = sumphi_ho / phi_factor;
456 occupancy_vs_ieta_HO->
Fill(
float(ieta), cnorm);
457 cnorm = sumphi_hf1 / phi_factor;
458 occupancy_vs_ieta_HF1->
Fill(
float(ieta), cnorm);
459 cnorm = sumphi_hf2 / phi_factor;
460 occupancy_vs_ieta_HF2->Fill(
float(ieta), cnorm);
465 RMS_vs_ieta_HB1->
Fill(ieta,cnorm);
466 cnorm = emean_vs_ieta_HB2->getBinError(
i);
467 RMS_vs_ieta_HB2->Fill(ieta,cnorm);
469 RMS_vs_ieta_HE1->
Fill(ieta,cnorm);
470 cnorm = emean_vs_ieta_HE2->getBinError(
i);
471 RMS_vs_ieta_HE2->Fill(ieta,cnorm);
473 RMS_vs_ieta_HE3->Fill(ieta,cnorm);
475 RMS_vs_ieta_HO->
Fill(ieta,cnorm);
477 RMS_vs_ieta_HF1->
Fill(ieta,cnorm);
479 RMS_vs_ieta_HF2->Fill(ieta,cnorm);
487 for(
int ibin = 1; ibin <= nx; ibin++ ){
490 cnorm = emean_seqHB2->getBinError(ibin);
491 RMS_seq_HB2->setBinContent(ibin, cnorm);
496 for(
int ibin = 1; ibin <= nx; ibin++ ){
499 cnorm = emean_seqHE2->getBinError(ibin);
500 RMS_seq_HE2->setBinContent(ibin, cnorm);
501 cnorm = emean_seqHE3->getBinError(ibin);
502 RMS_seq_HE3->setBinContent(ibin, cnorm);
505 for(
int ibin = 1; ibin <= nx; ibin++ ){
508 cnorm = emean_seqHF2->getBinError(ibin);
509 RMS_seq_HF2->setBinContent(ibin, cnorm);
514 for (
int ibin = 1; ibin <= nx; ibin++) {
515 cnorm = RecHit_StatusWord_HB->
getBinContent(ibin) / (fev * 2592.);
518 cnorm = RecHit_StatusWord_HE->getBinContent(ibin) / (fev * 2592.);
519 RecHit_StatusWord_HE->setBinContent(ibin,cnorm);
521 cnorm = RecHit_StatusWord_HO->getBinContent(ibin) / (fev * 2160.);
522 RecHit_StatusWord_HO->setBinContent(ibin,cnorm);
524 cnorm = RecHit_StatusWord_HF->getBinContent(ibin) / (fev * 1728.);
525 RecHit_StatusWord_HF->setBinContent(ibin,cnorm);
527 cnorm = RecHit_Aux_StatusWord_HB->
getBinContent(ibin) / (fev * 2592.);
530 cnorm = RecHit_Aux_StatusWord_HE->getBinContent(ibin) / (fev * 2592.);
531 RecHit_Aux_StatusWord_HE->setBinContent(ibin,cnorm);
533 cnorm = RecHit_Aux_StatusWord_HO->getBinContent(ibin) / (fev * 2160.);
534 RecHit_Aux_StatusWord_HO->setBinContent(ibin,cnorm);
536 cnorm = RecHit_Aux_StatusWord_HF->getBinContent(ibin) / (fev * 1728.);
537 RecHit_Aux_StatusWord_HF->setBinContent(ibin,cnorm);
541 nx = RecHit_StatusWord_HF67->getNbinsX();
542 for (
int ibin = 1; ibin <= nx; ibin++) {
543 cnorm = RecHit_StatusWord_HF67->getBinContent(ibin) / (fev * 1728.);
544 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