351 1.89, 5.07, 8.25, 11.43, 14.61, 17.78, 20.96, 24.14, 27.32, 30.50, 33.68,
352 36.86, 40.04, 43.22, 46.40, 49.58, 54.35, 60.71, 67.07, 73.43, 79.79, 86.15,
353 92.51, 98.87, 105.2, 111.6, 117.9, 124.3, 130.7, 137.0, 143.4, 149.7, 156.1,
354 162.5, 168.8, 175.2, 184.7, 197.4, 210.2, 222.9, 235.6, 248.3, 261.0, 273.7,
355 286.5, 299.2, 311.9, 324.6, 337.3, 350.1, 362.8, 375.5, 388.2, 400.9, 413.6,
356 426.4, 439.1, 458.2, 483.6, 509.0, 534.5, 559.9, 585.3, 610.8, 558.9, 584.2,
357 609.5, 634.7, 660.0, 685.3, 710.6, 735.9, 761.2, 786.5, 811.8, 837.1, 862.4,
358 887.7, 913.0, 938.3, 976.2, 1026.8, 1077.4, 1128.0, 1178.6, 1229.1, 1279.7, 1330.3,
359 1380.9, 1431.5, 1482.1, 1532.7, 1583.3, 1633.8, 1684.4, 1735.0, 1785.6, 1836.2, 1886.8,
360 1937.4, 2013.2, 2114.4, 2215.6, 2316.8, 2417.9, 2519.1, 2620.3, 2721.5, 2822.6, 2923.8,
361 3025.0, 3126.2, 3227.3, 3328.5, 3429.7, 3530.9, 3632.0, 3733.2, 3834.4, 3935.5, 4036.7,
362 4188.5, 4390.8, 4593.2, 4795.5, 4997.9, 5200.2, 5402.6, 5057.5, 5262.3, 5467.1, 5671.8,
363 5876.6, 6081.4, 6286.2, 6491.0, 6695.8, 6900.6, 7105.3, 7310.1, 7514.9, 7719.7, 7924.5,
364 8129.3, 8436.4, 8846.0, 9255.6, 9665.1, 10074.7, 10484.3, 10893.9, 11303.4, 11713.0, 12122.6,
365 12532.1, 12941.7, 13351.3, 13760.8, 14170.4, 14580.0, 14989.5, 15399.1, 15808.7, 16218.2, 16832.6,
366 17651.7, 18470.9, 19290.0, 20109.2, 20928.3, 21747.4, 22566.6, 23385.7, 24204.8, 25024.0, 25843.1,
367 26662.3, 27481.4, 28300.5, 29119.7, 29938.8, 30757.9, 31577.1, 32396.2, 33215.4, 34444.1, 36082.3,
368 37720.6, 39358.9, 40997.2, 42635.4, 44273.7, 40908.7, 42542.6, 44176.5, 45810.4, 47444.3, 49078.3,
369 50712.2, 52346.1, 53980.0, 55613.9, 57247.8, 58881.8, 60515.7, 62149.6, 63783.5, 65417.4, 67868.3,
370 71136.1, 74404.0, 77671.8, 80939.7, 84207.5, 87475.3, 90743.2, 94011.0, 97278.8, 100546.7, 103814.5,
371 107082.3, 110350.2, 113618.0, 116885.8, 120153.7, 123421.5, 126689.3, 129957.2, 134858.9, 141394.6, 147930.3,
372 154465.9, 161001.6, 167537.3, 174072.9, 180608.6, 187144.3, 193679.9, 200215.6, 206751.3, 213287.0, 219822.6,
373 226358.3, 232894.0, 239429.6, 245965.3, 252501.0, 259036.6, 265572.3, 275375.8, 288447.2, 301518.5, 314589.8,
374 327661.2, 340732.5, 353803.8};
377 9.56, 28.24, 46.91, 65.59, 84.27, 102.9, 121.6, 140.3, 159.0, 177.7,
378 196.3, 215.0, 233.7, 252.4, 271.0, 289.7, 317.7, 355.1, 392.4, 429.8,
379 467.1, 504.5, 541.9, 579.2, 616.6, 653.9, 691.3, 728.6, 766.0, 803.3,
380 840.7, 878.0, 915.4, 952.8, 990.1, 1027.5, 1083.5, 1158.2, 1232.9, 1307.6,
381 1382.3, 1457.0, 1531.7, 1606.4, 1681.2, 1755.9, 1830.6, 1905.3, 1980.0, 2054.7,
382 2129.4, 2204.1, 2278.8, 2353.5, 2428.2, 2502.9, 2577.7, 2689.7, 2839.1, 2988.6,
383 3138.0, 3287.4, 3436.8, 3586.2, 3263.0, 3411.3, 3559.6, 3707.9, 3856.2, 4004.5,
384 4152.9, 4301.2, 4449.5, 4597.8, 4746.1, 4894.4, 5042.7, 5191.0, 5339.4, 5487.7,
385 5710.1, 6006.8, 6303.4, 6600.0, 6896.6, 7193.3, 7489.9, 7786.5, 8083.1, 8379.8,
386 8676.4, 8973.0, 9269.6, 9566.3, 9862.9, 10159.5, 10456.2, 10752.8, 11049.4, 11346.0,
387 11791.0, 12384.2, 12977.5, 13570.7, 14164.0, 14757.2, 15350.5, 15943.7, 16537.0, 17130.2,
388 17723.5, 18316.7, 18910.0, 19503.2, 20096.5, 20689.7, 21283.0, 21876.2, 22469.5, 23062.8,
389 23656.0, 24545.9, 25732.4, 26918.9, 28105.4, 29291.9, 30478.4, 31664.9, 29399.4, 30590.1,
390 31780.9, 32971.7, 34162.4, 35353.2, 36544.0, 37734.7, 38925.5, 40116.3, 41307.0, 42497.8,
391 43688.5, 44879.3, 46070.1, 47260.8, 49047.0, 51428.5, 53810.1, 56191.6, 58573.1, 60954.6,
392 63336.2, 65717.7, 68099.2, 70480.8, 72862.3, 75243.8, 77625.4, 80006.9, 82388.4, 84769.9,
393 87151.5, 89533.0, 91914.5, 94296.1, 97868.4, 102631.4, 107394.5, 112157.5, 116920.6, 121683.7,
394 126446.7, 131209.8, 135972.8, 140735.9, 145499.0, 150262.0, 155025.1, 159788.2, 164551.2, 169314.3,
395 174077.3, 178840.4, 183603.5, 188366.5, 193129.6, 200274.2, 209800.3, 219326.4, 228852.5, 238378.7,
396 247904.8, 257430.9, 237822.7, 247326.7, 256830.7, 266334.8, 275838.8, 285342.9, 294846.9, 304351.0,
397 313855.0, 323359.1, 332863.1, 342367.1, 351871.2, 361375.2, 370879.3, 380383.3, 394639.4, 413647.5,
398 432655.6, 451663.6, 470671.7, 489679.8, 508687.9, 527696.0, 546704.1, 565712.2, 584720.3, 603728.3,
399 622736.4, 641744.5, 660752.6, 679760.7, 698768.8, 717776.9, 736785.0, 755793.0, 784305.2, 822321.3,
400 860337.5, 898353.7, 936369.9, 974386.0, 1012402.2, 1050418.4, 1088434.6, 1126450.7, 1164466.9, 1202483.1,
401 1240499.3, 1278515.4, 1316531.6, 1354547.8, 1392564.0, 1430580.1, 1468596.3, 1506612.5, 1544628.7, 1601652.9,
402 1677685.3, 1753717.6, 1829750.0, 1905782.3, 1981814.7, 2057847.0};
405 -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,
406 12.5, 13.5, 15., 17., 19., 21., 23., 25., 27., 29.5, 32.5, 35.5, 38.5,
407 42., 46., 50., 54.5, 59.5, 64.5, 59.5, 64.5, 69.5, 74.5, 79.5, 84.5, 89.5,
408 94.5, 99.5, 104.5, 109.5, 114.5, 119.5, 124.5, 129.5, 137., 147., 157., 167., 177.,
409 187., 197., 209.5, 224.5, 239.5, 254.5, 272., 292., 312., 334.5, 359.5, 384.5, 359.5,
410 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,
411 709.5, 747., 797., 847., 897., 947., 997., 1047., 1109.5, 1184.5, 1259.5, 1334.5, 1422.,
412 1522., 1622., 1734.5, 1859.5, 1984.5, 1859.5, 1984.5, 2109.5, 2234.5, 2359.5, 2484.5, 2609.5, 2734.5,
413 2859.5, 2984.5, 3109.5, 3234.5, 3359.5, 3484.5, 3609.5, 3797., 4047., 4297., 4547., 4797., 5047.,
414 5297., 5609.5, 5984.5, 6359.5, 6734.5, 7172., 7672., 8172., 8734.5, 9359.5, 9984.5};
438 double phi12(
double phi1,
double en1,
double phi2,
double en2);
439 double dR(
double eta1,
double phi1,
double eta2,
double phi2);
440 double dPhiWsign(
double phi1,
double phi2);
2453 std::cout <<
"============== endJob ===================================" << std::endl;
2455 std::cout <<
" --------------------------------------- " << std::endl;
2458 std::cout <<
" --------------------------------------------- " << std::endl;
2464 std::cout <<
"===== possible max number of events in ntuple(each 50th recorded) = " <<
nevent50 << std::endl;
2468 std::cout <<
"===== full number of events*HBHEdigis = " <<
nnnnnn << std::endl;
2473 std::cout <<
"===== Start writing user histograms =====" << std::endl;
2708 for (
int k1 = 0; k1 <
ndepth; k1++) {
2709 for (
int k2 = 0; k2 <
neta; k2++) {
2714 for (
int k3 = 0; k3 <
nphi; k3++) {
2773 int outabortgap = 1;
2783 std::cout <<
" --------------------------------------- " << std::endl;
2787 std::cout <<
" --------------------------------------------- " << std::endl;
2793 std::cout <<
" ---------***********************------------- " << std::endl;
2843 unsigned long int pcountall1 = 0;
2844 unsigned long int pcountall3 = 0;
2845 unsigned long int pcountall6 = 0;
2846 unsigned long int pcountall8 = 0;
2851 unsigned long int mcountall1 = 0;
2852 unsigned long int mcountall3 = 0;
2853 unsigned long int mcountall6 = 0;
2854 unsigned long int mcountall8 = 0;
2859 int pnnmin1 = 999999999;
2860 int pnnmin3 = 999999999;
2861 int pnnmin6 = 999999999;
2862 int pnnmin8 = 999999999;
2867 int mnnmin1 = 999999999;
2868 int mnnmin3 = 999999999;
2869 int mnnmin6 = 999999999;
2870 int mnnmin8 = 999999999;
2872 for (
int k1 = 0; k1 <
ndepth; k1++) {
2873 for (
int k3 = 0; k3 <
nphi; k3++) {
2874 for (
int k2 = 0; k2 <
neta; k2++) {
2983 int kkkk2 = (k2 - 1) / 4;
2988 int kkkk3 = (k3) / 4 + 1;
2992 ietaphi = ((kkkk2)-1) *
znphi + (kkkk3);
3060 if (bbbc / bbb1 > 25.) {
3075 if (bbbc / bbb1 > 25.) {
3160 if (bbbc / bbb1 > 15. && k3 % 2 == 0) {
3174 if (bbbc / bbb1 > 15. && k3 % 2 == 0) {
3314 if (bbbc / bbb1 > 20.) {
3329 if (bbbc / bbb1 > 20.) {
3421 if (bbbc / bbb1 > 80.) {
3436 if (bbbc / bbb1 > 80.) {
3991 if (
k0 == 0 && k1 == 0) {
3992 if (pcountmin1 > 0) {
3993 if (pcountmin1 < pnnmin1)
3994 pnnmin1 = pcountmin1;
3997 if (mcountmin1 > 0) {
3998 if (mcountmin1 < mnnmin1)
3999 mnnmin1 = mcountmin1;
4005 if (
k0 == 1 && k1 == 0) {
4006 if (pcountmin3 > 0) {
4007 if (pcountmin3 < pnnmin3)
4008 pnnmin3 = pcountmin3;
4011 if (mcountmin3 > 0) {
4012 if (mcountmin3 < mnnmin3)
4013 mnnmin3 = mcountmin3;
4019 if (
k0 == 2 && k1 == 3) {
4020 if (pcountmin8 > 0) {
4021 if (pcountmin8 < pnnmin8)
4022 pnnmin8 = pcountmin8;
4025 if (mcountmin8 > 0) {
4026 if (mcountmin8 < mnnmin8)
4027 mnnmin8 = mcountmin8;
4033 if (
k0 == 3 && k1 == 0) {
4034 if (pcountmin6 > 0) {
4035 if (pcountmin6 < pnnmin6)
4036 pnnmin6 = pcountmin6;
4039 if (mcountmin6 > 0) {
4040 if (mcountmin6 < mnnmin6)
4041 mnnmin6 = mcountmin6;
4056 float patiooccupancy1 = 0.;
4057 if (pcountall1 != 0)
4058 patiooccupancy1 = (
float)pnnmin1 * mnnbins1 / pcountall1;
4060 float matiooccupancy1 = 0.;
4061 if (mcountall1 != 0)
4062 matiooccupancy1 = (
float)mnnmin1 * mnnbins1 / mcountall1;
4065 float patiooccupancy3 = 0.;
4066 if (pcountall3 != 0)
4067 patiooccupancy3 = (
float)pnnmin3 * mnnbins3 / pcountall3;
4069 float matiooccupancy3 = 0.;
4070 if (mcountall3 != 0)
4071 matiooccupancy3 = (
float)mnnmin3 * mnnbins3 / mcountall3;
4074 float patiooccupancy6 = 0.;
4075 if (pcountall6 != 0)
4076 patiooccupancy6 = (
float)pnnmin6 * mnnbins6 / pcountall6;
4078 float matiooccupancy6 = 0.;
4079 if (mcountall6 != 0)
4080 matiooccupancy6 = (
float)mnnmin6 * mnnbins6 / mcountall6;
4083 float patiooccupancy8 = 0.;
4084 if (pcountall8 != 0)
4085 patiooccupancy8 = (
float)pnnmin8 * mnnbins8 / pcountall8;
4087 float matiooccupancy8 = 0.;
4088 if (mcountall8 != 0)
4089 matiooccupancy8 = (
float)mnnmin8 * mnnbins8 / mcountall8;
4093 for (
int k1 = 0; k1 <
ndepth; k1++) {
4094 for (
int k2 = 0; k2 <
neta; k2++) {
4095 for (
int k3 = 0; k3 <
nphi; k3++) {
4195 for (
int k1 = 0; k1 <
ndepth; k1++) {
4196 for (
int k2 = 0; k2 <
neta; k2++) {
4197 for (
int k3 = 0; k3 <
nphi; k3++) {
4219 for (
int k1 = 0; k1 <
ndepth; k1++) {
4220 for (
int k2 = 0; k2 <
neta; k2++) {
4221 for (
int k3 = 0; k3 <
nphi; k3++) {
4253 bool gotHFDigis =
true;
4257 if (!(
hf.isValid())) {
4261 std::cout <<
" ****************************** =========================== No HFDigiCollection found " 4266 eta = digi->id().ieta();
4267 phi = digi->id().iphi();
4268 depth = digi->id().depth();
4285 for (
int i = 0;
i <
nTS;
i++) {
4304 bool gotQIE10Digis =
true;
4306 gotQIE10Digis =
false;
4309 gotQIE10Digis =
false;
4311 if (!gotQIE10Digis) {
4312 std::cout <<
" No QIE10DigiCollection collection is found " << std::endl;
4315 double totalAmplitudeHF = 0.;
4316 for (
unsigned int j = 0;
j < qie10dc.
size();
j++) {
4340 double amplitudefullTSs = 0.;
4341 double nnnnnnTS = 0.;
4342 for (
int i = 0;
i <
nTS; ++
i) {
4345 int adc = qie10df[
i].adc();
4377 bool gotHBHEDigis =
true;
4379 gotHBHEDigis =
false;
4380 if (!(
hbhe.isValid()))
4381 gotHBHEDigis =
false;
4382 if (!gotHBHEDigis) {
4383 std::cout <<
" No HBHEDigiCollection collection is found " << std::endl;
4386 double totalAmplitudeHB = 0.;
4387 double totalAmplitudeHE = 0.;
4388 double nnnnnnTSHB = 0.;
4389 double nnnnnnTSHE = 0.;
4392 eta = digi->id().ieta();
4393 phi = digi->id().iphi();
4394 depth = digi->id().depth();
4434 double amplitudefullTSs = 0.;
4437 for (
int i = 0;
i <
nTS;
i++) {
4449 double amplitudefullTSs = 0.;
4452 for (
int i = 0;
i <
nTS;
i++) {
4467 if (totalAmplitudeHB != 0.) {
4472 if (totalAmplitudeHE != 0.) {
4488 bool gotQIE11Digis =
true;
4490 gotQIE11Digis =
false;
4492 gotQIE11Digis =
false;
4493 if (!gotQIE11Digis) {
4494 std::cout <<
" No QIE11DigiCollection collection is found " << std::endl;
4497 double totalAmplitudeHBQIE11 = 0.;
4498 double totalAmplitudeHEQIE11 = 0.;
4499 double nnnnnnTSHBQIE11 = 0.;
4500 double nnnnnnTSHEQIE11 = 0.;
4501 for (
unsigned int j = 0;
j < qie11dc.
size();
j++) {
4510 int sub = hcaldetid.
subdet();
4523 if (sub == 1 &&
depth == 1 && qwert1 == 0) {
4527 if (sub == 1 &&
depth == 2 && qwert2 == 0) {
4531 if (sub == 2 &&
depth == 1 && qwert3 == 0) {
4535 if (sub == 2 &&
depth == 2 && qwert4 == 0) {
4539 if (sub == 2 &&
depth == 3 && qwert5 == 0) {
4551 double amplitudefullTSs1 = 0.;
4552 double amplitudefullTSs6 = 0.;
4554 for (
int i = 0;
i <
nTS; ++
i) {
4555 int adc = qie11df[
i].adc();
4558 amplitudefullTSs1 += charge1;
4559 amplitudefullTSs6 += charge6;
4565 totalAmplitudeHBQIE11 +=
charge;
4572 double amplitudefullTSs1 = 0.;
4573 double amplitudefullTSs6 = 0.;
4575 for (
int i = 0;
i <
nTS;
i++) {
4576 int adc = qie11df[
i].adc();
4579 amplitudefullTSs1 += charge1;
4580 amplitudefullTSs6 += charge6;
4586 totalAmplitudeHEQIE11 +=
charge;
4595 if (totalAmplitudeHBQIE11 != 0.) {
4600 if (totalAmplitudeHEQIE11 != 0.) {
4611 bool gotHODigis =
true;
4614 if (!(
ho.isValid()))
4618 std::cout <<
" No HO collection is found " << std::endl;
4621 double totalAmplitudeHO = 0.;
4623 eta = digi->id().ieta();
4624 phi = digi->id().iphi();
4625 depth = digi->id().depth();
4646 double nnnnnnTS = 0.;
4647 double amplitudefullTSs = 0.;
4649 for (
int i = 0;
i <
nTS;
i++) {
4678 bool gotHBHERecHitsNoise =
true;
4680 gotHBHERecHitsNoise =
4683 gotHBHERecHitsNoise =
false;
4684 if (!gotHBHERecHitsNoise) {
4686 std::cout <<
" No RecHits HBHENoise collection is found " << std::endl;
4691 float icalconst = 1.;
4694 HBHERecHit aHit(hbheItr->id(), hbheItr->eraw() * icalconst, hbheItr->time());
4696 double energyhit = aHit.
energy();
4697 DetId id = (*hbheItr).detid();
4699 int sub = ((hid).
rawId() >> 25) & 0x7;
4708 if (hid.
depth() == 1 && sub == 1 && hid.
iphi() == 25) {
4710 std::cout <<
" Noise,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4711 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4713 if (hid.
depth() == 1 && sub == 2 && hid.
iphi() == 25) {
4715 std::cout <<
" Noise,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4716 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4722 int mdepth = hid.
depth();
4733 bool gotHFRecHitsNoise =
true;
4735 gotHFRecHitsNoise =
false;
4737 gotHFRecHitsNoise =
false;
4738 if (!gotHFRecHitsNoise) {
4740 std::cout <<
" No RecHits HFNoise collection is found " << std::endl;
4744 float icalconst = 1.;
4747 HFRecHit aHit(hfItr->id(), hfItr->energy() * icalconst, hfItr->time());
4748 double energyhit = aHit.
energy();
4749 DetId id = (*hfItr).detid();
4751 int sub = ((hid).
rawId() >> 25) & 0x7;
4757 if (hid.
iphi() == 25) {
4759 std::cout <<
"HF Noise,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4760 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4766 int mdepth = hid.
depth();
4780 bool gotHBHERecHitsSignal =
true;
4782 gotHBHERecHitsSignal =
4785 gotHBHERecHitsSignal =
false;
4786 if (!gotHBHERecHitsSignal) {
4788 std::cout <<
" No RecHits HBHESignal collection is found " << std::endl;
4793 float icalconst = 1.;
4796 HBHERecHit aHit(hbheItr->id(), hbheItr->eraw() * icalconst, hbheItr->time());
4798 double energyhit = aHit.
energy();
4799 DetId id = (*hbheItr).detid();
4801 int sub = ((hid).
rawId() >> 25) & 0x7;
4809 if (hid.
depth() == 1 && sub == 1 && hid.
iphi() == 25) {
4811 std::cout <<
"HBHE Signal,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4812 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4814 if (hid.
depth() == 1 && sub == 2 && hid.
iphi() == 25) {
4816 std::cout <<
"HBHE Signal,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4817 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4823 int mdepth = hid.
depth();
4834 bool gotHFRecHitsSignal =
true;
4836 gotHFRecHitsSignal =
false;
4838 gotHFRecHitsSignal =
false;
4839 if (!gotHFRecHitsSignal) {
4841 std::cout <<
" No RecHits HFSignal collection is found " << std::endl;
4845 float icalconst = 1.;
4848 HFRecHit aHit(hfItr->id(), hfItr->energy() * icalconst, hfItr->time());
4849 double energyhit = aHit.
energy();
4850 DetId id = (*hfItr).detid();
4852 int sub = ((hid).
rawId() >> 25) & 0x7;
4858 if (hid.
iphi() == 25) {
4860 std::cout <<
"HF Signal,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4861 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4867 int mdepth = hid.
depth();
4897 for (
int k1 = 0; k1 <
ndepth; k1++) {
4899 for (
int k2 = 0; k2 <
neta; k2++) {
4901 int k2plot = k2 - 41;
4904 double sumoverphi = 0;
4905 int nsumoverphi = 0;
4906 for (
int k3 = 0; k3 <
nphi; k3++) {
4911 std::cout <<
"==== nsumoverphi = " << nsumoverphi <<
" sumoverphi = " << sumoverphi
4912 <<
" k1 = " << k1 <<
" k2 = " << k2 <<
" kkk = " << kkk <<
" k3 = " << k3 << std::endl;
4916 for (
int k3 = 0; k3 <
nphi; k3++) {
4917 if (nsumoverphi != 0) {
4920 std::cout <<
"nsumoverphi= " << nsumoverphi <<
" sumoverphi= " << sumoverphi <<
" k1= " << k1
4921 <<
" k2= " << k2 <<
" kkk= " << kkk <<
" k3= " << k3
4922 <<
" maprphinorm= " <<
maprphinorm[
k0][k1][k2][k3] << std::endl;
4997 for (
int k1 = 0; k1 <
ndepth; k1++) {
4999 for (
int k2 = 0; k2 <
neta; k2++) {
5001 int k2plot = k2 - 41;
5003 for (
int k3 = 0; k3 <
nphi; k3++) {
5031 for (
int k1 = 0; k1 <
ndepth; k1++) {
5033 for (
int k2 = 0; k2 <
neta; k2++) {
5035 int k2plot = k2 - 41;
5037 for (
int k3 = 0; k3 <
nphi; k3++) {
5080 for (
int k1 = 0; k1 <
ndepth; k1++) {
5082 for (
int k2 = 0; k2 <
neta; k2++) {
5084 int k2plot = k2 - 41;
5086 for (
int k3 = 0; k3 <
nphi; k3++) {
5129 for (
int k1 = 0; k1 <
ndepth; k1++) {
5131 for (
int k2 = 0; k2 <
neta; k2++) {
5133 int k2plot = k2 - 41;
5135 for (
int k3 = 0; k3 <
nphi; k3++) {
5175 for (
int k1 = 0; k1 <
ndepth; k1++) {
5177 for (
int k2 = 0; k2 <
neta; k2++) {
5179 int k2plot = k2 - 41;
5181 for (
int k3 = 0; k3 <
nphi; k3++) {
5245 for (
int k1 = 0; k1 <
ndepth; k1++) {
5247 for (
int k2 = 0; k2 <
neta; k2++) {
5249 int k2plot = k2 - 41;
5251 for (
int k3 = 0; k3 <
nphi; k3++) {
5299 for (
int k1 = 0; k1 <
ndepth; k1++) {
5300 for (
int k2 = 0; k2 <
neta; k2++) {
5301 for (
int k3 = 0; k3 <
nphi; k3++) {
5313 for (
int k1 = 0; k1 <
ndepth; k1++) {
5314 for (
int k2 = 0; k2 <
neta; k2++) {
5317 if (k1 == 2 && k2 == 25) {
5321 }
else if (k1 == 2 && k2 == 56) {
5336 for (
int k1 = 0; k1 <
ndepth; k1++) {
5337 for (
int k2 = 0; k2 <
neta; k2++) {
5338 if (k1 == 2 && (k2 == 25 || k2 == 56)) {
5341 int k2plot = k2 - 41;
5351 int LLLLLL111111 = 0;
5352 if ((k1 == 0 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 30))
5355 int LLLLLL222222 = 0;
5356 if ((k1 == 1 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 27) ||
5357 (k1 == 2 && fabs(kkk - 41) > 26 && fabs(kkk - 41) < 29))
5360 if (LLLLLL111111 == 1) {
5364 if (LLLLLL222222 == 1) {
5377 for (
int k1 = 0; k1 <
ndepth; k1++) {
5378 for (
int k2 = 0; k2 <
neta; k2++) {
5379 if (k1 == 2 && (k2 == 25 || k2 == 56)) {
5381 int k2plot = k2 - 41;
5388 for (
int k3 = 0; k3 <
nphi; k3++) {
5391 int LLLLLL111111 = 0;
5392 if ((k1 == 0 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 30))
5394 int LLLLLL222222 = 0;
5395 if ((k1 == 1 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 27) ||
5396 (k1 == 2 && fabs(kkk - 41) > 26 && fabs(kkk - 41) < 29))
5398 if (LLLLLL111111 == 1) {
5402 if (LLLLLL222222 == 1) {
5413 for (
int k1 = 0; k1 <
ndepth; k1++) {
5414 for (
int k2 = 0; k2 <
neta; k2++) {
5417 for (
int k3 = 0; k3 <
nphi; k3++) {
5432 int sumofchannels = 0;
5433 double sumamplitudesubdet = 0.;
5434 int sumofchannels0 = 0;
5435 double sumamplitudesubdet0 = 0.;
5436 for (
int k2 = 0; k2 <
neta; k2++) {
5437 for (
int k3 = 0; k3 <
nphi; k3++) {
5442 if (sumamplitudechannel_HB > 80.) {
5443 sumamplitudesubdet += sumamplitudechannel_HB;
5446 if (sumamplitudechannel_HB > 0.) {
5447 sumamplitudesubdet0 += sumamplitudechannel_HB;
5455 double sumamplitudechannel_HE =
5458 if (sumamplitudechannel_HE > 200.) {
5459 sumamplitudesubdet += sumamplitudechannel_HE;
5462 if (sumamplitudechannel_HE > 0.) {
5463 sumamplitudesubdet0 += sumamplitudechannel_HE;
5473 if (sumamplitudechannel_HO > 1200.) {
5474 sumamplitudesubdet += sumamplitudechannel_HO;
5477 if (sumamplitudechannel_HO > 0.) {
5478 sumamplitudesubdet0 += sumamplitudechannel_HO;
5487 if (sumamplitudechannel_HF > 600.) {
5488 sumamplitudesubdet += sumamplitudechannel_HF;
5491 if (sumamplitudechannel_HF > 0.) {
5492 sumamplitudesubdet0 += sumamplitudechannel_HF;
5505 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM1)
5506 maxxSUM1 = sumamplitudesubdet + sumamplitudesubdet0;
5507 if ((sumofchannels + sumofchannels0) >
maxxOCCUP1)
5513 if ((sumamplitudesubdet + sumamplitudesubdet0) > 60000) {
5514 for (
int k2 = 0; k2 <
neta; k2++) {
5515 for (
int k3 = 0; k3 <
nphi; k3++) {
5534 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM2)
5535 maxxSUM2 = sumamplitudesubdet + sumamplitudesubdet0;
5536 if ((sumofchannels + sumofchannels0) >
maxxOCCUP2)
5546 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM3)
5547 maxxSUM3 = sumamplitudesubdet + sumamplitudesubdet0;
5548 if ((sumofchannels + sumofchannels0) >
maxxOCCUP3)
5558 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM4)
5559 maxxSUM4 = sumamplitudesubdet + sumamplitudesubdet0;
5560 if ((sumofchannels + sumofchannels0) >
maxxOCCUP4)
5581 for (
int k1 = 0; k1 <
ndepth; k1++) {
5583 int nbadchannels = 0;
5584 for (
int k2 = 0; k2 <
neta; k2++) {
5585 for (
int k3 = 0; k3 <
nphi; k3++) {
5598 if (nbadchannels != 0)
5611 if (nbadchannels != 0)
5622 if (nbadchannels > 0.)
5624 if (nbadchannels != 0)
5626 if (nbadchannels > 0.)
5633 if (nbadchannels > 0.)
5635 if (nbadchannels != 0)
5637 if (nbadchannels > 0.)
5649 if (nbadchannels != 0)
5662 if (nbadchannels != 0)
5675 if (nbadchannels != 0)
5686 if (nbadchannels > 0.)
5688 if (nbadchannels != 0)
5690 if (nbadchannels > 0.)
5697 if (nbadchannels > 0.)
5699 if (nbadchannels != 0)
5701 if (nbadchannels > 0.)
5708 if (nbadchannels > 0.)
5710 if (nbadchannels != 0)
5712 if (nbadchannels > 0.)
5719 if (nbadchannels > 0.)
5721 if (nbadchannels != 0)
5723 if (nbadchannels > 0.)
5735 if (nbadchannels != 0)
5751 if (nbadchannels != 0)
5764 if (nbadchannels != 0)
5775 if (nbadchannels > 0.)
5777 if (nbadchannels != 0)
5779 if (nbadchannels > 0.)
5786 if (nbadchannels > 0.)
5788 if (nbadchannels != 0)
5790 if (nbadchannels > 0.)
5806 bool gotCALIBDigis =
true;
5808 gotCALIBDigis =
false;
5810 if (!(
calib.isValid())) {
5811 gotCALIBDigis =
false;
5813 if (!gotCALIBDigis) {
5816 int cal_det = digi->id().hcalSubdet();
5817 int cal_phi = digi->id().iphi();
5818 int cal_eta = digi->id().ieta();
5819 int cal_cbox = digi->id().cboxChannel();
5823 if (cal_det > 0 && cal_det < 5 && cal_cbox == 0) {
5824 int iphi = cal_phi - 1;
5829 double max_signal = -100.;
5830 int ts_with_max_signal = -100;
5835 for (
int i = 0;
i <
nTS;
i++) {
5836 double ampldefault =
adc2fC[digi->sample(
i).adc() & 0xff];
5837 if (max_signal < ampldefault) {
5838 max_signal = ampldefault;
5839 ts_with_max_signal =
i;
5844 timew += (
i + 1) * ampldefault;
5849 double aveamplitude = -100.;
5852 double aveamplitude1 = aveamplitude - 1;
5857 double ampldefault =
adc2fC[digi->sample(
ii).adc() & 0xff];
5858 double aaaaaa = (
ii + 1) - aveamplitude;
5859 double aaaaaa2 = aaaaaa * aaaaaa;
5860 rmsamp += (aaaaaa2 * ampldefault);
5862 double rmsamplitude = -100.;
5870 if (ts_with_max_signal > -1 && ts_with_max_signal <
nTS)
5872 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 <
nTS)
5873 calib2[cal_det - 1][
ieta + 41][
iphi] +=
adc2fC[digi->sample(ts_with_max_signal + 1).adc() & 0xff];
5874 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 <
nTS)
5875 calib2[cal_det - 1][
ieta + 41][
iphi] +=
adc2fC[digi->sample(ts_with_max_signal + 2).adc() & 0xff];
5876 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 <
nTS)
5877 calib2[cal_det - 1][
ieta + 41][
iphi] +=
adc2fC[digi->sample(ts_with_max_signal - 1).adc() & 0xff];
5878 if (ts_with_max_signal - 2 > -1 && ts_with_max_signal - 2 <
nTS)
5879 calib2[cal_det - 1][
ieta + 41][
iphi] +=
adc2fC[digi->sample(ts_with_max_signal - 2).adc() & 0xff];
5881 bool anycapid =
true;
5884 int error1 = 0, error2 = 0, error3 = 0;
5885 int lastcapid = 0, capid = 0;
5886 for (
int ii = 0;
ii < (*digi).size();
ii++) {
5887 capid = (*digi)[
ii].capid();
5888 bool er = (*digi)[
ii].er();
5889 bool dv = (*digi)[
ii].dv();
5890 if (
ii != 0 && ((lastcapid + 1) % 4) != capid)
5904 if (error1 != 0 || error2 != 0 || error3 != 0)
5919 for (
int k1 = 0; k1 <
nsub; k1++) {
5920 for (
int k2 = 0; k2 <
neta; k2++) {
5921 for (
int k3 = 0; k3 <
nphi; k3++) {
5922 int k2plot = k2 - 41;
5926 if (
signal[k1][k2][k3] > 0.) {
5952 if (k20 > 0 || k20 == 0)
5957 int kk2 = 0, kk3 = 0;
5958 if (k1 == 0 || k1 == 1) {
5963 if (k30 == 71 || k30 ==
nphi || k30 == 1 || k30 == 2)
5966 kk3 = ((k30 - 3) / 4) * 4 + 3;
5967 }
else if (k1 == 2) {
5968 if (
abs(k20) <= 4) {
5970 if (k30 == 71 || k30 ==
nphi || k30 == 1 || k30 == 2 || k30 == 3 || k30 == 4)
5973 kk3 = ((k30 - 5) / 6) * 6 + 5;
5975 if (
abs(k20) > 4 &&
abs(k20) <= 10)
5977 if (
abs(k20) > 10 &&
abs(k20) <= 15)
5981 if (k30 == 71 || k30 ==
nphi || (k30 >= 1 && k30 <= 10))
5984 kk3 = ((k30 - 11) / 12) * 12 + 11;
5986 }
else if (k1 == 3) {
5991 if (k30 >= 1 && k30 <= 18)
5993 if (k30 >= 19 && k30 <= 36)
5995 if (k30 >= 37 && k30 <= 54)
5997 if (k30 >= 55 && k30 <=
nphi)
6001 int kkk2 = kk2 + 41;
6008 double GetRMSOverNormalizedSignal = -1.;
6009 if (
signal[k1][k2][k3] > 0. &&
calib0[k1][kkk2][kkk3] > 0.) {
6010 GetRMSOverNormalizedSignal =
signal[k1][k2][k3] /
calib0[k1][kkk2][kkk3];
6030 if (
signal[k1][k2][k3] > 0.) {
6033 if (
calib0[k1][kkk2][kkk3] > 0.)
6037 if (
calib0[k1][kkk2][kkk3] > 0.)
6038 ratio =
calib2[k1][kkk2][kkk3] /
calib0[k1][kkk2][kkk3];
6040 float calibtsmax =
calibt[k1][kkk2][kkk3];
6042 float calibtsmean =
caliba[k1][kkk2][kkk3];
6044 float calibwidth =
calibw[k1][kkk2][kkk3];
6046 float calibcap = -100.;
6063 if (calibtsmax > -0.5) {
6070 if (calibtsmean > -0.5) {
6077 if (calibwidth > -0.5) {
6087 h_map_HB->Fill(
double(k2plot),
double(k3), 1.);
6103 if (calibtsmax > -0.5) {
6110 if (calibtsmean > -0.5) {
6117 if (calibwidth > -0.5) {
6127 h_map_HE->Fill(
double(k2plot),
double(k3), 1.);
6143 if (calibtsmax > -0.5) {
6150 if (calibtsmean > -0.5) {
6157 if (calibwidth > -0.5) {
6167 h_map_HO->Fill(
double(k2plot),
double(k3), 1.);
6183 if (calibtsmax > -0.5) {
6190 if (calibtsmean > -0.5) {
6197 if (calibwidth > -0.5) {
6207 h_map_HF->Fill(
double(k2plot),
double(k3), 1.);
6239 std::cout <<
"======================== beignJob START +++++++++++++++++++++++++++" << std::endl;
6489 fs_->
make<TH1F>(
"h_ndepthesperamplitudebins_HE",
" ", 10, 0., 10.);
6499 "h_mapADCAmplsecondpeak0_HE",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
6856 fs_->
make<TH1F>(
"h_AamplitudewithPedSubtr_RADDAM_HEzoom0",
" ", 100, min80, 4000.);
6858 fs_->
make<TH1F>(
"h_AamplitudewithPedSubtr_RADDAM_HEzoom1",
" ", 100, min80, 1000.);
6872 fs_->
make<TH1F>(
"h_sumphiEta16Depth3RADDAM_HED2ALL",
" ", 100, min90, 70. * max90);
7020 float zaR = 10000000.;
7021 float zaR2 = 50000000.;
7055 fs_->
make<TH1F>(
"h_sumADCAmplEtaPhiLs_lscounterM1orbitNum",
" ", 600, 1., 601.);
7063 float bac2 = bac + 1.;
8093 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HB",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
8095 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HB",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
8097 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HE",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
8099 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HE",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
8101 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HF",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
8103 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HF",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
8105 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HO",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
8107 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HO",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
8156 float spls0 = 10000;
8166 float spls1 = 100000;
8202 std::cout <<
"======================== booking DONE +++++++++++++++++++++++++++" << std::endl;
8211 std::cout <<
"======================== beignJob finish +++++++++++++++++++++++++++" << std::endl;
8224 int mdepth = cell.
depth();
8225 int iphi = cell.iphi() - 1;
8226 int ieta = cell.ieta();
8229 int sub = cell.subdet();
8231 int errorGeneral = 0;
8240 bool anycapid =
true;
8245 int lastcapid = 0, capid = 0;
8246 int ERRORfiber = -10;
8247 int ERRORfiberChan = -10;
8248 int ERRORfiberAndChan = -10;
8249 int repetedcapid = 0;
8251 TSsize = digiItr->size();
8254 for (
int ii = 0;
ii < TSsize;
ii++) {
8255 capid = (*digiItr)[
ii].capid();
8256 bool er = (*digiItr)[
ii].er();
8257 bool dv = (*digiItr)[
ii].dv();
8258 int fiber = (*digiItr)[
ii].fiber();
8259 int fiberChan = (*digiItr)[
ii].fiberChan();
8260 int fiberAndChan = (*digiItr)[
ii].fiberAndChan();
8261 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8264 ERRORfiberChan = fiberChan;
8265 ERRORfiberAndChan = fiberAndChan;
8266 if (capid != lastcapid) {
8279 ERRORfiberChan = fiberChan;
8280 ERRORfiberAndChan = fiberAndChan;
8285 ERRORfiberChan = fiberChan;
8286 ERRORfiberAndChan = fiberAndChan;
8292 if (firstcapid == 0 && !anycapid)
8294 if (firstcapid == 1 && !anycapid)
8296 if (firstcapid == 2 && !anycapid)
8298 if (firstcapid == 3 && !anycapid)
8307 if (!anycapid && anyer)
8309 if (!anycapid && !anydv)
8311 if (!anycapid && anyer && !anydv)
8313 if (anyer && !anydv)
8319 for (
int ii = 0;
ii < TSsize;
ii++) {
8320 double ampldefault =
adc2fC[digiItr->sample(
ii).adc()];
8321 ampl += ampldefault;
8338 if (error1 != 0 || error2 != 0 || error3 != 0) {
8365 if (error1 != 0 || error2 != 0 || error3 != 0) {
8396 int mdepth = hcaldetid.
depth();
8397 int sub = hcaldetid.
subdet();
8400 for (
int ii = 0;
ii < 10;
ii++) {
8404 bool sticking =
false;
8406 bool anycapid =
true;
8409 int error2 = 0, error3 = 0;
8412 int lastcapid = 0, capid = 0;
8413 int repetedcapid = 0;
8417 if (mdepth == 0 || sub > 4)
8429 capid = qie11df[
ii].capid();
8434 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8436 if (capid != lastcapid) {
8457 int adc = qie11df[
ii].adc();
8468 ampl += ampldefault;
8470 if (sub == 1 || sub == 2) {
8471 if (
verbosity == 111111 && (ampts[0] == ampts[1]) && (ampts[1] == ampts[2]) && (ampts[2] == ampts[3]) &&
8472 (ampts[3] == ampts[4]) && (ampts[4] == ampts[5]) && (ampts[5] == ampts[6]) && (ampts[6] == ampts[7])) {
8473 std::cout <<
" HBHE DIGI A,fC = " << ampl <<
" sub = " << sub <<
" mdepth = " << mdepth <<
" nTS = " <<
nTS 8475 std::cout <<
" Ai,adc counts: " << std::endl;
8476 std::cout <<
" A1 = " << ampts[0] <<
" A2 = " << ampts[1] <<
" A3 = " << ampts[2] <<
" A4 = " << ampts[3]
8478 std::cout <<
" A5 = " << ampts[4] <<
" A6 = " << ampts[5] <<
" A7 = " << ampts[6] <<
" A8 = " << ampts[7]
8481 std::cout <<
" icapid1 = " << icapid[0] <<
" icapid2 = " << icapid[1] <<
" icapid3 = " << icapid[2]
8482 <<
" icapid4 = " << icapid[3] << std::endl;
8483 std::cout <<
" icapid5 = " << icapid[4] <<
" icapid6 = " << icapid[5] <<
" icapid7 = " << icapid[6]
8484 <<
" icapid8 = " << icapid[7] << std::endl;
8486 if (
verbosity == 222222 && (ampts[3] > 95 && (ampts[3] == ampts[4]) && (ampts[4] == ampts[5]))) {
8487 std::cout <<
" HBHE DIGI A,fC = " << ampl <<
" sub = " << sub <<
" mdepth = " << mdepth <<
" nTS = " <<
nTS 8489 std::cout <<
" Ai,adc counts: " << std::endl;
8490 std::cout <<
" A1 = " << ampts[0] <<
" A2 = " << ampts[1] <<
" A3 = " << ampts[2] <<
" A4 = " << ampts[3]
8492 std::cout <<
" A5 = " << ampts[4] <<
" A6 = " << ampts[5] <<
" A7 = " << ampts[6] <<
" A8 = " << ampts[7]
8495 std::cout <<
" icapid1 = " << icapid[0] <<
" icapid2 = " << icapid[1] <<
" icapid3 = " << icapid[2]
8496 <<
" icapid4 = " << icapid[3] << std::endl;
8497 std::cout <<
" icapid5 = " << icapid[4] <<
" icapid6 = " << icapid[5] <<
" icapid7 = " << icapid[6]
8498 <<
" icapid8 = " << icapid[7] << std::endl;
8500 if (
verbosity == 333333 && (ampts[0] > 14 && ampts[0] == ampts[1]) && (ampts[1] == ampts[2]) &&
8501 (ampts[2] == ampts[3]) && (ampts[3] == ampts[4]) && (ampts[4] == ampts[5]) && (ampts[5] == ampts[6]) &&
8502 (ampts[6] == ampts[7])) {
8503 std::cout <<
" HBHE DIGI A,fC = " << ampl <<
" sub = " << sub <<
" mdepth = " << mdepth <<
" nTS = " <<
nTS 8505 std::cout <<
" Ai,adc counts: " << std::endl;
8506 std::cout <<
" A1 = " << ampts[0] <<
" A2 = " << ampts[1] <<
" A3 = " << ampts[2] <<
" A4 = " << ampts[3]
8508 std::cout <<
" A5 = " << ampts[4] <<
" A6 = " << ampts[5] <<
" A7 = " << ampts[6] <<
" A8 = " << ampts[7]
8511 std::cout <<
" icapid1 = " << icapid[0] <<
" icapid2 = " << icapid[1] <<
" icapid3 = " << icapid[2]
8512 <<
" icapid4 = " << icapid[3] << std::endl;
8513 std::cout <<
" icapid5 = " << icapid[4] <<
" icapid6 = " << icapid[5] <<
" icapid7 = " << icapid[6]
8514 <<
" icapid8 = " << icapid[7] << std::endl;
8516 if ((ampts[0] == ampts[1]) && (ampts[1] == ampts[2]) && (ampts[2] == ampts[3]) && (ampts[3] == ampts[4]) &&
8517 (ampts[4] == ampts[5]) && (ampts[5] == ampts[6]) && (ampts[6] == ampts[7])) {
8527 if (((ampts[0] == ampts[1]) && (ampts[1] == ampts[2]) && (ampts[2] == ampts[3]) && (ampts[3] == ampts[4]) &&
8528 (ampts[4] == ampts[5])) ||
8529 ((ampts[1] == ampts[2]) && (ampts[2] == ampts[3]) && (ampts[3] == ampts[4]) && (ampts[4] == ampts[5]) &&
8530 (ampts[5] == ampts[6])) ||
8531 ((ampts[2] == ampts[3]) && (ampts[3] == ampts[4]) && (ampts[4] == ampts[5]) && (ampts[5] == ampts[6]) &&
8532 (ampts[6] == ampts[7]))) {
8540 if (((ampts[0] == ampts[1]) && (ampts[1] == ampts[2]) && (ampts[2] == ampts[3])) ||
8541 ((ampts[1] == ampts[2]) && (ampts[2] == ampts[3]) && (ampts[3] == ampts[4])) ||
8542 ((ampts[2] == ampts[3]) && (ampts[3] == ampts[4]) && (ampts[4] == ampts[5])) ||
8543 ((ampts[3] == ampts[4]) && (ampts[4] == ampts[5]) && (ampts[5] == ampts[6]))) {
8551 if (((ampts[4] == ampts[5]) && (ampts[5] == ampts[6]) && (ampts[6] == ampts[7]))) {
8559 if (((ampts[0] == ampts[1]) && (ampts[1] == ampts[2])) || ((ampts[1] == ampts[2]) && (ampts[2] == ampts[3])) ||
8560 ((ampts[2] == ampts[3]) && (ampts[3] == ampts[4]))) {
8568 if (((ampts[3] == ampts[4]) && (ampts[4] == ampts[5])) || ((ampts[4] == ampts[5]) && (ampts[5] == ampts[6])) ||
8569 ((ampts[5] == ampts[6]) && (ampts[6] == ampts[7]))) {
8577 if (ampts[0] == ampts[1]) {
8585 if (ampts[1] == ampts[2]) {
8593 if (ampts[2] == ampts[3]) {
8601 if (ampts[3] == ampts[4]) {
8609 if (ampts[4] == ampts[5]) {
8617 if (ampts[5] == ampts[6]) {
8625 if (ampts[6] == ampts[7]) {
8637 std::cout <<
" A>139, sticking = " << sticking <<
" sub = " << sub <<
" ampts0 = " << ampts[0]
8638 <<
" ampts1 = " << ampts[1] <<
" ampts2 = " << ampts[2] <<
" ampts3 = " << ampts[3]
8639 <<
" ampts4 = " << ampts[4] <<
" ampts5 = " << ampts[5] <<
" ampts6 = " << ampts[6]
8640 <<
" ampts7 = " << ampts[7] << std::endl;
8641 std::cout <<
" icapid0 = " << icapid[0] <<
" icapid1 = " << icapid[1] <<
" icapid2 = " << icapid[2]
8642 <<
" icapid3 = " << icapid[3] <<
" icapid4 = " << icapid[4] <<
" icapid5 = " << icapid[5]
8643 <<
" icapid6 = " << icapid[6] <<
" icapid7 = " << icapid[7] << std::endl;
8652 if (error1 != 0 || error2 != 0 || error3 != 0) {
8673 if (error1 != 0 || error2 != 0 || error3 != 0) {
8703 int mdepth = cell.
depth();
8704 int iphi = cell.iphi() - 1;
8705 int ieta = cell.ieta();
8708 int sub = cell.subdet();
8710 std::cout <<
" HF DIGI ?????????????? ERROR mdepth = " << mdepth << std::endl;
8712 int errorGeneral = 0;
8721 bool anycapid =
true;
8726 int lastcapid = 0, capid = 0;
8727 int ERRORfiber = -10;
8728 int ERRORfiberChan = -10;
8729 int ERRORfiberAndChan = -10;
8730 int repetedcapid = 0;
8733 TSsize = digiItr->size();
8735 for (
int ii = 0;
ii < TSsize;
ii++) {
8736 capid = (*digiItr)[
ii].capid();
8737 bool er = (*digiItr)[
ii].er();
8738 bool dv = (*digiItr)[
ii].dv();
8739 int fiber = (*digiItr)[
ii].fiber();
8740 int fiberChan = (*digiItr)[
ii].fiberChan();
8741 int fiberAndChan = (*digiItr)[
ii].fiberAndChan();
8742 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8745 ERRORfiberChan = fiberChan;
8746 ERRORfiberAndChan = fiberAndChan;
8747 if (capid != lastcapid) {
8758 ERRORfiberChan = fiberChan;
8759 ERRORfiberAndChan = fiberAndChan;
8764 ERRORfiberChan = fiberChan;
8765 ERRORfiberAndChan = fiberAndChan;
8769 if (firstcapid == 0 && !anycapid)
8771 if (firstcapid == 1 && !anycapid)
8773 if (firstcapid == 2 && !anycapid)
8775 if (firstcapid == 3 && !anycapid)
8783 if (!anycapid && anyer)
8785 if (!anycapid && !anydv)
8787 if (!anycapid && anyer && !anydv)
8789 if (anyer && !anydv)
8793 for (
int ii = 0;
ii < TSsize;
ii++) {
8794 double ampldefault =
adc2fC[digiItr->sample(
ii).adc()];
8795 ampl += ampldefault;
8809 if (error1 != 0 || error2 != 0 || error3 != 0) {
8839 int mdepth = hcaldetid.
depth();
8840 int sub = hcaldetid.
subdet();
8844 bool anycapid =
true;
8849 int lastcapid = 0, capid = 0;
8850 int repetedcapid = 0;
8855 capid = qie10df[
ii].capid();
8856 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8858 if (capid != lastcapid) {
8873 int adc = qie10df[
ii].adc();
8875 ampl += ampldefault;
8906 int mdepth = cell.
depth();
8907 int iphi = cell.iphi() - 1;
8908 int ieta = cell.ieta();
8911 int sub = cell.subdet();
8912 int errorGeneral = 0;
8921 bool anycapid =
true;
8926 int lastcapid = 0, capid = 0;
8927 int ERRORfiber = -10;
8928 int ERRORfiberChan = -10;
8929 int ERRORfiberAndChan = -10;
8930 int repetedcapid = 0;
8931 for (
int ii = 0;
ii < (*digiItr).size();
ii++) {
8932 capid = (*digiItr)[
ii].capid();
8933 bool er = (*digiItr)[
ii].er();
8934 bool dv = (*digiItr)[
ii].dv();
8935 int fiber = (*digiItr)[
ii].fiber();
8936 int fiberChan = (*digiItr)[
ii].fiberChan();
8937 int fiberAndChan = (*digiItr)[
ii].fiberAndChan();
8938 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8941 ERRORfiberChan = fiberChan;
8942 ERRORfiberAndChan = fiberAndChan;
8943 if (capid != lastcapid) {
8956 ERRORfiberChan = fiberChan;
8957 ERRORfiberAndChan = fiberAndChan;
8962 ERRORfiberChan = fiberChan;
8963 ERRORfiberAndChan = fiberAndChan;
8969 if (firstcapid == 0 && !anycapid)
8971 if (firstcapid == 1 && !anycapid)
8973 if (firstcapid == 2 && !anycapid)
8975 if (firstcapid == 3 && !anycapid)
8984 if (!anycapid && anyer)
8986 if (!anycapid && !anydv)
8988 if (!anycapid && anyer && !anydv)
8990 if (anyer && !anydv)
8995 for (
int ii = 0;
ii < (*digiItr).size();
ii++) {
8996 double ampldefault =
adc2fC[digiItr->sample(
ii).adc()];
8997 ampl += ampldefault;
9013 if (error1 != 0 || error2 != 0 || error3 != 0) {
9022 std::cout <<
" mdepth HO = " << mdepth << std::endl;
9035 double toolwithPedSubtr[100];
9036 double lintoolwithoutPedSubtr[100];
9038 int mdepth = cell.
depth();
9039 int iphi = cell.iphi() - 1;
9040 int ieta0 = cell.ieta();
9044 int sub = cell.subdet();
9056 coder.
adc2fC(*digiItr, toolOriginal);
9057 double pedestalaver9 = 0.;
9058 double pedestalaver4 = 0.;
9059 double pedestal0 = 0.;
9060 double pedestal1 = 0.;
9061 double pedestal2 = 0.;
9062 double pedestal3 = 0.;
9063 double pedestalwaver9 = 0.;
9064 double pedestalwaver4 = 0.;
9065 double pedestalw0 = 0.;
9066 double pedestalw1 = 0.;
9067 double pedestalw2 = 0.;
9068 double pedestalw3 = 0.;
9069 double difpedestal0 = 0.;
9070 double difpedestal1 = 0.;
9071 double difpedestal2 = 0.;
9072 double difpedestal3 = 0.;
9075 double absamplitude = 0.;
9076 double amplitude345 = 0.;
9078 double linamplitudewithoutPedSubtr = 0.;
9081 double max_signal = -100.;
9082 int ts_with_max_signal = -100;
9088 double errorBtype = 0.;
9091 TSsize = digiItr->size();
9092 if ((*digiItr).size() != 10)
9095 for (
int ii = 0;
ii < TSsize;
ii++) {
9097 double ampldefaultwithPedSubtr = 0.;
9098 double linampldefaultwithoutPedSubtr = 0.;
9099 double ampldefault = 0.;
9100 double ampldefault0 = 0.;
9101 double ampldefault1 = 0.;
9102 double ampldefault2 = 0.;
9103 ampldefault0 =
adc2fC[digiItr->sample(
ii).adc()];
9105 ampldefault1 = toolOriginal[
ii];
9106 ampldefault2 = (*digiItr)[
ii].adc();
9108 ampldefault = ampldefault0;
9111 ampldefault = ampldefault1;
9114 ampldefault = ampldefault2;
9116 ampldefaultwithPedSubtr = ampldefault0;
9117 linampldefaultwithoutPedSubtr = ampldefault2;
9119 int capid = ((*digiItr)[
ii]).capid();
9121 double pedestalINI = pedestal00->
getValue(capid);
9123 double pedestalw = pedw->
getSigma(capid, capid);
9124 ampldefaultwithPedSubtr -=
pedestal;
9128 tool[
ii] = ampldefault;
9129 toolwithPedSubtr[
ii] = ampldefaultwithPedSubtr;
9130 lintoolwithoutPedSubtr[
ii] = linampldefaultwithoutPedSubtr;
9133 pedestalwaver9 += pedestalw * pedestalw;
9135 if (capid == 0 &&
c0 == 0) {
9140 pedestalw0 = pedestalw;
9141 pedestalwaver4 += pedestalw * pedestalw;
9142 difpedestal0 =
pedestal - pedestalINI;
9145 if (capid == 1 &&
c1 == 0) {
9150 pedestalw1 = pedestalw;
9151 pedestalwaver4 += pedestalw * pedestalw;
9152 difpedestal1 =
pedestal - pedestalINI;
9154 if (capid == 2 && c2 == 0) {
9159 pedestalw2 = pedestalw;
9160 pedestalwaver4 += pedestalw * pedestalw;
9161 difpedestal2 =
pedestal - pedestalINI;
9163 if (capid == 3 && c3 == 0) {
9168 pedestalw3 = pedestalw;
9169 pedestalwaver4 += pedestalw * pedestalw;
9170 difpedestal3 =
pedestal - pedestalINI;
9173 if (max_signal < ampldefault) {
9174 max_signal = ampldefault;
9175 ts_with_max_signal =
ii;
9190 absamplitude +=
abs(ampldefault);
9192 if (
ii == 3 ||
ii == 4 ||
ii == 5)
9193 amplitude345 += ampldefault;
9197 timew += (
ii + 1) *
abs(ampldefault);
9198 timeww += (
ii + 1) * ampldefault;
9202 if (sub == 1 || sub == 2) {
9207 pedestalaver9 /= TSsize;
9208 pedestalaver4 /= c4;
9209 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
9210 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
9211 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
9212 ampl = tool[ts_with_max_signal];
9213 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
9214 ampl += tool[ts_with_max_signal + 2];
9215 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
9216 ampl += tool[ts_with_max_signal + 1];
9217 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
9218 ampl += tool[ts_with_max_signal - 1];
9221 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
9222 linamplitudewithoutPedSubtr = lintoolwithoutPedSubtr[ts_with_max_signal];
9223 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
9224 linamplitudewithoutPedSubtr += lintoolwithoutPedSubtr[ts_with_max_signal + 2];
9225 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
9226 linamplitudewithoutPedSubtr += lintoolwithoutPedSubtr[ts_with_max_signal + 1];
9227 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
9228 linamplitudewithoutPedSubtr += lintoolwithoutPedSubtr[ts_with_max_signal - 1];
9233 if (ratio < 0. || ratio > 1.02)
9235 double aveamplitude = 0.;
9236 double aveamplitudew = 0.;
9237 if (absamplitude > 0 && timew > 0)
9238 aveamplitude = timew / absamplitude;
9244 bool anycapid =
true;
9249 for (
int ii = 0;
ii < TSsize;
ii++) {
9250 double aaaaaa = (
ii + 1) - aveamplitudew;
9251 double aaaaaa2 = aaaaaa * aaaaaa;
9252 double ampldefault = tool[
ii];
9253 rmsamp += (aaaaaa2 * ampldefault);
9254 capid = ((*digiItr)[
ii]).capid();
9255 bool er = (*digiItr)[
ii].er();
9256 bool dv = (*digiItr)[
ii].dv();
9257 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
9269 if (!anycapid || anyer || !anydv)
9272 double rmsamplitude = 0.;
9275 double aveamplitude1 = aveamplitude - 1;
9281 if (sub == 1 && mdepth == 1)
9283 if (sub == 1 && mdepth == 2)
9285 if (sub == 2 && mdepth == 1)
9287 if (sub == 2 && mdepth == 2)
9289 if (sub == 2 && mdepth == 3)
9293 if (sub == 1 && mdepth == 1)
9295 if (sub == 1 && mdepth == 2)
9297 if (sub == 2 && mdepth == 1)
9299 if (sub == 2 && mdepth == 2)
9301 if (sub == 2 && mdepth == 3)
9305 for (
int ii = 0;
ii < TSsize;
ii++) {
9307 double ampldefault = tool[
ii];
9353 for (
int ii = 0;
ii < TSsize;
ii++) {
9354 double ampldefault = tool[
ii];
9361 for (
int ii = 0;
ii < TSsize;
ii++) {
9362 double ampldefault = tool[
ii];
9374 for (
int ii = 0;
ii < TSsize;
ii++) {
9375 double ampldefault = tool[
ii];
9382 for (
int ii = 0;
ii < TSsize;
ii++) {
9383 double ampldefault = tool[
ii];
9399 double mypedestal = pedestal0;
9400 double mypedestalw = pedestalw0;
9471 for (
int ii = 0;
ii < TSsize;
ii++) {
9504 if (amplitude < ADCAmplHBMin_ || amplitude >
ADCAmplHBMax_) {
9544 if (aveamplitude1 < TSmeanHBMin_ || aveamplitude1 >
TSmeanHBMax_) {
9565 if (ts_with_max_signal < TSpeakHBMin_ || ts_with_max_signal >
TSpeakHBMax_) {
9586 if (rmsamplitude < rmsHBMin_ || rmsamplitude >
rmsHBMax_) {
9645 double mypedestal = pedestal0;
9646 double mypedestalw = pedestalw0;
9717 for (
int ii = 0;
ii < TSsize;
ii++) {
9739 if (amplitude < ADCAmplHEMin_ || amplitude >
ADCAmplHEMax_) {
9800 if (aveamplitude1 < TSmeanHEMin_ || aveamplitude1 >
TSmeanHEMax_) {
9824 if (ts_with_max_signal < TSpeakHEMin_ || ts_with_max_signal >
TSpeakHEMax_) {
9848 if (rmsamplitude < rmsHEMin_ || rmsamplitude >
rmsHEMax_) {
9908 double amplitudewithPedSubtr = 0.;
9911 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
9912 amplitudewithPedSubtr = toolwithPedSubtr[ts_with_max_signal];
9913 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
9914 amplitudewithPedSubtr += toolwithPedSubtr[ts_with_max_signal + 2];
9915 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
9916 amplitudewithPedSubtr += toolwithPedSubtr[ts_with_max_signal + 1];
9917 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
9918 amplitudewithPedSubtr += toolwithPedSubtr[ts_with_max_signal - 1];
9924 if (amplitudewithPedSubtr > 50.) {
9947 int LLLLLL111111 = 0;
9948 if ((mdepth == 1 && fabs(ieta0) > 16 && fabs(ieta0) < 30))
9951 int LLLLLL222222 = 0;
9952 if ((mdepth == 2 && fabs(ieta0) > 16 && fabs(ieta0) < 27) ||
9953 (mdepth == 3 && fabs(ieta0) > 26 && fabs(ieta0) < 29))
9956 if (LLLLLL111111 == 1) {
9961 if (LLLLLL222222 == 1) {
9968 if (mdepth == 3 && fabs(ieta0) == 16) {
9997 int mdepth = hcaldetid.
depth();
9998 int sub = hcaldetid.
subdet();
10001 if (mdepth == 0 || sub > 4)
10026 coder.
adc2fC(qie11df, toolOriginal);
10027 double pedestalaver9 = 0.;
10028 double pedestalaver4 = 0.;
10029 double pedestal0 = 0.;
10030 double pedestal1 = 0.;
10031 double pedestal2 = 0.;
10032 double pedestal3 = 0.;
10033 double pedestalwaver9 = 0.;
10034 double pedestalwaver4 = 0.;
10035 double pedestalw0 = 0.;
10036 double pedestalw1 = 0.;
10037 double pedestalw2 = 0.;
10038 double pedestalw3 = 0.;
10039 double difpedestal0 = 0.;
10040 double difpedestal1 = 0.;
10041 double difpedestal2 = 0.;
10042 double difpedestal3 = 0.;
10045 double amplitude0 = 0.;
10046 double absamplitude = 0.;
10047 double tocampl = 0.;
10049 double amplitude345 = 0.;
10051 double ampl3ts = 0.;
10053 double timeww = 0.;
10054 double max_signal = -100.;
10055 int ts_with_max_signal = -100;
10061 double errorBtype = 0.;
10072 for (
int ii = 0;
ii < TSsize;
ii++) {
10073 double ampldefault = 0.;
10074 double tocdefault = 0.;
10075 double ampldefault0 = 0.;
10076 double ampldefault1 = 0.;
10077 double ampldefault2 = 0.;
10083 ampldefault1 = toolOriginal[
ii];
10084 ampldefault2 = qie11df[
ii].adc();
10087 ampldefault = ampldefault0;
10090 ampldefault = ampldefault1;
10093 ampldefault = ampldefault2;
10095 tocdefault = ampldefault;
10097 int capid = (qie11df[
ii]).capid();
10099 double pedestalw = pedw->
getSigma(capid, capid);
10100 double pedestalINI = pedestal00->
getValue(capid);
10104 tool[
ii] = ampldefault;
10106 pedestalwaver9 += pedestalw * pedestalw;
10108 if (capid == 0 &&
c0 == 0) {
10113 pedestalw0 = pedestalw;
10114 pedestalwaver4 += pedestalw * pedestalw;
10115 difpedestal0 =
pedestal - pedestalINI;
10118 if (capid == 1 &&
c1 == 0) {
10123 pedestalw1 = pedestalw;
10124 pedestalwaver4 += pedestalw * pedestalw;
10125 difpedestal1 =
pedestal - pedestalINI;
10127 if (capid == 2 && c2 == 0) {
10132 pedestalw2 = pedestalw;
10133 pedestalwaver4 += pedestalw * pedestalw;
10134 difpedestal2 =
pedestal - pedestalINI;
10136 if (capid == 3 && c3 == 0) {
10141 pedestalw3 = pedestalw;
10142 pedestalwaver4 += pedestalw * pedestalw;
10143 difpedestal3 =
pedestal - pedestalINI;
10146 if (max_signal < ampldefault) {
10147 max_signal = ampldefault;
10148 ts_with_max_signal =
ii;
10151 absamplitude +=
abs(ampldefault);
10152 tocampl += tocdefault;
10154 if (
ii == 1 ||
ii == 2 ||
ii == 3 ||
ii == 4 ||
ii == 5 ||
ii == 6 ||
ii == 7 ||
ii == 8)
10155 amplitude345 += ampldefault;
10160 timew += (
ii + 1) *
abs(ampldefault);
10161 timeww += (
ii + 1) * ampldefault;
10163 if (
ii == 2 && ampldefault > 0.)
10169 pedestalaver9 /= TSsize;
10170 pedestalaver4 /= c4;
10171 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
10172 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
10175 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize) {
10176 ampl = tool[ts_with_max_signal];
10177 ampl3ts = tool[ts_with_max_signal];
10179 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize) {
10180 ampl += tool[ts_with_max_signal - 1];
10181 ampl3ts += tool[ts_with_max_signal - 1];
10183 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize) {
10184 ampl += tool[ts_with_max_signal + 1];
10185 ampl3ts += tool[ts_with_max_signal + 1];
10187 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize) {
10188 ampl += tool[ts_with_max_signal + 2];
10194 double xb = ampl / 40.;
10195 double xc = amplitude345 / 40.;
10196 double xd = ampl3ts / 40.;
10197 double txa = tocampl / 40.;
10199 if (((
ieta == -16 ||
ieta == 15) && mdepth == 4) ||
10200 ((
ieta == -17 ||
ieta == 16) && (mdepth == 2 || mdepth == 3)) ||
10201 ((
ieta == -18 ||
ieta == 17) && mdepth == 5)) {
10202 double c0 = 1.000000;
10203 double b1 = 2.59096e-05;
10204 double a2 = 4.60721e-11;
10205 double corrforxa =
a2 * xa * xa +
b1 * xa +
c0;
10206 double corrforxb =
a2 * xb * xb +
b1 * xb +
c0;
10207 double corrforxc =
a2 * xc * xc +
b1 * xc +
c0;
10208 double corrforxd =
a2 * xd * xd +
b1 * xd +
c0;
10209 double corrfortxa =
a2 * txa * txa +
b1 * txa +
c0;
10214 amplitude345 *= corrforxc;
10215 ampl3ts *= corrforxd;
10216 tocampl *= corrfortxa;
10220 double c0 = 1.000000;
10221 double b1 = 2.71238e-05;
10222 double a2 = 1.32877e-10;
10223 double corrforxa =
a2 * xa * xa +
b1 * xa +
c0;
10224 double corrforxb =
a2 * xb * xb +
b1 * xb +
c0;
10225 double corrforxc =
a2 * xc * xc +
b1 * xc +
c0;
10226 double corrforxd =
a2 * xd * xd +
b1 * xd +
c0;
10227 double corrfortxa =
a2 * txa * txa +
b1 * txa +
c0;
10232 amplitude345 *= corrforxc;
10233 ampl3ts *= corrforxd;
10234 tocampl *= corrfortxa;
10240 if (sub == 1 || sub == 2) {
10251 if (ratio < 0. || ratio > 1.02)
10253 double aveamplitude = 0.;
10254 double aveamplitudew = 0.;
10255 if (absamplitude > 0 && timew > 0)
10256 aveamplitude = timew / absamplitude;
10257 if (amplitude0 > 0 && timeww > 0)
10258 aveamplitudew = timeww / amplitude0;
10259 double rmsamp = 0.;
10262 bool anycapid =
true;
10265 for (
int ii = 0;
ii < TSsize;
ii++) {
10266 double aaaaaa = (
ii + 1) - aveamplitudew;
10267 double aaaaaa2 = aaaaaa * aaaaaa;
10268 double ampldefault = tool[
ii];
10269 rmsamp += (aaaaaa2 * ampldefault);
10270 capid = (qie11df[
ii]).capid();
10271 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
10279 double rmsamplitude = 0.;
10280 if ((amplitude0 > 0 && rmsamp > 0) || (amplitude0 < 0 && rmsamp < 0))
10281 rmsamplitude =
sqrt(rmsamp / amplitude0);
10282 double aveamplitude1 = aveamplitude - 1;
10287 if (sub == 1 && mdepth == 1)
10289 if (sub == 1 && mdepth == 2)
10291 if (sub == 2 && mdepth == 1)
10293 if (sub == 2 && mdepth == 2)
10295 if (sub == 2 && mdepth == 3)
10299 if (sub == 1 && mdepth == 1)
10301 if (sub == 1 && mdepth == 2)
10303 if (sub == 2 && mdepth == 1)
10305 if (sub == 2 && mdepth == 2)
10307 if (sub == 2 && mdepth == 3)
10311 for (
int ii = 0;
ii < TSsize;
ii++) {
10313 double ampldefault = tool[
ii];
10316 if (amplitude0 > 120) {
10359 for (
int ii = 0;
ii < TSsize;
ii++) {
10360 double ampldefault = tool[
ii];
10367 for (
int ii = 0;
ii < TSsize;
ii++) {
10368 double ampldefault = tool[
ii];
10382 for (
int ii = 0;
ii < TSsize;
ii++) {
10383 double ampldefault = tool[
ii];
10390 for (
int ii = 0;
ii < TSsize;
ii++) {
10391 double ampldefault = tool[
ii];
10407 double mypedestal = pedestal0;
10408 double mypedestalw = pedestalw0;
10479 for (
int ii = 0;
ii < TSsize;
ii++) {
10512 if (amplitude < ADCAmplHBMin_ || amplitude >
ADCAmplHBMax_) {
10570 if (aveamplitude1 < TSmeanHBMin_ || aveamplitude1 >
TSmeanHBMax_) {
10598 h_TSmaxA_HB->Fill(
float(ts_with_max_signal), 1.);
10599 if (ts_with_max_signal < TSpeakHBMin_ || ts_with_max_signal >
TSpeakHBMax_) {
10628 if (rmsamplitude < rmsHBMin_ || rmsamplitude >
rmsHBMax_) {
10657 if (ratio < ratioHBMin_ || ratio >
ratioHBMax_) {
10706 if (flagTS2 == 1) {
10725 double mypedestal = pedestal0;
10726 double mypedestalw = pedestalw0;
10813 for (
int ii = 0;
ii < TSsize;
ii++) {
10831 if (amplitude < ADCAmplHEMin_ || amplitude >
ADCAmplHEMax_) {
10879 if (ampl3ts > 70 && ampl3ts < 110) {
10882 }
else if (ampl3ts > 110 && ampl3ts < 150) {
10886 if (ampl > 87 && ampl < 127) {
10889 }
else if (ampl > 127 && ampl < 167) {
10952 if (aveamplitude1 < TSmeanHEMin_ || aveamplitude1 >
TSmeanHEMax_) {
10991 h_TSmaxA_HE->Fill(
float(ts_with_max_signal), 1.);
10992 if (ts_with_max_signal < TSpeakHEMin_ || ts_with_max_signal >
TSpeakHEMax_) {
11032 if (rmsamplitude < rmsHEMin_ || rmsamplitude >
rmsHEMax_) {
11073 if (ratio < ratioHEMin_ || ratio >
ratioHEMax_) {
11143 if (flagTS2 == 1) {
11168 int mdepth = cell.
depth();
11169 int iphi = cell.iphi() - 1;
11170 int ieta = cell.ieta();
11173 int sub = cell.subdet();
11183 coder.
adc2fC(*digiItr, toolOriginal);
11184 double pedestalaver9 = 0.;
11185 double pedestalaver4 = 0.;
11186 double pedestal0 = 0.;
11187 double pedestal1 = 0.;
11188 double pedestal2 = 0.;
11189 double pedestal3 = 0.;
11190 double pedestalwaver9 = 0.;
11191 double pedestalwaver4 = 0.;
11192 double pedestalw0 = 0.;
11193 double pedestalw1 = 0.;
11194 double pedestalw2 = 0.;
11195 double pedestalw3 = 0.;
11197 double absamplitude = 0.;
11200 double timeww = 0.;
11201 double max_signal = -100.;
11202 int ts_with_max_signal = -100;
11208 double errorBtype = 0.;
11210 if ((*digiItr).size() != TSsize)
11212 TSsize = digiItr->size();
11213 for (
int ii = 0;
ii < TSsize;
ii++) {
11215 double ampldefault = 0.;
11216 double ampldefault0 = 0.;
11217 double ampldefault1 = 0.;
11218 double ampldefault2 = 0.;
11219 ampldefault0 =
adc2fC[digiItr->sample(
ii).adc()];
11221 ampldefault1 = toolOriginal[
ii];
11222 ampldefault2 = (*digiItr)[
ii].adc();
11224 ampldefault = ampldefault0;
11227 ampldefault = ampldefault1;
11230 ampldefault = ampldefault2;
11233 int capid = ((*digiItr)[
ii]).capid();
11236 double pedestalw = pedw->
getSigma(capid, capid);
11240 tool[
ii] = ampldefault;
11243 pedestalwaver9 += pedestalw * pedestalw;
11245 if (capid == 0 &&
c0 == 0) {
11250 pedestalw0 = pedestalw;
11251 pedestalwaver4 += pedestalw * pedestalw;
11254 if (capid == 1 &&
c1 == 0) {
11259 pedestalw1 = pedestalw;
11260 pedestalwaver4 += pedestalw * pedestalw;
11262 if (capid == 2 && c2 == 0) {
11267 pedestalw2 = pedestalw;
11268 pedestalwaver4 += pedestalw * pedestalw;
11270 if (capid == 3 && c3 == 0) {
11275 pedestalw3 = pedestalw;
11276 pedestalwaver4 += pedestalw * pedestalw;
11279 if (max_signal < ampldefault) {
11280 max_signal = ampldefault;
11281 ts_with_max_signal =
ii;
11284 absamplitude +=
abs(ampldefault);
11289 timew += (
ii + 1) *
abs(ampldefault);
11290 timeww += (
ii + 1) * ampldefault;
11300 pedestalaver9 /= TSsize;
11301 pedestalaver4 /= c4;
11302 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
11303 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
11306 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
11307 ampl = tool[ts_with_max_signal];
11308 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
11309 ampl += tool[ts_with_max_signal + 2];
11310 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
11311 ampl += tool[ts_with_max_signal + 1];
11312 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
11313 ampl += tool[ts_with_max_signal - 1];
11320 if (ratio < 0. || ratio > 1.02)
11323 double aveamplitude = 0.;
11324 double aveamplitudew = 0.;
11325 if (absamplitude > 0 && timew > 0)
11326 aveamplitude = timew / absamplitude;
11330 double rmsamp = 0.;
11333 bool anycapid =
true;
11334 bool anyer =
false;
11338 for (
int ii = 0;
ii < TSsize;
ii++) {
11339 double aaaaaa = (
ii + 1) - aveamplitudew;
11340 double aaaaaa2 = aaaaaa * aaaaaa;
11341 double ampldefault = tool[
ii];
11342 rmsamp += (aaaaaa2 * ampldefault);
11343 capid = ((*digiItr)[
ii]).capid();
11344 bool er = (*digiItr)[
ii].er();
11345 bool dv = (*digiItr)[
ii].dv();
11346 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
11359 if (!anycapid || anyer || !anydv)
11361 double rmsamplitude = 0.;
11364 double aveamplitude1 = aveamplitude - 1;
11366 if (sub == 4 && mdepth == 1)
11368 if (sub == 4 && mdepth == 2)
11372 if (sub == 4 && mdepth == 1)
11374 if (sub == 4 && mdepth == 2)
11386 for (
int ii = 0;
ii < TSsize;
ii++) {
11387 double ampldefault = tool[
ii];
11394 for (
int ii = 0;
ii < TSsize;
ii++) {
11395 double ampldefault = tool[
ii];
11409 double mypedestal = pedestal0;
11410 double mypedestalw = pedestalw0;
11470 for (
int ii = 0;
ii < TSsize;
ii++) {
11487 if (amplitude < ADCAmplHFMin_ || amplitude >
ADCAmplHFMax_) {
11529 if (aveamplitude1 < TSmeanHFMin_ || aveamplitude1 >
TSmeanHFMax_) {
11549 h_TSmaxA_HF->Fill(
float(ts_with_max_signal), 1.);
11550 if (ts_with_max_signal < TSpeakHFMin_ || ts_with_max_signal >
TSpeakHFMax_) {
11570 if (rmsamplitude < rmsHFMin_ || rmsamplitude >
rmsHFMax_) {
11590 if (ratio < ratioHFMin_ || ratio >
ratioHFMax_) {
11638 int mdepth = hcaldetid.
depth();
11639 int sub = hcaldetid.
subdet();
11654 if (mdepth == 0 || sub != 4)
11674 coder.
adc2fC(qie10df, toolOriginal);
11677 double pedestalaver9 = 0.;
11678 double pedestalaver4 = 0.;
11679 double pedestal0 = 0.;
11680 double pedestal1 = 0.;
11681 double pedestal2 = 0.;
11682 double pedestal3 = 0.;
11683 double pedestalwaver9 = 0.;
11684 double pedestalwaver4 = 0.;
11685 double pedestalw0 = 0.;
11686 double pedestalw1 = 0.;
11687 double pedestalw2 = 0.;
11688 double pedestalw3 = 0.;
11690 double absamplitude = 0.;
11693 double timeww = 0.;
11694 double max_signal = -100.;
11695 int ts_with_max_signal = -100;
11701 double errorBtype = 0.;
11708 for (
int ii = 0;
ii < TSsize;
ii++) {
11709 double ampldefault = 0.;
11710 double ampldefault0 = 0.;
11711 double ampldefault1 = 0.;
11712 double ampldefault2 = 0.;
11715 ampldefault1 = toolOriginal[
ii];
11716 ampldefault2 = qie10df[
ii].adc();
11718 ampldefault = ampldefault0;
11721 ampldefault = ampldefault1;
11724 ampldefault = ampldefault2;
11727 int capid = (qie10df[
ii]).capid();
11729 double pedestalw = pedw->
getSigma(capid, capid);
11734 tool[
ii] = ampldefault;
11737 pedestalwaver9 += pedestalw * pedestalw;
11739 if (capid == 0 &&
c0 == 0) {
11744 pedestalw0 = pedestalw;
11745 pedestalwaver4 += pedestalw * pedestalw;
11748 if (capid == 1 &&
c1 == 0) {
11753 pedestalw1 = pedestalw;
11754 pedestalwaver4 += pedestalw * pedestalw;
11756 if (capid == 2 && c2 == 0) {
11761 pedestalw2 = pedestalw;
11762 pedestalwaver4 += pedestalw * pedestalw;
11764 if (capid == 3 && c3 == 0) {
11769 pedestalw3 = pedestalw;
11770 pedestalwaver4 += pedestalw * pedestalw;
11773 if (max_signal < ampldefault) {
11774 max_signal = ampldefault;
11775 ts_with_max_signal =
ii;
11778 absamplitude +=
abs(ampldefault);
11780 timew += (
ii + 1) *
abs(ampldefault);
11781 timeww += (
ii + 1) * ampldefault;
11782 if (
ii == 1 && ampldefault > 0.)
11793 pedestalaver9 /= TSsize;
11794 pedestalaver4 /= c4;
11795 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
11796 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
11799 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
11800 ampl = tool[ts_with_max_signal];
11801 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
11802 ampl += tool[ts_with_max_signal + 2];
11803 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
11804 ampl += tool[ts_with_max_signal + 1];
11805 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
11806 ampl += tool[ts_with_max_signal - 1];
11812 if (ratio < 0. || ratio > 1.02)
11814 double aveamplitude = 0.;
11815 double aveamplitudew = 0.;
11816 if (absamplitude > 0 && timew > 0)
11817 aveamplitude = timew / absamplitude;
11821 double rmsamp = 0.;
11823 bool anycapid =
true;
11826 for (
int ii = 0;
ii < TSsize;
ii++) {
11827 double aaaaaa = (
ii + 1) - aveamplitudew;
11828 double aaaaaa2 = aaaaaa * aaaaaa;
11829 double ampldefault = tool[
ii];
11830 rmsamp += (aaaaaa2 * ampldefault);
11831 capid = (qie10df[
ii]).capid();
11832 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
11840 double rmsamplitude = 0.;
11843 double aveamplitude1 = aveamplitude - 1;
11848 if (sub == 4 && mdepth == 1)
11850 if (sub == 4 && mdepth == 2)
11854 if (sub == 4 && mdepth == 1)
11856 if (sub == 4 && mdepth == 2)
11868 for (
int ii = 0;
ii < TSsize;
ii++) {
11869 double ampldefault = tool[
ii];
11876 for (
int ii = 0;
ii < TSsize;
ii++) {
11877 double ampldefault = tool[
ii];
11892 double mypedestal = pedestal0;
11893 double mypedestalw = pedestalw0;
11961 for (
int ii = 0;
ii < TSsize;
ii++) {
11978 if (amplitude < ADCAmplHFMin_ || amplitude >
ADCAmplHFMax_) {
12034 if (aveamplitude1 < TSmeanHFMin_ || aveamplitude1 >
TSmeanHFMax_) {
12062 h_TSmaxA_HF->Fill(
float(ts_with_max_signal), 1.);
12063 if (ts_with_max_signal < TSpeakHFMin_ || ts_with_max_signal >
TSpeakHFMax_) {
12092 if (rmsamplitude < rmsHFMin_ || rmsamplitude >
rmsHFMax_) {
12121 if (ratio < ratioHFMin_ || ratio >
ratioHFMax_) {
12170 if (flagTS1 == 1) {
12191 int mdepth = cell.
depth();
12192 int iphi = cell.iphi() - 1;
12193 int ieta = cell.ieta();
12196 int sub = cell.subdet();
12206 coder.
adc2fC(*digiItr, toolOriginal);
12207 double pedestalaver9 = 0.;
12208 double pedestalaver4 = 0.;
12209 double pedestal0 = 0.;
12210 double pedestal1 = 0.;
12211 double pedestal2 = 0.;
12212 double pedestal3 = 0.;
12213 double pedestalwaver9 = 0.;
12214 double pedestalwaver4 = 0.;
12215 double pedestalw0 = 0.;
12216 double pedestalw1 = 0.;
12217 double pedestalw2 = 0.;
12218 double pedestalw3 = 0.;
12220 double absamplitude = 0.;
12223 double timeww = 0.;
12224 double max_signal = -100.;
12225 int ts_with_max_signal = -100;
12231 double errorBtype = 0.;
12233 if ((*digiItr).size() != TSsize)
12235 TSsize = digiItr->size();
12237 for (
int ii = 0;
ii < TSsize;
ii++) {
12238 double ampldefault = 0.;
12239 double ampldefault0 = 0.;
12240 double ampldefault1 = 0.;
12241 double ampldefault2 = 0.;
12242 ampldefault0 =
adc2fC[digiItr->sample(
ii).adc()];
12244 ampldefault1 = toolOriginal[
ii];
12245 ampldefault2 = (*digiItr)[
ii].adc();
12247 ampldefault = ampldefault0;
12250 ampldefault = ampldefault1;
12253 ampldefault = ampldefault2;
12255 int capid = ((*digiItr)[
ii]).capid();
12257 double pedestalw = pedw->
getSigma(capid, capid);
12260 tool[
ii] = ampldefault;
12262 pedestalwaver9 += pedestalw * pedestalw;
12263 if (capid == 0 &&
c0 == 0) {
12268 pedestalw0 = pedestalw;
12269 pedestalwaver4 += pedestalw * pedestalw;
12271 if (capid == 1 &&
c1 == 0) {
12276 pedestalw1 = pedestalw;
12277 pedestalwaver4 += pedestalw * pedestalw;
12279 if (capid == 2 && c2 == 0) {
12284 pedestalw2 = pedestalw;
12285 pedestalwaver4 += pedestalw * pedestalw;
12287 if (capid == 3 && c3 == 0) {
12292 pedestalw3 = pedestalw;
12293 pedestalwaver4 += pedestalw * pedestalw;
12296 if (max_signal < ampldefault) {
12297 max_signal = ampldefault;
12298 ts_with_max_signal =
ii;
12301 absamplitude +=
abs(ampldefault);
12305 timew += (
ii + 1) *
abs(ampldefault);
12306 timeww += (
ii + 1) * ampldefault;
12307 if (ii < 3 && ampldefault > 0.)
12312 pedestalaver9 /= TSsize;
12313 pedestalaver4 /= c4;
12314 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
12315 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
12316 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
12317 ampl = tool[ts_with_max_signal];
12318 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
12319 ampl += tool[ts_with_max_signal + 2];
12320 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
12321 ampl += tool[ts_with_max_signal + 1];
12322 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
12323 ampl += tool[ts_with_max_signal - 1];
12327 if (ratio < 0. || ratio > 1.04)
12329 double aveamplitude = 0.;
12330 double aveamplitudew = 0.;
12331 if (absamplitude > 0 && timew > 0)
12332 aveamplitude = timew / absamplitude;
12335 double rmsamp = 0.;
12337 bool anycapid =
true;
12338 bool anyer =
false;
12342 for (
int ii = 0;
ii < TSsize;
ii++) {
12343 double aaaaaa = (
ii + 1) - aveamplitudew;
12344 double aaaaaa2 = aaaaaa * aaaaaa;
12345 double ampldefault = tool[
ii];
12346 rmsamp += (aaaaaa2 * ampldefault);
12347 capid = ((*digiItr)[
ii]).capid();
12348 bool er = (*digiItr)[
ii].er();
12349 bool dv = (*digiItr)[
ii].dv();
12350 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
12362 if (!anycapid || anyer || !anydv)
12364 double rmsamplitude = 0.;
12367 double aveamplitude1 = aveamplitude - 1;
12369 if (sub == 3 && mdepth == 4)
12373 if (sub == 3 && mdepth == 4)
12384 for (
int ii = 0;
ii < TSsize;
ii++) {
12385 double ampldefault = tool[
ii];
12390 for (
int ii = 0;
ii < TSsize;
ii++) {
12391 double ampldefault = tool[
ii];
12405 double mypedestal = pedestal0;
12406 double mypedestalw = pedestalw0;
12449 for (
int ii = 0;
ii < TSsize;
ii++) {
12465 if (amplitude < ADCAmplHOMin_ || amplitude >
ADCAmplHOMax_) {
12497 if (aveamplitude1 < TSmeanHOMin_ || aveamplitude1 >
TSmeanHOMax_) {
12509 h_TSmaxA_HO->Fill(
float(ts_with_max_signal), 1.);
12510 if (ts_with_max_signal < TSpeakHOMin_ || ts_with_max_signal >
TSpeakHOMax_) {
12524 if (rmsamplitude < rmsHOMin_ || rmsamplitude >
rmsHOMax_) {
12537 if (ratio < ratioHOMin_ || ratio >
ratioHOMax_) {
12554 if (flagTS012 == 1)
12561 if (kdet == 1 || kdet == 2) {
12565 cal_RBX = (kphi + 1) / 4;
12566 cal_RBX = cal_RBX + 18 * (keta + 1) / 2;
12569 cal_RBX = (
int)(kphi / 18) + 1;
12576 cal_RBX = kphi / 12 + 1;
12582 cal_RBX = kphi / 12 + 1 + 6;
12588 cal_RBX = kphi / 6 + 1 + 12;
12594 cal_RBX = kphi / 12 + 1 + 24;
12600 cal_RBX = kphi / 12 + 1 + 30;
12621 for (
int k1 = 0; k1 <
ndepth; k1++) {
12622 for (
int k2 = 0; k2 <
neta; k2++) {
12623 for (
int k3 = 0; k3 <
nphi; k3++) {
12624 int ieta = k2 - 41;
12719 int kkkk2 = (k2 - 1) / 4;
12724 int kkkk3 = (k3) / 4 + 1;
12726 ietaphi = ((kkkk2)-1) *
znphi + (kkkk3);
12729 bbb3 = bbbc / bbb1;
13573 std::cout <<
" ==== Edn of run " << std::endl;
13591 MAPfile <<
"#define LogEleMapdb_h" << std::endl;
13592 MAPfile <<
"#include <algorithm>" << std::endl;
13593 MAPfile <<
"#include <iostream>" << std::endl;
13594 MAPfile <<
"#include <vector>" << std::endl;
13595 MAPfile <<
"#include <string>" << std::endl;
13596 MAPfile <<
"#include <sstream>" << std::endl;
13599 MAPfile <<
"struct Cell {" << std::endl;
13600 MAPfile <<
" std::string subdet;" << std::endl;
13601 MAPfile <<
" int Eta;" << std::endl;
13602 MAPfile <<
" int Phi;" << std::endl;
13603 MAPfile <<
" int Depth;" << std::endl;
13604 MAPfile <<
" std::string RBX;" << std::endl;
13605 MAPfile <<
" int RM;" << std::endl;
13606 MAPfile <<
" int Pixel;" << std::endl;
13607 MAPfile <<
" int RMfiber;" << std::endl;
13608 MAPfile <<
" int FiberCh;" << std::endl;
13609 MAPfile <<
" int QIE;" << std::endl;
13610 MAPfile <<
" int ADC;" << std::endl;
13611 MAPfile <<
" int VMECardID;" << std::endl;
13612 MAPfile <<
" int dccID;" << std::endl;
13613 MAPfile <<
" int Spigot;" << std::endl;
13614 MAPfile <<
" int FiberIndex;" << std::endl;
13615 MAPfile <<
" int HtrSlot;" << std::endl;
13616 MAPfile <<
" int HtrTB;" << std::endl;
13619 MAPfile <<
"// the function check, if \"par\" == \"val\" for this cell" << std::endl;
13620 MAPfile <<
" bool check(const std::string par, const int val) const " << std::endl;
13621 MAPfile <<
" {" << std::endl;
13622 MAPfile <<
" if (par == \"Eta\") return (val == Eta);" << std::endl;
13623 MAPfile <<
" else if (par == \"Phi\") return (val == Phi);" << std::endl;
13624 MAPfile <<
" else if (par == \"Depth\") return (val == Depth);" << std::endl;
13625 MAPfile <<
" else if (par == \"RM\") return (val == RM);" << std::endl;
13626 MAPfile <<
" else if (par == \"Pixel\") return (val == Pixel);" << std::endl;
13627 MAPfile <<
" else if (par == \"RMfiber\") return (val == RMfiber);" << std::endl;
13628 MAPfile <<
" else if (par == \"FiberCh\") return (val == FiberCh);" << std::endl;
13629 MAPfile <<
" else if (par == \"QIE\") return (val == QIE);" << std::endl;
13630 MAPfile <<
" else if (par == \"ADC\") return (val == ADC);" << std::endl;
13631 MAPfile <<
" else if (par == \"VMECardID\") return (val == VMECardID);" << std::endl;
13632 MAPfile <<
" else if (par == \"dccID\") return (val == dccID);" << std::endl;
13633 MAPfile <<
" else if (par == \"Spigot\") return (val == Spigot);" << std::endl;
13634 MAPfile <<
" else if (par == \"FiberIndex\") return (val == FiberIndex);" << std::endl;
13635 MAPfile <<
" else if (par == \"HtrSlot\") return (val == HtrSlot);" << std::endl;
13636 MAPfile <<
" else if (par == \"HtrTB\") return (val == HtrTB);" << std::endl;
13637 MAPfile <<
" else return false;" << std::endl;
13638 MAPfile <<
" }" << std::endl;
13641 MAPfile <<
" bool check(const std::string par, const std::string val) const" << std::endl;
13642 MAPfile <<
" {" << std::endl;
13643 MAPfile <<
" if (par == \"subdet\") return (val == subdet);" << std::endl;
13644 MAPfile <<
" else if (par == \"RBX\") return (val == RBX);" << std::endl;
13645 MAPfile <<
" else return false;" << std::endl;
13646 MAPfile <<
" }" << std::endl;
13648 MAPfile <<
"};" << std::endl;
13651 MAPfile <<
"const Cell AllCells[] = {" << std::endl;
13652 MAPfile <<
"//{ SD, Eta, Phi, Depth, RBX, RM, PIXEL, RMfiber, Fiber Ch., QIE, ADC, VMECrateId, dccid, spigot, " 13653 "fiberIndex, htrSlot, htrTopBottom }" 13666 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13667 << detid->
depth() <<
" ,";
13668 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13669 << lmap_entry.
rmFiber() <<
" , ";
13673 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13675 MAPfile <<
"}," << std::endl;
13688 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() - 1 <<
" , " << detid->
iphi() - 1 <<
" , " 13689 << detid->
depth() <<
" ,";
13690 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13691 << lmap_entry.
rmFiber() <<
" , ";
13695 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13697 MAPfile <<
"}," << std::endl;
13704 for (
int eta = -20;
eta <= -20;
eta++) {
13712 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13713 << detid->
depth() <<
" ,";
13714 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13715 << lmap_entry.
rmFiber() <<
" , ";
13719 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13721 MAPfile <<
"}," << std::endl;
13727 for (
int eta = -19;
eta <= -16;
eta++) {
13735 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13736 << detid->
depth() <<
" ,";
13737 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13738 << lmap_entry.
rmFiber() <<
" , ";
13742 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13744 MAPfile <<
"}," << std::endl;
13749 for (
int eta = -29;
eta <= -16;
eta++) {
13757 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13758 << detid->
depth() <<
" ,";
13759 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13760 << lmap_entry.
rmFiber() <<
" , ";
13764 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13766 MAPfile <<
"}," << std::endl;
13779 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13780 << detid->
depth() <<
" ,";
13781 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13782 << lmap_entry.
rmFiber() <<
" , ";
13786 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13788 MAPfile <<
"}," << std::endl;
13796 for (
int eta = -41;
eta <= -29;
eta++) {
13804 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13805 << detid->
depth() <<
" ,";
13806 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13807 << lmap_entry.
rmFiber() <<
" , ";
13811 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13813 MAPfile <<
"}," << std::endl;
13827 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13828 << detid->
depth() <<
" ,";
13829 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13830 << lmap_entry.
rmFiber() <<
" , ";
13834 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13836 MAPfile <<
"}," << std::endl;
13852 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13853 << detid->
depth() <<
" ,";
13854 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13855 << lmap_entry.
rmFiber() <<
" , ";
13859 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13861 MAPfile <<
"}," << std::endl;
13875 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() - 1 <<
" , " << detid->
iphi() - 1 <<
" , " 13876 << detid->
depth() <<
" ,";
13877 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13878 << lmap_entry.
rmFiber() <<
" , ";
13882 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13884 MAPfile <<
"}," << std::endl;
13889 MAPfile <<
"};" << std::endl;
13892 MAPfile <<
"// macro for array length calculation" << std::endl;
13893 MAPfile <<
"#define DIM(a) (sizeof(a)/sizeof(a[0]))" << std::endl;
13896 MAPfile <<
"// class for cells array managing" << std::endl;
13897 MAPfile <<
"class CellDB {" << std::endl;
13898 MAPfile <<
"public:" << std::endl;
13899 MAPfile <<
" CellDB()" << std::endl;
13900 MAPfile <<
" : cells(AllCells, AllCells + DIM(AllCells))" << std::endl;
13901 MAPfile <<
"{}" << std::endl;
13904 MAPfile <<
"// return i-th cell" << std::endl;
13905 MAPfile <<
"Cell operator [] (int i) const {return cells[i];}" << std::endl;
13907 MAPfile <<
"// number of cells in database" << std::endl;
13908 MAPfile <<
"int size() const {return cells.size();}" << std::endl;
13911 MAPfile <<
"// select cells for which \"par\" == \"val\"" << std::endl;
13912 MAPfile <<
"template<typename T>" << std::endl;
13913 MAPfile <<
"CellDB find(const std::string par, const T val) const" << std::endl;
13915 MAPfile <<
" std::vector<Cell> s;" << std::endl;
13916 MAPfile <<
" for (size_t i = 0; i < cells.size(); ++i)" << std::endl;
13917 MAPfile <<
" if (cells[i].check(par, val))" << std::endl;
13918 MAPfile <<
" s.push_back(cells[i]);" << std::endl;
13919 MAPfile <<
" return CellDB(s);" << std::endl;
13920 MAPfile <<
"} " << std::endl;
13922 MAPfile <<
"private:" << std::endl;
13923 MAPfile <<
" CellDB( const std::vector<Cell> s)" << std::endl;
13924 MAPfile <<
" : cells(s)" << std::endl;
13925 MAPfile <<
"{}" << std::endl;
13926 MAPfile <<
"std::vector<Cell> cells;" << std::endl;
13927 MAPfile <<
"};" << std::endl;
13929 std::cout <<
"===== Finish writing Channel MAP =====" << std::endl;
13933 double deltaphi = phi1 - phi2;
13935 deltaphi = phi2 - phi1;
13937 if (deltaphi >
M_PI) {
13938 deltaphi = 2. *
M_PI - deltaphi;
13941 double tmp =
sqrt(deltaeta * deltaeta + deltaphi * deltaphi);
13951 if (a1 > 0.5 *
M_PI &&
a2 < 0.)
13953 if (
a2 > 0.5 *
M_PI && a1 < 0.)
13955 double tmp = (a1 * en1 +
a2 * en2) / (en1 + en2);
13968 double tmp =
a2 - a1;
13969 if (a1 *
a2 < 0.) {
13970 if (a1 > 0.5 *
M_PI)
TH2F * h_mapDepth4Pedw1_HO
TH2F * h_mapDepth4AmplE34_HO
double lsdep_estimator2_HEdepth3_
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_runbadrate0_depth4_HF
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]
TH1F * h_runbadrate0_depth3_HF
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_
TH1F * h_shape0_good_channels_HO
TH1F * h_bcnbadrate0_depth4_HO
constexpr DetId detid() const
Get the detector id.
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
T getParameter(std::string const &) const
TH1F * h_bcnvsamplitude0_HB
TH2F * h2_TSnVsAyear2023_HO
int flagToUseDigiCollectionsORNot_
TH1F * h_runnbadchannels_depth6_HE
TH2F * h_recSignalEnergy1_HE5
constexpr int32_t htrSlot() const
get the htr slot
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
const HcalPedestal * getPedestal(const HcalGenericDetId &fId) const
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 * h1_TSnVsAyear2023_HF
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
TH1F * h_nbadchannels_depth6_HE
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
T const & getData(const ESGetToken< T, R > &iToken) const noexcept(false)
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
TH1F * h_runbadrate_depth3_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
TH1F * h_nbadchannels_depth3_HB
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
TH1F * h_runnbadchannels_depth4_HB
TH2F * h_2DsumPedestalLS4
TH2F * h_recSignalEnergy1_HE2
TH2F * h_2D0sumADCAmplLS1
TH1F * h_A_Depth2RADDAM_HE
float getValue(int fCapId) const
get value for capId = 0..3
TH1F * h_sumPedestalperLS4
TH2F * h_mapADCCalib047_HE
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
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
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
TH2F * h2_TSnVsAyear2023_HE
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
const DetId lookup(HcalElectronicsId fId) const
lookup the logical detid associated with the given electronics id
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
TH1F * h_runbadrate_depth4_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_runnbadchannels_depth5_HE
TH1F * h_gsmdifferencefit3_HE
const HcalQIEShape * shape
TH1F * h_sumamplitude_depth2_HE
double lsdep_estimator1_HFdepth4_
TH2F * h_mapDepth2linADCAmpl12_HE
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
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_
constexpr int32_t htrTopBottom() const
get the htr top/bottom (1=top/0=bottom), valid for VME
TH1F * h_Amplitude_forCapIdErrors_HF1
TH2F * h_2DsumAmplitudeLS3
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]
const HcalElectronicsMap * getHcalMapping() const
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_
TH1F * h_runbadrate_depth7_HE
TH1F * h_runnbadchannelsC_depth7_HE
TH2F * h_mapDepth1ADCAmpl225_HF
TH1F * h_nbadchannels_depth4_HB
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
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
TH2F * h_mapDepth2Pedw1_HE
TH1F * h_sum0TSmaxAperLS4
TH1F * h_runbadrateC_depth7_HE
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
const HcalPedestalWidth * getPedestalWidth(const HcalGenericDetId &fId) const
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 * h1_TSnVsAyear20230_HB
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
constexpr DetId detid() const
Get the detector id.
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
constexpr float energy() const
TH1F * h_sumCutAmplperLS4
TH1F * h_sumCutTSmeanAperLS3
TH2F * h_mapDepth2ADCAmplSiPM_HE
TH2F * h_mapDepth3pedestal_HB
TH2F * h_amplitudechannel1_HB1
TH1F * h_nbadchannels_depth3_HF
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
TH1F * h_runnbadchannelsC_depth5_HE
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
T getUntrackedParameter(std::string const &, T const &) const
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_
TH2F * h_mapDepth1TSmeanA225_HF
TH1F * h1_TSnVsAyear20230_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
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
TH2F * h_mapDepth1Pedw3_HB
TH1F * h_pedestalwaver9_HO
TH2F * h_mapDepth1TSmaxA_HE
TH2F * h_mapDepth3ADCAmpl12_HE
constexpr HcalSubdetector subdet() const
get the subdetector
void adc2fC(const HBHEDataFrame &df, CaloSamples &lf) const override
TH2F * h_maprphinorm0_HE2
TH2F * h_mapDepth4AmplE34_HB
TH2F * h_mapDepth2Ped0_HE
constexpr int32_t readoutVMECrateId() const
get the readout VME crate number
TH2F * h_mapDepth1pedestal_HB
TH1F * h_Amplitude_forCapIdErrors_HE2
TH1F * h_bcnvsamplitude_HO
TH2F * h_mapDepth3TSmaxA225_HB
TH1F * h_sumamplitudechannel_HE
double averNOSIGNALoccupancy_HB
TH2F * h_mapDepth4TSmeanA_HF
TH2F * h_mapADCAmplfirstpeak_HE
TH1F * h_sum0TSmaxAperLS2
TH2F * h_mapDepth4ADCAmpl225Copy_HF
TH1F * h_sum0PedestalperLS4
TH2F * h_recNoiseEnergy0_HE5
TH1F * h_runnbadchannelsC_depth4_HF
TH2F * h_mapDepth7TSmeanA_HE
TH2F * h_mapDepth4Ampl_HO
TH1F * h_nbadchannels_depth2_HF
TH1F * h_sum0AmplitudeperLS5
TH2F * h_mapDepth3ADCAmpl225_HF
TH2F * h_amplitudechannel0_HF4
TH1F * h_runbadrate_depth5_HE
TH1F * h_orbitNumvsamplitude0_HE
TH1F * h_sigLayer2RADDAM5_HED2
TH2F * h_recSignalEnergy0_HB2
TH1F * h_sumADCAmplperLS4
TH1F * h_eventoccupancy_HB
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
TH1F * h_runnbadchannels_depth4_HE
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
TH1F * h_runbadrate_depth4_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
TH1F * h_runbadrate0_depth5_HE
TH2F * h_mapDepth3Amplitude_HB
TH1F * h_energyhitNoise_HB
TH1F * h_sumADCAmplperLS1_M1
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
const HcalRespCorr * getHcalRespCorr(const HcalGenericDetId &fId) const
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)
constexpr int32_t spigot() const
get the spigot (input number on DCC, AMC card number for uTCA)
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
constexpr int ieta() const
get the cell ieta
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 * h1_TSnVsAyear2023_HB
TH1F * h_sumTSmeanAperLS2
const edm::ESGetToken< HcalTopology, HcalRecNumberingRecord > tokTopo_
TH2F * h_mapDepth4Amplitude225_HO
TH1F * h_sumamplitude_depth1_HB0
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
TH1F * h_runbadrate0_depth4_HB
TH2F * h_mapDepth1Pedw1_HB
bool usePedestalSubtraction_
TH2F * h_mapDepth2Pedw2_HF
TH1F * h1_TSnVsAyear20230_HE
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]
double const adc2fC_QIE10[256]
TH2F * h_2D0sumAmplitudeLS8
TH1F * h_sumCutADCAmplperLS3
double averNOSIGNALoccupancy_HE
TH2F * h_recSignalEnergy1_HB3
TH2F * h_recSignalEnergy0_HB3
TH2F * h_mapDepth6ADCAmpl225_HE
bool studyCalibCellsHist_
TH1F * h_nbadchannels_depth5_HE
TH2F * h_recSignalEnergy1_HB1
TH2F * h_recNoiseEnergy1_HF1
Abs< T >::type abs(const T &t)
TH1F * h_runnbadchannelsC_depth4_HE
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
static const float adc2fC[128]
TH1F * h_sumADCAmplperLS1_P2
TH2F * h_mapDepth2Ampl_HB
TH1F * h_runbadrate0_depth2_HF
TH2F * h_mapDepth2Amplitude_HE
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
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
#define DEFINE_FWK_MODULE(type)
double lsdep_estimator1_HBdepth2_
TH2F * h_2D0sumADCAmplLSdepth3HBu
TH1F * h_runbadrate_depth1_HE
TH1F * h_sumADCAmplperLS6_M1
TH2F * h_recNoiseEnergy2_HF3
const HcalLUTCorr * getHcalLUTCorr(const HcalGenericDetId &fId) const
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
TH1F * h_runnbadchannels_depth3_HB
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
const_iterator begin() const
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
float getSigma(int fCapId1, int fCapId2) const
get correlation element for capId1/2 = 0..3
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
TH1F * h_runbadrateC_depth3_HF
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
TH1F * h_runbadrateC_depth5_HE
double recSignalEnergy0[nsub][ndepth][neta][nphi]
TH2F * h_recNoiseEnergy2_HE1
def gen(fragment, howMuch)
Event to runs.
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
TH1F * h_runnbadchannelsC_depth3_HB
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
TH1F * h_nbadchannels_depth4_HE
TH2F * h_mapDepth3ADCAmpl225_HE
TH1F * h_sumCutAmplperLS1
TH2F * h_mapDepth4ADCAmpl12_HO
TH1F * h_runbadrate_depth3_HB
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
TH1F * h_runnbadchannelsC_depth6_HE
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
const_iterator end() const
TH2F * h_mapDepth6pedestalw_HE
TH2F * h_mapDepth3TSmeanA_HB
TH1F * h_runbadrate0_depth4_HE
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
TH1F * h_runbadrate_depth4_HF
double sumEstimator2[nsub][ndepth][neta][nphi]
TH2F * h_mapTSmaxCalib047_HF
TH1F * h_runbadrate0_depth3_HB
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
TH1F * h_runbadrateC_depth4_HF
TH2F * h_mapDepth1ADCAmpl225Copy_HB
TH2F * h_mapDepth2AmplE34_HB
double lsdep_estimator3_HEdepth1_
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
constexpr int samples() const
total number of samples in the digi
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
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
TH1F * h_runbadrate_depth6_HE
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
TH1F * h_runbadrate0_depth7_HE
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
const HcalGain * getGain(const HcalGenericDetId &fId) const
TH2F * h_2D0sumADCAmplLSdepth4HBu
TH2F * h_maprphinorm2_HE4
TH2F * h_recNoiseEnergy1_HE2
TH1F * h_sumCutTSmaxAperLS8
TH2F * h_mapDepth4TSmeanA_HO
TH2F * h_mapDepth3Ped0_HE
TH1F * h_averSIGNALsumamplitude_HE
TH2F * h_mapRatioCalib047_HF
double signal3[ndepth][neta][nphi]
TH2F * h_mapDepth2Ped1_HF
TH1F * h_bcnbadrate0_depth1_HB
TH2F * h_2D0sumAmplitudeLS7
constexpr int32_t fiberIndex() const
get the fiber index. For VME 1-8 (which of eight fibers carried by a spigot), for uTCA fibers are zer...
TH1F * h_bcnnbadchannels_depth2_HB
TH2F * h_mapDepth1Pedw2_HB
TH1F * h_runbadrate0_depth6_HE
TH1F * h_sumCutADCAmplperLS1
TH1F * h_sumCutAmplperLS2
TH2F * h_mapDepth6TSmaxA225_HE
TH2F * h_mapDepth4pedestalw_HO
TH1F * h_Eta16Depth3RADDAM_HED2
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
TH2F * h2_TSnVsAyear2023_HF
TH1F * h_sumphiEta16Depth3RADDAM_HED2P
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_runnbadchannelsC_depth3_HF
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
TH2F * h2_TSnVsAyear2023_HB
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_runnbadchannels_depth3_HF
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
const HcalQIECoder * getHcalCoder(const HcalGenericDetId &fId) const
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
TH1F * h_runnbadchannelsC_depth4_HB
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
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
TH1F * h1_TSnVsAyear2023_HO
TH2F * h_mapDepth2ADCAmpl225Copy_HF
TH1F * h_gsmdifferencefit2_HE
T * make(const Args &...args) const
make new ROOT object
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
TH1F * h1_TSnVsAyear2023_HE
double averSIGNALsumamplitude_HB
TH2F * h_maprphinorm0_HE6
TH2F * h_amplitudechannel2_HE4
TH1F * h_sumCutTSmaxAperLS7
double averSIGNALsumamplitude_HF
TH2F * h_mapTSmeanCalib_HO
TH2F * h_mapDepth2AmplE34_HF
TH1F * h_runbadrateC_depth4_HB
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_nbadchannels_depth7_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
TH1F * h_amplitudeaveragedbydepthes_HE
TH1F * h_sumADCAmplperLS3_P2
TH1F * h1_TSnVsAyear20230_HO
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_runbadrateC_depth4_HE
TH1F * h_sumCutTSmaxAperLS5
TH2F * h_mapDepth3ADCAmpl12SiPM_HE
TH2F * h_mapTSmaxCalib047_HO
const HcalTimeCorr * getHcalTimeCorr(const HcalGenericDetId &fId) const
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
TH1F * h_runnbadchannels_depth4_HF
TH1F * h_nbadchannels_depth4_HF
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
constexpr int iphi() const
get the cell iphi
TH1F * h_gsmdifferencefit1_HE
TH2F * h_recNoiseEnergy0_HE1
TH2F * h_mapenophinorm2_HE2
TH2F * h_mapDepth1ADCAmplSiPM_HE
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
TH1F * h_runnbadchannels_depth7_HE
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
constexpr int samples() const
total number of samples in the digi
TH1F * h_Amplitude_notCapIdErrors_HF1
double averSIGNALoccupancy_HF
TH2F * h_mapDepth3AmplE34_HE
constexpr int32_t dccid() const
get the (Hcal local) DCC id for VME, crate number for uTCA
TH2F * h_recNoiseEnergy2_HE2
TH1F * h_runbadrateC_depth6_HE
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
TH1F * h_runbadrateC_depth3_HB
TH1F * h_sumCutADCAmplperLS7
TH2F * h_2D0sumADCAmplLS4
TH1F * h_sumCutAmplitudeperLS4
constexpr int depth() const
get the tower depth
bool studyCapIDErrorsHist_
bool studyPedestalCorrelations_