355 1.89, 5.07, 8.25, 11.43, 14.61, 17.78, 20.96, 24.14, 27.32, 30.50, 33.68,
356 36.86, 40.04, 43.22, 46.40, 49.58, 54.35, 60.71, 67.07, 73.43, 79.79, 86.15,
357 92.51, 98.87, 105.2, 111.6, 117.9, 124.3, 130.7, 137.0, 143.4, 149.7, 156.1,
358 162.5, 168.8, 175.2, 184.7, 197.4, 210.2, 222.9, 235.6, 248.3, 261.0, 273.7,
359 286.5, 299.2, 311.9, 324.6, 337.3, 350.1, 362.8, 375.5, 388.2, 400.9, 413.6,
360 426.4, 439.1, 458.2, 483.6, 509.0, 534.5, 559.9, 585.3, 610.8, 558.9, 584.2,
361 609.5, 634.7, 660.0, 685.3, 710.6, 735.9, 761.2, 786.5, 811.8, 837.1, 862.4,
362 887.7, 913.0, 938.3, 976.2, 1026.8, 1077.4, 1128.0, 1178.6, 1229.1, 1279.7, 1330.3,
363 1380.9, 1431.5, 1482.1, 1532.7, 1583.3, 1633.8, 1684.4, 1735.0, 1785.6, 1836.2, 1886.8,
364 1937.4, 2013.2, 2114.4, 2215.6, 2316.8, 2417.9, 2519.1, 2620.3, 2721.5, 2822.6, 2923.8,
365 3025.0, 3126.2, 3227.3, 3328.5, 3429.7, 3530.9, 3632.0, 3733.2, 3834.4, 3935.5, 4036.7,
366 4188.5, 4390.8, 4593.2, 4795.5, 4997.9, 5200.2, 5402.6, 5057.5, 5262.3, 5467.1, 5671.8,
367 5876.6, 6081.4, 6286.2, 6491.0, 6695.8, 6900.6, 7105.3, 7310.1, 7514.9, 7719.7, 7924.5,
368 8129.3, 8436.4, 8846.0, 9255.6, 9665.1, 10074.7, 10484.3, 10893.9, 11303.4, 11713.0, 12122.6,
369 12532.1, 12941.7, 13351.3, 13760.8, 14170.4, 14580.0, 14989.5, 15399.1, 15808.7, 16218.2, 16832.6,
370 17651.7, 18470.9, 19290.0, 20109.2, 20928.3, 21747.4, 22566.6, 23385.7, 24204.8, 25024.0, 25843.1,
371 26662.3, 27481.4, 28300.5, 29119.7, 29938.8, 30757.9, 31577.1, 32396.2, 33215.4, 34444.1, 36082.3,
372 37720.6, 39358.9, 40997.2, 42635.4, 44273.7, 40908.7, 42542.6, 44176.5, 45810.4, 47444.3, 49078.3,
373 50712.2, 52346.1, 53980.0, 55613.9, 57247.8, 58881.8, 60515.7, 62149.6, 63783.5, 65417.4, 67868.3,
374 71136.1, 74404.0, 77671.8, 80939.7, 84207.5, 87475.3, 90743.2, 94011.0, 97278.8, 100546.7, 103814.5,
375 107082.3, 110350.2, 113618.0, 116885.8, 120153.7, 123421.5, 126689.3, 129957.2, 134858.9, 141394.6, 147930.3,
376 154465.9, 161001.6, 167537.3, 174072.9, 180608.6, 187144.3, 193679.9, 200215.6, 206751.3, 213287.0, 219822.6,
377 226358.3, 232894.0, 239429.6, 245965.3, 252501.0, 259036.6, 265572.3, 275375.8, 288447.2, 301518.5, 314589.8,
378 327661.2, 340732.5, 353803.8};
381 9.56, 28.24, 46.91, 65.59, 84.27, 102.9, 121.6, 140.3, 159.0, 177.7,
382 196.3, 215.0, 233.7, 252.4, 271.0, 289.7, 317.7, 355.1, 392.4, 429.8,
383 467.1, 504.5, 541.9, 579.2, 616.6, 653.9, 691.3, 728.6, 766.0, 803.3,
384 840.7, 878.0, 915.4, 952.8, 990.1, 1027.5, 1083.5, 1158.2, 1232.9, 1307.6,
385 1382.3, 1457.0, 1531.7, 1606.4, 1681.2, 1755.9, 1830.6, 1905.3, 1980.0, 2054.7,
386 2129.4, 2204.1, 2278.8, 2353.5, 2428.2, 2502.9, 2577.7, 2689.7, 2839.1, 2988.6,
387 3138.0, 3287.4, 3436.8, 3586.2, 3263.0, 3411.3, 3559.6, 3707.9, 3856.2, 4004.5,
388 4152.9, 4301.2, 4449.5, 4597.8, 4746.1, 4894.4, 5042.7, 5191.0, 5339.4, 5487.7,
389 5710.1, 6006.8, 6303.4, 6600.0, 6896.6, 7193.3, 7489.9, 7786.5, 8083.1, 8379.8,
390 8676.4, 8973.0, 9269.6, 9566.3, 9862.9, 10159.5, 10456.2, 10752.8, 11049.4, 11346.0,
391 11791.0, 12384.2, 12977.5, 13570.7, 14164.0, 14757.2, 15350.5, 15943.7, 16537.0, 17130.2,
392 17723.5, 18316.7, 18910.0, 19503.2, 20096.5, 20689.7, 21283.0, 21876.2, 22469.5, 23062.8,
393 23656.0, 24545.9, 25732.4, 26918.9, 28105.4, 29291.9, 30478.4, 31664.9, 29399.4, 30590.1,
394 31780.9, 32971.7, 34162.4, 35353.2, 36544.0, 37734.7, 38925.5, 40116.3, 41307.0, 42497.8,
395 43688.5, 44879.3, 46070.1, 47260.8, 49047.0, 51428.5, 53810.1, 56191.6, 58573.1, 60954.6,
396 63336.2, 65717.7, 68099.2, 70480.8, 72862.3, 75243.8, 77625.4, 80006.9, 82388.4, 84769.9,
397 87151.5, 89533.0, 91914.5, 94296.1, 97868.4, 102631.4, 107394.5, 112157.5, 116920.6, 121683.7,
398 126446.7, 131209.8, 135972.8, 140735.9, 145499.0, 150262.0, 155025.1, 159788.2, 164551.2, 169314.3,
399 174077.3, 178840.4, 183603.5, 188366.5, 193129.6, 200274.2, 209800.3, 219326.4, 228852.5, 238378.7,
400 247904.8, 257430.9, 237822.7, 247326.7, 256830.7, 266334.8, 275838.8, 285342.9, 294846.9, 304351.0,
401 313855.0, 323359.1, 332863.1, 342367.1, 351871.2, 361375.2, 370879.3, 380383.3, 394639.4, 413647.5,
402 432655.6, 451663.6, 470671.7, 489679.8, 508687.9, 527696.0, 546704.1, 565712.2, 584720.3, 603728.3,
403 622736.4, 641744.5, 660752.6, 679760.7, 698768.8, 717776.9, 736785.0, 755793.0, 784305.2, 822321.3,
404 860337.5, 898353.7, 936369.9, 974386.0, 1012402.2, 1050418.4, 1088434.6, 1126450.7, 1164466.9, 1202483.1,
405 1240499.3, 1278515.4, 1316531.6, 1354547.8, 1392564.0, 1430580.1, 1468596.3, 1506612.5, 1544628.7, 1601652.9,
406 1677685.3, 1753717.6, 1829750.0, 1905782.3, 1981814.7, 2057847.0};
409 -0.5, 0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5, 10.5, 11.5,
410 12.5, 13.5, 15., 17., 19., 21., 23., 25., 27., 29.5, 32.5, 35.5, 38.5,
411 42., 46., 50., 54.5, 59.5, 64.5, 59.5, 64.5, 69.5, 74.5, 79.5, 84.5, 89.5,
412 94.5, 99.5, 104.5, 109.5, 114.5, 119.5, 124.5, 129.5, 137., 147., 157., 167., 177.,
413 187., 197., 209.5, 224.5, 239.5, 254.5, 272., 292., 312., 334.5, 359.5, 384.5, 359.5,
414 384.5, 409.5, 434.5, 459.5, 484.5, 509.5, 534.5, 559.5, 584.5, 609.5, 634.5, 659.5, 684.5,
415 709.5, 747., 797., 847., 897., 947., 997., 1047., 1109.5, 1184.5, 1259.5, 1334.5, 1422.,
416 1522., 1622., 1734.5, 1859.5, 1984.5, 1859.5, 1984.5, 2109.5, 2234.5, 2359.5, 2484.5, 2609.5, 2734.5,
417 2859.5, 2984.5, 3109.5, 3234.5, 3359.5, 3484.5, 3609.5, 3797., 4047., 4297., 4547., 4797., 5047.,
418 5297., 5609.5, 5984.5, 6359.5, 6734.5, 7172., 7672., 8172., 8734.5, 9359.5, 9984.5};
442 double phi12(
double phi1,
double en1,
double phi2,
double en2);
443 double dR(
double eta1,
double phi1,
double eta2,
double phi2);
444 double dPhiWsign(
double phi1,
double phi2);
2375 std::cout <<
"============== endJob ===================================" << std::endl;
2377 std::cout <<
" --------------------------------------- " << std::endl;
2380 std::cout <<
" --------------------------------------------- " << std::endl;
2386 std::cout <<
"===== possible max number of events in ntuple(each 50th recorded) = " <<
nevent50 << std::endl;
2390 std::cout <<
"===== full number of events*HBHEdigis = " <<
nnnnnn << std::endl;
2395 std::cout <<
"===== Start writing user histograms =====" << std::endl;
2629 for (
int k1 = 0; k1 <
ndepth; k1++) {
2635 for (
int k3 = 0; k3 <
nphi; k3++) {
2694 int outabortgap = 1;
2704 std::cout <<
" --------------------------------------- " << std::endl;
2708 std::cout <<
" --------------------------------------------- " << std::endl;
2714 std::cout <<
" ---------***********************------------- " << std::endl;
2764 unsigned long int pcountall1 = 0;
2765 unsigned long int pcountall3 = 0;
2766 unsigned long int pcountall6 = 0;
2767 unsigned long int pcountall8 = 0;
2772 unsigned long int mcountall1 = 0;
2773 unsigned long int mcountall3 = 0;
2774 unsigned long int mcountall6 = 0;
2775 unsigned long int mcountall8 = 0;
2784 int pnnmin1 = 999999999;
2785 int pnnmin3 = 999999999;
2786 int pnnmin6 = 999999999;
2787 int pnnmin8 = 999999999;
2792 int mnnmin1 = 999999999;
2793 int mnnmin3 = 999999999;
2794 int mnnmin6 = 999999999;
2795 int mnnmin8 = 999999999;
2797 for (
int k1 = 0; k1 <
ndepth; k1++) {
2798 for (
int k3 = 0; k3 <
nphi; k3++) {
2908 int kkkk2 = (
k2 - 1) / 4;
2913 int kkkk3 = (k3) / 4 + 1;
2917 ietaphi = ((kkkk2)-1) *
znphi + (kkkk3);
2985 if (bbbc / bbb1 > 25.) {
3000 if (bbbc / bbb1 > 25.) {
3085 if (bbbc / bbb1 > 15. && k3 % 2 == 0) {
3099 if (bbbc / bbb1 > 15. && k3 % 2 == 0) {
3239 if (bbbc / bbb1 > 20.) {
3254 if (bbbc / bbb1 > 20.) {
3346 if (bbbc / bbb1 > 80.) {
3361 if (bbbc / bbb1 > 80.) {
3916 if (
k0 == 0 && k1 == 0) {
3917 if (pcountmin1 > 0) {
3918 if (pcountmin1 < pnnmin1)
3919 pnnmin1 = pcountmin1;
3923 if (mcountmin1 > 0) {
3924 if (mcountmin1 < mnnmin1)
3925 mnnmin1 = mcountmin1;
3931 if (
k0 == 1 && k1 == 0) {
3932 if (pcountmin3 > 0) {
3933 if (pcountmin3 < pnnmin3)
3934 pnnmin3 = pcountmin3;
3938 if (mcountmin3 > 0) {
3939 if (mcountmin3 < mnnmin3)
3940 mnnmin3 = mcountmin3;
3946 if (
k0 == 2 && k1 == 3) {
3947 if (pcountmin8 > 0) {
3948 if (pcountmin8 < pnnmin8)
3949 pnnmin8 = pcountmin8;
3953 if (mcountmin8 > 0) {
3954 if (mcountmin8 < mnnmin8)
3955 mnnmin8 = mcountmin8;
3961 if (
k0 == 3 && k1 == 0) {
3962 if (pcountmin6 > 0) {
3963 if (pcountmin6 < pnnmin6)
3964 pnnmin6 = pcountmin6;
3968 if (mcountmin6 > 0) {
3969 if (mcountmin6 < mnnmin6)
3970 mnnmin6 = mcountmin6;
3985 float patiooccupancy1 = 0.;
3986 if (pcountall1 != 0)
3987 patiooccupancy1 = (float)pnnmin1 * mnnbins1 / pcountall1;
3989 float matiooccupancy1 = 0.;
3990 if (mcountall1 != 0)
3991 matiooccupancy1 = (float)mnnmin1 * mnnbins1 / mcountall1;
3994 float patiooccupancy3 = 0.;
3995 if (pcountall3 != 0)
3996 patiooccupancy3 = (float)pnnmin3 * mnnbins3 / pcountall3;
3998 float matiooccupancy3 = 0.;
3999 if (mcountall3 != 0)
4000 matiooccupancy3 = (float)mnnmin3 * mnnbins3 / mcountall3;
4003 float patiooccupancy6 = 0.;
4004 if (pcountall6 != 0)
4005 patiooccupancy6 = (float)pnnmin6 * mnnbins6 / pcountall6;
4007 float matiooccupancy6 = 0.;
4008 if (mcountall6 != 0)
4009 matiooccupancy6 = (float)mnnmin6 * mnnbins6 / mcountall6;
4012 float patiooccupancy8 = 0.;
4013 if (pcountall8 != 0)
4014 patiooccupancy8 = (float)pnnmin8 * mnnbins8 / pcountall8;
4016 float matiooccupancy8 = 0.;
4017 if (mcountall8 != 0)
4018 matiooccupancy8 = (float)mnnmin8 * mnnbins8 / mcountall8;
4022 for (
int k1 = 0; k1 <
ndepth; k1++) {
4024 for (
int k3 = 0; k3 <
nphi; k3++) {
4124 for (
int k1 = 0; k1 <
ndepth; k1++) {
4126 for (
int k3 = 0; k3 <
nphi; k3++) {
4148 for (
int k1 = 0; k1 <
ndepth; k1++) {
4150 for (
int k3 = 0; k3 <
nphi; k3++) {
4179 bool gotHFDigis =
true;
4187 std::cout <<
" ****************************** =========================== No HFDigiCollection found "
4192 eta = digi->id().ieta();
4193 phi = digi->id().iphi();
4194 depth = digi->id().depth();
4211 for (
int i = 0;
i <
nTS;
i++) {
4230 bool gotQIE10Digis =
true;
4232 gotQIE10Digis =
false;
4235 gotQIE10Digis =
false;
4237 if (!gotQIE10Digis) {
4238 std::cout <<
" No QIE10DigiCollection collection is found " << std::endl;
4241 double totalAmplitudeHF = 0.;
4242 for (
unsigned int j = 0;
j < qie10dc.
size();
j++) {
4266 double amplitudefullTSs = 0.;
4267 double nnnnnnTS = 0.;
4268 for (
int i = 0;
i <
nTS; ++
i) {
4271 int adc = qie10df[
i].adc();
4303 bool gotHBHEDigis =
true;
4305 gotHBHEDigis =
false;
4307 gotHBHEDigis =
false;
4308 if (!gotHBHEDigis) {
4309 std::cout <<
" No HBHEDigiCollection collection is found " << std::endl;
4312 double totalAmplitudeHB = 0.;
4313 double totalAmplitudeHE = 0.;
4314 double nnnnnnTSHB = 0.;
4315 double nnnnnnTSHE = 0.;
4318 eta = digi->id().ieta();
4319 phi = digi->id().iphi();
4320 depth = digi->id().depth();
4330 if (digi->id().subdet() ==
HcalBarrel && depth == 2 && qwert2 == 0) {
4334 if (digi->id().subdet() ==
HcalEndcap && depth == 1 && qwert3 == 0) {
4338 if (digi->id().subdet() ==
HcalEndcap && depth == 2 && qwert4 == 0) {
4342 if (digi->id().subdet() ==
HcalEndcap && depth == 3 && qwert5 == 0) {
4360 double amplitudefullTSs = 0.;
4363 for (
int i = 0;
i <
nTS;
i++) {
4375 double amplitudefullTSs = 0.;
4378 for (
int i = 0;
i <
nTS;
i++) {
4393 if (totalAmplitudeHB != 0.) {
4398 if (totalAmplitudeHE != 0.) {
4413 bool gotQIE11Digis =
true;
4415 gotQIE11Digis =
false;
4417 gotQIE11Digis =
false;
4418 if (!gotQIE11Digis) {
4419 std::cout <<
" No QIE11DigiCollection collection is found " << std::endl;
4422 double totalAmplitudeHBQIE11 = 0.;
4423 double totalAmplitudeHEQIE11 = 0.;
4424 double nnnnnnTSHBQIE11 = 0.;
4425 double nnnnnnTSHEQIE11 = 0.;
4426 for (
unsigned int j = 0;
j < qie11dc.
size();
j++) {
4435 int sub = hcaldetid.
subdet();
4448 if (sub == 1 && depth == 1 && qwert1 == 0) {
4452 if (sub == 1 && depth == 2 && qwert2 == 0) {
4456 if (sub == 2 && depth == 1 && qwert3 == 0) {
4460 if (sub == 2 && depth == 2 && qwert4 == 0) {
4464 if (sub == 2 && depth == 3 && qwert5 == 0) {
4476 double amplitudefullTSs1 = 0.;
4477 double amplitudefullTSs6 = 0.;
4479 for (
int i = 0;
i <
nTS; ++
i) {
4480 int adc = qie11df[
i].adc();
4483 amplitudefullTSs1 += charge1;
4484 amplitudefullTSs6 += charge6;
4490 totalAmplitudeHBQIE11 +=
charge;
4497 double amplitudefullTSs1 = 0.;
4498 double amplitudefullTSs6 = 0.;
4500 for (
int i = 0;
i <
nTS;
i++) {
4501 int adc = qie11df[
i].adc();
4504 amplitudefullTSs1 += charge1;
4505 amplitudefullTSs6 += charge6;
4511 totalAmplitudeHEQIE11 +=
charge;
4520 if (totalAmplitudeHBQIE11 != 0.) {
4525 if (totalAmplitudeHEQIE11 != 0.) {
4536 bool gotHODigis =
true;
4543 std::cout <<
" No HO collection is found " << std::endl;
4546 double totalAmplitudeHO = 0.;
4548 eta = digi->id().ieta();
4549 phi = digi->id().iphi();
4550 depth = digi->id().depth();
4571 double nnnnnnTS = 0.;
4572 double amplitudefullTSs = 0.;
4574 for (
int i = 0;
i <
nTS;
i++) {
4602 bool gotHBHERecHitsNoise =
true;
4604 gotHBHERecHitsNoise =
4607 gotHBHERecHitsNoise =
false;
4608 if (!gotHBHERecHitsNoise) {
4610 std::cout <<
" No RecHits HBHENoise collection is found " << std::endl;
4615 float icalconst = 1.;
4618 HBHERecHit aHit(hbheItr->id(), hbheItr->eraw() * icalconst, hbheItr->time());
4620 double energyhit = aHit.
energy();
4621 DetId id = (*hbheItr).detid();
4623 int sub = ((hid).rawId() >> 25) & 0x7;
4632 if (hid.
depth() == 1 && sub == 1 && hid.
iphi() == 25) {
4634 std::cout <<
" Noise,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4635 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4637 if (hid.
depth() == 1 && sub == 2 && hid.
iphi() == 25) {
4639 std::cout <<
" Noise,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4640 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4642 int ieta = hid.
ieta();
4645 int iphi = hid.
iphi() - 1;
4646 int mdepth = hid.
depth();
4657 bool gotHFRecHitsNoise =
true;
4659 gotHFRecHitsNoise =
false;
4661 gotHFRecHitsNoise =
false;
4662 if (!gotHFRecHitsNoise) {
4664 std::cout <<
" No RecHits HFNoise collection is found " << std::endl;
4668 float icalconst = 1.;
4671 HFRecHit aHit(hfItr->id(), hfItr->energy() * icalconst, hfItr->time());
4672 double energyhit = aHit.
energy();
4673 DetId id = (*hfItr).detid();
4675 int sub = ((hid).rawId() >> 25) & 0x7;
4681 if (hid.
iphi() == 25) {
4683 std::cout <<
"HF Noise,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4684 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4686 int ieta = hid.
ieta();
4689 int iphi = hid.
iphi() - 1;
4690 int mdepth = hid.
depth();
4704 bool gotHBHERecHitsSignal =
true;
4706 gotHBHERecHitsSignal =
4709 gotHBHERecHitsSignal =
false;
4710 if (!gotHBHERecHitsSignal) {
4712 std::cout <<
" No RecHits HBHESignal collection is found " << std::endl;
4717 float icalconst = 1.;
4720 HBHERecHit aHit(hbheItr->id(), hbheItr->eraw() * icalconst, hbheItr->time());
4722 double energyhit = aHit.
energy();
4723 DetId id = (*hbheItr).detid();
4725 int sub = ((hid).rawId() >> 25) & 0x7;
4733 if (hid.
depth() == 1 && sub == 1 && hid.
iphi() == 25) {
4735 std::cout <<
"HBHE Signal,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4736 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4738 if (hid.
depth() == 1 && sub == 2 && hid.
iphi() == 25) {
4740 std::cout <<
"HBHE Signal,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4741 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4743 int ieta = hid.
ieta();
4746 int iphi = hid.
iphi() - 1;
4747 int mdepth = hid.
depth();
4758 bool gotHFRecHitsSignal =
true;
4760 gotHFRecHitsSignal =
false;
4762 gotHFRecHitsSignal =
false;
4763 if (!gotHFRecHitsSignal) {
4765 std::cout <<
" No RecHits HFSignal collection is found " << std::endl;
4769 float icalconst = 1.;
4772 HFRecHit aHit(hfItr->id(), hfItr->energy() * icalconst, hfItr->time());
4773 double energyhit = aHit.
energy();
4774 DetId id = (*hfItr).detid();
4776 int sub = ((hid).rawId() >> 25) & 0x7;
4782 if (hid.
iphi() == 25) {
4784 std::cout <<
"HF Signal,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4785 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4787 int ieta = hid.
ieta();
4790 int iphi = hid.
iphi() - 1;
4791 int mdepth = hid.
depth();
4821 for (
int k1 = 0; k1 <
ndepth; k1++) {
4825 int k2plot =
k2 - 41;
4828 double sumoverphi = 0;
4829 int nsumoverphi = 0;
4830 for (
int k3 = 0; k3 <
nphi; k3++) {
4835 std::cout <<
"==== nsumoverphi = " << nsumoverphi <<
" sumoverphi = " << sumoverphi
4836 <<
" k1 = " << k1 <<
" k2 = " <<
k2 <<
" kkk = " << kkk <<
" k3 = " << k3 << std::endl;
4840 for (
int k3 = 0; k3 <
nphi; k3++) {
4841 if (nsumoverphi != 0) {
4844 std::cout <<
"nsumoverphi= " << nsumoverphi <<
" sumoverphi= " << sumoverphi <<
" k1= " << k1
4845 <<
" k2= " <<
k2 <<
" kkk= " << kkk <<
" k3= " << k3
4921 for (
int k1 = 0; k1 <
ndepth; k1++) {
4925 int k2plot =
k2 - 41;
4927 for (
int k3 = 0; k3 <
nphi; k3++) {
4955 for (
int k1 = 0; k1 <
ndepth; k1++) {
4959 int k2plot =
k2 - 41;
4961 for (
int k3 = 0; k3 <
nphi; k3++) {
5004 for (
int k1 = 0; k1 <
ndepth; k1++) {
5008 int k2plot =
k2 - 41;
5010 for (
int k3 = 0; k3 <
nphi; k3++) {
5053 for (
int k1 = 0; k1 <
ndepth; k1++) {
5057 int k2plot =
k2 - 41;
5059 for (
int k3 = 0; k3 <
nphi; k3++) {
5099 for (
int k1 = 0; k1 <
ndepth; k1++) {
5103 int k2plot =
k2 - 41;
5105 for (
int k3 = 0; k3 <
nphi; k3++) {
5169 for (
int k1 = 0; k1 <
ndepth; k1++) {
5173 int k2plot =
k2 - 41;
5175 for (
int k3 = 0; k3 <
nphi; k3++) {
5223 for (
int k1 = 0; k1 <
ndepth; k1++) {
5225 for (
int k3 = 0; k3 <
nphi; k3++) {
5237 for (
int k1 = 0; k1 <
ndepth; k1++) {
5241 if (k1 == 2 &&
k2 == 25) {
5245 }
else if (k1 == 2 &&
k2 == 56) {
5260 for (
int k1 = 0; k1 <
ndepth; k1++) {
5262 if (k1 == 2 && (
k2 == 25 ||
k2 == 56)) {
5265 int k2plot =
k2 - 41;
5275 int LLLLLL111111 = 0;
5276 if ((k1 == 0 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 30))
5279 int LLLLLL222222 = 0;
5280 if ((k1 == 1 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 27) ||
5281 (k1 == 2 && fabs(kkk - 41) > 26 && fabs(kkk - 41) < 29))
5284 if (LLLLLL111111 == 1) {
5288 if (LLLLLL222222 == 1) {
5301 for (
int k1 = 0; k1 <
ndepth; k1++) {
5303 if (k1 == 2 && (
k2 == 25 ||
k2 == 56)) {
5305 int k2plot =
k2 - 41;
5312 for (
int k3 = 0; k3 <
nphi; k3++) {
5315 int LLLLLL111111 = 0;
5316 if ((k1 == 0 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 30))
5318 int LLLLLL222222 = 0;
5319 if ((k1 == 1 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 27) ||
5320 (k1 == 2 && fabs(kkk - 41) > 26 && fabs(kkk - 41) < 29))
5322 if (LLLLLL111111 == 1) {
5326 if (LLLLLL222222 == 1) {
5337 for (
int k1 = 0; k1 <
ndepth; k1++) {
5341 for (
int k3 = 0; k3 <
nphi; k3++) {
5359 int sumofchannels = 0;
5360 double sumamplitudesubdet = 0.;
5361 int sumofchannels0 = 0;
5362 double sumamplitudesubdet0 = 0.;
5364 for (
int k3 = 0; k3 <
nphi; k3++) {
5376 if (sumamplitudechannel_HB > 80.) {
5377 sumamplitudesubdet += sumamplitudechannel_HB;
5380 if (sumamplitudechannel_HB > 0.) {
5381 sumamplitudesubdet0 += sumamplitudechannel_HB;
5389 double sumamplitudechannel_HE =
5392 if (sumamplitudechannel_HE > 200.) {
5393 sumamplitudesubdet += sumamplitudechannel_HE;
5396 if (sumamplitudechannel_HE > 0.) {
5397 sumamplitudesubdet0 += sumamplitudechannel_HE;
5407 if (sumamplitudechannel_HO > 1200.) {
5408 sumamplitudesubdet += sumamplitudechannel_HO;
5411 if (sumamplitudechannel_HO > 0.) {
5412 sumamplitudesubdet0 += sumamplitudechannel_HO;
5421 if (sumamplitudechannel_HF > 600.) {
5422 sumamplitudesubdet += sumamplitudechannel_HF;
5425 if (sumamplitudechannel_HF > 0.) {
5426 sumamplitudesubdet0 += sumamplitudechannel_HF;
5439 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM1)
5440 maxxSUM1 = sumamplitudesubdet + sumamplitudesubdet0;
5441 if ((sumofchannels + sumofchannels0) >
maxxOCCUP1)
5447 if ((sumamplitudesubdet + sumamplitudesubdet0) > 60000) {
5449 for (
int k3 = 0; k3 <
nphi; k3++) {
5468 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM2)
5469 maxxSUM2 = sumamplitudesubdet + sumamplitudesubdet0;
5470 if ((sumofchannels + sumofchannels0) >
maxxOCCUP2)
5480 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM3)
5481 maxxSUM3 = sumamplitudesubdet + sumamplitudesubdet0;
5482 if ((sumofchannels + sumofchannels0) >
maxxOCCUP3)
5492 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM4)
5493 maxxSUM4 = sumamplitudesubdet + sumamplitudesubdet0;
5494 if ((sumofchannels + sumofchannels0) >
maxxOCCUP4)
5515 for (
int k1 = 0; k1 <
ndepth; k1++) {
5517 int nbadchannels = 0;
5519 for (
int k3 = 0; k3 <
nphi; k3++) {
5532 if (nbadchannels != 0)
5545 if (nbadchannels != 0)
5561 if (nbadchannels != 0)
5574 if (nbadchannels != 0)
5587 if (nbadchannels != 0)
5603 if (nbadchannels != 0)
5619 if (nbadchannels != 0)
5632 if (nbadchannels != 0)
5652 bool gotCALIBDigis =
true;
5654 gotCALIBDigis =
false;
5657 gotCALIBDigis =
false;
5659 if (!gotCALIBDigis) {
5662 int cal_det = digi->id().hcalSubdet();
5663 int cal_phi = digi->id().iphi();
5664 int cal_eta = digi->id().ieta();
5665 int cal_cbox = digi->id().cboxChannel();
5669 if (cal_det > 0 && cal_det < 5 && cal_cbox == 0) {
5670 int iphi = cal_phi - 1;
5675 double max_signal = -100.;
5676 int ts_with_max_signal = -100;
5681 for (
int i = 0;
i <
nTS;
i++) {
5682 double ampldefault =
adc2fC[digi->sample(
i).adc() & 0xff];
5683 if (max_signal < ampldefault) {
5684 max_signal = ampldefault;
5685 ts_with_max_signal =
i;
5688 calib3[cal_det - 1][ieta + 41][iphi] += ampldefault;
5689 calib0[cal_det - 1][ieta + 41][iphi] += ampldefault;
5690 timew += (
i + 1) * ampldefault;
5694 double amplitude =
calib0[cal_det - 1][ieta + 41][iphi];
5695 double aveamplitude = -100.;
5696 if (amplitude > 0 && timew > 0)
5697 aveamplitude = timew / amplitude;
5698 double aveamplitude1 = aveamplitude - 1;
5699 caliba[cal_det - 1][ieta + 41][iphi] = aveamplitude1;
5703 double ampldefault =
adc2fC[digi->sample(
ii).adc() & 0xff];
5704 double aaaaaa = (
ii + 1) - aveamplitude;
5705 double aaaaaa2 = aaaaaa * aaaaaa;
5706 rmsamp += (aaaaaa2 * ampldefault);
5708 double rmsamplitude = -100.;
5709 if ((amplitude > 0 && rmsamp > 0) || (amplitude < 0 && rmsamp < 0))
5710 rmsamplitude =
sqrt(rmsamp / amplitude);
5711 calibw[cal_det - 1][ieta + 41][iphi] = rmsamplitude;
5713 calibt[cal_det - 1][ieta + 41][iphi] = ts_with_max_signal;
5716 if (ts_with_max_signal > -1 && ts_with_max_signal < nTS)
5717 calib2[cal_det - 1][ieta + 41][iphi] =
adc2fC[digi->sample(ts_with_max_signal).adc() & 0xff];
5718 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < nTS)
5719 calib2[cal_det - 1][ieta + 41][iphi] +=
adc2fC[digi->sample(ts_with_max_signal + 1).adc() & 0xff];
5720 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < nTS)
5721 calib2[cal_det - 1][ieta + 41][iphi] +=
adc2fC[digi->sample(ts_with_max_signal + 2).adc() & 0xff];
5722 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < nTS)
5723 calib2[cal_det - 1][ieta + 41][iphi] +=
adc2fC[digi->sample(ts_with_max_signal - 1).adc() & 0xff];
5724 if (ts_with_max_signal - 2 > -1 && ts_with_max_signal - 2 < nTS)
5725 calib2[cal_det - 1][ieta + 41][iphi] +=
adc2fC[digi->sample(ts_with_max_signal - 2).adc() & 0xff];
5727 bool anycapid =
true;
5730 int error1 = 0, error2 = 0, error3 = 0;
5731 int lastcapid = 0, capid = 0;
5732 for (
int ii = 0;
ii < (*digi).size();
ii++) {
5733 capid = (*digi)[
ii].capid();
5734 bool er = (*digi)[
ii].er();
5735 bool dv = (*digi)[
ii].dv();
5736 if (
ii != 0 && ((lastcapid + 1) % 4) != capid)
5750 if (error1 != 0 || error2 != 0 || error3 != 0)
5765 for (
int k1 = 0; k1 <
nsub; k1++) {
5767 for (
int k3 = 0; k3 <
nphi; k3++) {
5768 int k2plot =
k2 - 41;
5798 if (k20 > 0 || k20 == 0)
5803 int kk2 = 0, kk3 = 0;
5804 if (k1 == 0 || k1 == 1) {
5809 if (k30 == 71 || k30 == nphi || k30 == 1 || k30 == 2)
5812 kk3 = ((k30 - 3) / 4) * 4 + 3;
5813 }
else if (k1 == 2) {
5814 if (
abs(k20) <= 4) {
5816 if (k30 == 71 || k30 == nphi || k30 == 1 || k30 == 2 || k30 == 3 || k30 == 4)
5819 kk3 = ((k30 - 5) / 6) * 6 + 5;
5821 if (
abs(k20) > 4 &&
abs(k20) <= 10)
5823 if (
abs(k20) > 10 &&
abs(k20) <= 15)
5827 if (k30 == 71 || k30 == nphi || (k30 >= 1 && k30 <= 10))
5830 kk3 = ((k30 - 11) / 12) * 12 + 11;
5832 }
else if (k1 == 3) {
5837 if (k30 >= 1 && k30 <= 18)
5839 if (k30 >= 19 && k30 <= 36)
5841 if (k30 >= 37 && k30 <= 54)
5843 if (k30 >= 55 && k30 <= nphi)
5847 int kkk2 = kk2 + 41;
5854 double GetRMSOverNormalizedSignal = -1.;
5856 GetRMSOverNormalizedSignal =
signal[k1][
k2][k3] /
calib0[k1][kkk2][kkk3];
5879 if (
calib0[k1][kkk2][kkk3] > 0.)
5880 adc =
calib0[k1][kkk2][kkk3];
5883 if (
calib0[k1][kkk2][kkk3] > 0.)
5884 ratio =
calib2[k1][kkk2][kkk3] /
calib0[k1][kkk2][kkk3];
5886 float calibtsmax =
calibt[k1][kkk2][kkk3];
5888 float calibtsmean =
caliba[k1][kkk2][kkk3];
5890 float calibwidth =
calibw[k1][kkk2][kkk3];
5892 float calibcap = -100.;
5909 if (calibtsmax > -0.5) {
5916 if (calibtsmean > -0.5) {
5923 if (calibwidth > -0.5) {
5933 h_map_HB->Fill(
double(k2plot),
double(k3), 1.);
5949 if (calibtsmax > -0.5) {
5956 if (calibtsmean > -0.5) {
5963 if (calibwidth > -0.5) {
5973 h_map_HE->Fill(
double(k2plot),
double(k3), 1.);
5989 if (calibtsmax > -0.5) {
5996 if (calibtsmean > -0.5) {
6003 if (calibwidth > -0.5) {
6013 h_map_HO->Fill(
double(k2plot),
double(k3), 1.);
6029 if (calibtsmax > -0.5) {
6036 if (calibtsmean > -0.5) {
6043 if (calibwidth > -0.5) {
6053 h_map_HF->Fill(
double(k2plot),
double(k3), 1.);
6085 std::cout <<
"======================== beignJob START +++++++++++++++++++++++++++" << std::endl;
6323 fs_->
make<TH1F>(
"h_ndepthesperamplitudebins_HE",
" ", 10, 0., 10.);
6333 "h_mapADCAmplsecondpeak0_HE",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
6690 fs_->
make<TH1F>(
"h_AamplitudewithPedSubtr_RADDAM_HEzoom0",
" ", 100, min80, 4000.);
6692 fs_->
make<TH1F>(
"h_AamplitudewithPedSubtr_RADDAM_HEzoom1",
" ", 100, min80, 1000.);
6706 fs_->
make<TH1F>(
"h_sumphiEta16Depth3RADDAM_HED2ALL",
" ", 100, min90, 70. * max90);
6763 h_Amplitude_HF =
fs_->
make<TH1F>(
"h_Amplitude_HF",
" ", 100, 0., 5.);
6854 float zaR = 10000000.;
6855 float zaR2 = 50000000.;
6889 fs_->
make<TH1F>(
"h_sumADCAmplEtaPhiLs_lscounterM1orbitNum",
" ", 600, 1., 601.);
6897 float bac2 = bac + 1.;
7857 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HB",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7859 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HB",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7861 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HE",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7863 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HE",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7865 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HF",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7867 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HF",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7869 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HO",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7871 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HO",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7920 float spls0 = 10000;
7930 float spls1 = 100000;
7966 std::cout <<
"======================== booking DONE +++++++++++++++++++++++++++" << std::endl;
7975 std::cout <<
"======================== beignJob finish +++++++++++++++++++++++++++" << std::endl;
7988 int mdepth = cell.
depth();
7989 int iphi = cell.iphi() - 1;
7990 int ieta = cell.ieta();
7993 int sub = cell.subdet();
7995 int errorGeneral = 0;
8004 bool anycapid =
true;
8010 int lastcapid = 0, capid = 0;
8011 int ERRORfiber = -10;
8012 int ERRORfiberChan = -10;
8013 int ERRORfiberAndChan = -10;
8014 int repetedcapid = 0;
8016 TSsize = digiItr->size();
8019 for (
int ii = 0;
ii < TSsize;
ii++) {
8020 capid = (*digiItr)[
ii].capid();
8021 bool er = (*digiItr)[
ii].er();
8022 bool dv = (*digiItr)[
ii].dv();
8023 int fiber = (*digiItr)[
ii].fiber();
8024 int fiberChan = (*digiItr)[
ii].fiberChan();
8025 int fiberAndChan = (*digiItr)[
ii].fiberAndChan();
8026 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8029 ERRORfiberChan = fiberChan;
8030 ERRORfiberAndChan = fiberAndChan;
8031 if (capid != lastcapid) {
8045 ERRORfiberChan = fiberChan;
8046 ERRORfiberAndChan = fiberAndChan;
8051 ERRORfiberChan = fiberChan;
8052 ERRORfiberAndChan = fiberAndChan;
8058 if (firstcapid == 0 && !anycapid)
8060 if (firstcapid == 1 && !anycapid)
8062 if (firstcapid == 2 && !anycapid)
8064 if (firstcapid == 3 && !anycapid)
8073 if (!anycapid && anyer)
8075 if (!anycapid && !anydv)
8077 if (!anycapid && anyer && !anydv)
8079 if (anyer && !anydv)
8085 for (
int ii = 0;
ii < TSsize;
ii++) {
8086 double ampldefault =
adc2fC[digiItr->sample(
ii).adc()];
8087 ampl += ampldefault;
8104 if (error1 != 0 || error2 != 0 || error3 != 0) {
8106 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
8131 if (error1 != 0 || error2 != 0 || error3 != 0) {
8133 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
8158 int ieta = hcaldetid.
ieta();
8161 int iphi = hcaldetid.
iphi() - 1;
8162 int mdepth = hcaldetid.
depth();
8163 int sub = hcaldetid.
subdet();
8167 bool anycapid =
true;
8173 int lastcapid = 0, capid = 0;
8174 int repetedcapid = 0;
8179 capid = qie11df[
ii].capid();
8180 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8182 if (capid != lastcapid) {
8198 int adc = qie11df[
ii].adc();
8204 ampl += ampldefault;
8214 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
8235 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
8263 int mdepth = cell.
depth();
8264 int iphi = cell.iphi() - 1;
8265 int ieta = cell.ieta();
8268 int sub = cell.subdet();
8270 std::cout <<
" HF DIGI ?????????????? ERROR mdepth = " << mdepth << std::endl;
8272 int errorGeneral = 0;
8281 bool anycapid =
true;
8287 int lastcapid = 0, capid = 0;
8288 int ERRORfiber = -10;
8289 int ERRORfiberChan = -10;
8290 int ERRORfiberAndChan = -10;
8291 int repetedcapid = 0;
8294 TSsize = digiItr->size();
8296 for (
int ii = 0;
ii < TSsize;
ii++) {
8297 capid = (*digiItr)[
ii].capid();
8298 bool er = (*digiItr)[
ii].er();
8299 bool dv = (*digiItr)[
ii].dv();
8300 int fiber = (*digiItr)[
ii].fiber();
8301 int fiberChan = (*digiItr)[
ii].fiberChan();
8302 int fiberAndChan = (*digiItr)[
ii].fiberAndChan();
8303 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8306 ERRORfiberChan = fiberChan;
8307 ERRORfiberAndChan = fiberAndChan;
8308 if (capid != lastcapid) {
8320 ERRORfiberChan = fiberChan;
8321 ERRORfiberAndChan = fiberAndChan;
8326 ERRORfiberChan = fiberChan;
8327 ERRORfiberAndChan = fiberAndChan;
8331 if (firstcapid == 0 && !anycapid)
8333 if (firstcapid == 1 && !anycapid)
8335 if (firstcapid == 2 && !anycapid)
8337 if (firstcapid == 3 && !anycapid)
8345 if (!anycapid && anyer)
8347 if (!anycapid && !anydv)
8349 if (!anycapid && anyer && !anydv)
8351 if (anyer && !anydv)
8355 for (
int ii = 0;
ii < TSsize;
ii++) {
8356 double ampldefault =
adc2fC[digiItr->sample(
ii).adc()];
8357 ampl += ampldefault;
8371 if (error1 != 0 || error2 != 0 || error3 != 0) {
8373 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
8397 int ieta = hcaldetid.
ieta();
8400 int iphi = hcaldetid.
iphi() - 1;
8401 int mdepth = hcaldetid.
depth();
8402 int sub = hcaldetid.
subdet();
8406 bool anycapid =
true;
8412 int lastcapid = 0, capid = 0;
8413 int repetedcapid = 0;
8418 capid = qie10df[
ii].capid();
8419 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8421 if (capid != lastcapid) {
8437 int adc = qie10df[
ii].adc();
8439 ampl += ampldefault;
8449 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
8470 int mdepth = cell.
depth();
8471 int iphi = cell.iphi() - 1;
8472 int ieta = cell.ieta();
8475 int sub = cell.subdet();
8476 int errorGeneral = 0;
8485 bool anycapid =
true;
8491 int lastcapid = 0, capid = 0;
8492 int ERRORfiber = -10;
8493 int ERRORfiberChan = -10;
8494 int ERRORfiberAndChan = -10;
8495 int repetedcapid = 0;
8496 for (
int ii = 0;
ii < (*digiItr).size();
ii++) {
8497 capid = (*digiItr)[
ii].capid();
8498 bool er = (*digiItr)[
ii].er();
8499 bool dv = (*digiItr)[
ii].dv();
8500 int fiber = (*digiItr)[
ii].fiber();
8501 int fiberChan = (*digiItr)[
ii].fiberChan();
8502 int fiberAndChan = (*digiItr)[
ii].fiberAndChan();
8503 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8506 ERRORfiberChan = fiberChan;
8507 ERRORfiberAndChan = fiberAndChan;
8508 if (capid != lastcapid) {
8522 ERRORfiberChan = fiberChan;
8523 ERRORfiberAndChan = fiberAndChan;
8528 ERRORfiberChan = fiberChan;
8529 ERRORfiberAndChan = fiberAndChan;
8535 if (firstcapid == 0 && !anycapid)
8537 if (firstcapid == 1 && !anycapid)
8539 if (firstcapid == 2 && !anycapid)
8541 if (firstcapid == 3 && !anycapid)
8550 if (!anycapid && anyer)
8552 if (!anycapid && !anydv)
8554 if (!anycapid && anyer && !anydv)
8556 if (anyer && !anydv)
8561 for (
int ii = 0;
ii < (*digiItr).size();
ii++) {
8562 double ampldefault =
adc2fC[digiItr->sample(
ii).adc()];
8563 ampl += ampldefault;
8579 if (error1 != 0 || error2 != 0 || error3 != 0) {
8581 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
8588 std::cout <<
" mdepth HO = " << mdepth << std::endl;
8601 double toolwithPedSubtr[100];
8602 double lintoolwithoutPedSubtr[100];
8604 int mdepth = cell.
depth();
8605 int iphi = cell.iphi() - 1;
8606 int ieta0 = cell.ieta();
8610 int sub = cell.subdet();
8622 coder.
adc2fC(*digiItr, toolOriginal);
8623 double pedestalaver9 = 0.;
8624 double pedestalaver4 = 0.;
8625 double pedestal0 = 0.;
8626 double pedestal1 = 0.;
8627 double pedestal2 = 0.;
8628 double pedestal3 = 0.;
8629 double pedestalwaver9 = 0.;
8630 double pedestalwaver4 = 0.;
8631 double pedestalw0 = 0.;
8632 double pedestalw1 = 0.;
8633 double pedestalw2 = 0.;
8634 double pedestalw3 = 0.;
8635 double difpedestal0 = 0.;
8636 double difpedestal1 = 0.;
8637 double difpedestal2 = 0.;
8638 double difpedestal3 = 0.;
8640 double amplitudewithPedSubtr1 = 0.;
8641 double amplitudewithPedSubtr2 = 0.;
8642 double amplitudewithPedSubtr3 = 0.;
8643 double amplitudewithPedSubtr4 = 0.;
8644 double amplitude = 0.;
8645 double absamplitude = 0.;
8646 double amplitude345 = 0.;
8648 double linamplitudewithoutPedSubtr = 0.;
8651 double max_signal = -100.;
8652 int ts_with_max_signal = -100;
8658 double errorBtype = 0.;
8661 TSsize = digiItr->size();
8662 if ((*digiItr).size() != 10)
8665 for (
int ii = 0;
ii < TSsize;
ii++) {
8667 double ampldefaultwithPedSubtr = 0.;
8668 double linampldefaultwithoutPedSubtr = 0.;
8669 double ampldefault = 0.;
8670 double ampldefault0 = 0.;
8671 double ampldefault1 = 0.;
8672 double ampldefault2 = 0.;
8673 ampldefault0 =
adc2fC[digiItr->sample(
ii).adc()];
8675 ampldefault1 = toolOriginal[
ii];
8676 ampldefault2 = (*digiItr)[
ii].adc();
8678 ampldefault = ampldefault0;
8681 ampldefault = ampldefault1;
8684 ampldefault = ampldefault2;
8686 ampldefaultwithPedSubtr = ampldefault0;
8687 linampldefaultwithoutPedSubtr = ampldefault2;
8689 int capid = ((*digiItr)[
ii]).capid();
8691 double pedestalINI = pedestal00->
getValue(capid);
8693 double pedestalw = pedw->
getSigma(capid, capid);
8694 ampldefaultwithPedSubtr -=
pedestal;
8698 tool[
ii] = ampldefault;
8699 toolwithPedSubtr[
ii] = ampldefaultwithPedSubtr;
8700 lintoolwithoutPedSubtr[
ii] = linampldefaultwithoutPedSubtr;
8703 pedestalwaver9 += pedestalw * pedestalw;
8705 if (capid == 0 && c0 == 0) {
8710 pedestalw0 = pedestalw;
8711 pedestalwaver4 += pedestalw * pedestalw;
8712 difpedestal0 = pedestal - pedestalINI;
8715 if (capid == 1 && c1 == 0) {
8720 pedestalw1 = pedestalw;
8721 pedestalwaver4 += pedestalw * pedestalw;
8722 difpedestal1 = pedestal - pedestalINI;
8724 if (capid == 2 && c2 == 0) {
8729 pedestalw2 = pedestalw;
8730 pedestalwaver4 += pedestalw * pedestalw;
8731 difpedestal2 = pedestal - pedestalINI;
8733 if (capid == 3 && c3 == 0) {
8738 pedestalw3 = pedestalw;
8739 pedestalwaver4 += pedestalw * pedestalw;
8740 difpedestal3 = pedestal - pedestalINI;
8743 if (max_signal < ampldefault) {
8744 max_signal = ampldefault;
8745 ts_with_max_signal =
ii;
8756 if (
ii > 0 &&
ii < 9)
8757 amplitudewithPedSubtr1 += ampldefaultwithPedSubtr;
8759 if (
ii > 1 &&
ii < 8)
8760 amplitudewithPedSubtr2 += ampldefaultwithPedSubtr;
8762 if (
ii > 2 &&
ii < 7)
8763 amplitudewithPedSubtr3 += ampldefaultwithPedSubtr;
8765 if (
ii > 2 &&
ii < 6)
8766 amplitudewithPedSubtr4 += ampldefaultwithPedSubtr;
8768 amplitude += ampldefault;
8769 absamplitude +=
abs(ampldefault);
8771 if (
ii == 3 ||
ii == 4 ||
ii == 5)
8772 amplitude345 += ampldefault;
8776 timew += (
ii + 1) *
abs(ampldefault);
8777 timeww += (
ii + 1) * ampldefault;
8781 if (sub == 1 || sub == 2) {
8786 pedestalaver9 /= TSsize;
8787 pedestalaver4 /= c4;
8788 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
8789 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
8790 if (ts_with_max_signal > -1 && ts_with_max_signal < 10)
8791 ampl = tool[ts_with_max_signal];
8792 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < 10)
8793 ampl += tool[ts_with_max_signal + 2];
8794 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < 10)
8795 ampl += tool[ts_with_max_signal + 1];
8796 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < 10)
8797 ampl += tool[ts_with_max_signal - 1];
8800 if (ts_with_max_signal > -1 && ts_with_max_signal < 10)
8801 linamplitudewithoutPedSubtr = lintoolwithoutPedSubtr[ts_with_max_signal];
8802 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < 10)
8803 linamplitudewithoutPedSubtr += lintoolwithoutPedSubtr[ts_with_max_signal + 2];
8804 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < 10)
8805 linamplitudewithoutPedSubtr += lintoolwithoutPedSubtr[ts_with_max_signal + 1];
8806 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < 10)
8807 linamplitudewithoutPedSubtr += lintoolwithoutPedSubtr[ts_with_max_signal - 1];
8810 if (amplitude != 0.)
8811 ratio = ampl / amplitude;
8812 if (ratio < 0. || ratio > 1.02)
8814 double aveamplitude = 0.;
8815 double aveamplitudew = 0.;
8816 if (absamplitude > 0 && timew > 0)
8817 aveamplitude = timew / absamplitude;
8818 if (amplitude > 0 && timeww > 0)
8819 aveamplitudew = timeww / amplitude;
8823 bool anycapid =
true;
8828 for (
int ii = 0;
ii < TSsize;
ii++) {
8829 double aaaaaa = (
ii + 1) - aveamplitudew;
8830 double aaaaaa2 = aaaaaa * aaaaaa;
8831 double ampldefault = tool[
ii];
8832 rmsamp += (aaaaaa2 * ampldefault);
8833 capid = ((*digiItr)[
ii]).capid();
8834 bool er = (*digiItr)[
ii].er();
8835 bool dv = (*digiItr)[
ii].dv();
8836 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8848 if (!anycapid || anyer || !anydv)
8851 double rmsamplitude = 0.;
8852 if ((amplitude > 0 && rmsamp > 0) || (amplitude < 0 && rmsamp < 0))
8853 rmsamplitude =
sqrt(rmsamp / amplitude);
8854 double aveamplitude1 = aveamplitude - 1;
8860 if (sub == 1 && mdepth == 1)
8862 if (sub == 1 && mdepth == 2)
8864 if (sub == 2 && mdepth == 1)
8866 if (sub == 2 && mdepth == 2)
8868 if (sub == 2 && mdepth == 3)
8872 if (sub == 1 && mdepth == 1)
8874 if (sub == 1 && mdepth == 2)
8876 if (sub == 2 && mdepth == 1)
8878 if (sub == 2 && mdepth == 2)
8880 if (sub == 2 && mdepth == 3)
8884 for (
int ii = 0;
ii < TSsize;
ii++) {
8886 double ampldefault = tool[
ii];
8889 if (amplitude > 120) {
8932 for (
int ii = 0;
ii < TSsize;
ii++) {
8933 double ampldefault = tool[
ii];
8940 for (
int ii = 0;
ii < TSsize;
ii++) {
8941 double ampldefault = tool[
ii];
8953 for (
int ii = 0;
ii < TSsize;
ii++) {
8954 double ampldefault = tool[
ii];
8961 for (
int ii = 0;
ii < TSsize;
ii++) {
8962 double ampldefault = tool[
ii];
8972 sumEstimator6[sub - 1][mdepth - 1][ieta + 41][iphi] += errorBtype;
8973 sumEstimator0[sub - 1][mdepth - 1][ieta + 41][iphi] += pedestal0;
8978 double mypedestal = pedestal0;
8979 double mypedestalw = pedestalw0;
9050 for (
int ii = 0;
ii < TSsize;
ii++) {
9083 if (amplitude < ADCAmplHBMin_ || amplitude >
ADCAmplHBMax_) {
9085 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
9092 if (amplitude < 35.) {
9112 sumEstimator1[sub - 1][mdepth - 1][ieta + 41][iphi] += amplitude;
9119 if (aveamplitude1 < TSmeanHBMin_ || aveamplitude1 >
TSmeanHBMax_) {
9121 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
9133 sumEstimator2[sub - 1][mdepth - 1][ieta + 41][iphi] += aveamplitude1;
9140 if (ts_with_max_signal < TSpeakHBMin_ || ts_with_max_signal >
TSpeakHBMax_) {
9142 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
9154 sumEstimator3[sub - 1][mdepth - 1][ieta + 41][iphi] += float(ts_with_max_signal);
9161 if (rmsamplitude < rmsHBMin_ || rmsamplitude >
rmsHBMax_) {
9163 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
9175 sumEstimator4[sub - 1][mdepth - 1][ieta + 41][iphi] += rmsamplitude;
9184 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
9197 sumEstimator5[sub - 1][mdepth - 1][ieta + 41][iphi] += ratio;
9220 double mypedestal = pedestal0;
9221 double mypedestalw = pedestalw0;
9292 for (
int ii = 0;
ii < TSsize;
ii++) {
9314 if (amplitude < ADCAmplHEMin_ || amplitude >
ADCAmplHEMax_) {
9316 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
9325 if (amplitude < 500.) {
9364 sumEstimator1[sub - 1][mdepth - 1][ieta + 41][iphi] += amplitude;
9371 if (aveamplitude1 < TSmeanHEMin_ || aveamplitude1 >
TSmeanHEMax_) {
9373 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
9389 sumEstimator2[sub - 1][mdepth - 1][ieta + 41][iphi] += aveamplitude1;
9395 if (ts_with_max_signal < TSpeakHEMin_ || ts_with_max_signal >
TSpeakHEMax_) {
9397 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
9413 sumEstimator3[sub - 1][mdepth - 1][ieta + 41][iphi] += float(ts_with_max_signal);
9419 if (rmsamplitude < rmsHEMin_ || rmsamplitude >
rmsHEMax_) {
9421 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
9437 sumEstimator4[sub - 1][mdepth - 1][ieta + 41][iphi] += rmsamplitude;
9446 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
9462 sumEstimator5[sub - 1][mdepth - 1][ieta + 41][iphi] += ratio;
9479 double amplitudewithPedSubtr = 0.;
9482 if (ts_with_max_signal > -1 && ts_with_max_signal < 10)
9483 amplitudewithPedSubtr = toolwithPedSubtr[ts_with_max_signal];
9484 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < 10)
9485 amplitudewithPedSubtr += toolwithPedSubtr[ts_with_max_signal + 2];
9486 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < 10)
9487 amplitudewithPedSubtr += toolwithPedSubtr[ts_with_max_signal + 1];
9488 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < 10)
9489 amplitudewithPedSubtr += toolwithPedSubtr[ts_with_max_signal - 1];
9495 if (amplitudewithPedSubtr > 50.) {
9497 mapRADDAM_HE[mdepth - 1][ieta + 41][iphi] += amplitudewithPedSubtr;
9518 int LLLLLL111111 = 0;
9519 if ((mdepth == 1 && fabs(ieta0) > 16 && fabs(ieta0) < 30))
9522 int LLLLLL222222 = 0;
9523 if ((mdepth == 2 && fabs(ieta0) > 16 && fabs(ieta0) < 27) ||
9524 (mdepth == 3 && fabs(ieta0) > 26 && fabs(ieta0) < 29))
9527 if (LLLLLL111111 == 1) {
9532 if (LLLLLL222222 == 1) {
9539 if (mdepth == 3 && fabs(ieta0) == 16) {
9564 int ieta = hcaldetid.
ieta();
9567 int iphi = hcaldetid.
iphi() - 1;
9568 int mdepth = hcaldetid.
depth();
9569 int sub = hcaldetid.
subdet();
9572 if (mdepth == 0 || sub > 4)
9597 coder.
adc2fC(qie11df, toolOriginal);
9598 double pedestalaver9 = 0.;
9599 double pedestalaver4 = 0.;
9600 double pedestal0 = 0.;
9601 double pedestal1 = 0.;
9602 double pedestal2 = 0.;
9603 double pedestal3 = 0.;
9604 double pedestalwaver9 = 0.;
9605 double pedestalwaver4 = 0.;
9606 double pedestalw0 = 0.;
9607 double pedestalw1 = 0.;
9608 double pedestalw2 = 0.;
9609 double pedestalw3 = 0.;
9610 double difpedestal0 = 0.;
9611 double difpedestal1 = 0.;
9612 double difpedestal2 = 0.;
9613 double difpedestal3 = 0.;
9615 double amplitude = 0.;
9616 double amplitude0 = 0.;
9617 double absamplitude = 0.;
9618 double tocampl = 0.;
9620 double amplitude345 = 0.;
9622 double ampl3ts = 0.;
9623 double amplmaxts = 0.;
9626 double max_signal = -100.;
9627 int ts_with_max_signal = -100;
9633 double errorBtype = 0.;
9644 for (
int ii = 0;
ii < TSsize;
ii++) {
9645 double ampldefault = 0.;
9646 double tocdefault = 0.;
9647 double ampldefault0 = 0.;
9648 double ampldefault1 = 0.;
9649 double ampldefault2 = 0.;
9655 ampldefault1 = toolOriginal[
ii];
9656 ampldefault2 = qie11df[
ii].adc();
9659 ampldefault = ampldefault0;
9662 ampldefault = ampldefault1;
9665 ampldefault = ampldefault2;
9667 tocdefault = ampldefault;
9669 int capid = (qie11df[
ii]).capid();
9671 double pedestalw = pedw->
getSigma(capid, capid);
9672 double pedestalINI = pedestal00->
getValue(capid);
9676 tool[
ii] = ampldefault;
9678 pedestalwaver9 += pedestalw * pedestalw;
9680 if (capid == 0 && c0 == 0) {
9685 pedestalw0 = pedestalw;
9686 pedestalwaver4 += pedestalw * pedestalw;
9687 difpedestal0 = pedestal - pedestalINI;
9690 if (capid == 1 && c1 == 0) {
9695 pedestalw1 = pedestalw;
9696 pedestalwaver4 += pedestalw * pedestalw;
9697 difpedestal1 = pedestal - pedestalINI;
9699 if (capid == 2 && c2 == 0) {
9704 pedestalw2 = pedestalw;
9705 pedestalwaver4 += pedestalw * pedestalw;
9706 difpedestal2 = pedestal - pedestalINI;
9708 if (capid == 3 && c3 == 0) {
9713 pedestalw3 = pedestalw;
9714 pedestalwaver4 += pedestalw * pedestalw;
9715 difpedestal3 = pedestal - pedestalINI;
9718 if (max_signal < ampldefault) {
9719 max_signal = ampldefault;
9720 ts_with_max_signal =
ii;
9722 amplitude += ampldefault;
9723 absamplitude +=
abs(ampldefault);
9724 tocampl += tocdefault;
9726 if (
ii == 1 ||
ii == 2 ||
ii == 3 ||
ii == 4 ||
ii == 5 ||
ii == 6 ||
ii == 7 ||
ii == 8)
9727 amplitude345 += ampldefault;
9732 timew += (
ii + 1) *
abs(ampldefault);
9733 timeww += (
ii + 1) * ampldefault;
9735 if (
ii == 2 && ampldefault > 0.)
9739 amplitude0 = amplitude;
9741 pedestalaver9 /= TSsize;
9742 pedestalaver4 /= c4;
9743 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
9744 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
9747 if (ts_with_max_signal > -1 && ts_with_max_signal < 10) {
9748 ampl = tool[ts_with_max_signal];
9749 ampl3ts = tool[ts_with_max_signal];
9750 amplmaxts = tool[ts_with_max_signal];
9752 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < 10) {
9753 ampl += tool[ts_with_max_signal - 1];
9754 ampl3ts += tool[ts_with_max_signal - 1];
9756 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < 10) {
9757 ampl += tool[ts_with_max_signal + 1];
9758 ampl3ts += tool[ts_with_max_signal + 1];
9760 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < 10) {
9761 ampl += tool[ts_with_max_signal + 2];
9766 double xa = amplitude / 40.;
9767 double xb = ampl / 40.;
9768 double xc = amplitude345 / 40.;
9769 double xd = ampl3ts / 40.;
9770 double xe = amplmaxts / 40.;
9771 double txa = tocampl / 40.;
9773 if (((ieta == -16 || ieta == 15) && mdepth == 4) ||
9774 ((ieta == -17 || ieta == 16) && (mdepth == 2 || mdepth == 3)) ||
9775 ((ieta == -18 || ieta == 17) && mdepth == 5)) {
9776 double c0 = 1.000000;
9777 double b1 = 2.59096e-05;
9778 double a2 = 4.60721e-11;
9779 double corrforxa = a2 * xa * xa + b1 * xa + c0;
9780 double corrforxb = a2 * xb * xb + b1 * xb + c0;
9781 double corrforxc = a2 * xc * xc + b1 * xc + c0;
9782 double corrforxd = a2 * xd * xd + b1 * xd + c0;
9783 double corrforxe = a2 * xe * xe + b1 * xe + c0;
9784 double corrfortxa = a2 * txa * txa + b1 * txa + c0;
9787 amplitude *= corrforxa;
9789 amplitude345 *= corrforxc;
9790 ampl3ts *= corrforxd;
9791 amplmaxts *= corrforxe;
9792 tocampl *= corrfortxa;
9796 double c0 = 1.000000;
9797 double b1 = 2.71238e-05;
9798 double a2 = 1.32877e-10;
9799 double corrforxa = a2 * xa * xa + b1 * xa + c0;
9800 double corrforxb = a2 * xb * xb + b1 * xb + c0;
9801 double corrforxc = a2 * xc * xc + b1 * xc + c0;
9802 double corrforxd = a2 * xd * xd + b1 * xd + c0;
9803 double corrforxe = a2 * xe * xe + b1 * xe + c0;
9804 double corrfortxa = a2 * txa * txa + b1 * txa + c0;
9807 amplitude *= corrforxa;
9809 amplitude345 *= corrforxc;
9810 ampl3ts *= corrforxd;
9811 amplmaxts *= corrforxe;
9812 tocampl *= corrfortxa;
9818 if (sub == 1 || sub == 2) {
9827 if (amplitude != 0.)
9828 ratio = ampl / amplitude;
9829 if (ratio < 0. || ratio > 1.02)
9831 double aveamplitude = 0.;
9832 double aveamplitudew = 0.;
9833 if (absamplitude > 0 && timew > 0)
9834 aveamplitude = timew / absamplitude;
9835 if (amplitude0 > 0 && timeww > 0)
9836 aveamplitudew = timeww / amplitude0;
9840 bool anycapid =
true;
9843 for (
int ii = 0;
ii < TSsize;
ii++) {
9844 double aaaaaa = (
ii + 1) - aveamplitudew;
9845 double aaaaaa2 = aaaaaa * aaaaaa;
9846 double ampldefault = tool[
ii];
9847 rmsamp += (aaaaaa2 * ampldefault);
9848 capid = (qie11df[
ii]).capid();
9849 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
9857 double rmsamplitude = 0.;
9858 if ((amplitude0 > 0 && rmsamp > 0) || (amplitude0 < 0 && rmsamp < 0))
9859 rmsamplitude =
sqrt(rmsamp / amplitude0);
9860 double aveamplitude1 = aveamplitude - 1;
9865 if (sub == 1 && mdepth == 1)
9867 if (sub == 1 && mdepth == 2)
9869 if (sub == 2 && mdepth == 1)
9871 if (sub == 2 && mdepth == 2)
9873 if (sub == 2 && mdepth == 3)
9877 if (sub == 1 && mdepth == 1)
9879 if (sub == 1 && mdepth == 2)
9881 if (sub == 2 && mdepth == 1)
9883 if (sub == 2 && mdepth == 2)
9885 if (sub == 2 && mdepth == 3)
9889 for (
int ii = 0;
ii < TSsize;
ii++) {
9891 double ampldefault = tool[
ii];
9894 if (amplitude0 > 120) {
9933 if (error == 1 || amplitude0 < ADCAmplHBMin_ || amplitude0 >
ADCAmplHBMax_ || rmsamplitude <
rmsHBMin_ ||
9937 for (
int ii = 0;
ii < TSsize;
ii++) {
9938 double ampldefault = tool[
ii];
9945 for (
int ii = 0;
ii < TSsize;
ii++) {
9946 double ampldefault = tool[
ii];
9956 if (error == 1 || amplitude0 < ADCAmplHEMin_ || amplitude0 >
ADCAmplHEMax_ || rmsamplitude <
rmsHEMin_ ||
9960 for (
int ii = 0;
ii < TSsize;
ii++) {
9961 double ampldefault = tool[
ii];
9968 for (
int ii = 0;
ii < TSsize;
ii++) {
9969 double ampldefault = tool[
ii];
9978 sumEstimator6[sub - 1][mdepth - 1][ieta + 41][iphi] += errorBtype;
9980 sumEstimator0[sub - 1][mdepth - 1][ieta + 41][iphi] += pedestal0;
9985 double mypedestal = pedestal0;
9986 double mypedestalw = pedestalw0;
10057 for (
int ii = 0;
ii < TSsize;
ii++) {
10090 if (amplitude < ADCAmplHBMin_ || amplitude >
ADCAmplHBMax_) {
10092 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
10103 if (amplitude < 35.) {
10137 sumEstimator1[sub - 1][mdepth - 1][ieta + 41][iphi] += amplitude;
10144 if (aveamplitude1 < TSmeanHBMin_ || aveamplitude1 >
TSmeanHBMax_) {
10146 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
10166 sumEstimator2[sub - 1][mdepth - 1][ieta + 41][iphi] += aveamplitude1;
10172 h_TSmaxA_HB->Fill(
float(ts_with_max_signal), 1.);
10173 if (ts_with_max_signal < TSpeakHBMin_ || ts_with_max_signal >
TSpeakHBMax_) {
10175 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
10195 sumEstimator3[sub - 1][mdepth - 1][ieta + 41][iphi] += float(ts_with_max_signal);
10202 if (rmsamplitude < rmsHBMin_ || rmsamplitude >
rmsHBMax_) {
10204 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
10224 sumEstimator4[sub - 1][mdepth - 1][ieta + 41][iphi] += rmsamplitude;
10231 if (ratio < ratioHBMin_ || ratio >
ratioHBMax_) {
10233 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
10254 sumEstimator5[sub - 1][mdepth - 1][ieta + 41][iphi] += ratio;
10280 if (flagTS2 == 1) {
10299 double mypedestal = pedestal0;
10300 double mypedestalw = pedestalw0;
10387 for (
int ii = 0;
ii < TSsize;
ii++) {
10405 if (amplitude < ADCAmplHEMin_ || amplitude >
ADCAmplHEMax_) {
10407 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
10424 if (amplitude < 500.) {
10445 if (amplitude > 110 && amplitude < 150) {
10448 }
else if (amplitude > 150 && amplitude < 190) {
10453 if (ampl3ts > 70 && ampl3ts < 110) {
10456 }
else if (ampl3ts > 110 && ampl3ts < 150) {
10460 if (ampl > 87 && ampl < 127) {
10463 }
else if (ampl > 127 && ampl < 167) {
10517 sumEstimator1[sub - 1][mdepth - 1][ieta + 41][iphi] += amplitude;
10522 if (aveamplitude1 < TSmeanHEMin_ || aveamplitude1 >
TSmeanHEMax_) {
10524 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
10556 sumEstimator2[sub - 1][mdepth - 1][ieta + 41][iphi] += aveamplitude1;
10561 h_TSmaxA_HE->Fill(
float(ts_with_max_signal), 1.);
10562 if (ts_with_max_signal < TSpeakHEMin_ || ts_with_max_signal >
TSpeakHEMax_) {
10564 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
10596 sumEstimator3[sub - 1][mdepth - 1][ieta + 41][iphi] += float(ts_with_max_signal);
10602 if (rmsamplitude < rmsHEMin_ || rmsamplitude >
rmsHEMax_) {
10604 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
10636 sumEstimator4[sub - 1][mdepth - 1][ieta + 41][iphi] += rmsamplitude;
10643 if (ratio < ratioHEMin_ || ratio >
ratioHEMax_) {
10645 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
10677 sumEstimator5[sub - 1][mdepth - 1][ieta + 41][iphi] += ratio;
10713 if (flagTS2 == 1) {
10738 int mdepth = cell.
depth();
10739 int iphi = cell.iphi() - 1;
10740 int ieta = cell.ieta();
10743 int sub = cell.subdet();
10753 coder.
adc2fC(*digiItr, toolOriginal);
10754 double pedestalaver9 = 0.;
10755 double pedestalaver4 = 0.;
10756 double pedestal0 = 0.;
10757 double pedestal1 = 0.;
10758 double pedestal2 = 0.;
10759 double pedestal3 = 0.;
10760 double pedestalwaver9 = 0.;
10761 double pedestalwaver4 = 0.;
10762 double pedestalw0 = 0.;
10763 double pedestalw1 = 0.;
10764 double pedestalw2 = 0.;
10765 double pedestalw3 = 0.;
10766 double amplitude = 0.;
10767 double absamplitude = 0.;
10770 double timeww = 0.;
10771 double max_signal = -100.;
10772 int ts_with_max_signal = -100;
10778 double errorBtype = 0.;
10780 if ((*digiItr).size() != TSsize)
10782 TSsize = digiItr->size();
10783 for (
int ii = 0;
ii < TSsize;
ii++) {
10785 double ampldefault = 0.;
10786 double ampldefault0 = 0.;
10787 double ampldefault1 = 0.;
10788 double ampldefault2 = 0.;
10789 ampldefault0 =
adc2fC[digiItr->sample(
ii).adc()];
10791 ampldefault1 = toolOriginal[
ii];
10792 ampldefault2 = (*digiItr)[
ii].adc();
10794 ampldefault = ampldefault0;
10797 ampldefault = ampldefault1;
10800 ampldefault = ampldefault2;
10803 int capid = ((*digiItr)[
ii]).capid();
10806 double pedestalw = pedw->
getSigma(capid, capid);
10810 tool[
ii] = ampldefault;
10813 pedestalwaver9 += pedestalw * pedestalw;
10815 if (capid == 0 && c0 == 0) {
10820 pedestalw0 = pedestalw;
10821 pedestalwaver4 += pedestalw * pedestalw;
10824 if (capid == 1 && c1 == 0) {
10829 pedestalw1 = pedestalw;
10830 pedestalwaver4 += pedestalw * pedestalw;
10832 if (capid == 2 && c2 == 0) {
10837 pedestalw2 = pedestalw;
10838 pedestalwaver4 += pedestalw * pedestalw;
10840 if (capid == 3 && c3 == 0) {
10845 pedestalw3 = pedestalw;
10846 pedestalwaver4 += pedestalw * pedestalw;
10849 if (max_signal < ampldefault) {
10850 max_signal = ampldefault;
10851 ts_with_max_signal =
ii;
10853 amplitude += ampldefault;
10854 absamplitude +=
abs(ampldefault);
10859 timew += (
ii + 1) *
abs(ampldefault);
10860 timeww += (
ii + 1) * ampldefault;
10870 pedestalaver9 /= TSsize;
10871 pedestalaver4 /= c4;
10872 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
10873 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
10876 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
10877 ampl = tool[ts_with_max_signal];
10878 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
10879 ampl += tool[ts_with_max_signal + 2];
10880 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
10881 ampl += tool[ts_with_max_signal + 1];
10882 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
10883 ampl += tool[ts_with_max_signal - 1];
10887 if (amplitude != 0.)
10888 ratio = ampl / amplitude;
10890 if (ratio < 0. || ratio > 1.02)
10893 double aveamplitude = 0.;
10894 double aveamplitudew = 0.;
10895 if (absamplitude > 0 && timew > 0)
10896 aveamplitude = timew / absamplitude;
10897 if (amplitude > 0 && timeww > 0)
10898 aveamplitudew = timeww / amplitude;
10900 double rmsamp = 0.;
10903 bool anycapid =
true;
10904 bool anyer =
false;
10908 for (
int ii = 0;
ii < TSsize;
ii++) {
10909 double aaaaaa = (
ii + 1) - aveamplitudew;
10910 double aaaaaa2 = aaaaaa * aaaaaa;
10911 double ampldefault = tool[
ii];
10912 rmsamp += (aaaaaa2 * ampldefault);
10913 capid = ((*digiItr)[
ii]).capid();
10914 bool er = (*digiItr)[
ii].er();
10915 bool dv = (*digiItr)[
ii].dv();
10916 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
10929 if (!anycapid || anyer || !anydv)
10931 double rmsamplitude = 0.;
10932 if ((amplitude > 0 && rmsamp > 0) || (amplitude < 0 && rmsamp < 0))
10933 rmsamplitude =
sqrt(rmsamp / amplitude);
10934 double aveamplitude1 = aveamplitude - 1;
10936 if (sub == 4 && mdepth == 1)
10938 if (sub == 4 && mdepth == 2)
10942 if (sub == 4 && mdepth == 1)
10944 if (sub == 4 && mdepth == 2)
10956 for (
int ii = 0;
ii < TSsize;
ii++) {
10957 double ampldefault = tool[
ii];
10964 for (
int ii = 0;
ii < TSsize;
ii++) {
10965 double ampldefault = tool[
ii];
10973 sumEstimator6[sub - 1][mdepth - 1][ieta + 41][iphi] += errorBtype;
10974 sumEstimator0[sub - 1][mdepth - 1][ieta + 41][iphi] += pedestal0;
10979 double mypedestal = pedestal0;
10980 double mypedestalw = pedestalw0;
11040 for (
int ii = 0;
ii < TSsize;
ii++) {
11057 if (amplitude < ADCAmplHFMin_ || amplitude >
ADCAmplHFMax_) {
11059 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
11066 if (amplitude < 20.) {
11088 sumEstimator1[sub - 1][mdepth - 1][ieta + 41][iphi] += amplitude;
11095 if (aveamplitude1 < TSmeanHFMin_ || aveamplitude1 >
TSmeanHFMax_) {
11097 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
11110 sumEstimator2[sub - 1][mdepth - 1][ieta + 41][iphi] += aveamplitude1;
11115 h_TSmaxA_HF->Fill(
float(ts_with_max_signal), 1.);
11116 if (ts_with_max_signal < TSpeakHFMin_ || ts_with_max_signal >
TSpeakHFMax_) {
11118 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
11130 sumEstimator3[sub - 1][mdepth - 1][ieta + 41][iphi] += float(ts_with_max_signal);
11136 if (rmsamplitude < rmsHFMin_ || rmsamplitude >
rmsHFMax_) {
11138 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
11150 sumEstimator4[sub - 1][mdepth - 1][ieta + 41][iphi] += rmsamplitude;
11156 if (ratio < ratioHFMin_ || ratio >
ratioHFMax_) {
11158 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
11171 sumEstimator5[sub - 1][mdepth - 1][ieta + 41][iphi] += ratio;
11200 int ieta = hcaldetid.
ieta();
11203 int iphi = hcaldetid.
iphi() - 1;
11204 int mdepth = hcaldetid.
depth();
11205 int sub = hcaldetid.
subdet();
11220 if (mdepth == 0 || sub != 4)
11240 coder.
adc2fC(qie10df, toolOriginal);
11243 double pedestalaver9 = 0.;
11244 double pedestalaver4 = 0.;
11245 double pedestal0 = 0.;
11246 double pedestal1 = 0.;
11247 double pedestal2 = 0.;
11248 double pedestal3 = 0.;
11249 double pedestalwaver9 = 0.;
11250 double pedestalwaver4 = 0.;
11251 double pedestalw0 = 0.;
11252 double pedestalw1 = 0.;
11253 double pedestalw2 = 0.;
11254 double pedestalw3 = 0.;
11255 double amplitude = 0.;
11256 double absamplitude = 0.;
11259 double timeww = 0.;
11260 double max_signal = -100.;
11261 int ts_with_max_signal = -100;
11267 double errorBtype = 0.;
11274 for (
int ii = 0;
ii < TSsize;
ii++) {
11275 double ampldefault = 0.;
11276 double ampldefault0 = 0.;
11277 double ampldefault1 = 0.;
11278 double ampldefault2 = 0.;
11281 ampldefault1 = toolOriginal[
ii];
11282 ampldefault2 = qie10df[
ii].adc();
11284 ampldefault = ampldefault0;
11287 ampldefault = ampldefault1;
11290 ampldefault = ampldefault2;
11293 int capid = (qie10df[
ii]).capid();
11295 double pedestalw = pedw->
getSigma(capid, capid);
11300 tool[
ii] = ampldefault;
11303 pedestalwaver9 += pedestalw * pedestalw;
11305 if (capid == 0 && c0 == 0) {
11310 pedestalw0 = pedestalw;
11311 pedestalwaver4 += pedestalw * pedestalw;
11314 if (capid == 1 && c1 == 0) {
11319 pedestalw1 = pedestalw;
11320 pedestalwaver4 += pedestalw * pedestalw;
11322 if (capid == 2 && c2 == 0) {
11327 pedestalw2 = pedestalw;
11328 pedestalwaver4 += pedestalw * pedestalw;
11330 if (capid == 3 && c3 == 0) {
11335 pedestalw3 = pedestalw;
11336 pedestalwaver4 += pedestalw * pedestalw;
11339 if (max_signal < ampldefault) {
11340 max_signal = ampldefault;
11341 ts_with_max_signal =
ii;
11343 amplitude += ampldefault;
11344 absamplitude +=
abs(ampldefault);
11346 timew += (
ii + 1) *
abs(ampldefault);
11347 timeww += (
ii + 1) * ampldefault;
11348 if (
ii == 1 && ampldefault > 0.)
11359 pedestalaver9 /= TSsize;
11360 pedestalaver4 /= c4;
11361 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
11362 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
11365 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
11366 ampl = tool[ts_with_max_signal];
11367 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
11368 ampl += tool[ts_with_max_signal + 2];
11369 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
11370 ampl += tool[ts_with_max_signal + 1];
11371 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
11372 ampl += tool[ts_with_max_signal - 1];
11376 if (amplitude != 0.)
11377 ratio = ampl / amplitude;
11378 if (ratio < 0. || ratio > 1.02)
11380 double aveamplitude = 0.;
11381 double aveamplitudew = 0.;
11382 if (absamplitude > 0 && timew > 0)
11383 aveamplitude = timew / absamplitude;
11384 if (amplitude > 0 && timeww > 0)
11385 aveamplitudew = timeww / amplitude;
11387 double rmsamp = 0.;
11389 bool anycapid =
true;
11392 for (
int ii = 0;
ii < TSsize;
ii++) {
11393 double aaaaaa = (
ii + 1) - aveamplitudew;
11394 double aaaaaa2 = aaaaaa * aaaaaa;
11395 double ampldefault = tool[
ii];
11396 rmsamp += (aaaaaa2 * ampldefault);
11397 capid = (qie10df[
ii]).capid();
11398 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
11406 double rmsamplitude = 0.;
11407 if ((amplitude > 0 && rmsamp > 0) || (amplitude < 0 && rmsamp < 0))
11408 rmsamplitude =
sqrt(rmsamp / amplitude);
11409 double aveamplitude1 = aveamplitude - 1;
11414 if (sub == 4 && mdepth == 1)
11416 if (sub == 4 && mdepth == 2)
11420 if (sub == 4 && mdepth == 1)
11422 if (sub == 4 && mdepth == 2)
11434 for (
int ii = 0;
ii < TSsize;
ii++) {
11435 double ampldefault = tool[
ii];
11442 for (
int ii = 0;
ii < TSsize;
ii++) {
11443 double ampldefault = tool[
ii];
11452 sumEstimator6[sub - 1][mdepth - 1][ieta + 41][iphi] += errorBtype;
11453 sumEstimator0[sub - 1][mdepth - 1][ieta + 41][iphi] += pedestal0;
11458 double mypedestal = pedestal0;
11459 double mypedestalw = pedestalw0;
11527 for (
int ii = 0;
ii < TSsize;
ii++) {
11544 if (amplitude < ADCAmplHFMin_ || amplitude >
ADCAmplHFMax_) {
11546 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
11557 if (amplitude < 20.) {
11591 sumEstimator1[sub - 1][mdepth - 1][ieta + 41][iphi] += amplitude;
11596 if (aveamplitude1 < TSmeanHFMin_ || aveamplitude1 >
TSmeanHFMax_) {
11598 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
11619 sumEstimator2[sub - 1][mdepth - 1][ieta + 41][iphi] += aveamplitude1;
11624 h_TSmaxA_HF->Fill(
float(ts_with_max_signal), 1.);
11625 if (ts_with_max_signal < TSpeakHFMin_ || ts_with_max_signal >
TSpeakHFMax_) {
11627 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
11648 sumEstimator3[sub - 1][mdepth - 1][ieta + 41][iphi] += float(ts_with_max_signal);
11654 if (rmsamplitude < rmsHFMin_ || rmsamplitude >
rmsHFMax_) {
11656 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
11677 sumEstimator4[sub - 1][mdepth - 1][ieta + 41][iphi] += rmsamplitude;
11683 if (ratio < ratioHFMin_ || ratio >
ratioHFMax_) {
11685 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
11706 sumEstimator5[sub - 1][mdepth - 1][ieta + 41][iphi] += ratio;
11732 if (flagTS1 == 1) {
11753 int mdepth = cell.
depth();
11754 int iphi = cell.iphi() - 1;
11755 int ieta = cell.ieta();
11758 int sub = cell.subdet();
11768 coder.
adc2fC(*digiItr, toolOriginal);
11769 double pedestalaver9 = 0.;
11770 double pedestalaver4 = 0.;
11771 double pedestal0 = 0.;
11772 double pedestal1 = 0.;
11773 double pedestal2 = 0.;
11774 double pedestal3 = 0.;
11775 double pedestalwaver9 = 0.;
11776 double pedestalwaver4 = 0.;
11777 double pedestalw0 = 0.;
11778 double pedestalw1 = 0.;
11779 double pedestalw2 = 0.;
11780 double pedestalw3 = 0.;
11781 double amplitude = 0.;
11782 double absamplitude = 0.;
11785 double timeww = 0.;
11786 double max_signal = -100.;
11787 int ts_with_max_signal = -100;
11793 double errorBtype = 0.;
11795 if ((*digiItr).size() != TSsize)
11797 TSsize = digiItr->size();
11799 for (
int ii = 0;
ii < TSsize;
ii++) {
11800 double ampldefault = 0.;
11801 double ampldefault0 = 0.;
11802 double ampldefault1 = 0.;
11803 double ampldefault2 = 0.;
11804 ampldefault0 =
adc2fC[digiItr->sample(
ii).adc()];
11806 ampldefault1 = toolOriginal[
ii];
11807 ampldefault2 = (*digiItr)[
ii].adc();
11809 ampldefault = ampldefault0;
11812 ampldefault = ampldefault1;
11815 ampldefault = ampldefault2;
11817 int capid = ((*digiItr)[
ii]).capid();
11819 double pedestalw = pedw->
getSigma(capid, capid);
11822 tool[
ii] = ampldefault;
11824 pedestalwaver9 += pedestalw * pedestalw;
11825 if (capid == 0 && c0 == 0) {
11830 pedestalw0 = pedestalw;
11831 pedestalwaver4 += pedestalw * pedestalw;
11833 if (capid == 1 && c1 == 0) {
11838 pedestalw1 = pedestalw;
11839 pedestalwaver4 += pedestalw * pedestalw;
11841 if (capid == 2 && c2 == 0) {
11846 pedestalw2 = pedestalw;
11847 pedestalwaver4 += pedestalw * pedestalw;
11849 if (capid == 3 && c3 == 0) {
11854 pedestalw3 = pedestalw;
11855 pedestalwaver4 += pedestalw * pedestalw;
11858 if (max_signal < ampldefault) {
11859 max_signal = ampldefault;
11860 ts_with_max_signal =
ii;
11862 amplitude += ampldefault;
11863 absamplitude +=
abs(ampldefault);
11867 timew += (
ii + 1) *
abs(ampldefault);
11868 timeww += (
ii + 1) * ampldefault;
11869 if (ii < 3 && ampldefault > 0.)
11874 pedestalaver9 /= TSsize;
11875 pedestalaver4 /= c4;
11876 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
11877 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
11878 if (ts_with_max_signal > -1 && ts_with_max_signal < 10)
11879 ampl = tool[ts_with_max_signal];
11880 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < 10)
11881 ampl += tool[ts_with_max_signal + 2];
11882 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < 10)
11883 ampl += tool[ts_with_max_signal + 1];
11884 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < 10)
11885 ampl += tool[ts_with_max_signal - 1];
11887 if (amplitude != 0.)
11888 ratio = ampl / amplitude;
11889 if (ratio < 0. || ratio > 1.04)
11891 double aveamplitude = 0.;
11892 double aveamplitudew = 0.;
11893 if (absamplitude > 0 && timew > 0)
11894 aveamplitude = timew / absamplitude;
11895 if (amplitude > 0 && timeww > 0)
11896 aveamplitudew = timeww / amplitude;
11897 double rmsamp = 0.;
11899 bool anycapid =
true;
11900 bool anyer =
false;
11904 for (
int ii = 0;
ii < TSsize;
ii++) {
11905 double aaaaaa = (
ii + 1) - aveamplitudew;
11906 double aaaaaa2 = aaaaaa * aaaaaa;
11907 double ampldefault = tool[
ii];
11908 rmsamp += (aaaaaa2 * ampldefault);
11909 capid = ((*digiItr)[
ii]).capid();
11910 bool er = (*digiItr)[
ii].er();
11911 bool dv = (*digiItr)[
ii].dv();
11912 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
11924 if (!anycapid || anyer || !anydv)
11926 double rmsamplitude = 0.;
11927 if ((amplitude > 0 && rmsamp > 0) || (amplitude < 0 && rmsamp < 0))
11928 rmsamplitude =
sqrt(rmsamp / amplitude);
11929 double aveamplitude1 = aveamplitude - 1;
11931 if (sub == 3 && mdepth == 4)
11935 if (sub == 3 && mdepth == 4)
11946 for (
int ii = 0;
ii < TSsize;
ii++) {
11947 double ampldefault = tool[
ii];
11952 for (
int ii = 0;
ii < TSsize;
ii++) {
11953 double ampldefault = tool[
ii];
11962 sumEstimator6[sub - 1][mdepth - 1][ieta + 41][iphi] += errorBtype;
11963 sumEstimator0[sub - 1][mdepth - 1][ieta + 41][iphi] += pedestal0;
11967 double mypedestal = pedestal0;
11968 double mypedestalw = pedestalw0;
12011 for (
int ii = 0;
ii < TSsize;
ii++) {
12027 if (amplitude < ADCAmplHOMin_ || amplitude >
ADCAmplHOMax_) {
12029 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
12035 if (amplitude < 100.) {
12051 sumEstimator1[sub - 1][mdepth - 1][ieta + 41][iphi] += amplitude;
12055 if (aveamplitude1 < TSmeanHOMin_ || aveamplitude1 >
TSmeanHOMax_) {
12057 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
12064 sumEstimator2[sub - 1][mdepth - 1][ieta + 41][iphi] += aveamplitude1;
12067 h_TSmaxA_HO->Fill(
float(ts_with_max_signal), 1.);
12068 if (ts_with_max_signal < TSpeakHOMin_ || ts_with_max_signal >
TSpeakHOMax_) {
12070 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
12078 sumEstimator3[sub - 1][mdepth - 1][ieta + 41][iphi] += float(ts_with_max_signal);
12082 if (rmsamplitude < rmsHOMin_ || rmsamplitude >
rmsHOMax_) {
12084 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
12091 sumEstimator4[sub - 1][mdepth - 1][ieta + 41][iphi] += rmsamplitude;
12095 if (ratio < ratioHOMin_ || ratio >
ratioHOMax_) {
12097 ++
badchannels[sub - 1][mdepth - 1][ieta + 41][iphi];
12104 sumEstimator5[sub - 1][mdepth - 1][ieta + 41][iphi] += ratio;
12112 if (flagTS012 == 1)
12119 if (kdet == 1 || kdet == 2) {
12123 cal_RBX = (kphi + 1) / 4;
12124 cal_RBX = cal_RBX + 18 * (keta + 1) / 2;
12127 cal_RBX = (int)(kphi / 18) + 1;
12134 cal_RBX = kphi / 12 + 1;
12140 cal_RBX = kphi / 12 + 1 + 6;
12146 cal_RBX = kphi / 6 + 1 + 12;
12152 cal_RBX = kphi / 12 + 1 + 24;
12158 cal_RBX = kphi / 12 + 1 + 30;
12179 for (
int k1 = 0; k1 <
ndepth; k1++) {
12181 for (
int k3 = 0; k3 <
nphi; k3++) {
12182 int ieta =
k2 - 41;
12277 int kkkk2 = (
k2 - 1) / 4;
12282 int kkkk3 = (k3) / 4 + 1;
12284 ietaphi = ((kkkk2)-1) *
znphi + (kkkk3);
12287 bbb3 = bbbc / bbb1;
13131 std::cout <<
" ==== Edn of run " << std::endl;
13149 MAPfile <<
"#define LogEleMapdb_h" << std::endl;
13150 MAPfile <<
"#include <algorithm>" << std::endl;
13151 MAPfile <<
"#include <iostream>" << std::endl;
13152 MAPfile <<
"#include <vector>" << std::endl;
13153 MAPfile <<
"#include <string>" << std::endl;
13154 MAPfile <<
"#include <sstream>" << std::endl;
13157 MAPfile <<
"struct Cell {" << std::endl;
13158 MAPfile <<
" std::string subdet;" << std::endl;
13159 MAPfile <<
" int Eta;" << std::endl;
13160 MAPfile <<
" int Phi;" << std::endl;
13161 MAPfile <<
" int Depth;" << std::endl;
13162 MAPfile <<
" std::string RBX;" << std::endl;
13163 MAPfile <<
" int RM;" << std::endl;
13164 MAPfile <<
" int Pixel;" << std::endl;
13165 MAPfile <<
" int RMfiber;" << std::endl;
13166 MAPfile <<
" int FiberCh;" << std::endl;
13167 MAPfile <<
" int QIE;" << std::endl;
13168 MAPfile <<
" int ADC;" << std::endl;
13169 MAPfile <<
" int VMECardID;" << std::endl;
13170 MAPfile <<
" int dccID;" << std::endl;
13171 MAPfile <<
" int Spigot;" << std::endl;
13172 MAPfile <<
" int FiberIndex;" << std::endl;
13173 MAPfile <<
" int HtrSlot;" << std::endl;
13174 MAPfile <<
" int HtrTB;" << std::endl;
13177 MAPfile <<
"// the function check, if \"par\" == \"val\" for this cell" << std::endl;
13178 MAPfile <<
" bool check(const std::string par, const int val) const " << std::endl;
13179 MAPfile <<
" {" << std::endl;
13180 MAPfile <<
" if (par == \"Eta\") return (val == Eta);" << std::endl;
13181 MAPfile <<
" else if (par == \"Phi\") return (val == Phi);" << std::endl;
13182 MAPfile <<
" else if (par == \"Depth\") return (val == Depth);" << std::endl;
13183 MAPfile <<
" else if (par == \"RM\") return (val == RM);" << std::endl;
13184 MAPfile <<
" else if (par == \"Pixel\") return (val == Pixel);" << std::endl;
13185 MAPfile <<
" else if (par == \"RMfiber\") return (val == RMfiber);" << std::endl;
13186 MAPfile <<
" else if (par == \"FiberCh\") return (val == FiberCh);" << std::endl;
13187 MAPfile <<
" else if (par == \"QIE\") return (val == QIE);" << std::endl;
13188 MAPfile <<
" else if (par == \"ADC\") return (val == ADC);" << std::endl;
13189 MAPfile <<
" else if (par == \"VMECardID\") return (val == VMECardID);" << std::endl;
13190 MAPfile <<
" else if (par == \"dccID\") return (val == dccID);" << std::endl;
13191 MAPfile <<
" else if (par == \"Spigot\") return (val == Spigot);" << std::endl;
13192 MAPfile <<
" else if (par == \"FiberIndex\") return (val == FiberIndex);" << std::endl;
13193 MAPfile <<
" else if (par == \"HtrSlot\") return (val == HtrSlot);" << std::endl;
13194 MAPfile <<
" else if (par == \"HtrTB\") return (val == HtrTB);" << std::endl;
13195 MAPfile <<
" else return false;" << std::endl;
13196 MAPfile <<
" }" << std::endl;
13199 MAPfile <<
" bool check(const std::string par, const std::string val) const" << std::endl;
13200 MAPfile <<
" {" << std::endl;
13201 MAPfile <<
" if (par == \"subdet\") return (val == subdet);" << std::endl;
13202 MAPfile <<
" else if (par == \"RBX\") return (val == RBX);" << std::endl;
13203 MAPfile <<
" else return false;" << std::endl;
13204 MAPfile <<
" }" << std::endl;
13206 MAPfile <<
"};" << std::endl;
13209 MAPfile <<
"const Cell AllCells[] = {" << std::endl;
13210 MAPfile <<
"//{ SD, Eta, Phi, Depth, RBX, RM, PIXEL, RMfiber, Fiber Ch., QIE, ADC, VMECrateId, dccid, spigot, "
13211 "fiberIndex, htrSlot, htrTopBottom }"
13224 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , "
13225 << detid->
depth() <<
" ,";
13226 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , "
13227 << lmap_entry.
rmFiber() <<
" , ";
13231 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13233 MAPfile <<
"}," << std::endl;
13246 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() - 1 <<
" , " << detid->
iphi() - 1 <<
" , "
13247 << detid->
depth() <<
" ,";
13248 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , "
13249 << lmap_entry.
rmFiber() <<
" , ";
13253 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13255 MAPfile <<
"}," << std::endl;
13262 for (
int eta = -20;
eta <= -20;
eta++) {
13270 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , "
13271 << detid->
depth() <<
" ,";
13272 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , "
13273 << lmap_entry.
rmFiber() <<
" , ";
13277 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13279 MAPfile <<
"}," << std::endl;
13285 for (
int eta = -19;
eta <= -16;
eta++) {
13293 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , "
13294 << detid->
depth() <<
" ,";
13295 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , "
13296 << lmap_entry.
rmFiber() <<
" , ";
13300 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13302 MAPfile <<
"}," << std::endl;
13307 for (
int eta = -29;
eta <= -16;
eta++) {
13315 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , "
13316 << detid->
depth() <<
" ,";
13317 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , "
13318 << lmap_entry.
rmFiber() <<
" , ";
13322 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13324 MAPfile <<
"}," << std::endl;
13337 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , "
13338 << detid->
depth() <<
" ,";
13339 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , "
13340 << lmap_entry.
rmFiber() <<
" , ";
13344 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13346 MAPfile <<
"}," << std::endl;
13354 for (
int eta = -41;
eta <= -29;
eta++) {
13362 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , "
13363 << detid->
depth() <<
" ,";
13364 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , "
13365 << lmap_entry.
rmFiber() <<
" , ";
13369 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13371 MAPfile <<
"}," << std::endl;
13385 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , "
13386 << detid->
depth() <<
" ,";
13387 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , "
13388 << lmap_entry.
rmFiber() <<
" , ";
13392 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13394 MAPfile <<
"}," << std::endl;
13410 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , "
13411 << detid->
depth() <<
" ,";
13412 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , "
13413 << lmap_entry.
rmFiber() <<
" , ";
13417 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13419 MAPfile <<
"}," << std::endl;
13433 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() - 1 <<
" , " << detid->
iphi() - 1 <<
" , "
13434 << detid->
depth() <<
" ,";
13435 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , "
13436 << lmap_entry.
rmFiber() <<
" , ";
13440 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13442 MAPfile <<
"}," << std::endl;
13447 MAPfile <<
"};" << std::endl;
13450 MAPfile <<
"// macro for array length calculation" << std::endl;
13451 MAPfile <<
"#define DIM(a) (sizeof(a)/sizeof(a[0]))" << std::endl;
13454 MAPfile <<
"// class for cells array managing" << std::endl;
13455 MAPfile <<
"class CellDB {" << std::endl;
13456 MAPfile <<
"public:" << std::endl;
13457 MAPfile <<
" CellDB()" << std::endl;
13458 MAPfile <<
" : cells(AllCells, AllCells + DIM(AllCells))" << std::endl;
13459 MAPfile <<
"{}" << std::endl;
13462 MAPfile <<
"// return i-th cell" << std::endl;
13463 MAPfile <<
"Cell operator [] (int i) const {return cells[i];}" << std::endl;
13465 MAPfile <<
"// number of cells in database" << std::endl;
13466 MAPfile <<
"int size() const {return cells.size();}" << std::endl;
13469 MAPfile <<
"// select cells for which \"par\" == \"val\"" << std::endl;
13470 MAPfile <<
"template<typename T>" << std::endl;
13471 MAPfile <<
"CellDB find(const std::string par, const T val) const" << std::endl;
13473 MAPfile <<
" std::vector<Cell> s;" << std::endl;
13474 MAPfile <<
" for (size_t i = 0; i < cells.size(); ++i)" << std::endl;
13475 MAPfile <<
" if (cells[i].check(par, val))" << std::endl;
13476 MAPfile <<
" s.push_back(cells[i]);" << std::endl;
13477 MAPfile <<
" return CellDB(s);" << std::endl;
13478 MAPfile <<
"} " << std::endl;
13480 MAPfile <<
"private:" << std::endl;
13481 MAPfile <<
" CellDB( const std::vector<Cell> s)" << std::endl;
13482 MAPfile <<
" : cells(s)" << std::endl;
13483 MAPfile <<
"{}" << std::endl;
13484 MAPfile <<
"std::vector<Cell> cells;" << std::endl;
13485 MAPfile <<
"};" << std::endl;
13487 std::cout <<
"===== Finish writing Channel MAP =====" << std::endl;
13491 double deltaphi = phi1 - phi2;
13493 deltaphi = phi2 - phi1;
13495 if (deltaphi >
M_PI) {
13496 deltaphi = 2. *
M_PI - deltaphi;
13498 double deltaeta = eta2 -
eta1;
13499 double tmp =
sqrt(deltaeta * deltaeta + deltaphi * deltaphi);
13509 if (a1 > 0.5 *
M_PI && a2 < 0.)
13511 if (a2 > 0.5 *
M_PI && a1 < 0.)
13513 double tmp = (a1 * en1 + a2 * en2) / (en1 + en2);
13526 double tmp = a2 -
a1;
13527 if (a1 * a2 < 0.) {
13528 if (a1 > 0.5 *
M_PI)
13530 if (a2 > 0.5 *
M_PI)
TH2F * h_mapDepth4Pedw1_HO
TH2F * h_mapDepth4AmplE34_HO
double lsdep_estimator2_HEdepth3_
constexpr float energy() const
TH1F * h_averSIGNALsumamplitude_HO
TH1F * h_runnbadchannels_depth1_HF
TH2F * h_recNoiseEnergy0_HF1
TH2F * h_2D0Ataildepth2_HB
TH1F * h_sumAmplitudeperLS8
double lsdep_estimator3_HFdepth1_
TH1F * h_nevents_per_LSzoom
TH1F * h_sum0PedestalperLS3
double amplitudechannel2[nsub][ndepth][neta][nphi]
static const std::string kSharedResource
TH1F * h_sumCutADCAmplperLSdepth6HEu
TH1F * h_bcnnbadchannels_depth4_HO
TH2F * h_recNoiseEnergy1_HB1
TH2F * h_2DsumAmplitudeLS5
TH1F * h_RatioOccupancy_HFM
TH2F * h_mapDepth6Amplitude225_HE
TH2F * h_2DsumADCAmplEtaPhiLs3
TH2F * h_mapWidthCalib047_HE
TH1F * h_averSIGNALsumamplitude_HF
TH2F * h_2D0sumADCAmplLS6
TH2F * h_recSignalEnergy1_HF3
bool studyRunDependenceHist_
TH2F * h_mapDepth4TSmeanA225_HE
TH2F * h_mapDepth1TSmaxA225_HE
TH1F * h_sum0ADCAmplperLS3_P1
double const adc2fC_QIE11_shunt6[256]
EventNumber_t event() const
TH2F * h_mapDepth3TSmaxA225_HE
TH2F * h_FullSignal3D0_HE
TH2F * h_mapDepth1Ped3_HB
TH2F * h_2D0sumTSmeanALS6
TH1F * h_sum0ADCAmplperLS5
TH1F * h_sumADCAmplperLS7u
int lsdep_cut1_peak_HFdepth2_
T getUntrackedParameter(std::string const &, T const &) const
TH1F * h_shape0_good_channels_HO
TH1F * h_bcnbadrate0_depth4_HO
TH2F * h_mapDepth3ADCAmpl225Copy_HF
TH2F * h_FullSignal3D0_HF
TH1F * h_bcnvsamplitude_HE
TH1F * h_RatioOccupancy_HOM
TH2F * h_amplitudechannel0_HE4
TH1F * h_sum0ADCAmplperLS8_M2
TH2F * h_mapDepth1linADCAmpl12_HE
TH1F * h_pedestalaver4_HF
TH2F * h_amplitudechannel1_HB4
TH2F * h_mapDepth2Amplitude_HB
double lsdep_estimator4_HEdepth3_
TH1F * h_runnbadchannels_depth4_HO
TH1F * h_RatioOccupancy_HOP
TH1F * h_runnbadchannels_depth1_HE
TH1F * h_sum0ErrorBperLS6
double calib3[ndepth][neta][nphi]
TH2F * h_mapDepth2Error_HB
TH2F * h_mapDepth2Ped0_HF
TH2F * h_mapDepth2Ampl047_HE
edm::EDGetTokenT< QIE10DigiCollection > tok_qie10_
TH1F * h_runbadrate0_depth3_HE
TH2F * h_amplitudechannel0_HB4
TH1F * h_orbitNumvsamplitude_HF
TH1F * h_pedestalaver4_HE
TH2F * h_mapDepth1Error_HF
TH2F * h_mapDepth5TSmeanA_HE
TH2F * h_mapDepth6Ampl_HE
TH2F * h_recNoiseEnergy1_HE5
TH1F * h_bcnvsamplitude0_HB
TH2F * h_recSignalEnergy1_HE5
TH2F * h_mapDepth2AmplE34_HE
TH2F * h_recSignalEnergy0_HF1
TH1F * h_ADCAmpl345Zoom_HB
TH2F * h_mapTSmeanCalib_HF
TH1F * h_averNOSIGNALsumamplitude_HO
TH2F * h_mapDepth1Ped2_HE
TH1F * h_runbadrate0_depth4_HO
TH1F * h_bcnbadrate0_depth1_HF
TH2F * h_mapGetRMSOverNormalizedSignal0_HE
TH2F * h_amplitudechannel0_HE3
TH2F * h_mapDepth5TSmaxA_HE
TH2F * h_2DsumPedestalLS8
TH1F * h_ndepthesperamplitudebins_HE
const HcalDbService * conditions
TH1F * h_AmplitudeHEtest6
TH1F * h_pedestalwaver9_HE
TH1F * h_sumTSmeanAperLS7
TH2F * h_2DsumADCAmplLS1_LSselected
TH1F * h_sumADCAmplperLS8_P2
TH2F * h_mapDepth1Amplitude225_HE
TH1F * h_averNOSIGNALsumamplitude_HE
TH1F * h_sum0ADCAmplperLS6_M2
TH1F * h_pedestalaver9_HB
TH1F * h_sum0TSmaxAperLS7
TH2F * h_mapDepth4ADCAmpl225_HB
double mapRADDAM_HE[ndepth][neta][nphi]
TH2F * h_recSignalEnergy1_HB4
TH2F * h_2DsumADCAmplLS8_LSselected
TH1F * h_sumADCAmplperLS6_P1
TH2F * h_mapDepth3Pedw2_HE
TH1F * h_AmplitudeHBrest1
TH1F * h_sumamplitude_depth2_HE0
TH2F * h_mapCapCalib047_HE
TH1F * h_sumCutTSmeanAperLS4
double lsdep_estimator2_HBdepth1_
TH2F * h_mapTSmaxCalib_HE
TH2F * h_mapTSmeanCalib047_HO
TH2F * h_amplitudechannel1_HF3
TH2F * h_mapDepth3Error_HB
TH1F * h_Amplitude_notCapIdErrors_HO4
TH2F * h_mapGetRMSOverNormalizedSignal_HF
TH1F * h_sumCutAmplitudeperLS1
TH2F * h_amplitudechannel1_HB2
double calib0[ndepth][neta][nphi]
TH2F * h_amplitudechannel0_HB3
TH1F * h_sumamplitudechannel_HB
edm::EDGetTokenT< HBHERecHitCollection > tok_hbheNoise_
TH2F * h_mapDepth4Pedw0_HO
double lsdep_estimator2_HOdepth4_
float calibt[ndepth][neta][nphi]
TH1F * h_sumamplitude_depth1_HE
TH1F * h_sumCutTSmaxAperLS1
TH1F * h_sumADCAmplperLSdepth4HBu
TH2F * h_mapDepth1Pedw2_HE
TH1F * h_sumCutAmplperLS7
TH1F * h_sumCutADCAmplperLSdepth4HEu
TH2F * h_mapDepth3pedestalw_HB
TH2F * h_recNoiseEnergy2_HF1
TH1F * h_pedestalwaver4_HB
TH1F * h_ADCAmpl345Zoom1_HE
TH2F * h_recSignalEnergy0_HF4
double lsdep_estimator3_HBdepth2_
TH2F * h_recNoiseEnergy1_HB4
CMTRawAnalyzer(const edm::ParameterSet &)
TH2F * h_mapDepth3Ped3_HE
TH2F * h_recSignalEnergy0_HE6
TH2F * h_amplitudechannel1_HF4
TH1F * h_averSIGNALoccupancy_HB
TH2F * h_mapDepth2TSmaxA225_HE
TH2F * h_mapDepth2Amplitude225_HF
TH1F * h_sumTSmeanAperLS4
TH2F * h_mapDepth4Amplitude_HF
TH2F * h_mapDepth2Error_HF
void fillDigiAmplitudeHO(HODigiCollection::const_iterator &digiItr)
TH1F * h_sumADCAmplperLS8
TH1F * h_sumADCAmplEtaPhiLs
TH2F * h_mapDepth4ADCAmpl12_HF
TH1F * h_runbadrate_depth2_HE
TH1F * h_Eta16Depth3RADDAM_HED2P
TH2F * h_amplitudechannel0_HE7
TH2F * h_mapDepth2Pedw3_HF
TH2F * h_mapDepth4ADCAmpl225Copy_HB
double averSIGNALsumamplitude_HO
TH2F * h_2DsumADCAmplLSdepth3HFu
void fillDigiErrors(HBHEDigiCollection::const_iterator &digiItr)
TH2F * h_recNoiseEnergy1_HE6
double lsdep_estimator3_HBdepth1_
TH1F * h_numberofhitsHEtest
double sum0Estimator[nsub][ndepth][neta][nphi]
TH2F * h_mapAiForLS_bad_HE
TH1F * h_AmplitudeHEtest1
double lsdep_estimator1_HFdepth3_
TH2F * h_2DsumADCAmplLSdepth5HEu
TH1F * h_sumamplitude_depth4_HO
TH1F * h_pedestalaver9_HO
TH2F * h_2D0sumTSmeanALS7
TH2F * h_2DsumAmplitudeLS1
TH1F * h_sum0ADCAmplperLS7u
TH1F * h_runnbadchannels_depth3_HE
TH2F * h_amplitudechannel2_HF1
TH2F * h_mapDepth1ADCAmpl12_HF
TH1F * h_AmplitudeHBrest6
TH2F * h_mapDepth4AmplE34_HE
TH2F * h_mapDepth3ADCAmpl12_HF
TH1F * h_sigLayer1RADDAM6_HED2
double mapRADDAM_HED2[ndepth][neta]
TH2F * h_maprphinorm0_HE4
TH2F * h_mapDepth1Pedw1_HE
TH2F * h_2D0sumPedestalLS1
TH2F * h_mapAiForLS_good_HE
TH2F * h_2DsumADCAmplLSdepth4HBu
TH2F * h_mapDepth2ADCAmpl225_HE
TH1F * h_bcnnbadchannels_depth2_HF
TH2F * h_mapDepth1Ped1_HE
TH2F * h_mapDepth2ADCAmpl225_HB
TH2F * h_2DsumPedestalLS4
TH2F * h_recSignalEnergy1_HE2
TH2F * h_2D0sumADCAmplLS1
TH1F * h_A_Depth2RADDAM_HE
TH1F * h_sumPedestalperLS4
TH2F * h_mapADCCalib047_HE
bool getByToken(EDGetToken token, Handle< PROD > &result) const
TH2F * h_mapDepth4Pedw3_HO
TH1F * h_shape_good_channels_HE
TH1F * h_sumADCAmplperLS1_P1
TH1F * h_sum0AmplitudeperLS1
TH2F * h_amplitudechannel2_HE6
int calibcapiderror[ndepth][neta][nphi]
TH2F * h_mapDepth1Ped3_HE
float getSigma(int fCapId1, int fCapId2) const
get correlation element for capId1/2 = 0..3
TH1F * h_sumphiEta16Depth3RADDAM_HED2
TH2F * h_2D0sumPedestalLS6
TH2F * h_mapDepth1Error_HB
TH2F * h_mapDepth2Ped1_HE
bool studyTSmaxShapeHist_
int flagIterativeMethodCalibrationGroupDigi_
TH2F * h_mapRatioCalib_HO
TH1F * h_sumCutAmplperLS5
double maprphinorm[nsub][ndepth][neta][nphi]
TH1F * h_gsmdifferencefit4_HE
TH2F * h_recSignalEnergy2_HB2
TH2F * h_mapDepth6ADCAmpl12_HE
TH2F * h_mapADCAmpl11firstpeak_HE
TH1F * h_eventoccupancy_HF
#define DEFINE_FWK_MODULE(type)
TH2F * h_2DsumAmplitudeLS6
int mapRADDAM0_HE[ndepth][neta][nphi]
TH1F * h_sum0ErrorBperLS4
TH1F * h_runbadrate_depth1_HB
double lsdep_estimator3_HOdepth4_
TH2F * h_recNoiseEnergy2_HE5
TH2F * h_recNoiseEnergy2_HE6
TH1F * h_gsmdifferencefit5_HE
TH2F * h_2DAtaildepth2_HB
TH2F * h_mapDepth2Amplitude225_HE
int lsdep_cut3_max_HEdepth2_
TH1F * h_Amplitude_notCapIdErrors_HB1
TH1F * h_runnbadchannelsC_depth1_HE
TH2F * h_recNoiseEnergy0_HF3
double amplitudechannel[nsub][ndepth][neta][nphi]
double averSIGNALoccupancy_HO
TH1F * h_Amplitude_forCapIdErrors_HB1
TH2F * h_2DsumPedestalLS7
TH1F * h_runbadrateC_depth2_HB
TH2F * h_mapDepth3Error_HE
TH2F * h_mapCapCalib047_HO
double lsdep_estimator4_HFdepth2_
TH2F * h_mapCapCalib047_HB
TH2F * h_mapDepth2Ampl_HE
TH1F * h_pedestalwaver4_HO
TH2F * h_mapDepth3Amplitude225_HB
TH2F * h_mapDepth4Amplitude225_HB
TH2F * h_recSignalEnergy1_HF2
TH2F * h_amplitudechannel0_HE2
TH2F * h_mapDepth2Pedw3_HE
int lsdep_cut3_max_HEdepth1_
TH2F * h_mapDepth4Ped1_HO
TH1F * h_runnbadchannelsC_depth1_HF
TH2F * h_mapDepth4Ampl047_HF
TH1F * h_sumamplitude_depth4_HO0
TH1F * h_shape0_bad_channels_HB
TH1F * h_RatioOccupancy_HEP
TH2F * h_mapWidthCalib_HO
TH1F * h_Eta16Depth3RADDAM_HED2ALL
TH2F * h_amplitudechannel1_HE5
TH2F * h_mapDepth1TSmaxA_HF
TH2F * h_recNoiseEnergy0_HE6
TH1F * h_sigLayer1RADDAM6_HE
TH1F * h_shape0_good_channels_HB
double calibrTSmeanHOMax_
TH2F * h_2D0sumPedestalLS4
TH2F * h_2DsumPedestalLS1
TH1F * h_sigLayer1RADDAM0_HE
TH1F * h_sum0ErrorBperLS5
TH2F * h_mapDepth7ADCAmpl_HE
TH2F * h_2DsumADCAmplEtaPhiLs0
double lsdep_estimator2_HEdepth1_
TH2F * h_mapDepth4ADCAmpl_HE
TH1F * h_sumamplitude_depth2_HB
TH1F * h_NphiForEta16Depth3RADDAM_HED2ALL
TH2F * h_2DsumAmplitudeLS4
TH2F * h_recNoiseEnergy1_HE3
TH2F * h_mapDepth4ADCAmpl12_HE
double phi12(double phi1, double en1, double phi2, double en2)
TH1F * h_energyhitNoise_HE
double lsdep_estimator1_HBdepth1_
bool studyRatioShapeHist_
edm::EDGetTokenT< HBHEDigiCollection > tok_hbhe_
TH2F * h_mapDepth4Error_HB
std::vector< T >::const_iterator const_iterator
TH2F * h_2DsumPedestalLS3
TH2F * h_FullSignal3D0_HO
TH1F * h_sumTSmeanAperLS3
TH2F * h_2D0Ataildepth1_HB
TH2F * h_mapDepth4Ampl_HF
int lsdep_cut3_max_HEdepth3_
TH1F * h_sum0ADCAmplperLSdepth5HEu
TH2F * h_mapDepth3ADCAmpl_HF
TH1F * h_gsmdifferencefit3_HE
const HcalQIEShape * shape
TH1F * h_sumamplitude_depth2_HE
double lsdep_estimator1_HFdepth4_
TH2F * h_mapDepth2linADCAmpl12_HE
int bunchCrossing() const
TH2F * h_mapGetRMSOverNormalizedSignal0_HO
TH2F * h_recNoiseEnergy1_HF4
TH2F * h_mapDepth1Ampl047_HB
TH2F * h_mapDepth2Amplitude225_HB
TH2F * h_mapTSmeanCalib047_HB
TH2F * h_recNoiseEnergy0_HB4
TH2F * h_mapDepth4Amplitude_HO
TH2F * h_recNoiseEnergy0_HE2
edm::EDGetTokenT< HFRecHitCollection > tok_hfSignal_
TH2F * h_mapDepth1AmplE34_HF
TH2F * h_mapDepth4pedestal_HB
double recSignalEnergy2[nsub][ndepth][neta][nphi]
TH2F * h_recSignalEnergy1_HB2
TH2F * h_amplitudechannel2_HE7
TH2F * h_mapDepth3Ampl047_HE
edm::LuminosityBlockNumber_t luminosityBlock() const
TH2F * h_mapDepth1ADCAmpl225_HE
TH2F * h_mapDepth7ADCAmpl225_HE
double averSIGNALoccupancy_HB
TH2F * h_mapDepth1Ped0_HB
TH2F * h_2D0sumAmplitudeLS2
TH2F * h_mapDepth1Ped2_HB
TH1F * h_nbadchannels_depth3_HE
TH1F * h_sum0ADCAmplperLS6_P2
TH1F * h_Amplitude_forCapIdErrors_HE3
TH2F * h_mapDepth4TSmeanA225_HB
TH1F * h_sumamplitude_depth2_HB0
TH1F * h_sumAmplperLS1_LSselected
TH1F * h_corrforxaMAIN0_HE
TH2F * h_mapDepth2ADCAmpl_HF
TH1F * h_sumamplitude_depth1_HF1
int lsdep_cut1_peak_HEdepth3_
TH2F * h_mapGetRMSOverNormalizedSignal_HE
double lsdep_estimator1_HEdepth4_
TH1F * h_Amplitude_forCapIdErrors_HF1
TH2F * h_2DsumAmplitudeLS3
T * make(const Args &...args) const
make new ROOT object
TH1F * h_sum0PedestalperLS1
TH1F * h_pedestalaver4_HO
TH2F * h_mapDepth2TSmaxA_HB
TH2F * h_FullSignal3D0_HB
TH1F * h_sum0ADCAmplperLSdepth3HBu
TH2F * h_mapDepth3Amplitude_HE
TH1F * h_ADCAmpl345Zoom_HE
TH2F * h_mapDepth5TSmeanA225_HE
TH2F * h_recNoiseEnergy0_HE7
TH2F * h_recNoiseEnergy2_HB1
TH2F * h_recNoiseEnergy0_HB2
TH2F * h_mapTSmaxCalib047_HB
TH1F * h_Amplitude_notCapIdErrors_HF2
TH1F * h_sumADCAmplperLS8_P1
TH2F * h_mapDepth1AmplE34_HB
TH2F * h_recSignalEnergy2_HE7
TH1F * h_sum0TSmeanAperLS7
TH2F * h_mapTSmaxCalib_HO
double lsdep_estimator2_HEdepth2_
TH2F * h_mapDepth4Error_HO
double sumEstimator3[nsub][ndepth][neta][nphi]
TH2F * h_recSignalEnergy0_HE2
TH2F * h_mapDepth2ADCAmpl_HB
TH2F * h_mapDepth1TSmaxA225_HF
TH2F * h_recNoiseEnergy1_HB3
double lsdep_estimator4_HBdepth2_
TH2F * h_mapDepth7TSmeanA225_HE
TH2F * h_mapRatioCalib_HB
double lsdep_estimator1_HEdepth3_
TH2F * h_mapDepth1ADCAmpl225_HF
TH1F * h_Amplitude_forCapIdErrors_HE1
int lsdep_cut1_peak_HFdepth1_
TH2F * h_mapDepth4Amplitude225_HF
TH2F * h_mapDepth2pedestal_HB
double lsdep_estimator3_HEdepth2_
TH2F * h_mapGetRMSOverNormalizedSignal_HO
TH2F * h_mapDepth1Pedw3_HE
bool studyTSmeanShapeHist_
void fillDigiAmplitudeQIE11(QIE11DataFrame qie11df)
TH2F * h_2D0sumADCAmplLSdepth4HFu
double calibrTSmeanHFMin_
TH1F * h_sumamplitude_depth1_HB
TH1F * h_corrforxaADDI_HE
TH1F * h_sumPedestalperLS1
TH1F * h_sumAmplitudeperLS1
constexpr int dccid() const
get the (Hcal local) DCC id for VME, crate number for uTCA
void fillDigiErrorsHF(HFDigiCollection::const_iterator &digiItr)
TH2F * h_mapDepth3ADCAmpl225_HB
TH1F * h_nbadchannels_depth2_HB
TH2F * h_2DAtaildepth1_HB
double lsdep_estimator5_HEdepth1_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbheSignal_
TH1F * h_shape_good_channels_HF
float getValue(int fCapId) const
get value for capId = 0..3
TH2F * h_mapDepth2Pedw1_HE
TH1F * h_sum0TSmaxAperLS4
TH2F * h_amplitudechannel2_HF3
double lsdep_estimator1_HFdepth2_
TH2F * h_mapDepth4TSmaxA225_HF
TH2F * h_recNoiseEnergy1_HE1
int lsdep_cut3_max_HOdepth4_
TH1F * h_sigLayer2RADDAM6_HE
TH2F * h_recSignalEnergy2_HE6
TH2F * h_mapDepth4TSmeanA_HB
TH2F * h_mapDepth1RADDAM0_HE
TH2F * h_2DsumADCAmplLSdepth3HBu
TH2F * h_2D0sumADCAmplLSdepth6HEu
TH2F * h_recSignalEnergy1_HE1
TH1F * h_sum0ADCAmplperLS8_P2
TH1F * h_nbadchannels_depth1_HB
TH2F * h_recNoiseEnergy1_HB2
TH2F * h_mapDepth3Ampl_HB
TH1F * h_sum0ADCAmplperLS3_P2
TH2F * h_mapDepth4ADCAmpl_HB
TH1F * h_energyhitSignal_HB
TH1F * h_runnbadchannelsC_depth2_HE
TH1F * h_bcnvsamplitude_HB
TH2F * h_amplitudechannel0_HF2
TH2F * h_mapDepth5ADCAmpl225_HE
TH1F * h_sumADCAmplperLS3_M2
TH1F * h_sum0ADCAmplperLS1_P2
TH2F * h_mapDepth1TSmeanA225_HE
TH2F * h_recSignalEnergy1_HE3
int flagfitshunt1pedorledlowintensity_
TH2F * h_mapDepth2Pedw3_HB
TH1F * h_sumamplitude_depth3_HE
TH1F * h_sumCutTSmeanAperLS8
TH2F * h_mapDepth1TSmeanA_HB
TH2F * h_mapDepth4pedestal_HE
TH1F * h_Amplitude_forCapIdErrors_HB2
TH2F * h_mapDepth4Ampl047_HB
TH1F * h_sumADCAmplperLS2
TH2F * h_mapDepth1ADCAmpl12_HE
TH2F * h_2DsumADCAmplEtaPhiLs1
TH2F * h_2DsumADCAmplLS4_LSselected
TH1F * h_nbadchannels_depth2_HE
TH1F * h_sumADCAmplLS1copy1
TH2F * h_mapDepth1Ampl_HE
double dPhiWsign(double phi1, double phi2)
TH1F * h_sum0ADCAmplperLS3_M1
TH2F * h_recNoiseEnergy0_HB3
double tocamplchannel[nsub][ndepth][neta][nphi]
double calib2[ndepth][neta][nphi]
TH2F * h_mapDepth4ADCAmpl225_HF
TH1F * h_sum0ErrorBperLS3
TH1F * h_sum0PedestalperLS8
TH1F * h_sumADCAmplperLS6u
TH1F * h_runbadrate_depth4_HO
double lsdep_estimator2_HFdepth2_
TH2F * h_mapDepth2Ped2_HB
TH2F * h_mapDepth1Ampl_HB
TH2F * h_maprphinorm2_HE7
double averNOSIGNALsumamplitude_HO
TH2F * h_mapDepth4Ampl047_HO
TH2F * h_2DsumPedestalLS5
TH1F * h_Amplitude_notCapIdErrors_HE3
TH2F * h_mapDepth2Pedw1_HB
TH2F * h_mapDepth2Ampl047_HB
TH2F * h_mapDepth1TSmaxA_HB
TH2F * h_mapDepth2Ampl_HF
TH2F * h_mapDepth2pedestalw_HF
TH2F * h_mapADCAmpl11firstpeak0_HE
TH1F * h_sumCutAmplitudeperLS7
edm::Service< TFileService > fs_
TH2F * h_2D0sumTSmeanALS8
TH2F * h_mapCapCalib047_HF
TH1F * h_sum0AmplitudeperLS2
TH1F * h_bcnvsamplitude0_HF
double lsdep_estimator1_HEdepth5_
double calibrTSmeanHEMax_
TH2F * h_mapDepth7ADCAmpl225Copy_HE
TH1F * h_sumCutTSmaxAperLS2
TH1F * h_sumCutAmplperLS4
TH1F * h_sumCutTSmeanAperLS3
TH2F * h_mapDepth2ADCAmplSiPM_HE
TH2F * h_mapDepth3pedestal_HB
TH2F * h_amplitudechannel1_HB1
TH2F * h_mapDepth4ADCAmpl225Copy_HO
TH2F * h_mapenophinorm_HE2
TH2F * h_2DsumADCAmplLSdepth4HEu
TH1F * h_sum0PedestalperLS6
TH2F * h_mapDepth4pedestalw_HF
TH1F * h_sum0ADCAmplperLS3_M2
TH2F * h_amplitudechannel0_HB1
int lsdep_cut1_peak_HBdepth1_
TH2F * h_recSignalEnergy2_HF2
TH1F * h_sumCutTSmeanAperLS7
TH2F * h_recSignalEnergy0_HF3
TH2F * h_amplitudechannel1_HE4
TH2F * h_recSignalEnergy1_HE7
TH1F * h_sumADCAmplperLS6_P2
TH2F * h_mapADCAmpl11secondpeak_HE
TH1F * h_nbadchannels_depth4_HO
TH2F * h_mapWidthCalib047_HB
TH2F * h_amplitudechannel2_HE5
TH1F * h_sumPedestalperLS3
TH1F * h_sumADCAmplLS1copy2
TH1F * h_sumCutADCAmplperLS5
TH1F * h_runbadrate0_depth2_HB
int lsdep_cut1_peak_HEdepth2_
double sumEstimator5[nsub][ndepth][neta][nphi]
TH1F * h_sumADCAmplperLSdepth6HEu
TH2F * h_recSignalEnergy2_HE3
TH1F * h_runbadrate_depth1_HF
TH1F * h_sum0ErrorBperLS2
TH2F * h_mapDepth4Ampl_HE
TH2F * h_mapDepth4TSmaxA_HB
TH1F * h_sum0ADCAmplperLS6_M1
TH2F * h_mapDepth4Error_HE
TH2F * h_2D0sumADCAmplLS3
TH1F * h_sumCutADCAmplperLSdepth7HEu
TH1F * h_pedestalwaver9_HF
TH1F * h_sumTSmeanAperLS8
TH2F * h_recSignalEnergy1_HE6
double lsdep_estimator2_HFdepth1_
edm::EDGetTokenT< HFRecHitCollection > tok_hfNoise_
bool getData(T &iHolder) const
TH2F * h_mapDepth1TSmeanA225_HF
TH2F * h_mapDepth5AmplE34_HE
double lsdep_estimator1_HBdepth3_
TH2F * h_mapAiForLS_good_HF
TH2F * h_map0AiForLS_bad_HB
TH2F * h_mapWidthCalib047_HO
void fillDigiErrorsQIE11(QIE11DataFrame qie11df)
TH1F * h_sigLayer1RADDAM5_HED2
TH2F * h_mapDepth1ADCAmpl225_HB
TH2F * h_mapDepth2Error_HE
TH1F * h_sumamplitude_depth2_HF
TH2F * h_amplitudechannel2_HB4
TH2F * h_mapenophinorm2_HE1
const HcalPedestalWidth * getPedestalWidth(const HcalGenericDetId &fId) const
TH2F * h_mapDepth2TSmeanA225_HE
TH2F * h_mapDepth4TSmaxA_HO
TH1F * h_sumCutAmplperLS6
TH2F * h_2D0sumTSmeanALS5
TH1F * h_shape0_bad_channels_HF
TH2F * h_amplitudechannel0_HE6
TH1F * h_sum0TSmaxAperLS3
constexpr DetId detid() const
Get the detector id.
TH2F * h_mapDepth1Pedw3_HB
TH1F * h_pedestalwaver9_HO
constexpr HcalSubdetector subdet() const
get the subdetector
TH2F * h_mapDepth1TSmaxA_HE
TH2F * h_mapDepth3ADCAmpl12_HE
const HcalRespCorr * getHcalRespCorr(const HcalGenericDetId &fId) const
void adc2fC(const HBHEDataFrame &df, CaloSamples &lf) const override
TH2F * h_maprphinorm0_HE2
TH2F * h_mapDepth4AmplE34_HB
TH2F * h_mapDepth2Ped0_HE
TH2F * h_mapDepth1pedestal_HB
TH1F * h_Amplitude_forCapIdErrors_HE2
TH1F * h_bcnvsamplitude_HO
TH2F * h_mapDepth3TSmaxA225_HB
TH1F * h_sumamplitudechannel_HE
double averNOSIGNALoccupancy_HB
constexpr int htrTopBottom() const
get the htr top/bottom (1=top/0=bottom), valid for VME
TH2F * h_mapDepth4TSmeanA_HF
TH2F * h_mapADCAmplfirstpeak_HE
TH1F * h_sum0TSmaxAperLS2
TH2F * h_mapDepth4ADCAmpl225Copy_HF
TH1F * h_sum0PedestalperLS4
TH2F * h_recNoiseEnergy0_HE5
TH2F * h_mapDepth7TSmeanA_HE
TH2F * h_mapDepth4Ampl_HO
TH1F * h_nbadchannels_depth2_HF
TH1F * h_sum0AmplitudeperLS5
TH2F * h_mapDepth3ADCAmpl225_HF
TH2F * h_amplitudechannel0_HF4
const HcalLUTCorr * getHcalLUTCorr(const HcalGenericDetId &fId) const
TH1F * h_orbitNumvsamplitude0_HE
TH1F * h_sigLayer2RADDAM5_HED2
TH2F * h_recSignalEnergy0_HB2
TH1F * h_sumADCAmplperLS4
TH1F * h_eventoccupancy_HB
const HcalTimeCorr * getHcalTimeCorr(const HcalGenericDetId &fId) const
TH2F * h_2D0sumAmplitudeLS1
TH1F * h_sumADCAmplperLS3_M1
TH2F * h_recSignalEnergy2_HE4
TH1F * h_shape0_good_channels_HE
TH2F * h_mapADCCalib047_HO
TH1F * h_eventoccupancy_HE
double lsdep_estimator1_HOdepth4_
TH2F * h_mapDepth2Ped2_HE
TH1F * h_sumADCAmplEtaPhiLs_ietaphi
double lsdep_estimator4_HOdepth4_
TH1F * h_mapDepth3RADDAM16_HE
TH1F * h_sum0ADCAmplperLSdepth4HBu
TH1F * h_shape_bad_channels_HO
TH1F * h_Amplitude_forCapIdErrors_HF2
TH1F * h_RatioOccupancy_HBP
double lsdep_estimator4_HBdepth1_
TH1F * h_sumCutADCAmplperLS6
TH2F * h_2D0sumPedestalLS8
TH2F * h_recSignalEnergy0_HB1
TH1F * h_runbadrateC_depth2_HE
TH2F * h_recSignalEnergy1_HE4
TH1F * h_sumADCAmplperLS8_M2
TH1F * h_sumamplitudechannel_HF
TH2F * h_mapDepth2ADCAmpl12SiPM_HE
int lsdep_cut3_max_HFdepth1_
TH2F * h_mapDepth1pedestalw_HB
TH2F * h_mapDepth2RADDAM_HE
TH2F * h_mapDepth6pedestal_HE
TH2F * h_mapDepth3TSmeanA225_HF
TH2F * h_mapDepth4ADCAmpl12_HB
TH2F * h_2D0sumPedestalLS7
TH2F * h_mapDepth1Amplitude_HF
TH1F * h_sumADCAmplLS1copy5
TH2F * h_mapDepth2Ped1_HB
TH2F * h_2D0sumADCAmplLSdepth3HFu
const HcalTopology * topo
TH2F * h_amplitudechannel2_HF4
TH1F * h_runbadrateC_depth2_HF
double forallestimators_amplitude_bigger_
TH1F * h_sum0AmplitudeperLS8
TH2F * h_mapDepth5Ampl047_HE
TH1F * h_sum0AmplitudeperLS3
TH1F * h_ADCAmpl345Zoom1_HB
int lsdep_cut3_max_HBdepth2_
TH1F * h_sumADCAmplperLS7
double lsdep_estimator5_HFdepth2_
TH1F * h_runbadrate_depth2_HF
TH1F * h_sumTSmeanAperLS1
TH2F * h_mapDepth3Amplitude_HB
TH1F * h_energyhitNoise_HB
TH1F * h_sumADCAmplperLS1_M1
constexpr int iphi() const
get the cell iphi
TH1F * h_totalAmplitudeHEperEvent
TH2F * h_mapDepth3RADDAM_HE
TH2F * h_mapTSmaxCalib_HF
TH2F * h_mapDepth3ADCAmpl225Copy_HE
TH2F * h_2DsumPedestalLS2
TH1F * h_sumADCAmplperLSdepth3HBu
TH2F * h_mapDepth2Ped3_HF
TH1F * h_energyhitSignal_HF
TH1F * h_sumADCAmplperLS6
TH2F * h_mapDepth5pedestalw_HE
TH1F * h_runbadrateC_depth1_HF
int lsdep_cut3_max_HFdepth2_
TH2F * h_mapWidthCalib047_HF
TH2F * h_recSignalEnergy2_HE2
TH2F * h_amplitudechannel0_HE1
TH1F * h_runnbadchannelsC_depth2_HF
int getRBX(int &i, int &j, int &k)
TH1F * h_sumCutAmplitudeperLS8
TH2F * h_mapAiForLS_bad_HF
TH2F * h_mapDepth7Ampl047_HE
TH1F * h_sumCutTSmeanAperLS1
TH1F * h_sum0ADCAmplperLS1
TH1F * h_sumamplitude_depth2_HE1
TH1F * h_sumADCAmplperLS8_M1
TH2F * h_maprphinorm2_HE5
TH2F * h_mapDepth2TSmaxA_HE
TH1F * h_bcnnbadchannels_depth1_HB
TH1F * h_sumADCAmplperLS1
TH2F * h_mapDepth2Amplitude_HF
TH2F * h_mapDepth1Ampl_HF
TH2F * h_mapADCAmpl12secondpeak_HE
TH2F * h_mapDepth1TSmeanA225_HB
TH1F * h_sum0ErrorBperLS1
TH2F * h_recNoiseEnergy0_HE4
TH2F * h_amplitudechannel2_HB3
TH1F * h_runnbadchannels_depth2_HB
TH2F * h_mapDepth3Amplitude_HF
TH1F * h_sumTSmeanAperLS2
const edm::ESGetToken< HcalTopology, HcalRecNumberingRecord > tokTopo_
TH2F * h_mapDepth4Amplitude225_HO
TH1F * h_sumamplitude_depth1_HB0
constexpr int htrSlot() const
get the htr slot
TH2F * h_2D0sumTSmeanALS2
TH1F * h_sumamplitude_depth3_HE1
TH1F * h_corrforxaADDI0_HE
TH2F * h_mapDepth5ADCAmpl12_HE
TH1F * h_eventamplitude_HO
double lsdep_estimator4_HFdepth1_
TH2F * h_mapDepth5ADCAmpl_HE
TH1F * h_sum0ADCAmplperLS8_P1
TH1F * h_totalAmplitudeHOperEvent
TH2F * h_mapDepth1Ampl047_HE
TH2F * h_mapDepth2Ped2_HF
TH1F * h_sumCutADCAmplperLS2
TH2F * h_mapDepth3Pedw3_HE
TH2F * h_mapDepth4Ampl_HB
TH1F * h_sumAmplitudeperLS6
TH1F * h_AamplitudewithPedSubtr_RADDAM_HEzoom0
TH2F * h_maprphinorm0_HE3
TH2F * h_mapDepth6Amplitude_HE
TH2F * h_mapenophinorm_HE5
TH2F * h_map0AiForLS_bad_HF
TH2F * h_recSignalEnergy2_HF3
TH1F * h_sumCutTSmeanAperLS5
TH1F * h_sum0ADCAmplperLS2
TH1F * h_shape_good_channels_HO
double averNOSIGNALsumamplitude_HB
TH2F * h_mapDepth1Pedw1_HB
bool usePedestalSubtraction_
TH2F * h_mapDepth2Pedw2_HF
TH2F * h_mapTSmeanCalib047_HF
TH1F * h_totalAmplitudeHO
TH1F * h_numberofhitsHFtest
TH2F * h_mapTSmaxCalib047_HE
TH1F * h_runbadrate0_depth1_HF
TH2F * h_mapDepth2Pedw0_HF
TH2F * h_maprphinorm2_HE1
TH2F * h_amplitudechannel0_HF1
double calibrTSmeanHOMin_
TH1F * h_NphiForEta16Depth3RADDAM_HED2P
TH1F * h_sigLayer2RADDAM_HE
TH2F * h_mapDepth4ADCAmpl_HO
TH1F * h_sum0ADCAmplperLS6_P1
TH2F * h_recNoiseEnergy0_HB1
TH2F * h_amplitudechannel0_HF3
TH2F * h_mapDepth2ADCAmpl225Copy_HB
TH2F * h_mapDepth7pedestalw_HE
TH1F * h_sumCutADCAmplperLSdepth3HBu
TH2F * h_amplitudechannel1_HF1
TH1F * h_ADCAmpl_HBCapIdNoError
void analyze(const edm::Event &, const edm::EventSetup &) override
double caliba[ndepth][neta][nphi]
TH2F * h_2D0sumAmplitudeLS8
TH1F * h_sumCutADCAmplperLS3
double averNOSIGNALoccupancy_HE
TH2F * h_recSignalEnergy1_HB3
TH2F * h_recSignalEnergy0_HB3
TH2F * h_mapDepth6ADCAmpl225_HE
bool studyCalibCellsHist_
TH2F * h_recSignalEnergy1_HB1
TH2F * h_recNoiseEnergy1_HF1
Abs< T >::type abs(const T &t)
TH1F * h_Amplitude_notCapIdErrors_HB2
TH1F * h_energyhitSignal_HE
TH1F * h_runbadrateC_depth1_HE
TH1F * h_orbitNumvsamplitude0_HF
double sumEstimator1[nsub][ndepth][neta][nphi]
TH2F * h_2DsumAmplitudeLS8
TH2F * h_mapDepth3Amplitude225_HF
TH2F * h_mapDepth3TSmeanA225_HE
constexpr int ieta() const
get the cell ieta
static const float adc2fC[128]
TH1F * h_sumADCAmplperLS1_P2
TH2F * h_mapDepth2Ampl_HB
TH1F * h_runbadrate0_depth2_HF
TH2F * h_mapDepth2Amplitude_HE
constexpr int readoutVMECrateId() const
get the readout VME crate number
TH1F * h_sum0ADCAmplperLS6
double calibw[ndepth][neta][nphi]
TH1F * h_nevents_per_eachLS
TH2F * h_recSignalEnergy2_HF4
TH2F * h_recSignalEnergy0_HF2
edm::EDGetTokenT< QIE11DigiCollection > tok_qie11_
TH2F * h_2DsumADCAmplEtaPhiLs20
TH1F * h_ADCAmpl345_HBCapIdNoError
float getValue(int fCapId) const
get value for capId = 0..3
TH2F * h_amplitudechannel1_HE6
TH2F * h_amplitudechannel2_HE3
TH2F * h_mapDepth7Amplitude225_HE
TH2F * h_mapDepth3ADCAmpl12_HB
int lsdep_cut1_peak_HBdepth2_
TH2F * h_mapDepth4Ped0_HO
TH1F * h_sum0TSmeanAperLS6
TH1F * h_sum0TSmaxAperLS1
double lsdep_estimator1_HBdepth2_
TH2F * h_2D0sumADCAmplLSdepth3HBu
TH1F * h_runbadrate_depth1_HE
TH1F * h_sumADCAmplperLS6_M1
TH2F * h_recNoiseEnergy2_HF3
constexpr int spigot() const
get the spigot (input number on DCC, AMC card number for uTCA)
TH2F * h_mapTSmaxCalib_HB
TH2F * h_amplitudechannel2_HE1
TH2F * h_2D0sumTSmeanALS4
TH2F * h_mapDepth2TSmeanA225_HB
TH2F * h_mapDepth1ADCAmpl12_HB
TH2F * h_mapDepth1Pedw2_HF
double averNOSIGNALsumamplitude_HF
TH2F * h_mapDepth5Amplitude225_HE
TH2F * h_mapDepth2ADCAmpl12_HB
TH2F * h_map0AiForLS_good_HF
TH2F * h_2D0sumADCAmplLS5
TH2F * h_mapenophinorm2_HE7
TH2F * h_mapDepth2TSmeanA225_HF
TH2F * h_mapDepth5TSmaxA225_HE
TH1F * h_averNOSIGNALoccupancy_HO
TH1F * h_runnbadchannelsC_depth2_HB
TH1F * h_sumphiEta16Depth3RADDAM_HED2ALL
TH2F * h_recNoiseEnergy2_HE7
TH1F * h_sum0ADCAmplperLS8
TH2F * h_mapenophinorm_HE1
TH1F * h_sumAmplitudeperLS1_LSselected
TH2F * h_mapDepth1Pedw0_HF
TH2F * h_mapDepth2TSmaxA_HF
TH1F * h_sumCutAmplitudeperLS3
double lsdep_estimator5_HBdepth2_
TH1F * h_sumamplitude_depth3_HE0
TH1F * h_shape0_Ahigh_HB2
bool usecontinuousnumbering_
TH1F * h_AmplitudeHBtest1
TH2F * h_mapDepth3pedestalw_HE
TH1F * h_averNOSIGNALoccupancy_HB
TH1F * h_ADCAmpl_HBCapIdError
TH2F * h_2DsumADCAmplLS3_LSselected
TH1F * h_runnbadchannels_depth2_HE
TH1F * h_sumADCAmplperLSdepth5HEu
TH2F * h_mapDepth2Ped3_HB
TH1F * h_averSIGNALoccupancy_HO
TH1F * h_totalAmplitudeHE
TH2F * h_maprphinorm2_HE2
TH2F * h_mapDepth7TSmaxA225_HE
TH2F * h_mapDepth3Ped2_HE
TH2F * h_mapDepth1Ped3_HF
TH1F * h_sumADCAmplLS1copy4
TH2F * h_mapDepth7pedestal_HE
TH1F * h_sumCutADCAmplperLS8
TH1F * h_lsnumber_per_eachLS
TH2F * h_mapADCCalib047_HF
double amplitudechannel0[nsub][ndepth][neta][nphi]
TH2F * h_recNoiseEnergy1_HF2
TH1F * h_Amplitude_forCapIdErrors_HO4
int flagestimatornormalization_
TH1F * h_orbitNumvsamplitude0_HO
TH2F * h_mapRatioCalib_HE
TH2F * h_mapADCAmpl12firstpeak_HE
TH2F * h_mapADCAmplfirstpeak0_HE
TH1F * h_sum0TSmaxAperLS6
TH1F * h_pedestalwaver4_HE
double recSignalEnergy0[nsub][ndepth][neta][nphi]
TH2F * h_recNoiseEnergy2_HE1
TH2F * h_recNoiseEnergy1_HF3
TH2F * h_mapDepth2TSmaxA225_HB
TH1F * h_sumAmplitudeperLS3
TH2F * h_mapDepth4TSmaxA225_HB
TH2F * h_mapDepth6TSmeanA_HE
TH1F * h_sumCutADCAmplperLS6u
TH2F * h_map0AiForLS_bad_HO
double recSignalEnergy1[nsub][ndepth][neta][nphi]
TH1F * h_sumADCAmplperLS1_M2
TH1F * h_sumTSmeanAperLS1_LSselected
double lsdep_estimator1_HBdepth4_
TH1F * h_sumADCAmplperLS3_P1
TH1F * h_sumADCAmplLS1copy3
TH2F * h_mapDepth2Pedw2_HB
TH1F * h_sum0ADCAmplperLSdepth7HEu
TH1F * h_runbadrate_depth3_HE
TH2F * h_mapAiForLS_good_HB
TH1F * h_runbadrate0_depth2_HE
TH2F * h_mapDepth5Ampl_HE
double averNOSIGNALoccupancy_HO
TH1F * h_sigLayer1RADDAM_HE
TH2F * h_mapADCAmpl12secondpeak0_HE
TH1F * h_sumCutADCAmplperLSdepth4HBu
TH1F * h_sumPedestalperLS7
TH2F * h_mapDepth5Error_HE
TH1F * h_shape_bad_channels_HF
int flagtodefinebadchannel_
double recNoiseEnergy0[nsub][ndepth][neta][nphi]
TH1F * h_shape0_Ahigh_HB3
TH2F * h_mapDepth4TS012_HO
TH1F * h_sum0ADCAmplperLS7
TH2F * h_mapDepth3ADCAmpl225_HE
TH1F * h_sumCutAmplperLS1
HcalLogicalMap createMap(const HcalTopology *topo, unsigned int mapIOV=4)
TH2F * h_mapDepth4ADCAmpl12_HO
TH2F * h_mapDepth1AmplE34_HE
TH2F * h_mapAiForLS_good_HO
TH1F * h_bcnbadrate0_depth2_HF
TH2F * h_mapDepth4AmplE34_HF
TH1F * h_sum0PedestalperLS5
TH2F * h_2D0sumAmplitudeLS5
TH2F * h_mapDepth2TSmaxA225_HF
TH2F * h_mapDepth3Ampl047_HB
TH1F * h_sumCutTSmeanAperLS2
TH2F * h_maprphinorm0_HE1
TH1F * h_sumamplitude_depth2_HB1
TH2F * h_mapADCAmplsecondpeak_HE
TH2F * h_mapDepth3pedestalw_HF
TH2F * h_mapDepth2TSmeanA_HE
TH1F * h_sumamplitude_depth1_HE1
double lsdep_estimator4_HEdepth2_
TH2F * h_amplitudechannel2_HB1
TH2F * h_mapDepth1ADCAmpl_HB
TH1F * h_bcnbadrate0_depth2_HB
TH2F * h_2DsumADCAmplEtaPhiLs10
TH2F * h_recSignalEnergy2_HB3
TH1F * h_runnbadchannelsC_depth3_HE
TH2F * h_mapDepth6pedestalw_HE
TH2F * h_mapDepth3TSmeanA_HB
const HcalFrontEndId getHcalFrontEndId(const DetId &)
double sumEstimator0[nsub][ndepth][neta][nphi]
TH2F * h_2DsumAmplitudeLS2
TH2F * h_amplitudechannel1_HE2
TH2F * h_mapDepth3Ped1_HE
TH2F * h_map0AiForLS_good_HO
int lsdep_cut1_peak_HOdepth4_
TH2F * h_2D0sumPedestalLS5
double lsdep_estimator5_HBdepth1_
TH2F * h_mapADCAmpl12firstpeak0_HE
TH2F * h_mapDepth1Error_HE
double sumEstimator2[nsub][ndepth][neta][nphi]
TH2F * h_mapTSmaxCalib047_HF
TH1F * h_sum0TSmeanAperLS4
TH2F * h_mapDepth3pedestal_HF
TH2F * h_mapDepth3ADCAmpl_HB
TH2F * h_recSignalEnergy2_HF1
TH1F * h_sumADCAmplEtaPhiLs_bbbc
TH1F * h_pedestalwaver4_HF
TH2F * h_recNoiseEnergy2_HE3
TH2F * h_mapADCCalib047_HB
TH2F * h_mapDepth1Amplitude_HE
TH1F * h_sum0TSmaxAperLS8
TH2F * h_mapDepth1ADCAmpl225Copy_HB
TH2F * h_mapDepth2AmplE34_HB
double lsdep_estimator3_HEdepth1_
constexpr int fiberIndex() const
get the fiber index. For VME 1-8 (which of eight fibers carried by a spigot), for uTCA fibers are zer...
TH2F * h_mapDepth1Pedw0_HB
TH2F * h_2D0sumPedestalLS3
TH1F * h_sumAmplitudeperLS4
TH2F * h_2DsumADCAmplEtaPhiLs00
TH2F * h_2D0sumADCAmplLS8
TH2F * h_mapDepth1Ped2_HF
TH1F * h_energyhitNoise_HF
TH1F * h_sumADCAmplperLS5
TH2F * h_mapDepth3Error_HF
TH2F * h_amplitudechannel2_HB2
TH2F * h_mapDepth3TSmaxA_HF
TH2F * h_2D0sumADCAmplLS2
TH2F * h_maprphinorm0_HE7
TH2F * h_maprphinorm0_HE5
TH1F * h_runbadrate_depth2_HB
TH2F * h_mapDepth2Ped0_HB
TH2F * h_2D0sumPedestalLS2
TH1F * h_Amplitude_notCapIdErrors_HE1
TH2F * h_mapADCAmpl11secondpeak0_HE
TH2F * h_recNoiseEnergy0_HE3
TH1F * h_sum0AmplitudeperLS4
constexpr int samples() const
total number of samples in the digi
TH2F * h_mapDepth4TSmaxA_HF
TH1F * h_pedestalaver9_HE
double lsdep_estimator5_HFdepth1_
TH2F * h_maprphinorm2_HE6
TH2F * h_2DsumADCAmplEtaPhiLs2
TH2F * h_recNoiseEnergy2_HE4
TH2F * h_mapDepth1Pedw1_HF
TH1F * h_eventamplitude_HB
TH1F * h_pedestalwaver9_HB
TH1F * h_sumamplitude_depth2_HF0
TH2F * h_recSignalEnergy0_HE4
~CMTRawAnalyzer() override
TH1F * h_shape0_Ahigh_HB0
TH2F * h_mapenophinorm2_HE5
TH2F * h_mapDepth4Amplitude_HE
TH1F * h_sum0ADCAmplperLS3
TH2F * h_2D0sumAmplitudeLS6
void fillDigiAmplitudeHF(HFDigiCollection::const_iterator &digiItr)
TH2F * h_mapDepth3Amplitude225_HE
TH1F * h_sumamplitude_depth1_HF0
TH1F * h_sum0TSmeanAperLS3
TH1F * h_sumCutADCAmplperLS4
TH2F * h_map0AiForLS_bad_HE
TH1F * h_sumCutADCAmplperLS7u
double lsdep_estimator5_HOdepth4_
TH2F * h_recNoiseEnergy2_HF4
TH1F * h_sumAmplitudeperLS5
TH2F * h_recSignalEnergy0_HE1
TH2F * h_mapDepth2ADCAmpl225Copy_HE
TH2F * h_mapenophinorm_HE4
TH2F * h_mapDepth1Pedw3_HF
TH2F * h_mapDepth2TSmeanA_HB
TH2F * h_mapTSmeanCalib_HB
int flagIterativeMethodCalibrationGroupReco_
TH2F * h_mapDepth7AmplE34_HE
TH2F * h_amplitudechannel1_HB3
void beginRun(const edm::Run &r, const edm::EventSetup &iSetup) override
TH2F * h_mapWidthCalib_HF
TH2F * h_mapenophinorm_HE6
TH1F * h_nbadchannels_depth1_HE
TH2F * h_mapDepth3TSmeanA_HF
TH1F * h_sumPedestalperLS2
TH2F * h_mapTSmeanCalib047_HE
TH1F * h_sumTSmaxAperLS1_LSselected
TH1F * h_runbadrateC_depth4_HO
TH1F * h_sumamplitude_depth4_HO1
TH1F * h_runnbadchannels_depth1_HB
TH2F * h_mapDepth3Ampl047_HF
TH1F * h_runnbadchannelsC_depth1_HB
TH1F * h_sum0TSmeanAperLS2
TH2F * h_mapenophinorm_HE7
TH1F * h_bcnbadrate0_depth1_HE
TH1F * h_averNOSIGNALsumamplitude_HF
TH2F * h_mapDepth4Ped2_HO
double sumEstimator6[nsub][ndepth][neta][nphi]
TH2F * h_2DsumADCAmplLS6_LSselected
double lsdep_estimator1_HEdepth6_
TH1F * h_sumCutTSmaxAperLS6
TH2F * h_mapDepth7Error_HE
TH2F * h_2D0sumADCAmplLSdepth4HBu
TH2F * h_maprphinorm2_HE4
TH2F * h_recNoiseEnergy1_HE2
TH1F * h_sumCutTSmaxAperLS8
TH2F * h_mapDepth4TSmeanA_HO
TH2F * h_mapDepth3Ped0_HE
T getParameter(std::string const &) const
TH1F * h_averSIGNALsumamplitude_HE
TH2F * h_mapRatioCalib047_HF
double signal3[ndepth][neta][nphi]
TH2F * h_mapDepth2Ped1_HF
TH1F * h_bcnbadrate0_depth1_HB
const HcalGain * getGain(const HcalGenericDetId &fId) const
TH2F * h_2D0sumAmplitudeLS7
TH1F * h_bcnnbadchannels_depth2_HB
TH2F * h_mapDepth1Pedw2_HB
TH1F * h_sumCutADCAmplperLS1
TH1F * h_sumCutAmplperLS2
TH2F * h_mapDepth6TSmaxA225_HE
TH2F * h_mapDepth4pedestalw_HO
TH1F * h_Eta16Depth3RADDAM_HED2
const HcalQIECoder * getHcalCoder(const HcalGenericDetId &fId) const
TH1F * h_sumamplitude_depth1_HB1
TH2F * h_mapDepth2TSmeanA_HF
TH1F * h_pedestalaver9_HF
double averSIGNALoccupancy_HE
TH2F * h_amplitudechannel2_HF2
double lsdep_estimator1_HEdepth2_
TH2F * h_mapDepth6ADCAmpl_HE
TH2F * h_mapDepth5Amplitude_HE
TH2F * h_mapDepth3TSmaxA225_HF
TH2F * h_mapDepth4Error_HF
TH2F * h_mapDepth1ADCAmpl_HF
TH1F * h_sum0PedestalperLS7
TH1F * h_runnbadchannelsC_depth4_HO
double lsdep_estimator1_HFdepth1_
TH2F * h_mapDepth4pedestalw_HB
TH1F * h_sigLayer1RADDAM5_HE
TH2F * h_mapDepth4pedestalw_HE
int badchannels[nsub][ndepth][neta][nphi]
TH1F * h_sum0ADCAmplperLS1_P1
TH1F * h_eventoccupancy_HO
TH1F * h_orbitNumvsamplitude_HB
TH2F * h_amplitudechannel1_HE7
TH1F * h_shape0_good_channels_HF
TH2F * h_mapRatioCalib047_HO
TH2F * h_mapDepth3Ampl_HF
TH2F * h_mapDepth1ADCAmpl12SiPM_HE
TH1F * h_sumphiEta16Depth3RADDAM_HED2P
constexpr int samples() const
total number of samples in the digi
TH2F * h_2D0sumADCAmplLSdepth5HEu
TH2F * h_mapDepth2pedestal_HF
TH2F * h_mapDepth7TSmaxA_HE
TH2F * h_mapDepth6ADCAmpl225Copy_HE
TH2F * h_mapDepth3TSmaxA_HB
TH1F * h_RatioOccupancy_HEM
TH1F * h_sumCutAmplitudeperLS6
TH1F * h_bcnnbadchannels_depth1_HF
TH2F * h_recNoiseEnergy2_HB3
TH2F * h_mapDepth1Amplitude225_HB
TH2F * h_2DsumADCAmplLS5_LSselected
TH1F * h_averNOSIGNALsumamplitude_HB
TH2F * h_mapDepth2Pedw2_HE
TH1F * h_sumADCAmplperLS3
TH1F * h_sumCutTSmaxAperLS3
TH1F * h_sumPedestalperLS6
TH1F * h_sum0TSmaxAperLS5
TH1F * h_sumCutAmplperLS8
TH2F * h_recNoiseEnergy0_HF2
TH2F * h_mapDepth1pedestalw_HE
TH1F * h_bcnnbadchannels_depth1_HE
TH2F * h_2DsumADCAmplLSdepth6HEu
double averNOSIGNALoccupancy_HF
TH2F * h_mapTSmeanCalib_HE
TH2F * h_mapDepth1ADCAmpl225Copy_HE
TH1F * h_totalAmplitudeHF
TH2F * h_mapDepth7Ampl_HE
TH1F * h_sum0TSmeanAperLS1
TH1F * h_averNOSIGNALoccupancy_HE
edm::EDGetTokenT< HFDigiCollection > tok_hf_
TH2F * h_mapDepth2ADCAmpl225_HF
TH1F * h_runbadrateC_depth3_HE
TH1F * h_sum0ErrorBperLS8
TH2F * h_recSignalEnergy2_HB4
double lsdep_estimator5_HEdepth2_
double signal[ndepth][neta][nphi]
TH1F * h_sumADCAmplEtaPhiLs_lscounterM1orbitNum
TH1F * h_RatioOccupancy_HFP
TH1F * h_totalAmplitudeHB
TH2F * h_mapenophinorm_HE3
double lsdep_estimator3_HFdepth2_
TH1F * h_shape_bad_channels_HE
TH2F * h_mapDepth4Ampl047_HE
TH2F * h_mapDepth1Pedw0_HE
TH1F * h_numberofhitsHBtest
TH1F * h_eventamplitude_HE
TH2F * h_mapGetRMSOverNormalizedSignal_HB
TH2F * h_mapDepth1Ampl047_HF
double lsdep_estimator1_HEdepth1_
TH2F * h_2D0sumTSmeanALS3
TH2F * h_mapDepth6TSmeanA225_HE
TH1F * h_sumTSmeanAperLS5
TH1F * h_corrforxaMAIN_HE
TH1F * h_averSIGNALoccupancy_HF
TH1F * h_A_Depth3RADDAM_HE
TH1F * h_bcnvsamplitude_HF
TH2F * h_mapDepth3Pedw1_HE
TH2F * h_recSignalEnergy1_HF4
TH2F * h_recNoiseEnergy0_HF4
TH1F * h_sigLayer2RADDAM6_HED2
TH1F * h_sumCutTSmeanAperLS6
TH2F * h_mapDepth4Pedw2_HO
TH2F * h_mapDepth6TSmaxA_HE
double calibrTSmeanHEMin_
TH2F * h_mapDepth2ADCAmpl12_HF
double calibrTSmeanHBMin_
TH2F * h_2D0sumAmplitudeLS3
TH2F * h_maprphinorm2_HE3
TH2F * h_2DsumADCAmplLS7_LSselected
int mapRADDAM_HED20[ndepth][neta]
TH2F * h_mapDepth2pedestalw_HE
TH1F * h_RatioOccupancy_HBM
TH1F * h_sum0ADCAmplperLSdepth4HEu
long int gsmdepth1sipm[npfit][neta][nphi][ndepth]
int flagabortgaprejected_
TH2F * h_mapDepth2RADDAM0_HE
TH2F * h_mapDepth7Amplitude_HE
TH2F * h_mapDepth1pedestalw_HF
TH1F * h_sumamplitude_depth1_HE0
TH2F * h_mapDepth1pedestal_HF
TH2F * h_mapDepth5pedestal_HE
TH2F * h_mapRatioCalib_HF
void fillDigiErrorsHO(HODigiCollection::const_iterator &digiItr)
TH1F * h_A_Depth1RADDAM_HE
TH1F * h_averNOSIGNALoccupancy_HF
TH1F * h_sum0ErrorBperLS7
TH1F * h_sumamplitude_depth1_HF
TH1F * h_runbadrate0_depth1_HE
TH2F * h_mapDepth3Ampl_HE
const HcalElectronicsMap * getHcalMapping() const
TH2F * h_mapenophinorm2_HE4
TH1F * h_sum0ADCAmplperLS8_M1
TH2F * h_2DsumADCAmplLS2_LSselected
int lsdep_cut1_peak_HEdepth1_
TH2F * h_mapDepth3TSmeanA225_HB
TH2F * h_mapDepth4pedestal_HF
TH1F * h_sum0AmplitudeperLS6
double lsdep_estimator2_HBdepth2_
TH2F * h_mapDepth2pedestalw_HB
TH2F * h_2DsumADCAmplEtaPhiLs30
TH2F * h_mapDepth6Ampl047_HE
TH1F * h_sum0ADCAmplperLS6u
TH2F * h_mapDepth1Ped1_HF
TH2F * h_mapDepth4ADCAmpl225_HE
TH2F * h_recNoiseEnergy2_HB2
edm::EDGetTokenT< HODigiCollection > tok_ho_
TH2F * h_mapAiForLS_bad_HO
TH1F * h_bcnbadrate0_depth3_HE
TH1F * h_sumADCAmplEtaPhiLs_orbitNum
TH1F * h_AmplitudeHBtest6
TH2F * h_mapDepth4TSmaxA225_HE
TH2F * h_recNoiseEnergy2_HB4
TH2F * h_mapADCAmplsecondpeak0_HE
TH2F * h_mapDepth2ADCAmpl225Copy_HF
TH1F * h_gsmdifferencefit2_HE
TH1F * h_averSIGNALsumamplitude_HB
TH1F * h_sumCutAmplitudeperLS2
TH2F * h_mapDepth2Ampl047_HF
TH1F * h_sum0ADCAmplperLS4
TH2F * h_mapDepth3ADCAmplSiPM_HE
TH1F * h_ADCAmpl345_HBCapIdError
void fillDigiErrorsHFQIE10(QIE10DataFrame qie10df)
TH1F * h_bcnnbadchannels_depth2_HE
void endRun(const edm::Run &r, const edm::EventSetup &iSetup) override
double calibrTSmeanHBMax_
TH2F * h_mapDepth2ADCAmpl_HE
TH2F * h_mapDepth3AmplE34_HF
TH2F * h_mapDepth4ADCAmpl225Copy_HE
TH1F * h_sumADCAmplEtaPhiLs_bbb1
TH2F * h_mapDepth1pedestal_HE
TH1F * h_sumAmplitudeperLS2
TH2F * h_amplitudechannel1_HE3
TH1F * h_AamplitudewithPedSubtr_RADDAM_HE
TH1F * h_bcnnbadchannels_depth3_HE
TH2F * h_mapDepth1TSmeanA_HF
double lsdep_estimator3_HEdepth3_
TH1F * h_gsmdifferencefit6_HE
double lsdep_estimator4_HEdepth1_
TH2F * h_mapDepth4Ped3_HO
TH2F * h_mapDepth2Ped3_HE
TH1F * h_sumCutAmplitudeperLS5
TH1F * h_sum0ADCAmplperLS1_M2
TH1F * h_AamplitudewithPedSubtr_RADDAM_HEzoom1
double averNOSIGNALsumamplitude_HE
TH1F * h_sigLayer2RADDAM0_HE
TH1F * h_sumamplitudechannel_HO
TH1F * h_pedestalaver4_HB
TH2F * h_mapWidthCalib_HB
TH1F * h_sumAmplitudeperLS7
TH2F * h_recSignalEnergy1_HF1
TH1F * h_sum0ADCAmplperLS1_M1
TH2F * h_mapDepth1TSmaxA225_HB
TH1F * h_orbitNumvsamplitude0_HB
TH2F * h_mapDepth2ADCAmpl12_HE
int lsdep_cut3_max_HBdepth1_
double recNoiseEnergy1[nsub][ndepth][neta][nphi]
TH2F * h_recNoiseEnergy1_HE7
double sumEstimator4[nsub][ndepth][neta][nphi]
TH2F * h_mapDepth2Pedw0_HB
double averSIGNALsumamplitude_HB
TH2F * h_maprphinorm0_HE6
constexpr int depth() const
get the tower depth
TH2F * h_amplitudechannel2_HE4
TH1F * h_sumCutTSmaxAperLS7
double averSIGNALsumamplitude_HF
TH2F * h_mapTSmeanCalib_HO
TH2F * h_mapDepth2AmplE34_HF
void fillDigiAmplitudeHFQIE10(QIE10DataFrame qie10df)
TH1F * h_eventamplitude_HF
TH1F * h_sum0ADCAmplperLSdepth6HEu
TH2F * h_mapDepth3pedestal_HE
TH1F * h_numberofhitsHOtest
TH1F * h_sumamplitude_depth2_HF1
TH2F * h_2D0sumTSmeanALS1
TH1F * h_totalAmplitudeHFperEvent
TH1F * h_orbitNumvsamplitude_HO
TH2F * h_mapDepth3ADCAmpl_HE
TH1F * h_totalAmplitudeHBperEvent
TH2F * h_2DsumAmplitudeLS7
TH2F * h_mapGetRMSOverNormalizedSignal0_HB
TH2F * h_mapDepth1ADCAmpl225Copy_HF
TH2F * h_recSignalEnergy2_HB1
TH2F * h_mapenophinorm2_HE3
TH2F * h_recSignalEnergy2_HE1
double lsdep_estimator5_HEdepth3_
TH2F * h_recSignalEnergy0_HB4
double const adc2fC_QIE11_shunt1[256]
TH1F * h_runbadrate0_depth1_HB
TH1F * h_sum0AmplitudeperLS7
TH2F * h_mapDepth4ADCAmpl_HF
TH2F * h_mapDepth3linADCAmpl12_HE
TH2F * h_recSignalEnergy0_HE5
TH1F * h_shape_good_channels_HB
TH2F * h_map0AiForLS_good_HE
TH2F * h_mapDepth1Ped1_HB
TH2F * h_recSignalEnergy0_HE7
TH2F * h_2DsumPedestalLS6
TH1F * h_nevents_per_eachRealLS
double lsdep_estimator1_HEdepth7_
TH1F * h_shape0_Ahigh_HB1
const HcalPedestal * getPedestal(const HcalGenericDetId &fId) const
TH1F * h_amplitudeaveragedbydepthes_HE
TH1F * h_sumADCAmplperLS3_P2
TH2F * h_mapDepth4Amplitude_HB
TH2F * h_mapRatioCalib047_HB
TH2F * h_mapDepth1ADCAmpl_HE
TH1F * h_shape0_bad_channels_HE
TH1F * h_shape0_bad_channels_HO
TH2F * h_amplitudechannel1_HF2
TH2F * h_mapDepth2Pedw0_HE
TH2F * h_2D0sumADCAmplLSdepth7HEu
TH2F * h_mapDepth6Error_HE
TH1F * h_sumTSmeanAperLS6
TH1F * h_sumADCAmplperLSdepth4HEu
TH2F * h_2D0sumADCAmplLSdepth4HEu
TH1F * h_sumPedestalperLS5
TH2F * h_mapDepth1Amplitude225_HF
TH2F * h_mapWidthCalib_HE
TH1F * h_sumCutTSmaxAperLS5
TH2F * h_mapDepth3ADCAmpl12SiPM_HE
TH2F * h_mapTSmaxCalib047_HO
TH1F * h_sigLayer2RADDAM5_HE
double calibrTSmeanHFMax_
TH2F * h_amplitudechannel0_HB2
TH2F * h_mapDepth4TSmeanA225_HF
TH2F * h_mapDepth4TSmaxA_HE
TH2F * h_mapGetRMSOverNormalizedSignal0_HF
double averSIGNALsumamplitude_HE
TH2F * h_mapDepth4ADCAmpl225_HO
TH2F * h_2DsumADCAmplLSdepth7HEu
TH2F * h_mapDepth3Pedw0_HE
Readout chain identification for Hcal.
TH1F * h_sumADCAmplperLS6_M2
TH1F * h_bcnvsamplitude0_HE
TH1F * h_sum0TSmeanAperLS5
TH1F * h_bcnbadrate0_depth2_HE
TH2F * h_mapAiForLS_bad_HB
TH2F * h_2D0sumADCAmplLS7
TH2F * h_mapDepth2Pedw1_HF
TH2F * h_mapDepth1Ped0_HF
TH2F * h_mapDepth1RADDAM_HE
TH1F * h_NphiForEta16Depth3RADDAM_HED2
TH2F * h_mapDepth4TSmeanA225_HO
TH2F * h_recNoiseEnergy2_HF2
TH1F * h_sumCutTSmaxAperLS4
TH2F * h_mapDepth4TSmaxA225_HO
TH1F * h_Amplitude_notCapIdErrors_HE2
TH2F * h_mapDepth5ADCAmpl225Copy_HE
TH2F * h_mapDepth7ADCAmpl12_HE
TH1F * h_gsmdifferencefit1_HE
TH2F * h_recNoiseEnergy0_HE1
TH2F * h_mapenophinorm2_HE2
TH2F * h_mapDepth1ADCAmplSiPM_HE
const DetId lookup(HcalElectronicsId fId) const
lookup the logical detid associated with the given electronics id
TH1F * h_sum0PedestalperLS2
TH2F * h_mapDepth3AmplE34_HB
TH2F * h_2D0sumAmplitudeLS4
TH1F * h_sumADCAmplEtaPhiLs_lscounterM1
TH2F * h_mapDepth6AmplE34_HE
Power< A, B >::type pow(const A &a, const B &b)
TH1F * h_sum0TSmeanAperLS8
TH2F * h_mapDepth1TSmeanA_HE
TH2F * h_mapDepth3ADCAmpl225Copy_HB
TH2F * h_recNoiseEnergy1_HE4
TH2F * h_mapenophinorm2_HE6
TH1F * h_runbadrateC_depth1_HB
TH2F * h_amplitudechannel1_HE1
TH1F * h_shape_bad_channels_HB
TH1F * h_Amplitude_notCapIdErrors_HF1
double averSIGNALoccupancy_HF
TH2F * h_mapDepth3AmplE34_HE
TH2F * h_recNoiseEnergy2_HE2
static constexpr float b1
TH2F * h_mapDepth3TSmeanA_HE
const edm::ESGetToken< HcalDbService, HcalDbRecord > tokDB_
TH2F * h_mapDepth4pedestal_HO
TH2F * h_mapDepth4TSmeanA_HE
TH1F * h_sumCutAmplperLS3
TH2F * h_mapDepth1Amplitude_HB
double dR(double eta1, double phi1, double eta2, double phi2)
uint16_t *__restrict__ uint16_t const *__restrict__ adc
TH1F * h_sumPedestalperLS8
TH2F * h_mapDepth4Amplitude225_HE
edm::EDGetTokenT< HcalCalibDigiCollection > tok_calib_
TH2F * h_recSignalEnergy2_HE5
TH1F * h_orbitNumvsamplitude_HE
TH1F * h_sumADCAmplperLSdepth7HEu
double recNoiseEnergy2[nsub][ndepth][neta][nphi]
TH2F * h_mapDepth3RADDAM0_HE
TH2F * h_mapDepth2pedestal_HE
void fillDigiAmplitude(HBHEDigiCollection::const_iterator &digiItr)
TH2F * h_amplitudechannel0_HE5
TH1F * h_averSIGNALoccupancy_HE
TH1F * h_runnbadchannels_depth2_HF
TH2F * h_2DsumADCAmplLSdepth4HFu
TH1F * h_bcnvsamplitude0_HO
TH1F * h_sumCutADCAmplperLSdepth5HEu
TH1F * h_nbadchannels_depth1_HF
TH2F * h_recSignalEnergy0_HE3
TH2F * h_mapDepth1Ped0_HE
TH2F * h_map0AiForLS_good_HB
TH2F * h_mapRatioCalib047_HE
TH2F * h_amplitudechannel2_HE2
TH2F * h_mapDepth3TSmaxA_HE
constexpr DetId detid() const
Get the detector id.
TH1F * h_sumCutADCAmplperLS7
TH2F * h_2D0sumADCAmplLS4
TH1F * h_sumCutAmplitudeperLS4
bool studyCapIDErrorsHist_
bool studyPedestalCorrelations_