351 1.89, 5.07, 8.25, 11.43, 14.61, 17.78, 20.96, 24.14, 27.32, 30.50, 33.68,
352 36.86, 40.04, 43.22, 46.40, 49.58, 54.35, 60.71, 67.07, 73.43, 79.79, 86.15,
353 92.51, 98.87, 105.2, 111.6, 117.9, 124.3, 130.7, 137.0, 143.4, 149.7, 156.1,
354 162.5, 168.8, 175.2, 184.7, 197.4, 210.2, 222.9, 235.6, 248.3, 261.0, 273.7,
355 286.5, 299.2, 311.9, 324.6, 337.3, 350.1, 362.8, 375.5, 388.2, 400.9, 413.6,
356 426.4, 439.1, 458.2, 483.6, 509.0, 534.5, 559.9, 585.3, 610.8, 558.9, 584.2,
357 609.5, 634.7, 660.0, 685.3, 710.6, 735.9, 761.2, 786.5, 811.8, 837.1, 862.4,
358 887.7, 913.0, 938.3, 976.2, 1026.8, 1077.4, 1128.0, 1178.6, 1229.1, 1279.7, 1330.3,
359 1380.9, 1431.5, 1482.1, 1532.7, 1583.3, 1633.8, 1684.4, 1735.0, 1785.6, 1836.2, 1886.8,
360 1937.4, 2013.2, 2114.4, 2215.6, 2316.8, 2417.9, 2519.1, 2620.3, 2721.5, 2822.6, 2923.8,
361 3025.0, 3126.2, 3227.3, 3328.5, 3429.7, 3530.9, 3632.0, 3733.2, 3834.4, 3935.5, 4036.7,
362 4188.5, 4390.8, 4593.2, 4795.5, 4997.9, 5200.2, 5402.6, 5057.5, 5262.3, 5467.1, 5671.8,
363 5876.6, 6081.4, 6286.2, 6491.0, 6695.8, 6900.6, 7105.3, 7310.1, 7514.9, 7719.7, 7924.5,
364 8129.3, 8436.4, 8846.0, 9255.6, 9665.1, 10074.7, 10484.3, 10893.9, 11303.4, 11713.0, 12122.6,
365 12532.1, 12941.7, 13351.3, 13760.8, 14170.4, 14580.0, 14989.5, 15399.1, 15808.7, 16218.2, 16832.6,
366 17651.7, 18470.9, 19290.0, 20109.2, 20928.3, 21747.4, 22566.6, 23385.7, 24204.8, 25024.0, 25843.1,
367 26662.3, 27481.4, 28300.5, 29119.7, 29938.8, 30757.9, 31577.1, 32396.2, 33215.4, 34444.1, 36082.3,
368 37720.6, 39358.9, 40997.2, 42635.4, 44273.7, 40908.7, 42542.6, 44176.5, 45810.4, 47444.3, 49078.3,
369 50712.2, 52346.1, 53980.0, 55613.9, 57247.8, 58881.8, 60515.7, 62149.6, 63783.5, 65417.4, 67868.3,
370 71136.1, 74404.0, 77671.8, 80939.7, 84207.5, 87475.3, 90743.2, 94011.0, 97278.8, 100546.7, 103814.5,
371 107082.3, 110350.2, 113618.0, 116885.8, 120153.7, 123421.5, 126689.3, 129957.2, 134858.9, 141394.6, 147930.3,
372 154465.9, 161001.6, 167537.3, 174072.9, 180608.6, 187144.3, 193679.9, 200215.6, 206751.3, 213287.0, 219822.6,
373 226358.3, 232894.0, 239429.6, 245965.3, 252501.0, 259036.6, 265572.3, 275375.8, 288447.2, 301518.5, 314589.8,
374 327661.2, 340732.5, 353803.8};
377 9.56, 28.24, 46.91, 65.59, 84.27, 102.9, 121.6, 140.3, 159.0, 177.7,
378 196.3, 215.0, 233.7, 252.4, 271.0, 289.7, 317.7, 355.1, 392.4, 429.8,
379 467.1, 504.5, 541.9, 579.2, 616.6, 653.9, 691.3, 728.6, 766.0, 803.3,
380 840.7, 878.0, 915.4, 952.8, 990.1, 1027.5, 1083.5, 1158.2, 1232.9, 1307.6,
381 1382.3, 1457.0, 1531.7, 1606.4, 1681.2, 1755.9, 1830.6, 1905.3, 1980.0, 2054.7,
382 2129.4, 2204.1, 2278.8, 2353.5, 2428.2, 2502.9, 2577.7, 2689.7, 2839.1, 2988.6,
383 3138.0, 3287.4, 3436.8, 3586.2, 3263.0, 3411.3, 3559.6, 3707.9, 3856.2, 4004.5,
384 4152.9, 4301.2, 4449.5, 4597.8, 4746.1, 4894.4, 5042.7, 5191.0, 5339.4, 5487.7,
385 5710.1, 6006.8, 6303.4, 6600.0, 6896.6, 7193.3, 7489.9, 7786.5, 8083.1, 8379.8,
386 8676.4, 8973.0, 9269.6, 9566.3, 9862.9, 10159.5, 10456.2, 10752.8, 11049.4, 11346.0,
387 11791.0, 12384.2, 12977.5, 13570.7, 14164.0, 14757.2, 15350.5, 15943.7, 16537.0, 17130.2,
388 17723.5, 18316.7, 18910.0, 19503.2, 20096.5, 20689.7, 21283.0, 21876.2, 22469.5, 23062.8,
389 23656.0, 24545.9, 25732.4, 26918.9, 28105.4, 29291.9, 30478.4, 31664.9, 29399.4, 30590.1,
390 31780.9, 32971.7, 34162.4, 35353.2, 36544.0, 37734.7, 38925.5, 40116.3, 41307.0, 42497.8,
391 43688.5, 44879.3, 46070.1, 47260.8, 49047.0, 51428.5, 53810.1, 56191.6, 58573.1, 60954.6,
392 63336.2, 65717.7, 68099.2, 70480.8, 72862.3, 75243.8, 77625.4, 80006.9, 82388.4, 84769.9,
393 87151.5, 89533.0, 91914.5, 94296.1, 97868.4, 102631.4, 107394.5, 112157.5, 116920.6, 121683.7,
394 126446.7, 131209.8, 135972.8, 140735.9, 145499.0, 150262.0, 155025.1, 159788.2, 164551.2, 169314.3,
395 174077.3, 178840.4, 183603.5, 188366.5, 193129.6, 200274.2, 209800.3, 219326.4, 228852.5, 238378.7,
396 247904.8, 257430.9, 237822.7, 247326.7, 256830.7, 266334.8, 275838.8, 285342.9, 294846.9, 304351.0,
397 313855.0, 323359.1, 332863.1, 342367.1, 351871.2, 361375.2, 370879.3, 380383.3, 394639.4, 413647.5,
398 432655.6, 451663.6, 470671.7, 489679.8, 508687.9, 527696.0, 546704.1, 565712.2, 584720.3, 603728.3,
399 622736.4, 641744.5, 660752.6, 679760.7, 698768.8, 717776.9, 736785.0, 755793.0, 784305.2, 822321.3,
400 860337.5, 898353.7, 936369.9, 974386.0, 1012402.2, 1050418.4, 1088434.6, 1126450.7, 1164466.9, 1202483.1,
401 1240499.3, 1278515.4, 1316531.6, 1354547.8, 1392564.0, 1430580.1, 1468596.3, 1506612.5, 1544628.7, 1601652.9,
402 1677685.3, 1753717.6, 1829750.0, 1905782.3, 1981814.7, 2057847.0};
405 -0.5, 0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5, 10.5, 11.5,
406 12.5, 13.5, 15., 17., 19., 21., 23., 25., 27., 29.5, 32.5, 35.5, 38.5,
407 42., 46., 50., 54.5, 59.5, 64.5, 59.5, 64.5, 69.5, 74.5, 79.5, 84.5, 89.5,
408 94.5, 99.5, 104.5, 109.5, 114.5, 119.5, 124.5, 129.5, 137., 147., 157., 167., 177.,
409 187., 197., 209.5, 224.5, 239.5, 254.5, 272., 292., 312., 334.5, 359.5, 384.5, 359.5,
410 384.5, 409.5, 434.5, 459.5, 484.5, 509.5, 534.5, 559.5, 584.5, 609.5, 634.5, 659.5, 684.5,
411 709.5, 747., 797., 847., 897., 947., 997., 1047., 1109.5, 1184.5, 1259.5, 1334.5, 1422.,
412 1522., 1622., 1734.5, 1859.5, 1984.5, 1859.5, 1984.5, 2109.5, 2234.5, 2359.5, 2484.5, 2609.5, 2734.5,
413 2859.5, 2984.5, 3109.5, 3234.5, 3359.5, 3484.5, 3609.5, 3797., 4047., 4297., 4547., 4797., 5047.,
414 5297., 5609.5, 5984.5, 6359.5, 6734.5, 7172., 7672., 8172., 8734.5, 9359.5, 9984.5};
438 double phi12(
double phi1,
double en1,
double phi2,
double en2);
439 double dR(
double eta1,
double phi1,
double eta2,
double phi2);
440 double dPhiWsign(
double phi1,
double phi2);
2385 std::cout <<
"============== endJob ===================================" << std::endl;
2387 std::cout <<
" --------------------------------------- " << std::endl;
2390 std::cout <<
" --------------------------------------------- " << std::endl;
2396 std::cout <<
"===== possible max number of events in ntuple(each 50th recorded) = " <<
nevent50 << std::endl;
2400 std::cout <<
"===== full number of events*HBHEdigis = " <<
nnnnnn << std::endl;
2405 std::cout <<
"===== Start writing user histograms =====" << std::endl;
2640 for (
int k1 = 0; k1 <
ndepth; k1++) {
2641 for (
int k2 = 0; k2 <
neta; k2++) {
2646 for (
int k3 = 0; k3 <
nphi; k3++) {
2705 int outabortgap = 1;
2715 std::cout <<
" --------------------------------------- " << std::endl;
2719 std::cout <<
" --------------------------------------------- " << std::endl;
2725 std::cout <<
" ---------***********************------------- " << std::endl;
2775 unsigned long int pcountall1 = 0;
2776 unsigned long int pcountall3 = 0;
2777 unsigned long int pcountall6 = 0;
2778 unsigned long int pcountall8 = 0;
2783 unsigned long int mcountall1 = 0;
2784 unsigned long int mcountall3 = 0;
2785 unsigned long int mcountall6 = 0;
2786 unsigned long int mcountall8 = 0;
2791 int pnnmin1 = 999999999;
2792 int pnnmin3 = 999999999;
2793 int pnnmin6 = 999999999;
2794 int pnnmin8 = 999999999;
2799 int mnnmin1 = 999999999;
2800 int mnnmin3 = 999999999;
2801 int mnnmin6 = 999999999;
2802 int mnnmin8 = 999999999;
2804 for (
int k1 = 0; k1 <
ndepth; k1++) {
2805 for (
int k3 = 0; k3 <
nphi; k3++) {
2806 for (
int k2 = 0; k2 <
neta; k2++) {
2915 int kkkk2 = (k2 - 1) / 4;
2920 int kkkk3 = (k3) / 4 + 1;
2924 ietaphi = ((kkkk2)-1) *
znphi + (kkkk3);
2992 if (bbbc / bbb1 > 25.) {
3007 if (bbbc / bbb1 > 25.) {
3092 if (bbbc / bbb1 > 15. && k3 % 2 == 0) {
3106 if (bbbc / bbb1 > 15. && k3 % 2 == 0) {
3246 if (bbbc / bbb1 > 20.) {
3261 if (bbbc / bbb1 > 20.) {
3353 if (bbbc / bbb1 > 80.) {
3368 if (bbbc / bbb1 > 80.) {
3923 if (
k0 == 0 && k1 == 0) {
3924 if (pcountmin1 > 0) {
3925 if (pcountmin1 < pnnmin1)
3926 pnnmin1 = pcountmin1;
3929 if (mcountmin1 > 0) {
3930 if (mcountmin1 < mnnmin1)
3931 mnnmin1 = mcountmin1;
3937 if (
k0 == 1 && k1 == 0) {
3938 if (pcountmin3 > 0) {
3939 if (pcountmin3 < pnnmin3)
3940 pnnmin3 = pcountmin3;
3943 if (mcountmin3 > 0) {
3944 if (mcountmin3 < mnnmin3)
3945 mnnmin3 = mcountmin3;
3951 if (
k0 == 2 && k1 == 3) {
3952 if (pcountmin8 > 0) {
3953 if (pcountmin8 < pnnmin8)
3954 pnnmin8 = pcountmin8;
3957 if (mcountmin8 > 0) {
3958 if (mcountmin8 < mnnmin8)
3959 mnnmin8 = mcountmin8;
3965 if (
k0 == 3 && k1 == 0) {
3966 if (pcountmin6 > 0) {
3967 if (pcountmin6 < pnnmin6)
3968 pnnmin6 = pcountmin6;
3971 if (mcountmin6 > 0) {
3972 if (mcountmin6 < mnnmin6)
3973 mnnmin6 = mcountmin6;
3988 float patiooccupancy1 = 0.;
3989 if (pcountall1 != 0)
3990 patiooccupancy1 = (
float)pnnmin1 * mnnbins1 / pcountall1;
3992 float matiooccupancy1 = 0.;
3993 if (mcountall1 != 0)
3994 matiooccupancy1 = (
float)mnnmin1 * mnnbins1 / mcountall1;
3997 float patiooccupancy3 = 0.;
3998 if (pcountall3 != 0)
3999 patiooccupancy3 = (
float)pnnmin3 * mnnbins3 / pcountall3;
4001 float matiooccupancy3 = 0.;
4002 if (mcountall3 != 0)
4003 matiooccupancy3 = (
float)mnnmin3 * mnnbins3 / mcountall3;
4006 float patiooccupancy6 = 0.;
4007 if (pcountall6 != 0)
4008 patiooccupancy6 = (
float)pnnmin6 * mnnbins6 / pcountall6;
4010 float matiooccupancy6 = 0.;
4011 if (mcountall6 != 0)
4012 matiooccupancy6 = (
float)mnnmin6 * mnnbins6 / mcountall6;
4015 float patiooccupancy8 = 0.;
4016 if (pcountall8 != 0)
4017 patiooccupancy8 = (
float)pnnmin8 * mnnbins8 / pcountall8;
4019 float matiooccupancy8 = 0.;
4020 if (mcountall8 != 0)
4021 matiooccupancy8 = (
float)mnnmin8 * mnnbins8 / mcountall8;
4025 for (
int k1 = 0; k1 <
ndepth; k1++) {
4026 for (
int k2 = 0; k2 <
neta; k2++) {
4027 for (
int k3 = 0; k3 <
nphi; k3++) {
4127 for (
int k1 = 0; k1 <
ndepth; k1++) {
4128 for (
int k2 = 0; k2 <
neta; k2++) {
4129 for (
int k3 = 0; k3 <
nphi; k3++) {
4151 for (
int k1 = 0; k1 <
ndepth; k1++) {
4152 for (
int k2 = 0; k2 <
neta; k2++) {
4153 for (
int k3 = 0; k3 <
nphi; k3++) {
4185 bool gotHFDigis =
true;
4189 if (!(
hf.isValid())) {
4193 std::cout <<
" ****************************** =========================== No HFDigiCollection found " 4198 eta = digi->id().ieta();
4199 phi = digi->id().iphi();
4200 depth = digi->id().depth();
4217 for (
int i = 0;
i <
nTS;
i++) {
4236 bool gotQIE10Digis =
true;
4238 gotQIE10Digis =
false;
4241 gotQIE10Digis =
false;
4243 if (!gotQIE10Digis) {
4244 std::cout <<
" No QIE10DigiCollection collection is found " << std::endl;
4247 double totalAmplitudeHF = 0.;
4248 for (
unsigned int j = 0;
j < qie10dc.
size();
j++) {
4272 double amplitudefullTSs = 0.;
4273 double nnnnnnTS = 0.;
4274 for (
int i = 0;
i <
nTS; ++
i) {
4277 int adc = qie10df[
i].adc();
4309 bool gotHBHEDigis =
true;
4311 gotHBHEDigis =
false;
4312 if (!(
hbhe.isValid()))
4313 gotHBHEDigis =
false;
4314 if (!gotHBHEDigis) {
4315 std::cout <<
" No HBHEDigiCollection collection is found " << std::endl;
4318 double totalAmplitudeHB = 0.;
4319 double totalAmplitudeHE = 0.;
4320 double nnnnnnTSHB = 0.;
4321 double nnnnnnTSHE = 0.;
4324 eta = digi->id().ieta();
4325 phi = digi->id().iphi();
4326 depth = digi->id().depth();
4366 double amplitudefullTSs = 0.;
4369 for (
int i = 0;
i <
nTS;
i++) {
4381 double amplitudefullTSs = 0.;
4384 for (
int i = 0;
i <
nTS;
i++) {
4399 if (totalAmplitudeHB != 0.) {
4404 if (totalAmplitudeHE != 0.) {
4420 bool gotQIE11Digis =
true;
4422 gotQIE11Digis =
false;
4424 gotQIE11Digis =
false;
4425 if (!gotQIE11Digis) {
4426 std::cout <<
" No QIE11DigiCollection collection is found " << std::endl;
4429 double totalAmplitudeHBQIE11 = 0.;
4430 double totalAmplitudeHEQIE11 = 0.;
4431 double nnnnnnTSHBQIE11 = 0.;
4432 double nnnnnnTSHEQIE11 = 0.;
4433 for (
unsigned int j = 0;
j < qie11dc.
size();
j++) {
4442 int sub = hcaldetid.
subdet();
4455 if (sub == 1 &&
depth == 1 && qwert1 == 0) {
4459 if (sub == 1 &&
depth == 2 && qwert2 == 0) {
4463 if (sub == 2 &&
depth == 1 && qwert3 == 0) {
4467 if (sub == 2 &&
depth == 2 && qwert4 == 0) {
4471 if (sub == 2 &&
depth == 3 && qwert5 == 0) {
4483 double amplitudefullTSs1 = 0.;
4484 double amplitudefullTSs6 = 0.;
4486 for (
int i = 0;
i <
nTS; ++
i) {
4487 int adc = qie11df[
i].adc();
4490 amplitudefullTSs1 += charge1;
4491 amplitudefullTSs6 += charge6;
4497 totalAmplitudeHBQIE11 +=
charge;
4504 double amplitudefullTSs1 = 0.;
4505 double amplitudefullTSs6 = 0.;
4507 for (
int i = 0;
i <
nTS;
i++) {
4508 int adc = qie11df[
i].adc();
4511 amplitudefullTSs1 += charge1;
4512 amplitudefullTSs6 += charge6;
4518 totalAmplitudeHEQIE11 +=
charge;
4527 if (totalAmplitudeHBQIE11 != 0.) {
4532 if (totalAmplitudeHEQIE11 != 0.) {
4543 bool gotHODigis =
true;
4546 if (!(
ho.isValid()))
4550 std::cout <<
" No HO collection is found " << std::endl;
4553 double totalAmplitudeHO = 0.;
4555 eta = digi->id().ieta();
4556 phi = digi->id().iphi();
4557 depth = digi->id().depth();
4578 double nnnnnnTS = 0.;
4579 double amplitudefullTSs = 0.;
4581 for (
int i = 0;
i <
nTS;
i++) {
4610 bool gotHBHERecHitsNoise =
true;
4612 gotHBHERecHitsNoise =
4615 gotHBHERecHitsNoise =
false;
4616 if (!gotHBHERecHitsNoise) {
4618 std::cout <<
" No RecHits HBHENoise collection is found " << std::endl;
4623 float icalconst = 1.;
4626 HBHERecHit aHit(hbheItr->id(), hbheItr->eraw() * icalconst, hbheItr->time());
4628 double energyhit = aHit.
energy();
4629 DetId id = (*hbheItr).detid();
4631 int sub = ((hid).
rawId() >> 25) & 0x7;
4640 if (hid.
depth() == 1 && sub == 1 && hid.
iphi() == 25) {
4642 std::cout <<
" Noise,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4643 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4645 if (hid.
depth() == 1 && sub == 2 && hid.
iphi() == 25) {
4647 std::cout <<
" Noise,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4648 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4654 int mdepth = hid.
depth();
4665 bool gotHFRecHitsNoise =
true;
4667 gotHFRecHitsNoise =
false;
4669 gotHFRecHitsNoise =
false;
4670 if (!gotHFRecHitsNoise) {
4672 std::cout <<
" No RecHits HFNoise collection is found " << std::endl;
4676 float icalconst = 1.;
4679 HFRecHit aHit(hfItr->id(), hfItr->energy() * icalconst, hfItr->time());
4680 double energyhit = aHit.
energy();
4681 DetId id = (*hfItr).detid();
4683 int sub = ((hid).
rawId() >> 25) & 0x7;
4689 if (hid.
iphi() == 25) {
4691 std::cout <<
"HF Noise,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4692 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4698 int mdepth = hid.
depth();
4712 bool gotHBHERecHitsSignal =
true;
4714 gotHBHERecHitsSignal =
4717 gotHBHERecHitsSignal =
false;
4718 if (!gotHBHERecHitsSignal) {
4720 std::cout <<
" No RecHits HBHESignal collection is found " << std::endl;
4725 float icalconst = 1.;
4728 HBHERecHit aHit(hbheItr->id(), hbheItr->eraw() * icalconst, hbheItr->time());
4730 double energyhit = aHit.
energy();
4731 DetId id = (*hbheItr).detid();
4733 int sub = ((hid).
rawId() >> 25) & 0x7;
4741 if (hid.
depth() == 1 && sub == 1 && hid.
iphi() == 25) {
4743 std::cout <<
"HBHE Signal,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4744 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4746 if (hid.
depth() == 1 && sub == 2 && hid.
iphi() == 25) {
4748 std::cout <<
"HBHE Signal,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4749 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4755 int mdepth = hid.
depth();
4766 bool gotHFRecHitsSignal =
true;
4768 gotHFRecHitsSignal =
false;
4770 gotHFRecHitsSignal =
false;
4771 if (!gotHFRecHitsSignal) {
4773 std::cout <<
" No RecHits HFSignal collection is found " << std::endl;
4777 float icalconst = 1.;
4780 HFRecHit aHit(hfItr->id(), hfItr->energy() * icalconst, hfItr->time());
4781 double energyhit = aHit.
energy();
4782 DetId id = (*hfItr).detid();
4784 int sub = ((hid).
rawId() >> 25) & 0x7;
4790 if (hid.
iphi() == 25) {
4792 std::cout <<
"HF Signal,sub = " << sub <<
" mdepth = " << hid.
depth() <<
" ieta= " << hid.
ieta()
4793 <<
" iphi= " << hid.
iphi() <<
" energyhit= " << energyhit << std::endl;
4799 int mdepth = hid.
depth();
4829 for (
int k1 = 0; k1 <
ndepth; k1++) {
4831 for (
int k2 = 0; k2 <
neta; k2++) {
4833 int k2plot = k2 - 41;
4836 double sumoverphi = 0;
4837 int nsumoverphi = 0;
4838 for (
int k3 = 0; k3 <
nphi; k3++) {
4843 std::cout <<
"==== nsumoverphi = " << nsumoverphi <<
" sumoverphi = " << sumoverphi
4844 <<
" k1 = " << k1 <<
" k2 = " << k2 <<
" kkk = " << kkk <<
" k3 = " << k3 << std::endl;
4848 for (
int k3 = 0; k3 <
nphi; k3++) {
4849 if (nsumoverphi != 0) {
4852 std::cout <<
"nsumoverphi= " << nsumoverphi <<
" sumoverphi= " << sumoverphi <<
" k1= " << k1
4853 <<
" k2= " << k2 <<
" kkk= " << kkk <<
" k3= " << k3
4854 <<
" maprphinorm= " <<
maprphinorm[
k0][k1][k2][k3] << std::endl;
4929 for (
int k1 = 0; k1 <
ndepth; k1++) {
4931 for (
int k2 = 0; k2 <
neta; k2++) {
4933 int k2plot = k2 - 41;
4935 for (
int k3 = 0; k3 <
nphi; k3++) {
4963 for (
int k1 = 0; k1 <
ndepth; k1++) {
4965 for (
int k2 = 0; k2 <
neta; k2++) {
4967 int k2plot = k2 - 41;
4969 for (
int k3 = 0; k3 <
nphi; k3++) {
5012 for (
int k1 = 0; k1 <
ndepth; k1++) {
5014 for (
int k2 = 0; k2 <
neta; k2++) {
5016 int k2plot = k2 - 41;
5018 for (
int k3 = 0; k3 <
nphi; k3++) {
5061 for (
int k1 = 0; k1 <
ndepth; k1++) {
5063 for (
int k2 = 0; k2 <
neta; k2++) {
5065 int k2plot = k2 - 41;
5067 for (
int k3 = 0; k3 <
nphi; k3++) {
5107 for (
int k1 = 0; k1 <
ndepth; k1++) {
5109 for (
int k2 = 0; k2 <
neta; k2++) {
5111 int k2plot = k2 - 41;
5113 for (
int k3 = 0; k3 <
nphi; k3++) {
5177 for (
int k1 = 0; k1 <
ndepth; k1++) {
5179 for (
int k2 = 0; k2 <
neta; k2++) {
5181 int k2plot = k2 - 41;
5183 for (
int k3 = 0; k3 <
nphi; k3++) {
5231 for (
int k1 = 0; k1 <
ndepth; k1++) {
5232 for (
int k2 = 0; k2 <
neta; k2++) {
5233 for (
int k3 = 0; k3 <
nphi; k3++) {
5245 for (
int k1 = 0; k1 <
ndepth; k1++) {
5246 for (
int k2 = 0; k2 <
neta; k2++) {
5249 if (k1 == 2 && k2 == 25) {
5253 }
else if (k1 == 2 && k2 == 56) {
5268 for (
int k1 = 0; k1 <
ndepth; k1++) {
5269 for (
int k2 = 0; k2 <
neta; k2++) {
5270 if (k1 == 2 && (k2 == 25 || k2 == 56)) {
5273 int k2plot = k2 - 41;
5283 int LLLLLL111111 = 0;
5284 if ((k1 == 0 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 30))
5287 int LLLLLL222222 = 0;
5288 if ((k1 == 1 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 27) ||
5289 (k1 == 2 && fabs(kkk - 41) > 26 && fabs(kkk - 41) < 29))
5292 if (LLLLLL111111 == 1) {
5296 if (LLLLLL222222 == 1) {
5309 for (
int k1 = 0; k1 <
ndepth; k1++) {
5310 for (
int k2 = 0; k2 <
neta; k2++) {
5311 if (k1 == 2 && (k2 == 25 || k2 == 56)) {
5313 int k2plot = k2 - 41;
5320 for (
int k3 = 0; k3 <
nphi; k3++) {
5323 int LLLLLL111111 = 0;
5324 if ((k1 == 0 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 30))
5326 int LLLLLL222222 = 0;
5327 if ((k1 == 1 && fabs(kkk - 41) > 16 && fabs(kkk - 41) < 27) ||
5328 (k1 == 2 && fabs(kkk - 41) > 26 && fabs(kkk - 41) < 29))
5330 if (LLLLLL111111 == 1) {
5334 if (LLLLLL222222 == 1) {
5345 for (
int k1 = 0; k1 <
ndepth; k1++) {
5346 for (
int k2 = 0; k2 <
neta; k2++) {
5349 for (
int k3 = 0; k3 <
nphi; k3++) {
5364 int sumofchannels = 0;
5365 double sumamplitudesubdet = 0.;
5366 int sumofchannels0 = 0;
5367 double sumamplitudesubdet0 = 0.;
5368 for (
int k2 = 0; k2 <
neta; k2++) {
5369 for (
int k3 = 0; k3 <
nphi; k3++) {
5374 if (sumamplitudechannel_HB > 80.) {
5375 sumamplitudesubdet += sumamplitudechannel_HB;
5378 if (sumamplitudechannel_HB > 0.) {
5379 sumamplitudesubdet0 += sumamplitudechannel_HB;
5387 double sumamplitudechannel_HE =
5390 if (sumamplitudechannel_HE > 200.) {
5391 sumamplitudesubdet += sumamplitudechannel_HE;
5394 if (sumamplitudechannel_HE > 0.) {
5395 sumamplitudesubdet0 += sumamplitudechannel_HE;
5405 if (sumamplitudechannel_HO > 1200.) {
5406 sumamplitudesubdet += sumamplitudechannel_HO;
5409 if (sumamplitudechannel_HO > 0.) {
5410 sumamplitudesubdet0 += sumamplitudechannel_HO;
5419 if (sumamplitudechannel_HF > 600.) {
5420 sumamplitudesubdet += sumamplitudechannel_HF;
5423 if (sumamplitudechannel_HF > 0.) {
5424 sumamplitudesubdet0 += sumamplitudechannel_HF;
5437 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM1)
5438 maxxSUM1 = sumamplitudesubdet + sumamplitudesubdet0;
5439 if ((sumofchannels + sumofchannels0) >
maxxOCCUP1)
5445 if ((sumamplitudesubdet + sumamplitudesubdet0) > 60000) {
5446 for (
int k2 = 0; k2 <
neta; k2++) {
5447 for (
int k3 = 0; k3 <
nphi; k3++) {
5466 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM2)
5467 maxxSUM2 = sumamplitudesubdet + sumamplitudesubdet0;
5468 if ((sumofchannels + sumofchannels0) >
maxxOCCUP2)
5478 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM3)
5479 maxxSUM3 = sumamplitudesubdet + sumamplitudesubdet0;
5480 if ((sumofchannels + sumofchannels0) >
maxxOCCUP3)
5490 if ((sumamplitudesubdet + sumamplitudesubdet0) >
maxxSUM4)
5491 maxxSUM4 = sumamplitudesubdet + sumamplitudesubdet0;
5492 if ((sumofchannels + sumofchannels0) >
maxxOCCUP4)
5513 for (
int k1 = 0; k1 <
ndepth; k1++) {
5515 int nbadchannels = 0;
5516 for (
int k2 = 0; k2 <
neta; k2++) {
5517 for (
int k3 = 0; k3 <
nphi; k3++) {
5530 if (nbadchannels != 0)
5543 if (nbadchannels != 0)
5559 if (nbadchannels != 0)
5572 if (nbadchannels != 0)
5585 if (nbadchannels != 0)
5601 if (nbadchannels != 0)
5617 if (nbadchannels != 0)
5630 if (nbadchannels != 0)
5650 bool gotCALIBDigis =
true;
5652 gotCALIBDigis =
false;
5654 if (!(
calib.isValid())) {
5655 gotCALIBDigis =
false;
5657 if (!gotCALIBDigis) {
5660 int cal_det = digi->id().hcalSubdet();
5661 int cal_phi = digi->id().iphi();
5662 int cal_eta = digi->id().ieta();
5663 int cal_cbox = digi->id().cboxChannel();
5667 if (cal_det > 0 && cal_det < 5 && cal_cbox == 0) {
5668 int iphi = cal_phi - 1;
5673 double max_signal = -100.;
5674 int ts_with_max_signal = -100;
5679 for (
int i = 0;
i <
nTS;
i++) {
5680 double ampldefault =
adc2fC[digi->sample(
i).adc() & 0xff];
5681 if (max_signal < ampldefault) {
5682 max_signal = ampldefault;
5683 ts_with_max_signal =
i;
5688 timew += (
i + 1) * ampldefault;
5693 double aveamplitude = -100.;
5696 double aveamplitude1 = aveamplitude - 1;
5701 double ampldefault =
adc2fC[digi->sample(
ii).adc() & 0xff];
5702 double aaaaaa = (
ii + 1) - aveamplitude;
5703 double aaaaaa2 = aaaaaa * aaaaaa;
5704 rmsamp += (aaaaaa2 * ampldefault);
5706 double rmsamplitude = -100.;
5714 if (ts_with_max_signal > -1 && ts_with_max_signal <
nTS)
5716 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 <
nTS)
5717 calib2[cal_det - 1][
ieta + 41][
iphi] +=
adc2fC[digi->sample(ts_with_max_signal + 1).adc() & 0xff];
5718 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 <
nTS)
5719 calib2[cal_det - 1][
ieta + 41][
iphi] +=
adc2fC[digi->sample(ts_with_max_signal + 2).adc() & 0xff];
5720 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 <
nTS)
5721 calib2[cal_det - 1][
ieta + 41][
iphi] +=
adc2fC[digi->sample(ts_with_max_signal - 1).adc() & 0xff];
5722 if (ts_with_max_signal - 2 > -1 && ts_with_max_signal - 2 <
nTS)
5723 calib2[cal_det - 1][
ieta + 41][
iphi] +=
adc2fC[digi->sample(ts_with_max_signal - 2).adc() & 0xff];
5725 bool anycapid =
true;
5728 int error1 = 0, error2 = 0, error3 = 0;
5729 int lastcapid = 0, capid = 0;
5730 for (
int ii = 0;
ii < (*digi).size();
ii++) {
5731 capid = (*digi)[
ii].capid();
5732 bool er = (*digi)[
ii].er();
5733 bool dv = (*digi)[
ii].dv();
5734 if (
ii != 0 && ((lastcapid + 1) % 4) != capid)
5748 if (error1 != 0 || error2 != 0 || error3 != 0)
5763 for (
int k1 = 0; k1 <
nsub; k1++) {
5764 for (
int k2 = 0; k2 <
neta; k2++) {
5765 for (
int k3 = 0; k3 <
nphi; k3++) {
5766 int k2plot = k2 - 41;
5770 if (
signal[k1][k2][k3] > 0.) {
5796 if (k20 > 0 || k20 == 0)
5801 int kk2 = 0, kk3 = 0;
5802 if (k1 == 0 || k1 == 1) {
5807 if (k30 == 71 || k30 ==
nphi || k30 == 1 || k30 == 2)
5810 kk3 = ((k30 - 3) / 4) * 4 + 3;
5811 }
else if (k1 == 2) {
5812 if (
abs(k20) <= 4) {
5814 if (k30 == 71 || k30 ==
nphi || k30 == 1 || k30 == 2 || k30 == 3 || k30 == 4)
5817 kk3 = ((k30 - 5) / 6) * 6 + 5;
5819 if (
abs(k20) > 4 &&
abs(k20) <= 10)
5821 if (
abs(k20) > 10 &&
abs(k20) <= 15)
5825 if (k30 == 71 || k30 ==
nphi || (k30 >= 1 && k30 <= 10))
5828 kk3 = ((k30 - 11) / 12) * 12 + 11;
5830 }
else if (k1 == 3) {
5835 if (k30 >= 1 && k30 <= 18)
5837 if (k30 >= 19 && k30 <= 36)
5839 if (k30 >= 37 && k30 <= 54)
5841 if (k30 >= 55 && k30 <=
nphi)
5845 int kkk2 = kk2 + 41;
5852 double GetRMSOverNormalizedSignal = -1.;
5853 if (
signal[k1][k2][k3] > 0. &&
calib0[k1][kkk2][kkk3] > 0.) {
5854 GetRMSOverNormalizedSignal =
signal[k1][k2][k3] /
calib0[k1][kkk2][kkk3];
5874 if (
signal[k1][k2][k3] > 0.) {
5877 if (
calib0[k1][kkk2][kkk3] > 0.)
5881 if (
calib0[k1][kkk2][kkk3] > 0.)
5884 float calibtsmax =
calibt[k1][kkk2][kkk3];
5886 float calibtsmean =
caliba[k1][kkk2][kkk3];
5888 float calibwidth =
calibw[k1][kkk2][kkk3];
5890 float calibcap = -100.;
5907 if (calibtsmax > -0.5) {
5914 if (calibtsmean > -0.5) {
5921 if (calibwidth > -0.5) {
5931 h_map_HB->Fill(
double(k2plot),
double(k3), 1.);
5947 if (calibtsmax > -0.5) {
5954 if (calibtsmean > -0.5) {
5961 if (calibwidth > -0.5) {
5971 h_map_HE->Fill(
double(k2plot),
double(k3), 1.);
5987 if (calibtsmax > -0.5) {
5994 if (calibtsmean > -0.5) {
6001 if (calibwidth > -0.5) {
6011 h_map_HO->Fill(
double(k2plot),
double(k3), 1.);
6027 if (calibtsmax > -0.5) {
6034 if (calibtsmean > -0.5) {
6041 if (calibwidth > -0.5) {
6051 h_map_HF->Fill(
double(k2plot),
double(k3), 1.);
6083 std::cout <<
"======================== beignJob START +++++++++++++++++++++++++++" << std::endl;
6321 fs_->
make<TH1F>(
"h_ndepthesperamplitudebins_HE",
" ", 10, 0., 10.);
6331 "h_mapADCAmplsecondpeak0_HE",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
6688 fs_->
make<TH1F>(
"h_AamplitudewithPedSubtr_RADDAM_HEzoom0",
" ", 100, min80, 4000.);
6690 fs_->
make<TH1F>(
"h_AamplitudewithPedSubtr_RADDAM_HEzoom1",
" ", 100, min80, 1000.);
6704 fs_->
make<TH1F>(
"h_sumphiEta16Depth3RADDAM_HED2ALL",
" ", 100, min90, 70. * max90);
6852 float zaR = 10000000.;
6853 float zaR2 = 50000000.;
6887 fs_->
make<TH1F>(
"h_sumADCAmplEtaPhiLs_lscounterM1orbitNum",
" ", 600, 1., 601.);
6895 float bac2 = bac + 1.;
7869 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HB",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7871 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HB",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7873 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HE",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7875 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HE",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7877 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HF",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7879 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HF",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7881 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal_HO",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7883 fs_->
make<TH2F>(
"h_mapGetRMSOverNormalizedSignal0_HO",
" ",
neta, -41., 41.,
nphi, 0.,
bphi);
7932 float spls0 = 10000;
7942 float spls1 = 100000;
7978 std::cout <<
"======================== booking DONE +++++++++++++++++++++++++++" << std::endl;
7987 std::cout <<
"======================== beignJob finish +++++++++++++++++++++++++++" << std::endl;
8000 int mdepth = cell.
depth();
8001 int iphi = cell.iphi() - 1;
8002 int ieta = cell.ieta();
8005 int sub = cell.subdet();
8007 int errorGeneral = 0;
8016 bool anycapid =
true;
8021 int lastcapid = 0, capid = 0;
8022 int ERRORfiber = -10;
8023 int ERRORfiberChan = -10;
8024 int ERRORfiberAndChan = -10;
8025 int repetedcapid = 0;
8027 TSsize = digiItr->size();
8030 for (
int ii = 0;
ii < TSsize;
ii++) {
8031 capid = (*digiItr)[
ii].capid();
8032 bool er = (*digiItr)[
ii].er();
8033 bool dv = (*digiItr)[
ii].dv();
8034 int fiber = (*digiItr)[
ii].fiber();
8035 int fiberChan = (*digiItr)[
ii].fiberChan();
8036 int fiberAndChan = (*digiItr)[
ii].fiberAndChan();
8037 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8040 ERRORfiberChan = fiberChan;
8041 ERRORfiberAndChan = fiberAndChan;
8042 if (capid != lastcapid) {
8055 ERRORfiberChan = fiberChan;
8056 ERRORfiberAndChan = fiberAndChan;
8061 ERRORfiberChan = fiberChan;
8062 ERRORfiberAndChan = fiberAndChan;
8068 if (firstcapid == 0 && !anycapid)
8070 if (firstcapid == 1 && !anycapid)
8072 if (firstcapid == 2 && !anycapid)
8074 if (firstcapid == 3 && !anycapid)
8083 if (!anycapid && anyer)
8085 if (!anycapid && !anydv)
8087 if (!anycapid && anyer && !anydv)
8089 if (anyer && !anydv)
8095 for (
int ii = 0;
ii < TSsize;
ii++) {
8096 double ampldefault =
adc2fC[digiItr->sample(
ii).adc()];
8097 ampl += ampldefault;
8114 if (error1 != 0 || error2 != 0 || error3 != 0) {
8141 if (error1 != 0 || error2 != 0 || error3 != 0) {
8172 int mdepth = hcaldetid.
depth();
8173 int sub = hcaldetid.
subdet();
8177 bool anycapid =
true;
8182 int lastcapid = 0, capid = 0;
8183 int repetedcapid = 0;
8188 capid = qie11df[
ii].capid();
8189 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8191 if (capid != lastcapid) {
8206 int adc = qie11df[
ii].adc();
8212 ampl += ampldefault;
8271 int mdepth = cell.
depth();
8272 int iphi = cell.iphi() - 1;
8273 int ieta = cell.ieta();
8276 int sub = cell.subdet();
8278 std::cout <<
" HF DIGI ?????????????? ERROR mdepth = " << mdepth << std::endl;
8280 int errorGeneral = 0;
8289 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) {
8326 ERRORfiberChan = fiberChan;
8327 ERRORfiberAndChan = fiberAndChan;
8332 ERRORfiberChan = fiberChan;
8333 ERRORfiberAndChan = fiberAndChan;
8337 if (firstcapid == 0 && !anycapid)
8339 if (firstcapid == 1 && !anycapid)
8341 if (firstcapid == 2 && !anycapid)
8343 if (firstcapid == 3 && !anycapid)
8351 if (!anycapid && anyer)
8353 if (!anycapid && !anydv)
8355 if (!anycapid && anyer && !anydv)
8357 if (anyer && !anydv)
8361 for (
int ii = 0;
ii < TSsize;
ii++) {
8362 double ampldefault =
adc2fC[digiItr->sample(
ii).adc()];
8363 ampl += ampldefault;
8377 if (error1 != 0 || error2 != 0 || error3 != 0) {
8407 int mdepth = hcaldetid.
depth();
8408 int sub = hcaldetid.
subdet();
8412 bool anycapid =
true;
8417 int lastcapid = 0, capid = 0;
8418 int repetedcapid = 0;
8423 capid = qie10df[
ii].capid();
8424 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8426 if (capid != lastcapid) {
8441 int adc = qie10df[
ii].adc();
8443 ampl += ampldefault;
8474 int mdepth = cell.
depth();
8475 int iphi = cell.iphi() - 1;
8476 int ieta = cell.ieta();
8479 int sub = cell.subdet();
8480 int errorGeneral = 0;
8489 bool anycapid =
true;
8494 int lastcapid = 0, capid = 0;
8495 int ERRORfiber = -10;
8496 int ERRORfiberChan = -10;
8497 int ERRORfiberAndChan = -10;
8498 int repetedcapid = 0;
8499 for (
int ii = 0;
ii < (*digiItr).size();
ii++) {
8500 capid = (*digiItr)[
ii].capid();
8501 bool er = (*digiItr)[
ii].er();
8502 bool dv = (*digiItr)[
ii].dv();
8503 int fiber = (*digiItr)[
ii].fiber();
8504 int fiberChan = (*digiItr)[
ii].fiberChan();
8505 int fiberAndChan = (*digiItr)[
ii].fiberAndChan();
8506 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8509 ERRORfiberChan = fiberChan;
8510 ERRORfiberAndChan = fiberAndChan;
8511 if (capid != lastcapid) {
8524 ERRORfiberChan = fiberChan;
8525 ERRORfiberAndChan = fiberAndChan;
8530 ERRORfiberChan = fiberChan;
8531 ERRORfiberAndChan = fiberAndChan;
8537 if (firstcapid == 0 && !anycapid)
8539 if (firstcapid == 1 && !anycapid)
8541 if (firstcapid == 2 && !anycapid)
8543 if (firstcapid == 3 && !anycapid)
8552 if (!anycapid && anyer)
8554 if (!anycapid && !anydv)
8556 if (!anycapid && anyer && !anydv)
8558 if (anyer && !anydv)
8563 for (
int ii = 0;
ii < (*digiItr).size();
ii++) {
8564 double ampldefault =
adc2fC[digiItr->sample(
ii).adc()];
8565 ampl += ampldefault;
8581 if (error1 != 0 || error2 != 0 || error3 != 0) {
8590 std::cout <<
" mdepth HO = " << mdepth << std::endl;
8603 double toolwithPedSubtr[100];
8604 double lintoolwithoutPedSubtr[100];
8606 int mdepth = cell.
depth();
8607 int iphi = cell.iphi() - 1;
8608 int ieta0 = cell.ieta();
8612 int sub = cell.subdet();
8624 coder.
adc2fC(*digiItr, toolOriginal);
8625 double pedestalaver9 = 0.;
8626 double pedestalaver4 = 0.;
8627 double pedestal0 = 0.;
8628 double pedestal1 = 0.;
8629 double pedestal2 = 0.;
8630 double pedestal3 = 0.;
8631 double pedestalwaver9 = 0.;
8632 double pedestalwaver4 = 0.;
8633 double pedestalw0 = 0.;
8634 double pedestalw1 = 0.;
8635 double pedestalw2 = 0.;
8636 double pedestalw3 = 0.;
8637 double difpedestal0 = 0.;
8638 double difpedestal1 = 0.;
8639 double difpedestal2 = 0.;
8640 double difpedestal3 = 0.;
8643 double absamplitude = 0.;
8644 double amplitude345 = 0.;
8646 double linamplitudewithoutPedSubtr = 0.;
8649 double max_signal = -100.;
8650 int ts_with_max_signal = -100;
8656 double errorBtype = 0.;
8659 TSsize = digiItr->size();
8660 if ((*digiItr).size() != 10)
8663 for (
int ii = 0;
ii < TSsize;
ii++) {
8665 double ampldefaultwithPedSubtr = 0.;
8666 double linampldefaultwithoutPedSubtr = 0.;
8667 double ampldefault = 0.;
8668 double ampldefault0 = 0.;
8669 double ampldefault1 = 0.;
8670 double ampldefault2 = 0.;
8671 ampldefault0 =
adc2fC[digiItr->sample(
ii).adc()];
8673 ampldefault1 = toolOriginal[
ii];
8674 ampldefault2 = (*digiItr)[
ii].adc();
8676 ampldefault = ampldefault0;
8679 ampldefault = ampldefault1;
8682 ampldefault = ampldefault2;
8684 ampldefaultwithPedSubtr = ampldefault0;
8685 linampldefaultwithoutPedSubtr = ampldefault2;
8687 int capid = ((*digiItr)[
ii]).capid();
8689 double pedestalINI = pedestal00->
getValue(capid);
8691 double pedestalw = pedw->
getSigma(capid, capid);
8692 ampldefaultwithPedSubtr -=
pedestal;
8696 tool[
ii] = ampldefault;
8697 toolwithPedSubtr[
ii] = ampldefaultwithPedSubtr;
8698 lintoolwithoutPedSubtr[
ii] = linampldefaultwithoutPedSubtr;
8701 pedestalwaver9 += pedestalw * pedestalw;
8703 if (capid == 0 &&
c0 == 0) {
8708 pedestalw0 = pedestalw;
8709 pedestalwaver4 += pedestalw * pedestalw;
8710 difpedestal0 =
pedestal - pedestalINI;
8713 if (capid == 1 &&
c1 == 0) {
8718 pedestalw1 = pedestalw;
8719 pedestalwaver4 += pedestalw * pedestalw;
8720 difpedestal1 =
pedestal - pedestalINI;
8722 if (capid == 2 && c2 == 0) {
8727 pedestalw2 = pedestalw;
8728 pedestalwaver4 += pedestalw * pedestalw;
8729 difpedestal2 =
pedestal - pedestalINI;
8731 if (capid == 3 && c3 == 0) {
8736 pedestalw3 = pedestalw;
8737 pedestalwaver4 += pedestalw * pedestalw;
8738 difpedestal3 =
pedestal - pedestalINI;
8741 if (max_signal < ampldefault) {
8742 max_signal = ampldefault;
8743 ts_with_max_signal =
ii;
8758 absamplitude +=
abs(ampldefault);
8760 if (
ii == 3 ||
ii == 4 ||
ii == 5)
8761 amplitude345 += ampldefault;
8765 timew += (
ii + 1) *
abs(ampldefault);
8766 timeww += (
ii + 1) * ampldefault;
8770 if (sub == 1 || sub == 2) {
8775 pedestalaver9 /= TSsize;
8776 pedestalaver4 /= c4;
8777 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
8778 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
8779 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
8780 ampl = tool[ts_with_max_signal];
8781 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
8782 ampl += tool[ts_with_max_signal + 2];
8783 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
8784 ampl += tool[ts_with_max_signal + 1];
8785 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
8786 ampl += tool[ts_with_max_signal - 1];
8789 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
8790 linamplitudewithoutPedSubtr = lintoolwithoutPedSubtr[ts_with_max_signal];
8791 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
8792 linamplitudewithoutPedSubtr += lintoolwithoutPedSubtr[ts_with_max_signal + 2];
8793 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
8794 linamplitudewithoutPedSubtr += lintoolwithoutPedSubtr[ts_with_max_signal + 1];
8795 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
8796 linamplitudewithoutPedSubtr += lintoolwithoutPedSubtr[ts_with_max_signal - 1];
8801 if (ratio < 0. || ratio > 1.02)
8803 double aveamplitude = 0.;
8804 double aveamplitudew = 0.;
8805 if (absamplitude > 0 && timew > 0)
8806 aveamplitude = timew / absamplitude;
8812 bool anycapid =
true;
8817 for (
int ii = 0;
ii < TSsize;
ii++) {
8818 double aaaaaa = (
ii + 1) - aveamplitudew;
8819 double aaaaaa2 = aaaaaa * aaaaaa;
8820 double ampldefault = tool[
ii];
8821 rmsamp += (aaaaaa2 * ampldefault);
8822 capid = ((*digiItr)[
ii]).capid();
8823 bool er = (*digiItr)[
ii].er();
8824 bool dv = (*digiItr)[
ii].dv();
8825 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
8837 if (!anycapid || anyer || !anydv)
8840 double rmsamplitude = 0.;
8843 double aveamplitude1 = aveamplitude - 1;
8849 if (sub == 1 && mdepth == 1)
8851 if (sub == 1 && mdepth == 2)
8853 if (sub == 2 && mdepth == 1)
8855 if (sub == 2 && mdepth == 2)
8857 if (sub == 2 && mdepth == 3)
8861 if (sub == 1 && mdepth == 1)
8863 if (sub == 1 && mdepth == 2)
8865 if (sub == 2 && mdepth == 1)
8867 if (sub == 2 && mdepth == 2)
8869 if (sub == 2 && mdepth == 3)
8873 for (
int ii = 0;
ii < TSsize;
ii++) {
8875 double ampldefault = tool[
ii];
8921 for (
int ii = 0;
ii < TSsize;
ii++) {
8922 double ampldefault = tool[
ii];
8929 for (
int ii = 0;
ii < TSsize;
ii++) {
8930 double ampldefault = tool[
ii];
8942 for (
int ii = 0;
ii < TSsize;
ii++) {
8943 double ampldefault = tool[
ii];
8950 for (
int ii = 0;
ii < TSsize;
ii++) {
8951 double ampldefault = tool[
ii];
8967 double mypedestal = pedestal0;
8968 double mypedestalw = pedestalw0;
9039 for (
int ii = 0;
ii < TSsize;
ii++) {
9072 if (amplitude < ADCAmplHBMin_ || amplitude >
ADCAmplHBMax_) {
9112 if (aveamplitude1 < TSmeanHBMin_ || aveamplitude1 >
TSmeanHBMax_) {
9133 if (ts_with_max_signal < TSpeakHBMin_ || ts_with_max_signal >
TSpeakHBMax_) {
9154 if (rmsamplitude < rmsHBMin_ || rmsamplitude >
rmsHBMax_) {
9213 double mypedestal = pedestal0;
9214 double mypedestalw = pedestalw0;
9285 for (
int ii = 0;
ii < TSsize;
ii++) {
9307 if (amplitude < ADCAmplHEMin_ || amplitude >
ADCAmplHEMax_) {
9368 if (aveamplitude1 < TSmeanHEMin_ || aveamplitude1 >
TSmeanHEMax_) {
9392 if (ts_with_max_signal < TSpeakHEMin_ || ts_with_max_signal >
TSpeakHEMax_) {
9416 if (rmsamplitude < rmsHEMin_ || rmsamplitude >
rmsHEMax_) {
9476 double amplitudewithPedSubtr = 0.;
9479 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
9480 amplitudewithPedSubtr = toolwithPedSubtr[ts_with_max_signal];
9481 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
9482 amplitudewithPedSubtr += toolwithPedSubtr[ts_with_max_signal + 2];
9483 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
9484 amplitudewithPedSubtr += toolwithPedSubtr[ts_with_max_signal + 1];
9485 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
9486 amplitudewithPedSubtr += toolwithPedSubtr[ts_with_max_signal - 1];
9492 if (amplitudewithPedSubtr > 50.) {
9515 int LLLLLL111111 = 0;
9516 if ((mdepth == 1 && fabs(ieta0) > 16 && fabs(ieta0) < 30))
9519 int LLLLLL222222 = 0;
9520 if ((mdepth == 2 && fabs(ieta0) > 16 && fabs(ieta0) < 27) ||
9521 (mdepth == 3 && fabs(ieta0) > 26 && fabs(ieta0) < 29))
9524 if (LLLLLL111111 == 1) {
9529 if (LLLLLL222222 == 1) {
9536 if (mdepth == 3 && fabs(ieta0) == 16) {
9565 int mdepth = hcaldetid.
depth();
9566 int sub = hcaldetid.
subdet();
9569 if (mdepth == 0 || sub > 4)
9594 coder.
adc2fC(qie11df, toolOriginal);
9595 double pedestalaver9 = 0.;
9596 double pedestalaver4 = 0.;
9597 double pedestal0 = 0.;
9598 double pedestal1 = 0.;
9599 double pedestal2 = 0.;
9600 double pedestal3 = 0.;
9601 double pedestalwaver9 = 0.;
9602 double pedestalwaver4 = 0.;
9603 double pedestalw0 = 0.;
9604 double pedestalw1 = 0.;
9605 double pedestalw2 = 0.;
9606 double pedestalw3 = 0.;
9607 double difpedestal0 = 0.;
9608 double difpedestal1 = 0.;
9609 double difpedestal2 = 0.;
9610 double difpedestal3 = 0.;
9613 double amplitude0 = 0.;
9614 double absamplitude = 0.;
9615 double tocampl = 0.;
9617 double amplitude345 = 0.;
9619 double ampl3ts = 0.;
9622 double max_signal = -100.;
9623 int ts_with_max_signal = -100;
9629 double errorBtype = 0.;
9640 for (
int ii = 0;
ii < TSsize;
ii++) {
9641 double ampldefault = 0.;
9642 double tocdefault = 0.;
9643 double ampldefault0 = 0.;
9644 double ampldefault1 = 0.;
9645 double ampldefault2 = 0.;
9651 ampldefault1 = toolOriginal[
ii];
9652 ampldefault2 = qie11df[
ii].adc();
9655 ampldefault = ampldefault0;
9658 ampldefault = ampldefault1;
9661 ampldefault = ampldefault2;
9663 tocdefault = ampldefault;
9665 int capid = (qie11df[
ii]).capid();
9667 double pedestalw = pedw->
getSigma(capid, capid);
9668 double pedestalINI = pedestal00->
getValue(capid);
9672 tool[
ii] = ampldefault;
9674 pedestalwaver9 += pedestalw * pedestalw;
9676 if (capid == 0 &&
c0 == 0) {
9681 pedestalw0 = pedestalw;
9682 pedestalwaver4 += pedestalw * pedestalw;
9683 difpedestal0 =
pedestal - pedestalINI;
9686 if (capid == 1 &&
c1 == 0) {
9691 pedestalw1 = pedestalw;
9692 pedestalwaver4 += pedestalw * pedestalw;
9693 difpedestal1 =
pedestal - pedestalINI;
9695 if (capid == 2 && c2 == 0) {
9700 pedestalw2 = pedestalw;
9701 pedestalwaver4 += pedestalw * pedestalw;
9702 difpedestal2 =
pedestal - pedestalINI;
9704 if (capid == 3 && c3 == 0) {
9709 pedestalw3 = pedestalw;
9710 pedestalwaver4 += pedestalw * pedestalw;
9711 difpedestal3 =
pedestal - pedestalINI;
9714 if (max_signal < ampldefault) {
9715 max_signal = ampldefault;
9716 ts_with_max_signal =
ii;
9719 absamplitude +=
abs(ampldefault);
9720 tocampl += tocdefault;
9722 if (
ii == 1 ||
ii == 2 ||
ii == 3 ||
ii == 4 ||
ii == 5 ||
ii == 6 ||
ii == 7 ||
ii == 8)
9723 amplitude345 += ampldefault;
9728 timew += (
ii + 1) *
abs(ampldefault);
9729 timeww += (
ii + 1) * ampldefault;
9731 if (
ii == 2 && ampldefault > 0.)
9737 pedestalaver9 /= TSsize;
9738 pedestalaver4 /= c4;
9739 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
9740 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
9743 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize) {
9744 ampl = tool[ts_with_max_signal];
9745 ampl3ts = tool[ts_with_max_signal];
9747 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize) {
9748 ampl += tool[ts_with_max_signal - 1];
9749 ampl3ts += tool[ts_with_max_signal - 1];
9751 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize) {
9752 ampl += tool[ts_with_max_signal + 1];
9753 ampl3ts += tool[ts_with_max_signal + 1];
9755 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize) {
9756 ampl += tool[ts_with_max_signal + 2];
9762 double xb = ampl / 40.;
9763 double xc = amplitude345 / 40.;
9764 double xd = ampl3ts / 40.;
9765 double txa = tocampl / 40.;
9767 if (((
ieta == -16 ||
ieta == 15) && mdepth == 4) ||
9768 ((
ieta == -17 ||
ieta == 16) && (mdepth == 2 || mdepth == 3)) ||
9769 ((
ieta == -18 ||
ieta == 17) && mdepth == 5)) {
9770 double c0 = 1.000000;
9771 double b1 = 2.59096e-05;
9772 double a2 = 4.60721e-11;
9773 double corrforxa =
a2 * xa * xa +
b1 * xa +
c0;
9774 double corrforxb =
a2 * xb * xb +
b1 * xb +
c0;
9775 double corrforxc =
a2 * xc * xc +
b1 * xc +
c0;
9776 double corrforxd =
a2 * xd * xd +
b1 * xd +
c0;
9777 double corrfortxa =
a2 * txa * txa +
b1 * txa +
c0;
9782 amplitude345 *= corrforxc;
9783 ampl3ts *= corrforxd;
9784 tocampl *= corrfortxa;
9788 double c0 = 1.000000;
9789 double b1 = 2.71238e-05;
9790 double a2 = 1.32877e-10;
9791 double corrforxa =
a2 * xa * xa +
b1 * xa +
c0;
9792 double corrforxb =
a2 * xb * xb +
b1 * xb +
c0;
9793 double corrforxc =
a2 * xc * xc +
b1 * xc +
c0;
9794 double corrforxd =
a2 * xd * xd +
b1 * xd +
c0;
9795 double corrfortxa =
a2 * txa * txa +
b1 * txa +
c0;
9800 amplitude345 *= corrforxc;
9801 ampl3ts *= corrforxd;
9802 tocampl *= corrfortxa;
9808 if (sub == 1 || sub == 2) {
9819 if (ratio < 0. || ratio > 1.02)
9821 double aveamplitude = 0.;
9822 double aveamplitudew = 0.;
9823 if (absamplitude > 0 && timew > 0)
9824 aveamplitude = timew / absamplitude;
9825 if (amplitude0 > 0 && timeww > 0)
9826 aveamplitudew = timeww / amplitude0;
9830 bool anycapid =
true;
9833 for (
int ii = 0;
ii < TSsize;
ii++) {
9834 double aaaaaa = (
ii + 1) - aveamplitudew;
9835 double aaaaaa2 = aaaaaa * aaaaaa;
9836 double ampldefault = tool[
ii];
9837 rmsamp += (aaaaaa2 * ampldefault);
9838 capid = (qie11df[
ii]).capid();
9839 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
9847 double rmsamplitude = 0.;
9848 if ((amplitude0 > 0 && rmsamp > 0) || (amplitude0 < 0 && rmsamp < 0))
9849 rmsamplitude =
sqrt(rmsamp / amplitude0);
9850 double aveamplitude1 = aveamplitude - 1;
9855 if (sub == 1 && mdepth == 1)
9857 if (sub == 1 && mdepth == 2)
9859 if (sub == 2 && mdepth == 1)
9861 if (sub == 2 && mdepth == 2)
9863 if (sub == 2 && mdepth == 3)
9867 if (sub == 1 && mdepth == 1)
9869 if (sub == 1 && mdepth == 2)
9871 if (sub == 2 && mdepth == 1)
9873 if (sub == 2 && mdepth == 2)
9875 if (sub == 2 && mdepth == 3)
9879 for (
int ii = 0;
ii < TSsize;
ii++) {
9881 double ampldefault = tool[
ii];
9884 if (amplitude0 > 120) {
9927 for (
int ii = 0;
ii < TSsize;
ii++) {
9928 double ampldefault = tool[
ii];
9935 for (
int ii = 0;
ii < TSsize;
ii++) {
9936 double ampldefault = tool[
ii];
9950 for (
int ii = 0;
ii < TSsize;
ii++) {
9951 double ampldefault = tool[
ii];
9958 for (
int ii = 0;
ii < TSsize;
ii++) {
9959 double ampldefault = tool[
ii];
9975 double mypedestal = pedestal0;
9976 double mypedestalw = pedestalw0;
10047 for (
int ii = 0;
ii < TSsize;
ii++) {
10080 if (amplitude < ADCAmplHBMin_ || amplitude >
ADCAmplHBMax_) {
10138 if (aveamplitude1 < TSmeanHBMin_ || aveamplitude1 >
TSmeanHBMax_) {
10166 h_TSmaxA_HB->Fill(
float(ts_with_max_signal), 1.);
10167 if (ts_with_max_signal < TSpeakHBMin_ || ts_with_max_signal >
TSpeakHBMax_) {
10196 if (rmsamplitude < rmsHBMin_ || rmsamplitude >
rmsHBMax_) {
10225 if (ratio < ratioHBMin_ || ratio >
ratioHBMax_) {
10274 if (flagTS2 == 1) {
10293 double mypedestal = pedestal0;
10294 double mypedestalw = pedestalw0;
10381 for (
int ii = 0;
ii < TSsize;
ii++) {
10399 if (amplitude < ADCAmplHEMin_ || amplitude >
ADCAmplHEMax_) {
10447 if (ampl3ts > 70 && ampl3ts < 110) {
10450 }
else if (ampl3ts > 110 && ampl3ts < 150) {
10454 if (ampl > 87 && ampl < 127) {
10457 }
else if (ampl > 127 && ampl < 167) {
10520 if (aveamplitude1 < TSmeanHEMin_ || aveamplitude1 >
TSmeanHEMax_) {
10559 h_TSmaxA_HE->Fill(
float(ts_with_max_signal), 1.);
10560 if (ts_with_max_signal < TSpeakHEMin_ || ts_with_max_signal >
TSpeakHEMax_) {
10600 if (rmsamplitude < rmsHEMin_ || rmsamplitude >
rmsHEMax_) {
10641 if (ratio < ratioHEMin_ || ratio >
ratioHEMax_) {
10711 if (flagTS2 == 1) {
10736 int mdepth = cell.
depth();
10737 int iphi = cell.iphi() - 1;
10738 int ieta = cell.ieta();
10741 int sub = cell.subdet();
10751 coder.
adc2fC(*digiItr, toolOriginal);
10752 double pedestalaver9 = 0.;
10753 double pedestalaver4 = 0.;
10754 double pedestal0 = 0.;
10755 double pedestal1 = 0.;
10756 double pedestal2 = 0.;
10757 double pedestal3 = 0.;
10758 double pedestalwaver9 = 0.;
10759 double pedestalwaver4 = 0.;
10760 double pedestalw0 = 0.;
10761 double pedestalw1 = 0.;
10762 double pedestalw2 = 0.;
10763 double pedestalw3 = 0.;
10765 double absamplitude = 0.;
10768 double timeww = 0.;
10769 double max_signal = -100.;
10770 int ts_with_max_signal = -100;
10776 double errorBtype = 0.;
10778 if ((*digiItr).size() != TSsize)
10780 TSsize = digiItr->size();
10781 for (
int ii = 0;
ii < TSsize;
ii++) {
10783 double ampldefault = 0.;
10784 double ampldefault0 = 0.;
10785 double ampldefault1 = 0.;
10786 double ampldefault2 = 0.;
10787 ampldefault0 =
adc2fC[digiItr->sample(
ii).adc()];
10789 ampldefault1 = toolOriginal[
ii];
10790 ampldefault2 = (*digiItr)[
ii].adc();
10792 ampldefault = ampldefault0;
10795 ampldefault = ampldefault1;
10798 ampldefault = ampldefault2;
10801 int capid = ((*digiItr)[
ii]).capid();
10804 double pedestalw = pedw->
getSigma(capid, capid);
10808 tool[
ii] = ampldefault;
10811 pedestalwaver9 += pedestalw * pedestalw;
10813 if (capid == 0 &&
c0 == 0) {
10818 pedestalw0 = pedestalw;
10819 pedestalwaver4 += pedestalw * pedestalw;
10822 if (capid == 1 &&
c1 == 0) {
10827 pedestalw1 = pedestalw;
10828 pedestalwaver4 += pedestalw * pedestalw;
10830 if (capid == 2 && c2 == 0) {
10835 pedestalw2 = pedestalw;
10836 pedestalwaver4 += pedestalw * pedestalw;
10838 if (capid == 3 && c3 == 0) {
10843 pedestalw3 = pedestalw;
10844 pedestalwaver4 += pedestalw * pedestalw;
10847 if (max_signal < ampldefault) {
10848 max_signal = ampldefault;
10849 ts_with_max_signal =
ii;
10852 absamplitude +=
abs(ampldefault);
10857 timew += (
ii + 1) *
abs(ampldefault);
10858 timeww += (
ii + 1) * ampldefault;
10868 pedestalaver9 /= TSsize;
10869 pedestalaver4 /= c4;
10870 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
10871 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
10874 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
10875 ampl = tool[ts_with_max_signal];
10876 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
10877 ampl += tool[ts_with_max_signal + 2];
10878 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
10879 ampl += tool[ts_with_max_signal + 1];
10880 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
10881 ampl += tool[ts_with_max_signal - 1];
10888 if (ratio < 0. || ratio > 1.02)
10891 double aveamplitude = 0.;
10892 double aveamplitudew = 0.;
10893 if (absamplitude > 0 && timew > 0)
10894 aveamplitude = timew / absamplitude;
10898 double rmsamp = 0.;
10901 bool anycapid =
true;
10902 bool anyer =
false;
10906 for (
int ii = 0;
ii < TSsize;
ii++) {
10907 double aaaaaa = (
ii + 1) - aveamplitudew;
10908 double aaaaaa2 = aaaaaa * aaaaaa;
10909 double ampldefault = tool[
ii];
10910 rmsamp += (aaaaaa2 * ampldefault);
10911 capid = ((*digiItr)[
ii]).capid();
10912 bool er = (*digiItr)[
ii].er();
10913 bool dv = (*digiItr)[
ii].dv();
10914 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
10927 if (!anycapid || anyer || !anydv)
10929 double rmsamplitude = 0.;
10932 double aveamplitude1 = aveamplitude - 1;
10934 if (sub == 4 && mdepth == 1)
10936 if (sub == 4 && mdepth == 2)
10940 if (sub == 4 && mdepth == 1)
10942 if (sub == 4 && mdepth == 2)
10954 for (
int ii = 0;
ii < TSsize;
ii++) {
10955 double ampldefault = tool[
ii];
10962 for (
int ii = 0;
ii < TSsize;
ii++) {
10963 double ampldefault = tool[
ii];
10977 double mypedestal = pedestal0;
10978 double mypedestalw = pedestalw0;
11038 for (
int ii = 0;
ii < TSsize;
ii++) {
11055 if (amplitude < ADCAmplHFMin_ || amplitude >
ADCAmplHFMax_) {
11097 if (aveamplitude1 < TSmeanHFMin_ || aveamplitude1 >
TSmeanHFMax_) {
11117 h_TSmaxA_HF->Fill(
float(ts_with_max_signal), 1.);
11118 if (ts_with_max_signal < TSpeakHFMin_ || ts_with_max_signal >
TSpeakHFMax_) {
11138 if (rmsamplitude < rmsHFMin_ || rmsamplitude >
rmsHFMax_) {
11158 if (ratio < ratioHFMin_ || ratio >
ratioHFMax_) {
11206 int mdepth = hcaldetid.
depth();
11207 int sub = hcaldetid.
subdet();
11222 if (mdepth == 0 || sub != 4)
11242 coder.
adc2fC(qie10df, toolOriginal);
11245 double pedestalaver9 = 0.;
11246 double pedestalaver4 = 0.;
11247 double pedestal0 = 0.;
11248 double pedestal1 = 0.;
11249 double pedestal2 = 0.;
11250 double pedestal3 = 0.;
11251 double pedestalwaver9 = 0.;
11252 double pedestalwaver4 = 0.;
11253 double pedestalw0 = 0.;
11254 double pedestalw1 = 0.;
11255 double pedestalw2 = 0.;
11256 double pedestalw3 = 0.;
11258 double absamplitude = 0.;
11261 double timeww = 0.;
11262 double max_signal = -100.;
11263 int ts_with_max_signal = -100;
11269 double errorBtype = 0.;
11276 for (
int ii = 0;
ii < TSsize;
ii++) {
11277 double ampldefault = 0.;
11278 double ampldefault0 = 0.;
11279 double ampldefault1 = 0.;
11280 double ampldefault2 = 0.;
11283 ampldefault1 = toolOriginal[
ii];
11284 ampldefault2 = qie10df[
ii].adc();
11286 ampldefault = ampldefault0;
11289 ampldefault = ampldefault1;
11292 ampldefault = ampldefault2;
11295 int capid = (qie10df[
ii]).capid();
11297 double pedestalw = pedw->
getSigma(capid, capid);
11302 tool[
ii] = ampldefault;
11305 pedestalwaver9 += pedestalw * pedestalw;
11307 if (capid == 0 &&
c0 == 0) {
11312 pedestalw0 = pedestalw;
11313 pedestalwaver4 += pedestalw * pedestalw;
11316 if (capid == 1 &&
c1 == 0) {
11321 pedestalw1 = pedestalw;
11322 pedestalwaver4 += pedestalw * pedestalw;
11324 if (capid == 2 && c2 == 0) {
11329 pedestalw2 = pedestalw;
11330 pedestalwaver4 += pedestalw * pedestalw;
11332 if (capid == 3 && c3 == 0) {
11337 pedestalw3 = pedestalw;
11338 pedestalwaver4 += pedestalw * pedestalw;
11341 if (max_signal < ampldefault) {
11342 max_signal = ampldefault;
11343 ts_with_max_signal =
ii;
11346 absamplitude +=
abs(ampldefault);
11348 timew += (
ii + 1) *
abs(ampldefault);
11349 timeww += (
ii + 1) * ampldefault;
11350 if (
ii == 1 && ampldefault > 0.)
11361 pedestalaver9 /= TSsize;
11362 pedestalaver4 /= c4;
11363 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
11364 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
11367 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
11368 ampl = tool[ts_with_max_signal];
11369 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
11370 ampl += tool[ts_with_max_signal + 2];
11371 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
11372 ampl += tool[ts_with_max_signal + 1];
11373 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
11374 ampl += tool[ts_with_max_signal - 1];
11380 if (ratio < 0. || ratio > 1.02)
11382 double aveamplitude = 0.;
11383 double aveamplitudew = 0.;
11384 if (absamplitude > 0 && timew > 0)
11385 aveamplitude = timew / absamplitude;
11389 double rmsamp = 0.;
11391 bool anycapid =
true;
11394 for (
int ii = 0;
ii < TSsize;
ii++) {
11395 double aaaaaa = (
ii + 1) - aveamplitudew;
11396 double aaaaaa2 = aaaaaa * aaaaaa;
11397 double ampldefault = tool[
ii];
11398 rmsamp += (aaaaaa2 * ampldefault);
11399 capid = (qie10df[
ii]).capid();
11400 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
11408 double rmsamplitude = 0.;
11411 double aveamplitude1 = aveamplitude - 1;
11416 if (sub == 4 && mdepth == 1)
11418 if (sub == 4 && mdepth == 2)
11422 if (sub == 4 && mdepth == 1)
11424 if (sub == 4 && mdepth == 2)
11436 for (
int ii = 0;
ii < TSsize;
ii++) {
11437 double ampldefault = tool[
ii];
11444 for (
int ii = 0;
ii < TSsize;
ii++) {
11445 double ampldefault = tool[
ii];
11460 double mypedestal = pedestal0;
11461 double mypedestalw = pedestalw0;
11529 for (
int ii = 0;
ii < TSsize;
ii++) {
11546 if (amplitude < ADCAmplHFMin_ || amplitude >
ADCAmplHFMax_) {
11602 if (aveamplitude1 < TSmeanHFMin_ || aveamplitude1 >
TSmeanHFMax_) {
11630 h_TSmaxA_HF->Fill(
float(ts_with_max_signal), 1.);
11631 if (ts_with_max_signal < TSpeakHFMin_ || ts_with_max_signal >
TSpeakHFMax_) {
11660 if (rmsamplitude < rmsHFMin_ || rmsamplitude >
rmsHFMax_) {
11689 if (ratio < ratioHFMin_ || ratio >
ratioHFMax_) {
11738 if (flagTS1 == 1) {
11759 int mdepth = cell.
depth();
11760 int iphi = cell.iphi() - 1;
11761 int ieta = cell.ieta();
11764 int sub = cell.subdet();
11774 coder.
adc2fC(*digiItr, toolOriginal);
11775 double pedestalaver9 = 0.;
11776 double pedestalaver4 = 0.;
11777 double pedestal0 = 0.;
11778 double pedestal1 = 0.;
11779 double pedestal2 = 0.;
11780 double pedestal3 = 0.;
11781 double pedestalwaver9 = 0.;
11782 double pedestalwaver4 = 0.;
11783 double pedestalw0 = 0.;
11784 double pedestalw1 = 0.;
11785 double pedestalw2 = 0.;
11786 double pedestalw3 = 0.;
11788 double absamplitude = 0.;
11791 double timeww = 0.;
11792 double max_signal = -100.;
11793 int ts_with_max_signal = -100;
11799 double errorBtype = 0.;
11801 if ((*digiItr).size() != TSsize)
11803 TSsize = digiItr->size();
11805 for (
int ii = 0;
ii < TSsize;
ii++) {
11806 double ampldefault = 0.;
11807 double ampldefault0 = 0.;
11808 double ampldefault1 = 0.;
11809 double ampldefault2 = 0.;
11810 ampldefault0 =
adc2fC[digiItr->sample(
ii).adc()];
11812 ampldefault1 = toolOriginal[
ii];
11813 ampldefault2 = (*digiItr)[
ii].adc();
11815 ampldefault = ampldefault0;
11818 ampldefault = ampldefault1;
11821 ampldefault = ampldefault2;
11823 int capid = ((*digiItr)[
ii]).capid();
11825 double pedestalw = pedw->
getSigma(capid, capid);
11828 tool[
ii] = ampldefault;
11830 pedestalwaver9 += pedestalw * pedestalw;
11831 if (capid == 0 &&
c0 == 0) {
11836 pedestalw0 = pedestalw;
11837 pedestalwaver4 += pedestalw * pedestalw;
11839 if (capid == 1 &&
c1 == 0) {
11844 pedestalw1 = pedestalw;
11845 pedestalwaver4 += pedestalw * pedestalw;
11847 if (capid == 2 && c2 == 0) {
11852 pedestalw2 = pedestalw;
11853 pedestalwaver4 += pedestalw * pedestalw;
11855 if (capid == 3 && c3 == 0) {
11860 pedestalw3 = pedestalw;
11861 pedestalwaver4 += pedestalw * pedestalw;
11864 if (max_signal < ampldefault) {
11865 max_signal = ampldefault;
11866 ts_with_max_signal =
ii;
11869 absamplitude +=
abs(ampldefault);
11873 timew += (
ii + 1) *
abs(ampldefault);
11874 timeww += (
ii + 1) * ampldefault;
11875 if (ii < 3 && ampldefault > 0.)
11880 pedestalaver9 /= TSsize;
11881 pedestalaver4 /= c4;
11882 pedestalwaver9 =
sqrt(pedestalwaver9 / TSsize);
11883 pedestalwaver4 =
sqrt(pedestalwaver4 / c4);
11884 if (ts_with_max_signal > -1 && ts_with_max_signal < TSsize)
11885 ampl = tool[ts_with_max_signal];
11886 if (ts_with_max_signal + 2 > -1 && ts_with_max_signal + 2 < TSsize)
11887 ampl += tool[ts_with_max_signal + 2];
11888 if (ts_with_max_signal + 1 > -1 && ts_with_max_signal + 1 < TSsize)
11889 ampl += tool[ts_with_max_signal + 1];
11890 if (ts_with_max_signal - 1 > -1 && ts_with_max_signal - 1 < TSsize)
11891 ampl += tool[ts_with_max_signal - 1];
11895 if (ratio < 0. || ratio > 1.04)
11897 double aveamplitude = 0.;
11898 double aveamplitudew = 0.;
11899 if (absamplitude > 0 && timew > 0)
11900 aveamplitude = timew / absamplitude;
11903 double rmsamp = 0.;
11905 bool anycapid =
true;
11906 bool anyer =
false;
11910 for (
int ii = 0;
ii < TSsize;
ii++) {
11911 double aaaaaa = (
ii + 1) - aveamplitudew;
11912 double aaaaaa2 = aaaaaa * aaaaaa;
11913 double ampldefault = tool[
ii];
11914 rmsamp += (aaaaaa2 * ampldefault);
11915 capid = ((*digiItr)[
ii]).capid();
11916 bool er = (*digiItr)[
ii].er();
11917 bool dv = (*digiItr)[
ii].dv();
11918 if (
ii != 0 && ((lastcapid + 1) % 4) != capid) {
11930 if (!anycapid || anyer || !anydv)
11932 double rmsamplitude = 0.;
11935 double aveamplitude1 = aveamplitude - 1;
11937 if (sub == 3 && mdepth == 4)
11941 if (sub == 3 && mdepth == 4)
11952 for (
int ii = 0;
ii < TSsize;
ii++) {
11953 double ampldefault = tool[
ii];
11958 for (
int ii = 0;
ii < TSsize;
ii++) {
11959 double ampldefault = tool[
ii];
11973 double mypedestal = pedestal0;
11974 double mypedestalw = pedestalw0;
12017 for (
int ii = 0;
ii < TSsize;
ii++) {
12033 if (amplitude < ADCAmplHOMin_ || amplitude >
ADCAmplHOMax_) {
12065 if (aveamplitude1 < TSmeanHOMin_ || aveamplitude1 >
TSmeanHOMax_) {
12077 h_TSmaxA_HO->Fill(
float(ts_with_max_signal), 1.);
12078 if (ts_with_max_signal < TSpeakHOMin_ || ts_with_max_signal >
TSpeakHOMax_) {
12092 if (rmsamplitude < rmsHOMin_ || rmsamplitude >
rmsHOMax_) {
12105 if (ratio < ratioHOMin_ || ratio >
ratioHOMax_) {
12122 if (flagTS012 == 1)
12129 if (kdet == 1 || kdet == 2) {
12133 cal_RBX = (kphi + 1) / 4;
12134 cal_RBX = cal_RBX + 18 * (keta + 1) / 2;
12137 cal_RBX = (
int)(kphi / 18) + 1;
12144 cal_RBX = kphi / 12 + 1;
12150 cal_RBX = kphi / 12 + 1 + 6;
12156 cal_RBX = kphi / 6 + 1 + 12;
12162 cal_RBX = kphi / 12 + 1 + 24;
12168 cal_RBX = kphi / 12 + 1 + 30;
12189 for (
int k1 = 0; k1 <
ndepth; k1++) {
12190 for (
int k2 = 0; k2 <
neta; k2++) {
12191 for (
int k3 = 0; k3 <
nphi; k3++) {
12192 int ieta = k2 - 41;
12287 int kkkk2 = (k2 - 1) / 4;
12292 int kkkk3 = (k3) / 4 + 1;
12294 ietaphi = ((kkkk2)-1) *
znphi + (kkkk3);
12297 bbb3 = bbbc / bbb1;
13141 std::cout <<
" ==== Edn of run " << std::endl;
13159 MAPfile <<
"#define LogEleMapdb_h" << std::endl;
13160 MAPfile <<
"#include <algorithm>" << std::endl;
13161 MAPfile <<
"#include <iostream>" << std::endl;
13162 MAPfile <<
"#include <vector>" << std::endl;
13163 MAPfile <<
"#include <string>" << std::endl;
13164 MAPfile <<
"#include <sstream>" << std::endl;
13167 MAPfile <<
"struct Cell {" << std::endl;
13168 MAPfile <<
" std::string subdet;" << std::endl;
13169 MAPfile <<
" int Eta;" << std::endl;
13170 MAPfile <<
" int Phi;" << std::endl;
13171 MAPfile <<
" int Depth;" << std::endl;
13172 MAPfile <<
" std::string RBX;" << std::endl;
13173 MAPfile <<
" int RM;" << std::endl;
13174 MAPfile <<
" int Pixel;" << std::endl;
13175 MAPfile <<
" int RMfiber;" << std::endl;
13176 MAPfile <<
" int FiberCh;" << std::endl;
13177 MAPfile <<
" int QIE;" << std::endl;
13178 MAPfile <<
" int ADC;" << std::endl;
13179 MAPfile <<
" int VMECardID;" << std::endl;
13180 MAPfile <<
" int dccID;" << std::endl;
13181 MAPfile <<
" int Spigot;" << std::endl;
13182 MAPfile <<
" int FiberIndex;" << std::endl;
13183 MAPfile <<
" int HtrSlot;" << std::endl;
13184 MAPfile <<
" int HtrTB;" << std::endl;
13187 MAPfile <<
"// the function check, if \"par\" == \"val\" for this cell" << std::endl;
13188 MAPfile <<
" bool check(const std::string par, const int val) const " << std::endl;
13189 MAPfile <<
" {" << std::endl;
13190 MAPfile <<
" if (par == \"Eta\") return (val == Eta);" << std::endl;
13191 MAPfile <<
" else if (par == \"Phi\") return (val == Phi);" << std::endl;
13192 MAPfile <<
" else if (par == \"Depth\") return (val == Depth);" << std::endl;
13193 MAPfile <<
" else if (par == \"RM\") return (val == RM);" << std::endl;
13194 MAPfile <<
" else if (par == \"Pixel\") return (val == Pixel);" << std::endl;
13195 MAPfile <<
" else if (par == \"RMfiber\") return (val == RMfiber);" << std::endl;
13196 MAPfile <<
" else if (par == \"FiberCh\") return (val == FiberCh);" << std::endl;
13197 MAPfile <<
" else if (par == \"QIE\") return (val == QIE);" << std::endl;
13198 MAPfile <<
" else if (par == \"ADC\") return (val == ADC);" << std::endl;
13199 MAPfile <<
" else if (par == \"VMECardID\") return (val == VMECardID);" << std::endl;
13200 MAPfile <<
" else if (par == \"dccID\") return (val == dccID);" << std::endl;
13201 MAPfile <<
" else if (par == \"Spigot\") return (val == Spigot);" << std::endl;
13202 MAPfile <<
" else if (par == \"FiberIndex\") return (val == FiberIndex);" << std::endl;
13203 MAPfile <<
" else if (par == \"HtrSlot\") return (val == HtrSlot);" << std::endl;
13204 MAPfile <<
" else if (par == \"HtrTB\") return (val == HtrTB);" << std::endl;
13205 MAPfile <<
" else return false;" << std::endl;
13206 MAPfile <<
" }" << std::endl;
13209 MAPfile <<
" bool check(const std::string par, const std::string val) const" << std::endl;
13210 MAPfile <<
" {" << std::endl;
13211 MAPfile <<
" if (par == \"subdet\") return (val == subdet);" << std::endl;
13212 MAPfile <<
" else if (par == \"RBX\") return (val == RBX);" << std::endl;
13213 MAPfile <<
" else return false;" << std::endl;
13214 MAPfile <<
" }" << std::endl;
13216 MAPfile <<
"};" << std::endl;
13219 MAPfile <<
"const Cell AllCells[] = {" << std::endl;
13220 MAPfile <<
"//{ SD, Eta, Phi, Depth, RBX, RM, PIXEL, RMfiber, Fiber Ch., QIE, ADC, VMECrateId, dccid, spigot, " 13221 "fiberIndex, htrSlot, htrTopBottom }" 13234 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13235 << detid->
depth() <<
" ,";
13236 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13237 << lmap_entry.
rmFiber() <<
" , ";
13241 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13243 MAPfile <<
"}," << std::endl;
13256 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() - 1 <<
" , " << detid->
iphi() - 1 <<
" , " 13257 << detid->
depth() <<
" ,";
13258 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13259 << lmap_entry.
rmFiber() <<
" , ";
13263 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13265 MAPfile <<
"}," << std::endl;
13272 for (
int eta = -20;
eta <= -20;
eta++) {
13280 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13281 << detid->
depth() <<
" ,";
13282 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13283 << lmap_entry.
rmFiber() <<
" , ";
13287 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13289 MAPfile <<
"}," << std::endl;
13295 for (
int eta = -19;
eta <= -16;
eta++) {
13303 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13304 << detid->
depth() <<
" ,";
13305 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13306 << lmap_entry.
rmFiber() <<
" , ";
13310 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13312 MAPfile <<
"}," << std::endl;
13317 for (
int eta = -29;
eta <= -16;
eta++) {
13325 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13326 << detid->
depth() <<
" ,";
13327 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13328 << lmap_entry.
rmFiber() <<
" , ";
13332 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13334 MAPfile <<
"}," << std::endl;
13347 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13348 << detid->
depth() <<
" ,";
13349 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13350 << lmap_entry.
rmFiber() <<
" , ";
13354 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13356 MAPfile <<
"}," << std::endl;
13364 for (
int eta = -41;
eta <= -29;
eta++) {
13372 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13373 << detid->
depth() <<
" ,";
13374 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13375 << lmap_entry.
rmFiber() <<
" , ";
13379 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13381 MAPfile <<
"}," << std::endl;
13395 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13396 << detid->
depth() <<
" ,";
13397 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13398 << lmap_entry.
rmFiber() <<
" , ";
13402 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13404 MAPfile <<
"}," << std::endl;
13420 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() <<
" , " << detid->
iphi() - 1 <<
" , " 13421 << detid->
depth() <<
" ,";
13422 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13423 << lmap_entry.
rmFiber() <<
" , ";
13427 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13429 MAPfile <<
"}," << std::endl;
13443 MAPfile <<
" {\"" << subdet <<
"\" , " << detid->
ieta() - 1 <<
" , " << detid->
iphi() - 1 <<
" , " 13444 << detid->
depth() <<
" ,";
13445 MAPfile <<
"\"" << lmap_entry.
rbx() <<
"\" , " << lmap_entry.
rm() <<
" , " << lmap_entry.
pixel() <<
" , " 13446 << lmap_entry.
rmFiber() <<
" , ";
13450 <<
" , " << emap_entry.
fiberIndex() <<
" , ";
13452 MAPfile <<
"}," << std::endl;
13457 MAPfile <<
"};" << std::endl;
13460 MAPfile <<
"// macro for array length calculation" << std::endl;
13461 MAPfile <<
"#define DIM(a) (sizeof(a)/sizeof(a[0]))" << std::endl;
13464 MAPfile <<
"// class for cells array managing" << std::endl;
13465 MAPfile <<
"class CellDB {" << std::endl;
13466 MAPfile <<
"public:" << std::endl;
13467 MAPfile <<
" CellDB()" << std::endl;
13468 MAPfile <<
" : cells(AllCells, AllCells + DIM(AllCells))" << std::endl;
13469 MAPfile <<
"{}" << std::endl;
13472 MAPfile <<
"// return i-th cell" << std::endl;
13473 MAPfile <<
"Cell operator [] (int i) const {return cells[i];}" << std::endl;
13475 MAPfile <<
"// number of cells in database" << std::endl;
13476 MAPfile <<
"int size() const {return cells.size();}" << std::endl;
13479 MAPfile <<
"// select cells for which \"par\" == \"val\"" << std::endl;
13480 MAPfile <<
"template<typename T>" << std::endl;
13481 MAPfile <<
"CellDB find(const std::string par, const T val) const" << std::endl;
13483 MAPfile <<
" std::vector<Cell> s;" << std::endl;
13484 MAPfile <<
" for (size_t i = 0; i < cells.size(); ++i)" << std::endl;
13485 MAPfile <<
" if (cells[i].check(par, val))" << std::endl;
13486 MAPfile <<
" s.push_back(cells[i]);" << std::endl;
13487 MAPfile <<
" return CellDB(s);" << std::endl;
13488 MAPfile <<
"} " << std::endl;
13490 MAPfile <<
"private:" << std::endl;
13491 MAPfile <<
" CellDB( const std::vector<Cell> s)" << std::endl;
13492 MAPfile <<
" : cells(s)" << std::endl;
13493 MAPfile <<
"{}" << std::endl;
13494 MAPfile <<
"std::vector<Cell> cells;" << std::endl;
13495 MAPfile <<
"};" << std::endl;
13497 std::cout <<
"===== Finish writing Channel MAP =====" << std::endl;
13501 double deltaphi = phi1 - phi2;
13503 deltaphi = phi2 - phi1;
13505 if (deltaphi >
M_PI) {
13506 deltaphi = 2. *
M_PI - deltaphi;
13509 double tmp =
sqrt(deltaeta * deltaeta + deltaphi * deltaphi);
13519 if (a1 > 0.5 *
M_PI &&
a2 < 0.)
13521 if (
a2 > 0.5 *
M_PI && a1 < 0.)
13523 double tmp = (a1 * en1 +
a2 * en2) / (en1 + en2);
13536 double tmp =
a2 - a1;
13537 if (a1 *
a2 < 0.) {
13538 if (a1 > 0.5 *
M_PI)
TH2F * h_mapDepth4Pedw1_HO
TH2F * h_mapDepth4AmplE34_HO
double lsdep_estimator2_HEdepth3_
TH1F * h_averSIGNALsumamplitude_HO
TH1F * h_runnbadchannels_depth1_HF
TH2F * h_recNoiseEnergy0_HF1
TH2F * h_2D0Ataildepth2_HB
TH1F * h_sumAmplitudeperLS8
double lsdep_estimator3_HFdepth1_
TH1F * h_nevents_per_LSzoom
TH1F * h_sum0PedestalperLS3
double amplitudechannel2[nsub][ndepth][neta][nphi]
static const std::string kSharedResource
TH1F * h_sumCutADCAmplperLSdepth6HEu
TH1F * h_bcnnbadchannels_depth4_HO
TH2F * h_recNoiseEnergy1_HB1
TH2F * h_2DsumAmplitudeLS5
TH1F * h_RatioOccupancy_HFM
TH2F * h_mapDepth6Amplitude225_HE
TH2F * h_2DsumADCAmplEtaPhiLs3
TH2F * h_mapWidthCalib047_HE
TH1F * h_averSIGNALsumamplitude_HF
TH2F * h_2D0sumADCAmplLS6
TH2F * h_recSignalEnergy1_HF3
bool studyRunDependenceHist_
TH2F * h_mapDepth4TSmeanA225_HE
TH2F * h_mapDepth1TSmaxA225_HE
TH1F * h_sum0ADCAmplperLS3_P1
double const adc2fC_QIE11_shunt6[256]
TH2F * h_mapDepth3TSmaxA225_HE
TH2F * h_FullSignal3D0_HE
TH2F * h_mapDepth1Ped3_HB
TH2F * h_2D0sumTSmeanALS6
TH1F * h_sum0ADCAmplperLS5
TH1F * h_sumADCAmplperLS7u
int lsdep_cut1_peak_HFdepth2_
TH1F * h_shape0_good_channels_HO
TH1F * h_bcnbadrate0_depth4_HO
constexpr DetId detid() const
Get the detector id.
TH2F * h_mapDepth3ADCAmpl225Copy_HF
TH2F * h_FullSignal3D0_HF
TH1F * h_bcnvsamplitude_HE
TH1F * h_RatioOccupancy_HOM
TH2F * h_amplitudechannel0_HE4
TH1F * h_sum0ADCAmplperLS8_M2
TH2F * h_mapDepth1linADCAmpl12_HE
TH1F * h_pedestalaver4_HF
TH2F * h_amplitudechannel1_HB4
TH2F * h_mapDepth2Amplitude_HB
double lsdep_estimator4_HEdepth3_
TH1F * h_runnbadchannels_depth4_HO
TH1F * h_RatioOccupancy_HOP
TH1F * h_runnbadchannels_depth1_HE
TH1F * h_sum0ErrorBperLS6
double calib3[ndepth][neta][nphi]
TH2F * h_mapDepth2Error_HB
constexpr int readoutVMECrateId() const
get the readout VME crate number
TH2F * h_mapDepth2Ped0_HF
TH2F * h_mapDepth2Ampl047_HE
edm::EDGetTokenT< QIE10DigiCollection > tok_qie10_
TH1F * h_runbadrate0_depth3_HE
TH2F * h_amplitudechannel0_HB4
TH1F * h_orbitNumvsamplitude_HF
TH1F * h_pedestalaver4_HE
TH2F * h_mapDepth1Error_HF
TH2F * h_mapDepth5TSmeanA_HE
TH2F * h_mapDepth6Ampl_HE
TH2F * h_recNoiseEnergy1_HE5
T getParameter(std::string const &) const
TH1F * h_bcnvsamplitude0_HB
TH2F * h2_TSnVsAyear2023_HO
int flagToUseDigiCollectionsORNot_
TH2F * h_recSignalEnergy1_HE5
TH2F * h_mapDepth2AmplE34_HE
TH2F * h_recSignalEnergy0_HF1
TH1F * h_ADCAmpl345Zoom_HB
TH2F * h_mapTSmeanCalib_HF
TH1F * h_averNOSIGNALsumamplitude_HO
TH2F * h_mapDepth1Ped2_HE
TH1F * h_runbadrate0_depth4_HO
TH1F * h_bcnbadrate0_depth1_HF
TH2F * h_mapGetRMSOverNormalizedSignal0_HE
TH2F * h_amplitudechannel0_HE3
TH2F * h_mapDepth5TSmaxA_HE
TH2F * h_2DsumPedestalLS8
TH1F * h_ndepthesperamplitudebins_HE
const HcalDbService * conditions
TH1F * h_AmplitudeHEtest6
TH1F * h_pedestalwaver9_HE
TH1F * h_sumTSmeanAperLS7
TH2F * h_2DsumADCAmplLS1_LSselected
const HcalPedestal * getPedestal(const HcalGenericDetId &fId) const
TH1F * h_sumADCAmplperLS8_P2
TH2F * h_mapDepth1Amplitude225_HE
TH1F * h_averNOSIGNALsumamplitude_HE
TH1F * h_sum0ADCAmplperLS6_M2
TH1F * h_pedestalaver9_HB
TH1F * h_sum0TSmaxAperLS7
TH2F * h_mapDepth4ADCAmpl225_HB
double mapRADDAM_HE[ndepth][neta][nphi]
TH2F * h_recSignalEnergy1_HB4
TH2F * h_2DsumADCAmplLS8_LSselected
TH1F * h_sumADCAmplperLS6_P1
TH2F * h_mapDepth3Pedw2_HE
TH1F * h1_TSnVsAyear2023_HF
TH1F * h_AmplitudeHBrest1
TH1F * h_sumamplitude_depth2_HE0
TH2F * h_mapCapCalib047_HE
TH1F * h_sumCutTSmeanAperLS4
double lsdep_estimator2_HBdepth1_
TH2F * h_mapTSmaxCalib_HE
TH2F * h_mapTSmeanCalib047_HO
TH2F * h_amplitudechannel1_HF3
TH2F * h_mapDepth3Error_HB
TH1F * h_Amplitude_notCapIdErrors_HO4
TH2F * h_mapGetRMSOverNormalizedSignal_HF
TH1F * h_sumCutAmplitudeperLS1
TH2F * h_amplitudechannel1_HB2
double calib0[ndepth][neta][nphi]
TH2F * h_amplitudechannel0_HB3
TH1F * h_sumamplitudechannel_HB
edm::EDGetTokenT< HBHERecHitCollection > tok_hbheNoise_
TH2F * h_mapDepth4Pedw0_HO
double lsdep_estimator2_HOdepth4_
float calibt[ndepth][neta][nphi]
TH1F * h_sumamplitude_depth1_HE
TH1F * h_sumCutTSmaxAperLS1
TH1F * h_sumADCAmplperLSdepth4HBu
T const & getData(const ESGetToken< T, R > &iToken) const noexcept(false)
TH2F * h_mapDepth1Pedw2_HE
TH1F * h_sumCutAmplperLS7
TH1F * h_sumCutADCAmplperLSdepth4HEu
TH2F * h_mapDepth3pedestalw_HB
TH2F * h_recNoiseEnergy2_HF1
TH1F * h_pedestalwaver4_HB
TH1F * h_ADCAmpl345Zoom1_HE
TH2F * h_recSignalEnergy0_HF4
double lsdep_estimator3_HBdepth2_
TH2F * h_recNoiseEnergy1_HB4
CMTRawAnalyzer(const edm::ParameterSet &)
TH2F * h_mapDepth3Ped3_HE
TH2F * h_recSignalEnergy0_HE6
TH2F * h_amplitudechannel1_HF4
TH1F * h_averSIGNALoccupancy_HB
TH2F * h_mapDepth2TSmaxA225_HE
TH2F * h_mapDepth2Amplitude225_HF
TH1F * h_sumTSmeanAperLS4
TH2F * h_mapDepth4Amplitude_HF
TH2F * h_mapDepth2Error_HF
void fillDigiAmplitudeHO(HODigiCollection::const_iterator &digiItr)
TH1F * h_sumADCAmplperLS8
TH1F * h_sumADCAmplEtaPhiLs
TH2F * h_mapDepth4ADCAmpl12_HF
TH1F * h_runbadrate_depth2_HE
TH1F * h_Eta16Depth3RADDAM_HED2P
TH2F * h_amplitudechannel0_HE7
TH2F * h_mapDepth2Pedw3_HF
TH2F * h_mapDepth4ADCAmpl225Copy_HB
double averSIGNALsumamplitude_HO
TH2F * h_2DsumADCAmplLSdepth3HFu
void fillDigiErrors(HBHEDigiCollection::const_iterator &digiItr)
TH2F * h_recNoiseEnergy1_HE6
double lsdep_estimator3_HBdepth1_
TH1F * h_numberofhitsHEtest
double sum0Estimator[nsub][ndepth][neta][nphi]
TH2F * h_mapAiForLS_bad_HE
TH1F * h_AmplitudeHEtest1
double lsdep_estimator1_HFdepth3_
TH2F * h_2DsumADCAmplLSdepth5HEu
TH1F * h_sumamplitude_depth4_HO
TH1F * h_pedestalaver9_HO
TH2F * h_2D0sumTSmeanALS7
TH2F * h_2DsumAmplitudeLS1
TH1F * h_sum0ADCAmplperLS7u
TH1F * h_runnbadchannels_depth3_HE
TH2F * h_amplitudechannel2_HF1
TH2F * h_mapDepth1ADCAmpl12_HF
TH1F * h_AmplitudeHBrest6
TH2F * h_mapDepth4AmplE34_HE
TH2F * h_mapDepth3ADCAmpl12_HF
TH1F * h_sigLayer1RADDAM6_HED2
double mapRADDAM_HED2[ndepth][neta]
TH2F * h_maprphinorm0_HE4
TH2F * h_mapDepth1Pedw1_HE
TH2F * h_2D0sumPedestalLS1
TH2F * h_mapAiForLS_good_HE
TH2F * h_2DsumADCAmplLSdepth4HBu
TH2F * h_mapDepth2ADCAmpl225_HE
TH1F * h_bcnnbadchannels_depth2_HF
TH2F * h_mapDepth1Ped1_HE
TH2F * h_mapDepth2ADCAmpl225_HB
TH2F * h_2DsumPedestalLS4
TH2F * h_recSignalEnergy1_HE2
TH2F * h_2D0sumADCAmplLS1
TH1F * h_A_Depth2RADDAM_HE
float getValue(int fCapId) const
get value for capId = 0..3
TH1F * h_sumPedestalperLS4
TH2F * h_mapADCCalib047_HE
TH2F * h_mapDepth4Pedw3_HO
TH1F * h_shape_good_channels_HE
TH1F * h_sumADCAmplperLS1_P1
TH1F * h_sum0AmplitudeperLS1
TH2F * h_amplitudechannel2_HE6
int calibcapiderror[ndepth][neta][nphi]
TH2F * h_mapDepth1Ped3_HE
TH1F * h_sumphiEta16Depth3RADDAM_HED2
TH2F * h_2D0sumPedestalLS6
TH2F * h_mapDepth1Error_HB
TH2F * h_mapDepth2Ped1_HE
bool studyTSmaxShapeHist_
int flagIterativeMethodCalibrationGroupDigi_
TH2F * h_mapRatioCalib_HO
TH1F * h_sumCutAmplperLS5
double maprphinorm[nsub][ndepth][neta][nphi]
TH1F * h_gsmdifferencefit4_HE
TH2F * h_recSignalEnergy2_HB2
TH2F * h_mapDepth6ADCAmpl12_HE
TH2F * h_mapADCAmpl11firstpeak_HE
TH1F * h_eventoccupancy_HF
TH2F * h_2DsumAmplitudeLS6
int mapRADDAM0_HE[ndepth][neta][nphi]
TH1F * h_sum0ErrorBperLS4
TH1F * h_runbadrate_depth1_HB
double lsdep_estimator3_HOdepth4_
TH2F * h_recNoiseEnergy2_HE5
TH2F * h_recNoiseEnergy2_HE6
TH1F * h_gsmdifferencefit5_HE
TH2F * h_2DAtaildepth2_HB
TH2F * h_mapDepth2Amplitude225_HE
int lsdep_cut3_max_HEdepth2_
TH1F * h_Amplitude_notCapIdErrors_HB1
TH1F * h_runnbadchannelsC_depth1_HE
TH2F * h_recNoiseEnergy0_HF3
double amplitudechannel[nsub][ndepth][neta][nphi]
double averSIGNALoccupancy_HO
TH1F * h_Amplitude_forCapIdErrors_HB1
TH2F * h_2DsumPedestalLS7
TH1F * h_runbadrateC_depth2_HB
TH2F * h_mapDepth3Error_HE
TH2F * h_mapCapCalib047_HO
double lsdep_estimator4_HFdepth2_
TH2F * h_mapCapCalib047_HB
TH2F * h_mapDepth2Ampl_HE
TH1F * h_pedestalwaver4_HO
TH2F * h_mapDepth3Amplitude225_HB
TH2F * h_mapDepth4Amplitude225_HB
TH2F * h_recSignalEnergy1_HF2
TH2F * h_amplitudechannel0_HE2
TH2F * h_mapDepth2Pedw3_HE
int lsdep_cut3_max_HEdepth1_
TH2F * h_mapDepth4Ped1_HO
TH1F * h_runnbadchannelsC_depth1_HF
TH2F * h_mapDepth4Ampl047_HF
TH1F * h_sumamplitude_depth4_HO0
TH1F * h_shape0_bad_channels_HB
TH2F * h2_TSnVsAyear2023_HE
TH1F * h_RatioOccupancy_HEP
TH2F * h_mapWidthCalib_HO
TH1F * h_Eta16Depth3RADDAM_HED2ALL
TH2F * h_amplitudechannel1_HE5
TH2F * h_mapDepth1TSmaxA_HF
TH2F * h_recNoiseEnergy0_HE6
TH1F * h_sigLayer1RADDAM6_HE
TH1F * h_shape0_good_channels_HB
double calibrTSmeanHOMax_
TH2F * h_2D0sumPedestalLS4
TH2F * h_2DsumPedestalLS1
TH1F * h_sigLayer1RADDAM0_HE
TH1F * h_sum0ErrorBperLS5
TH2F * h_mapDepth7ADCAmpl_HE
const DetId lookup(HcalElectronicsId fId) const
lookup the logical detid associated with the given electronics id
TH2F * h_2DsumADCAmplEtaPhiLs0
double lsdep_estimator2_HEdepth1_
TH2F * h_mapDepth4ADCAmpl_HE
TH1F * h_sumamplitude_depth2_HB
TH1F * h_NphiForEta16Depth3RADDAM_HED2ALL
TH2F * h_2DsumAmplitudeLS4
TH2F * h_recNoiseEnergy1_HE3
TH2F * h_mapDepth4ADCAmpl12_HE
double phi12(double phi1, double en1, double phi2, double en2)
TH1F * h_energyhitNoise_HE
double lsdep_estimator1_HBdepth1_
bool studyRatioShapeHist_
edm::EDGetTokenT< HBHEDigiCollection > tok_hbhe_
TH2F * h_mapDepth4Error_HB
std::vector< T >::const_iterator const_iterator
TH2F * h_2DsumPedestalLS3
TH2F * h_FullSignal3D0_HO
TH1F * h_sumTSmeanAperLS3
TH2F * h_2D0Ataildepth1_HB
TH2F * h_mapDepth4Ampl_HF
int lsdep_cut3_max_HEdepth3_
TH1F * h_sum0ADCAmplperLSdepth5HEu
TH2F * h_mapDepth3ADCAmpl_HF
TH1F * h_gsmdifferencefit3_HE
const HcalQIEShape * shape
TH1F * h_sumamplitude_depth2_HE
double lsdep_estimator1_HFdepth4_
TH2F * h_mapDepth2linADCAmpl12_HE
TH2F * h_mapGetRMSOverNormalizedSignal0_HO
TH2F * h_recNoiseEnergy1_HF4
TH2F * h_mapDepth1Ampl047_HB
TH2F * h_mapDepth2Amplitude225_HB
TH2F * h_mapTSmeanCalib047_HB
TH2F * h_recNoiseEnergy0_HB4
TH2F * h_mapDepth4Amplitude_HO
TH2F * h_recNoiseEnergy0_HE2
edm::EDGetTokenT< HFRecHitCollection > tok_hfSignal_
TH2F * h_mapDepth1AmplE34_HF
TH2F * h_mapDepth4pedestal_HB
double recSignalEnergy2[nsub][ndepth][neta][nphi]
TH2F * h_recSignalEnergy1_HB2
TH2F * h_amplitudechannel2_HE7
TH2F * h_mapDepth3Ampl047_HE
TH2F * h_mapDepth1ADCAmpl225_HE
TH2F * h_mapDepth7ADCAmpl225_HE
double averSIGNALoccupancy_HB
TH2F * h_mapDepth1Ped0_HB
TH2F * h_2D0sumAmplitudeLS2
TH2F * h_mapDepth1Ped2_HB
TH1F * h_nbadchannels_depth3_HE
TH1F * h_sum0ADCAmplperLS6_P2
TH1F * h_Amplitude_forCapIdErrors_HE3
TH2F * h_mapDepth4TSmeanA225_HB
TH1F * h_sumamplitude_depth2_HB0
TH1F * h_sumAmplperLS1_LSselected
TH1F * h_corrforxaMAIN0_HE
TH2F * h_mapDepth2ADCAmpl_HF
TH1F * h_sumamplitude_depth1_HF1
int lsdep_cut1_peak_HEdepth3_
TH2F * h_mapGetRMSOverNormalizedSignal_HE
double lsdep_estimator1_HEdepth4_
TH1F * h_Amplitude_forCapIdErrors_HF1
TH2F * h_2DsumAmplitudeLS3
TH1F * h_sum0PedestalperLS1
TH1F * h_pedestalaver4_HO
TH2F * h_mapDepth2TSmaxA_HB
TH2F * h_FullSignal3D0_HB
TH1F * h_sum0ADCAmplperLSdepth3HBu
TH2F * h_mapDepth3Amplitude_HE
TH1F * h_ADCAmpl345Zoom_HE
TH2F * h_mapDepth5TSmeanA225_HE
TH2F * h_recNoiseEnergy0_HE7
TH2F * h_recNoiseEnergy2_HB1
TH2F * h_recNoiseEnergy0_HB2
TH2F * h_mapTSmaxCalib047_HB
TH1F * h_Amplitude_notCapIdErrors_HF2
TH1F * h_sumADCAmplperLS8_P1
TH2F * h_mapDepth1AmplE34_HB
TH2F * h_recSignalEnergy2_HE7
TH1F * h_sum0TSmeanAperLS7
TH2F * h_mapTSmaxCalib_HO
double lsdep_estimator2_HEdepth2_
TH2F * h_mapDepth4Error_HO
double sumEstimator3[nsub][ndepth][neta][nphi]
const HcalElectronicsMap * getHcalMapping() const
TH2F * h_recSignalEnergy0_HE2
TH2F * h_mapDepth2ADCAmpl_HB
TH2F * h_mapDepth1TSmaxA225_HF
TH2F * h_recNoiseEnergy1_HB3
double lsdep_estimator4_HBdepth2_
TH2F * h_mapDepth7TSmeanA225_HE
TH2F * h_mapRatioCalib_HB
double lsdep_estimator1_HEdepth3_
TH2F * h_mapDepth1ADCAmpl225_HF
TH1F * h_Amplitude_forCapIdErrors_HE1
int lsdep_cut1_peak_HFdepth1_
TH2F * h_mapDepth4Amplitude225_HF
TH2F * h_mapDepth2pedestal_HB
double lsdep_estimator3_HEdepth2_
TH2F * h_mapGetRMSOverNormalizedSignal_HO
TH2F * h_mapDepth1Pedw3_HE
bool studyTSmeanShapeHist_
void fillDigiAmplitudeQIE11(QIE11DataFrame qie11df)
TH2F * h_2D0sumADCAmplLSdepth4HFu
double calibrTSmeanHFMin_
TH1F * h_sumamplitude_depth1_HB
constexpr int htrSlot() const
get the htr slot
TH1F * h_corrforxaADDI_HE
TH1F * h_sumPedestalperLS1
TH1F * h_sumAmplitudeperLS1
void fillDigiErrorsHF(HFDigiCollection::const_iterator &digiItr)
TH2F * h_mapDepth3ADCAmpl225_HB
TH1F * h_nbadchannels_depth2_HB
TH2F * h_2DAtaildepth1_HB
double lsdep_estimator5_HEdepth1_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbheSignal_
TH1F * h_shape_good_channels_HF
TH2F * h_mapDepth2Pedw1_HE
TH1F * h_sum0TSmaxAperLS4
TH2F * h_amplitudechannel2_HF3
double lsdep_estimator1_HFdepth2_
TH2F * h_mapDepth4TSmaxA225_HF
TH2F * h_recNoiseEnergy1_HE1
int lsdep_cut3_max_HOdepth4_
TH1F * h_sigLayer2RADDAM6_HE
TH2F * h_recSignalEnergy2_HE6
TH2F * h_mapDepth4TSmeanA_HB
TH2F * h_mapDepth1RADDAM0_HE
TH2F * h_2DsumADCAmplLSdepth3HBu
TH2F * h_2D0sumADCAmplLSdepth6HEu
TH2F * h_recSignalEnergy1_HE1
TH1F * h_sum0ADCAmplperLS8_P2
TH1F * h_nbadchannels_depth1_HB
TH2F * h_recNoiseEnergy1_HB2
TH2F * h_mapDepth3Ampl_HB
const HcalPedestalWidth * getPedestalWidth(const HcalGenericDetId &fId) const
TH1F * h_sum0ADCAmplperLS3_P2
TH2F * h_mapDepth4ADCAmpl_HB
TH1F * h_energyhitSignal_HB
TH1F * h_runnbadchannelsC_depth2_HE
TH1F * h_bcnvsamplitude_HB
TH2F * h_amplitudechannel0_HF2
TH2F * h_mapDepth5ADCAmpl225_HE
TH1F * h_sumADCAmplperLS3_M2
TH1F * h_sum0ADCAmplperLS1_P2
TH2F * h_mapDepth1TSmeanA225_HE
TH2F * h_recSignalEnergy1_HE3
int flagfitshunt1pedorledlowintensity_
TH2F * h_mapDepth2Pedw3_HB
TH1F * h_sumamplitude_depth3_HE
TH1F * h_sumCutTSmeanAperLS8
TH2F * h_mapDepth1TSmeanA_HB
TH2F * h_mapDepth4pedestal_HE
TH1F * h_Amplitude_forCapIdErrors_HB2
TH2F * h_mapDepth4Ampl047_HB
TH1F * h_sumADCAmplperLS2
TH2F * h_mapDepth1ADCAmpl12_HE
TH2F * h_2DsumADCAmplEtaPhiLs1
TH2F * h_2DsumADCAmplLS4_LSselected
TH1F * h_nbadchannels_depth2_HE
TH1F * h_sumADCAmplLS1copy1
TH2F * h_mapDepth1Ampl_HE
double dPhiWsign(double phi1, double phi2)
TH1F * h_sum0ADCAmplperLS3_M1
TH2F * h_recNoiseEnergy0_HB3
double tocamplchannel[nsub][ndepth][neta][nphi]
double calib2[ndepth][neta][nphi]
TH2F * h_mapDepth4ADCAmpl225_HF
TH1F * h1_TSnVsAyear20230_HB
TH1F * h_sum0ErrorBperLS3
TH1F * h_sum0PedestalperLS8
TH1F * h_sumADCAmplperLS6u
TH1F * h_runbadrate_depth4_HO
double lsdep_estimator2_HFdepth2_
TH2F * h_mapDepth2Ped2_HB
TH2F * h_mapDepth1Ampl_HB
TH2F * h_maprphinorm2_HE7
double averNOSIGNALsumamplitude_HO
TH2F * h_mapDepth4Ampl047_HO
TH2F * h_2DsumPedestalLS5
TH1F * h_Amplitude_notCapIdErrors_HE3
TH2F * h_mapDepth2Pedw1_HB
constexpr DetId detid() const
Get the detector id.
TH2F * h_mapDepth2Ampl047_HB
TH2F * h_mapDepth1TSmaxA_HB
TH2F * h_mapDepth2Ampl_HF
TH2F * h_mapDepth2pedestalw_HF
TH2F * h_mapADCAmpl11firstpeak0_HE
TH1F * h_sumCutAmplitudeperLS7
edm::Service< TFileService > fs_
TH2F * h_2D0sumTSmeanALS8
TH2F * h_mapCapCalib047_HF
TH1F * h_sum0AmplitudeperLS2
TH1F * h_bcnvsamplitude0_HF
double lsdep_estimator1_HEdepth5_
double calibrTSmeanHEMax_
TH2F * h_mapDepth7ADCAmpl225Copy_HE
TH1F * h_sumCutTSmaxAperLS2
constexpr float energy() const
TH1F * h_sumCutAmplperLS4
TH1F * h_sumCutTSmeanAperLS3
TH2F * h_mapDepth2ADCAmplSiPM_HE
TH2F * h_mapDepth3pedestal_HB
TH2F * h_amplitudechannel1_HB1
TH2F * h_mapDepth4ADCAmpl225Copy_HO
TH2F * h_mapenophinorm_HE2
TH2F * h_2DsumADCAmplLSdepth4HEu
TH1F * h_sum0PedestalperLS6
TH2F * h_mapDepth4pedestalw_HF
TH1F * h_sum0ADCAmplperLS3_M2
TH2F * h_amplitudechannel0_HB1
int lsdep_cut1_peak_HBdepth1_
TH2F * h_recSignalEnergy2_HF2
TH1F * h_sumCutTSmeanAperLS7
TH2F * h_recSignalEnergy0_HF3
TH2F * h_amplitudechannel1_HE4
TH2F * h_recSignalEnergy1_HE7
TH1F * h_sumADCAmplperLS6_P2
TH2F * h_mapADCAmpl11secondpeak_HE
TH1F * h_nbadchannels_depth4_HO
TH2F * h_mapWidthCalib047_HB
TH2F * h_amplitudechannel2_HE5
TH1F * h_sumPedestalperLS3
TH1F * h_sumADCAmplLS1copy2
TH1F * h_sumCutADCAmplperLS5
TH1F * h_runbadrate0_depth2_HB
int lsdep_cut1_peak_HEdepth2_
double sumEstimator5[nsub][ndepth][neta][nphi]
TH1F * h_sumADCAmplperLSdepth6HEu
TH2F * h_recSignalEnergy2_HE3
TH1F * h_runbadrate_depth1_HF
T getUntrackedParameter(std::string const &, T const &) const
TH1F * h_sum0ErrorBperLS2
TH2F * h_mapDepth4Ampl_HE
TH2F * h_mapDepth4TSmaxA_HB
TH1F * h_sum0ADCAmplperLS6_M1
TH2F * h_mapDepth4Error_HE
TH2F * h_2D0sumADCAmplLS3
TH1F * h_sumCutADCAmplperLSdepth7HEu
TH1F * h_pedestalwaver9_HF
TH1F * h_sumTSmeanAperLS8
TH2F * h_recSignalEnergy1_HE6
double lsdep_estimator2_HFdepth1_
edm::EDGetTokenT< HFRecHitCollection > tok_hfNoise_
TH2F * h_mapDepth1TSmeanA225_HF
TH1F * h1_TSnVsAyear20230_HF
TH2F * h_mapDepth5AmplE34_HE
double lsdep_estimator1_HBdepth3_
TH2F * h_mapAiForLS_good_HF
TH2F * h_map0AiForLS_bad_HB
TH2F * h_mapWidthCalib047_HO
void fillDigiErrorsQIE11(QIE11DataFrame qie11df)
TH1F * h_sigLayer1RADDAM5_HED2
TH2F * h_mapDepth1ADCAmpl225_HB
TH2F * h_mapDepth2Error_HE
TH1F * h_sumamplitude_depth2_HF
TH2F * h_amplitudechannel2_HB4
TH2F * h_mapenophinorm2_HE1
TH2F * h_mapDepth2TSmeanA225_HE
TH2F * h_mapDepth4TSmaxA_HO
TH1F * h_sumCutAmplperLS6
TH2F * h_2D0sumTSmeanALS5
TH1F * h_shape0_bad_channels_HF
TH2F * h_amplitudechannel0_HE6
TH1F * h_sum0TSmaxAperLS3
TH2F * h_mapDepth1Pedw3_HB
TH1F * h_pedestalwaver9_HO
TH2F * h_mapDepth1TSmaxA_HE
TH2F * h_mapDepth3ADCAmpl12_HE
constexpr HcalSubdetector subdet() const
get the subdetector
void adc2fC(const HBHEDataFrame &df, CaloSamples &lf) const override
TH2F * h_maprphinorm0_HE2
TH2F * h_mapDepth4AmplE34_HB
TH2F * h_mapDepth2Ped0_HE
TH2F * h_mapDepth1pedestal_HB
TH1F * h_Amplitude_forCapIdErrors_HE2
TH1F * h_bcnvsamplitude_HO
TH2F * h_mapDepth3TSmaxA225_HB
TH1F * h_sumamplitudechannel_HE
double averNOSIGNALoccupancy_HB
TH2F * h_mapDepth4TSmeanA_HF
TH2F * h_mapADCAmplfirstpeak_HE
TH1F * h_sum0TSmaxAperLS2
TH2F * h_mapDepth4ADCAmpl225Copy_HF
TH1F * h_sum0PedestalperLS4
TH2F * h_recNoiseEnergy0_HE5
TH2F * h_mapDepth7TSmeanA_HE
TH2F * h_mapDepth4Ampl_HO
TH1F * h_nbadchannels_depth2_HF
TH1F * h_sum0AmplitudeperLS5
TH2F * h_mapDepth3ADCAmpl225_HF
TH2F * h_amplitudechannel0_HF4
TH1F * h_orbitNumvsamplitude0_HE
TH1F * h_sigLayer2RADDAM5_HED2
TH2F * h_recSignalEnergy0_HB2
TH1F * h_sumADCAmplperLS4
TH1F * h_eventoccupancy_HB
TH2F * h_2D0sumAmplitudeLS1
TH1F * h_sumADCAmplperLS3_M1
TH2F * h_recSignalEnergy2_HE4
TH1F * h_shape0_good_channels_HE
TH2F * h_mapADCCalib047_HO
TH1F * h_eventoccupancy_HE
double lsdep_estimator1_HOdepth4_
TH2F * h_mapDepth2Ped2_HE
TH1F * h_sumADCAmplEtaPhiLs_ietaphi
double lsdep_estimator4_HOdepth4_
TH1F * h_mapDepth3RADDAM16_HE
TH1F * h_sum0ADCAmplperLSdepth4HBu
TH1F * h_shape_bad_channels_HO
TH1F * h_Amplitude_forCapIdErrors_HF2
TH1F * h_RatioOccupancy_HBP
double lsdep_estimator4_HBdepth1_
TH1F * h_sumCutADCAmplperLS6
TH2F * h_2D0sumPedestalLS8
TH2F * h_recSignalEnergy0_HB1
TH1F * h_runbadrateC_depth2_HE
TH2F * h_recSignalEnergy1_HE4
TH1F * h_sumADCAmplperLS8_M2
TH1F * h_sumamplitudechannel_HF
TH2F * h_mapDepth2ADCAmpl12SiPM_HE
int lsdep_cut3_max_HFdepth1_
TH2F * h_mapDepth1pedestalw_HB
TH2F * h_mapDepth2RADDAM_HE
TH2F * h_mapDepth6pedestal_HE
TH2F * h_mapDepth3TSmeanA225_HF
TH2F * h_mapDepth4ADCAmpl12_HB
TH2F * h_2D0sumPedestalLS7
TH2F * h_mapDepth1Amplitude_HF
TH1F * h_sumADCAmplLS1copy5
TH2F * h_mapDepth2Ped1_HB
TH2F * h_2D0sumADCAmplLSdepth3HFu
const HcalTopology * topo
TH2F * h_amplitudechannel2_HF4
TH1F * h_runbadrateC_depth2_HF
double forallestimators_amplitude_bigger_
TH1F * h_sum0AmplitudeperLS8
TH2F * h_mapDepth5Ampl047_HE
TH1F * h_sum0AmplitudeperLS3
TH1F * h_ADCAmpl345Zoom1_HB
int lsdep_cut3_max_HBdepth2_
TH1F * h_sumADCAmplperLS7
double lsdep_estimator5_HFdepth2_
TH1F * h_runbadrate_depth2_HF
TH1F * h_sumTSmeanAperLS1
TH2F * h_mapDepth3Amplitude_HB
TH1F * h_energyhitNoise_HB
TH1F * h_sumADCAmplperLS1_M1
TH1F * h_totalAmplitudeHEperEvent
TH2F * h_mapDepth3RADDAM_HE
TH2F * h_mapTSmaxCalib_HF
TH2F * h_mapDepth3ADCAmpl225Copy_HE
TH2F * h_2DsumPedestalLS2
TH1F * h_sumADCAmplperLSdepth3HBu
TH2F * h_mapDepth2Ped3_HF
const HcalRespCorr * getHcalRespCorr(const HcalGenericDetId &fId) const
TH1F * h_energyhitSignal_HF
TH1F * h_sumADCAmplperLS6
TH2F * h_mapDepth5pedestalw_HE
TH1F * h_runbadrateC_depth1_HF
int lsdep_cut3_max_HFdepth2_
TH2F * h_mapWidthCalib047_HF
TH2F * h_recSignalEnergy2_HE2
TH2F * h_amplitudechannel0_HE1
TH1F * h_runnbadchannelsC_depth2_HF
int getRBX(int &i, int &j, int &k)
TH1F * h_sumCutAmplitudeperLS8
TH2F * h_mapAiForLS_bad_HF
TH2F * h_mapDepth7Ampl047_HE
TH1F * h_sumCutTSmeanAperLS1
TH1F * h_sum0ADCAmplperLS1
TH1F * h_sumamplitude_depth2_HE1
TH1F * h_sumADCAmplperLS8_M1
TH2F * h_maprphinorm2_HE5
TH2F * h_mapDepth2TSmaxA_HE
TH1F * h_bcnnbadchannels_depth1_HB
TH1F * h_sumADCAmplperLS1
constexpr int fiberIndex() const
get the fiber index. For VME 1-8 (which of eight fibers carried by a spigot), for uTCA fibers are zer...
constexpr int ieta() const
get the cell ieta
TH2F * h_mapDepth2Amplitude_HF
TH2F * h_mapDepth1Ampl_HF
TH2F * h_mapADCAmpl12secondpeak_HE
TH2F * h_mapDepth1TSmeanA225_HB
TH1F * h_sum0ErrorBperLS1
TH2F * h_recNoiseEnergy0_HE4
TH2F * h_amplitudechannel2_HB3
TH1F * h_runnbadchannels_depth2_HB
TH2F * h_mapDepth3Amplitude_HF
TH1F * h1_TSnVsAyear2023_HB
TH1F * h_sumTSmeanAperLS2
const edm::ESGetToken< HcalTopology, HcalRecNumberingRecord > tokTopo_
TH2F * h_mapDepth4Amplitude225_HO
TH1F * h_sumamplitude_depth1_HB0
TH2F * h_2D0sumTSmeanALS2
TH1F * h_sumamplitude_depth3_HE1
TH1F * h_corrforxaADDI0_HE
TH2F * h_mapDepth5ADCAmpl12_HE
TH1F * h_eventamplitude_HO
double lsdep_estimator4_HFdepth1_
TH2F * h_mapDepth5ADCAmpl_HE
TH1F * h_sum0ADCAmplperLS8_P1
TH1F * h_totalAmplitudeHOperEvent
TH2F * h_mapDepth1Ampl047_HE
TH2F * h_mapDepth2Ped2_HF
TH1F * h_sumCutADCAmplperLS2
TH2F * h_mapDepth3Pedw3_HE
TH2F * h_mapDepth4Ampl_HB
TH1F * h_sumAmplitudeperLS6
TH1F * h_AamplitudewithPedSubtr_RADDAM_HEzoom0
TH2F * h_maprphinorm0_HE3
TH2F * h_mapDepth6Amplitude_HE
TH2F * h_mapenophinorm_HE5
TH2F * h_map0AiForLS_bad_HF
TH2F * h_recSignalEnergy2_HF3
TH1F * h_sumCutTSmeanAperLS5
TH1F * h_sum0ADCAmplperLS2
TH1F * h_shape_good_channels_HO
double averNOSIGNALsumamplitude_HB
TH2F * h_mapDepth1Pedw1_HB
bool usePedestalSubtraction_
TH2F * h_mapDepth2Pedw2_HF
TH1F * h1_TSnVsAyear20230_HE
TH2F * h_mapTSmeanCalib047_HF
TH1F * h_totalAmplitudeHO
TH1F * h_numberofhitsHFtest
TH2F * h_mapTSmaxCalib047_HE
TH1F * h_runbadrate0_depth1_HF
TH2F * h_mapDepth2Pedw0_HF
TH2F * h_maprphinorm2_HE1
TH2F * h_amplitudechannel0_HF1
double calibrTSmeanHOMin_
TH1F * h_NphiForEta16Depth3RADDAM_HED2P
TH1F * h_sigLayer2RADDAM_HE
TH2F * h_mapDepth4ADCAmpl_HO
TH1F * h_sum0ADCAmplperLS6_P1
TH2F * h_recNoiseEnergy0_HB1
TH2F * h_amplitudechannel0_HF3
TH2F * h_mapDepth2ADCAmpl225Copy_HB
TH2F * h_mapDepth7pedestalw_HE
TH1F * h_sumCutADCAmplperLSdepth3HBu
TH2F * h_amplitudechannel1_HF1
TH1F * h_ADCAmpl_HBCapIdNoError
void analyze(const edm::Event &, const edm::EventSetup &) override
double caliba[ndepth][neta][nphi]
double const adc2fC_QIE10[256]
TH2F * h_2D0sumAmplitudeLS8
TH1F * h_sumCutADCAmplperLS3
double averNOSIGNALoccupancy_HE
TH2F * h_recSignalEnergy1_HB3
TH2F * h_recSignalEnergy0_HB3
TH2F * h_mapDepth6ADCAmpl225_HE
bool studyCalibCellsHist_
TH2F * h_recSignalEnergy1_HB1
TH2F * h_recNoiseEnergy1_HF1
Abs< T >::type abs(const T &t)
TH1F * h_Amplitude_notCapIdErrors_HB2
TH1F * h_energyhitSignal_HE
TH1F * h_runbadrateC_depth1_HE
TH1F * h_orbitNumvsamplitude0_HF
double sumEstimator1[nsub][ndepth][neta][nphi]
TH2F * h_2DsumAmplitudeLS8
TH2F * h_mapDepth3Amplitude225_HF
TH2F * h_mapDepth3TSmeanA225_HE
static const float adc2fC[128]
TH1F * h_sumADCAmplperLS1_P2
TH2F * h_mapDepth2Ampl_HB
TH1F * h_runbadrate0_depth2_HF
TH2F * h_mapDepth2Amplitude_HE
TH1F * h_sum0ADCAmplperLS6
double calibw[ndepth][neta][nphi]
TH1F * h_nevents_per_eachLS
TH2F * h_recSignalEnergy2_HF4
TH2F * h_recSignalEnergy0_HF2
edm::EDGetTokenT< QIE11DigiCollection > tok_qie11_
TH2F * h_2DsumADCAmplEtaPhiLs20
TH1F * h_ADCAmpl345_HBCapIdNoError
TH2F * h_amplitudechannel1_HE6
TH2F * h_amplitudechannel2_HE3
TH2F * h_mapDepth7Amplitude225_HE
TH2F * h_mapDepth3ADCAmpl12_HB
int lsdep_cut1_peak_HBdepth2_
TH2F * h_mapDepth4Ped0_HO
TH1F * h_sum0TSmeanAperLS6
TH1F * h_sum0TSmaxAperLS1
#define DEFINE_FWK_MODULE(type)
double lsdep_estimator1_HBdepth2_
TH2F * h_2D0sumADCAmplLSdepth3HBu
TH1F * h_runbadrate_depth1_HE
TH1F * h_sumADCAmplperLS6_M1
TH2F * h_recNoiseEnergy2_HF3
const HcalLUTCorr * getHcalLUTCorr(const HcalGenericDetId &fId) const
TH2F * h_mapTSmaxCalib_HB
TH2F * h_amplitudechannel2_HE1
TH2F * h_2D0sumTSmeanALS4
TH2F * h_mapDepth2TSmeanA225_HB
TH2F * h_mapDepth1ADCAmpl12_HB
TH2F * h_mapDepth1Pedw2_HF
double averNOSIGNALsumamplitude_HF
TH2F * h_mapDepth5Amplitude225_HE
TH2F * h_mapDepth2ADCAmpl12_HB
TH2F * h_map0AiForLS_good_HF
TH2F * h_2D0sumADCAmplLS5
TH2F * h_mapenophinorm2_HE7
TH2F * h_mapDepth2TSmeanA225_HF
TH2F * h_mapDepth5TSmaxA225_HE
TH1F * h_averNOSIGNALoccupancy_HO
TH1F * h_runnbadchannelsC_depth2_HB
const_iterator begin() const
TH1F * h_sumphiEta16Depth3RADDAM_HED2ALL
TH2F * h_recNoiseEnergy2_HE7
TH1F * h_sum0ADCAmplperLS8
TH2F * h_mapenophinorm_HE1
TH1F * h_sumAmplitudeperLS1_LSselected
TH2F * h_mapDepth1Pedw0_HF
TH2F * h_mapDepth2TSmaxA_HF
TH1F * h_sumCutAmplitudeperLS3
double lsdep_estimator5_HBdepth2_
TH1F * h_sumamplitude_depth3_HE0
TH1F * h_shape0_Ahigh_HB2
bool usecontinuousnumbering_
TH1F * h_AmplitudeHBtest1
TH2F * h_mapDepth3pedestalw_HE
TH1F * h_averNOSIGNALoccupancy_HB
float getSigma(int fCapId1, int fCapId2) const
get correlation element for capId1/2 = 0..3
TH1F * h_ADCAmpl_HBCapIdError
TH2F * h_2DsumADCAmplLS3_LSselected
TH1F * h_runnbadchannels_depth2_HE
TH1F * h_sumADCAmplperLSdepth5HEu
TH2F * h_mapDepth2Ped3_HB
TH1F * h_averSIGNALoccupancy_HO
TH1F * h_totalAmplitudeHE
TH2F * h_maprphinorm2_HE2
TH2F * h_mapDepth7TSmaxA225_HE
TH2F * h_mapDepth3Ped2_HE
TH2F * h_mapDepth1Ped3_HF
TH1F * h_sumADCAmplLS1copy4
TH2F * h_mapDepth7pedestal_HE
TH1F * h_sumCutADCAmplperLS8
TH1F * h_lsnumber_per_eachLS
TH2F * h_mapADCCalib047_HF
double amplitudechannel0[nsub][ndepth][neta][nphi]
TH2F * h_recNoiseEnergy1_HF2
TH1F * h_Amplitude_forCapIdErrors_HO4
int flagestimatornormalization_
TH1F * h_orbitNumvsamplitude0_HO
TH2F * h_mapRatioCalib_HE
TH2F * h_mapADCAmpl12firstpeak_HE
TH2F * h_mapADCAmplfirstpeak0_HE
TH1F * h_sum0TSmaxAperLS6
TH1F * h_pedestalwaver4_HE
double recSignalEnergy0[nsub][ndepth][neta][nphi]
TH2F * h_recNoiseEnergy2_HE1
def gen(fragment, howMuch)
Production test section ####.
TH2F * h_recNoiseEnergy1_HF3
TH2F * h_mapDepth2TSmaxA225_HB
TH1F * h_sumAmplitudeperLS3
TH2F * h_mapDepth4TSmaxA225_HB
TH2F * h_mapDepth6TSmeanA_HE
TH1F * h_sumCutADCAmplperLS6u
TH2F * h_map0AiForLS_bad_HO
double recSignalEnergy1[nsub][ndepth][neta][nphi]
TH1F * h_sumADCAmplperLS1_M2
TH1F * h_sumTSmeanAperLS1_LSselected
double lsdep_estimator1_HBdepth4_
TH1F * h_sumADCAmplperLS3_P1
TH1F * h_sumADCAmplLS1copy3
TH2F * h_mapDepth2Pedw2_HB
TH1F * h_sum0ADCAmplperLSdepth7HEu
TH1F * h_runbadrate_depth3_HE
TH2F * h_mapAiForLS_good_HB
TH1F * h_runbadrate0_depth2_HE
TH2F * h_mapDepth5Ampl_HE
double averNOSIGNALoccupancy_HO
TH1F * h_sigLayer1RADDAM_HE
TH2F * h_mapADCAmpl12secondpeak0_HE
TH1F * h_sumCutADCAmplperLSdepth4HBu
TH1F * h_sumPedestalperLS7
TH2F * h_mapDepth5Error_HE
TH1F * h_shape_bad_channels_HF
int flagtodefinebadchannel_
double recNoiseEnergy0[nsub][ndepth][neta][nphi]
TH1F * h_shape0_Ahigh_HB3
TH2F * h_mapDepth4TS012_HO
TH1F * h_sum0ADCAmplperLS7
TH2F * h_mapDepth3ADCAmpl225_HE
TH1F * h_sumCutAmplperLS1
TH2F * h_mapDepth4ADCAmpl12_HO
TH2F * h_mapDepth1AmplE34_HE
TH2F * h_mapAiForLS_good_HO
TH1F * h_bcnbadrate0_depth2_HF
TH2F * h_mapDepth4AmplE34_HF
TH1F * h_sum0PedestalperLS5
TH2F * h_2D0sumAmplitudeLS5
TH2F * h_mapDepth2TSmaxA225_HF
TH2F * h_mapDepth3Ampl047_HB
TH1F * h_sumCutTSmeanAperLS2
TH2F * h_maprphinorm0_HE1
TH1F * h_sumamplitude_depth2_HB1
TH2F * h_mapADCAmplsecondpeak_HE
TH2F * h_mapDepth3pedestalw_HF
TH2F * h_mapDepth2TSmeanA_HE
TH1F * h_sumamplitude_depth1_HE1
double lsdep_estimator4_HEdepth2_
TH2F * h_amplitudechannel2_HB1
TH2F * h_mapDepth1ADCAmpl_HB
TH1F * h_bcnbadrate0_depth2_HB
TH2F * h_2DsumADCAmplEtaPhiLs10
TH2F * h_recSignalEnergy2_HB3
TH1F * h_runnbadchannelsC_depth3_HE
const_iterator end() const
TH2F * h_mapDepth6pedestalw_HE
TH2F * h_mapDepth3TSmeanA_HB
const HcalFrontEndId getHcalFrontEndId(const DetId &)
double sumEstimator0[nsub][ndepth][neta][nphi]
TH2F * h_2DsumAmplitudeLS2
TH2F * h_amplitudechannel1_HE2
TH2F * h_mapDepth3Ped1_HE
TH2F * h_map0AiForLS_good_HO
int lsdep_cut1_peak_HOdepth4_
TH2F * h_2D0sumPedestalLS5
double lsdep_estimator5_HBdepth1_
TH2F * h_mapADCAmpl12firstpeak0_HE
TH2F * h_mapDepth1Error_HE
double sumEstimator2[nsub][ndepth][neta][nphi]
TH2F * h_mapTSmaxCalib047_HF
TH1F * h_sum0TSmeanAperLS4
TH2F * h_mapDepth3pedestal_HF
TH2F * h_mapDepth3ADCAmpl_HB
TH2F * h_recSignalEnergy2_HF1
TH1F * h_sumADCAmplEtaPhiLs_bbbc
TH1F * h_pedestalwaver4_HF
TH2F * h_recNoiseEnergy2_HE3
TH2F * h_mapADCCalib047_HB
TH2F * h_mapDepth1Amplitude_HE
TH1F * h_sum0TSmaxAperLS8
TH2F * h_mapDepth1ADCAmpl225Copy_HB
TH2F * h_mapDepth2AmplE34_HB
double lsdep_estimator3_HEdepth1_
TH2F * h_mapDepth1Pedw0_HB
TH2F * h_2D0sumPedestalLS3
TH1F * h_sumAmplitudeperLS4
TH2F * h_2DsumADCAmplEtaPhiLs00
TH2F * h_2D0sumADCAmplLS8
TH2F * h_mapDepth1Ped2_HF
TH1F * h_energyhitNoise_HF
TH1F * h_sumADCAmplperLS5
constexpr int samples() const
total number of samples in the digi
TH2F * h_mapDepth3Error_HF
TH2F * h_amplitudechannel2_HB2
TH2F * h_mapDepth3TSmaxA_HF
TH2F * h_2D0sumADCAmplLS2
TH2F * h_maprphinorm0_HE7
TH2F * h_maprphinorm0_HE5
TH1F * h_runbadrate_depth2_HB
TH2F * h_mapDepth2Ped0_HB
TH2F * h_2D0sumPedestalLS2
TH1F * h_Amplitude_notCapIdErrors_HE1
TH2F * h_mapADCAmpl11secondpeak0_HE
TH2F * h_recNoiseEnergy0_HE3
TH1F * h_sum0AmplitudeperLS4
TH2F * h_mapDepth4TSmaxA_HF
TH1F * h_pedestalaver9_HE
double lsdep_estimator5_HFdepth1_
TH2F * h_maprphinorm2_HE6
TH2F * h_2DsumADCAmplEtaPhiLs2
TH2F * h_recNoiseEnergy2_HE4
TH2F * h_mapDepth1Pedw1_HF
TH1F * h_eventamplitude_HB
TH1F * h_pedestalwaver9_HB
TH1F * h_sumamplitude_depth2_HF0
TH2F * h_recSignalEnergy0_HE4
~CMTRawAnalyzer() override
TH1F * h_shape0_Ahigh_HB0
TH2F * h_mapenophinorm2_HE5
TH2F * h_mapDepth4Amplitude_HE
TH1F * h_sum0ADCAmplperLS3
TH2F * h_2D0sumAmplitudeLS6
void fillDigiAmplitudeHF(HFDigiCollection::const_iterator &digiItr)
TH2F * h_mapDepth3Amplitude225_HE
TH1F * h_sumamplitude_depth1_HF0
TH1F * h_sum0TSmeanAperLS3
TH1F * h_sumCutADCAmplperLS4
TH2F * h_map0AiForLS_bad_HE
TH1F * h_sumCutADCAmplperLS7u
double lsdep_estimator5_HOdepth4_
TH2F * h_recNoiseEnergy2_HF4
TH1F * h_sumAmplitudeperLS5
TH2F * h_recSignalEnergy0_HE1
TH2F * h_mapDepth2ADCAmpl225Copy_HE
TH2F * h_mapenophinorm_HE4
TH2F * h_mapDepth1Pedw3_HF
TH2F * h_mapDepth2TSmeanA_HB
TH2F * h_mapTSmeanCalib_HB
int flagIterativeMethodCalibrationGroupReco_
TH2F * h_mapDepth7AmplE34_HE
TH2F * h_amplitudechannel1_HB3
void beginRun(const edm::Run &r, const edm::EventSetup &iSetup) override
TH2F * h_mapWidthCalib_HF
TH2F * h_mapenophinorm_HE6
TH1F * h_nbadchannels_depth1_HE
TH2F * h_mapDepth3TSmeanA_HF
TH1F * h_sumPedestalperLS2
TH2F * h_mapTSmeanCalib047_HE
TH1F * h_sumTSmaxAperLS1_LSselected
TH1F * h_runbadrateC_depth4_HO
TH1F * h_sumamplitude_depth4_HO1
TH1F * h_runnbadchannels_depth1_HB
TH2F * h_mapDepth3Ampl047_HF
TH1F * h_runnbadchannelsC_depth1_HB
TH1F * h_sum0TSmeanAperLS2
TH2F * h_mapenophinorm_HE7
TH1F * h_bcnbadrate0_depth1_HE
TH1F * h_averNOSIGNALsumamplitude_HF
TH2F * h_mapDepth4Ped2_HO
double sumEstimator6[nsub][ndepth][neta][nphi]
TH2F * h_2DsumADCAmplLS6_LSselected
double lsdep_estimator1_HEdepth6_
TH1F * h_sumCutTSmaxAperLS6
TH2F * h_mapDepth7Error_HE
const HcalGain * getGain(const HcalGenericDetId &fId) const
TH2F * h_2D0sumADCAmplLSdepth4HBu
TH2F * h_maprphinorm2_HE4
TH2F * h_recNoiseEnergy1_HE2
TH1F * h_sumCutTSmaxAperLS8
TH2F * h_mapDepth4TSmeanA_HO
TH2F * h_mapDepth3Ped0_HE
constexpr int dccid() const
get the (Hcal local) DCC id for VME, crate number for uTCA
TH1F * h_averSIGNALsumamplitude_HE
TH2F * h_mapRatioCalib047_HF
double signal3[ndepth][neta][nphi]
TH2F * h_mapDepth2Ped1_HF
TH1F * h_bcnbadrate0_depth1_HB
TH2F * h_2D0sumAmplitudeLS7
TH1F * h_bcnnbadchannels_depth2_HB
TH2F * h_mapDepth1Pedw2_HB
TH1F * h_sumCutADCAmplperLS1
TH1F * h_sumCutAmplperLS2
TH2F * h_mapDepth6TSmaxA225_HE
TH2F * h_mapDepth4pedestalw_HO
TH1F * h_Eta16Depth3RADDAM_HED2
TH1F * h_sumamplitude_depth1_HB1
TH2F * h_mapDepth2TSmeanA_HF
TH1F * h_pedestalaver9_HF
double averSIGNALoccupancy_HE
TH2F * h_amplitudechannel2_HF2
double lsdep_estimator1_HEdepth2_
TH2F * h_mapDepth6ADCAmpl_HE
TH2F * h_mapDepth5Amplitude_HE
TH2F * h_mapDepth3TSmaxA225_HF
TH2F * h_mapDepth4Error_HF
TH2F * h_mapDepth1ADCAmpl_HF
TH1F * h_sum0PedestalperLS7
TH1F * h_runnbadchannelsC_depth4_HO
double lsdep_estimator1_HFdepth1_
TH2F * h_mapDepth4pedestalw_HB
TH1F * h_sigLayer1RADDAM5_HE
TH2F * h_mapDepth4pedestalw_HE
int badchannels[nsub][ndepth][neta][nphi]
TH1F * h_sum0ADCAmplperLS1_P1
TH1F * h_eventoccupancy_HO
TH1F * h_orbitNumvsamplitude_HB
TH2F * h_amplitudechannel1_HE7
TH1F * h_shape0_good_channels_HF
TH2F * h_mapRatioCalib047_HO
TH2F * h_mapDepth3Ampl_HF
TH2F * h_mapDepth1ADCAmpl12SiPM_HE
TH2F * h2_TSnVsAyear2023_HF
TH1F * h_sumphiEta16Depth3RADDAM_HED2P
TH2F * h_2D0sumADCAmplLSdepth5HEu
TH2F * h_mapDepth2pedestal_HF
TH2F * h_mapDepth7TSmaxA_HE
TH2F * h_mapDepth6ADCAmpl225Copy_HE
TH2F * h_mapDepth3TSmaxA_HB
TH1F * h_RatioOccupancy_HEM
TH1F * h_sumCutAmplitudeperLS6
TH1F * h_bcnnbadchannels_depth1_HF
TH2F * h_recNoiseEnergy2_HB3
TH2F * h_mapDepth1Amplitude225_HB
TH2F * h_2DsumADCAmplLS5_LSselected
TH1F * h_averNOSIGNALsumamplitude_HB
TH2F * h_mapDepth2Pedw2_HE
TH1F * h_sumADCAmplperLS3
TH1F * h_sumCutTSmaxAperLS3
TH1F * h_sumPedestalperLS6
TH1F * h_sum0TSmaxAperLS5
TH1F * h_sumCutAmplperLS8
TH2F * h_recNoiseEnergy0_HF2
TH2F * h_mapDepth1pedestalw_HE
TH1F * h_bcnnbadchannels_depth1_HE
TH2F * h_2DsumADCAmplLSdepth6HEu
double averNOSIGNALoccupancy_HF
TH2F * h_mapTSmeanCalib_HE
TH2F * h_mapDepth1ADCAmpl225Copy_HE
TH1F * h_totalAmplitudeHF
TH2F * h_mapDepth7Ampl_HE
TH1F * h_sum0TSmeanAperLS1
TH2F * h2_TSnVsAyear2023_HB
TH1F * h_averNOSIGNALoccupancy_HE
edm::EDGetTokenT< HFDigiCollection > tok_hf_
TH2F * h_mapDepth2ADCAmpl225_HF
TH1F * h_runbadrateC_depth3_HE
TH1F * h_sum0ErrorBperLS8
TH2F * h_recSignalEnergy2_HB4
double lsdep_estimator5_HEdepth2_
double signal[ndepth][neta][nphi]
TH1F * h_sumADCAmplEtaPhiLs_lscounterM1orbitNum
constexpr int htrTopBottom() const
get the htr top/bottom (1=top/0=bottom), valid for VME
TH1F * h_RatioOccupancy_HFP
TH1F * h_totalAmplitudeHB
TH2F * h_mapenophinorm_HE3
double lsdep_estimator3_HFdepth2_
TH1F * h_shape_bad_channels_HE
TH2F * h_mapDepth4Ampl047_HE
TH2F * h_mapDepth1Pedw0_HE
TH1F * h_numberofhitsHBtest
TH1F * h_eventamplitude_HE
TH2F * h_mapGetRMSOverNormalizedSignal_HB
TH2F * h_mapDepth1Ampl047_HF
double lsdep_estimator1_HEdepth1_
TH2F * h_2D0sumTSmeanALS3
TH2F * h_mapDepth6TSmeanA225_HE
TH1F * h_sumTSmeanAperLS5
TH1F * h_corrforxaMAIN_HE
TH1F * h_averSIGNALoccupancy_HF
TH1F * h_A_Depth3RADDAM_HE
TH1F * h_bcnvsamplitude_HF
TH2F * h_mapDepth3Pedw1_HE
TH2F * h_recSignalEnergy1_HF4
const HcalQIECoder * getHcalCoder(const HcalGenericDetId &fId) const
TH2F * h_recNoiseEnergy0_HF4
TH1F * h_sigLayer2RADDAM6_HED2
TH1F * h_sumCutTSmeanAperLS6
TH2F * h_mapDepth4Pedw2_HO
TH2F * h_mapDepth6TSmaxA_HE
double calibrTSmeanHEMin_
TH2F * h_mapDepth2ADCAmpl12_HF
double calibrTSmeanHBMin_
TH2F * h_2D0sumAmplitudeLS3
TH2F * h_maprphinorm2_HE3
TH2F * h_2DsumADCAmplLS7_LSselected
int mapRADDAM_HED20[ndepth][neta]
TH2F * h_mapDepth2pedestalw_HE
TH1F * h_RatioOccupancy_HBM
TH1F * h_sum0ADCAmplperLSdepth4HEu
long int gsmdepth1sipm[npfit][neta][nphi][ndepth]
int flagabortgaprejected_
TH2F * h_mapDepth2RADDAM0_HE
TH2F * h_mapDepth7Amplitude_HE
TH2F * h_mapDepth1pedestalw_HF
constexpr int spigot() const
get the spigot (input number on DCC, AMC card number for uTCA)
TH1F * h_sumamplitude_depth1_HE0
TH2F * h_mapDepth1pedestal_HF
TH2F * h_mapDepth5pedestal_HE
TH2F * h_mapRatioCalib_HF
void fillDigiErrorsHO(HODigiCollection::const_iterator &digiItr)
TH1F * h_A_Depth1RADDAM_HE
TH1F * h_averNOSIGNALoccupancy_HF
TH1F * h_sum0ErrorBperLS7
TH1F * h_sumamplitude_depth1_HF
TH1F * h_runbadrate0_depth1_HE
TH2F * h_mapDepth3Ampl_HE
TH2F * h_mapenophinorm2_HE4
TH1F * h_sum0ADCAmplperLS8_M1
TH2F * h_2DsumADCAmplLS2_LSselected
int lsdep_cut1_peak_HEdepth1_
TH2F * h_mapDepth3TSmeanA225_HB
TH2F * h_mapDepth4pedestal_HF
TH1F * h_sum0AmplitudeperLS6
double lsdep_estimator2_HBdepth2_
TH2F * h_mapDepth2pedestalw_HB
TH2F * h_2DsumADCAmplEtaPhiLs30
TH2F * h_mapDepth6Ampl047_HE
TH1F * h_sum0ADCAmplperLS6u
TH2F * h_mapDepth1Ped1_HF
TH2F * h_mapDepth4ADCAmpl225_HE
TH2F * h_recNoiseEnergy2_HB2
edm::EDGetTokenT< HODigiCollection > tok_ho_
TH2F * h_mapAiForLS_bad_HO
TH1F * h_bcnbadrate0_depth3_HE
TH1F * h_sumADCAmplEtaPhiLs_orbitNum
TH1F * h_AmplitudeHBtest6
TH2F * h_mapDepth4TSmaxA225_HE
TH2F * h_recNoiseEnergy2_HB4
TH2F * h_mapADCAmplsecondpeak0_HE
TH1F * h1_TSnVsAyear2023_HO
TH2F * h_mapDepth2ADCAmpl225Copy_HF
TH1F * h_gsmdifferencefit2_HE
T * make(const Args &...args) const
make new ROOT object
TH1F * h_averSIGNALsumamplitude_HB
TH1F * h_sumCutAmplitudeperLS2
TH2F * h_mapDepth2Ampl047_HF
TH1F * h_sum0ADCAmplperLS4
TH2F * h_mapDepth3ADCAmplSiPM_HE
TH1F * h_ADCAmpl345_HBCapIdError
void fillDigiErrorsHFQIE10(QIE10DataFrame qie10df)
TH1F * h_bcnnbadchannels_depth2_HE
void endRun(const edm::Run &r, const edm::EventSetup &iSetup) override
double calibrTSmeanHBMax_
TH2F * h_mapDepth2ADCAmpl_HE
TH2F * h_mapDepth3AmplE34_HF
TH2F * h_mapDepth4ADCAmpl225Copy_HE
TH1F * h_sumADCAmplEtaPhiLs_bbb1
TH2F * h_mapDepth1pedestal_HE
TH1F * h_sumAmplitudeperLS2
TH2F * h_amplitudechannel1_HE3
TH1F * h_AamplitudewithPedSubtr_RADDAM_HE
TH1F * h_bcnnbadchannels_depth3_HE
TH2F * h_mapDepth1TSmeanA_HF
double lsdep_estimator3_HEdepth3_
TH1F * h_gsmdifferencefit6_HE
double lsdep_estimator4_HEdepth1_
TH2F * h_mapDepth4Ped3_HO
TH2F * h_mapDepth2Ped3_HE
TH1F * h_sumCutAmplitudeperLS5
TH1F * h_sum0ADCAmplperLS1_M2
TH1F * h_AamplitudewithPedSubtr_RADDAM_HEzoom1
double averNOSIGNALsumamplitude_HE
TH1F * h_sigLayer2RADDAM0_HE
TH1F * h_sumamplitudechannel_HO
TH1F * h_pedestalaver4_HB
TH2F * h_mapWidthCalib_HB
TH1F * h_sumAmplitudeperLS7
TH2F * h_recSignalEnergy1_HF1
TH1F * h_sum0ADCAmplperLS1_M1
TH2F * h_mapDepth1TSmaxA225_HB
TH1F * h_orbitNumvsamplitude0_HB
TH2F * h_mapDepth2ADCAmpl12_HE
int lsdep_cut3_max_HBdepth1_
double recNoiseEnergy1[nsub][ndepth][neta][nphi]
TH2F * h_recNoiseEnergy1_HE7
double sumEstimator4[nsub][ndepth][neta][nphi]
TH2F * h_mapDepth2Pedw0_HB
TH1F * h1_TSnVsAyear2023_HE
double averSIGNALsumamplitude_HB
TH2F * h_maprphinorm0_HE6
TH2F * h_amplitudechannel2_HE4
TH1F * h_sumCutTSmaxAperLS7
double averSIGNALsumamplitude_HF
TH2F * h_mapTSmeanCalib_HO
TH2F * h_mapDepth2AmplE34_HF
void fillDigiAmplitudeHFQIE10(QIE10DataFrame qie10df)
TH1F * h_eventamplitude_HF
TH1F * h_sum0ADCAmplperLSdepth6HEu
TH2F * h_mapDepth3pedestal_HE
TH1F * h_numberofhitsHOtest
TH1F * h_sumamplitude_depth2_HF1
TH2F * h_2D0sumTSmeanALS1
TH1F * h_totalAmplitudeHFperEvent
TH1F * h_orbitNumvsamplitude_HO
TH2F * h_mapDepth3ADCAmpl_HE
TH1F * h_totalAmplitudeHBperEvent
TH2F * h_2DsumAmplitudeLS7
TH2F * h_mapGetRMSOverNormalizedSignal0_HB
TH2F * h_mapDepth1ADCAmpl225Copy_HF
TH2F * h_recSignalEnergy2_HB1
TH2F * h_mapenophinorm2_HE3
TH2F * h_recSignalEnergy2_HE1
double lsdep_estimator5_HEdepth3_
TH2F * h_recSignalEnergy0_HB4
double const adc2fC_QIE11_shunt1[256]
TH1F * h_runbadrate0_depth1_HB
TH1F * h_sum0AmplitudeperLS7
TH2F * h_mapDepth4ADCAmpl_HF
TH2F * h_mapDepth3linADCAmpl12_HE
TH2F * h_recSignalEnergy0_HE5
TH1F * h_shape_good_channels_HB
TH2F * h_map0AiForLS_good_HE
TH2F * h_mapDepth1Ped1_HB
TH2F * h_recSignalEnergy0_HE7
TH2F * h_2DsumPedestalLS6
TH1F * h_nevents_per_eachRealLS
double lsdep_estimator1_HEdepth7_
TH1F * h_shape0_Ahigh_HB1
TH1F * h_amplitudeaveragedbydepthes_HE
TH1F * h_sumADCAmplperLS3_P2
TH1F * h1_TSnVsAyear20230_HO
TH2F * h_mapDepth4Amplitude_HB
TH2F * h_mapRatioCalib047_HB
TH2F * h_mapDepth1ADCAmpl_HE
TH1F * h_shape0_bad_channels_HE
TH1F * h_shape0_bad_channels_HO
TH2F * h_amplitudechannel1_HF2
TH2F * h_mapDepth2Pedw0_HE
TH2F * h_2D0sumADCAmplLSdepth7HEu
TH2F * h_mapDepth6Error_HE
TH1F * h_sumTSmeanAperLS6
TH1F * h_sumADCAmplperLSdepth4HEu
TH2F * h_2D0sumADCAmplLSdepth4HEu
TH1F * h_sumPedestalperLS5
TH2F * h_mapDepth1Amplitude225_HF
TH2F * h_mapWidthCalib_HE
TH1F * h_sumCutTSmaxAperLS5
TH2F * h_mapDepth3ADCAmpl12SiPM_HE
TH2F * h_mapTSmaxCalib047_HO
const HcalTimeCorr * getHcalTimeCorr(const HcalGenericDetId &fId) const
TH1F * h_sigLayer2RADDAM5_HE
double calibrTSmeanHFMax_
TH2F * h_amplitudechannel0_HB2
TH2F * h_mapDepth4TSmeanA225_HF
TH2F * h_mapDepth4TSmaxA_HE
TH2F * h_mapGetRMSOverNormalizedSignal0_HF
double averSIGNALsumamplitude_HE
TH2F * h_mapDepth4ADCAmpl225_HO
TH2F * h_2DsumADCAmplLSdepth7HEu
TH2F * h_mapDepth3Pedw0_HE
Readout chain identification for Hcal.
TH1F * h_sumADCAmplperLS6_M2
TH1F * h_bcnvsamplitude0_HE
TH1F * h_sum0TSmeanAperLS5
TH1F * h_bcnbadrate0_depth2_HE
TH2F * h_mapAiForLS_bad_HB
TH2F * h_2D0sumADCAmplLS7
TH2F * h_mapDepth2Pedw1_HF
TH2F * h_mapDepth1Ped0_HF
TH2F * h_mapDepth1RADDAM_HE
TH1F * h_NphiForEta16Depth3RADDAM_HED2
TH2F * h_mapDepth4TSmeanA225_HO
TH2F * h_recNoiseEnergy2_HF2
TH1F * h_sumCutTSmaxAperLS4
TH2F * h_mapDepth4TSmaxA225_HO
TH1F * h_Amplitude_notCapIdErrors_HE2
TH2F * h_mapDepth5ADCAmpl225Copy_HE
TH2F * h_mapDepth7ADCAmpl12_HE
constexpr int iphi() const
get the cell iphi
TH1F * h_gsmdifferencefit1_HE
TH2F * h_recNoiseEnergy0_HE1
TH2F * h_mapenophinorm2_HE2
TH2F * h_mapDepth1ADCAmplSiPM_HE
TH1F * h_sum0PedestalperLS2
TH2F * h_mapDepth3AmplE34_HB
TH2F * h_2D0sumAmplitudeLS4
TH1F * h_sumADCAmplEtaPhiLs_lscounterM1
TH2F * h_mapDepth6AmplE34_HE
TH1F * h_sum0TSmeanAperLS8
TH2F * h_mapDepth1TSmeanA_HE
TH2F * h_mapDepth3ADCAmpl225Copy_HB
TH2F * h_recNoiseEnergy1_HE4
TH2F * h_mapenophinorm2_HE6
TH1F * h_runbadrateC_depth1_HB
TH2F * h_amplitudechannel1_HE1
TH1F * h_shape_bad_channels_HB
constexpr int samples() const
total number of samples in the digi
TH1F * h_Amplitude_notCapIdErrors_HF1
double averSIGNALoccupancy_HF
TH2F * h_mapDepth3AmplE34_HE
TH2F * h_recNoiseEnergy2_HE2
TH2F * h_mapDepth3TSmeanA_HE
const edm::ESGetToken< HcalDbService, HcalDbRecord > tokDB_
TH2F * h_mapDepth4pedestal_HO
TH2F * h_mapDepth4TSmeanA_HE
TH1F * h_sumCutAmplperLS3
TH2F * h_mapDepth1Amplitude_HB
double dR(double eta1, double phi1, double eta2, double phi2)
uint16_t *__restrict__ uint16_t const *__restrict__ adc
TH1F * h_sumPedestalperLS8
TH2F * h_mapDepth4Amplitude225_HE
edm::EDGetTokenT< HcalCalibDigiCollection > tok_calib_
TH2F * h_recSignalEnergy2_HE5
TH1F * h_orbitNumvsamplitude_HE
TH1F * h_sumADCAmplperLSdepth7HEu
double recNoiseEnergy2[nsub][ndepth][neta][nphi]
TH2F * h_mapDepth3RADDAM0_HE
TH2F * h_mapDepth2pedestal_HE
void fillDigiAmplitude(HBHEDigiCollection::const_iterator &digiItr)
TH2F * h_amplitudechannel0_HE5
TH1F * h_averSIGNALoccupancy_HE
TH1F * h_runnbadchannels_depth2_HF
TH2F * h_2DsumADCAmplLSdepth4HFu
TH1F * h_bcnvsamplitude0_HO
TH1F * h_sumCutADCAmplperLSdepth5HEu
TH1F * h_nbadchannels_depth1_HF
TH2F * h_recSignalEnergy0_HE3
TH2F * h_mapDepth1Ped0_HE
TH2F * h_map0AiForLS_good_HB
TH2F * h_mapRatioCalib047_HE
TH2F * h_amplitudechannel2_HE2
TH2F * h_mapDepth3TSmaxA_HE
TH1F * h_sumCutADCAmplperLS7
TH2F * h_2D0sumADCAmplLS4
TH1F * h_sumCutAmplitudeperLS4
constexpr int depth() const
get the tower depth
bool studyCapIDErrorsHist_
bool studyPedestalCorrelations_