355 1.89, 5.07, 8.25, 11.43, 14.61, 17.78, 20.96, 24.14, 27.32, 30.50, 33.68,
356 36.86, 40.04, 43.22, 46.40, 49.58, 54.35, 60.71, 67.07, 73.43, 79.79, 86.15,
357 92.51, 98.87, 105.2, 111.6, 117.9, 124.3, 130.7, 137.0, 143.4, 149.7, 156.1,
358 162.5, 168.8, 175.2, 184.7, 197.4, 210.2, 222.9, 235.6, 248.3, 261.0, 273.7,
359 286.5, 299.2, 311.9, 324.6, 337.3, 350.1, 362.8, 375.5, 388.2, 400.9, 413.6,
360 426.4, 439.1, 458.2, 483.6, 509.0, 534.5, 559.9, 585.3, 610.8, 558.9, 584.2,
361 609.5, 634.7, 660.0, 685.3, 710.6, 735.9, 761.2, 786.5, 811.8, 837.1, 862.4,
362 887.7, 913.0, 938.3, 976.2, 1026.8, 1077.4, 1128.0, 1178.6, 1229.1, 1279.7, 1330.3,
363 1380.9, 1431.5, 1482.1, 1532.7, 1583.3, 1633.8, 1684.4, 1735.0, 1785.6, 1836.2, 1886.8,
364 1937.4, 2013.2, 2114.4, 2215.6, 2316.8, 2417.9, 2519.1, 2620.3, 2721.5, 2822.6, 2923.8,
365 3025.0, 3126.2, 3227.3, 3328.5, 3429.7, 3530.9, 3632.0, 3733.2, 3834.4, 3935.5, 4036.7,
366 4188.5, 4390.8, 4593.2, 4795.5, 4997.9, 5200.2, 5402.6, 5057.5, 5262.3, 5467.1, 5671.8,
367 5876.6, 6081.4, 6286.2, 6491.0, 6695.8, 6900.6, 7105.3, 7310.1, 7514.9, 7719.7, 7924.5,
368 8129.3, 8436.4, 8846.0, 9255.6, 9665.1, 10074.7, 10484.3, 10893.9, 11303.4, 11713.0, 12122.6,
369 12532.1, 12941.7, 13351.3, 13760.8, 14170.4, 14580.0, 14989.5, 15399.1, 15808.7, 16218.2, 16832.6,
370 17651.7, 18470.9, 19290.0, 20109.2, 20928.3, 21747.4, 22566.6, 23385.7, 24204.8, 25024.0, 25843.1,
371 26662.3, 27481.4, 28300.5, 29119.7, 29938.8, 30757.9, 31577.1, 32396.2, 33215.4, 34444.1, 36082.3,
372 37720.6, 39358.9, 40997.2, 42635.4, 44273.7, 40908.7, 42542.6, 44176.5, 45810.4, 47444.3, 49078.3,
373 50712.2, 52346.1, 53980.0, 55613.9, 57247.8, 58881.8, 60515.7, 62149.6, 63783.5, 65417.4, 67868.3,
374 71136.1, 74404.0, 77671.8, 80939.7, 84207.5, 87475.3, 90743.2, 94011.0, 97278.8, 100546.7, 103814.5,
375 107082.3, 110350.2, 113618.0, 116885.8, 120153.7, 123421.5, 126689.3, 129957.2, 134858.9, 141394.6, 147930.3,
376 154465.9, 161001.6, 167537.3, 174072.9, 180608.6, 187144.3, 193679.9, 200215.6, 206751.3, 213287.0, 219822.6,
377 226358.3, 232894.0, 239429.6, 245965.3, 252501.0, 259036.6, 265572.3, 275375.8, 288447.2, 301518.5, 314589.8,
378 327661.2, 340732.5, 353803.8};
381 9.56, 28.24, 46.91, 65.59, 84.27, 102.9, 121.6, 140.3, 159.0, 177.7,
382 196.3, 215.0, 233.7, 252.4, 271.0, 289.7, 317.7, 355.1, 392.4, 429.8,
383 467.1, 504.5, 541.9, 579.2, 616.6, 653.9, 691.3, 728.6, 766.0, 803.3,
384 840.7, 878.0, 915.4, 952.8, 990.1, 1027.5, 1083.5, 1158.2, 1232.9, 1307.6,
385 1382.3, 1457.0, 1531.7, 1606.4, 1681.2, 1755.9, 1830.6, 1905.3, 1980.0, 2054.7,
386 2129.4, 2204.1, 2278.8, 2353.5, 2428.2, 2502.9, 2577.7, 2689.7, 2839.1, 2988.6,
387 3138.0, 3287.4, 3436.8, 3586.2, 3263.0, 3411.3, 3559.6, 3707.9, 3856.2, 4004.5,
388 4152.9, 4301.2, 4449.5, 4597.8, 4746.1, 4894.4, 5042.7, 5191.0, 5339.4, 5487.7,
389 5710.1, 6006.8, 6303.4, 6600.0, 6896.6, 7193.3, 7489.9, 7786.5, 8083.1, 8379.8,
390 8676.4, 8973.0, 9269.6, 9566.3, 9862.9, 10159.5, 10456.2, 10752.8, 11049.4, 11346.0,
391 11791.0, 12384.2, 12977.5, 13570.7, 14164.0, 14757.2, 15350.5, 15943.7, 16537.0, 17130.2,
392 17723.5, 18316.7, 18910.0, 19503.2, 20096.5, 20689.7, 21283.0, 21876.2, 22469.5, 23062.8,
393 23656.0, 24545.9, 25732.4, 26918.9, 28105.4, 29291.9, 30478.4, 31664.9, 29399.4, 30590.1,
394 31780.9, 32971.7, 34162.4, 35353.2, 36544.0, 37734.7, 38925.5, 40116.3, 41307.0, 42497.8,
395 43688.5, 44879.3, 46070.1, 47260.8, 49047.0, 51428.5, 53810.1, 56191.6, 58573.1, 60954.6,
396 63336.2, 65717.7, 68099.2, 70480.8, 72862.3, 75243.8, 77625.4, 80006.9, 82388.4, 84769.9,
397 87151.5, 89533.0, 91914.5, 94296.1, 97868.4, 102631.4, 107394.5, 112157.5, 116920.6, 121683.7,
398 126446.7, 131209.8, 135972.8, 140735.9, 145499.0, 150262.0, 155025.1, 159788.2, 164551.2, 169314.3,
399 174077.3, 178840.4, 183603.5, 188366.5, 193129.6, 200274.2, 209800.3, 219326.4, 228852.5, 238378.7,
400 247904.8, 257430.9, 237822.7, 247326.7, 256830.7, 266334.8, 275838.8, 285342.9, 294846.9, 304351.0,
401 313855.0, 323359.1, 332863.1, 342367.1, 351871.2, 361375.2, 370879.3, 380383.3, 394639.4, 413647.5,
402 432655.6, 451663.6, 470671.7, 489679.8, 508687.9, 527696.0, 546704.1, 565712.2, 584720.3, 603728.3,
403 622736.4, 641744.5, 660752.6, 679760.7, 698768.8, 717776.9, 736785.0, 755793.0, 784305.2, 822321.3,
404 860337.5, 898353.7, 936369.9, 974386.0, 1012402.2, 1050418.4, 1088434.6, 1126450.7, 1164466.9, 1202483.1,
405 1240499.3, 1278515.4, 1316531.6, 1354547.8, 1392564.0, 1430580.1, 1468596.3, 1506612.5, 1544628.7, 1601652.9,
406 1677685.3, 1753717.6, 1829750.0, 1905782.3, 1981814.7, 2057847.0};
409 -0.5, 0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5, 10.5, 11.5,
410 12.5, 13.5, 15., 17., 19., 21., 23., 25., 27., 29.5, 32.5, 35.5, 38.5,
411 42., 46., 50., 54.5, 59.5, 64.5, 59.5, 64.5, 69.5, 74.5, 79.5, 84.5, 89.5,
412 94.5, 99.5, 104.5, 109.5, 114.5, 119.5, 124.5, 129.5, 137., 147., 157., 167., 177.,
413 187., 197., 209.5, 224.5, 239.5, 254.5, 272., 292., 312., 334.5, 359.5, 384.5, 359.5,
414 384.5, 409.5, 434.5, 459.5, 484.5, 509.5, 534.5, 559.5, 584.5, 609.5, 634.5, 659.5, 684.5,
415 709.5, 747., 797., 847., 897., 947., 997., 1047., 1109.5, 1184.5, 1259.5, 1334.5, 1422.,
416 1522., 1622., 1734.5, 1859.5, 1984.5, 1859.5, 1984.5, 2109.5, 2234.5, 2359.5, 2484.5, 2609.5, 2734.5,
417 2859.5, 2984.5, 3109.5, 3234.5, 3359.5, 3484.5, 3609.5, 3797., 4047., 4297., 4547., 4797., 5047.,
418 5297., 5609.5, 5984.5, 6359.5, 6734.5, 7172., 7672., 8172., 8734.5, 9359.5, 9984.5};
442 double phi12(
double phi1,
double en1,
double phi2,
double en2);
443 double dR(
double eta1,
double phi1,
double eta2,
double phi2);
444 double dPhiWsign(
double phi1,
double phi2);
2376 std::cout <<
"============== endJob ===================================" << std::endl;
2378 std::cout <<
" --------------------------------------- " << std::endl;
2381 std::cout <<
" --------------------------------------------- " << std::endl;
2387 std::cout <<
"===== possible max number of events in ntuple(each 50th recorded) = " <<
nevent50 << std::endl;
2391 std::cout <<
"===== full number of events*HBHEdigis = " <<
nnnnnn << std::endl;
2396 std::cout <<
"===== Start writing user histograms =====" << std::endl;
2631 for (
int k1 = 0; k1 <
ndepth; k1++) {
2632 for (
int k2 = 0; k2 <
neta; k2++) {
2637 for (
int k3 = 0; k3 <
nphi; k3++) {
2696 int outabortgap = 1;
2706 std::cout <<
" --------------------------------------- " << std::endl;
2710 std::cout <<
" --------------------------------------------- " << std::endl;
2716 std::cout <<
" ---------***********************------------- " << std::endl;
2766 unsigned long int pcountall1 = 0;
2767 unsigned long int pcountall3 = 0;
2768 unsigned long int pcountall6 = 0;
2769 unsigned long int pcountall8 = 0;
2774 unsigned long int mcountall1 = 0;
2775 unsigned long int mcountall3 = 0;
2776 unsigned long int mcountall6 = 0;
2777 unsigned long int mcountall8 = 0;
2786 int pnnmin1 = 999999999;
2787 int pnnmin3 = 999999999;
2788 int pnnmin6 = 999999999;
2789 int pnnmin8 = 999999999;
2794 int mnnmin1 = 999999999;
2795 int mnnmin3 = 999999999;
2796 int mnnmin6 = 999999999;
2797 int mnnmin8 = 999999999;
2799 for (
int k1 = 0; k1 <
ndepth; k1++) {
2800 for (
int k3 = 0; k3 <
nphi; k3++) {
2801 for (
int k2 = 0; k2 <
neta; k2++) {
2910 int kkkk2 = (k2 - 1) / 4;
2915 int kkkk3 = (k3) / 4 + 1;
2919 ietaphi = ((kkkk2)-1) *
znphi + (kkkk3);
2987 if (bbbc / bbb1 > 25.) {
3002 if (bbbc / bbb1 > 25.) {
3087 if (bbbc / bbb1 > 15. && k3 % 2 == 0) {
3101 if (bbbc / bbb1 > 15. && k3 % 2 == 0) {
3241 if (bbbc / bbb1 > 20.) {
3256 if (bbbc / bbb1 > 20.) {
3348 if (bbbc / bbb1 > 80.) {
3363 if (bbbc / bbb1 > 80.) {
3918 if (
k0 == 0 && k1 == 0) {
3919 if (pcountmin1 > 0) {
3920 if (pcountmin1 < pnnmin1)
3921 pnnmin1 = pcountmin1;
3925 if (mcountmin1 > 0) {
3926 if (mcountmin1 < mnnmin1)
3927 mnnmin1 = mcountmin1;
3933 if (
k0 == 1 && k1 == 0) {
3934 if (pcountmin3 > 0) {
3935 if (pcountmin3 < pnnmin3)
3936 pnnmin3 = pcountmin3;
3940 if (mcountmin3 > 0) {
3941 if (mcountmin3 < mnnmin3)
3942 mnnmin3 = mcountmin3;
3948 if (
k0 == 2 && k1 == 3) {
3949 if (pcountmin8 > 0) {
3950 if (pcountmin8 < pnnmin8)
3951 pnnmin8 = pcountmin8;
3955 if (mcountmin8 > 0) {
3956 if (mcountmin8 < mnnmin8)
3957 mnnmin8 = mcountmin8;
3963 if (
k0 == 3 && k1 == 0) {
3964 if (pcountmin6 > 0) {
3965 if (pcountmin6 < pnnmin6)
3966 pnnmin6 = pcountmin6;
3970 if (mcountmin6 > 0) {
3971 if (mcountmin6 < mnnmin6)
3972 mnnmin6 = mcountmin6;
3987 float patiooccupancy1 = 0.;
3988 if (pcountall1 != 0)
3989 patiooccupancy1 = (
float)pnnmin1 * mnnbins1 / pcountall1;
3991 float matiooccupancy1 = 0.;
3992 if (mcountall1 != 0)
3993 matiooccupancy1 = (
float)mnnmin1 * mnnbins1 / mcountall1;
3996 float patiooccupancy3 = 0.;
3997 if (pcountall3 != 0)
3998 patiooccupancy3 = (
float)pnnmin3 * mnnbins3 / pcountall3;
4000 float matiooccupancy3 = 0.;
4001 if (mcountall3 != 0)
4002 matiooccupancy3 = (
float)mnnmin3 * mnnbins3 / mcountall3;
4005 float patiooccupancy6 = 0.;
4006 if (pcountall6 != 0)
4007 patiooccupancy6 = (
float)pnnmin6 * mnnbins6 / pcountall6;
4009 float matiooccupancy6 = 0.;
4010 if (mcountall6 != 0)
4011 matiooccupancy6 = (
float)mnnmin6 * mnnbins6 / mcountall6;
4014 float patiooccupancy8 = 0.;
4015 if (pcountall8 != 0)
4016 patiooccupancy8 = (
float)pnnmin8 * mnnbins8 / pcountall8;
4018 float matiooccupancy8 = 0.;
4019 if (mcountall8 != 0)
4020 matiooccupancy8 = (
float)mnnmin8 * mnnbins8 / mcountall8;
4024 for (
int k1 = 0; k1 <
ndepth; k1++) {
4025 for (
int k2 = 0; k2 <
neta; k2++) {
4026 for (
int k3 = 0; k3 <
nphi; k3++) {
4126 for (
int k1 = 0; k1 <
ndepth; k1++) {
4127 for (
int k2 = 0; k2 <
neta; k2++) {
4128 for (
int k3 = 0; k3 <
nphi; k3++) {
4150 for (
int k1 = 0; k1 <
ndepth; k1++) {
4151 for (
int k2 = 0; k2 <
neta; k2++) {
4152 for (
int k3 = 0; k3 <
nphi; k3++) {
4184 bool gotHFDigis =
true;
4188 if (!(
hf.isValid())) {
4192 std::cout <<
" ****************************** =========================== No HFDigiCollection found " 4197 eta = digi->id().ieta();
4198 phi = digi->id().iphi();
4199 depth = digi->id().depth();
4216 for (
int i = 0;
i <
nTS;
i++) {
4235 bool gotQIE10Digis =
true;
4237 gotQIE10Digis =
false;
4240 gotQIE10Digis =
false;
4242 if (!gotQIE10Digis) {
4243 std::cout <<
" No QIE10DigiCollection collection is found " << std::endl;
4246 double totalAmplitudeHF = 0.;
4247 for (
unsigned int j = 0;
j < qie10dc.
size();
j++) {
4271 double amplitudefullTSs = 0.;
4272 double nnnnnnTS = 0.;
4273 for (
int i = 0;
i <
nTS; ++
i) {
4276 int adc = qie10df[
i].adc();
4308 bool gotHBHEDigis =
true;
4310 gotHBHEDigis =
false;
4311 if (!(
hbhe.isValid()))
4312 gotHBHEDigis =
false;
4313 if (!gotHBHEDigis) {
4314 std::cout <<
" No HBHEDigiCollection collection is found " << std::endl;
4317 double totalAmplitudeHB = 0.;
4318 double totalAmplitudeHE = 0.;
4319 double nnnnnnTSHB = 0.;
4320 double nnnnnnTSHE = 0.;
4323 eta = digi->id().ieta();
4324 phi = digi->id().iphi();
4325 depth = digi->id().depth();
4365 double amplitudefullTSs = 0.;
4368 for (
int i = 0;
i <
nTS;
i++) {
4380 double amplitudefullTSs = 0.;
4383 for (
int i = 0;
i <
nTS;
i++) {
4398 if (totalAmplitudeHB != 0.) {
4403 if (totalAmplitudeHE != 0.) {
4419 bool gotQIE11Digis =
true;
4421 gotQIE11Digis =
false;
4423 gotQIE11Digis =
false;
4424 if (!gotQIE11Digis) {
4425 std::cout <<
" No QIE11DigiCollection collection is found " << std::endl;
4428 double totalAmplitudeHBQIE11 = 0.;
4429 double totalAmplitudeHEQIE11 = 0.;
4430 double nnnnnnTSHBQIE11 = 0.;
4431 double nnnnnnTSHEQIE11 = 0.;
4432 for (
unsigned int j = 0;
j < qie11dc.
size();
j++) {
4441 int sub = hcaldetid.
subdet();
4454 if (sub == 1 &&
depth == 1 && qwert1 == 0) {
4458 if (sub == 1 &&
depth == 2 && qwert2 == 0) {
4462 if (sub == 2 &&
depth == 1 && qwert3 == 0) {
4466 if (sub == 2 &&
depth == 2 && qwert4 == 0) {
4470 if (sub == 2 &&
depth == 3 && qwert5 == 0) {
4482 double amplitudefullTSs1 = 0.;
4483 double amplitudefullTSs6 = 0.;
4485 for (
int i = 0;
i <
nTS; ++
i) {
4486 int adc = qie11df[
i].adc();
4489 amplitudefullTSs1 += charge1;
4490 amplitudefullTSs6 += charge6;
4496 totalAmplitudeHBQIE11 +=
charge;
4503 double amplitudefullTSs1 = 0.;
4504 double amplitudefullTSs6 = 0.;
4506 for (
int i = 0;
i <
nTS;
i++) {
4507 int adc = qie11df[
i].adc();
4510 amplitudefullTSs1 += charge1;
4511 amplitudefullTSs6 += charge6;
4517 totalAmplitudeHEQIE11 +=
charge;
4526 if (totalAmplitudeHBQIE11 != 0.) {
4531 if (totalAmplitudeHEQIE11 != 0.) {
4542 bool gotHODigis =
true;
4545 if (!(
ho.isValid()))
4549 std::cout <<
" No HO collection is found " << std::endl;
4552 double totalAmplitudeHO = 0.;
4554 eta = digi->id().ieta();
4555 phi = digi->id().iphi();
4556 depth = digi->id().depth();
4577 double nnnnnnTS = 0.;
4578 double amplitudefullTSs = 0.;
4580 for (
int i = 0;
i <
nTS;
i++) {
4609 bool gotHBHERecHitsNoise =
true;
4611 gotHBHERecHitsNoise =
4614 gotHBHERecHitsNoise =
false;
4615 if (!gotHBHERecHitsNoise) {
4617 std::cout <<
" No RecHits HBHENoise collection is found " << std::endl;
4622 float icalconst = 1.;
4625 HBHERecHit aHit(hbheItr->id(), hbheItr->eraw() * icalconst, hbheItr->time());
4627 double energyhit = aHit.
energy();
4628 DetId id = (*hbheItr).detid();
4630 int sub = ((hid).rawId() >> 25) & 0x7;
4639 if (hid.
depth() == 1 && sub == 1 && hid.
iphi() == 25) {
4641 std::cout <<
" Noise,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4642 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4644 if (hid.
depth() == 1 && sub == 2 && hid.
iphi() == 25) {
4646 std::cout <<
" Noise,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4647 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4653 int mdepth = hid.
depth();
4664 bool gotHFRecHitsNoise =
true;
4666 gotHFRecHitsNoise =
false;
4668 gotHFRecHitsNoise =
false;
4669 if (!gotHFRecHitsNoise) {
4671 std::cout <<
" No RecHits HFNoise collection is found " << std::endl;
4675 float icalconst = 1.;
4678 HFRecHit aHit(hfItr->id(), hfItr->energy() * icalconst, hfItr->time());
4679 double energyhit = aHit.
energy();
4680 DetId id = (*hfItr).detid();
4682 int sub = ((hid).rawId() >> 25) & 0x7;
4688 if (hid.
iphi() == 25) {
4690 std::cout <<
"HF Noise,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4691 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4697 int mdepth = hid.
depth();
4711 bool gotHBHERecHitsSignal =
true;
4713 gotHBHERecHitsSignal =
4716 gotHBHERecHitsSignal =
false;
4717 if (!gotHBHERecHitsSignal) {
4719 std::cout <<
" No RecHits HBHESignal collection is found " << std::endl;
4724 float icalconst = 1.;
4727 HBHERecHit aHit(hbheItr->id(), hbheItr->eraw() * icalconst, hbheItr->time());
4729 double energyhit = aHit.
energy();
4730 DetId id = (*hbheItr).detid();
4732 int sub = ((hid).rawId() >> 25) & 0x7;
4740 if (hid.
depth() == 1 && sub == 1 && hid.
iphi() == 25) {
4742 std::cout <<
"HBHE Signal,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4743 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4745 if (hid.
depth() == 1 && sub == 2 && hid.
iphi() == 25) {
4747 std::cout <<
"HBHE Signal,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4748 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4754 int mdepth = hid.
depth();
4765 bool gotHFRecHitsSignal =
true;
4767 gotHFRecHitsSignal =
false;
4769 gotHFRecHitsSignal =
false;
4770 if (!gotHFRecHitsSignal) {
4772 std::cout <<
" No RecHits HFSignal collection is found " << std::endl;
4776 float icalconst = 1.;
4779 HFRecHit aHit(hfItr->id(), hfItr->energy() * icalconst, hfItr->time());
4780 double energyhit = aHit.
energy();
4781 DetId id = (*hfItr).detid();
4783 int sub = ((hid).rawId() >> 25) & 0x7;
4789 if (hid.
iphi() == 25) {
4791 std::cout <<
"HF Signal,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4792 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4798 int mdepth = hid.
depth();
4828 for (
int k1 = 0; k1 <
ndepth; k1++) {
4830 for (
int k2 = 0; k2 <
neta; k2++) {
4832 int k2plot = k2 - 41;
4835 double sumoverphi = 0;
4836 int nsumoverphi = 0;
4837 for (
int k3 = 0; k3 <
nphi; k3++) {
4842 std::cout <<
"==== nsumoverphi = " << nsumoverphi <<
" sumoverphi = " << sumoverphi
4843 <<
" k1 = " << k1 <<
" k2 = " << k2 <<
" kkk = " << kkk <<
" k3 = " << k3 << std::endl;
4847 for (
int k3 = 0; k3 <
nphi; k3++) {
4848 if (nsumoverphi != 0) {
4851 std::cout <<
"nsumoverphi= " << nsumoverphi <<
" sumoverphi= " << sumoverphi <<
" k1= " << k1
4852 <<
" k2= " << k2 <<
" kkk= " << kkk <<
" k3= " << k3
4853 <<
" maprphinorm= " <<
maprphinorm[
k0][k1][k2][k3] << std::endl;
4928 for (
int k1 = 0; k1 <
ndepth; k1++) {
4930 for (
int k2 = 0; k2 <
neta; k2++) {
4932 int k2plot = k2 - 41;
4934 for (
int k3 = 0; k3 <
nphi; k3++) {
4962 for (
int k1 = 0; k1 <
ndepth; k1++) {
4964 for (
int k2 = 0; k2 <
neta; k2++) {
4966 int k2plot = k2 - 41;
4968 for (
int k3 = 0; k3 <
nphi; k3++) {
5011 for (
int k1 = 0; k1 <
ndepth; k1++) {
5013 for (
int k2 = 0; k2 <
neta; k2++) {
5015 int k2plot = k2 - 41;
5017 for (
int k3 = 0; k3 <
nphi; k3++) {
5060 for (
int k1 = 0; k1 <
ndepth; k1++) {
5062 for (
int k2 = 0; k2 <
neta; k2++) {
5064 int k2plot = k2 - 41;
5066 for (
int k3 = 0; k3 <
nphi; k3++) {
5106 for (
int k1 = 0; k1 <
ndepth; k1++) {
5108 for (
int k2 = 0; k2 <
neta; k2++) {
5110 int k2plot = k2 - 41;
5112 for (
int k3 = 0; k3 <
nphi; k3++) {
5176 for (
int k1 = 0; k1 <
ndepth; k1++) {
5178 for (
int k2 = 0; k2 <
neta; k2++) {
5180 int k2plot = k2 - 41;
5182 for (
int k3 = 0; k3 <
nphi; k3++) {
5230 for (
int k1 = 0; k1 <
ndepth; k1++) {
5231 for (
int k2 = 0; k2 <
neta; k2++) {
5232 for (
int k3 = 0; k3 <
nphi; k3++) {
5244 for (
int k1 = 0; k1 <
ndepth; k1++) {
5245 for (
int k2 = 0; k2 <
neta; k2++) {
5248 if (k1 == 2 && k2 == 25) {
5252 }
else if (k1 == 2 && k2 == 56) {
5267 for (
int k1 = 0; k1 <
ndepth; k1++) {
5268 for (
int k2 = 0; k2 <
neta; k2++) {
5269 if (k1 == 2 && (k2 == 25 || k2 == 56)) {
5272 int k2plot = k2 - 41;
5282 int LLLLLL111111 = 0;
5283 if ((k1 == 0 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 30))
5286 int LLLLLL222222 = 0;
5287 if ((k1 == 1 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 27) ||
5288 (k1 == 2 && fabs(kkk - 41) > 26 && fabs(kkk - 41) < 29))
5291 if (LLLLLL111111 == 1) {
5295 if (LLLLLL222222 == 1) {
5308 for (
int k1 = 0; k1 <
ndepth; k1++) {
5309 for (
int k2 = 0; k2 <
neta; k2++) {
5310 if (k1 == 2 && (k2 == 25 || k2 == 56)) {
5312 int k2plot = k2 - 41;
5319 for (
int k3 = 0; k3 <
nphi; k3++) {
5322 int LLLLLL111111 = 0;
5323 if ((k1 == 0 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 30))
5325 int LLLLLL222222 = 0;
5326 if ((k1 == 1 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 27) ||
5327 (k1 == 2 && fabs(kkk - 41) > 26 && fabs(kkk - 41) < 29))
5329 if (LLLLLL111111 == 1) {
5333 if (LLLLLL222222 == 1) {
5344 for (
int k1 = 0; k1 <
ndepth; k1++) {
5345 for (
int k2 = 0; k2 <
neta; k2++) {
5348 for (
int k3 = 0; k3 <
nphi; k3++) {
5366 int sumofchannels = 0;
5367 double sumamplitudesubdet = 0.;
5368 int sumofchannels0 = 0;
5369 double sumamplitudesubdet0 = 0.;
5370 for (
int k2 = 0; k2 <
neta; k2++) {
5371 for (
int k3 = 0; k3 <
nphi; k3++) {
5383 if (sumamplitudechannel_HB > 80.) {
5384 sumamplitudesubdet += sumamplitudechannel_HB;
5387 if (sumamplitudechannel_HB > 0.) {
5388 sumamplitudesubdet0 += sumamplitudechannel_HB;
5396 double sumamplitudechannel_HE =
5399 if (sumamplitudechannel_HE > 200.) {
5400 sumamplitudesubdet += sumamplitudechannel_HE;
5403 if (sumamplitudechannel_HE > 0.) {
5404 sumamplitudesubdet0 += sumamplitudechannel_HE;
5414 if (sumamplitudechannel_HO > 1200.) {
5415 sumamplitudesubdet += sumamplitudechannel_HO;
5418 if (sumamplitudechannel_HO > 0.) {
5419 sumamplitudesubdet0 += sumamplitudechannel_HO;
5428 if (sumamplitudechannel_HF > 600.) {
5429 sumamplitudesubdet += sumamplitudechannel_HF;
5432 if (sumamplitudechannel_HF > 0.) {
5433 sumamplitudesubdet0 += sumamplitudechannel_HF;
5446 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM1)
5447 maxxSUM1 = sumamplitudesubdet + sumamplitudesubdet0;
5448 if ((sumofchannels + sumofchannels0) >
maxxOCCUP1)
5454 if ((sumamplitudesubdet + sumamplitudesubdet0) > 60000) {
5455 for (
int k2 = 0; k2 <
neta; k2++) {
5456 for (
int k3 = 0; k3 <
nphi; k3++) {
5475 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM2)
5476 maxxSUM2 = sumamplitudesubdet + sumamplitudesubdet0;
5477 if ((sumofchannels + sumofchannels0) >
maxxOCCUP2)
5487 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM3)
5488 maxxSUM3 = sumamplitudesubdet + sumamplitudesubdet0;
5489 if ((sumofchannels + sumofchannels0) >
maxxOCCUP3)
5499 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM4)
5500 maxxSUM4 = sumamplitudesubdet + sumamplitudesubdet0;
5501 if ((sumofchannels + sumofchannels0) >
maxxOCCUP4)
5522 for (
int k1 = 0; k1 <
ndepth; k1++) {
5524 int nbadchannels = 0;
5525 for (
int k2 = 0; k2 <
neta; k2++) {
5526 for (
int k3 = 0; k3 <
nphi; k3++) {
5539 if (nbadchannels != 0)
5552 if (nbadchannels != 0)
5568 if (nbadchannels != 0)
5581 if (nbadchannels != 0)
5594 if (nbadchannels != 0)
5610 if (nbadchannels != 0)
5626 if (nbadchannels != 0)
5639 if (nbadchannels != 0)
5659 bool gotCALIBDigis =
true;
5661 gotCALIBDigis =
false;
5663 if (!(
calib.isValid())) {
5664 gotCALIBDigis =
false;
5666 if (!gotCALIBDigis) {
5669 int cal_det = digi->id().hcalSubdet();
5670 int cal_phi = digi->id().iphi();
5671 int cal_eta = digi->id().ieta();
5672 int cal_cbox = digi->id().cboxChannel();
5676 if (cal_det > 0 && cal_det < 5 && cal_cbox == 0) {
5677 int iphi = cal_phi - 1;
5682 double max_signal = -100.;
5683 int ts_with_max_signal = -100;
5688 for (
int i = 0;
i <
nTS;
i++) {
5689 double ampldefault =
adc2fC[digi->sample(
i).adc() & 0xff];
5690 if (max_signal < ampldefault) {
5691 max_signal = ampldefault;
5692 ts_with_max_signal =
i;
5697 timew += (
i + 1) * ampldefault;
5702 double aveamplitude = -100.;
5705 double aveamplitude1 = aveamplitude - 1;
5710 double ampldefault =
adc2fC[digi->sample(
ii).adc() & 0xff];
5711 double aaaaaa = (
ii + 1) - aveamplitude;
5712 double aaaaaa2 = aaaaaa * aaaaaa;
5713 rmsamp += (aaaaaa2 * ampldefault);
5715 double rmsamplitude = -100.;
5723 if (ts_with_max_signal > -1 && ts_with_max_signal <
nTS)
5725 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 <
nTS)
5726 calib2[cal_det - 1][
ieta + 41][
iphi] +=
adc2fC[digi->sample(ts_with_max_signal + 1).adc() & 0xff];
5727 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 <
nTS)
5728 calib2[cal_det - 1][
ieta + 41][
iphi] +=
adc2fC[digi->sample(ts_with_max_signal + 2).adc() & 0xff];
5729 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 <
nTS)
5730 calib2[cal_det - 1][
ieta + 41][
iphi] +=
adc2fC[digi->sample(ts_with_max_signal - 1).adc() & 0xff];
5731 if (ts_with_max_signal - 2 > -1 && ts_with_max_signal - 2 <
nTS)
5732 calib2[cal_det - 1][
ieta + 41][
iphi] +=
adc2fC[digi->sample(ts_with_max_signal - 2).adc() & 0xff];
5734 bool anycapid =
true;
5737 int error1 = 0, error2 = 0, error3 = 0;
5738 int lastcapid = 0, capid = 0;
5739 for (
int ii = 0;
ii < (*digi).size();
ii++) {
5740 capid = (*digi)[
ii].capid();
5741 bool er = (*digi)[
ii].er();
5742 bool dv = (*digi)[
ii].dv();
5743 if (
ii != 0 && ((lastcapid + 1) % 4) != capid)
5757 if (error1 != 0 || error2 != 0 || error3 != 0)
5772 for (
int k1 = 0; k1 <
nsub; k1++) {
5773 for (
int k2 = 0; k2 <
neta; k2++) {
5774 for (
int k3 = 0; k3 <
nphi; k3++) {
5775 int k2plot = k2 - 41;
5779 if (
signal[k1][k2][k3] > 0.) {
5805 if (k20 > 0 || k20 == 0)
5810 int kk2 = 0, kk3 = 0;
5811 if (k1 == 0 || k1 == 1) {
5816 if (k30 == 71 || k30 ==
nphi || k30 == 1 || k30 == 2)
5819 kk3 = ((k30 - 3) / 4) * 4 + 3;
5820 }
else if (k1 == 2) {
5821 if (
abs(k20) <= 4) {
5823 if (k30 == 71 || k30 ==
nphi || k30 == 1 || k30 == 2 || k30 == 3 || k30 == 4)
5826 kk3 = ((k30 - 5) / 6) * 6 + 5;
5828 if (
abs(k20) > 4 &&
abs(k20) <= 10)
5830 if (
abs(k20) > 10 &&
abs(k20) <= 15)
5834 if (k30 == 71 || k30 ==
nphi || (k30 >= 1 && k30 <= 10))
5837 kk3 = ((k30 - 11) / 12) * 12 + 11;
5839 }
else if (k1 == 3) {
5844 if (k30 >= 1 && k30 <= 18)
5846 if (k30 >= 19 && k30 <= 36)
5848 if (k30 >= 37 && k30 <= 54)
5850 if (k30 >= 55 && k30 <=
nphi)
5854 int kkk2 = kk2 + 41;
5861 double GetRMSOverNormalizedSignal = -1.;
5862 if (
signal[k1][k2][k3] > 0. &&
calib0[k1][kkk2][kkk3] > 0.) {
5863 GetRMSOverNormalizedSignal =
signal[k1][k2][k3] /
calib0[k1][kkk2][kkk3];
5883 if (
signal[k1][k2][k3] > 0.) {
5886 if (
calib0[k1][kkk2][kkk3] > 0.)
5890 if (
calib0[k1][kkk2][kkk3] > 0.)
5893 float calibtsmax =
calibt[k1][kkk2][kkk3];
5895 float calibtsmean =
caliba[k1][kkk2][kkk3];
5897 float calibwidth =
calibw[k1][kkk2][kkk3];
5899 float calibcap = -100.;
5916 if (calibtsmax > -0.5) {
5923 if (calibtsmean > -0.5) {
5930 if (calibwidth > -0.5) {
5940 h_map_HB->Fill(
double(k2plot),
double(k3), 1.);
5956 if (calibtsmax > -0.5) {
5963 if (calibtsmean > -0.5) {
5970 if (calibwidth > -0.5) {
5980 h_map_HE->Fill(
double(k2plot),
double(k3), 1.);
5996 if (calibtsmax > -0.5) {
6003 if (calibtsmean > -0.5) {
6010 if (calibwidth > -0.5) {
6020 h_map_HO->Fill(
double(k2plot),
double(k3), 1.);
6036 if (calibtsmax > -0.5) {
6043 if (calibtsmean > -0.5) {
6050 if (calibwidth > -0.5) {
6060 h_map_HF->Fill(
double(k2plot),
double(k3), 1.);
6092 std::cout <<
"======================== beignJob START +++++++++++++++++++++++++++" << std::endl;
6330 fs_->
make<TH1F>(
"h_ndepthesperamplitudebins_HE",
" ", 10, 0., 10.);
6340 "h_mapADCAmplsecondpeak0_HE",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
6697 fs_->
make<TH1F>(
"h_AamplitudewithPedSubtr_RADDAM_HEzoom0",
" ", 100, min80, 4000.);
6699 fs_->
make<TH1F>(
"h_AamplitudewithPedSubtr_RADDAM_HEzoom1",
" ", 100, min80, 1000.);
6713 fs_->
make<TH1F>(
"h_sumphiEta16Depth3RADDAM_HED2ALL",
" ", 100, min90, 70. * max90);
6861 float zaR = 10000000.;
6862 float zaR2 = 50000000.;
6896 fs_->
make<TH1F>(
"h_sumADCAmplEtaPhiLs_lscounterM1orbitNum",
" ", 600, 1., 601.);
6904 float bac2 = bac + 1.;
7864 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HB",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7866 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HB",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7868 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HE",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7870 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HE",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7872 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HF",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7874 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HF",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7876 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HO",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7878 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HO",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7927 float spls0 = 10000;
7937 float spls1 = 100000;
7973 std::cout <<
"======================== booking DONE +++++++++++++++++++++++++++" << std::endl;
7982 std::cout <<
"======================== beignJob finish +++++++++++++++++++++++++++" << std::endl;
7995 int mdepth = cell.
depth();
7996 int iphi = cell.iphi() - 1;
7997 int ieta = cell.ieta();
8000 int sub = cell.subdet();
8002 int errorGeneral = 0;
8011 bool anycapid =
true;
8017 int lastcapid = 0, capid = 0;
8018 int ERRORfiber = -10;
8019 int ERRORfiberChan = -10;
8020 int ERRORfiberAndChan = -10;
8021 int repetedcapid = 0;
8023 TSsize = digiItr->size();
8026 for (
int ii = 0;
ii < TSsize;
ii++) {
8027 capid = (*digiItr)[
ii].capid();
8028 bool er = (*digiItr)[
ii].er();
8029 bool dv = (*digiItr)[
ii].dv();
8030 int fiber = (*digiItr)[
ii].fiber();
8031 int fiberChan = (*digiItr)[
ii].fiberChan();
8032 int fiberAndChan = (*digiItr)[
ii].fiberAndChan();
8033 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8036 ERRORfiberChan = fiberChan;
8037 ERRORfiberAndChan = fiberAndChan;
8038 if (capid != lastcapid) {
8052 ERRORfiberChan = fiberChan;
8053 ERRORfiberAndChan = fiberAndChan;
8058 ERRORfiberChan = fiberChan;
8059 ERRORfiberAndChan = fiberAndChan;
8065 if (firstcapid == 0 && !anycapid)
8067 if (firstcapid == 1 && !anycapid)
8069 if (firstcapid == 2 && !anycapid)
8071 if (firstcapid == 3 && !anycapid)
8080 if (!anycapid && anyer)
8082 if (!anycapid && !anydv)
8084 if (!anycapid && anyer && !anydv)
8086 if (anyer && !anydv)
8092 for (
int ii = 0;
ii < TSsize;
ii++) {
8093 double ampldefault =
adc2fC[digiItr->sample(
ii).adc()];
8094 ampl += ampldefault;
8111 if (error1 != 0 || error2 != 0 || error3 != 0) {
8138 if (error1 != 0 || error2 != 0 || error3 != 0) {
8169 int mdepth = hcaldetid.
depth();
8170 int sub = hcaldetid.
subdet();
8174 bool anycapid =
true;
8180 int lastcapid = 0, capid = 0;
8181 int repetedcapid = 0;
8186 capid = qie11df[
ii].capid();
8187 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8189 if (capid != lastcapid) {
8205 int adc = qie11df[
ii].adc();
8211 ampl += ampldefault;
8270 int mdepth = cell.
depth();
8271 int iphi = cell.iphi() - 1;
8272 int ieta = cell.ieta();
8275 int sub = cell.subdet();
8277 std::cout <<
" HF DIGI ?????????????? ERROR mdepth = " << mdepth << std::endl;
8279 int errorGeneral = 0;
8288 bool anycapid =
true;
8294 int lastcapid = 0, capid = 0;
8295 int ERRORfiber = -10;
8296 int ERRORfiberChan = -10;
8297 int ERRORfiberAndChan = -10;
8298 int repetedcapid = 0;
8301 TSsize = digiItr->size();
8303 for (
int ii = 0;
ii < TSsize;
ii++) {
8304 capid = (*digiItr)[
ii].capid();
8305 bool er = (*digiItr)[
ii].er();
8306 bool dv = (*digiItr)[
ii].dv();
8307 int fiber = (*digiItr)[
ii].fiber();
8308 int fiberChan = (*digiItr)[
ii].fiberChan();
8309 int fiberAndChan = (*digiItr)[
ii].fiberAndChan();
8310 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8313 ERRORfiberChan = fiberChan;
8314 ERRORfiberAndChan = fiberAndChan;
8315 if (capid != lastcapid) {
8327 ERRORfiberChan = fiberChan;
8328 ERRORfiberAndChan = fiberAndChan;
8333 ERRORfiberChan = fiberChan;
8334 ERRORfiberAndChan = fiberAndChan;
8338 if (firstcapid == 0 && !anycapid)
8340 if (firstcapid == 1 && !anycapid)
8342 if (firstcapid == 2 && !anycapid)
8344 if (firstcapid == 3 && !anycapid)
8352 if (!anycapid && anyer)
8354 if (!anycapid && !anydv)
8356 if (!anycapid && anyer && !anydv)
8358 if (anyer && !anydv)
8362 for (
int ii = 0;
ii < TSsize;
ii++) {
8363 double ampldefault =
adc2fC[digiItr->sample(
ii).adc()];
8364 ampl += ampldefault;
8378 if (error1 != 0 || error2 != 0 || error3 != 0) {
8408 int mdepth = hcaldetid.
depth();
8409 int sub = hcaldetid.
subdet();
8413 bool anycapid =
true;
8419 int lastcapid = 0, capid = 0;
8420 int repetedcapid = 0;
8425 capid = qie10df[
ii].capid();
8426 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8428 if (capid != lastcapid) {
8444 int adc = qie10df[
ii].adc();
8446 ampl += ampldefault;
8477 int mdepth = cell.
depth();
8478 int iphi = cell.iphi() - 1;
8479 int ieta = cell.ieta();
8482 int sub = cell.subdet();
8483 int errorGeneral = 0;
8492 bool anycapid =
true;
8498 int lastcapid = 0, capid = 0;
8499 int ERRORfiber = -10;
8500 int ERRORfiberChan = -10;
8501 int ERRORfiberAndChan = -10;
8502 int repetedcapid = 0;
8503 for (
int ii = 0;
ii < (*digiItr).size();
ii++) {
8504 capid = (*digiItr)[
ii].capid();
8505 bool er = (*digiItr)[
ii].er();
8506 bool dv = (*digiItr)[
ii].dv();
8507 int fiber = (*digiItr)[
ii].fiber();
8508 int fiberChan = (*digiItr)[
ii].fiberChan();
8509 int fiberAndChan = (*digiItr)[
ii].fiberAndChan();
8510 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8513 ERRORfiberChan = fiberChan;
8514 ERRORfiberAndChan = fiberAndChan;
8515 if (capid != lastcapid) {
8529 ERRORfiberChan = fiberChan;
8530 ERRORfiberAndChan = fiberAndChan;
8535 ERRORfiberChan = fiberChan;
8536 ERRORfiberAndChan = fiberAndChan;
8542 if (firstcapid == 0 && !anycapid)
8544 if (firstcapid == 1 && !anycapid)
8546 if (firstcapid == 2 && !anycapid)
8548 if (firstcapid == 3 && !anycapid)
8557 if (!anycapid && anyer)
8559 if (!anycapid && !anydv)
8561 if (!anycapid && anyer && !anydv)
8563 if (anyer && !anydv)
8568 for (
int ii = 0;
ii < (*digiItr).size();
ii++) {
8569 double ampldefault =
adc2fC[digiItr->sample(
ii).adc()];
8570 ampl += ampldefault;
8586 if (error1 != 0 || error2 != 0 || error3 != 0) {
8595 std::cout <<
" mdepth HO = " << mdepth << std::endl;
8608 double toolwithPedSubtr[100];
8609 double lintoolwithoutPedSubtr[100];
8611 int mdepth = cell.
depth();
8612 int iphi = cell.iphi() - 1;
8613 int ieta0 = cell.ieta();
8617 int sub = cell.subdet();
8629 coder.
adc2fC(*digiItr, toolOriginal);
8630 double pedestalaver9 = 0.;
8631 double pedestalaver4 = 0.;
8632 double pedestal0 = 0.;
8633 double pedestal1 = 0.;
8634 double pedestal2 = 0.;
8635 double pedestal3 = 0.;
8636 double pedestalwaver9 = 0.;
8637 double pedestalwaver4 = 0.;
8638 double pedestalw0 = 0.;
8639 double pedestalw1 = 0.;
8640 double pedestalw2 = 0.;
8641 double pedestalw3 = 0.;
8642 double difpedestal0 = 0.;
8643 double difpedestal1 = 0.;
8644 double difpedestal2 = 0.;
8645 double difpedestal3 = 0.;
8647 double amplitudewithPedSubtr1 = 0.;
8648 double amplitudewithPedSubtr2 = 0.;
8649 double amplitudewithPedSubtr3 = 0.;
8650 double amplitudewithPedSubtr4 = 0.;
8652 double absamplitude = 0.;
8653 double amplitude345 = 0.;
8655 double linamplitudewithoutPedSubtr = 0.;
8658 double max_signal = -100.;
8659 int ts_with_max_signal = -100;
8665 double errorBtype = 0.;
8668 TSsize = digiItr->size();
8669 if ((*digiItr).size() != 10)
8672 for (
int ii = 0;
ii < TSsize;
ii++) {
8674 double ampldefaultwithPedSubtr = 0.;
8675 double linampldefaultwithoutPedSubtr = 0.;
8676 double ampldefault = 0.;
8677 double ampldefault0 = 0.;
8678 double ampldefault1 = 0.;
8679 double ampldefault2 = 0.;
8680 ampldefault0 =
adc2fC[digiItr->sample(
ii).adc()];
8682 ampldefault1 = toolOriginal[
ii];
8683 ampldefault2 = (*digiItr)[
ii].adc();
8685 ampldefault = ampldefault0;
8688 ampldefault = ampldefault1;
8691 ampldefault = ampldefault2;
8693 ampldefaultwithPedSubtr = ampldefault0;
8694 linampldefaultwithoutPedSubtr = ampldefault2;
8696 int capid = ((*digiItr)[
ii]).capid();
8698 double pedestalINI = pedestal00->
getValue(capid);
8700 double pedestalw = pedw->
getSigma(capid, capid);
8701 ampldefaultwithPedSubtr -=
pedestal;
8705 tool[
ii] = ampldefault;
8706 toolwithPedSubtr[
ii] = ampldefaultwithPedSubtr;
8707 lintoolwithoutPedSubtr[
ii] = linampldefaultwithoutPedSubtr;
8710 pedestalwaver9 += pedestalw * pedestalw;
8712 if (capid == 0 &&
c0 == 0) {
8717 pedestalw0 = pedestalw;
8718 pedestalwaver4 += pedestalw * pedestalw;
8719 difpedestal0 =
pedestal - pedestalINI;
8722 if (capid == 1 &&
c1 == 0) {
8727 pedestalw1 = pedestalw;
8728 pedestalwaver4 += pedestalw * pedestalw;
8729 difpedestal1 =
pedestal - pedestalINI;
8731 if (capid == 2 && c2 == 0) {
8736 pedestalw2 = pedestalw;
8737 pedestalwaver4 += pedestalw * pedestalw;
8738 difpedestal2 =
pedestal - pedestalINI;
8740 if (capid == 3 && c3 == 0) {
8745 pedestalw3 = pedestalw;
8746 pedestalwaver4 += pedestalw * pedestalw;
8747 difpedestal3 =
pedestal - pedestalINI;
8750 if (max_signal < ampldefault) {
8751 max_signal = ampldefault;
8752 ts_with_max_signal =
ii;
8763 if (
ii > 0 &&
ii < 9)
8764 amplitudewithPedSubtr1 += ampldefaultwithPedSubtr;
8766 if (
ii > 1 &&
ii < 8)
8767 amplitudewithPedSubtr2 += ampldefaultwithPedSubtr;
8769 if (
ii > 2 &&
ii < 7)
8770 amplitudewithPedSubtr3 += ampldefaultwithPedSubtr;
8772 if (
ii > 2 &&
ii < 6)
8773 amplitudewithPedSubtr4 += ampldefaultwithPedSubtr;
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 < 10)
8798 ampl = tool[ts_with_max_signal];
8799 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < 10)
8800 ampl += tool[ts_with_max_signal + 2];
8801 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < 10)
8802 ampl += tool[ts_with_max_signal + 1];
8803 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < 10)
8804 ampl += tool[ts_with_max_signal - 1];
8807 if (ts_with_max_signal > -1 && ts_with_max_signal < 10)
8808 linamplitudewithoutPedSubtr = lintoolwithoutPedSubtr[ts_with_max_signal];
8809 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < 10)
8810 linamplitudewithoutPedSubtr += lintoolwithoutPedSubtr[ts_with_max_signal + 2];
8811 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < 10)
8812 linamplitudewithoutPedSubtr += lintoolwithoutPedSubtr[ts_with_max_signal + 1];
8813 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < 10)
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_) {
9126 if (aveamplitude1 < TSmeanHBMin_ || aveamplitude1 >
TSmeanHBMax_) {
9147 if (ts_with_max_signal < TSpeakHBMin_ || ts_with_max_signal >
TSpeakHBMax_) {
9168 if (rmsamplitude < rmsHBMin_ || rmsamplitude >
rmsHBMax_) {
9227 double mypedestal = pedestal0;
9228 double mypedestalw = pedestalw0;
9299 for (
int ii = 0;
ii < TSsize;
ii++) {
9321 if (amplitude < ADCAmplHEMin_ || amplitude >
ADCAmplHEMax_) {
9378 if (aveamplitude1 < TSmeanHEMin_ || aveamplitude1 >
TSmeanHEMax_) {
9402 if (ts_with_max_signal < TSpeakHEMin_ || ts_with_max_signal >
TSpeakHEMax_) {
9426 if (rmsamplitude < rmsHEMin_ || rmsamplitude >
rmsHEMax_) {
9486 double amplitudewithPedSubtr = 0.;
9489 if (ts_with_max_signal > -1 && ts_with_max_signal < 10)
9490 amplitudewithPedSubtr = toolwithPedSubtr[ts_with_max_signal];
9491 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < 10)
9492 amplitudewithPedSubtr += toolwithPedSubtr[ts_with_max_signal + 2];
9493 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < 10)
9494 amplitudewithPedSubtr += toolwithPedSubtr[ts_with_max_signal + 1];
9495 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < 10)
9496 amplitudewithPedSubtr += toolwithPedSubtr[ts_with_max_signal - 1];
9502 if (amplitudewithPedSubtr > 50.) {
9525 int LLLLLL111111 = 0;
9526 if ((mdepth == 1 && fabs(ieta0) > 16 && fabs(ieta0) < 30))
9529 int LLLLLL222222 = 0;
9530 if ((mdepth == 2 && fabs(ieta0) > 16 && fabs(ieta0) < 27) ||
9531 (mdepth == 3 && fabs(ieta0) > 26 && fabs(ieta0) < 29))
9534 if (LLLLLL111111 == 1) {
9539 if (LLLLLL222222 == 1) {
9546 if (mdepth == 3 && fabs(ieta0) == 16) {
9575 int mdepth = hcaldetid.
depth();
9576 int sub = hcaldetid.
subdet();
9579 if (mdepth == 0 || sub > 4)
9604 coder.
adc2fC(qie11df, toolOriginal);
9605 double pedestalaver9 = 0.;
9606 double pedestalaver4 = 0.;
9607 double pedestal0 = 0.;
9608 double pedestal1 = 0.;
9609 double pedestal2 = 0.;
9610 double pedestal3 = 0.;
9611 double pedestalwaver9 = 0.;
9612 double pedestalwaver4 = 0.;
9613 double pedestalw0 = 0.;
9614 double pedestalw1 = 0.;
9615 double pedestalw2 = 0.;
9616 double pedestalw3 = 0.;
9617 double difpedestal0 = 0.;
9618 double difpedestal1 = 0.;
9619 double difpedestal2 = 0.;
9620 double difpedestal3 = 0.;
9623 double amplitude0 = 0.;
9624 double absamplitude = 0.;
9625 double tocampl = 0.;
9627 double amplitude345 = 0.;
9629 double ampl3ts = 0.;
9630 double amplmaxts = 0.;
9633 double max_signal = -100.;
9634 int ts_with_max_signal = -100;
9640 double errorBtype = 0.;
9651 for (
int ii = 0;
ii < TSsize;
ii++) {
9652 double ampldefault = 0.;
9653 double tocdefault = 0.;
9654 double ampldefault0 = 0.;
9655 double ampldefault1 = 0.;
9656 double ampldefault2 = 0.;
9662 ampldefault1 = toolOriginal[
ii];
9663 ampldefault2 = qie11df[
ii].adc();
9666 ampldefault = ampldefault0;
9669 ampldefault = ampldefault1;
9672 ampldefault = ampldefault2;
9674 tocdefault = ampldefault;
9676 int capid = (qie11df[
ii]).capid();
9678 double pedestalw = pedw->
getSigma(capid, capid);
9679 double pedestalINI = pedestal00->
getValue(capid);
9683 tool[
ii] = ampldefault;
9685 pedestalwaver9 += pedestalw * pedestalw;
9687 if (capid == 0 &&
c0 == 0) {
9692 pedestalw0 = pedestalw;
9693 pedestalwaver4 += pedestalw * pedestalw;
9694 difpedestal0 =
pedestal - pedestalINI;
9697 if (capid == 1 &&
c1 == 0) {
9702 pedestalw1 = pedestalw;
9703 pedestalwaver4 += pedestalw * pedestalw;
9704 difpedestal1 =
pedestal - pedestalINI;
9706 if (capid == 2 && c2 == 0) {
9711 pedestalw2 = pedestalw;
9712 pedestalwaver4 += pedestalw * pedestalw;
9713 difpedestal2 =
pedestal - pedestalINI;
9715 if (capid == 3 && c3 == 0) {
9720 pedestalw3 = pedestalw;
9721 pedestalwaver4 += pedestalw * pedestalw;
9722 difpedestal3 =
pedestal - pedestalINI;
9725 if (max_signal < ampldefault) {
9726 max_signal = ampldefault;
9727 ts_with_max_signal =
ii;
9730 absamplitude +=
abs(ampldefault);
9731 tocampl += tocdefault;
9733 if (
ii == 1 ||
ii == 2 ||
ii == 3 ||
ii == 4 ||
ii == 5 ||
ii == 6 ||
ii == 7 ||
ii == 8)
9734 amplitude345 += ampldefault;
9739 timew += (
ii + 1) *
abs(ampldefault);
9740 timeww += (
ii + 1) * ampldefault;
9742 if (
ii == 2 && ampldefault > 0.)
9748 pedestalaver9 /= TSsize;
9749 pedestalaver4 /= c4;
9750 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
9751 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
9754 if (ts_with_max_signal > -1 && ts_with_max_signal < 10) {
9755 ampl = tool[ts_with_max_signal];
9756 ampl3ts = tool[ts_with_max_signal];
9757 amplmaxts = tool[ts_with_max_signal];
9759 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < 10) {
9760 ampl += tool[ts_with_max_signal - 1];
9761 ampl3ts += tool[ts_with_max_signal - 1];
9763 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < 10) {
9764 ampl += tool[ts_with_max_signal + 1];
9765 ampl3ts += tool[ts_with_max_signal + 1];
9767 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < 10) {
9768 ampl += tool[ts_with_max_signal + 2];
9774 double xb = ampl / 40.;
9775 double xc = amplitude345 / 40.;
9776 double xd = ampl3ts / 40.;
9777 double xe = amplmaxts / 40.;
9778 double txa = tocampl / 40.;
9780 if (((
ieta == -16 ||
ieta == 15) && mdepth == 4) ||
9781 ((
ieta == -17 ||
ieta == 16) && (mdepth == 2 || mdepth == 3)) ||
9782 ((
ieta == -18 ||
ieta == 17) && mdepth == 5)) {
9783 double c0 = 1.000000;
9784 double b1 = 2.59096e-05;
9785 double a2 = 4.60721e-11;
9786 double corrforxa =
a2 * xa * xa +
b1 * xa +
c0;
9787 double corrforxb =
a2 * xb * xb +
b1 * xb +
c0;
9788 double corrforxc =
a2 * xc * xc +
b1 * xc +
c0;
9789 double corrforxd =
a2 * xd * xd +
b1 * xd +
c0;
9790 double corrforxe =
a2 * xe * xe +
b1 * xe +
c0;
9791 double corrfortxa =
a2 * txa * txa +
b1 * txa +
c0;
9796 amplitude345 *= corrforxc;
9797 ampl3ts *= corrforxd;
9798 amplmaxts *= corrforxe;
9799 tocampl *= corrfortxa;
9803 double c0 = 1.000000;
9804 double b1 = 2.71238e-05;
9805 double a2 = 1.32877e-10;
9806 double corrforxa =
a2 * xa * xa +
b1 * xa +
c0;
9807 double corrforxb =
a2 * xb * xb +
b1 * xb +
c0;
9808 double corrforxc =
a2 * xc * xc +
b1 * xc +
c0;
9809 double corrforxd =
a2 * xd * xd +
b1 * xd +
c0;
9810 double corrforxe =
a2 * xe * xe +
b1 * xe +
c0;
9811 double corrfortxa =
a2 * txa * txa +
b1 * txa +
c0;
9816 amplitude345 *= corrforxc;
9817 ampl3ts *= corrforxd;
9818 amplmaxts *= corrforxe;
9819 tocampl *= corrfortxa;
9825 if (sub == 1 || sub == 2) {
9836 if (ratio < 0. || ratio > 1.02)
9838 double aveamplitude = 0.;
9839 double aveamplitudew = 0.;
9840 if (absamplitude > 0 && timew > 0)
9841 aveamplitude = timew / absamplitude;
9842 if (amplitude0 > 0 && timeww > 0)
9843 aveamplitudew = timeww / amplitude0;
9847 bool anycapid =
true;
9850 for (
int ii = 0;
ii < TSsize;
ii++) {
9851 double aaaaaa = (
ii + 1) - aveamplitudew;
9852 double aaaaaa2 = aaaaaa * aaaaaa;
9853 double ampldefault = tool[
ii];
9854 rmsamp += (aaaaaa2 * ampldefault);
9855 capid = (qie11df[
ii]).capid();
9856 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
9864 double rmsamplitude = 0.;
9865 if ((amplitude0 > 0 && rmsamp > 0) || (amplitude0 < 0 && rmsamp < 0))
9866 rmsamplitude =
sqrt(rmsamp / amplitude0);
9867 double aveamplitude1 = aveamplitude - 1;
9872 if (sub == 1 && mdepth == 1)
9874 if (sub == 1 && mdepth == 2)
9876 if (sub == 2 && mdepth == 1)
9878 if (sub == 2 && mdepth == 2)
9880 if (sub == 2 && mdepth == 3)
9884 if (sub == 1 && mdepth == 1)
9886 if (sub == 1 && mdepth == 2)
9888 if (sub == 2 && mdepth == 1)
9890 if (sub == 2 && mdepth == 2)
9892 if (sub == 2 && mdepth == 3)
9896 for (
int ii = 0;
ii < TSsize;
ii++) {
9898 double ampldefault = tool[
ii];
9901 if (amplitude0 > 120) {
9944 for (
int ii = 0;
ii < TSsize;
ii++) {
9945 double ampldefault = tool[
ii];
9952 for (
int ii = 0;
ii < TSsize;
ii++) {
9953 double ampldefault = tool[
ii];
9967 for (
int ii = 0;
ii < TSsize;
ii++) {
9968 double ampldefault = tool[
ii];
9975 for (
int ii = 0;
ii < TSsize;
ii++) {
9976 double ampldefault = tool[
ii];
9992 double mypedestal = pedestal0;
9993 double mypedestalw = pedestalw0;
10064 for (
int ii = 0;
ii < TSsize;
ii++) {
10097 if (amplitude < ADCAmplHBMin_ || amplitude >
ADCAmplHBMax_) {
10151 if (aveamplitude1 < TSmeanHBMin_ || aveamplitude1 >
TSmeanHBMax_) {
10179 h_TSmaxA_HB->Fill(
float(ts_with_max_signal), 1.);
10180 if (ts_with_max_signal < TSpeakHBMin_ || ts_with_max_signal >
TSpeakHBMax_) {
10209 if (rmsamplitude < rmsHBMin_ || rmsamplitude >
rmsHBMax_) {
10238 if (ratio < ratioHBMin_ || ratio >
ratioHBMax_) {
10287 if (flagTS2 == 1) {
10306 double mypedestal = pedestal0;
10307 double mypedestalw = pedestalw0;
10394 for (
int ii = 0;
ii < TSsize;
ii++) {
10412 if (amplitude < ADCAmplHEMin_ || amplitude >
ADCAmplHEMax_) {
10460 if (ampl3ts > 70 && ampl3ts < 110) {
10463 }
else if (ampl3ts > 110 && ampl3ts < 150) {
10467 if (ampl > 87 && ampl < 127) {
10470 }
else if (ampl > 127 && ampl < 167) {
10529 if (aveamplitude1 < TSmeanHEMin_ || aveamplitude1 >
TSmeanHEMax_) {
10568 h_TSmaxA_HE->Fill(
float(ts_with_max_signal), 1.);
10569 if (ts_with_max_signal < TSpeakHEMin_ || ts_with_max_signal >
TSpeakHEMax_) {
10609 if (rmsamplitude < rmsHEMin_ || rmsamplitude >
rmsHEMax_) {
10650 if (ratio < ratioHEMin_ || ratio >
ratioHEMax_) {
10720 if (flagTS2 == 1) {
10745 int mdepth = cell.
depth();
10746 int iphi = cell.iphi() - 1;
10747 int ieta = cell.ieta();
10750 int sub = cell.subdet();
10760 coder.
adc2fC(*digiItr, toolOriginal);
10761 double pedestalaver9 = 0.;
10762 double pedestalaver4 = 0.;
10763 double pedestal0 = 0.;
10764 double pedestal1 = 0.;
10765 double pedestal2 = 0.;
10766 double pedestal3 = 0.;
10767 double pedestalwaver9 = 0.;
10768 double pedestalwaver4 = 0.;
10769 double pedestalw0 = 0.;
10770 double pedestalw1 = 0.;
10771 double pedestalw2 = 0.;
10772 double pedestalw3 = 0.;
10774 double absamplitude = 0.;
10777 double timeww = 0.;
10778 double max_signal = -100.;
10779 int ts_with_max_signal = -100;
10785 double errorBtype = 0.;
10787 if ((*digiItr).size() != TSsize)
10789 TSsize = digiItr->size();
10790 for (
int ii = 0;
ii < TSsize;
ii++) {
10792 double ampldefault = 0.;
10793 double ampldefault0 = 0.;
10794 double ampldefault1 = 0.;
10795 double ampldefault2 = 0.;
10796 ampldefault0 =
adc2fC[digiItr->sample(
ii).adc()];
10798 ampldefault1 = toolOriginal[
ii];
10799 ampldefault2 = (*digiItr)[
ii].adc();
10801 ampldefault = ampldefault0;
10804 ampldefault = ampldefault1;
10807 ampldefault = ampldefault2;
10810 int capid = ((*digiItr)[
ii]).capid();
10813 double pedestalw = pedw->
getSigma(capid, capid);
10817 tool[
ii] = ampldefault;
10820 pedestalwaver9 += pedestalw * pedestalw;
10822 if (capid == 0 &&
c0 == 0) {
10827 pedestalw0 = pedestalw;
10828 pedestalwaver4 += pedestalw * pedestalw;
10831 if (capid == 1 &&
c1 == 0) {
10836 pedestalw1 = pedestalw;
10837 pedestalwaver4 += pedestalw * pedestalw;
10839 if (capid == 2 && c2 == 0) {
10844 pedestalw2 = pedestalw;
10845 pedestalwaver4 += pedestalw * pedestalw;
10847 if (capid == 3 && c3 == 0) {
10852 pedestalw3 = pedestalw;
10853 pedestalwaver4 += pedestalw * pedestalw;
10856 if (max_signal < ampldefault) {
10857 max_signal = ampldefault;
10858 ts_with_max_signal =
ii;
10861 absamplitude +=
abs(ampldefault);
10866 timew += (
ii + 1) *
abs(ampldefault);
10867 timeww += (
ii + 1) * ampldefault;
10877 pedestalaver9 /= TSsize;
10878 pedestalaver4 /= c4;
10879 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
10880 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
10883 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
10884 ampl = tool[ts_with_max_signal];
10885 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
10886 ampl += tool[ts_with_max_signal + 2];
10887 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
10888 ampl += tool[ts_with_max_signal + 1];
10889 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
10890 ampl += tool[ts_with_max_signal - 1];
10897 if (ratio < 0. || ratio > 1.02)
10900 double aveamplitude = 0.;
10901 double aveamplitudew = 0.;
10902 if (absamplitude > 0 && timew > 0)
10903 aveamplitude = timew / absamplitude;
10907 double rmsamp = 0.;
10910 bool anycapid =
true;
10911 bool anyer =
false;
10915 for (
int ii = 0;
ii < TSsize;
ii++) {
10916 double aaaaaa = (
ii + 1) - aveamplitudew;
10917 double aaaaaa2 = aaaaaa * aaaaaa;
10918 double ampldefault = tool[
ii];
10919 rmsamp += (aaaaaa2 * ampldefault);
10920 capid = ((*digiItr)[
ii]).capid();
10921 bool er = (*digiItr)[
ii].er();
10922 bool dv = (*digiItr)[
ii].dv();
10923 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
10936 if (!anycapid || anyer || !anydv)
10938 double rmsamplitude = 0.;
10941 double aveamplitude1 = aveamplitude - 1;
10943 if (sub == 4 && mdepth == 1)
10945 if (sub == 4 && mdepth == 2)
10949 if (sub == 4 && mdepth == 1)
10951 if (sub == 4 && mdepth == 2)
10963 for (
int ii = 0;
ii < TSsize;
ii++) {
10964 double ampldefault = tool[
ii];
10971 for (
int ii = 0;
ii < TSsize;
ii++) {
10972 double ampldefault = tool[
ii];
10986 double mypedestal = pedestal0;
10987 double mypedestalw = pedestalw0;
11047 for (
int ii = 0;
ii < TSsize;
ii++) {
11064 if (amplitude < ADCAmplHFMin_ || amplitude >
ADCAmplHFMax_) {
11102 if (aveamplitude1 < TSmeanHFMin_ || aveamplitude1 >
TSmeanHFMax_) {
11122 h_TSmaxA_HF->Fill(
float(ts_with_max_signal), 1.);
11123 if (ts_with_max_signal < TSpeakHFMin_ || ts_with_max_signal >
TSpeakHFMax_) {
11143 if (rmsamplitude < rmsHFMin_ || rmsamplitude >
rmsHFMax_) {
11163 if (ratio < ratioHFMin_ || ratio >
ratioHFMax_) {
11211 int mdepth = hcaldetid.
depth();
11212 int sub = hcaldetid.
subdet();
11227 if (mdepth == 0 || sub != 4)
11247 coder.
adc2fC(qie10df, toolOriginal);
11250 double pedestalaver9 = 0.;
11251 double pedestalaver4 = 0.;
11252 double pedestal0 = 0.;
11253 double pedestal1 = 0.;
11254 double pedestal2 = 0.;
11255 double pedestal3 = 0.;
11256 double pedestalwaver9 = 0.;
11257 double pedestalwaver4 = 0.;
11258 double pedestalw0 = 0.;
11259 double pedestalw1 = 0.;
11260 double pedestalw2 = 0.;
11261 double pedestalw3 = 0.;
11263 double absamplitude = 0.;
11266 double timeww = 0.;
11267 double max_signal = -100.;
11268 int ts_with_max_signal = -100;
11274 double errorBtype = 0.;
11281 for (
int ii = 0;
ii < TSsize;
ii++) {
11282 double ampldefault = 0.;
11283 double ampldefault0 = 0.;
11284 double ampldefault1 = 0.;
11285 double ampldefault2 = 0.;
11288 ampldefault1 = toolOriginal[
ii];
11289 ampldefault2 = qie10df[
ii].adc();
11291 ampldefault = ampldefault0;
11294 ampldefault = ampldefault1;
11297 ampldefault = ampldefault2;
11300 int capid = (qie10df[
ii]).capid();
11302 double pedestalw = pedw->
getSigma(capid, capid);
11307 tool[
ii] = ampldefault;
11310 pedestalwaver9 += pedestalw * pedestalw;
11312 if (capid == 0 &&
c0 == 0) {
11317 pedestalw0 = pedestalw;
11318 pedestalwaver4 += pedestalw * pedestalw;
11321 if (capid == 1 &&
c1 == 0) {
11326 pedestalw1 = pedestalw;
11327 pedestalwaver4 += pedestalw * pedestalw;
11329 if (capid == 2 && c2 == 0) {
11334 pedestalw2 = pedestalw;
11335 pedestalwaver4 += pedestalw * pedestalw;
11337 if (capid == 3 && c3 == 0) {
11342 pedestalw3 = pedestalw;
11343 pedestalwaver4 += pedestalw * pedestalw;
11346 if (max_signal < ampldefault) {
11347 max_signal = ampldefault;
11348 ts_with_max_signal =
ii;
11351 absamplitude +=
abs(ampldefault);
11353 timew += (
ii + 1) *
abs(ampldefault);
11354 timeww += (
ii + 1) * ampldefault;
11355 if (
ii == 1 && ampldefault > 0.)
11366 pedestalaver9 /= TSsize;
11367 pedestalaver4 /= c4;
11368 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
11369 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
11372 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
11373 ampl = tool[ts_with_max_signal];
11374 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
11375 ampl += tool[ts_with_max_signal + 2];
11376 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
11377 ampl += tool[ts_with_max_signal + 1];
11378 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
11379 ampl += tool[ts_with_max_signal - 1];
11385 if (ratio < 0. || ratio > 1.02)
11387 double aveamplitude = 0.;
11388 double aveamplitudew = 0.;
11389 if (absamplitude > 0 && timew > 0)
11390 aveamplitude = timew / absamplitude;
11394 double rmsamp = 0.;
11396 bool anycapid =
true;
11399 for (
int ii = 0;
ii < TSsize;
ii++) {
11400 double aaaaaa = (
ii + 1) - aveamplitudew;
11401 double aaaaaa2 = aaaaaa * aaaaaa;
11402 double ampldefault = tool[
ii];
11403 rmsamp += (aaaaaa2 * ampldefault);
11404 capid = (qie10df[
ii]).capid();
11405 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
11413 double rmsamplitude = 0.;
11416 double aveamplitude1 = aveamplitude - 1;
11421 if (sub == 4 && mdepth == 1)
11423 if (sub == 4 && mdepth == 2)
11427 if (sub == 4 && mdepth == 1)
11429 if (sub == 4 && mdepth == 2)
11441 for (
int ii = 0;
ii < TSsize;
ii++) {
11442 double ampldefault = tool[
ii];
11449 for (
int ii = 0;
ii < TSsize;
ii++) {
11450 double ampldefault = tool[
ii];
11465 double mypedestal = pedestal0;
11466 double mypedestalw = pedestalw0;
11534 for (
int ii = 0;
ii < TSsize;
ii++) {
11551 if (amplitude < ADCAmplHFMin_ || amplitude >
ADCAmplHFMax_) {
11603 if (aveamplitude1 < TSmeanHFMin_ || aveamplitude1 >
TSmeanHFMax_) {
11631 h_TSmaxA_HF->Fill(
float(ts_with_max_signal), 1.);
11632 if (ts_with_max_signal < TSpeakHFMin_ || ts_with_max_signal >
TSpeakHFMax_) {
11661 if (rmsamplitude < rmsHFMin_ || rmsamplitude >
rmsHFMax_) {
11690 if (ratio < ratioHFMin_ || ratio >
ratioHFMax_) {
11739 if (flagTS1 == 1) {
11760 int mdepth = cell.
depth();
11761 int iphi = cell.iphi() - 1;
11762 int ieta = cell.ieta();
11765 int sub = cell.subdet();
11775 coder.
adc2fC(*digiItr, toolOriginal);
11776 double pedestalaver9 = 0.;
11777 double pedestalaver4 = 0.;
11778 double pedestal0 = 0.;
11779 double pedestal1 = 0.;
11780 double pedestal2 = 0.;
11781 double pedestal3 = 0.;
11782 double pedestalwaver9 = 0.;
11783 double pedestalwaver4 = 0.;
11784 double pedestalw0 = 0.;
11785 double pedestalw1 = 0.;
11786 double pedestalw2 = 0.;
11787 double pedestalw3 = 0.;
11789 double absamplitude = 0.;
11792 double timeww = 0.;
11793 double max_signal = -100.;
11794 int ts_with_max_signal = -100;
11800 double errorBtype = 0.;
11802 if ((*digiItr).size() != TSsize)
11804 TSsize = digiItr->size();
11806 for (
int ii = 0;
ii < TSsize;
ii++) {
11807 double ampldefault = 0.;
11808 double ampldefault0 = 0.;
11809 double ampldefault1 = 0.;
11810 double ampldefault2 = 0.;
11811 ampldefault0 =
adc2fC[digiItr->sample(
ii).adc()];
11813 ampldefault1 = toolOriginal[
ii];
11814 ampldefault2 = (*digiItr)[
ii].adc();
11816 ampldefault = ampldefault0;
11819 ampldefault = ampldefault1;
11822 ampldefault = ampldefault2;
11824 int capid = ((*digiItr)[
ii]).capid();
11826 double pedestalw = pedw->
getSigma(capid, capid);
11829 tool[
ii] = ampldefault;
11831 pedestalwaver9 += pedestalw * pedestalw;
11832 if (capid == 0 &&
c0 == 0) {
11837 pedestalw0 = pedestalw;
11838 pedestalwaver4 += pedestalw * pedestalw;
11840 if (capid == 1 &&
c1 == 0) {
11845 pedestalw1 = pedestalw;
11846 pedestalwaver4 += pedestalw * pedestalw;
11848 if (capid == 2 && c2 == 0) {
11853 pedestalw2 = pedestalw;
11854 pedestalwaver4 += pedestalw * pedestalw;
11856 if (capid == 3 && c3 == 0) {
11861 pedestalw3 = pedestalw;
11862 pedestalwaver4 += pedestalw * pedestalw;
11865 if (max_signal < ampldefault) {
11866 max_signal = ampldefault;
11867 ts_with_max_signal =
ii;
11870 absamplitude +=
abs(ampldefault);
11874 timew += (
ii + 1) *
abs(ampldefault);
11875 timeww += (
ii + 1) * ampldefault;
11876 if (ii < 3 && ampldefault > 0.)
11881 pedestalaver9 /= TSsize;
11882 pedestalaver4 /= c4;
11883 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
11884 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
11885 if (ts_with_max_signal > -1 && ts_with_max_signal < 10)
11886 ampl = tool[ts_with_max_signal];
11887 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < 10)
11888 ampl += tool[ts_with_max_signal + 2];
11889 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < 10)
11890 ampl += tool[ts_with_max_signal + 1];
11891 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < 10)
11892 ampl += tool[ts_with_max_signal - 1];
11896 if (ratio < 0. || ratio > 1.04)
11898 double aveamplitude = 0.;
11899 double aveamplitudew = 0.;
11900 if (absamplitude > 0 && timew > 0)
11901 aveamplitude = timew / absamplitude;
11904 double rmsamp = 0.;
11906 bool anycapid =
true;
11907 bool anyer =
false;
11911 for (
int ii = 0;
ii < TSsize;
ii++) {
11912 double aaaaaa = (
ii + 1) - aveamplitudew;
11913 double aaaaaa2 = aaaaaa * aaaaaa;
11914 double ampldefault = tool[
ii];
11915 rmsamp += (aaaaaa2 * ampldefault);
11916 capid = ((*digiItr)[
ii]).capid();
11917 bool er = (*digiItr)[
ii].er();
11918 bool dv = (*digiItr)[
ii].dv();
11919 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
11931 if (!anycapid || anyer || !anydv)
11933 double rmsamplitude = 0.;
11936 double aveamplitude1 = aveamplitude - 1;
11938 if (sub == 3 && mdepth == 4)
11942 if (sub == 3 && mdepth == 4)
11953 for (
int ii = 0;
ii < TSsize;
ii++) {
11954 double ampldefault = tool[
ii];
11959 for (
int ii = 0;
ii < TSsize;
ii++) {
11960 double ampldefault = tool[
ii];
11974 double mypedestal = pedestal0;
11975 double mypedestalw = pedestalw0;
12018 for (
int ii = 0;
ii < TSsize;
ii++) {
12034 if (amplitude < ADCAmplHOMin_ || amplitude >
ADCAmplHOMax_) {
12062 if (aveamplitude1 < TSmeanHOMin_ || aveamplitude1 >
TSmeanHOMax_) {
12074 h_TSmaxA_HO->Fill(
float(ts_with_max_signal), 1.);
12075 if (ts_with_max_signal < TSpeakHOMin_ || ts_with_max_signal >
TSpeakHOMax_) {
12089 if (rmsamplitude < rmsHOMin_ || rmsamplitude >
rmsHOMax_) {
12102 if (ratio < ratioHOMin_ || ratio >
ratioHOMax_) {
12119 if (flagTS012 == 1)
12126 if (kdet == 1 || kdet == 2) {
12130 cal_RBX = (kphi + 1) / 4;
12131 cal_RBX = cal_RBX + 18 * (keta + 1) / 2;
12134 cal_RBX = (
int)(kphi / 18) + 1;
12141 cal_RBX = kphi / 12 + 1;
12147 cal_RBX = kphi / 12 + 1 + 6;
12153 cal_RBX = kphi / 6 + 1 + 12;
12159 cal_RBX = kphi / 12 + 1 + 24;
12165 cal_RBX = kphi / 12 + 1 + 30;
12186 for (
int k1 = 0; k1 <
ndepth; k1++) {
12187 for (
int k2 = 0; k2 <
neta; k2++) {
12188 for (
int k3 = 0; k3 <
nphi; k3++) {
12189 int ieta = k2 - 41;
12284 int kkkk2 = (k2 - 1) / 4;
12289 int kkkk3 = (k3) / 4 + 1;
12291 ietaphi = ((kkkk2)-1) *
znphi + (kkkk3);
12294 bbb3 = bbbc / bbb1;
13138 std::cout <<
" ==== Edn of run " << std::endl;
13156 MAPfile <<
"#define LogEleMapdb_h" << std::endl;
13157 MAPfile <<
"#include <algorithm>" << std::endl;
13158 MAPfile <<
"#include <iostream>" << std::endl;
13159 MAPfile <<
"#include <vector>" << std::endl;
13160 MAPfile <<
"#include <string>" << std::endl;
13161 MAPfile <<
"#include <sstream>" << std::endl;
13164 MAPfile <<
"struct Cell {" << std::endl;
13165 MAPfile <<
" std::string subdet;" << std::endl;
13166 MAPfile <<
" int Eta;" << std::endl;
13167 MAPfile <<
" int Phi;" << std::endl;
13168 MAPfile <<
" int Depth;" << std::endl;
13169 MAPfile <<
" std::string RBX;" << std::endl;
13170 MAPfile <<
" int RM;" << std::endl;
13171 MAPfile <<
" int Pixel;" << std::endl;
13172 MAPfile <<
" int RMfiber;" << std::endl;
13173 MAPfile <<
" int FiberCh;" << std::endl;
13174 MAPfile <<
" int QIE;" << std::endl;
13175 MAPfile <<
" int ADC;" << std::endl;
13176 MAPfile <<
" int VMECardID;" << std::endl;
13177 MAPfile <<
" int dccID;" << std::endl;
13178 MAPfile <<
" int Spigot;" << std::endl;
13179 MAPfile <<
" int FiberIndex;" << std::endl;
13180 MAPfile <<
" int HtrSlot;" << std::endl;
13181 MAPfile <<
" int HtrTB;" << std::endl;
13184 MAPfile <<
"// the function check, if \"par\" == \"val\" for this cell" << std::endl;
13185 MAPfile <<
" bool check(const std::string par, const int val) const " << std::endl;
13186 MAPfile <<
" {" << std::endl;
13187 MAPfile <<
" if (par == \"Eta\") return (val == Eta);" << std::endl;
13188 MAPfile <<
" else if (par == \"Phi\") return (val == Phi);" << std::endl;
13189 MAPfile <<
" else if (par == \"Depth\") return (val == Depth);" << std::endl;
13190 MAPfile <<
" else if (par == \"RM\") return (val == RM);" << std::endl;
13191 MAPfile <<
" else if (par == \"Pixel\") return (val == Pixel);" << std::endl;
13192 MAPfile <<
" else if (par == \"RMfiber\") return (val == RMfiber);" << std::endl;
13193 MAPfile <<
" else if (par == \"FiberCh\") return (val == FiberCh);" << std::endl;
13194 MAPfile <<
" else if (par == \"QIE\") return (val == QIE);" << std::endl;
13195 MAPfile <<
" else if (par == \"ADC\") return (val == ADC);" << std::endl;
13196 MAPfile <<
" else if (par == \"VMECardID\") return (val == VMECardID);" << std::endl;
13197 MAPfile <<
" else if (par == \"dccID\") return (val == dccID);" << std::endl;
13198 MAPfile <<
" else if (par == \"Spigot\") return (val == Spigot);" << std::endl;
13199 MAPfile <<
" else if (par == \"FiberIndex\") return (val == FiberIndex);" << std::endl;
13200 MAPfile <<
" else if (par == \"HtrSlot\") return (val == HtrSlot);" << std::endl;
13201 MAPfile <<
" else if (par == \"HtrTB\") return (val == HtrTB);" << std::endl;
13202 MAPfile <<
" else return false;" << std::endl;
13203 MAPfile <<
" }" << std::endl;
13206 MAPfile <<
" bool check(const std::string par, const std::string val) const" << std::endl;
13207 MAPfile <<
" {" << std::endl;
13208 MAPfile <<
" if (par == \"subdet\") return (val == subdet);" << std::endl;
13209 MAPfile <<
" else if (par == \"RBX\") return (val == RBX);" << std::endl;
13210 MAPfile <<
" else return false;" << std::endl;
13211 MAPfile <<
" }" << std::endl;
13213 MAPfile <<
"};" << std::endl;
13216 MAPfile <<
"const Cell AllCells[] = {" << std::endl;
13217 MAPfile <<
"//{ SD, Eta, Phi, Depth, RBX, RM, PIXEL, RMfiber, Fiber Ch., QIE, ADC, VMECrateId, dccid, spigot, " 13218 "fiberIndex, htrSlot, htrTopBottom }" 13231 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13232 << detid->
depth() <<
" ,";
13233 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13234 << lmap_entry.
rmFiber() <<
" , ";
13238 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13240 MAPfile <<
"}," << std::endl;
13253 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() - 1 <<
" , " << detid->
iphi() - 1 <<
" , " 13254 << detid->
depth() <<
" ,";
13255 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13256 << lmap_entry.
rmFiber() <<
" , ";
13260 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13262 MAPfile <<
"}," << std::endl;
13269 for (
int eta = -20;
eta <= -20;
eta++) {
13277 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13278 << detid->
depth() <<
" ,";
13279 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13280 << lmap_entry.
rmFiber() <<
" , ";
13284 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13286 MAPfile <<
"}," << std::endl;
13292 for (
int eta = -19;
eta <= -16;
eta++) {
13300 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13301 << detid->
depth() <<
" ,";
13302 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13303 << lmap_entry.
rmFiber() <<
" , ";
13307 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13309 MAPfile <<
"}," << std::endl;
13314 for (
int eta = -29;
eta <= -16;
eta++) {
13322 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13323 << detid->
depth() <<
" ,";
13324 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13325 << lmap_entry.
rmFiber() <<
" , ";
13329 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13331 MAPfile <<
"}," << std::endl;
13344 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13345 << detid->
depth() <<
" ,";
13346 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13347 << lmap_entry.
rmFiber() <<
" , ";
13351 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13353 MAPfile <<
"}," << std::endl;
13361 for (
int eta = -41;
eta <= -29;
eta++) {
13369 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13370 << detid->
depth() <<
" ,";
13371 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13372 << lmap_entry.
rmFiber() <<
" , ";
13376 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13378 MAPfile <<
"}," << std::endl;
13392 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13393 << detid->
depth() <<
" ,";
13394 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13395 << lmap_entry.
rmFiber() <<
" , ";
13399 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13401 MAPfile <<
"}," << std::endl;
13417 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13418 << detid->
depth() <<
" ,";
13419 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13420 << lmap_entry.
rmFiber() <<
" , ";
13424 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13426 MAPfile <<
"}," << std::endl;
13440 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() - 1 <<
" , " << detid->
iphi() - 1 <<
" , " 13441 << detid->
depth() <<
" ,";
13442 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13443 << lmap_entry.
rmFiber() <<
" , ";
13447 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13449 MAPfile <<
"}," << std::endl;
13454 MAPfile <<
"};" << std::endl;
13457 MAPfile <<
"// macro for array length calculation" << std::endl;
13458 MAPfile <<
"#define DIM(a) (sizeof(a)/sizeof(a[0]))" << std::endl;
13461 MAPfile <<
"// class for cells array managing" << std::endl;
13462 MAPfile <<
"class CellDB {" << std::endl;
13463 MAPfile <<
"public:" << std::endl;
13464 MAPfile <<
" CellDB()" << std::endl;
13465 MAPfile <<
" : cells(AllCells, AllCells + DIM(AllCells))" << std::endl;
13466 MAPfile <<
"{}" << std::endl;
13469 MAPfile <<
"// return i-th cell" << std::endl;
13470 MAPfile <<
"Cell operator [] (int i) const {return cells[i];}" << std::endl;
13472 MAPfile <<
"// number of cells in database" << std::endl;
13473 MAPfile <<
"int size() const {return cells.size();}" << std::endl;
13476 MAPfile <<
"// select cells for which \"par\" == \"val\"" << std::endl;
13477 MAPfile <<
"template<typename T>" << std::endl;
13478 MAPfile <<
"CellDB find(const std::string par, const T val) const" << std::endl;
13480 MAPfile <<
" std::vector<Cell> s;" << std::endl;
13481 MAPfile <<
" for (size_t i = 0; i < cells.size(); ++i)" << std::endl;
13482 MAPfile <<
" if (cells[i].check(par, val))" << std::endl;
13483 MAPfile <<
" s.push_back(cells[i]);" << std::endl;
13484 MAPfile <<
" return CellDB(s);" << std::endl;
13485 MAPfile <<
"} " << std::endl;
13487 MAPfile <<
"private:" << std::endl;
13488 MAPfile <<
" CellDB( const std::vector<Cell> s)" << std::endl;
13489 MAPfile <<
" : cells(s)" << std::endl;
13490 MAPfile <<
"{}" << std::endl;
13491 MAPfile <<
"std::vector<Cell> cells;" << std::endl;
13492 MAPfile <<
"};" << std::endl;
13494 std::cout <<
"===== Finish writing Channel MAP =====" << std::endl;
13498 double deltaphi = phi1 - phi2;
13500 deltaphi = phi2 - phi1;
13502 if (deltaphi >
M_PI) {
13503 deltaphi = 2. *
M_PI - deltaphi;
13506 double tmp =
sqrt(deltaeta * deltaeta + deltaphi * deltaphi);
13516 if (a1 > 0.5 *
M_PI &&
a2 < 0.)
13518 if (
a2 > 0.5 *
M_PI && a1 < 0.)
13520 double tmp = (a1 * en1 +
a2 * en2) / (en1 + en2);
13533 double tmp =
a2 - a1;
13534 if (a1 *
a2 < 0.) {
13535 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
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 * h_AmplitudeHBrest1
TH1F * h_sumamplitude_depth2_HE0
TH2F * h_mapCapCalib047_HE
TH1F * h_sumCutTSmeanAperLS4
double lsdep_estimator2_HBdepth1_
TH2F * h_mapTSmaxCalib_HE
TH2F * h_mapTSmeanCalib047_HO
TH2F * h_amplitudechannel1_HF3
TH2F * h_mapDepth3Error_HB
TH1F * h_Amplitude_notCapIdErrors_HO4
TH2F * h_mapGetRMSOverNormalizedSignal_HF
TH1F * h_sumCutAmplitudeperLS1
TH2F * h_amplitudechannel1_HB2
double calib0[ndepth][neta][nphi]
TH2F * h_amplitudechannel0_HB3
TH1F * h_sumamplitudechannel_HB
edm::EDGetTokenT< HBHERecHitCollection > tok_hbheNoise_
TH2F * h_mapDepth4Pedw0_HO
double lsdep_estimator2_HOdepth4_
float calibt[ndepth][neta][nphi]
TH1F * h_sumamplitude_depth1_HE
TH1F * h_sumCutTSmaxAperLS1
TH1F * h_sumADCAmplperLSdepth4HBu
TH2F * h_mapDepth1Pedw2_HE
TH1F * h_sumCutAmplperLS7
TH1F * h_sumCutADCAmplperLSdepth4HEu
TH2F * h_mapDepth3pedestalw_HB
TH2F * h_recNoiseEnergy2_HF1
TH1F * h_pedestalwaver4_HB
TH1F * h_ADCAmpl345Zoom1_HE
TH2F * h_recSignalEnergy0_HF4
double lsdep_estimator3_HBdepth2_
TH2F * h_recNoiseEnergy1_HB4
CMTRawAnalyzer(const edm::ParameterSet &)
TH2F * h_mapDepth3Ped3_HE
TH2F * h_recSignalEnergy0_HE6
TH2F * h_amplitudechannel1_HF4
TH1F * h_averSIGNALoccupancy_HB
TH2F * h_mapDepth2TSmaxA225_HE
TH2F * h_mapDepth2Amplitude225_HF
TH1F * h_sumTSmeanAperLS4
TH2F * h_mapDepth4Amplitude_HF
TH2F * h_mapDepth2Error_HF
void fillDigiAmplitudeHO(HODigiCollection::const_iterator &digiItr)
TH1F * h_sumADCAmplperLS8
TH1F * h_sumADCAmplEtaPhiLs
TH2F * h_mapDepth4ADCAmpl12_HF
TH1F * h_runbadrate_depth2_HE
TH1F * h_Eta16Depth3RADDAM_HED2P
TH2F * h_amplitudechannel0_HE7
TH2F * h_mapDepth2Pedw3_HF
TH2F * h_mapDepth4ADCAmpl225Copy_HB
double averSIGNALsumamplitude_HO
TH2F * h_2DsumADCAmplLSdepth3HFu
void fillDigiErrors(HBHEDigiCollection::const_iterator &digiItr)
TH2F * h_recNoiseEnergy1_HE6
double lsdep_estimator3_HBdepth1_
TH1F * h_numberofhitsHEtest
double sum0Estimator[nsub][ndepth][neta][nphi]
TH2F * h_mapAiForLS_bad_HE
TH1F * h_AmplitudeHEtest1
double lsdep_estimator1_HFdepth3_
TH2F * h_2DsumADCAmplLSdepth5HEu
TH1F * h_sumamplitude_depth4_HO
TH1F * h_pedestalaver9_HO
TH2F * h_2D0sumTSmeanALS7
TH2F * h_2DsumAmplitudeLS1
TH1F * h_sum0ADCAmplperLS7u
TH1F * h_runnbadchannels_depth3_HE
TH2F * h_amplitudechannel2_HF1
TH2F * h_mapDepth1ADCAmpl12_HF
TH1F * h_AmplitudeHBrest6
TH2F * h_mapDepth4AmplE34_HE
TH2F * h_mapDepth3ADCAmpl12_HF
TH1F * h_sigLayer1RADDAM6_HED2
double mapRADDAM_HED2[ndepth][neta]
TH2F * h_maprphinorm0_HE4
TH2F * h_mapDepth1Pedw1_HE
TH2F * h_2D0sumPedestalLS1
TH2F * h_mapAiForLS_good_HE
TH2F * h_2DsumADCAmplLSdepth4HBu
TH2F * h_mapDepth2ADCAmpl225_HE
TH1F * h_bcnnbadchannels_depth2_HF
TH2F * h_mapDepth1Ped1_HE
TH2F * h_mapDepth2ADCAmpl225_HB
TH2F * h_2DsumPedestalLS4
TH2F * h_recSignalEnergy1_HE2
TH2F * h_2D0sumADCAmplLS1
TH1F * h_A_Depth2RADDAM_HE
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
#define DEFINE_FWK_MODULE(type)
TH2F * h_2DsumAmplitudeLS6
int mapRADDAM0_HE[ndepth][neta][nphi]
TH1F * h_sum0ErrorBperLS4
TH1F * h_runbadrate_depth1_HB
double lsdep_estimator3_HOdepth4_
TH2F * h_recNoiseEnergy2_HE5
TH2F * h_recNoiseEnergy2_HE6
TH1F * h_gsmdifferencefit5_HE
TH2F * h_2DAtaildepth2_HB
TH2F * h_mapDepth2Amplitude225_HE
int lsdep_cut3_max_HEdepth2_
TH1F * h_Amplitude_notCapIdErrors_HB1
TH1F * h_runnbadchannelsC_depth1_HE
TH2F * h_recNoiseEnergy0_HF3
double amplitudechannel[nsub][ndepth][neta][nphi]
double averSIGNALoccupancy_HO
TH1F * h_Amplitude_forCapIdErrors_HB1
TH2F * h_2DsumPedestalLS7
TH1F * h_runbadrateC_depth2_HB
TH2F * h_mapDepth3Error_HE
TH2F * h_mapCapCalib047_HO
double lsdep_estimator4_HFdepth2_
TH2F * h_mapCapCalib047_HB
TH2F * h_mapDepth2Ampl_HE
TH1F * h_pedestalwaver4_HO
TH2F * h_mapDepth3Amplitude225_HB
TH2F * h_mapDepth4Amplitude225_HB
TH2F * h_recSignalEnergy1_HF2
TH2F * h_amplitudechannel0_HE2
TH2F * h_mapDepth2Pedw3_HE
int lsdep_cut3_max_HEdepth1_
TH2F * h_mapDepth4Ped1_HO
TH1F * h_runnbadchannelsC_depth1_HF
TH2F * h_mapDepth4Ampl047_HF
TH1F * h_sumamplitude_depth4_HO0
TH1F * h_shape0_bad_channels_HB
TH1F * h_RatioOccupancy_HEP
TH2F * h_mapWidthCalib_HO
TH1F * h_Eta16Depth3RADDAM_HED2ALL
TH2F * h_amplitudechannel1_HE5
TH2F * h_mapDepth1TSmaxA_HF
TH2F * h_recNoiseEnergy0_HE6
TH1F * h_sigLayer1RADDAM6_HE
TH1F * h_shape0_good_channels_HB
double calibrTSmeanHOMax_
TH2F * h_2D0sumPedestalLS4
TH2F * h_2DsumPedestalLS1
TH1F * h_sigLayer1RADDAM0_HE
TH1F * h_sum0ErrorBperLS5
TH2F * h_mapDepth7ADCAmpl_HE
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 * 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
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 * 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
TH2F * h_mapTSmeanCalib047_HF
TH1F * h_totalAmplitudeHO
TH1F * h_numberofhitsHFtest
TH2F * h_mapTSmaxCalib047_HE
TH1F * h_runbadrate0_depth1_HF
TH2F * h_mapDepth2Pedw0_HF
TH2F * h_maprphinorm2_HE1
TH2F * h_amplitudechannel0_HF1
double calibrTSmeanHOMin_
TH1F * h_NphiForEta16Depth3RADDAM_HED2P
TH1F * h_sigLayer2RADDAM_HE
TH2F * h_mapDepth4ADCAmpl_HO
TH1F * h_sum0ADCAmplperLS6_P1
TH2F * h_recNoiseEnergy0_HB1
TH2F * h_amplitudechannel0_HF3
TH2F * h_mapDepth2ADCAmpl225Copy_HB
TH2F * h_mapDepth7pedestalw_HE
TH1F * h_sumCutADCAmplperLSdepth3HBu
TH2F * h_amplitudechannel1_HF1
TH1F * h_ADCAmpl_HBCapIdNoError
void analyze(const edm::Event &, const edm::EventSetup &) override
double caliba[ndepth][neta][nphi]
TH2F * h_2D0sumAmplitudeLS8
TH1F * h_sumCutADCAmplperLS3
double averNOSIGNALoccupancy_HE
TH2F * h_recSignalEnergy1_HB3
TH2F * h_recSignalEnergy0_HB3
TH2F * h_mapDepth6ADCAmpl225_HE
bool studyCalibCellsHist_
TH2F * h_recSignalEnergy1_HB1
TH2F * h_recNoiseEnergy1_HF1
Abs< T >::type abs(const T &t)
TH1F * h_Amplitude_notCapIdErrors_HB2
TH1F * h_energyhitSignal_HE
TH1F * h_runbadrateC_depth1_HE
TH1F * h_orbitNumvsamplitude0_HF
double sumEstimator1[nsub][ndepth][neta][nphi]
TH2F * h_2DsumAmplitudeLS8
TH2F * h_mapDepth3Amplitude225_HF
TH2F * h_mapDepth3TSmeanA225_HE
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
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
bool getData(T &iHolder) const
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
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
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
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
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
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
static constexpr float b1
TH2F * h_mapDepth3TSmeanA_HE
const edm::ESGetToken< HcalDbService, HcalDbRecord > tokDB_
TH2F * h_mapDepth4pedestal_HO
TH2F * h_mapDepth4TSmeanA_HE
TH1F * h_sumCutAmplperLS3
TH2F * h_mapDepth1Amplitude_HB
double dR(double eta1, double phi1, double eta2, double phi2)
uint16_t *__restrict__ uint16_t const *__restrict__ adc
TH1F * h_sumPedestalperLS8
TH2F * h_mapDepth4Amplitude225_HE
edm::EDGetTokenT< HcalCalibDigiCollection > tok_calib_
TH2F * h_recSignalEnergy2_HE5
TH1F * h_orbitNumvsamplitude_HE
TH1F * h_sumADCAmplperLSdepth7HEu
double recNoiseEnergy2[nsub][ndepth][neta][nphi]
TH2F * h_mapDepth3RADDAM0_HE
TH2F * h_mapDepth2pedestal_HE
void fillDigiAmplitude(HBHEDigiCollection::const_iterator &digiItr)
TH2F * h_amplitudechannel0_HE5
TH1F * h_averSIGNALoccupancy_HE
TH1F * h_runnbadchannels_depth2_HF
TH2F * h_2DsumADCAmplLSdepth4HFu
TH1F * h_bcnvsamplitude0_HO
TH1F * h_sumCutADCAmplperLSdepth5HEu
TH1F * h_nbadchannels_depth1_HF
TH2F * h_recSignalEnergy0_HE3
TH2F * h_mapDepth1Ped0_HE
TH2F * h_map0AiForLS_good_HB
TH2F * h_mapRatioCalib047_HE
TH2F * h_amplitudechannel2_HE2
TH2F * h_mapDepth3TSmaxA_HE
TH1F * h_sumCutADCAmplperLS7
TH2F * h_2D0sumADCAmplLS4
TH1F * h_sumCutAmplitudeperLS4
constexpr int depth() const
get the tower depth
bool studyCapIDErrorsHist_
bool studyPedestalCorrelations_