Skip to content

Commit d69766c

Browse files
author
Sabrina Hernandez
committed
added pZDC correlations with ft0c and glb tracks as cent estimators
1 parent eaf70c7 commit d69766c

1 file changed

Lines changed: 47 additions & 4 deletions

File tree

‎PWGCF/Flow/Tasks/flowZdcTask.cxx‎

Lines changed: 47 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -54,6 +54,9 @@ using namespace o2::aod::evsel;
5454

5555
struct FlowZdcTask {
5656
SliceCache cache;
57+
Configurable<bool> useFT0M{"useFT0M", true, "FT0M for centrality"};
58+
Configurable<bool> useFT0C{"useFT0C", false, "FT0C for centrality"};
59+
Configurable<bool> useGlbTracks{"useGlbTracks", false, "GlbTracks for centrality"};
5760
Configurable<float> cfgCutVertex{"cfgCutVertex", 10.0, "Accepted z-vertex range"};
5861
Configurable<float> maxZem{"maxZem", 3099.5, "Max ZEM signal"};
5962
// for ZDC info and analysis
@@ -65,6 +68,8 @@ struct FlowZdcTask {
6568
Configurable<float> maxAmpFT0{"maxAmpFT0", 2500, "Max FT0 amp"};
6669
Configurable<float> maxAmpFT0M{"maxAmpFT0M", 2500, "Max FT0M amp"};
6770
Configurable<int> nBinsAmpFT0M{"nBinsAmpFT0M", 100, "N bins FT0M amp"};
71+
Configurable<int> nBinsAmpFT0C{"nBinsAmpFT0C", 100, "N bins FT0C amp"};
72+
Configurable<int> nBinsGlbTracks{"nBinsGlbTracks", 100, "N bins global tracks"};
6873
Configurable<int> nBinsZDC{"nBinsZDC", 400, "nBinsZDC"};
6974
Configurable<int> nBinsZP{"nBinsZP", 50, "nBinsZP"};
7075
Configurable<int> nBinsZN{"nBinsZN", 50, "nBinsZN"};
@@ -199,6 +204,18 @@ struct FlowZdcTask {
199204
xAxis->SetBinLabel(16, "Within TDC cut?");
200205

201206
if (doprocessQA) {
207+
if (useFT0M) {
208+
histos.add("ZPAZNAVsFT0M", ";T0A+T0C (#times 1/100);ZPA Amplitude;ZNA Amplitude;", kTH3F, {{{nBinsAmpFT0M, 0., maxAmpFT0M}, {nBinsZP, -0.5, maxZp}, {nBinsZN, -0.5, maxZn}}});
209+
histos.add("ZPCZNCVsFT0M", ";T0A+T0C (#times 1/100);ZPC Amplitude;ZNC Amplitude;", kTH3F, {{{nBinsAmpFT0M, 0., maxAmpFT0M}, {nBinsZP, -0.5, maxZp}, {nBinsZN, -0.5, maxZn}}});
210+
}
211+
if (useFT0C) {
212+
histos.add("ZPAZNAVsFT0C", ";T0C (#times 1/100);ZPA Amplitude;ZNA Amplitude;", kTH3F, {{{nBinsAmpFT0C, 0., maxAmpFT0}, {nBinsZP, -0.5, maxZp}, {nBinsZN, -0.5, maxZn}}});
213+
histos.add("ZPCZNCVsFT0C", ";T0C (#times 1/100);ZPC Amplitude;ZNC Amplitude;", kTH3F, {{{nBinsAmpFT0C, 0., maxAmpFT0}, {nBinsZP, -0.5, maxZp}, {nBinsZN, -0.5, maxZn}}});
214+
}
215+
if (useGlbTracks) {
216+
histos.add("ZPAZNAVsGlbTracks", ";N_{global tracks} (|#eta|<0.8);ZPA Amplitude;ZNA Amplitude;", kTH3F, {{{nBinsGlbTracks, minNch, maxNch}, {nBinsZP, -0.5, maxZp}, {nBinsZN, -0.5, maxZn}}});
217+
histos.add("ZPCZNCVsGlbTracks", ";N_{global tracks} (|#eta|<0.8);ZPC Amplitude;ZNC Amplitude;", kTH3F, {{{nBinsGlbTracks, minNch, maxNch}, {nBinsZP, -0.5, maxZp}, {nBinsZN, -0.5, maxZn}}});
218+
}
202219
histos.add("ZNAVsFT0C", ";T0C (#times 1/100);ZNA Amplitude;", kTH2F, {{{nBinsAmpFT0, 0., maxAmpFT0}, {nBinsZDC, -0.5, maxZn}}});
203220
histos.add("ZNAVsFT0M", ";T0A+T0C (#times 1/100);ZNA Amplitude;", kTH2F, {{{nBinsAmpFT0, 0., maxAmpFT0M}, {nBinsZDC, -0.5, maxZn}}});
204221
histos.add("ZNCVsFT0C", ";T0C (#times 1/100);ZNC Amplitude;", kTH2F, {{{nBinsAmpFT0, 0., maxAmpFT0}, {nBinsZDC, -0.5, maxZn}}});
@@ -545,12 +562,28 @@ struct FlowZdcTask {
545562
if (((tZNA >= minTdcZn) && (tZNA <= maxTdcZn)) && ((tZPA >= minTdcZp) && (tZPA <= maxTdcZp))) {
546563
histos.fill(HIST("ZNAVsZPA"), zpA, znA);
547564
histos.fill(HIST("CommonZNAVsZPA"), commonSumZpa, commonSumZna);
548-
histos.fill(HIST("ZPAZNAVsFT0M"), (aT0A + aT0C) / 100., zpA, znA);
565+
if (useFT0M) {
566+
histos.fill(HIST("ZPAZNAVsFT0M"), (aT0A + aT0C) / 100., zpA, znA);
567+
}
568+
if (useFT0C) {
569+
histos.fill(HIST("ZPAZNAVsFT0C"), aT0C / 100., zpA, znA);
570+
}
571+
if (useGlbTracks) {
572+
histos.fill(HIST("ZPAZNAVsGlbTracks"), glbTracks, zpA, znA);
573+
}
549574
}
550575
if (((tZNC >= minTdcZn) && (tZNC <= maxTdcZn)) && ((tZPC >= minTdcZp) && (tZPC <= maxTdcZp))) {
551576
histos.fill(HIST("ZNCVsZPC"), zpC, znC);
552577
histos.fill(HIST("CommonZNCVsZPC"), commonSumZpc, commonSumZnc);
553-
histos.fill(HIST("ZPCZNCVsFT0M"), (aT0A + aT0C) / 100., zpC, znC);
578+
if (useFT0M) {
579+
histos.fill(HIST("ZPCZNCVsFT0M"), (aT0A + aT0C) / 100., zpC, znC);
580+
}
581+
if (useFT0C) {
582+
histos.fill(HIST("ZPCZNCVsFT0C"), aT0C / 100., zpC, znC);
583+
}
584+
if (useGlbTracks) {
585+
histos.fill(HIST("ZPCZNCVsGlbTracks"), glbTracks, zpC, znC);
586+
}
554587
}
555588
} else {
556589
histos.fill(HIST("ZNA"), znA);
@@ -583,8 +616,18 @@ struct FlowZdcTask {
583616
histos.fill(HIST("ZPCVsNch"), glbTracks, zpC);
584617
histos.fill(HIST("ZNCVsNch"), glbTracks, znC);
585618
histos.fill(HIST("ZNAVsNch"), glbTracks, znA);
586-
histos.fill(HIST("ZPAZNAVsFT0M"), (aT0A + aT0C) / 100., zpA, znA);
587-
histos.fill(HIST("ZPCZNCVsFT0M"), (aT0A + aT0C) / 100., zpC, znC);
619+
if (useFT0M) {
620+
histos.fill(HIST("ZPAZNAVsFT0M"), (aT0A + aT0C) / 100., zpA, znA);
621+
histos.fill(HIST("ZPCZNCVsFT0M"), (aT0A + aT0C) / 100., zpC, znC);
622+
}
623+
if (useFT0C) {
624+
histos.fill(HIST("ZPAZNAVsFT0C"), aT0C / 100., zpA, znA);
625+
histos.fill(HIST("ZPCZNCVsFT0C"), aT0C / 100., zpC, znC);
626+
}
627+
if (useGlbTracks) {
628+
histos.fill(HIST("ZPAZNAVsGlbTracks"), glbTracks, zpA, znA);
629+
histos.fill(HIST("ZPCZNCVsGlbTracks"), glbTracks, zpC, znC);
630+
}
588631
}
589632
histos.fill(HIST("ZEM1"), aZEM1);
590633
histos.fill(HIST("ZEM2"), aZEM2);

0 commit comments

Comments
 (0)