Skip to content

Commit e3b4880

Browse files
omvazqueomvazque
authored
[PWGLF] Additional histograms to calculate ZP over ZN v3 (#11430)
1 parent 675e6ba commit e3b4880

File tree

1 file changed

+55
-29
lines changed

1 file changed

+55
-29
lines changed

‎PWGLF/Tasks/GlobalEventProperties/uccZdc.cxx‎

Lines changed: 55 additions & 29 deletions
Original file line numberDiff line numberDiff line change
@@ -65,8 +65,8 @@ using SimTracks = soa::Join<aod::Tracks, aod::TrackSelection, aod::TracksExtra,
6565

6666
struct UccZdc {
6767

68-
static constexpr float collEnergy{2.68};
69-
static constexpr float zEro{0.};
68+
static constexpr float kCollEnergy{2.68};
69+
static constexpr float kZero{0.};
7070

7171
// Configurables Event Selection
7272
Configurable<bool> isNoCollInTimeRangeStrict{"isNoCollInTimeRangeStrict", true, "use isNoCollInTimeRangeStrict?"};
@@ -101,13 +101,19 @@ struct UccZdc {
101101
Configurable<float> maxEta{"maxEta", +0.8, "maximum eta"};
102102

103103
// Configurables, binning
104+
Configurable<int> nBinsITSTrack{"nBinsITSTrack", 2000, "N bins ITS tracks"};
105+
Configurable<float> minITSTrack{"minITSTrack", 0., "Min ITS tracks"};
106+
Configurable<float> maxITSTrack{"maxITSTrack", 6000., "Min ITS tracks"};
104107
Configurable<float> maxAmpFV0{"maxAmpFV0", 2000, "Max FV0 amp"};
105108
Configurable<int> nBinsAmpFT0{"nBinsAmpFT0", 100, "N bins FT0 amp"};
109+
Configurable<int> nBinsAmpFT0Fine{"nBinsAmpFT0Fine", 1000, "N bins FT0 amp"};
106110
Configurable<float> maxAmpFT0{"maxAmpFT0", 2500, "Max FT0 amp"};
107111
Configurable<int> nBinsNch{"nBinsNch", 2501, "N bins Nch (|eta|<0.8)"};
112+
Configurable<int> nBinsNchFine{"nBinsNchFine", 3000, "N bins Nch (|eta|<0.8)"};
108113
Configurable<float> minNch{"minNch", 0, "Min Nch (|eta|<0.8)"};
109-
Configurable<float> maxNch{"maxNch", 2500, "Max Nch (|eta|<0.8)"};
110-
Configurable<int> nBinsZDC{"nBinsZDC", 400, "N bins ZDC"};
114+
Configurable<float> maxNch{"maxNch", 3000, "Max Nch (|eta|<0.8)"};
115+
Configurable<int> nBinsZN{"nBinsZN", 400, "N bins ZN"};
116+
Configurable<int> nBinsZP{"nBinsZP", 160, "N bins ZP"};
111117
Configurable<float> minZN{"minZN", 0, "Min ZN signal"};
112118
Configurable<float> maxZN{"maxZN", 150, "Max ZN signal"};
113119
Configurable<float> maxZP{"maxZP", 60, "Max ZP signal"};
@@ -171,7 +177,7 @@ struct UccZdc {
171177
registry.add("T0Ccent", ";;Entries", kTH1F, {axisCent});
172178
registry.add("NchUncorrected", ";#it{N}_{ch} (|#eta| < 0.8);Entries;", kTH1F, {{300, 0., 3000.}});
173179
registry.add("hEventCounter", ";;Events", kTH1F, {axisEvent});
174-
registry.add("ZNamp", ";ZNA+ZNC;Entries;", kTH1F, {{nBinsZDC, -0.5, maxZN}});
180+
registry.add("ZNamp", ";ZNA+ZNC;Entries;", kTH1F, {{nBinsZN, -0.5, maxZN}});
175181
registry.add("ExcludedEvtVsFT0M", ";T0A+T0C (#times 1/100, -3.3 < #eta < -2.1 and 3.5 < #eta < 4.9);Entries;", kTH1F, {{nBinsAmpFT0, 0., maxAmpFT0}});
176182
registry.add("ExcludedEvtVsNch", ";#it{N}_{ch} (|#eta|<0.8);Entries;", kTH1F, {{300, 0, 3000}});
177183
registry.add("Nch", ";#it{N}_{ch} (|#eta| < 0.8, Corrected);", kTH1F, {{nBinsNch, minNch, maxNch}});
@@ -202,11 +208,17 @@ struct UccZdc {
202208
x->SetBinLabel(17, "Within ZEM cut?");
203209

204210
if (doprocessZdcCollAss) {
205-
registry.add("NchVsPt", ";#it{N}_{ch} (|#eta| < 0.8, Corrected);;", kTH2F, {{{nBinsNch, minNch, maxNch}, {axisPt}}});
206-
registry.add("NchVsOneParCorrVsZN", ";#it{N}_{ch} (|#eta| < 0.8, Corrected); ZNA+ZNC; #LT[#it{p}_{T}^{(1)}]#GT", kTProfile2D, {{{nBinsNch, minNch, maxNch}, {nBinsZDC, -0.5, maxZN}}});
207-
registry.add("NchVsTwoParCorrVsZN", ";#it{N}_{ch} (|#eta| < 0.8, Corrected);ZNA+ZNC;#LT[#it{p}_{T}^{(2)}]#GT", kTProfile2D, {{{nBinsNch, minNch, maxNch}, {nBinsZDC, -0.5, maxZN}}});
208-
registry.add("NchVsThreeParCorrVsZN", ";#it{N}_{ch} (|#eta| < 0.8, Corrected);ZNA+ZNC;#LT[#it{p}_{T}^{(3)}]#GT", kTProfile2D, {{{nBinsNch, minNch, maxNch}, {nBinsZDC, -0.5, maxZN}}});
209-
registry.add("NchVsFourParCorrVsZN", ";#it{N}_{ch} (|#eta| < 0.8, Corrected);ZNA+ZNC;#LT[#it{p}_{T}^{(4)}]#GT", kTProfile2D, {{{nBinsNch, minNch, maxNch}, {nBinsZDC, -0.5, maxZN}}});
211+
registry.add("NchVsZN", ";#it{N}_{ch} (|#eta| < 0.8); ZNA+ZNC;", kTH2F, {{{nBinsNchFine, minNch, maxNch}, {nBinsZN, -0.5, maxZN}}});
212+
registry.add("NchVsZP", ";#it{N}_{ch} (|#eta| < 0.8); ZPA+ZPC;", kTH2F, {{{nBinsNchFine, minNch, maxNch}, {nBinsZP, -0.5, maxZP}}});
213+
registry.add("NITSTacksVsZN", ";ITS tracks; ZNA+ZNC;", kTH2F, {{{nBinsITSTrack, minITSTrack, maxITSTrack}, {nBinsZN, -0.5, maxZN}}});
214+
registry.add("NITSTacksVsZP", ";ITS tracks; ZPA+ZPC;", kTH2F, {{{nBinsITSTrack, minITSTrack, maxITSTrack}, {nBinsZP, -0.5, maxZP}}});
215+
registry.add("T0MVsZN", ";T0A+T0C amp (#times 1/100); ZNA+ZNC;", kTH2F, {{{nBinsAmpFT0Fine, 0., maxAmpFT0}, {nBinsZN, -0.5, maxZN}}});
216+
registry.add("T0MVsZP", ";T0A+T0C amp (#times 1/100); ZPA+ZPC;", kTH2F, {{{nBinsAmpFT0Fine, 0., maxAmpFT0}, {nBinsZP, -0.5, maxZP}}});
217+
registry.add("NchVsZNVsPt", ";#it{N}_{ch} (|#eta| < 0.8); ZNA+ZNC;#it{p}_{T} (GeV/#it{c})", kTH3F, {{{nBinsNch, minNch, maxNch}, {nBinsZN, -0.5, maxZN}, {axisPt}}});
218+
registry.add("NchVsOneParCorrVsZN", ";#it{N}_{ch} (|#eta| < 0.8, Corrected); ZNA+ZNC; #LT[#it{p}_{T}^{(1)}]#G (GeV/#it{c})T", kTProfile2D, {{{nBinsNch, minNch, maxNch}, {nBinsZN, -0.5, maxZN}}});
219+
registry.add("NchVsTwoParCorrVsZN", ";#it{N}_{ch} (|#eta| < 0.8, Corrected);ZNA+ZNC;#LT[#it{p}_{T}^{(2)}]#GT", kTProfile2D, {{{nBinsNch, minNch, maxNch}, {nBinsZN, -0.5, maxZN}}});
220+
registry.add("NchVsThreeParCorrVsZN", ";#it{N}_{ch} (|#eta| < 0.8, Corrected);ZNA+ZNC;#LT[#it{p}_{T}^{(3)}]#GT", kTProfile2D, {{{nBinsNch, minNch, maxNch}, {nBinsZN, -0.5, maxZN}}});
221+
registry.add("NchVsFourParCorrVsZN", ";#it{N}_{ch} (|#eta| < 0.8, Corrected);ZNA+ZNC;#LT[#it{p}_{T}^{(4)}]#GT", kTProfile2D, {{{nBinsNch, minNch, maxNch}, {nBinsZN, -0.5, maxZN}}});
210222
}
211223

212224
if (doprocessMCclosure) {
@@ -260,21 +272,21 @@ struct UccZdc {
260272
registry.add("NchVsEt", ";#it{E}_{T} (|#eta|<0.8);#LTITS+TPC tracks#GT (|#eta|<0.8);", kTH2F, {{{nBinsNch, minNch, maxNch}, {nBinsNch, minNch, maxNch}}});
261273
registry.add("NchVsNPV", ";#it{N}_{PV} (|#eta|<1);ITS+TPC tracks (|#eta|<0.8);", kTH2F, {{{300, -0.5, 5999.5}, {nBinsNch, minNch, maxNch}}});
262274
registry.add("NchVsITStracks", ";ITS tracks nCls >= 5;TITS+TPC tracks (|#eta|<0.8);", kTH2F, {{{300, -0.5, 5999.5}, {nBinsNch, minNch, maxNch}}});
263-
registry.add("ZNVsFT0A", ";T0A (#times 1/100);ZNA+ZNC;", kTH2F, {{{80, 0., 1800.}, {nBinsZDC, -0.5, maxZN}}});
264-
registry.add("ZNVsFT0C", ";T0C (#times 1/100);ZNA+ZNC;", kTH2F, {{{80, 0., 600.}, {nBinsZDC, -0.5, maxZN}}});
265-
registry.add("ZNVsFT0M", ";T0A+T0C (#times 1/100);ZNA+ZNC;", kTH2F, {{{nBinsAmpFT0, 0., maxAmpFT0}, {nBinsZDC, -0.5, maxZN}}});
266-
registry.add("ZNAamp", ";ZNA;Entries;", kTH1F, {{nBinsZDC, -0.5, maxZN}});
267-
registry.add("ZPAamp", ";ZPA;Entries;", kTH1F, {{nBinsZDC, -0.5, maxZP}});
268-
registry.add("ZNCamp", ";ZNC;Entries;", kTH1F, {{nBinsZDC, -0.5, maxZN}});
269-
registry.add("ZPCamp", ";ZPC;Entries;", kTH1F, {{nBinsZDC, -0.5, maxZP}});
275+
registry.add("ZNVsFT0A", ";T0A (#times 1/100);ZNA+ZNC;", kTH2F, {{{80, 0., 1800.}, {nBinsZN, -0.5, maxZN}}});
276+
registry.add("ZNVsFT0C", ";T0C (#times 1/100);ZNA+ZNC;", kTH2F, {{{80, 0., 600.}, {nBinsZN, -0.5, maxZN}}});
277+
registry.add("ZNVsFT0M", ";T0A+T0C (#times 1/100);ZNA+ZNC;", kTH2F, {{{nBinsAmpFT0, 0., maxAmpFT0}, {nBinsZN, -0.5, maxZN}}});
278+
registry.add("ZNAamp", ";ZNA;Entries;", kTH1F, {{nBinsZN, -0.5, maxZN}});
279+
registry.add("ZPAamp", ";ZPA;Entries;", kTH1F, {{nBinsZP, -0.5, maxZP}});
280+
registry.add("ZNCamp", ";ZNC;Entries;", kTH1F, {{nBinsZN, -0.5, maxZN}});
281+
registry.add("ZPCamp", ";ZPC;Entries;", kTH1F, {{nBinsZP, -0.5, maxZP}});
270282
registry.add("ZNAVsZNC", ";ZNC;ZNA", kTH2F, {{{30, -0.5, maxZN}, {30, -0.5, maxZN}}});
271283
registry.add("ZPAVsZPC", ";ZPC;ZPA;", kTH2F, {{{100, -0.5, maxZP}, {100, -0.5, maxZP}}});
272284
registry.add("ZNAVsZPA", ";ZPA;ZNA;", kTH2F, {{{20, -0.5, maxZP}, {30, -0.5, maxZN}}});
273285
registry.add("ZNCVsZPC", ";ZPC;ZNC;", kTH2F, {{{20, -0.5, maxZP}, {30, -0.5, maxZN}}});
274286
registry.add("ZNVsZEM", ";ZEM;ZNA+ZNC;", kTH2F, {{{60, -0.5, maxZEM}, {60, -0.5, maxZN}}});
275-
registry.add("ZNCVsNch", ";#it{N}_{ch} (|#eta|<0.8);ZNC;", kTH2F, {{{nBinsNch, minNch, maxNch}, {nBinsZDC, minZN, maxZN}}});
276-
registry.add("ZNAVsNch", ";#it{N}_{ch} (|#eta|<0.8);ZNA;", kTH2F, {{{nBinsNch, minNch, maxNch}, {nBinsZDC, minZN, maxZN}}});
277-
registry.add("ZNVsNch", ";#it{N}_{ch} (|#eta|<0.8);ZNA+ZNC;", kTH2F, {{{nBinsNch, minNch, maxNch}, {nBinsZDC, minZN, maxZN}}});
287+
registry.add("ZNCVsNch", ";#it{N}_{ch} (|#eta|<0.8);ZNC;", kTH2F, {{{nBinsNch, minNch, maxNch}, {nBinsZN, minZN, maxZN}}});
288+
registry.add("ZNAVsNch", ";#it{N}_{ch} (|#eta|<0.8);ZNA;", kTH2F, {{{nBinsNch, minNch, maxNch}, {nBinsZN, minZN, maxZN}}});
289+
registry.add("ZNVsNch", ";#it{N}_{ch} (|#eta|<0.8);ZNA+ZNC;", kTH2F, {{{nBinsNch, minNch, maxNch}, {nBinsZN, minZN, maxZN}}});
278290
registry.add("ZNDifVsNch", ";#it{N}_{ch} (|#eta|<0.8);ZNA-ZNC;", kTH2F, {{{nBinsNch, minNch, maxNch}, {100, -50., 50.}}});
279291
}
280292

@@ -429,10 +441,10 @@ struct UccZdc {
429441
float tZPC{zdc.timeZPC()};
430442
float tZDCdif{tZNC + tZPC - tZNA - tZPA};
431443
float tZDCsum{tZNC + tZPC + tZNA + tZPA};
432-
znA /= collEnergy;
433-
znC /= collEnergy;
434-
zpA /= collEnergy;
435-
zpC /= collEnergy;
444+
znA /= kCollEnergy;
445+
znC /= kCollEnergy;
446+
zpA /= kCollEnergy;
447+
zpC /= kCollEnergy;
436448
float sumZNs{znA + znC};
437449
float sumZEMs{aZEM1 + aZEM2};
438450

@@ -576,6 +588,8 @@ struct UccZdc {
576588
const double normT0M{(aT0A + aT0C) / 100.};
577589
float znA{foundBC.zdc().amplitudeZNA()};
578590
float znC{foundBC.zdc().amplitudeZNC()};
591+
float zpA{foundBC.zdc().amplitudeZPA()};
592+
float zpC{foundBC.zdc().amplitudeZPC()};
579593
float aZEM1{foundBC.zdc().amplitudeZEM1()};
580594
float aZEM2{foundBC.zdc().amplitudeZEM2()};
581595
float tZNA{foundBC.zdc().timeZNA()};
@@ -584,9 +598,12 @@ struct UccZdc {
584598
float tZPC{foundBC.zdc().timeZPC()};
585599
float tZDCdif{tZNC + tZPC - tZNA - tZPA};
586600
float tZDCsum{tZNC + tZPC + tZNA + tZPA};
587-
znA /= collEnergy;
588-
znC /= collEnergy;
601+
znA /= kCollEnergy;
602+
znC /= kCollEnergy;
603+
zpA /= kCollEnergy;
604+
zpC /= kCollEnergy;
589605
float sumZNs{znA + znC};
606+
float sumZPs{zpA + zpC};
590607
float sumZEMs{aZEM1 + aZEM2};
591608

592609
// TDC cut
@@ -609,9 +626,12 @@ struct UccZdc {
609626
registry.fill(HIST("T0Ccent"), collision.centFT0C());
610627

611628
// Nch-based selection
612-
int glbTracks{0};
629+
int itsTracks = 0, glbTracks = 0;
613630
for (const auto& track : tracks) {
614631
// Track Selection
632+
if (track.hasITS()) {
633+
itsTracks++;
634+
}
615635
if (!track.isGlobalTrack()) {
616636
continue;
617637
}
@@ -693,7 +713,7 @@ struct UccZdc {
693713
if (!track.isGlobalTrack()) {
694714
continue;
695715
}
696-
registry.fill(HIST("NchVsPt"), w1, track.pt());
716+
registry.fill(HIST("NchVsZNVsPt"), w1, sumZNs, track.pt());
697717
}
698718

699719
// EbE one-particle pT correlation
@@ -716,6 +736,12 @@ struct UccZdc {
716736

717737
registry.fill(HIST("Nch"), w1);
718738
registry.fill(HIST("ZNamp"), sumZNs);
739+
registry.fill(HIST("NchVsZN"), w1, sumZNs);
740+
registry.fill(HIST("NchVsZP"), w1, sumZPs);
741+
registry.fill(HIST("NITSTacksVsZN"), itsTracks, sumZNs);
742+
registry.fill(HIST("NITSTacksVsZP"), itsTracks, sumZPs);
743+
registry.fill(HIST("T0MVsZN"), normT0M, sumZNs);
744+
registry.fill(HIST("T0MVsZP"), normT0M, sumZPs);
719745
registry.fill(HIST("NchUncorrected"), glbTracks);
720746
registry.fill(HIST("NchVsOneParCorr"), w1, oneParCorr, w1);
721747
registry.fill(HIST("NchVsOneParCorrVsZN"), w1, sumZNs, oneParCorr, w1);
@@ -760,7 +786,7 @@ struct UccZdc {
760786
registry.fill(HIST("T0Ccent"), cent);
761787

762788
// Half of the statistics for MC closure
763-
if (rndNum >= zEro && rndNum < evtFracMCcl) {
789+
if (rndNum >= kZero && rndNum < evtFracMCcl) {
764790
registry.fill(HIST("EvtsDivided"), 0);
765791

766792
// To use run-by-run efficiency

0 commit comments

Comments
 (0)