Skip to content

Commit d970d35

Browse files
fuchuncuifuchuncui
authored
PWGCF: Add variable mass axis in flowGFWOmegaXi.cxx (AliceO2Group#7977)
* Add files via upload * Add files via upload * Add files via upload
1 parent d1237bb commit d970d35

File tree

1 file changed

+85
-57
lines changed

1 file changed

+85
-57
lines changed

‎PWGCF/Flow/Tasks/flowGFWOmegaXi.cxx‎

Lines changed: 85 additions & 57 deletions
Original file line numberDiff line numberDiff line change
@@ -40,6 +40,7 @@
4040
#include <TRandom3.h>
4141
#include <TF1.h>
4242
#include <TF2.h>
43+
#include <TPDGCode.h>
4344

4445
using namespace o2;
4546
using namespace o2::framework;
@@ -88,6 +89,8 @@ struct FlowGFWOmegaXi {
8889
ConfigurableAxis cfgaxisVertex{"axisVertex", {20, -10, 10}, "vertex axis for histograms"};
8990
ConfigurableAxis cfgaxisPhi{"axisPhi", {60, 0.0, constants::math::TwoPI}, "phi axis for histograms"};
9091
ConfigurableAxis cfgaxisEta{"axisEta", {40, -1., 1.}, "eta axis for histograms"};
92+
ConfigurableAxis cfgaxisPtXi{"axisPtXi", {VARIABLE_WIDTH, 0.9, 1.1, 1.3, 1.5, 1.7, 1.9, 2.1, 2.3, 2.5, 2.7, 2.9, 3.9, 4.9, 5.9, 9.9}, "pt (GeV)"};
93+
ConfigurableAxis cfgaxisPtV0{"axisPtV0", {VARIABLE_WIDTH, 0.9, 1.1, 1.3, 1.5, 1.7, 1.9, 2.1, 2.3, 2.5, 2.7, 2.9, 3.9, 4.9, 5.9, 9.9}, "pt (GeV)"};
9194

9295
Configurable<std::vector<float>> cfgvecwacc{"vecwacc", std::vector<float>{0.879543, 0.893808, 0.993375, 1.09663, 0.983883, 0.984094, 1.11362, 0.963896, 0.911212, 1.02934, 1.00295, 0.950711, 0.996856, 1.11934, 0.993665, 0.99087, 1.11915, 1.0198, 0.966849, 1.03237, 0.989367, 0.948312, 0.970883, 0.984305, 0.920335, 0.929722, 1.07467, 1.00862, 0.977185, 0.870868, 1.06552, 0.962393, 1.01025, 1.09959, 0.984226, 0.986361, 1.0931, 0.994377, 0.976051, 1.05249, 0.995538, 0.886452, 0.936763, 0.993613, 0.94491, 0.966559, 1.10829, 1.01998, 0.991503, 1.07918, 1.05655, 0.973784, 1.00914, 1.11678, 1.00092, 0.95232, 1.09814, 1.02322, 0.958543, 0.947231}, "wacc in phi bins"};
9396
AxisSpec axisPt{{0.20, 0.25, 0.30, 0.35, 0.40, 0.45, 0.50, 0.55, 0.60, 0.65, 0.70, 0.75, 0.80, 0.85, 0.90, 0.95, 1.00, 1.10, 1.20, 1.30, 1.40, 1.50, 1.60, 1.70, 1.80, 1.90, 2.00, 2.20, 2.40, 2.60, 2.80, 3.00, 3.50, 4.00, 4.50, 5.00, 5.50, 6.00, 7.00, 8.00, 9.00, 10.0}, "pt(GeV)"};
@@ -96,6 +99,10 @@ struct FlowGFWOmegaXi {
9699
AxisSpec axisXiminusMass = {70, 1.3f, 1.37f, "Inv. Mass (GeV)"};
97100
AxisSpec axisK0sMass = {400, 0.4f, 0.6f, "Inv. Mass (GeV)"};
98101
AxisSpec axisLambdaMass = {160, 1.08f, 1.16f, "Inv. Mass (GeV)"};
102+
AxisSpec axisOmegaminusMassforflow = {cfgOmegaMassbins, 1.63f, 1.71f, "Inv. Mass (GeV)"};
103+
AxisSpec axisXiminusMassforflow = {cfgXiMassbins, 1.3f, 1.37f, "Inv. Mass (GeV)"};
104+
AxisSpec axisK0sMassforflow = {cfgK0sMassbins, 0.4f, 0.6f, "Inv. Mass (GeV)"};
105+
AxisSpec axisLambdaMassforflow = {cfgLambdaMassbins, 1.08f, 1.16f, "Inv. Mass (GeV)"};
99106

100107
Configurable<bool> cfgcheckDauTPC{"checkDauTPC", false, "check if daughter tracks have TPC match"};
101108
Configurable<float> cfgCasc_rapidity{"Casc_rapidity", 0.5, "rapidity"};
@@ -138,11 +145,20 @@ struct FlowGFWOmegaXi {
138145
using DaughterTracks = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::pidTPCPi, aod::pidTPCPr, aod::pidTPCKa>;
139146

140147
// Set the pt, mult and phi Axis;
141-
o2::framework::AxisSpec axis = axisPt;
142-
int nPtBins = axis.binEdges.size() - 1;
143-
double* PtBins = &(axis.binEdges)[0];
148+
int nPtBins = axisPt.binEdges.size() - 1;
149+
double* PtBins = &(axisPt.binEdges)[0];
144150
TAxis* fPtAxis = new TAxis(nPtBins, PtBins);
145151

152+
o2::framework::AxisSpec axisXiPt = cfgaxisPtXi;
153+
int nXiPtBins = axisXiPt.binEdges.size() - 1;
154+
double* XiPtBins = &(axisXiPt.binEdges)[0];
155+
TAxis* fXiPtAxis = new TAxis(nXiPtBins, XiPtBins);
156+
157+
o2::framework::AxisSpec axisV0Pt = cfgaxisPtV0;
158+
int nV0PtBins = axisV0Pt.binEdges.size() - 1;
159+
double* V0PtBins = &(axisV0Pt.binEdges)[0];
160+
TAxis* fV0PtAxis = new TAxis(nV0PtBins, V0PtBins);
161+
146162
o2::framework::AxisSpec axisMult = axisMultiplicity;
147163
int nMultBins = axisMult.binEdges.size() - 1;
148164
double* MultBins = &(axisMult.binEdges)[0];
@@ -176,40 +192,40 @@ struct FlowGFWOmegaXi {
176192
registry.add("hMult", "", {HistType::kTH1D, {{3000, 0.5, 3000.5}}});
177193
registry.add("hCent", "", {HistType::kTH1D, {{90, 0, 90}}});
178194
registry.add("hPt", "", {HistType::kTH1D, {axisPt}});
179-
registry.add("hEtaPhiREF", "", {HistType::kTH2D, {cfgaxisEta, cfgaxisPhi}});
180-
registry.add("hEtaPhiPOIXi", "", {HistType::kTH2D, {cfgaxisEta, cfgaxisPhi}});
181-
registry.add("hEtaPhiPOIOmega", "", {HistType::kTH2D, {cfgaxisEta, cfgaxisPhi}});
182-
registry.add("hEtaPhiPOIK0s", "", {HistType::kTH2D, {cfgaxisEta, cfgaxisPhi}});
183-
registry.add("hEtaPhiPOILambda", "", {HistType::kTH2D, {cfgaxisEta, cfgaxisPhi}});
195+
registry.add("hEtaPhiVtxzREF", "", {HistType::kTH3D, {cfgaxisEta, cfgaxisPhi, {20, -10, 10}}});
196+
registry.add("hEtaPhiVtxzPOIXi", "", {HistType::kTH3D, {cfgaxisEta, cfgaxisPhi, {20, -10, 10}}});
197+
registry.add("hEtaPhiVtxzPOIOmega", "", {HistType::kTH3D, {cfgaxisEta, cfgaxisPhi, {20, -10, 10}}});
198+
registry.add("hEtaPhiVtxzPOIK0s", "", {HistType::kTH3D, {cfgaxisEta, cfgaxisPhi, {20, -10, 10}}});
199+
registry.add("hEtaPhiVtxzPOILambda", "", {HistType::kTH3D, {cfgaxisEta, cfgaxisPhi, {20, -10, 10}}});
184200
// cumulant of flow
185201
registry.add("c22", ";Centrality (%) ; C_{2}{2}", {HistType::kTProfile, {axisMultiplicity}});
186202
registry.add("c24", ";Centrality (%) ; C_{2}{4}", {HistType::kTProfile, {axisMultiplicity}});
187203
// pt-diff cumulant of flow
188-
registry.add("Xic22dpt", ";pt ; C_{2}{2} ", {HistType::kTProfile3D, {axisPt, axisXiminusMass, axisMultiplicity}});
189-
registry.add("Omegac22dpt", ";pt ; C_{2}{2} ", {HistType::kTProfile3D, {axisPt, axisOmegaminusMass, axisMultiplicity}});
190-
registry.add("K0sc22dpt", ";pt ; C_{2}{2} ", {HistType::kTProfile3D, {axisPt, axisK0sMass, axisMultiplicity}});
191-
registry.add("Lambdac22dpt", ";pt ; C_{2}{2} ", {HistType::kTProfile3D, {axisPt, axisLambdaMass, axisMultiplicity}});
204+
registry.add("Xic22dpt", ";pt ; C_{2}{2} ", {HistType::kTProfile3D, {cfgaxisPtXi, axisXiminusMassforflow, axisMultiplicity}});
205+
registry.add("Omegac22dpt", ";pt ; C_{2}{2} ", {HistType::kTProfile3D, {cfgaxisPtXi, axisOmegaminusMassforflow, axisMultiplicity}});
206+
registry.add("K0sc22dpt", ";pt ; C_{2}{2} ", {HistType::kTProfile3D, {cfgaxisPtV0, axisK0sMassforflow, axisMultiplicity}});
207+
registry.add("Lambdac22dpt", ";pt ; C_{2}{2} ", {HistType::kTProfile3D, {cfgaxisPtV0, axisLambdaMassforflow, axisMultiplicity}});
192208
// InvMass(GeV) of casc and v0
193-
registry.add("InvMassXiMinus_all", "", {HistType::kTHnSparseF, {axisPt, axisXiminusMass, cfgaxisEta, axisMultiplicity}});
194-
registry.add("InvMassOmegaMinus_all", "", {HistType::kTHnSparseF, {axisPt, axisOmegaminusMass, cfgaxisEta, axisMultiplicity}});
195-
registry.add("InvMassOmegaMinus", "", {HistType::kTHnSparseF, {axisPt, axisOmegaminusMass, cfgaxisEta, axisMultiplicity}});
196-
registry.add("InvMassXiMinus", "", {HistType::kTHnSparseF, {axisPt, axisXiminusMass, cfgaxisEta, axisMultiplicity}});
197-
registry.add("InvMassK0s_all", "", {HistType::kTHnSparseF, {axisPt, axisK0sMass, cfgaxisEta, axisMultiplicity}});
198-
registry.add("InvMassLambda_all", "", {HistType::kTHnSparseF, {axisPt, axisLambdaMass, cfgaxisEta, axisMultiplicity}});
199-
registry.add("InvMassK0s", "", {HistType::kTHnSparseF, {axisPt, axisK0sMass, cfgaxisEta, axisMultiplicity}});
200-
registry.add("InvMassLambda", "", {HistType::kTHnSparseF, {axisPt, axisLambdaMass, cfgaxisEta, axisMultiplicity}});
209+
registry.add("InvMassXiMinus_all", "", {HistType::kTHnSparseF, {cfgaxisPtXi, axisXiminusMass, cfgaxisEta, axisMultiplicity}});
210+
registry.add("InvMassOmegaMinus_all", "", {HistType::kTHnSparseF, {cfgaxisPtXi, axisOmegaminusMass, cfgaxisEta, axisMultiplicity}});
211+
registry.add("InvMassOmegaMinus", "", {HistType::kTHnSparseF, {cfgaxisPtXi, axisOmegaminusMass, cfgaxisEta, axisMultiplicity}});
212+
registry.add("InvMassXiMinus", "", {HistType::kTHnSparseF, {cfgaxisPtXi, axisXiminusMass, cfgaxisEta, axisMultiplicity}});
213+
registry.add("InvMassK0s_all", "", {HistType::kTHnSparseF, {cfgaxisPtV0, axisK0sMass, cfgaxisEta, axisMultiplicity}});
214+
registry.add("InvMassLambda_all", "", {HistType::kTHnSparseF, {cfgaxisPtV0, axisLambdaMass, cfgaxisEta, axisMultiplicity}});
215+
registry.add("InvMassK0s", "", {HistType::kTHnSparseF, {cfgaxisPtV0, axisK0sMass, cfgaxisEta, axisMultiplicity}});
216+
registry.add("InvMassLambda", "", {HistType::kTHnSparseF, {cfgaxisPtV0, axisLambdaMass, cfgaxisEta, axisMultiplicity}});
201217

202218
fGFW->AddRegion("full", -0.8, 0.8, 1, 1); // ("name", etamin, etamax, ptbinnum, bitmask)eta region -0.8 to 0.8
203219
// with (-0.5, 0.5) eta gap
204220
fGFW->AddRegion("refN10", -0.8, -0.5, 1, 1);
205221
fGFW->AddRegion("refP10", 0.5, 0.8, 1, 1);
206-
int nXiptMassBins = nPtBins * nXiMassBins;
222+
int nXiptMassBins = nXiPtBins * nXiMassBins;
207223
fGFW->AddRegion("poiXiP", 0.5, 0.8, nXiptMassBins, 2);
208-
int nOmegaptMassBins = nPtBins * nOmegaMassBins;
224+
int nOmegaptMassBins = nXiPtBins * nOmegaMassBins;
209225
fGFW->AddRegion("poiOmegaP", 0.5, 0.8, nOmegaptMassBins, 4);
210-
int nK0sptMassBins = nPtBins * nK0sMassBins;
226+
int nK0sptMassBins = nV0PtBins * nK0sMassBins;
211227
fGFW->AddRegion("poiK0sP", 0.5, 0.8, nK0sptMassBins, 8);
212-
int nLambdaptMassBins = nPtBins * nLambdaMassBins;
228+
int nLambdaptMassBins = nV0PtBins * nLambdaMassBins;
213229
fGFW->AddRegion("poiLambdaP", 0.5, 0.8, nLambdaptMassBins, 16);
214230
// pushback
215231
corrconfigs.push_back(fGFW->GetCorrelatorConfig("refP10 {2} refN10 {-2}", "ChFull220", kFALSE));
@@ -255,35 +271,45 @@ struct FlowGFWOmegaXi {
255271
}
256272

257273
template <char... chars>
258-
void FillProfilepT(const GFW::CorrConfig& corrconf, const ConstStr<chars...>& tarName, const int& ptbin, const int& partical, const float& cent)
274+
void FillProfilepT(const GFW::CorrConfig& corrconf, const ConstStr<chars...>& tarName, const int& ptbin, const int& PDGCode, const float& cent)
259275
{
260276
int nMassBins = 0;
277+
int nptbins = nPtBins;
261278
TAxis* fMass = nullptr;
262-
if (partical == 3312) {
279+
TAxis* fpt = nullptr;
280+
if (PDGCode == kXiMinus) {
263281
nMassBins = nXiMassBins;
282+
nptbins = nXiPtBins;
283+
fpt = fXiPtAxis;
264284
fMass = fXiMass;
265-
} else if (partical == 3334) {
285+
} else if (PDGCode == kOmegaMinus) {
266286
nMassBins = nOmegaMassBins;
287+
nptbins = nXiPtBins;
288+
fpt = fXiPtAxis;
267289
fMass = fOmegaMass;
268-
} else if (partical == 310) {
290+
} else if (PDGCode == kK0Short) {
269291
nMassBins = nK0sMassBins;
292+
nptbins = nV0PtBins;
293+
fpt = fV0PtAxis;
270294
fMass = fK0sMass;
271-
} else if (partical == 3122) {
295+
} else if (PDGCode == kLambda0) {
272296
nMassBins = nLambdaMassBins;
297+
nptbins = nV0PtBins;
298+
fpt = fV0PtAxis;
273299
fMass = fLambdaMass;
274300
} else {
275-
LOGF(error, "Error, please put in correct PDG particalID of K0s, Lambda, Xi or Omega");
301+
LOGF(error, "Error, please put in correct PDGCode of K0s, Lambda, Xi or Omega");
276302
return;
277303
}
278304
for (int massbin = 1; massbin <= nMassBins; massbin++) {
279305
float dnx = 0;
280306
float val = 0;
281-
dnx = fGFW->Calculate(corrconf, (ptbin - 1) + ((massbin - 1) * nPtBins), kTRUE).real();
307+
dnx = fGFW->Calculate(corrconf, (ptbin - 1) + ((massbin - 1) * nptbins), kTRUE).real();
282308
if (dnx == 0)
283309
continue;
284-
val = fGFW->Calculate(corrconf, (ptbin - 1) + ((massbin - 1) * nPtBins), kFALSE).real() / dnx;
310+
val = fGFW->Calculate(corrconf, (ptbin - 1) + ((massbin - 1) * nptbins), kFALSE).real() / dnx;
285311
if (TMath::Abs(val) < 1) {
286-
registry.fill(tarName, fPtAxis->GetBinCenter(ptbin), fMass->GetBinCenter(massbin), cent, val, dnx);
312+
registry.fill(tarName, fpt->GetBinCenter(ptbin), fMass->GetBinCenter(massbin), cent, val, dnx);
287313
}
288314
}
289315
return;
@@ -409,7 +435,7 @@ struct FlowGFWOmegaXi {
409435
wacc = 1 / vecwa[phibin];
410436
registry.fill(HIST("hPhi"), track.phi());
411437
registry.fill(HIST("hEta"), track.eta());
412-
registry.fill(HIST("hEtaPhiREF"), track.eta(), track.phi());
438+
registry.fill(HIST("hEtaPhiVtxzREF"), track.eta(), track.phi(), vtxz);
413439
registry.fill(HIST("hPt"), track.pt());
414440
int ptbin = fPtAxis->FindBin(track.pt()) - 1;
415441
if ((track.pt() > cfgCutPtMin) && (track.pt() < cfgCutPtMax)) {
@@ -421,13 +447,13 @@ struct FlowGFWOmegaXi {
421447
auto v0posdau = v0.posTrack_as<DaughterTracks>();
422448
auto v0negdau = v0.negTrack_as<DaughterTracks>();
423449
// check tpc
424-
int partical = 0;
450+
int PDGCode = 0;
425451
if (v0.qtarm() / TMath::Abs(v0.alpha()) > cfgv0_ArmPodocut && TMath::Abs(v0posdau.tpcNSigmaPi()) < cfgNSigmaCascPion && TMath::Abs(v0negdau.tpcNSigmaPi()) < cfgNSigmaCascPion) {
426452
registry.fill(HIST("InvMassK0s_all"), v0.pt(), v0.mK0Short(), v0.eta(), cent);
427-
partical = 310;
453+
PDGCode = kK0Short;
428454
} else if (v0.qtarm() / TMath::Abs(v0.alpha()) < cfgv0_ArmPodocut && TMath::Abs(v0posdau.tpcNSigmaPr()) < cfgNSigmaCascProton && TMath::Abs(v0negdau.tpcNSigmaPi()) < cfgNSigmaCascPion) {
429455
registry.fill(HIST("InvMassLambda_all"), v0.pt(), v0.mLambda(), v0.eta(), cent);
430-
partical = 3122;
456+
PDGCode = kLambda0;
431457
}
432458
// track quality check
433459
if (v0posdau.tpcNClsFound() < cfgtpcclusters)
@@ -445,17 +471,17 @@ struct FlowGFWOmegaXi {
445471
continue;
446472
if (v0.dcaV0daughters() > cfgv0_dcav0dau)
447473
continue;
448-
if (partical == 310) {
474+
if (PDGCode == kK0Short) {
449475
if (TMath::Abs(v0.mK0Short() - 0.49761) < cfgv0_mk0swindow) {
450476
registry.fill(HIST("InvMassK0s"), v0.pt(), v0.mK0Short(), v0.eta(), cent);
451-
registry.fill(HIST("hEtaPhiPOIK0s"), v0.eta(), v0.phi());
452-
fGFW->Fill(v0.eta(), fPtAxis->FindBin(v0.pt()) - 1 + ((fK0sMass->FindBin(v0.mK0Short()) - 1) * nPtBins), v0.phi(), wacc * weff, 8);
477+
registry.fill(HIST("hEtaPhiVtxzPOIK0s"), v0.eta(), v0.phi(), vtxz);
478+
fGFW->Fill(v0.eta(), fV0PtAxis->FindBin(v0.pt()) - 1 + ((fK0sMass->FindBin(v0.mK0Short()) - 1) * nV0PtBins), v0.phi(), wacc * weff, 8);
453479
}
454-
} else if (partical == 3122) {
480+
} else if (PDGCode == kLambda0) {
455481
if (TMath::Abs(v0.mLambda() - 1.115683) < cfgv0_mlambdawindow) {
456482
registry.fill(HIST("InvMassLambda"), v0.pt(), v0.mLambda(), v0.eta(), cent);
457-
registry.fill(HIST("hEtaPhiPOILambda"), v0.eta(), v0.phi());
458-
fGFW->Fill(v0.eta(), fPtAxis->FindBin(v0.pt()) - 1 + ((fLambdaMass->FindBin(v0.mLambda()) - 1) * nPtBins), v0.phi(), wacc * weff, 16);
483+
registry.fill(HIST("hEtaPhiVtxzPOILambda"), v0.eta(), v0.phi(), vtxz);
484+
fGFW->Fill(v0.eta(), fV0PtAxis->FindBin(v0.pt()) - 1 + ((fLambdaMass->FindBin(v0.mLambda()) - 1) * nV0PtBins), v0.phi(), wacc * weff, 16);
459485
}
460486
}
461487
}
@@ -468,13 +494,13 @@ struct FlowGFWOmegaXi {
468494
if (cfgcheckDauTPC && (!posdau.hasTPC() || !negdau.hasTPC() || !bachelor.hasTPC())) {
469495
continue;
470496
}
471-
int partical = 0;
497+
int PDGCode = 0;
472498
if (casc.sign() < 0 && TMath::Abs(casc.yOmega()) < cfgCasc_rapidity && TMath::Abs(bachelor.tpcNSigmaKa()) < cfgNSigmaCascKaon && TMath::Abs(posdau.tpcNSigmaPr()) < cfgNSigmaCascProton && TMath::Abs(negdau.tpcNSigmaPi()) < cfgNSigmaCascPion) {
473499
registry.fill(HIST("InvMassOmegaMinus_all"), casc.pt(), casc.mOmega(), casc.eta(), cent);
474-
partical = 3334;
500+
PDGCode = kOmegaMinus;
475501
} else if (casc.sign() < 0 && TMath::Abs(casc.yXi()) < cfgCasc_rapidity && TMath::Abs(bachelor.tpcNSigmaPi()) < cfgNSigmaCascPion && TMath::Abs(posdau.tpcNSigmaPr()) < cfgNSigmaCascProton && TMath::Abs(negdau.tpcNSigmaPi()) < cfgNSigmaCascPion) {
476502
registry.fill(HIST("InvMassXiMinus_all"), casc.pt(), casc.mXi(), casc.eta(), cent);
477-
partical = 3312;
503+
PDGCode = kXiMinus;
478504
} else {
479505
continue;
480506
}
@@ -508,30 +534,32 @@ struct FlowGFWOmegaXi {
508534
continue;
509535
if (negdau.itsNCls() < cfgitsclusters)
510536
continue;
511-
if (partical == 3334) {
512-
registry.fill(HIST("hEtaPhiPOIOmega"), casc.eta(), casc.phi());
537+
if (PDGCode == kOmegaMinus) {
538+
registry.fill(HIST("hEtaPhiVtxzPOIOmega"), casc.eta(), casc.phi(), vtxz);
513539
registry.fill(HIST("InvMassOmegaMinus"), casc.pt(), casc.mOmega(), casc.eta(), cent);
514540
if ((casc.pt() < cfgCutPtPOIMax) && (casc.pt() > cfgCutPtPOIMin) && (casc.mOmega() > 1.63) && (casc.mOmega() < 1.71)) {
515-
fGFW->Fill(casc.eta(), fPtAxis->FindBin(casc.pt()) - 1 + ((fOmegaMass->FindBin(casc.mOmega()) - 1) * nPtBins), casc.phi(), wacc * weff, 4);
541+
fGFW->Fill(casc.eta(), fXiPtAxis->FindBin(casc.pt()) - 1 + ((fOmegaMass->FindBin(casc.mOmega()) - 1) * nXiPtBins), casc.phi(), wacc * weff, 4);
516542
}
517543
}
518-
if (partical == 3312) {
519-
registry.fill(HIST("hEtaPhiPOIXi"), casc.eta(), casc.phi());
544+
if (PDGCode == kXiMinus) {
545+
registry.fill(HIST("hEtaPhiVtxzPOIXi"), casc.eta(), casc.phi(), vtxz);
520546
registry.fill(HIST("InvMassXiMinus"), casc.pt(), casc.mXi(), casc.eta(), cent);
521547
if ((casc.pt() < cfgCutPtPOIMax) && (casc.pt() > cfgCutPtPOIMin) && (casc.mXi() > 1.30) && (casc.mXi() < 1.37)) {
522-
fGFW->Fill(casc.eta(), fPtAxis->FindBin(casc.pt()) - 1 + ((fXiMass->FindBin(casc.mXi()) - 1) * nPtBins), casc.phi(), wacc * weff, 2);
548+
fGFW->Fill(casc.eta(), fXiPtAxis->FindBin(casc.pt()) - 1 + ((fXiMass->FindBin(casc.mXi()) - 1) * nXiPtBins), casc.phi(), wacc * weff, 2);
523549
}
524550
}
525551
}
526552
// Filling cumulant with ROOT TProfile
527553
FillProfile(corrconfigs.at(0), HIST("c22"), cent);
528554
FillProfile(corrconfigs.at(1), HIST("c24"), cent);
529-
for (int i = 1; i <= nPtBins; i++) // loop for all ptBins
555+
for (int i = 1; i <= nV0PtBins; i++) // loop for all ptBins
530556
{
531-
FillProfilepT(corrconfigs.at(2), HIST("Xic22dpt"), i, 3312, cent);
532-
FillProfilepT(corrconfigs.at(3), HIST("Omegac22dpt"), i, 3334, cent);
533-
FillProfilepT(corrconfigs.at(4), HIST("K0sc22dpt"), i, 310, cent);
534-
FillProfilepT(corrconfigs.at(5), HIST("Lambdac22dpt"), i, 3122, cent);
557+
FillProfilepT(corrconfigs.at(4), HIST("K0sc22dpt"), i, kK0Short, cent);
558+
FillProfilepT(corrconfigs.at(5), HIST("Lambdac22dpt"), i, kLambda0, cent);
559+
}
560+
for (int i = 1; i <= nXiPtBins; i++) {
561+
FillProfilepT(corrconfigs.at(2), HIST("Xic22dpt"), i, kXiMinus, cent);
562+
FillProfilepT(corrconfigs.at(3), HIST("Omegac22dpt"), i, kOmegaMinus, cent);
535563
}
536564
}
537565
};

0 commit comments

Comments
 (0)