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);
2385 std::cout <<
"============== endJob ===================================" << std::endl;
2387 std::cout <<
" --------------------------------------- " << std::endl;
2390 std::cout <<
" --------------------------------------------- " << std::endl;
2396 std::cout <<
"===== possible max number of events in ntuple(each 50th recorded) = " <<
nevent50 << std::endl;
2400 std::cout <<
"===== full number of events*HBHEdigis = " <<
nnnnnn << std::endl;
2405 std::cout <<
"===== Start writing user histograms =====" << std::endl;
2640 for (
int k1 = 0; k1 <
ndepth; k1++) {
2641 for (
int k2 = 0; k2 <
neta; k2++) {
2646 for (
int k3 = 0; k3 <
nphi; k3++) {
2705 int outabortgap = 1;
2715 std::cout <<
" --------------------------------------- " << std::endl;
2719 std::cout <<
" --------------------------------------------- " << std::endl;
2725 std::cout <<
" ---------***********************------------- " << std::endl;
2775 unsigned long int pcountall1 = 0;
2776 unsigned long int pcountall3 = 0;
2777 unsigned long int pcountall6 = 0;
2778 unsigned long int pcountall8 = 0;
2783 unsigned long int mcountall1 = 0;
2784 unsigned long int mcountall3 = 0;
2785 unsigned long int mcountall6 = 0;
2786 unsigned long int mcountall8 = 0;
2795 int pnnmin1 = 999999999;
2796 int pnnmin3 = 999999999;
2797 int pnnmin6 = 999999999;
2798 int pnnmin8 = 999999999;
2803 int mnnmin1 = 999999999;
2804 int mnnmin3 = 999999999;
2805 int mnnmin6 = 999999999;
2806 int mnnmin8 = 999999999;
2808 for (
int k1 = 0; k1 <
ndepth; k1++) {
2809 for (
int k3 = 0; k3 <
nphi; k3++) {
2810 for (
int k2 = 0; k2 <
neta; k2++) {
2919 int kkkk2 = (k2 - 1) / 4;
2924 int kkkk3 = (k3) / 4 + 1;
2928 ietaphi = ((kkkk2)-1) *
znphi + (kkkk3);
2996 if (bbbc / bbb1 > 25.) {
3011 if (bbbc / bbb1 > 25.) {
3096 if (bbbc / bbb1 > 15. && k3 % 2 == 0) {
3110 if (bbbc / bbb1 > 15. && k3 % 2 == 0) {
3250 if (bbbc / bbb1 > 20.) {
3265 if (bbbc / bbb1 > 20.) {
3357 if (bbbc / bbb1 > 80.) {
3372 if (bbbc / bbb1 > 80.) {
3927 if (
k0 == 0 && k1 == 0) {
3928 if (pcountmin1 > 0) {
3929 if (pcountmin1 < pnnmin1)
3930 pnnmin1 = pcountmin1;
3934 if (mcountmin1 > 0) {
3935 if (mcountmin1 < mnnmin1)
3936 mnnmin1 = mcountmin1;
3942 if (
k0 == 1 && k1 == 0) {
3943 if (pcountmin3 > 0) {
3944 if (pcountmin3 < pnnmin3)
3945 pnnmin3 = pcountmin3;
3949 if (mcountmin3 > 0) {
3950 if (mcountmin3 < mnnmin3)
3951 mnnmin3 = mcountmin3;
3957 if (
k0 == 2 && k1 == 3) {
3958 if (pcountmin8 > 0) {
3959 if (pcountmin8 < pnnmin8)
3960 pnnmin8 = pcountmin8;
3964 if (mcountmin8 > 0) {
3965 if (mcountmin8 < mnnmin8)
3966 mnnmin8 = mcountmin8;
3972 if (
k0 == 3 && k1 == 0) {
3973 if (pcountmin6 > 0) {
3974 if (pcountmin6 < pnnmin6)
3975 pnnmin6 = pcountmin6;
3979 if (mcountmin6 > 0) {
3980 if (mcountmin6 < mnnmin6)
3981 mnnmin6 = mcountmin6;
3996 float patiooccupancy1 = 0.;
3997 if (pcountall1 != 0)
3998 patiooccupancy1 = (
float)pnnmin1 * mnnbins1 / pcountall1;
4000 float matiooccupancy1 = 0.;
4001 if (mcountall1 != 0)
4002 matiooccupancy1 = (
float)mnnmin1 * mnnbins1 / mcountall1;
4005 float patiooccupancy3 = 0.;
4006 if (pcountall3 != 0)
4007 patiooccupancy3 = (
float)pnnmin3 * mnnbins3 / pcountall3;
4009 float matiooccupancy3 = 0.;
4010 if (mcountall3 != 0)
4011 matiooccupancy3 = (
float)mnnmin3 * mnnbins3 / mcountall3;
4014 float patiooccupancy6 = 0.;
4015 if (pcountall6 != 0)
4016 patiooccupancy6 = (
float)pnnmin6 * mnnbins6 / pcountall6;
4018 float matiooccupancy6 = 0.;
4019 if (mcountall6 != 0)
4020 matiooccupancy6 = (
float)mnnmin6 * mnnbins6 / mcountall6;
4023 float patiooccupancy8 = 0.;
4024 if (pcountall8 != 0)
4025 patiooccupancy8 = (
float)pnnmin8 * mnnbins8 / pcountall8;
4027 float matiooccupancy8 = 0.;
4028 if (mcountall8 != 0)
4029 matiooccupancy8 = (
float)mnnmin8 * mnnbins8 / mcountall8;
4033 for (
int k1 = 0; k1 <
ndepth; k1++) {
4034 for (
int k2 = 0; k2 <
neta; k2++) {
4035 for (
int k3 = 0; k3 <
nphi; k3++) {
4135 for (
int k1 = 0; k1 <
ndepth; k1++) {
4136 for (
int k2 = 0; k2 <
neta; k2++) {
4137 for (
int k3 = 0; k3 <
nphi; k3++) {
4159 for (
int k1 = 0; k1 <
ndepth; k1++) {
4160 for (
int k2 = 0; k2 <
neta; k2++) {
4161 for (
int k3 = 0; k3 <
nphi; k3++) {
4193 bool gotHFDigis =
true;
4197 if (!(
hf.isValid())) {
4201 std::cout <<
" ****************************** =========================== No HFDigiCollection found " 4206 eta = digi->id().ieta();
4207 phi = digi->id().iphi();
4208 depth = digi->id().depth();
4225 for (
int i = 0;
i <
nTS;
i++) {
4244 bool gotQIE10Digis =
true;
4246 gotQIE10Digis =
false;
4249 gotQIE10Digis =
false;
4251 if (!gotQIE10Digis) {
4252 std::cout <<
" No QIE10DigiCollection collection is found " << std::endl;
4255 double totalAmplitudeHF = 0.;
4256 for (
unsigned int j = 0;
j < qie10dc.
size();
j++) {
4280 double amplitudefullTSs = 0.;
4281 double nnnnnnTS = 0.;
4282 for (
int i = 0;
i <
nTS; ++
i) {
4285 int adc = qie10df[
i].adc();
4317 bool gotHBHEDigis =
true;
4319 gotHBHEDigis =
false;
4320 if (!(
hbhe.isValid()))
4321 gotHBHEDigis =
false;
4322 if (!gotHBHEDigis) {
4323 std::cout <<
" No HBHEDigiCollection collection is found " << std::endl;
4326 double totalAmplitudeHB = 0.;
4327 double totalAmplitudeHE = 0.;
4328 double nnnnnnTSHB = 0.;
4329 double nnnnnnTSHE = 0.;
4332 eta = digi->id().ieta();
4333 phi = digi->id().iphi();
4334 depth = digi->id().depth();
4374 double amplitudefullTSs = 0.;
4377 for (
int i = 0;
i <
nTS;
i++) {
4389 double amplitudefullTSs = 0.;
4392 for (
int i = 0;
i <
nTS;
i++) {
4407 if (totalAmplitudeHB != 0.) {
4412 if (totalAmplitudeHE != 0.) {
4428 bool gotQIE11Digis =
true;
4430 gotQIE11Digis =
false;
4432 gotQIE11Digis =
false;
4433 if (!gotQIE11Digis) {
4434 std::cout <<
" No QIE11DigiCollection collection is found " << std::endl;
4437 double totalAmplitudeHBQIE11 = 0.;
4438 double totalAmplitudeHEQIE11 = 0.;
4439 double nnnnnnTSHBQIE11 = 0.;
4440 double nnnnnnTSHEQIE11 = 0.;
4441 for (
unsigned int j = 0;
j < qie11dc.
size();
j++) {
4450 int sub = hcaldetid.
subdet();
4463 if (sub == 1 &&
depth == 1 && qwert1 == 0) {
4467 if (sub == 1 &&
depth == 2 && qwert2 == 0) {
4471 if (sub == 2 &&
depth == 1 && qwert3 == 0) {
4475 if (sub == 2 &&
depth == 2 && qwert4 == 0) {
4479 if (sub == 2 &&
depth == 3 && qwert5 == 0) {
4491 double amplitudefullTSs1 = 0.;
4492 double amplitudefullTSs6 = 0.;
4494 for (
int i = 0;
i <
nTS; ++
i) {
4495 int adc = qie11df[
i].adc();
4498 amplitudefullTSs1 += charge1;
4499 amplitudefullTSs6 += charge6;
4505 totalAmplitudeHBQIE11 +=
charge;
4512 double amplitudefullTSs1 = 0.;
4513 double amplitudefullTSs6 = 0.;
4515 for (
int i = 0;
i <
nTS;
i++) {
4516 int adc = qie11df[
i].adc();
4519 amplitudefullTSs1 += charge1;
4520 amplitudefullTSs6 += charge6;
4526 totalAmplitudeHEQIE11 +=
charge;
4535 if (totalAmplitudeHBQIE11 != 0.) {
4540 if (totalAmplitudeHEQIE11 != 0.) {
4551 bool gotHODigis =
true;
4554 if (!(
ho.isValid()))
4558 std::cout <<
" No HO collection is found " << std::endl;
4561 double totalAmplitudeHO = 0.;
4563 eta = digi->id().ieta();
4564 phi = digi->id().iphi();
4565 depth = digi->id().depth();
4586 double nnnnnnTS = 0.;
4587 double amplitudefullTSs = 0.;
4589 for (
int i = 0;
i <
nTS;
i++) {
4618 bool gotHBHERecHitsNoise =
true;
4620 gotHBHERecHitsNoise =
4623 gotHBHERecHitsNoise =
false;
4624 if (!gotHBHERecHitsNoise) {
4626 std::cout <<
" No RecHits HBHENoise collection is found " << std::endl;
4631 float icalconst = 1.;
4634 HBHERecHit aHit(hbheItr->id(), hbheItr->eraw() * icalconst, hbheItr->time());
4636 double energyhit = aHit.
energy();
4637 DetId id = (*hbheItr).detid();
4639 int sub = ((hid).rawId() >> 25) & 0x7;
4648 if (hid.
depth() == 1 && sub == 1 && hid.
iphi() == 25) {
4650 std::cout <<
" Noise,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4651 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4653 if (hid.
depth() == 1 && sub == 2 && hid.
iphi() == 25) {
4655 std::cout <<
" Noise,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4656 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4662 int mdepth = hid.
depth();
4673 bool gotHFRecHitsNoise =
true;
4675 gotHFRecHitsNoise =
false;
4677 gotHFRecHitsNoise =
false;
4678 if (!gotHFRecHitsNoise) {
4680 std::cout <<
" No RecHits HFNoise collection is found " << std::endl;
4684 float icalconst = 1.;
4687 HFRecHit aHit(hfItr->id(), hfItr->energy() * icalconst, hfItr->time());
4688 double energyhit = aHit.
energy();
4689 DetId id = (*hfItr).detid();
4691 int sub = ((hid).rawId() >> 25) & 0x7;
4697 if (hid.
iphi() == 25) {
4699 std::cout <<
"HF Noise,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4700 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4706 int mdepth = hid.
depth();
4720 bool gotHBHERecHitsSignal =
true;
4722 gotHBHERecHitsSignal =
4725 gotHBHERecHitsSignal =
false;
4726 if (!gotHBHERecHitsSignal) {
4728 std::cout <<
" No RecHits HBHESignal collection is found " << std::endl;
4733 float icalconst = 1.;
4736 HBHERecHit aHit(hbheItr->id(), hbheItr->eraw() * icalconst, hbheItr->time());
4738 double energyhit = aHit.
energy();
4739 DetId id = (*hbheItr).detid();
4741 int sub = ((hid).rawId() >> 25) & 0x7;
4749 if (hid.
depth() == 1 && sub == 1 && hid.
iphi() == 25) {
4751 std::cout <<
"HBHE Signal,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4752 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4754 if (hid.
depth() == 1 && sub == 2 && hid.
iphi() == 25) {
4756 std::cout <<
"HBHE Signal,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4757 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4763 int mdepth = hid.
depth();
4774 bool gotHFRecHitsSignal =
true;
4776 gotHFRecHitsSignal =
false;
4778 gotHFRecHitsSignal =
false;
4779 if (!gotHFRecHitsSignal) {
4781 std::cout <<
" No RecHits HFSignal collection is found " << std::endl;
4785 float icalconst = 1.;
4788 HFRecHit aHit(hfItr->id(), hfItr->energy() * icalconst, hfItr->time());
4789 double energyhit = aHit.
energy();
4790 DetId id = (*hfItr).detid();
4792 int sub = ((hid).rawId() >> 25) & 0x7;
4798 if (hid.
iphi() == 25) {
4800 std::cout <<
"HF Signal,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4801 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4807 int mdepth = hid.
depth();
4837 for (
int k1 = 0; k1 <
ndepth; k1++) {
4839 for (
int k2 = 0; k2 <
neta; k2++) {
4841 int k2plot = k2 - 41;
4844 double sumoverphi = 0;
4845 int nsumoverphi = 0;
4846 for (
int k3 = 0; k3 <
nphi; k3++) {
4851 std::cout <<
"==== nsumoverphi = " << nsumoverphi <<
" sumoverphi = " << sumoverphi
4852 <<
" k1 = " << k1 <<
" k2 = " << k2 <<
" kkk = " << kkk <<
" k3 = " << k3 << std::endl;
4856 for (
int k3 = 0; k3 <
nphi; k3++) {
4857 if (nsumoverphi != 0) {
4860 std::cout <<
"nsumoverphi= " << nsumoverphi <<
" sumoverphi= " << sumoverphi <<
" k1= " << k1
4861 <<
" k2= " << k2 <<
" kkk= " << kkk <<
" k3= " << k3
4862 <<
" maprphinorm= " <<
maprphinorm[
k0][k1][k2][k3] << std::endl;
4937 for (
int k1 = 0; k1 <
ndepth; k1++) {
4939 for (
int k2 = 0; k2 <
neta; k2++) {
4941 int k2plot = k2 - 41;
4943 for (
int k3 = 0; k3 <
nphi; k3++) {
4971 for (
int k1 = 0; k1 <
ndepth; k1++) {
4973 for (
int k2 = 0; k2 <
neta; k2++) {
4975 int k2plot = k2 - 41;
4977 for (
int k3 = 0; k3 <
nphi; k3++) {
5020 for (
int k1 = 0; k1 <
ndepth; k1++) {
5022 for (
int k2 = 0; k2 <
neta; k2++) {
5024 int k2plot = k2 - 41;
5026 for (
int k3 = 0; k3 <
nphi; k3++) {
5069 for (
int k1 = 0; k1 <
ndepth; k1++) {
5071 for (
int k2 = 0; k2 <
neta; k2++) {
5073 int k2plot = k2 - 41;
5075 for (
int k3 = 0; k3 <
nphi; k3++) {
5115 for (
int k1 = 0; k1 <
ndepth; k1++) {
5117 for (
int k2 = 0; k2 <
neta; k2++) {
5119 int k2plot = k2 - 41;
5121 for (
int k3 = 0; k3 <
nphi; k3++) {
5185 for (
int k1 = 0; k1 <
ndepth; k1++) {
5187 for (
int k2 = 0; k2 <
neta; k2++) {
5189 int k2plot = k2 - 41;
5191 for (
int k3 = 0; k3 <
nphi; k3++) {
5239 for (
int k1 = 0; k1 <
ndepth; k1++) {
5240 for (
int k2 = 0; k2 <
neta; k2++) {
5241 for (
int k3 = 0; k3 <
nphi; k3++) {
5253 for (
int k1 = 0; k1 <
ndepth; k1++) {
5254 for (
int k2 = 0; k2 <
neta; k2++) {
5257 if (k1 == 2 && k2 == 25) {
5261 }
else if (k1 == 2 && k2 == 56) {
5276 for (
int k1 = 0; k1 <
ndepth; k1++) {
5277 for (
int k2 = 0; k2 <
neta; k2++) {
5278 if (k1 == 2 && (k2 == 25 || k2 == 56)) {
5281 int k2plot = k2 - 41;
5291 int LLLLLL111111 = 0;
5292 if ((k1 == 0 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 30))
5295 int LLLLLL222222 = 0;
5296 if ((k1 == 1 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 27) ||
5297 (k1 == 2 && fabs(kkk - 41) > 26 && fabs(kkk - 41) < 29))
5300 if (LLLLLL111111 == 1) {
5304 if (LLLLLL222222 == 1) {
5317 for (
int k1 = 0; k1 <
ndepth; k1++) {
5318 for (
int k2 = 0; k2 <
neta; k2++) {
5319 if (k1 == 2 && (k2 == 25 || k2 == 56)) {
5321 int k2plot = k2 - 41;
5328 for (
int k3 = 0; k3 <
nphi; k3++) {
5331 int LLLLLL111111 = 0;
5332 if ((k1 == 0 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 30))
5334 int LLLLLL222222 = 0;
5335 if ((k1 == 1 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 27) ||
5336 (k1 == 2 && fabs(kkk - 41) > 26 && fabs(kkk - 41) < 29))
5338 if (LLLLLL111111 == 1) {
5342 if (LLLLLL222222 == 1) {
5353 for (
int k1 = 0; k1 <
ndepth; k1++) {
5354 for (
int k2 = 0; k2 <
neta; k2++) {
5357 for (
int k3 = 0; k3 <
nphi; k3++) {
5375 int sumofchannels = 0;
5376 double sumamplitudesubdet = 0.;
5377 int sumofchannels0 = 0;
5378 double sumamplitudesubdet0 = 0.;
5379 for (
int k2 = 0; k2 <
neta; k2++) {
5380 for (
int k3 = 0; k3 <
nphi; k3++) {
5392 if (sumamplitudechannel_HB > 80.) {
5393 sumamplitudesubdet += sumamplitudechannel_HB;
5396 if (sumamplitudechannel_HB > 0.) {
5397 sumamplitudesubdet0 += sumamplitudechannel_HB;
5405 double sumamplitudechannel_HE =
5408 if (sumamplitudechannel_HE > 200.) {
5409 sumamplitudesubdet += sumamplitudechannel_HE;
5412 if (sumamplitudechannel_HE > 0.) {
5413 sumamplitudesubdet0 += sumamplitudechannel_HE;
5423 if (sumamplitudechannel_HO > 1200.) {
5424 sumamplitudesubdet += sumamplitudechannel_HO;
5427 if (sumamplitudechannel_HO > 0.) {
5428 sumamplitudesubdet0 += sumamplitudechannel_HO;
5437 if (sumamplitudechannel_HF > 600.) {
5438 sumamplitudesubdet += sumamplitudechannel_HF;
5441 if (sumamplitudechannel_HF > 0.) {
5442 sumamplitudesubdet0 += sumamplitudechannel_HF;
5455 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM1)
5456 maxxSUM1 = sumamplitudesubdet + sumamplitudesubdet0;
5457 if ((sumofchannels + sumofchannels0) >
maxxOCCUP1)
5463 if ((sumamplitudesubdet + sumamplitudesubdet0) > 60000) {
5464 for (
int k2 = 0; k2 <
neta; k2++) {
5465 for (
int k3 = 0; k3 <
nphi; k3++) {
5484 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM2)
5485 maxxSUM2 = sumamplitudesubdet + sumamplitudesubdet0;
5486 if ((sumofchannels + sumofchannels0) >
maxxOCCUP2)
5496 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM3)
5497 maxxSUM3 = sumamplitudesubdet + sumamplitudesubdet0;
5498 if ((sumofchannels + sumofchannels0) >
maxxOCCUP3)
5508 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM4)
5509 maxxSUM4 = sumamplitudesubdet + sumamplitudesubdet0;
5510 if ((sumofchannels + sumofchannels0) >
maxxOCCUP4)
5531 for (
int k1 = 0; k1 <
ndepth; k1++) {
5533 int nbadchannels = 0;
5534 for (
int k2 = 0; k2 <
neta; k2++) {
5535 for (
int k3 = 0; k3 <
nphi; k3++) {
5548 if (nbadchannels != 0)
5561 if (nbadchannels != 0)
5577 if (nbadchannels != 0)
5590 if (nbadchannels != 0)
5603 if (nbadchannels != 0)
5619 if (nbadchannels != 0)
5635 if (nbadchannels != 0)
5648 if (nbadchannels != 0)
5668 bool gotCALIBDigis =
true;
5670 gotCALIBDigis =
false;
5672 if (!(
calib.isValid())) {
5673 gotCALIBDigis =
false;
5675 if (!gotCALIBDigis) {
5678 int cal_det = digi->id().hcalSubdet();
5679 int cal_phi = digi->id().iphi();
5680 int cal_eta = digi->id().ieta();
5681 int cal_cbox = digi->id().cboxChannel();
5685 if (cal_det > 0 && cal_det < 5 && cal_cbox == 0) {
5686 int iphi = cal_phi - 1;
5691 double max_signal = -100.;
5692 int ts_with_max_signal = -100;
5697 for (
int i = 0;
i <
nTS;
i++) {
5698 double ampldefault =
adc2fC[digi->sample(
i).adc() & 0xff];
5699 if (max_signal < ampldefault) {
5700 max_signal = ampldefault;
5701 ts_with_max_signal =
i;
5706 timew += (
i + 1) * ampldefault;
5711 double aveamplitude = -100.;
5714 double aveamplitude1 = aveamplitude - 1;
5719 double ampldefault =
adc2fC[digi->sample(
ii).adc() & 0xff];
5720 double aaaaaa = (
ii + 1) - aveamplitude;
5721 double aaaaaa2 = aaaaaa * aaaaaa;
5722 rmsamp += (aaaaaa2 * ampldefault);
5724 double rmsamplitude = -100.;
5732 if (ts_with_max_signal > -1 && ts_with_max_signal <
nTS)
5734 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 <
nTS)
5735 calib2[cal_det - 1][
ieta + 41][
iphi] +=
adc2fC[digi->sample(ts_with_max_signal + 1).adc() & 0xff];
5736 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 <
nTS)
5737 calib2[cal_det - 1][
ieta + 41][
iphi] +=
adc2fC[digi->sample(ts_with_max_signal + 2).adc() & 0xff];
5738 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 <
nTS)
5739 calib2[cal_det - 1][
ieta + 41][
iphi] +=
adc2fC[digi->sample(ts_with_max_signal - 1).adc() & 0xff];
5740 if (ts_with_max_signal - 2 > -1 && ts_with_max_signal - 2 <
nTS)
5741 calib2[cal_det - 1][
ieta + 41][
iphi] +=
adc2fC[digi->sample(ts_with_max_signal - 2).adc() & 0xff];
5743 bool anycapid =
true;
5746 int error1 = 0, error2 = 0, error3 = 0;
5747 int lastcapid = 0, capid = 0;
5748 for (
int ii = 0;
ii < (*digi).size();
ii++) {
5749 capid = (*digi)[
ii].capid();
5750 bool er = (*digi)[
ii].er();
5751 bool dv = (*digi)[
ii].dv();
5752 if (
ii != 0 && ((lastcapid + 1) % 4) != capid)
5766 if (error1 != 0 || error2 != 0 || error3 != 0)
5781 for (
int k1 = 0; k1 <
nsub; k1++) {
5782 for (
int k2 = 0; k2 <
neta; k2++) {
5783 for (
int k3 = 0; k3 <
nphi; k3++) {
5784 int k2plot = k2 - 41;
5788 if (
signal[k1][k2][k3] > 0.) {
5814 if (k20 > 0 || k20 == 0)
5819 int kk2 = 0, kk3 = 0;
5820 if (k1 == 0 || k1 == 1) {
5825 if (k30 == 71 || k30 ==
nphi || k30 == 1 || k30 == 2)
5828 kk3 = ((k30 - 3) / 4) * 4 + 3;
5829 }
else if (k1 == 2) {
5830 if (
abs(k20) <= 4) {
5832 if (k30 == 71 || k30 ==
nphi || k30 == 1 || k30 == 2 || k30 == 3 || k30 == 4)
5835 kk3 = ((k30 - 5) / 6) * 6 + 5;
5837 if (
abs(k20) > 4 &&
abs(k20) <= 10)
5839 if (
abs(k20) > 10 &&
abs(k20) <= 15)
5843 if (k30 == 71 || k30 ==
nphi || (k30 >= 1 && k30 <= 10))
5846 kk3 = ((k30 - 11) / 12) * 12 + 11;
5848 }
else if (k1 == 3) {
5853 if (k30 >= 1 && k30 <= 18)
5855 if (k30 >= 19 && k30 <= 36)
5857 if (k30 >= 37 && k30 <= 54)
5859 if (k30 >= 55 && k30 <=
nphi)
5863 int kkk2 = kk2 + 41;
5870 double GetRMSOverNormalizedSignal = -1.;
5871 if (
signal[k1][k2][k3] > 0. &&
calib0[k1][kkk2][kkk3] > 0.) {
5872 GetRMSOverNormalizedSignal =
signal[k1][k2][k3] /
calib0[k1][kkk2][kkk3];
5892 if (
signal[k1][k2][k3] > 0.) {
5895 if (
calib0[k1][kkk2][kkk3] > 0.)
5899 if (
calib0[k1][kkk2][kkk3] > 0.)
5902 float calibtsmax =
calibt[k1][kkk2][kkk3];
5904 float calibtsmean =
caliba[k1][kkk2][kkk3];
5906 float calibwidth =
calibw[k1][kkk2][kkk3];
5908 float calibcap = -100.;
5925 if (calibtsmax > -0.5) {
5932 if (calibtsmean > -0.5) {
5939 if (calibwidth > -0.5) {
5949 h_map_HB->Fill(
double(k2plot),
double(k3), 1.);
5965 if (calibtsmax > -0.5) {
5972 if (calibtsmean > -0.5) {
5979 if (calibwidth > -0.5) {
5989 h_map_HE->Fill(
double(k2plot),
double(k3), 1.);
6005 if (calibtsmax > -0.5) {
6012 if (calibtsmean > -0.5) {
6019 if (calibwidth > -0.5) {
6029 h_map_HO->Fill(
double(k2plot),
double(k3), 1.);
6045 if (calibtsmax > -0.5) {
6052 if (calibtsmean > -0.5) {
6059 if (calibwidth > -0.5) {
6069 h_map_HF->Fill(
double(k2plot),
double(k3), 1.);
6101 std::cout <<
"======================== beignJob START +++++++++++++++++++++++++++" << std::endl;
6339 fs_->
make<TH1F>(
"h_ndepthesperamplitudebins_HE",
" ", 10, 0., 10.);
6349 "h_mapADCAmplsecondpeak0_HE",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
6706 fs_->
make<TH1F>(
"h_AamplitudewithPedSubtr_RADDAM_HEzoom0",
" ", 100, min80, 4000.);
6708 fs_->
make<TH1F>(
"h_AamplitudewithPedSubtr_RADDAM_HEzoom1",
" ", 100, min80, 1000.);
6722 fs_->
make<TH1F>(
"h_sumphiEta16Depth3RADDAM_HED2ALL",
" ", 100, min90, 70. * max90);
6870 float zaR = 10000000.;
6871 float zaR2 = 50000000.;
6905 fs_->
make<TH1F>(
"h_sumADCAmplEtaPhiLs_lscounterM1orbitNum",
" ", 600, 1., 601.);
6913 float bac2 = bac + 1.;
7887 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HB",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7889 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HB",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7891 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HE",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7893 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HE",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7895 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HF",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7897 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HF",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7899 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HO",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7901 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HO",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7950 float spls0 = 10000;
7960 float spls1 = 100000;
7996 std::cout <<
"======================== booking DONE +++++++++++++++++++++++++++" << std::endl;
8005 std::cout <<
"======================== beignJob finish +++++++++++++++++++++++++++" << std::endl;
8018 int mdepth = cell.
depth();
8019 int iphi = cell.iphi() - 1;
8020 int ieta = cell.ieta();
8023 int sub = cell.subdet();
8025 int errorGeneral = 0;
8034 bool anycapid =
true;
8039 int lastcapid = 0, capid = 0;
8040 int ERRORfiber = -10;
8041 int ERRORfiberChan = -10;
8042 int ERRORfiberAndChan = -10;
8043 int repetedcapid = 0;
8045 TSsize = digiItr->size();
8048 for (
int ii = 0;
ii < TSsize;
ii++) {
8049 capid = (*digiItr)[
ii].capid();
8050 bool er = (*digiItr)[
ii].er();
8051 bool dv = (*digiItr)[
ii].dv();
8052 int fiber = (*digiItr)[
ii].fiber();
8053 int fiberChan = (*digiItr)[
ii].fiberChan();
8054 int fiberAndChan = (*digiItr)[
ii].fiberAndChan();
8055 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8058 ERRORfiberChan = fiberChan;
8059 ERRORfiberAndChan = fiberAndChan;
8060 if (capid != lastcapid) {
8073 ERRORfiberChan = fiberChan;
8074 ERRORfiberAndChan = fiberAndChan;
8079 ERRORfiberChan = fiberChan;
8080 ERRORfiberAndChan = fiberAndChan;
8086 if (firstcapid == 0 && !anycapid)
8088 if (firstcapid == 1 && !anycapid)
8090 if (firstcapid == 2 && !anycapid)
8092 if (firstcapid == 3 && !anycapid)
8101 if (!anycapid && anyer)
8103 if (!anycapid && !anydv)
8105 if (!anycapid && anyer && !anydv)
8107 if (anyer && !anydv)
8113 for (
int ii = 0;
ii < TSsize;
ii++) {
8114 double ampldefault =
adc2fC[digiItr->sample(
ii).adc()];
8115 ampl += ampldefault;
8132 if (error1 != 0 || error2 != 0 || error3 != 0) {
8159 if (error1 != 0 || error2 != 0 || error3 != 0) {
8190 int mdepth = hcaldetid.
depth();
8191 int sub = hcaldetid.
subdet();
8195 bool anycapid =
true;
8200 int lastcapid = 0, capid = 0;
8201 int repetedcapid = 0;
8206 capid = qie11df[
ii].capid();
8207 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8209 if (capid != lastcapid) {
8224 int adc = qie11df[
ii].adc();
8230 ampl += ampldefault;
8289 int mdepth = cell.
depth();
8290 int iphi = cell.iphi() - 1;
8291 int ieta = cell.ieta();
8294 int sub = cell.subdet();
8296 std::cout <<
" HF DIGI ?????????????? ERROR mdepth = " << mdepth << std::endl;
8298 int errorGeneral = 0;
8307 bool anycapid =
true;
8312 int lastcapid = 0, capid = 0;
8313 int ERRORfiber = -10;
8314 int ERRORfiberChan = -10;
8315 int ERRORfiberAndChan = -10;
8316 int repetedcapid = 0;
8319 TSsize = digiItr->size();
8321 for (
int ii = 0;
ii < TSsize;
ii++) {
8322 capid = (*digiItr)[
ii].capid();
8323 bool er = (*digiItr)[
ii].er();
8324 bool dv = (*digiItr)[
ii].dv();
8325 int fiber = (*digiItr)[
ii].fiber();
8326 int fiberChan = (*digiItr)[
ii].fiberChan();
8327 int fiberAndChan = (*digiItr)[
ii].fiberAndChan();
8328 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8331 ERRORfiberChan = fiberChan;
8332 ERRORfiberAndChan = fiberAndChan;
8333 if (capid != lastcapid) {
8344 ERRORfiberChan = fiberChan;
8345 ERRORfiberAndChan = fiberAndChan;
8350 ERRORfiberChan = fiberChan;
8351 ERRORfiberAndChan = fiberAndChan;
8355 if (firstcapid == 0 && !anycapid)
8357 if (firstcapid == 1 && !anycapid)
8359 if (firstcapid == 2 && !anycapid)
8361 if (firstcapid == 3 && !anycapid)
8369 if (!anycapid && anyer)
8371 if (!anycapid && !anydv)
8373 if (!anycapid && anyer && !anydv)
8375 if (anyer && !anydv)
8379 for (
int ii = 0;
ii < TSsize;
ii++) {
8380 double ampldefault =
adc2fC[digiItr->sample(
ii).adc()];
8381 ampl += ampldefault;
8395 if (error1 != 0 || error2 != 0 || error3 != 0) {
8425 int mdepth = hcaldetid.
depth();
8426 int sub = hcaldetid.
subdet();
8430 bool anycapid =
true;
8435 int lastcapid = 0, capid = 0;
8436 int repetedcapid = 0;
8441 capid = qie10df[
ii].capid();
8442 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8444 if (capid != lastcapid) {
8459 int adc = qie10df[
ii].adc();
8461 ampl += ampldefault;
8492 int mdepth = cell.
depth();
8493 int iphi = cell.iphi() - 1;
8494 int ieta = cell.ieta();
8497 int sub = cell.subdet();
8498 int errorGeneral = 0;
8507 bool anycapid =
true;
8512 int lastcapid = 0, capid = 0;
8513 int ERRORfiber = -10;
8514 int ERRORfiberChan = -10;
8515 int ERRORfiberAndChan = -10;
8516 int repetedcapid = 0;
8517 for (
int ii = 0;
ii < (*digiItr).size();
ii++) {
8518 capid = (*digiItr)[
ii].capid();
8519 bool er = (*digiItr)[
ii].er();
8520 bool dv = (*digiItr)[
ii].dv();
8521 int fiber = (*digiItr)[
ii].fiber();
8522 int fiberChan = (*digiItr)[
ii].fiberChan();
8523 int fiberAndChan = (*digiItr)[
ii].fiberAndChan();
8524 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8527 ERRORfiberChan = fiberChan;
8528 ERRORfiberAndChan = fiberAndChan;
8529 if (capid != lastcapid) {
8542 ERRORfiberChan = fiberChan;
8543 ERRORfiberAndChan = fiberAndChan;
8548 ERRORfiberChan = fiberChan;
8549 ERRORfiberAndChan = fiberAndChan;
8555 if (firstcapid == 0 && !anycapid)
8557 if (firstcapid == 1 && !anycapid)
8559 if (firstcapid == 2 && !anycapid)
8561 if (firstcapid == 3 && !anycapid)
8570 if (!anycapid && anyer)
8572 if (!anycapid && !anydv)
8574 if (!anycapid && anyer && !anydv)
8576 if (anyer && !anydv)
8581 for (
int ii = 0;
ii < (*digiItr).size();
ii++) {
8582 double ampldefault =
adc2fC[digiItr->sample(
ii).adc()];
8583 ampl += ampldefault;
8599 if (error1 != 0 || error2 != 0 || error3 != 0) {
8608 std::cout <<
" mdepth HO = " << mdepth << std::endl;
8621 double toolwithPedSubtr[100];
8622 double lintoolwithoutPedSubtr[100];
8624 int mdepth = cell.
depth();
8625 int iphi = cell.iphi() - 1;
8626 int ieta0 = cell.ieta();
8630 int sub = cell.subdet();
8642 coder.
adc2fC(*digiItr, toolOriginal);
8643 double pedestalaver9 = 0.;
8644 double pedestalaver4 = 0.;
8645 double pedestal0 = 0.;
8646 double pedestal1 = 0.;
8647 double pedestal2 = 0.;
8648 double pedestal3 = 0.;
8649 double pedestalwaver9 = 0.;
8650 double pedestalwaver4 = 0.;
8651 double pedestalw0 = 0.;
8652 double pedestalw1 = 0.;
8653 double pedestalw2 = 0.;
8654 double pedestalw3 = 0.;
8655 double difpedestal0 = 0.;
8656 double difpedestal1 = 0.;
8657 double difpedestal2 = 0.;
8658 double difpedestal3 = 0.;
8661 double absamplitude = 0.;
8662 double amplitude345 = 0.;
8664 double linamplitudewithoutPedSubtr = 0.;
8667 double max_signal = -100.;
8668 int ts_with_max_signal = -100;
8674 double errorBtype = 0.;
8677 TSsize = digiItr->size();
8678 if ((*digiItr).size() != 10)
8681 for (
int ii = 0;
ii < TSsize;
ii++) {
8683 double ampldefaultwithPedSubtr = 0.;
8684 double linampldefaultwithoutPedSubtr = 0.;
8685 double ampldefault = 0.;
8686 double ampldefault0 = 0.;
8687 double ampldefault1 = 0.;
8688 double ampldefault2 = 0.;
8689 ampldefault0 =
adc2fC[digiItr->sample(
ii).adc()];
8691 ampldefault1 = toolOriginal[
ii];
8692 ampldefault2 = (*digiItr)[
ii].adc();
8694 ampldefault = ampldefault0;
8697 ampldefault = ampldefault1;
8700 ampldefault = ampldefault2;
8702 ampldefaultwithPedSubtr = ampldefault0;
8703 linampldefaultwithoutPedSubtr = ampldefault2;
8705 int capid = ((*digiItr)[
ii]).capid();
8707 double pedestalINI = pedestal00->
getValue(capid);
8709 double pedestalw = pedw->
getSigma(capid, capid);
8710 ampldefaultwithPedSubtr -=
pedestal;
8714 tool[
ii] = ampldefault;
8715 toolwithPedSubtr[
ii] = ampldefaultwithPedSubtr;
8716 lintoolwithoutPedSubtr[
ii] = linampldefaultwithoutPedSubtr;
8719 pedestalwaver9 += pedestalw * pedestalw;
8721 if (capid == 0 &&
c0 == 0) {
8726 pedestalw0 = pedestalw;
8727 pedestalwaver4 += pedestalw * pedestalw;
8728 difpedestal0 =
pedestal - pedestalINI;
8731 if (capid == 1 &&
c1 == 0) {
8736 pedestalw1 = pedestalw;
8737 pedestalwaver4 += pedestalw * pedestalw;
8738 difpedestal1 =
pedestal - pedestalINI;
8740 if (capid == 2 && c2 == 0) {
8745 pedestalw2 = pedestalw;
8746 pedestalwaver4 += pedestalw * pedestalw;
8747 difpedestal2 =
pedestal - pedestalINI;
8749 if (capid == 3 && c3 == 0) {
8754 pedestalw3 = pedestalw;
8755 pedestalwaver4 += pedestalw * pedestalw;
8756 difpedestal3 =
pedestal - pedestalINI;
8759 if (max_signal < ampldefault) {
8760 max_signal = ampldefault;
8761 ts_with_max_signal =
ii;
8776 absamplitude +=
abs(ampldefault);
8778 if (
ii == 3 ||
ii == 4 ||
ii == 5)
8779 amplitude345 += ampldefault;
8783 timew += (
ii + 1) *
abs(ampldefault);
8784 timeww += (
ii + 1) * ampldefault;
8788 if (sub == 1 || sub == 2) {
8793 pedestalaver9 /= TSsize;
8794 pedestalaver4 /= c4;
8795 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
8796 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
8797 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
8798 ampl = tool[ts_with_max_signal];
8799 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
8800 ampl += tool[ts_with_max_signal + 2];
8801 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
8802 ampl += tool[ts_with_max_signal + 1];
8803 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
8804 ampl += tool[ts_with_max_signal - 1];
8807 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
8808 linamplitudewithoutPedSubtr = lintoolwithoutPedSubtr[ts_with_max_signal];
8809 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
8810 linamplitudewithoutPedSubtr += lintoolwithoutPedSubtr[ts_with_max_signal + 2];
8811 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
8812 linamplitudewithoutPedSubtr += lintoolwithoutPedSubtr[ts_with_max_signal + 1];
8813 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
8814 linamplitudewithoutPedSubtr += lintoolwithoutPedSubtr[ts_with_max_signal - 1];
8819 if (ratio < 0. || ratio > 1.02)
8821 double aveamplitude = 0.;
8822 double aveamplitudew = 0.;
8823 if (absamplitude > 0 && timew > 0)
8824 aveamplitude = timew / absamplitude;
8830 bool anycapid =
true;
8835 for (
int ii = 0;
ii < TSsize;
ii++) {
8836 double aaaaaa = (
ii + 1) - aveamplitudew;
8837 double aaaaaa2 = aaaaaa * aaaaaa;
8838 double ampldefault = tool[
ii];
8839 rmsamp += (aaaaaa2 * ampldefault);
8840 capid = ((*digiItr)[
ii]).capid();
8841 bool er = (*digiItr)[
ii].er();
8842 bool dv = (*digiItr)[
ii].dv();
8843 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8855 if (!anycapid || anyer || !anydv)
8858 double rmsamplitude = 0.;
8861 double aveamplitude1 = aveamplitude - 1;
8867 if (sub == 1 && mdepth == 1)
8869 if (sub == 1 && mdepth == 2)
8871 if (sub == 2 && mdepth == 1)
8873 if (sub == 2 && mdepth == 2)
8875 if (sub == 2 && mdepth == 3)
8879 if (sub == 1 && mdepth == 1)
8881 if (sub == 1 && mdepth == 2)
8883 if (sub == 2 && mdepth == 1)
8885 if (sub == 2 && mdepth == 2)
8887 if (sub == 2 && mdepth == 3)
8891 for (
int ii = 0;
ii < TSsize;
ii++) {
8893 double ampldefault = tool[
ii];
8939 for (
int ii = 0;
ii < TSsize;
ii++) {
8940 double ampldefault = tool[
ii];
8947 for (
int ii = 0;
ii < TSsize;
ii++) {
8948 double ampldefault = tool[
ii];
8960 for (
int ii = 0;
ii < TSsize;
ii++) {
8961 double ampldefault = tool[
ii];
8968 for (
int ii = 0;
ii < TSsize;
ii++) {
8969 double ampldefault = tool[
ii];
8985 double mypedestal = pedestal0;
8986 double mypedestalw = pedestalw0;
9057 for (
int ii = 0;
ii < TSsize;
ii++) {
9090 if (amplitude < ADCAmplHBMin_ || amplitude >
ADCAmplHBMax_) {
9130 if (aveamplitude1 < TSmeanHBMin_ || aveamplitude1 >
TSmeanHBMax_) {
9151 if (ts_with_max_signal < TSpeakHBMin_ || ts_with_max_signal >
TSpeakHBMax_) {
9172 if (rmsamplitude < rmsHBMin_ || rmsamplitude >
rmsHBMax_) {
9231 double mypedestal = pedestal0;
9232 double mypedestalw = pedestalw0;
9303 for (
int ii = 0;
ii < TSsize;
ii++) {
9325 if (amplitude < ADCAmplHEMin_ || amplitude >
ADCAmplHEMax_) {
9386 if (aveamplitude1 < TSmeanHEMin_ || aveamplitude1 >
TSmeanHEMax_) {
9410 if (ts_with_max_signal < TSpeakHEMin_ || ts_with_max_signal >
TSpeakHEMax_) {
9434 if (rmsamplitude < rmsHEMin_ || rmsamplitude >
rmsHEMax_) {
9494 double amplitudewithPedSubtr = 0.;
9497 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
9498 amplitudewithPedSubtr = toolwithPedSubtr[ts_with_max_signal];
9499 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
9500 amplitudewithPedSubtr += toolwithPedSubtr[ts_with_max_signal + 2];
9501 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
9502 amplitudewithPedSubtr += toolwithPedSubtr[ts_with_max_signal + 1];
9503 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
9504 amplitudewithPedSubtr += toolwithPedSubtr[ts_with_max_signal - 1];
9510 if (amplitudewithPedSubtr > 50.) {
9533 int LLLLLL111111 = 0;
9534 if ((mdepth == 1 && fabs(ieta0) > 16 && fabs(ieta0) < 30))
9537 int LLLLLL222222 = 0;
9538 if ((mdepth == 2 && fabs(ieta0) > 16 && fabs(ieta0) < 27) ||
9539 (mdepth == 3 && fabs(ieta0) > 26 && fabs(ieta0) < 29))
9542 if (LLLLLL111111 == 1) {
9547 if (LLLLLL222222 == 1) {
9554 if (mdepth == 3 && fabs(ieta0) == 16) {
9583 int mdepth = hcaldetid.
depth();
9584 int sub = hcaldetid.
subdet();
9587 if (mdepth == 0 || sub > 4)
9612 coder.
adc2fC(qie11df, toolOriginal);
9613 double pedestalaver9 = 0.;
9614 double pedestalaver4 = 0.;
9615 double pedestal0 = 0.;
9616 double pedestal1 = 0.;
9617 double pedestal2 = 0.;
9618 double pedestal3 = 0.;
9619 double pedestalwaver9 = 0.;
9620 double pedestalwaver4 = 0.;
9621 double pedestalw0 = 0.;
9622 double pedestalw1 = 0.;
9623 double pedestalw2 = 0.;
9624 double pedestalw3 = 0.;
9625 double difpedestal0 = 0.;
9626 double difpedestal1 = 0.;
9627 double difpedestal2 = 0.;
9628 double difpedestal3 = 0.;
9631 double amplitude0 = 0.;
9632 double absamplitude = 0.;
9633 double tocampl = 0.;
9635 double amplitude345 = 0.;
9637 double ampl3ts = 0.;
9640 double max_signal = -100.;
9641 int ts_with_max_signal = -100;
9647 double errorBtype = 0.;
9658 for (
int ii = 0;
ii < TSsize;
ii++) {
9659 double ampldefault = 0.;
9660 double tocdefault = 0.;
9661 double ampldefault0 = 0.;
9662 double ampldefault1 = 0.;
9663 double ampldefault2 = 0.;
9669 ampldefault1 = toolOriginal[
ii];
9670 ampldefault2 = qie11df[
ii].adc();
9673 ampldefault = ampldefault0;
9676 ampldefault = ampldefault1;
9679 ampldefault = ampldefault2;
9681 tocdefault = ampldefault;
9683 int capid = (qie11df[
ii]).capid();
9685 double pedestalw = pedw->
getSigma(capid, capid);
9686 double pedestalINI = pedestal00->
getValue(capid);
9690 tool[
ii] = ampldefault;
9692 pedestalwaver9 += pedestalw * pedestalw;
9694 if (capid == 0 &&
c0 == 0) {
9699 pedestalw0 = pedestalw;
9700 pedestalwaver4 += pedestalw * pedestalw;
9701 difpedestal0 =
pedestal - pedestalINI;
9704 if (capid == 1 &&
c1 == 0) {
9709 pedestalw1 = pedestalw;
9710 pedestalwaver4 += pedestalw * pedestalw;
9711 difpedestal1 =
pedestal - pedestalINI;
9713 if (capid == 2 && c2 == 0) {
9718 pedestalw2 = pedestalw;
9719 pedestalwaver4 += pedestalw * pedestalw;
9720 difpedestal2 =
pedestal - pedestalINI;
9722 if (capid == 3 && c3 == 0) {
9727 pedestalw3 = pedestalw;
9728 pedestalwaver4 += pedestalw * pedestalw;
9729 difpedestal3 =
pedestal - pedestalINI;
9732 if (max_signal < ampldefault) {
9733 max_signal = ampldefault;
9734 ts_with_max_signal =
ii;
9737 absamplitude +=
abs(ampldefault);
9738 tocampl += tocdefault;
9740 if (
ii == 1 ||
ii == 2 ||
ii == 3 ||
ii == 4 ||
ii == 5 ||
ii == 6 ||
ii == 7 ||
ii == 8)
9741 amplitude345 += ampldefault;
9746 timew += (
ii + 1) *
abs(ampldefault);
9747 timeww += (
ii + 1) * ampldefault;
9749 if (
ii == 2 && ampldefault > 0.)
9755 pedestalaver9 /= TSsize;
9756 pedestalaver4 /= c4;
9757 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
9758 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
9761 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize) {
9762 ampl = tool[ts_with_max_signal];
9763 ampl3ts = tool[ts_with_max_signal];
9765 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize) {
9766 ampl += tool[ts_with_max_signal - 1];
9767 ampl3ts += tool[ts_with_max_signal - 1];
9769 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize) {
9770 ampl += tool[ts_with_max_signal + 1];
9771 ampl3ts += tool[ts_with_max_signal + 1];
9773 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize) {
9774 ampl += tool[ts_with_max_signal + 2];
9780 double xb = ampl / 40.;
9781 double xc = amplitude345 / 40.;
9782 double xd = ampl3ts / 40.;
9783 double txa = tocampl / 40.;
9785 if (((
ieta == -16 ||
ieta == 15) && mdepth == 4) ||
9786 ((
ieta == -17 ||
ieta == 16) && (mdepth == 2 || mdepth == 3)) ||
9787 ((
ieta == -18 ||
ieta == 17) && mdepth == 5)) {
9788 double c0 = 1.000000;
9789 double b1 = 2.59096e-05;
9790 double a2 = 4.60721e-11;
9791 double corrforxa =
a2 * xa * xa +
b1 * xa +
c0;
9792 double corrforxb =
a2 * xb * xb +
b1 * xb +
c0;
9793 double corrforxc =
a2 * xc * xc +
b1 * xc +
c0;
9794 double corrforxd =
a2 * xd * xd +
b1 * xd +
c0;
9795 double corrfortxa =
a2 * txa * txa +
b1 * txa +
c0;
9800 amplitude345 *= corrforxc;
9801 ampl3ts *= corrforxd;
9802 tocampl *= corrfortxa;
9806 double c0 = 1.000000;
9807 double b1 = 2.71238e-05;
9808 double a2 = 1.32877e-10;
9809 double corrforxa =
a2 * xa * xa +
b1 * xa +
c0;
9810 double corrforxb =
a2 * xb * xb +
b1 * xb +
c0;
9811 double corrforxc =
a2 * xc * xc +
b1 * xc +
c0;
9812 double corrforxd =
a2 * xd * xd +
b1 * xd +
c0;
9813 double corrfortxa =
a2 * txa * txa +
b1 * txa +
c0;
9818 amplitude345 *= corrforxc;
9819 ampl3ts *= corrforxd;
9820 tocampl *= corrfortxa;
9826 if (sub == 1 || sub == 2) {
9837 if (ratio < 0. || ratio > 1.02)
9839 double aveamplitude = 0.;
9840 double aveamplitudew = 0.;
9841 if (absamplitude > 0 && timew > 0)
9842 aveamplitude = timew / absamplitude;
9843 if (amplitude0 > 0 && timeww > 0)
9844 aveamplitudew = timeww / amplitude0;
9848 bool anycapid =
true;
9851 for (
int ii = 0;
ii < TSsize;
ii++) {
9852 double aaaaaa = (
ii + 1) - aveamplitudew;
9853 double aaaaaa2 = aaaaaa * aaaaaa;
9854 double ampldefault = tool[
ii];
9855 rmsamp += (aaaaaa2 * ampldefault);
9856 capid = (qie11df[
ii]).capid();
9857 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
9865 double rmsamplitude = 0.;
9866 if ((amplitude0 > 0 && rmsamp > 0) || (amplitude0 < 0 && rmsamp < 0))
9867 rmsamplitude =
sqrt(rmsamp / amplitude0);
9868 double aveamplitude1 = aveamplitude - 1;
9873 if (sub == 1 && mdepth == 1)
9875 if (sub == 1 && mdepth == 2)
9877 if (sub == 2 && mdepth == 1)
9879 if (sub == 2 && mdepth == 2)
9881 if (sub == 2 && mdepth == 3)
9885 if (sub == 1 && mdepth == 1)
9887 if (sub == 1 && mdepth == 2)
9889 if (sub == 2 && mdepth == 1)
9891 if (sub == 2 && mdepth == 2)
9893 if (sub == 2 && mdepth == 3)
9897 for (
int ii = 0;
ii < TSsize;
ii++) {
9899 double ampldefault = tool[
ii];
9902 if (amplitude0 > 120) {
9945 for (
int ii = 0;
ii < TSsize;
ii++) {
9946 double ampldefault = tool[
ii];
9953 for (
int ii = 0;
ii < TSsize;
ii++) {
9954 double ampldefault = tool[
ii];
9968 for (
int ii = 0;
ii < TSsize;
ii++) {
9969 double ampldefault = tool[
ii];
9976 for (
int ii = 0;
ii < TSsize;
ii++) {
9977 double ampldefault = tool[
ii];
9993 double mypedestal = pedestal0;
9994 double mypedestalw = pedestalw0;
10065 for (
int ii = 0;
ii < TSsize;
ii++) {
10098 if (amplitude < ADCAmplHBMin_ || amplitude >
ADCAmplHBMax_) {
10156 if (aveamplitude1 < TSmeanHBMin_ || aveamplitude1 >
TSmeanHBMax_) {
10184 h_TSmaxA_HB->Fill(
float(ts_with_max_signal), 1.);
10185 if (ts_with_max_signal < TSpeakHBMin_ || ts_with_max_signal >
TSpeakHBMax_) {
10214 if (rmsamplitude < rmsHBMin_ || rmsamplitude >
rmsHBMax_) {
10243 if (ratio < ratioHBMin_ || ratio >
ratioHBMax_) {
10292 if (flagTS2 == 1) {
10311 double mypedestal = pedestal0;
10312 double mypedestalw = pedestalw0;
10399 for (
int ii = 0;
ii < TSsize;
ii++) {
10417 if (amplitude < ADCAmplHEMin_ || amplitude >
ADCAmplHEMax_) {
10465 if (ampl3ts > 70 && ampl3ts < 110) {
10468 }
else if (ampl3ts > 110 && ampl3ts < 150) {
10472 if (ampl > 87 && ampl < 127) {
10475 }
else if (ampl > 127 && ampl < 167) {
10538 if (aveamplitude1 < TSmeanHEMin_ || aveamplitude1 >
TSmeanHEMax_) {
10577 h_TSmaxA_HE->Fill(
float(ts_with_max_signal), 1.);
10578 if (ts_with_max_signal < TSpeakHEMin_ || ts_with_max_signal >
TSpeakHEMax_) {
10618 if (rmsamplitude < rmsHEMin_ || rmsamplitude >
rmsHEMax_) {
10659 if (ratio < ratioHEMin_ || ratio >
ratioHEMax_) {
10729 if (flagTS2 == 1) {
10754 int mdepth = cell.
depth();
10755 int iphi = cell.iphi() - 1;
10756 int ieta = cell.ieta();
10759 int sub = cell.subdet();
10769 coder.
adc2fC(*digiItr, toolOriginal);
10770 double pedestalaver9 = 0.;
10771 double pedestalaver4 = 0.;
10772 double pedestal0 = 0.;
10773 double pedestal1 = 0.;
10774 double pedestal2 = 0.;
10775 double pedestal3 = 0.;
10776 double pedestalwaver9 = 0.;
10777 double pedestalwaver4 = 0.;
10778 double pedestalw0 = 0.;
10779 double pedestalw1 = 0.;
10780 double pedestalw2 = 0.;
10781 double pedestalw3 = 0.;
10783 double absamplitude = 0.;
10786 double timeww = 0.;
10787 double max_signal = -100.;
10788 int ts_with_max_signal = -100;
10794 double errorBtype = 0.;
10796 if ((*digiItr).size() != TSsize)
10798 TSsize = digiItr->size();
10799 for (
int ii = 0;
ii < TSsize;
ii++) {
10801 double ampldefault = 0.;
10802 double ampldefault0 = 0.;
10803 double ampldefault1 = 0.;
10804 double ampldefault2 = 0.;
10805 ampldefault0 =
adc2fC[digiItr->sample(
ii).adc()];
10807 ampldefault1 = toolOriginal[
ii];
10808 ampldefault2 = (*digiItr)[
ii].adc();
10810 ampldefault = ampldefault0;
10813 ampldefault = ampldefault1;
10816 ampldefault = ampldefault2;
10819 int capid = ((*digiItr)[
ii]).capid();
10822 double pedestalw = pedw->
getSigma(capid, capid);
10826 tool[
ii] = ampldefault;
10829 pedestalwaver9 += pedestalw * pedestalw;
10831 if (capid == 0 &&
c0 == 0) {
10836 pedestalw0 = pedestalw;
10837 pedestalwaver4 += pedestalw * pedestalw;
10840 if (capid == 1 &&
c1 == 0) {
10845 pedestalw1 = pedestalw;
10846 pedestalwaver4 += pedestalw * pedestalw;
10848 if (capid == 2 && c2 == 0) {
10853 pedestalw2 = pedestalw;
10854 pedestalwaver4 += pedestalw * pedestalw;
10856 if (capid == 3 && c3 == 0) {
10861 pedestalw3 = pedestalw;
10862 pedestalwaver4 += pedestalw * pedestalw;
10865 if (max_signal < ampldefault) {
10866 max_signal = ampldefault;
10867 ts_with_max_signal =
ii;
10870 absamplitude +=
abs(ampldefault);
10875 timew += (
ii + 1) *
abs(ampldefault);
10876 timeww += (
ii + 1) * ampldefault;
10886 pedestalaver9 /= TSsize;
10887 pedestalaver4 /= c4;
10888 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
10889 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
10892 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
10893 ampl = tool[ts_with_max_signal];
10894 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
10895 ampl += tool[ts_with_max_signal + 2];
10896 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
10897 ampl += tool[ts_with_max_signal + 1];
10898 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
10899 ampl += tool[ts_with_max_signal - 1];
10906 if (ratio < 0. || ratio > 1.02)
10909 double aveamplitude = 0.;
10910 double aveamplitudew = 0.;
10911 if (absamplitude > 0 && timew > 0)
10912 aveamplitude = timew / absamplitude;
10916 double rmsamp = 0.;
10919 bool anycapid =
true;
10920 bool anyer =
false;
10924 for (
int ii = 0;
ii < TSsize;
ii++) {
10925 double aaaaaa = (
ii + 1) - aveamplitudew;
10926 double aaaaaa2 = aaaaaa * aaaaaa;
10927 double ampldefault = tool[
ii];
10928 rmsamp += (aaaaaa2 * ampldefault);
10929 capid = ((*digiItr)[
ii]).capid();
10930 bool er = (*digiItr)[
ii].er();
10931 bool dv = (*digiItr)[
ii].dv();
10932 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
10945 if (!anycapid || anyer || !anydv)
10947 double rmsamplitude = 0.;
10950 double aveamplitude1 = aveamplitude - 1;
10952 if (sub == 4 && mdepth == 1)
10954 if (sub == 4 && mdepth == 2)
10958 if (sub == 4 && mdepth == 1)
10960 if (sub == 4 && mdepth == 2)
10972 for (
int ii = 0;
ii < TSsize;
ii++) {
10973 double ampldefault = tool[
ii];
10980 for (
int ii = 0;
ii < TSsize;
ii++) {
10981 double ampldefault = tool[
ii];
10995 double mypedestal = pedestal0;
10996 double mypedestalw = pedestalw0;
11056 for (
int ii = 0;
ii < TSsize;
ii++) {
11073 if (amplitude < ADCAmplHFMin_ || amplitude >
ADCAmplHFMax_) {
11115 if (aveamplitude1 < TSmeanHFMin_ || aveamplitude1 >
TSmeanHFMax_) {
11135 h_TSmaxA_HF->Fill(
float(ts_with_max_signal), 1.);
11136 if (ts_with_max_signal < TSpeakHFMin_ || ts_with_max_signal >
TSpeakHFMax_) {
11156 if (rmsamplitude < rmsHFMin_ || rmsamplitude >
rmsHFMax_) {
11176 if (ratio < ratioHFMin_ || ratio >
ratioHFMax_) {
11224 int mdepth = hcaldetid.
depth();
11225 int sub = hcaldetid.
subdet();
11240 if (mdepth == 0 || sub != 4)
11260 coder.
adc2fC(qie10df, toolOriginal);
11263 double pedestalaver9 = 0.;
11264 double pedestalaver4 = 0.;
11265 double pedestal0 = 0.;
11266 double pedestal1 = 0.;
11267 double pedestal2 = 0.;
11268 double pedestal3 = 0.;
11269 double pedestalwaver9 = 0.;
11270 double pedestalwaver4 = 0.;
11271 double pedestalw0 = 0.;
11272 double pedestalw1 = 0.;
11273 double pedestalw2 = 0.;
11274 double pedestalw3 = 0.;
11276 double absamplitude = 0.;
11279 double timeww = 0.;
11280 double max_signal = -100.;
11281 int ts_with_max_signal = -100;
11287 double errorBtype = 0.;
11294 for (
int ii = 0;
ii < TSsize;
ii++) {
11295 double ampldefault = 0.;
11296 double ampldefault0 = 0.;
11297 double ampldefault1 = 0.;
11298 double ampldefault2 = 0.;
11301 ampldefault1 = toolOriginal[
ii];
11302 ampldefault2 = qie10df[
ii].adc();
11304 ampldefault = ampldefault0;
11307 ampldefault = ampldefault1;
11310 ampldefault = ampldefault2;
11313 int capid = (qie10df[
ii]).capid();
11315 double pedestalw = pedw->
getSigma(capid, capid);
11320 tool[
ii] = ampldefault;
11323 pedestalwaver9 += pedestalw * pedestalw;
11325 if (capid == 0 &&
c0 == 0) {
11330 pedestalw0 = pedestalw;
11331 pedestalwaver4 += pedestalw * pedestalw;
11334 if (capid == 1 &&
c1 == 0) {
11339 pedestalw1 = pedestalw;
11340 pedestalwaver4 += pedestalw * pedestalw;
11342 if (capid == 2 && c2 == 0) {
11347 pedestalw2 = pedestalw;
11348 pedestalwaver4 += pedestalw * pedestalw;
11350 if (capid == 3 && c3 == 0) {
11355 pedestalw3 = pedestalw;
11356 pedestalwaver4 += pedestalw * pedestalw;
11359 if (max_signal < ampldefault) {
11360 max_signal = ampldefault;
11361 ts_with_max_signal =
ii;
11364 absamplitude +=
abs(ampldefault);
11366 timew += (
ii + 1) *
abs(ampldefault);
11367 timeww += (
ii + 1) * ampldefault;
11368 if (
ii == 1 && ampldefault > 0.)
11379 pedestalaver9 /= TSsize;
11380 pedestalaver4 /= c4;
11381 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
11382 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
11385 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
11386 ampl = tool[ts_with_max_signal];
11387 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
11388 ampl += tool[ts_with_max_signal + 2];
11389 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
11390 ampl += tool[ts_with_max_signal + 1];
11391 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
11392 ampl += tool[ts_with_max_signal - 1];
11398 if (ratio < 0. || ratio > 1.02)
11400 double aveamplitude = 0.;
11401 double aveamplitudew = 0.;
11402 if (absamplitude > 0 && timew > 0)
11403 aveamplitude = timew / absamplitude;
11407 double rmsamp = 0.;
11409 bool anycapid =
true;
11412 for (
int ii = 0;
ii < TSsize;
ii++) {
11413 double aaaaaa = (
ii + 1) - aveamplitudew;
11414 double aaaaaa2 = aaaaaa * aaaaaa;
11415 double ampldefault = tool[
ii];
11416 rmsamp += (aaaaaa2 * ampldefault);
11417 capid = (qie10df[
ii]).capid();
11418 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
11426 double rmsamplitude = 0.;
11429 double aveamplitude1 = aveamplitude - 1;
11434 if (sub == 4 && mdepth == 1)
11436 if (sub == 4 && mdepth == 2)
11440 if (sub == 4 && mdepth == 1)
11442 if (sub == 4 && mdepth == 2)
11454 for (
int ii = 0;
ii < TSsize;
ii++) {
11455 double ampldefault = tool[
ii];
11462 for (
int ii = 0;
ii < TSsize;
ii++) {
11463 double ampldefault = tool[
ii];
11478 double mypedestal = pedestal0;
11479 double mypedestalw = pedestalw0;
11547 for (
int ii = 0;
ii < TSsize;
ii++) {
11564 if (amplitude < ADCAmplHFMin_ || amplitude >
ADCAmplHFMax_) {
11620 if (aveamplitude1 < TSmeanHFMin_ || aveamplitude1 >
TSmeanHFMax_) {
11648 h_TSmaxA_HF->Fill(
float(ts_with_max_signal), 1.);
11649 if (ts_with_max_signal < TSpeakHFMin_ || ts_with_max_signal >
TSpeakHFMax_) {
11678 if (rmsamplitude < rmsHFMin_ || rmsamplitude >
rmsHFMax_) {
11707 if (ratio < ratioHFMin_ || ratio >
ratioHFMax_) {
11756 if (flagTS1 == 1) {
11777 int mdepth = cell.
depth();
11778 int iphi = cell.iphi() - 1;
11779 int ieta = cell.ieta();
11782 int sub = cell.subdet();
11792 coder.
adc2fC(*digiItr, toolOriginal);
11793 double pedestalaver9 = 0.;
11794 double pedestalaver4 = 0.;
11795 double pedestal0 = 0.;
11796 double pedestal1 = 0.;
11797 double pedestal2 = 0.;
11798 double pedestal3 = 0.;
11799 double pedestalwaver9 = 0.;
11800 double pedestalwaver4 = 0.;
11801 double pedestalw0 = 0.;
11802 double pedestalw1 = 0.;
11803 double pedestalw2 = 0.;
11804 double pedestalw3 = 0.;
11806 double absamplitude = 0.;
11809 double timeww = 0.;
11810 double max_signal = -100.;
11811 int ts_with_max_signal = -100;
11817 double errorBtype = 0.;
11819 if ((*digiItr).size() != TSsize)
11821 TSsize = digiItr->size();
11823 for (
int ii = 0;
ii < TSsize;
ii++) {
11824 double ampldefault = 0.;
11825 double ampldefault0 = 0.;
11826 double ampldefault1 = 0.;
11827 double ampldefault2 = 0.;
11828 ampldefault0 =
adc2fC[digiItr->sample(
ii).adc()];
11830 ampldefault1 = toolOriginal[
ii];
11831 ampldefault2 = (*digiItr)[
ii].adc();
11833 ampldefault = ampldefault0;
11836 ampldefault = ampldefault1;
11839 ampldefault = ampldefault2;
11841 int capid = ((*digiItr)[
ii]).capid();
11843 double pedestalw = pedw->
getSigma(capid, capid);
11846 tool[
ii] = ampldefault;
11848 pedestalwaver9 += pedestalw * pedestalw;
11849 if (capid == 0 &&
c0 == 0) {
11854 pedestalw0 = pedestalw;
11855 pedestalwaver4 += pedestalw * pedestalw;
11857 if (capid == 1 &&
c1 == 0) {
11862 pedestalw1 = pedestalw;
11863 pedestalwaver4 += pedestalw * pedestalw;
11865 if (capid == 2 && c2 == 0) {
11870 pedestalw2 = pedestalw;
11871 pedestalwaver4 += pedestalw * pedestalw;
11873 if (capid == 3 && c3 == 0) {
11878 pedestalw3 = pedestalw;
11879 pedestalwaver4 += pedestalw * pedestalw;
11882 if (max_signal < ampldefault) {
11883 max_signal = ampldefault;
11884 ts_with_max_signal =
ii;
11887 absamplitude +=
abs(ampldefault);
11891 timew += (
ii + 1) *
abs(ampldefault);
11892 timeww += (
ii + 1) * ampldefault;
11893 if (ii < 3 && ampldefault > 0.)
11898 pedestalaver9 /= TSsize;
11899 pedestalaver4 /= c4;
11900 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
11901 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
11902 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
11903 ampl = tool[ts_with_max_signal];
11904 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
11905 ampl += tool[ts_with_max_signal + 2];
11906 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
11907 ampl += tool[ts_with_max_signal + 1];
11908 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
11909 ampl += tool[ts_with_max_signal - 1];
11913 if (ratio < 0. || ratio > 1.04)
11915 double aveamplitude = 0.;
11916 double aveamplitudew = 0.;
11917 if (absamplitude > 0 && timew > 0)
11918 aveamplitude = timew / absamplitude;
11921 double rmsamp = 0.;
11923 bool anycapid =
true;
11924 bool anyer =
false;
11928 for (
int ii = 0;
ii < TSsize;
ii++) {
11929 double aaaaaa = (
ii + 1) - aveamplitudew;
11930 double aaaaaa2 = aaaaaa * aaaaaa;
11931 double ampldefault = tool[
ii];
11932 rmsamp += (aaaaaa2 * ampldefault);
11933 capid = ((*digiItr)[
ii]).capid();
11934 bool er = (*digiItr)[
ii].er();
11935 bool dv = (*digiItr)[
ii].dv();
11936 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
11948 if (!anycapid || anyer || !anydv)
11950 double rmsamplitude = 0.;
11953 double aveamplitude1 = aveamplitude - 1;
11955 if (sub == 3 && mdepth == 4)
11959 if (sub == 3 && mdepth == 4)
11970 for (
int ii = 0;
ii < TSsize;
ii++) {
11971 double ampldefault = tool[
ii];
11976 for (
int ii = 0;
ii < TSsize;
ii++) {
11977 double ampldefault = tool[
ii];
11991 double mypedestal = pedestal0;
11992 double mypedestalw = pedestalw0;
12035 for (
int ii = 0;
ii < TSsize;
ii++) {
12051 if (amplitude < ADCAmplHOMin_ || amplitude >
ADCAmplHOMax_) {
12083 if (aveamplitude1 < TSmeanHOMin_ || aveamplitude1 >
TSmeanHOMax_) {
12095 h_TSmaxA_HO->Fill(
float(ts_with_max_signal), 1.);
12096 if (ts_with_max_signal < TSpeakHOMin_ || ts_with_max_signal >
TSpeakHOMax_) {
12110 if (rmsamplitude < rmsHOMin_ || rmsamplitude >
rmsHOMax_) {
12123 if (ratio < ratioHOMin_ || ratio >
ratioHOMax_) {
12140 if (flagTS012 == 1)
12147 if (kdet == 1 || kdet == 2) {
12151 cal_RBX = (kphi + 1) / 4;
12152 cal_RBX = cal_RBX + 18 * (keta + 1) / 2;
12155 cal_RBX = (
int)(kphi / 18) + 1;
12162 cal_RBX = kphi / 12 + 1;
12168 cal_RBX = kphi / 12 + 1 + 6;
12174 cal_RBX = kphi / 6 + 1 + 12;
12180 cal_RBX = kphi / 12 + 1 + 24;
12186 cal_RBX = kphi / 12 + 1 + 30;
12207 for (
int k1 = 0; k1 <
ndepth; k1++) {
12208 for (
int k2 = 0; k2 <
neta; k2++) {
12209 for (
int k3 = 0; k3 <
nphi; k3++) {
12210 int ieta = k2 - 41;
12305 int kkkk2 = (k2 - 1) / 4;
12310 int kkkk3 = (k3) / 4 + 1;
12312 ietaphi = ((kkkk2)-1) *
znphi + (kkkk3);
12315 bbb3 = bbbc / bbb1;
13159 std::cout <<
" ==== Edn of run " << std::endl;
13177 MAPfile <<
"#define LogEleMapdb_h" << std::endl;
13178 MAPfile <<
"#include <algorithm>" << std::endl;
13179 MAPfile <<
"#include <iostream>" << std::endl;
13180 MAPfile <<
"#include <vector>" << std::endl;
13181 MAPfile <<
"#include <string>" << std::endl;
13182 MAPfile <<
"#include <sstream>" << std::endl;
13185 MAPfile <<
"struct Cell {" << std::endl;
13186 MAPfile <<
" std::string subdet;" << std::endl;
13187 MAPfile <<
" int Eta;" << std::endl;
13188 MAPfile <<
" int Phi;" << std::endl;
13189 MAPfile <<
" int Depth;" << std::endl;
13190 MAPfile <<
" std::string RBX;" << std::endl;
13191 MAPfile <<
" int RM;" << std::endl;
13192 MAPfile <<
" int Pixel;" << std::endl;
13193 MAPfile <<
" int RMfiber;" << std::endl;
13194 MAPfile <<
" int FiberCh;" << std::endl;
13195 MAPfile <<
" int QIE;" << std::endl;
13196 MAPfile <<
" int ADC;" << std::endl;
13197 MAPfile <<
" int VMECardID;" << std::endl;
13198 MAPfile <<
" int dccID;" << std::endl;
13199 MAPfile <<
" int Spigot;" << std::endl;
13200 MAPfile <<
" int FiberIndex;" << std::endl;
13201 MAPfile <<
" int HtrSlot;" << std::endl;
13202 MAPfile <<
" int HtrTB;" << std::endl;
13205 MAPfile <<
"// the function check, if \"par\" == \"val\" for this cell" << std::endl;
13206 MAPfile <<
" bool check(const std::string par, const int val) const " << std::endl;
13207 MAPfile <<
" {" << std::endl;
13208 MAPfile <<
" if (par == \"Eta\") return (val == Eta);" << std::endl;
13209 MAPfile <<
" else if (par == \"Phi\") return (val == Phi);" << std::endl;
13210 MAPfile <<
" else if (par == \"Depth\") return (val == Depth);" << std::endl;
13211 MAPfile <<
" else if (par == \"RM\") return (val == RM);" << std::endl;
13212 MAPfile <<
" else if (par == \"Pixel\") return (val == Pixel);" << std::endl;
13213 MAPfile <<
" else if (par == \"RMfiber\") return (val == RMfiber);" << std::endl;
13214 MAPfile <<
" else if (par == \"FiberCh\") return (val == FiberCh);" << std::endl;
13215 MAPfile <<
" else if (par == \"QIE\") return (val == QIE);" << std::endl;
13216 MAPfile <<
" else if (par == \"ADC\") return (val == ADC);" << std::endl;
13217 MAPfile <<
" else if (par == \"VMECardID\") return (val == VMECardID);" << std::endl;
13218 MAPfile <<
" else if (par == \"dccID\") return (val == dccID);" << std::endl;
13219 MAPfile <<
" else if (par == \"Spigot\") return (val == Spigot);" << std::endl;
13220 MAPfile <<
" else if (par == \"FiberIndex\") return (val == FiberIndex);" << std::endl;
13221 MAPfile <<
" else if (par == \"HtrSlot\") return (val == HtrSlot);" << std::endl;
13222 MAPfile <<
" else if (par == \"HtrTB\") return (val == HtrTB);" << std::endl;
13223 MAPfile <<
" else return false;" << std::endl;
13224 MAPfile <<
" }" << std::endl;
13227 MAPfile <<
" bool check(const std::string par, const std::string val) const" << std::endl;
13228 MAPfile <<
" {" << std::endl;
13229 MAPfile <<
" if (par == \"subdet\") return (val == subdet);" << std::endl;
13230 MAPfile <<
" else if (par == \"RBX\") return (val == RBX);" << std::endl;
13231 MAPfile <<
" else return false;" << std::endl;
13232 MAPfile <<
" }" << std::endl;
13234 MAPfile <<
"};" << std::endl;
13237 MAPfile <<
"const Cell AllCells[] = {" << std::endl;
13238 MAPfile <<
"//{ SD, Eta, Phi, Depth, RBX, RM, PIXEL, RMfiber, Fiber Ch., QIE, ADC, VMECrateId, dccid, spigot, " 13239 "fiberIndex, htrSlot, htrTopBottom }" 13252 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13253 << detid->
depth() <<
" ,";
13254 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13255 << lmap_entry.
rmFiber() <<
" , ";
13259 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13261 MAPfile <<
"}," << std::endl;
13274 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() - 1 <<
" , " << detid->
iphi() - 1 <<
" , " 13275 << detid->
depth() <<
" ,";
13276 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13277 << lmap_entry.
rmFiber() <<
" , ";
13281 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13283 MAPfile <<
"}," << std::endl;
13290 for (
int eta = -20;
eta <= -20;
eta++) {
13298 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13299 << detid->
depth() <<
" ,";
13300 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13301 << lmap_entry.
rmFiber() <<
" , ";
13305 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13307 MAPfile <<
"}," << std::endl;
13313 for (
int eta = -19;
eta <= -16;
eta++) {
13321 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13322 << detid->
depth() <<
" ,";
13323 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13324 << lmap_entry.
rmFiber() <<
" , ";
13328 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13330 MAPfile <<
"}," << std::endl;
13335 for (
int eta = -29;
eta <= -16;
eta++) {
13343 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13344 << detid->
depth() <<
" ,";
13345 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13346 << lmap_entry.
rmFiber() <<
" , ";
13350 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13352 MAPfile <<
"}," << std::endl;
13365 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13366 << detid->
depth() <<
" ,";
13367 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13368 << lmap_entry.
rmFiber() <<
" , ";
13372 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13374 MAPfile <<
"}," << std::endl;
13382 for (
int eta = -41;
eta <= -29;
eta++) {
13390 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13391 << detid->
depth() <<
" ,";
13392 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13393 << lmap_entry.
rmFiber() <<
" , ";
13397 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13399 MAPfile <<
"}," << std::endl;
13413 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13414 << detid->
depth() <<
" ,";
13415 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13416 << lmap_entry.
rmFiber() <<
" , ";
13420 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13422 MAPfile <<
"}," << std::endl;
13438 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13439 << detid->
depth() <<
" ,";
13440 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13441 << lmap_entry.
rmFiber() <<
" , ";
13445 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13447 MAPfile <<
"}," << std::endl;
13461 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() - 1 <<
" , " << detid->
iphi() - 1 <<
" , " 13462 << detid->
depth() <<
" ,";
13463 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13464 << lmap_entry.
rmFiber() <<
" , ";
13468 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13470 MAPfile <<
"}," << std::endl;
13475 MAPfile <<
"};" << std::endl;
13478 MAPfile <<
"// macro for array length calculation" << std::endl;
13479 MAPfile <<
"#define DIM(a) (sizeof(a)/sizeof(a[0]))" << std::endl;
13482 MAPfile <<
"// class for cells array managing" << std::endl;
13483 MAPfile <<
"class CellDB {" << std::endl;
13484 MAPfile <<
"public:" << std::endl;
13485 MAPfile <<
" CellDB()" << std::endl;
13486 MAPfile <<
" : cells(AllCells, AllCells + DIM(AllCells))" << std::endl;
13487 MAPfile <<
"{}" << std::endl;
13490 MAPfile <<
"// return i-th cell" << std::endl;
13491 MAPfile <<
"Cell operator [] (int i) const {return cells[i];}" << std::endl;
13493 MAPfile <<
"// number of cells in database" << std::endl;
13494 MAPfile <<
"int size() const {return cells.size();}" << std::endl;
13497 MAPfile <<
"// select cells for which \"par\" == \"val\"" << std::endl;
13498 MAPfile <<
"template<typename T>" << std::endl;
13499 MAPfile <<
"CellDB find(const std::string par, const T val) const" << std::endl;
13501 MAPfile <<
" std::vector<Cell> s;" << std::endl;
13502 MAPfile <<
" for (size_t i = 0; i < cells.size(); ++i)" << std::endl;
13503 MAPfile <<
" if (cells[i].check(par, val))" << std::endl;
13504 MAPfile <<
" s.push_back(cells[i]);" << std::endl;
13505 MAPfile <<
" return CellDB(s);" << std::endl;
13506 MAPfile <<
"} " << std::endl;
13508 MAPfile <<
"private:" << std::endl;
13509 MAPfile <<
" CellDB( const std::vector<Cell> s)" << std::endl;
13510 MAPfile <<
" : cells(s)" << std::endl;
13511 MAPfile <<
"{}" << std::endl;
13512 MAPfile <<
"std::vector<Cell> cells;" << std::endl;
13513 MAPfile <<
"};" << std::endl;
13515 std::cout <<
"===== Finish writing Channel MAP =====" << std::endl;
13519 double deltaphi = phi1 - phi2;
13521 deltaphi = phi2 - phi1;
13523 if (deltaphi >
M_PI) {
13524 deltaphi = 2. *
M_PI - deltaphi;
13527 double tmp =
sqrt(deltaeta * deltaeta + deltaphi * deltaphi);
13537 if (a1 > 0.5 *
M_PI &&
a2 < 0.)
13539 if (
a2 > 0.5 *
M_PI && a1 < 0.)
13541 double tmp = (a1 * en1 +
a2 * en2) / (en1 + en2);
13554 double tmp =
a2 - a1;
13555 if (a1 *
a2 < 0.) {
13556 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_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]
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
constexpr int readoutVMECrateId() const
get the readout VME crate number
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_
TH2F * h_recSignalEnergy1_HE5
TH2F * h_mapDepth2AmplE34_HE
TH2F * h_recSignalEnergy0_HF1
TH1F * h_ADCAmpl345Zoom_HB
TH2F * h_mapTSmeanCalib_HF
TH1F * h_averNOSIGNALsumamplitude_HO
TH2F * h_mapDepth1Ped2_HE
TH1F * h_runbadrate0_depth4_HO
TH1F * h_bcnbadrate0_depth1_HF
TH2F * h_mapGetRMSOverNormalizedSignal0_HE
TH2F * h_amplitudechannel0_HE3
TH2F * h_mapDepth5TSmaxA_HE
TH2F * h_2DsumPedestalLS8
TH1F * h_ndepthesperamplitudebins_HE
const HcalDbService * conditions
TH1F * h_AmplitudeHEtest6
TH1F * h_pedestalwaver9_HE
TH1F * h_sumTSmeanAperLS7
TH2F * h_2DsumADCAmplLS1_LSselected
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
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
TH2F * h_mapDepth4ADCAmpl225Copy_HB
double averSIGNALsumamplitude_HO
TH2F * h_2DsumADCAmplLSdepth3HFu
void fillDigiErrors(HBHEDigiCollection::const_iterator &digiItr)
TH2F * h_recNoiseEnergy1_HE6
double lsdep_estimator3_HBdepth1_
TH1F * h_numberofhitsHEtest
double sum0Estimator[nsub][ndepth][neta][nphi]
TH2F * h_mapAiForLS_bad_HE
TH1F * h_AmplitudeHEtest1
double lsdep_estimator1_HFdepth3_
TH2F * h_2DsumADCAmplLSdepth5HEu
TH1F * h_sumamplitude_depth4_HO
TH1F * h_pedestalaver9_HO
TH2F * h_2D0sumTSmeanALS7
TH2F * h_2DsumAmplitudeLS1
TH1F * h_sum0ADCAmplperLS7u
TH1F * h_runnbadchannels_depth3_HE
TH2F * h_amplitudechannel2_HF1
TH2F * h_mapDepth1ADCAmpl12_HF
TH1F * h_AmplitudeHBrest6
TH2F * h_mapDepth4AmplE34_HE
TH2F * h_mapDepth3ADCAmpl12_HF
TH1F * h_sigLayer1RADDAM6_HED2
double mapRADDAM_HED2[ndepth][neta]
TH2F * h_maprphinorm0_HE4
TH2F * h_mapDepth1Pedw1_HE
TH2F * h_2D0sumPedestalLS1
TH2F * h_mapAiForLS_good_HE
TH2F * h_2DsumADCAmplLSdepth4HBu
TH2F * h_mapDepth2ADCAmpl225_HE
TH1F * h_bcnnbadchannels_depth2_HF
TH2F * h_mapDepth1Ped1_HE
TH2F * h_mapDepth2ADCAmpl225_HB
TH2F * h_2DsumPedestalLS4
TH2F * h_recSignalEnergy1_HE2
TH2F * h_2D0sumADCAmplLS1
TH1F * h_A_Depth2RADDAM_HE
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
double phi12(double phi1, double en1, double phi2, double en2)
TH1F * h_energyhitNoise_HE
double lsdep_estimator1_HBdepth1_
bool studyRatioShapeHist_
edm::EDGetTokenT< HBHEDigiCollection > tok_hbhe_
TH2F * h_mapDepth4Error_HB
std::vector< T >::const_iterator const_iterator
TH2F * h_2DsumPedestalLS3
TH2F * h_FullSignal3D0_HO
TH1F * h_sumTSmeanAperLS3
TH2F * h_2D0Ataildepth1_HB
TH2F * h_mapDepth4Ampl_HF
int lsdep_cut3_max_HEdepth3_
TH1F * h_sum0ADCAmplperLSdepth5HEu
TH2F * h_mapDepth3ADCAmpl_HF
TH1F * h_gsmdifferencefit3_HE
const HcalQIEShape * shape
TH1F * h_sumamplitude_depth2_HE
double lsdep_estimator1_HFdepth4_
TH2F * h_mapDepth2linADCAmpl12_HE
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_
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_
TH2F * h_mapDepth1ADCAmpl225_HF
TH1F * h_Amplitude_forCapIdErrors_HE1
int lsdep_cut1_peak_HFdepth1_
TH2F * h_mapDepth4Amplitude225_HF
TH2F * h_mapDepth2pedestal_HB
double lsdep_estimator3_HEdepth2_
TH2F * h_mapGetRMSOverNormalizedSignal_HO
TH2F * h_mapDepth1Pedw3_HE
bool studyTSmeanShapeHist_
void fillDigiAmplitudeQIE11(QIE11DataFrame qie11df)
TH2F * h_2D0sumADCAmplLSdepth4HFu
double calibrTSmeanHFMin_
TH1F * h_sumamplitude_depth1_HB
constexpr int htrSlot() const
get the htr slot
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
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
TH2F * h_mapDepth4ADCAmpl225Copy_HO
TH2F * h_mapenophinorm_HE2
TH2F * h_2DsumADCAmplLSdepth4HEu
TH1F * h_sum0PedestalperLS6
TH2F * h_mapDepth4pedestalw_HF
TH1F * h_sum0ADCAmplperLS3_M2
TH2F * h_amplitudechannel0_HB1
int lsdep_cut1_peak_HBdepth1_
TH2F * h_recSignalEnergy2_HF2
TH1F * h_sumCutTSmeanAperLS7
TH2F * h_recSignalEnergy0_HF3
TH2F * h_amplitudechannel1_HE4
TH2F * h_recSignalEnergy1_HE7
TH1F * h_sumADCAmplperLS6_P2
TH2F * h_mapADCAmpl11secondpeak_HE
TH1F * h_nbadchannels_depth4_HO
TH2F * h_mapWidthCalib047_HB
TH2F * h_amplitudechannel2_HE5
TH1F * h_sumPedestalperLS3
TH1F * h_sumADCAmplLS1copy2
TH1F * h_sumCutADCAmplperLS5
TH1F * h_runbadrate0_depth2_HB
int lsdep_cut1_peak_HEdepth2_
double sumEstimator5[nsub][ndepth][neta][nphi]
TH1F * h_sumADCAmplperLSdepth6HEu
TH2F * h_recSignalEnergy2_HE3
TH1F * h_runbadrate_depth1_HF
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
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
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_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
TH2F * h_recSignalEnergy0_HB1
TH1F * h_runbadrateC_depth2_HE
TH2F * h_recSignalEnergy1_HE4
TH1F * h_sumADCAmplperLS8_M2
TH1F * h_sumamplitudechannel_HF
TH2F * h_mapDepth2ADCAmpl12SiPM_HE
int lsdep_cut3_max_HFdepth1_
TH2F * h_mapDepth1pedestalw_HB
TH2F * h_mapDepth2RADDAM_HE
TH2F * h_mapDepth6pedestal_HE
TH2F * h_mapDepth3TSmeanA225_HF
TH2F * h_mapDepth4ADCAmpl12_HB
TH2F * h_2D0sumPedestalLS7
TH2F * h_mapDepth1Amplitude_HF
TH1F * h_sumADCAmplLS1copy5
TH2F * h_mapDepth2Ped1_HB
TH2F * h_2D0sumADCAmplLSdepth3HFu
const HcalTopology * topo
TH2F * h_amplitudechannel2_HF4
TH1F * h_runbadrateC_depth2_HF
double forallestimators_amplitude_bigger_
TH1F * h_sum0AmplitudeperLS8
TH2F * h_mapDepth5Ampl047_HE
TH1F * h_sum0AmplitudeperLS3
TH1F * h_ADCAmpl345Zoom1_HB
int lsdep_cut3_max_HBdepth2_
TH1F * h_sumADCAmplperLS7
double lsdep_estimator5_HFdepth2_
TH1F * h_runbadrate_depth2_HF
TH1F * h_sumTSmeanAperLS1
TH2F * h_mapDepth3Amplitude_HB
TH1F * h_energyhitNoise_HB
TH1F * h_sumADCAmplperLS1_M1
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)
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 fiberIndex() const
get the fiber index. For VME 1-8 (which of eight fibers carried by a spigot), for uTCA fibers are zer...
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
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_
TH2F * h_recSignalEnergy1_HB1
TH2F * h_recNoiseEnergy1_HF1
Abs< T >::type abs(const T &t)
TH1F * h_Amplitude_notCapIdErrors_HB2
TH1F * h_energyhitSignal_HE
TH1F * h_runbadrateC_depth1_HE
TH1F * h_orbitNumvsamplitude0_HF
double sumEstimator1[nsub][ndepth][neta][nphi]
TH2F * h_2DsumAmplitudeLS8
TH2F * h_mapDepth3Amplitude225_HF
TH2F * h_mapDepth3TSmeanA225_HE
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
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
TH2F * h_mapDepth3Ped2_HE
TH2F * h_mapDepth1Ped3_HF
TH1F * h_sumADCAmplLS1copy4
TH2F * h_mapDepth7pedestal_HE
TH1F * h_sumCutADCAmplperLS8
TH1F * h_lsnumber_per_eachLS
TH2F * h_mapADCCalib047_HF
double amplitudechannel0[nsub][ndepth][neta][nphi]
TH2F * h_recNoiseEnergy1_HF2
TH1F * h_Amplitude_forCapIdErrors_HO4
int flagestimatornormalization_
TH1F * h_orbitNumvsamplitude0_HO
TH2F * h_mapRatioCalib_HE
TH2F * h_mapADCAmpl12firstpeak_HE
TH2F * h_mapADCAmplfirstpeak0_HE
TH1F * h_sum0TSmaxAperLS6
TH1F * h_pedestalwaver4_HE
double recSignalEnergy0[nsub][ndepth][neta][nphi]
TH2F * h_recNoiseEnergy2_HE1
def gen(fragment, howMuch)
Production test section ####.
TH2F * h_recNoiseEnergy1_HF3
TH2F * h_mapDepth2TSmaxA225_HB
TH1F * h_sumAmplitudeperLS3
TH2F * h_mapDepth4TSmaxA225_HB
TH2F * h_mapDepth6TSmeanA_HE
TH1F * h_sumCutADCAmplperLS6u
TH2F * h_map0AiForLS_bad_HO
double recSignalEnergy1[nsub][ndepth][neta][nphi]
TH1F * h_sumADCAmplperLS1_M2
TH1F * h_sumTSmeanAperLS1_LSselected
double lsdep_estimator1_HBdepth4_
TH1F * h_sumADCAmplperLS3_P1
TH1F * h_sumADCAmplLS1copy3
TH2F * h_mapDepth2Pedw2_HB
TH1F * h_sum0ADCAmplperLSdepth7HEu
TH1F * h_runbadrate_depth3_HE
TH2F * h_mapAiForLS_good_HB
TH1F * h_runbadrate0_depth2_HE
TH2F * h_mapDepth5Ampl_HE
double averNOSIGNALoccupancy_HO
TH1F * h_sigLayer1RADDAM_HE
TH2F * h_mapADCAmpl12secondpeak0_HE
TH1F * h_sumCutADCAmplperLSdepth4HBu
TH1F * h_sumPedestalperLS7
TH2F * h_mapDepth5Error_HE
TH1F * h_shape_bad_channels_HF
int flagtodefinebadchannel_
double recNoiseEnergy0[nsub][ndepth][neta][nphi]
TH1F * h_shape0_Ahigh_HB3
TH2F * h_mapDepth4TS012_HO
TH1F * h_sum0ADCAmplperLS7
TH2F * h_mapDepth3ADCAmpl225_HE
TH1F * h_sumCutAmplperLS1
TH2F * h_mapDepth4ADCAmpl12_HO
TH2F * h_mapDepth1AmplE34_HE
TH2F * h_mapAiForLS_good_HO
TH1F * h_bcnbadrate0_depth2_HF
TH2F * h_mapDepth4AmplE34_HF
TH1F * h_sum0PedestalperLS5
TH2F * h_2D0sumAmplitudeLS5
TH2F * h_mapDepth2TSmaxA225_HF
TH2F * h_mapDepth3Ampl047_HB
TH1F * h_sumCutTSmeanAperLS2
TH2F * h_maprphinorm0_HE1
TH1F * h_sumamplitude_depth2_HB1
TH2F * h_mapADCAmplsecondpeak_HE
TH2F * h_mapDepth3pedestalw_HF
TH2F * h_mapDepth2TSmeanA_HE
TH1F * h_sumamplitude_depth1_HE1
double lsdep_estimator4_HEdepth2_
TH2F * h_amplitudechannel2_HB1
TH2F * h_mapDepth1ADCAmpl_HB
TH1F * h_bcnbadrate0_depth2_HB
TH2F * h_2DsumADCAmplEtaPhiLs10
TH2F * h_recSignalEnergy2_HB3
TH1F * h_runnbadchannelsC_depth3_HE
const_iterator end() const
TH2F * h_mapDepth6pedestalw_HE
TH2F * h_mapDepth3TSmeanA_HB
const HcalFrontEndId getHcalFrontEndId(const DetId &)
double sumEstimator0[nsub][ndepth][neta][nphi]
TH2F * h_2DsumAmplitudeLS2
TH2F * h_amplitudechannel1_HE2
TH2F * h_mapDepth3Ped1_HE
TH2F * h_map0AiForLS_good_HO
int lsdep_cut1_peak_HOdepth4_
TH2F * h_2D0sumPedestalLS5
double lsdep_estimator5_HBdepth1_
TH2F * h_mapADCAmpl12firstpeak0_HE
TH2F * h_mapDepth1Error_HE
double sumEstimator2[nsub][ndepth][neta][nphi]
TH2F * h_mapTSmaxCalib047_HF
TH1F * h_sum0TSmeanAperLS4
TH2F * h_mapDepth3pedestal_HF
TH2F * h_mapDepth3ADCAmpl_HB
TH2F * h_recSignalEnergy2_HF1
TH1F * h_sumADCAmplEtaPhiLs_bbbc
TH1F * h_pedestalwaver4_HF
TH2F * h_recNoiseEnergy2_HE3
TH2F * h_mapADCCalib047_HB
TH2F * h_mapDepth1Amplitude_HE
TH1F * h_sum0TSmaxAperLS8
TH2F * h_mapDepth1ADCAmpl225Copy_HB
TH2F * h_mapDepth2AmplE34_HB
double lsdep_estimator3_HEdepth1_
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
TH2F * h_mapDepth2ADCAmpl225Copy_HE
TH2F * h_mapenophinorm_HE4
TH2F * h_mapDepth1Pedw3_HF
TH2F * h_mapDepth2TSmeanA_HB
TH2F * h_mapTSmeanCalib_HB
int flagIterativeMethodCalibrationGroupReco_
TH2F * h_mapDepth7AmplE34_HE
TH2F * h_amplitudechannel1_HB3
void beginRun(const edm::Run &r, const edm::EventSetup &iSetup) override
TH2F * h_mapWidthCalib_HF
TH2F * h_mapenophinorm_HE6
TH1F * h_nbadchannels_depth1_HE
TH2F * h_mapDepth3TSmeanA_HF
TH1F * h_sumPedestalperLS2
TH2F * h_mapTSmeanCalib047_HE
TH1F * h_sumTSmaxAperLS1_LSselected
TH1F * h_runbadrateC_depth4_HO
TH1F * h_sumamplitude_depth4_HO1
TH1F * h_runnbadchannels_depth1_HB
TH2F * h_mapDepth3Ampl047_HF
TH1F * h_runnbadchannelsC_depth1_HB
TH1F * h_sum0TSmeanAperLS2
TH2F * h_mapenophinorm_HE7
TH1F * h_bcnbadrate0_depth1_HE
TH1F * h_averNOSIGNALsumamplitude_HF
TH2F * h_mapDepth4Ped2_HO
double sumEstimator6[nsub][ndepth][neta][nphi]
TH2F * h_2DsumADCAmplLS6_LSselected
double lsdep_estimator1_HEdepth6_
TH1F * h_sumCutTSmaxAperLS6
TH2F * h_mapDepth7Error_HE
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
constexpr int dccid() const
get the (Hcal local) DCC id for VME, crate number for uTCA
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
TH1F * h_bcnnbadchannels_depth2_HB
TH2F * h_mapDepth1Pedw2_HB
TH1F * h_sumCutADCAmplperLS1
TH1F * h_sumCutAmplperLS2
TH2F * h_mapDepth6TSmaxA225_HE
TH2F * h_mapDepth4pedestalw_HO
TH1F * h_Eta16Depth3RADDAM_HED2
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_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
constexpr int htrTopBottom() const
get the htr top/bottom (1=top/0=bottom), valid for VME
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
constexpr int spigot() const
get the spigot (input number on DCC, AMC card number for uTCA)
TH1F * h_sumamplitude_depth1_HE0
TH2F * h_mapDepth1pedestal_HF
TH2F * h_mapDepth5pedestal_HE
TH2F * h_mapRatioCalib_HF
void fillDigiErrorsHO(HODigiCollection::const_iterator &digiItr)
TH1F * h_A_Depth1RADDAM_HE
TH1F * h_averNOSIGNALoccupancy_HF
TH1F * h_sum0ErrorBperLS7
TH1F * h_sumamplitude_depth1_HF
TH1F * h_runbadrate0_depth1_HE
TH2F * h_mapDepth3Ampl_HE
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
void fillDigiAmplitudeHFQIE10(QIE10DataFrame qie10df)
TH1F * h_eventamplitude_HF
TH1F * h_sum0ADCAmplperLSdepth6HEu
TH2F * h_mapDepth3pedestal_HE
TH1F * h_numberofhitsHOtest
TH1F * h_sumamplitude_depth2_HF1
TH2F * h_2D0sumTSmeanALS1
TH1F * h_totalAmplitudeHFperEvent
TH1F * h_orbitNumvsamplitude_HO
TH2F * h_mapDepth3ADCAmpl_HE
TH1F * h_totalAmplitudeHBperEvent
TH2F * h_2DsumAmplitudeLS7
TH2F * h_mapGetRMSOverNormalizedSignal0_HB
TH2F * h_mapDepth1ADCAmpl225Copy_HF
TH2F * h_recSignalEnergy2_HB1
TH2F * h_mapenophinorm2_HE3
TH2F * h_recSignalEnergy2_HE1
double lsdep_estimator5_HEdepth3_
TH2F * h_recSignalEnergy0_HB4
double const adc2fC_QIE11_shunt1[256]
TH1F * h_runbadrate0_depth1_HB
TH1F * h_sum0AmplitudeperLS7
TH2F * h_mapDepth4ADCAmpl_HF
TH2F * h_mapDepth3linADCAmpl12_HE
TH2F * h_recSignalEnergy0_HE5
TH1F * h_shape_good_channels_HB
TH2F * h_map0AiForLS_good_HE
TH2F * h_mapDepth1Ped1_HB
TH2F * h_recSignalEnergy0_HE7
TH2F * h_2DsumPedestalLS6
TH1F * h_nevents_per_eachRealLS
double lsdep_estimator1_HEdepth7_
TH1F * h_shape0_Ahigh_HB1
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_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
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
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
TH2F * h_recNoiseEnergy2_HE2
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_sumCutADCAmplperLS7
TH2F * h_2D0sumADCAmplLS4
TH1F * h_sumCutAmplitudeperLS4
constexpr int depth() const
get the tower depth
bool studyCapIDErrorsHist_
bool studyPedestalCorrelations_