Skip to content

Commit 3ccf0c4

Browse files
author
Sabrina Hernandez
committed
added histograms
1 parent f3d6d5c commit 3ccf0c4

File tree

1 file changed

+34
-10
lines changed

1 file changed

+34
-10
lines changed

‎PWGCF/Flow/Tasks/flowZdcTask.cxx‎

Lines changed: 34 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -202,9 +202,11 @@ struct FlowZdcTask {
202202
histos.add("ZNVsFT0A", ";T0A (#times 1/100);ZNA+ZNC Amplitude;", kTH2F, {{{nBinsAmpFT0, 0., maxAmpFT0}, {nBinsZDC, -0.5, maxZn}}});
203203
histos.add("ZNVsFT0C", ";T0C (#times 1/100);ZNA+ZNC Amplitude;", kTH2F, {{{nBinsAmpFT0, 0., maxAmpFT0}, {nBinsZDC, -0.5, maxZn}}});
204204
histos.add("ZNVsFT0M", ";T0A+T0C (#times 1/100);ZNA+ZNC Amplitude;", kTH2F, {{{nBinsAmpFT0, 0., maxAmpFT0M}, {nBinsZDC, -0.5, maxZn}}});
205+
histos.add("CommonZNVsFT0M", ";T0A+T0C (#times 1/100);ZPA+ZPC Amplitude;", kTH2F, {{{nBinsAmpFT0, 0., maxAmpFT0M}, {nBinsZDC, -0.5, maxZn}}});
205206
histos.add("ZPVsFT0A", ";T0A (#times 1/100);ZPA+ZPC Amplitude;", kTH2F, {{{nBinsAmpFT0, 0., maxAmpFT0}, {nBinsZDC, -0.5, maxZp}}});
206207
histos.add("ZPVsFT0C", ";T0C (#times 1/100);ZPA+ZPC Amplitude;", kTH2F, {{{nBinsAmpFT0, 0., maxAmpFT0}, {nBinsZDC, -0.5, maxZp}}});
207208
histos.add("ZPVsFT0M", ";T0A+T0C (#times 1/100);ZPA+ZPC Amplitude;", kTH2F, {{{nBinsAmpFT0, 0., maxAmpFT0M}, {nBinsZDC, -0.5, maxZp}}});
209+
histos.add("CommonZPVsFT0M", ";T0A+T0C (#times 1/100);ZPA+ZPC Amplitude;", kTH2F, {{{nBinsAmpFT0, 0., maxAmpFT0M}, {nBinsZDC, -0.5, maxZp}}});
208210
histos.add("ZNAVsFT0A", ";T0A (#times 1/100);ZNA Amplitude;", kTH2F, {{{nBinsAmpFT0, 0., maxAmpFT0}, {nBinsZDC, -0.5, maxZn}}});
209211
histos.add("ZNAVsFT0C", ";T0C (#times 1/100);ZNA Amplitude;", kTH2F, {{{nBinsAmpFT0, 0., maxAmpFT0}, {nBinsZDC, -0.5, maxZn}}});
210212
histos.add("ZNAVsFT0M", ";T0A+T0C (#times 1/100);ZNA Amplitude;", kTH2F, {{{nBinsAmpFT0, 0., maxAmpFT0M}, {nBinsZDC, -0.5, maxZn}}});
@@ -279,6 +281,10 @@ struct FlowZdcTask {
279281
histos.add("ampZpa", ";ZPA Amplitude;Entries;", kTH1F, {{nBinsZDC, -0.5, maxZp}});
280282
histos.add("ampZnc", ";ZNC Amplitude;Entries;", kTH1F, {{nBinsZDC, -0.5, maxZn}});
281283
histos.add("ampZpc", ";ZPC Amplitude;Entries;", kTH1F, {{nBinsZDC, -0.5, maxZp}});
284+
histos.add("commonZna", ";ZNA Amplitude;Entries;", kTH1F, {{nBinsZDC, -0.5, maxZn}});
285+
histos.add("commonZpa", ";ZPA Amplitude;Entries;", kTH1F, {{nBinsZDC, -0.5, maxZp}});
286+
histos.add("commonZnc", ";ZNC Amplitude;Entries;", kTH1F, {{nBinsZDC, -0.5, maxZn}});
287+
histos.add("commonZpc", ";ZPC Amplitude;Entries;", kTH1F, {{nBinsZDC, -0.5, maxZp}});
282288
histos.add("ampZEM1", ";ZEM1 Amplitude;Entries;", kTH1F, {{nBinsZDC, -0.5, maxZem}});
283289
histos.add("ampZEM2", ";ZEM2 Amplitude;Entries;", kTH1F, {{nBinsZDC, -0.5, maxZem}});
284290
histos.add("ZnVsZem", "ZnVsZEM; ZEM; ZNA + ZNC", kTH2F, {{{nBinsZDC, -0.5, maxZem}, {nBinsZDC, -0.5, maxZn}}});
@@ -428,16 +434,18 @@ struct FlowZdcTask {
428434
float znC = zdc.amplitudeZNC() / cfgCollisionEnergy;
429435
float zpA = zdc.amplitudeZPA() / cfgCollisionEnergy;
430436
float zpC = zdc.amplitudeZPC() / cfgCollisionEnergy;
437+
float commonSumZnc = zdc.energyCommonZNC() / cfgCollisionEnergy;
438+
float commonSumZna = zdc.energyCommonZNA() / cfgCollisionEnergy;
439+
float commonSumZpc = zdc.energyCommonZPC() / cfgCollisionEnergy;
440+
float commonSumZpa = zdc.energyCommonZPA() / cfgCollisionEnergy;
431441
if (applyZdcCorrection) {
432442
const float a = zpaCoeff;
433443
const float b = zpcCoeff;
434444
zpA = zpA - a * znA;
445+
commonSumZpa = commonSumZpa - a * commonSumZna;
435446
zpC = zpC - b * znC;
436-
}
437-
float commonSumZnc = zdc.energyCommonZNC() / cfgCollisionEnergy;
438-
float commonSumZna = zdc.energyCommonZNA() / cfgCollisionEnergy;
439-
float commonSumZpc = zdc.energyCommonZPC() / cfgCollisionEnergy;
440-
float commonSumZpa = zdc.energyCommonZPA() / cfgCollisionEnergy;
447+
commonSumZpc = commonSumZpc - b * commonSumZnc;
448+
}
441449
float aZEM1{zdc.amplitudeZEM1()};
442450
float aZEM2{zdc.amplitudeZEM2()};
443451
float sumZEMs{aZEM1 + aZEM2};
@@ -492,12 +500,14 @@ struct FlowZdcTask {
492500
histos.fill(HIST("ZN"), znA + znC);
493501
histos.fill(HIST("ZNVsFT0C"), aT0C / 100., znA + znC);
494502
histos.fill(HIST("ZNVsFT0M"), (aT0A + aT0C) / 100., znA + znC);
503+
histos.fill(HIST("CommonZNVsFT0M"), (aT0A + aT0C) / 100., commonSumZna + commonSumZnc);
495504
}
496-
if (((tZPA >= minTdcZn) && (tZPA <= maxTdcZn)) && ((tZPC >= minTdcZn) && (tZPC <= maxTdcZn))) {
505+
if (((tZPA >= minTdcZp) && (tZPA <= maxTdcZp)) && ((tZPC >= minTdcZp) && (tZPC <= maxTdcZp))) {
497506
histos.fill(HIST("ZPAVsZPC"), zpC, zpA);
498507
histos.fill(HIST("ZPVsFT0A"), aT0A / 100., zpA + zpC);
499508
histos.fill(HIST("ZPVsFT0C"), aT0C / 100., zpA + zpC);
500509
histos.fill(HIST("ZPVsFT0M"), (aT0A + aT0C) / 100., zpA + zpC);
510+
histos.fill(HIST("CommonZPVsFT0M"), (aT0A + aT0C) / 100., commonSumZpa + commonSumZpc);
501511
}
502512
if (((tZNA >= minTdcZn) && (tZNA <= maxTdcZn)) && ((tZPA >= minTdcZp) && (tZPA <= maxTdcZp)))
503513
histos.fill(HIST("ZNAVsZPA"), zpA, znA);
@@ -525,11 +535,13 @@ struct FlowZdcTask {
525535
histos.fill(HIST("ZPVsFT0A"), aT0A / 100., zpA + zpC);
526536
histos.fill(HIST("ZPVsFT0C"), aT0C / 100., zpA + zpC);
527537
histos.fill(HIST("ZPVsFT0M"), (aT0A + aT0C) / 100., zpA + zpC);
538+
histos.fill(HIST("CommonZPVsFT0M"), (aT0A + aT0C) / 100., commonSumZpa + commonSumZpc);
528539
histos.fill(HIST("ZPAVsFT0A"), aT0A / 100., zpA);
529540
histos.fill(HIST("ZPAVsFT0C"), aT0C / 100., zpA);
530541
histos.fill(HIST("ZPAVsFT0M"), (aT0A + aT0C) / 100., zpA);
531542
histos.fill(HIST("ZNVsFT0C"), aT0C / 100., znA + znC);
532543
histos.fill(HIST("ZNVsFT0M"), (aT0A + aT0C) / 100., znA + znC);
544+
histos.fill(HIST("CommonZNVsFT0M"), (aT0A + aT0C) / 100., commonSumZna + commonSumZnc);
533545
histos.fill(HIST("ZPCVsFT0A"), aT0A / 100., zpC);
534546
histos.fill(HIST("ZPCVsFT0C"), aT0C / 100., zpC);
535547
histos.fill(HIST("ZPCVsFT0M"), (aT0A + aT0C) / 100., zpC);
@@ -690,14 +702,22 @@ struct FlowZdcTask {
690702
auto tZPA = zdc.timeZPA();
691703
auto tZPC = zdc.timeZPC();
692704
if (isTDCcut) {
693-
if ((tZNA >= minTdcZn) && (tZNA <= maxTdcZn))
705+
if ((tZNA >= minTdcZn) && (tZNA <= maxTdcZn)) {
694706
histos.fill(HIST("ampZna"), znA);
695-
if ((tZNC >= minTdcZn) && (tZNC <= maxTdcZn))
707+
histos.fill(HIST("commonZna"), commonSumZna);
708+
}
709+
if ((tZNC >= minTdcZn) && (tZNC <= maxTdcZn)) {
696710
histos.fill(HIST("ampZnc"), znC);
697-
if ((tZPA >= minTdcZp) && (tZPA <= maxTdcZp))
711+
histos.fill(HIST("commonZnc"), commonSumZnc);
712+
}
713+
if ((tZPA >= minTdcZp) && (tZPA <= maxTdcZp)) {
698714
histos.fill(HIST("ampZpa"), zpA);
699-
if ((tZPC >= minTdcZp) && (tZPC <= maxTdcZp))
715+
histos.fill(HIST("commonZpa"), commonSumZpa);
716+
}
717+
if ((tZPC >= minTdcZp) && (tZPC <= maxTdcZp)) {
700718
histos.fill(HIST("ampZpc"), zpC);
719+
histos.fill(HIST("commonZpc"), commonSumZpc);
720+
}
701721
if (((tZNC >= minTdcZn) && (tZNC <= maxTdcZn)) && ((tZNA >= minTdcZn) && (tZNA <= maxTdcZn)))
702722
histos.fill(HIST("ZnVsZem"), sumZEMs, znC + znA);
703723
if (((tZNC >= minTdcZn) && (tZNC <= maxTdcZn)) && ((tZNA >= minTdcZn) && (tZNA <= maxTdcZn)))
@@ -713,6 +733,10 @@ struct FlowZdcTask {
713733
histos.fill(HIST("ampZnc"), znC);
714734
histos.fill(HIST("ampZpa"), zpA);
715735
histos.fill(HIST("ampZpc"), zpC);
736+
histos.fill(HIST("commonZna"), commonSumZna);
737+
histos.fill(HIST("commonZnc"), commonSumZnc);
738+
histos.fill(HIST("commonZpa"), commonSumZpa);
739+
histos.fill(HIST("commonZpc"), commonSumZpc);
716740
histos.fill(HIST("ZnVsZem"), sumZEMs, znC + znA);
717741
histos.fill(HIST("ZnaVsZnc"), znA, znC);
718742
histos.fill(HIST("ZpaVsZpc"), zpA, zpC);

0 commit comments

Comments
 (0)