Skip to content

Commit 656dd81

Browse files
FDUEnrichFDUEnrich
committed
Modified according to Mattia's suggestions
1 parent fa7cc4e commit 656dd81

File tree

2 files changed

+474
-469
lines changed

2 files changed

+474
-469
lines changed

‎PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx‎

Lines changed: 66 additions & 65 deletions
Original file line numberDiff line numberDiff line change
@@ -138,10 +138,10 @@ struct HfCandidateCreatorXic0Omegac0 {
138138
struct {
139139
float chi2GeoV0;
140140
float ldlV0;
141-
float chi2TopoV0ToPv;
141+
float chi2NdfTopoV0ToPv;
142142
float chi2GeoCasc;
143143
float ldlCasc;
144-
float chi2TopoCascToPv;
144+
float chi2NdfTopoCascToPv;
145145
float decayLenXYLambda;
146146
float decayLenXYCasc;
147147
float cosPaV0ToCasc;
@@ -159,17 +159,17 @@ struct HfCandidateCreatorXic0Omegac0 {
159159
float rapOmegac;
160160
float massOmegac;
161161
float cosThetaStarPiFromOmegac;
162-
float chi2TopoPiFromOmegacToPv;
162+
float chi2NdfTopoPiFromOmegacToPv;
163163
float kfDcaXYPiFromOmegac;
164-
float chi2TopoV0ToCasc;
165-
float chi2TopoCascToOmegac;
164+
float chi2NdfTopoV0ToCasc;
165+
float chi2NdfTopoCascToOmegac;
166166
float decayLenXYOmegac;
167167
float chi2GeoOmegac;
168168
float kfDcaV0Dau;
169169
float kfDcaCascDau;
170170
float kfDcaOmegacDau;
171171
float kfDcaXYCascToPv;
172-
float chi2TopoOmegacToPv;
172+
float chi2NdfTopoOmegacToPv;
173173
float cosPaOmegacToPv;
174174
float cosPaXYOmegacToPv;
175175
float ldlOmegac;
@@ -278,7 +278,7 @@ struct HfCandidateCreatorXic0Omegac0 {
278278
hCascadesCounterToOmegaPi = registry.add<TH1>("hCascadesCounterToOmegaPi", "Cascades counter wrt derived data - #Omega #pi decay;status;entries", {HistType::kTH1D, {{2, -0.5, 1.5}}}); // 0 --> cascades in derived data table (and stored in AOD table), 1 --> cascades in derived data table and also accessible in cascData table
279279
hCascadesCounterToOmegaK = registry.add<TH1>("hCascadesCounterToOmegaK", "Cascades counter wrt derived data - #Omega K decay;status;entries", {HistType::kTH1D, {{2, -0.5, 1.5}}}); // 0 --> cascades in derived data table (and stored in AOD table), 1 --> cascades in derived data table and also accessible in cascData table
280280

281-
// KFparticle variables hist
281+
// KFParticle Variables Histograms
282282
registry.add("hKFParticleV0TopoChi2", "hKFParticleV0TopoChi2", kTH1D, {{1000, -0.10f, 100.0f}});
283283
registry.add("hKFParticleCascTopoChi2", "hKFParticleCascTopoChi2", kTH1D, {{1000, -0.1f, 100.0f}});
284284
registry.add("hKFParticleCascBachTopoChi2", "hKFParticleCascBachTopoChi2", kTH1D, {{1000, -0.1f, 100.0f}});
@@ -293,37 +293,39 @@ struct HfCandidateCreatorXic0Omegac0 {
293293
registry.add("hInvMassOmegaMinus", "hInvMassOmegaMinus", kTH1D, {{1000, 1.6f, 2.0f}});
294294
registry.add("hInvMassXiMinus", "hInvMassXiMinus", kTH1D, {{1000, 1.25f, 1.65f}});
295295
registry.add("hInvMassXiMinus_rej", "hInvMassXiMinus_rej", kTH1D, {{1000, 1.25f, 1.65f}});
296-
// Additional KF hist
297-
registry.add("hKFParticlechi2TopoOmegacToPv", "hKFParticlechi2TopoOmegacToPv", kTH1D, {{1000, -0.1f, 100.0f}}); // new
298-
registry.add("hKFParticlechi2TopoCascToPv", "hKFParticlechi2TopoCascToPv", kTH1D, {{1000, -0.1f, 100.0f}}); // new
299-
registry.add("hKFParticleDcaXYV0DauPosToPv", "hKFParticleDcaXYV0DauPosToPv", kTH1D, {{1000, -0.1f, 30.0f}}); // new
300-
registry.add("hKFParticleDcaXYV0DauNegToPv", "hKFParticleDcaXYV0DauNegToPv", kTH1D, {{1000, -0.1f, 30.0f}}); // new
301-
registry.add("hEtaV0PosDau", "hEtaV0PosDau", kTH1D, {{1000, -5.0f, 5.0f}}); // new
302-
registry.add("hEtaV0NegDau", "hEtaV0NegDau", kTH1D, {{1000, -5.0f, 5.0f}}); // new
303-
registry.add("hEtaKaFromCasc", "hEtaKaFromCasc", kTH1D, {{1000, -5.0f, 5.0f}}); // new
304-
registry.add("hEtaPiFromCharmBaryon", "hEtaPiFromCharmBaryon", kTH1D, {{1000, -5.0f, 5.0f}}); // new
305-
registry.add("hCascradius", "hCascradius", kTH1D, {{1000, 0.0f, 50.0f}}); // new
306-
registry.add("hV0radius", "hV0radius", kTH1D, {{1000, 0.0f, 50.0f}}); // new
307-
registry.add("hCosPACasc", "hCosPACasc", kTH1D, {{5000, 0.8f, 1.0f}}); // new
308-
registry.add("hCosPAV0", "hCosPAV0", kTH1D, {{5000, 0.8f, 1.0f}}); // new
309-
registry.add("hDcaCascDau", "hDcaCascDau", kTH1D, {{1000, -0.1f, 10.0f}}); // new
310-
registry.add("hDcaV0Dau", "hDcaV0Dau", kTH1D, {{1000, -0.1f, 10.0f}}); // new
311-
registry.add("hDcaXYToPvKa", "hDcaXYToPvKa", kTH1D, {{1000, -0.1f, 10.0f}}); // new
312-
registry.add("hImpactParBachFromCharmBaryonXY", "hImpactParBachFromCharmBaryonXY", kTH1D, {{1000, -1.0f, 1.0f}}); // new
313-
registry.add("hImpactParBachFromCharmBaryonZ", "hImpactParBachFromCharmBaryonZ", kTH1D, {{1000, -1.0f, 1.0f}}); // new
314-
registry.add("hImpactParCascXY", "hImpactParCascXY", kTH1D, {{1000, -4.0f, 4.0f}}); // new
315-
registry.add("hImpactParCascZ", "hImpactParCascZ", kTH1D, {{1000, -5.0f, 5.0f}}); // new
316-
registry.add("hPtKaFromCasc", "hPtKaFromCasc", kTH1D, {{1000, 0.0f, 5.0f}}); // new
317-
registry.add("hPtPiFromCharmBaryon", "hPtPiFromCharmBaryon", kTH1D, {{1000, 0.0f, 5.0f}}); // new
318-
registry.add("hCTauOmegac", "hCTauOmegac", kTH1D, {{1000, 0.0f, 0.1f}}); // new
319-
registry.add("hKFGeoV0Chi2OverNdf", "hKFGeoV0Chi2OverNdf", kTH1D, {{1000, 0.0f, 100.0f}}); // new
320-
registry.add("hKFGeoCascChi2OverNdf", "hKFGeoCascChi2OverNdf", kTH1D, {{1000, 0.0f, 100.0f}}); // new
321-
registry.add("hKFGeoCharmbaryonChi2OverNdf", "hKFGeoCharmbaryonChi2OverNdf", kTH1D, {{1000, 0.0f, 100.0f}}); // new
322-
registry.add("hKFdecayLenXYLambda", "hKFdecayLenXYLambda", kTH1D, {{1000, 0.0f, 50.0f}}); // new
323-
registry.add("hKFdecayLenXYCasc", "hKFdecayLenXYCasc", kTH1D, {{1000, 0.0f, 50.0f}}); // new
324-
registry.add("hKFdecayLenXYOmegac", "hKFdecayLenXYOmegac", kTH1D, {{1000, 0.0f, 50.0f}}); // new
325-
registry.add("hKFcosPaV0ToCasc", "hKFcosPaV0ToCasc", kTH1D, {{5000, 0.8f, 1.0f}}); // new
326-
registry.add("hKFcosPaCascToOmegac", "hKFcosPaCascToOmegac", kTH1D, {{5000, 0.8f, 1.0f}}); // new
296+
297+
// Additional KFParticle Histograms
298+
registry.add("hKFParticlechi2TopoOmegacToPv", "hKFParticlechi2TopoOmegacToPv", kTH1D, {{1000, -0.1f, 100.0f}});
299+
registry.add("hKFParticlechi2TopoCascToPv", "hKFParticlechi2TopoCascToPv", kTH1D, {{1000, -0.1f, 100.0f}});
300+
registry.add("hKFParticleDcaXYV0DauPosToPv", "hKFParticleDcaXYV0DauPosToPv", kTH1D, {{1000, -0.1f, 30.0f}});
301+
registry.add("hKFParticleDcaXYV0DauNegToPv", "hKFParticleDcaXYV0DauNegToPv", kTH1D, {{1000, -0.1f, 30.0f}});
302+
registry.add("hEtaV0PosDau", "hEtaV0PosDau", kTH1D, {{1000, -5.0f, 5.0f}});
303+
registry.add("hEtaV0NegDau", "hEtaV0NegDau", kTH1D, {{1000, -5.0f, 5.0f}});
304+
registry.add("hEtaKaFromCasc", "hEtaKaFromCasc", kTH1D, {{1000, -5.0f, 5.0f}});
305+
registry.add("hEtaPiFromCharmBaryon", "hEtaPiFromCharmBaryon", kTH1D, {{1000, -5.0f, 5.0f}});
306+
registry.add("hCascradius", "hCascradius", kTH1D, {{1000, 0.0f, 50.0f}});
307+
registry.add("hV0radius", "hV0radius", kTH1D, {{1000, 0.0f, 50.0f}});
308+
registry.add("hCosPACasc", "hCosPACasc", kTH1D, {{5000, 0.8f, 1.1f}});
309+
registry.add("hCosPAV0", "hCosPAV0", kTH1D, {{5000, 0.8f, 1.1f}});
310+
registry.add("hDcaCascDau", "hDcaCascDau", kTH1D, {{1000, -0.1f, 10.0f}});
311+
registry.add("hDcaV0Dau", "hDcaV0Dau", kTH1D, {{1000, -0.1f, 10.0f}});
312+
registry.add("hDcaXYToPvKa", "hDcaXYToPvKa", kTH1D, {{1000, -0.1f, 10.0f}});
313+
registry.add("hImpactParBachFromCharmBaryonXY", "hImpactParBachFromCharmBaryonXY", kTH1D, {{1000, -1.0f, 1.0f}});
314+
registry.add("hImpactParBachFromCharmBaryonZ", "hImpactParBachFromCharmBaryonZ", kTH1D, {{1000, -2.0f, 2.0f}});
315+
registry.add("hImpactParCascXY", "hImpactParCascXY", kTH1D, {{1000, -4.0f, 4.0f}});
316+
registry.add("hImpactParCascZ", "hImpactParCascZ", kTH1D, {{1000, -5.0f, 5.0f}});
317+
registry.add("hPtKaFromCasc", "hPtKaFromCasc", kTH1D, {{1000, 0.0f, 5.0f}});
318+
registry.add("hPtPiFromCharmBaryon", "hPtPiFromCharmBaryon", kTH1D, {{1000, 0.0f, 5.0f}});
319+
registry.add("hCTauOmegac", "hCTauOmegac", kTH1D, {{1000, 0.0f, 0.1f}});
320+
registry.add("hKFGeoV0Chi2OverNdf", "hKFGeoV0Chi2OverNdf", kTH1D, {{1000, 0.0f, 100.0f}});
321+
registry.add("hKFGeoCascChi2OverNdf", "hKFGeoCascChi2OverNdf", kTH1D, {{1000, 0.0f, 100.0f}});
322+
registry.add("hKFGeoCharmbaryonChi2OverNdf", "hKFGeoCharmbaryonChi2OverNdf", kTH1D, {{1000, 0.0f, 100.0f}});
323+
registry.add("hKFdecayLenXYLambda", "hKFdecayLenXYLambda", kTH1D, {{1000, 0.0f, 50.0f}});
324+
registry.add("hKFdecayLenXYCasc", "hKFdecayLenXYCasc", kTH1D, {{1000, 0.0f, 50.0f}});
325+
registry.add("hKFdecayLenXYOmegac", "hKFdecayLenXYOmegac", kTH1D, {{1000, 0.0f, 50.0f}});
326+
registry.add("hKFcosPaV0ToCasc", "hKFcosPaV0ToCasc", kTH1D, {{5000, 0.8f, 1.1f}});
327+
registry.add("hKFcosPaCascToOmegac", "hKFcosPaCascToOmegac", kTH1D, {{5000, 0.8f, 1.1f}});
328+
327329

328330
hfEvSel.addHistograms(registry); // collision monitoring
329331

@@ -758,26 +760,26 @@ struct HfCandidateCreatorXic0Omegac0 {
758760
KFParticle kfPosPr(kfTrack0, kProton);
759761
KFParticle kfNegPi(kfTrack1, kPiMinus);
760762
KFParticle kfNegKa(kfTrackBach, kKMinus);
761-
KFParticle kfNegPirej(kfTrackBach, kPiMinus); // rej
763+
KFParticle kfNegPirRej(kfTrackBach, kPiMinus); // rej
762764
KFParticle kfPosPi(kfTrack0, kPiPlus);
763765
KFParticle kfNegPr(kfTrack1, kProton);
764766
KFParticle kfPosKa(kfTrackBach, kKPlus);
765-
KFParticle kfPosPirej(kfTrackBach, kPiPlus); // rej
767+
KFParticle kfPosPiRej(kfTrackBach, kPiPlus); // rej
766768

767769
KFParticle kfBachKaon;
768770
KFParticle kfPos;
769771
KFParticle kfNeg;
770-
KFParticle kfBachrej; // rej
772+
KFParticle kfBachPionRej; // rej
771773
if (bachCharge < 0) {
772774
kfPos = kfPosPr;
773775
kfNeg = kfNegPi;
774776
kfBachKaon = kfNegKa;
775-
kfBachrej = kfNegPirej; // rej
777+
kfBachPionRej = kfNegPirRej; // rej
776778
} else {
777779
kfPos = kfPosPi;
778780
kfNeg = kfNegPr;
779781
kfBachKaon = kfPosKa;
780-
kfBachrej = kfPosPirej; // rej
782+
kfBachPionRej = kfPosPiRej; // rej
781783
}
782784

783785
//__________________________________________
@@ -815,15 +817,15 @@ struct HfCandidateCreatorXic0Omegac0 {
815817
//__________________________________________
816818
//*>~<* step 2 : reconstruct cascade(Omega) with KF
817819
const KFParticle* omegaDaugthers[2] = {&kfBachKaon, &kfV0};
818-
const KFParticle* omegaDaugthersRej[2] = {&kfBachRej, &kfV0}; // rej
820+
const KFParticle* omegaDaugthersRej[2] = {&kfBachPionRej, &kfV0}; // rej
819821
// construct cascade
820822
KFParticle kfOmega;
821823
KFParticle kfOmegarej; // rej
822824
kfOmega.SetConstructMethod(kfConstructMethod);
823825
kfOmegarej.SetConstructMethod(kfConstructMethod); // rej
824826
try {
825827
kfOmega.Construct(omegaDaugthers, 2);
826-
kfOmegarej.Construct(omegaDaugthersrej, 2); // rej
828+
kfOmegarej.Construct(omegaDaugthersRej, 2); // rej
827829
} catch (std::runtime_error& e) {
828830
LOG(debug) << "Failed to construct Omega or Omega_rej from V0 and bachelor track: " << e.what();
829831
continue;
@@ -851,6 +853,7 @@ struct HfCandidateCreatorXic0Omegac0 {
851853
}
852854
registry.fill(HIST("hInvMassOmegaMinus"), massCasc);
853855
kfOmega.TransportToDecayVertex();
856+
//rej: Add competing rejection to minimize misidentified Xi impact. Reject if kfBachPionRej is Pion and the constructed cascade has Xi's invariant mass.
854857

855858
//__________________________________________
856859
//*>~<* step 3 : reconstruc Omegac0 with KF
@@ -1013,14 +1016,14 @@ struct HfCandidateCreatorXic0Omegac0 {
10131016
auto cascChi2OverNdfm = kfOmegac0Candidate.chi2MassCasc / cascNdfm;
10141017

10151018
// KF topo Chi2 over NDF
1016-
kfOmegac0Candidate.chi2TopoV0ToPv = kfV0ToPv.GetChi2() / kfV0ToPv.GetNDF();
1017-
kfOmegac0Candidate.chi2TopoCascToPv = kfOmegaToPv.GetChi2() / kfOmegaToPv.GetNDF();
1018-
kfOmegac0Candidate.chi2TopoPiFromOmegacToPv = kfPiFromOmegacToPv.GetChi2() / kfPiFromOmegacToPv.GetNDF();
1019-
kfOmegac0Candidate.chi2TopoOmegacToPv = kfOmegac0ToPv.GetChi2() / kfOmegac0ToPv.GetNDF();
1019+
kfOmegac0Candidate.chi2NdfTopoV0ToPv = kfV0ToPv.GetChi2() / kfV0ToPv.GetNDF();
1020+
kfOmegac0Candidate.chi2NdfTopoCascToPv = kfOmegaToPv.GetChi2() / kfOmegaToPv.GetNDF();
1021+
kfOmegac0Candidate.chi2NdfTopoPiFromOmegacToPv = kfPiFromOmegacToPv.GetChi2() / kfPiFromOmegacToPv.GetNDF();
1022+
kfOmegac0Candidate.chi2NdfTopoOmegacToPv = kfOmegac0ToPv.GetChi2() / kfOmegac0ToPv.GetNDF();
10201023

1021-
auto cascBachTopoChi2 = kfBachKaonToOmega.GetChi2() / kfBachKaonToOmega.GetNDF();
1022-
kfOmegac0Candidate.chi2TopoV0ToCasc = kfV0ToCasc.GetChi2() / kfV0ToCasc.GetNDF();
1023-
kfOmegac0Candidate.chi2TopoCascToOmegac = kfOmegaToOmegaC.GetChi2() / kfOmegaToOmegaC.GetNDF();
1024+
auto cascBachTopoChi2Ndf = kfBachKaonToOmega.GetChi2() / kfBachKaonToOmega.GetNDF();
1025+
kfOmegac0Candidate.chi2NdfTopoV0ToCasc = kfV0ToCasc.GetChi2() / kfV0ToCasc.GetNDF();
1026+
kfOmegac0Candidate.chi2NdfTopoCascToOmegac = kfOmegaToOmegaC.GetChi2() / kfOmegaToOmegaC.GetNDF();
10241027

10251028
// KF ldl
10261029
kfOmegac0Candidate.ldlV0 = ldlFromKF(kfV0, kfPV);
@@ -1083,15 +1086,15 @@ struct HfCandidateCreatorXic0Omegac0 {
10831086
kfOmegac0Candidate.etaOmegac = kfOmegaC0.GetEta();
10841087

10851088
// fill KF hist
1086-
registry.fill(HIST("hKFParticleCascBachTopoChi2"), cascBachTopoChi2);
1087-
registry.fill(HIST("hKFParticleV0TopoChi2"), kfOmegac0Candidate.chi2TopoV0ToCasc);
1088-
registry.fill(HIST("hKFParticleCascTopoChi2"), kfOmegac0Candidate.chi2TopoCascToOmegac);
1089-
registry.fill(HIST("hKFParticlechi2TopoOmegacToPv"), kfOmegac0Candidate.chi2TopoOmegacToPv); // new
1090-
registry.fill(HIST("hKFParticlechi2TopoCascToPv"), kfOmegac0Candidate.chi2TopoCascToPv); // new
1089+
registry.fill(HIST("hKFParticleCascBachTopoChi2"), cascBachTopoChi2Ndf);
1090+
registry.fill(HIST("hKFParticleV0TopoChi2"), kfOmegac0Candidate.chi2NdfTopoV0ToCasc);
1091+
registry.fill(HIST("hKFParticleCascTopoChi2"), kfOmegac0Candidate.chi2NdfTopoCascToOmegac);
1092+
registry.fill(HIST("hKFParticlechi2TopoOmegacToPv"), kfOmegac0Candidate.chi2NdfTopoOmegacToPv);
1093+
registry.fill(HIST("hKFParticlechi2TopoCascToPv"), kfOmegac0Candidate.chi2NdfTopoCascToPv);
10911094
registry.fill(HIST("hKFParticleDcaCharmBaryonDau"), kfOmegac0Candidate.kfDcaOmegacDau);
10921095
registry.fill(HIST("hKFParticleDcaXYCascBachToPv"), dcaxyCascBachelor);
1093-
registry.fill(HIST("hKFParticleDcaXYV0DauPosToPv"), dcaxyV0Dau0); // new
1094-
registry.fill(HIST("hKFParticleDcaXYV0DauNegToPv"), dcaxyV0Dau1); // new
1096+
registry.fill(HIST("hKFParticleDcaXYV0DauPosToPv"), dcaxyV0Dau0);
1097+
registry.fill(HIST("hKFParticleDcaXYV0DauNegToPv"), dcaxyV0Dau1);
10951098
registry.fill(HIST("hKfLambda_ldl"), kfOmegac0Candidate.ldlV0);
10961099
registry.fill(HIST("hKfOmega_ldl"), kfOmegac0Candidate.ldlCasc);
10971100
registry.fill(HIST("hKfOmegaC0_ldl"), kfOmegac0Candidate.ldlOmegac);
@@ -1159,8 +1162,8 @@ struct HfCandidateCreatorXic0Omegac0 {
11591162
kfCandidateData(kfOmegac0Candidate.kfDcaXYPiFromOmegac, kfOmegac0Candidate.kfDcaXYCascToPv,
11601163
kfOmegac0Candidate.chi2GeoV0, kfOmegac0Candidate.chi2GeoCasc, kfOmegac0Candidate.chi2GeoOmegac, kfOmegac0Candidate.chi2MassV0, kfOmegac0Candidate.chi2MassCasc,
11611164
kfOmegac0Candidate.ldlV0, kfOmegac0Candidate.ldlCasc, kfOmegac0Candidate.ldlOmegac,
1162-
kfOmegac0Candidate.chi2TopoV0ToPv, kfOmegac0Candidate.chi2TopoCascToPv, kfOmegac0Candidate.chi2TopoPiFromOmegacToPv, kfOmegac0Candidate.chi2TopoOmegacToPv,
1163-
kfOmegac0Candidate.chi2TopoV0ToCasc, kfOmegac0Candidate.chi2TopoCascToOmegac,
1165+
kfOmegac0Candidate.chi2NdfTopoV0ToPv, kfOmegac0Candidate.chi2NdfTopoCascToPv, kfOmegac0Candidate.chi2NdfTopoPiFromOmegacToPv, kfOmegac0Candidate.chi2NdfTopoOmegacToPv,
1166+
kfOmegac0Candidate.chi2NdfTopoV0ToCasc, kfOmegac0Candidate.chi2NdfTopoCascToOmegac,
11641167
kfOmegac0Candidate.decayLenXYLambda, kfOmegac0Candidate.decayLenXYCasc, kfOmegac0Candidate.decayLenXYOmegac,
11651168
kfOmegac0Candidate.cosPaV0ToCasc, kfOmegac0Candidate.cosPaCascToOmegac, kfOmegac0Candidate.cosPaXYV0ToCasc, kfOmegac0Candidate.cosPaXYCascToOmegac,
11661169
kfOmegac0Candidate.rapOmegac, kfOmegac0Candidate.ptPiFromOmegac, kfOmegac0Candidate.ptOmegac,
@@ -2118,19 +2121,17 @@ struct HfCandidateCreatorXic0Omegac0Mc {
21182121
// Slice the collisions table to get the collision info for the current MC collision
21192122
float centrality{-1.f};
21202123
uint16_t rejectionMask{0};
2121-
int nSplitColl = 0;
21222124
if constexpr (centEstimator == CentralityEstimator::FT0C) {
21232125
const auto collSlice = collsWithMcLabels.sliceBy(colPerMcCollisionFT0C, mcCollision.globalIndex());
21242126
rejectionMask = hfEvSelMc.getHfMcCollisionRejectionMask<BCsInfo, centEstimator>(mcCollision, collSlice, centrality);
21252127
} else if constexpr (centEstimator == CentralityEstimator::FT0M) {
21262128
const auto collSlice = collsWithMcLabels.sliceBy(colPerMcCollisionFT0M, mcCollision.globalIndex());
2127-
nSplitColl = collSlice.size();
21282129
rejectionMask = hfEvSelMc.getHfMcCollisionRejectionMask<BCsInfo, centEstimator>(mcCollision, collSlice, centrality);
21292130
} else if constexpr (centEstimator == CentralityEstimator::None) {
21302131
const auto collSlice = collsWithMcLabels.sliceBy(colPerMcCollision, mcCollision.globalIndex());
21312132
rejectionMask = hfEvSelMc.getHfMcCollisionRejectionMask<BCsInfo, centEstimator>(mcCollision, collSlice, centrality);
21322133
}
2133-
hfEvSelMc.fillHistograms<centEstimator>(mcCollision, rejectionMask, nSplitColl);
2134+
hfEvSelMc.fillHistograms<centEstimator>(mcCollision, rejectionMask);
21342135
if (rejectionMask != 0) {
21352136
/// at least one event selection not satisfied --> reject all particles from this collision
21362137
for (unsigned int i = 0; i < mcParticlesPerMcColl.size(); ++i) {
@@ -2501,4 +2502,4 @@ WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)
25012502
return WorkflowSpec{
25022503
adaptAnalysisTask<HfCandidateCreatorXic0Omegac0>(cfgc),
25032504
adaptAnalysisTask<HfCandidateCreatorXic0Omegac0Mc>(cfgc)};
2504-
}
2505+
}

0 commit comments

Comments
 (0)