Skip to content

Commit a8c7c4f

Browse files
axti98axti98
authored
Merge branch 'master' into Alice3MultDev
2 parents 9bff573 + 814256a commit a8c7c4f

61 files changed

Lines changed: 6054 additions & 1855 deletions

File tree

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

‎ALICE3/DataModel/OTFStrangeness.h‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -314,7 +314,7 @@ DECLARE_SOA_COLUMN(DcaV0Daughters, dcaV0Daughters, float);
314314
DECLARE_SOA_COLUMN(V0Radius, v0Radius, float);
315315
DECLARE_SOA_COLUMN(MLambda, mLambda, float);
316316
DECLARE_SOA_COLUMN(MAntiLambda, mAntiLambda, float);
317-
DECLARE_SOA_COLUMN(MK0, mK0, float);
317+
DECLARE_SOA_COLUMN(MK0Short, mK0Short, float);
318318

319319
// kinematics
320320
DECLARE_SOA_COLUMN(Pt, pt, float);
@@ -330,7 +330,7 @@ DECLARE_SOA_TABLE(UpgradeV0s, "AOD", "UPGRADEV0S",
330330
otfv0::V0Radius,
331331
otfv0::MLambda,
332332
otfv0::MAntiLambda,
333-
otfv0::MK0,
333+
otfv0::MK0Short,
334334
otfv0::Pt);
335335

336336
using UpgradeV0 = UpgradeV0s::iterator;

‎ALICE3/TableProducer/OTF/onTheFlyTracker.cxx‎

Lines changed: 284 additions & 268 deletions
Large diffs are not rendered by default.

‎ALICE3/TableProducer/alice3TrackingTranslator.cxx‎

Lines changed: 473 additions & 235 deletions
Large diffs are not rendered by default.

‎ALICE3/TableProducer/alice3strangenessFinder.cxx‎

Lines changed: 26 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -309,6 +309,17 @@ struct Alice3strangenessFinder {
309309
histos.add("CascadeBuilding/hXiMass", "", kTH1D, {axisXiMass});
310310
histos.add("CascadeBuilding/hOmegaMass", "", kTH1D, {axisOmegaMass});
311311
}
312+
histos.add("hPx", "", kTH1D, {{200, -10.f, 10.f}});
313+
histos.add("hPy", "", kTH1D, {{200, -10.f, 10.f}});
314+
histos.add("hPz", "", kTH1D, {{200, -10.f, 10.f}});
315+
316+
histos.add("hPxPairs", "", kTH1D, {{200, -10.f, 10.f}});
317+
histos.add("hPyPairs", "", kTH1D, {{200, -10.f, 10.f}});
318+
histos.add("hPzPairs", "", kTH1D, {{200, -10.f, 10.f}});
319+
320+
histos.add("hPxV0CandDaughter", "", kTH1D, {{200, -10.f, 10.f}});
321+
histos.add("hPyV0CandDaughter", "", kTH1D, {{200, -10.f, 10.f}});
322+
histos.add("hPzV0CandDaughter", "", kTH1D, {{200, -10.f, 10.f}});
312323

313324
histos.print();
314325
}
@@ -516,6 +527,14 @@ struct Alice3strangenessFinder {
516527
}
517528

518529
o2::track::TrackParCov pos = getTrackParCov(posTrack);
530+
// if constexpr (requires { posTrack.pdgCode(); }) {
531+
// if (posTrack.pdgCode() == kPiPlus) {
532+
histos.fill(HIST("hPx"), posTrack.px());
533+
histos.fill(HIST("hPy"), posTrack.py());
534+
histos.fill(HIST("hPz"), posTrack.pz());
535+
// }
536+
// }
537+
519538
for (auto const& negTrack : negTracksGrouped) {
520539
if (!negTrack.isReconstructed()) {
521540
continue; // no ghost tracks
@@ -571,6 +590,13 @@ struct Alice3strangenessFinder {
571590
// } else {
572591
// histos.fill(HIST("hV0Counter"), 3.5);
573592
// }
593+
histos.fill(HIST("hPxPairs"), posTrack.px());
594+
histos.fill(HIST("hPyPairs"), posTrack.py());
595+
histos.fill(HIST("hPzPairs"), posTrack.pz());
596+
597+
histos.fill(HIST("hPxV0CandDaughter"), v0Cand.pDau0[0]);
598+
histos.fill(HIST("hPyV0CandDaughter"), v0Cand.pDau0[1]);
599+
histos.fill(HIST("hPzV0CandDaughter"), v0Cand.pDau0[2]);
574600

575601
v0CandidateIndices(collision.globalIndex(),
576602
posTrack.globalIndex(),

‎ALICE3/Tasks/alice3Strangeness.cxx‎

Lines changed: 30 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -90,7 +90,7 @@ struct Alice3Strangeness {
9090
Configurable<bool> applyLifetimeSelection{"applyLifetimeSelection", true, "apply lifetime selection"};
9191
Configurable<bool> applyEtaDaughterSelection{"applyEtaDaughterSelection", true, "apply eta daughter selection"};
9292
Configurable<bool> doQAforSelectionVariables{"doQAforSelectionVariables", false, "enable QA plots"};
93-
Configurable<bool> analyseOnlyTrueV0s{"analyseOnlyTrueV0s", false, "analyse only true V0s from MC"};
93+
Configurable<bool> analyseRecoV0Candidates{"analyseRecoV0Candidates", false, "analyse reconstructed V0 candidates"};
9494
} v0SelectionFlags;
9595

9696
struct : ConfigurableGroup {
@@ -299,6 +299,8 @@ struct Alice3Strangeness {
299299
histos.add("reconstructedCandidates/Lambda/hMass1D", "hMass1D", kTH1D, {histAxes.axisLambdaMass});
300300
histos.add("reconstructedCandidates/hArmeterosBeforeAllSelections", "hArmeterosBeforeAllSelections", kTH2D, {{100, -1.0f, 1.0f}, {200, 0.0f, 0.5f}});
301301
histos.add("reconstructedCandidates/hArmeterosAfterAllSelections", "hArmeterosAfterAllSelections", kTH2D, {{100, -1.0f, 1.0f}, {200, 0.0f, 0.5f}});
302+
histos.add("reconstructedCandidates/hPtEtaPosDaughter", "hPtEtaPosDaughter", kTH2D, {histAxes.axisPt, histAxes.axisEta});
303+
histos.add("reconstructedCandidates/hPtEtaNegDaughter", "hPtEtaNegDaughter", kTH2D, {histAxes.axisPt, histAxes.axisEta});
302304

303305
if (doprocessFoundCascadeCandidates) {
304306
analysedCascade.setCascadeType(static_cast<Cascade::Type>(cascadeFlags.analyseCascade.value));
@@ -432,9 +434,9 @@ struct Alice3Strangeness {
432434
auto negV0Daughter = v0Cand.negTrack_as<Alice3Tracks>(); // de-reference negative track
433435
auto posV0Daughter = v0Cand.posTrack_as<Alice3Tracks>(); // de-reference positive track
434436

435-
bool isK0 = v0Cand.mK0() > 0;
437+
bool isK0 = v0Cand.mK0Short() > 0;
436438
if (isK0) {
437-
histos.fill(HIST("K0/hMassAllCandidates"), v0Cand.mK0(), v0Cand.pt());
439+
histos.fill(HIST("K0/hMassAllCandidates"), v0Cand.mK0Short(), v0Cand.pt());
438440
histos.fill(HIST("K0/hSelections"), 0); // all candidates
439441
histos.fill(HIST("K0/hDCANegDaughter"), negV0Daughter.dcaXY());
440442
histos.fill(HIST("K0/hDCAPosDaughter"), posV0Daughter.dcaXY());
@@ -453,7 +455,7 @@ struct Alice3Strangeness {
453455
if (std::abs(negV0Daughter.eta()) > v0SelectionValues.etaDaughterSelection || std::abs(posV0Daughter.eta()) > v0SelectionValues.etaDaughterSelection)
454456
continue;
455457
histos.fill(HIST("K0/hSelections"), 5); // eta cut
456-
histos.fill(HIST("K0/hMassSelected"), v0Cand.mK0(), v0Cand.pt());
458+
histos.fill(HIST("K0/hMassSelected"), v0Cand.mK0Short(), v0Cand.pt());
457459
}
458460
}
459461

@@ -470,7 +472,7 @@ struct Alice3Strangeness {
470472
}
471473
}
472474

473-
void processFoundV0Candidates(aod::Collision const& collision, FullV0Candidates const& v0Candidates, Alice3Tracks const&, aod::McParticles const&)
475+
void processFoundV0Candidates(aod::Collision const& collision, FullV0Candidates const& v0Candidates, Alice3Tracks const&, aod::McParticles const& mcParticles)
474476
{
475477
// if(collision.lutConfigId()!=idGeometry)
476478
// return;
@@ -550,15 +552,29 @@ struct Alice3Strangeness {
550552
selectionCheck = v0.distOverTotMom(collision.posX(), collision.posY(), collision.posZ()) * o2::constants::physics::MassLambda0;
551553
}
552554
histos.fill(HIST("hV0CandidateCounter"), 8.5);
553-
auto posTrack = v0.template posTrack_as<Alice3Tracks>();
554-
auto negTrack = v0.template negTrack_as<Alice3Tracks>();
555-
if (v0SelectionFlags.applyEtaDaughterSelection) {
556-
if (std::abs(posTrack.eta()) > v0SelectionValues.etaDaughterSelection || std::abs(negTrack.eta()) > v0SelectionValues.etaDaughterSelection)
557-
continue;
555+
if (v0SelectionFlags.analyseRecoV0Candidates) {
556+
auto posTrack = v0.template posTrack_as<Alice3Tracks>();
557+
auto negTrack = v0.template negTrack_as<Alice3Tracks>();
558+
if (v0SelectionFlags.applyEtaDaughterSelection) {
559+
if (std::abs(posTrack.eta()) > v0SelectionValues.etaDaughterSelection || std::abs(negTrack.eta()) > v0SelectionValues.etaDaughterSelection)
560+
continue;
561+
}
562+
histos.fill(HIST("reconstructedCandidates/hEtaDaughters"), posTrack.eta());
563+
histos.fill(HIST("reconstructedCandidates/hEtaDaughters"), negTrack.eta());
564+
histos.fill(HIST("hV0CandidateCounter"), 9.5);
565+
histos.fill(HIST("reconstructedCandidates/hPtEtaPosDaughter"), posTrack.pt(), posTrack.eta());
566+
histos.fill(HIST("reconstructedCandidates/hPtEtaNegDaughter"), negTrack.pt(), negTrack.eta());
567+
} else {
568+
int posDaugID = v0.posTrackId();
569+
int negDaugID = v0.negTrackId();
570+
auto posMcParticle = mcParticles.rawIteratorAt(posDaugID);
571+
auto negMcParticle = mcParticles.rawIteratorAt(negDaugID);
572+
// auto posTrack = v0.posTrack_as<o2::aod::McParticle>();
573+
// auto negTrack = v0.negTrack_as<o2::aod::McParticle>();
574+
575+
histos.fill(HIST("reconstructedCandidates/hPtEtaPosDaughter"), posMcParticle.pt(), posMcParticle.eta());
576+
histos.fill(HIST("reconstructedCandidates/hPtEtaNegDaughter"), negMcParticle.pt(), negMcParticle.eta());
558577
}
559-
histos.fill(HIST("reconstructedCandidates/hEtaDaughters"), posTrack.eta());
560-
histos.fill(HIST("reconstructedCandidates/hEtaDaughters"), negTrack.eta());
561-
histos.fill(HIST("hV0CandidateCounter"), 9.5);
562578

563579
histos.fill(HIST("reconstructedCandidates/hArmeterosAfterAllSelections"), v0.alpha(), v0.qtArm());
564580
if (v0SelectionFlags.doQAforSelectionVariables) {
@@ -576,6 +592,7 @@ struct Alice3Strangeness {
576592
}
577593
});
578594
}
595+
579596
if (isK0) {
580597
histos.fill(HIST("reconstructedCandidates/K0/hMass"), v0.mK0Short(), v0.pt(), v0.eta());
581598
histos.fill(HIST("reconstructedCandidates/K0/hMass1D"), v0.mK0Short());

‎Common/Tasks/centralityQa.cxx‎

Lines changed: 44 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -66,6 +66,7 @@ struct CentralityQa {
6666
ConfigurableAxis axisMultiplicity{"axisMultiplicity", {1000, 0, 1000}, "Multiplicity"};
6767
ConfigurableAxis axisMultiplicityPV{"axisMultiplicityPV", {1000, 0, 1000}, "Multiplicity PV"};
6868
ConfigurableAxis axisChannelAmplitude{"axisChannelAmplitude", {5000, 0, 5000}, "Channel Amplitude"};
69+
ConfigurableAxis axisCentrality{"axisCentrality", {101, 0.0f, 101.0f}, "Centrality (%)"};
6970

7071
struct : ConfigurableGroup {
7172
std::string prefix = "eventSelections"; // JSON group name
@@ -257,6 +258,28 @@ struct CentralityQa {
257258
histos.add("hMultEta05VsGenMultNTPV", ";Multiplicity NTPV; Multiplicity PV contributors (|#it{#eta}| < 0.5)", kTH2D, {axisMultiplicity, axisMultiplicityPV});
258259
histos.add("hMultEta05VsGenMultNGlobal", ";Multiplicity NGlobal; Multiplicity PV contributors (|#it{#eta}| < 0.5)", kTH2D, {axisMultiplicity, axisMultiplicityPV});
259260
histos.add("hMultEta05VsGenMultMFT", ";Multiplicity MFT; Multiplicity PV contributors (|#it{#eta}| < 0.5)", kTH2D, {axisMultiplicity, axisMultiplicityPV});
261+
262+
histos.add("hGenMultEta05VsCentralityFV0A", ";FV0A Centrality (%); Generated multiplicity (|#it{#eta}| < 0.5)", kTH2D, {axisCentrality, axisMultiplicityPV});
263+
histos.add("hGenMultEta05VsCentralityFT0M", ";FT0M Centrality (%); Generated multiplicity (|#it{#eta}| < 0.5)", kTH2D, {axisCentrality, axisMultiplicityPV});
264+
histos.add("hGenMultEta05VsCentralityFT0A", ";FT0A Centrality (%); Generated multiplicity (|#it{#eta}| < 0.5)", kTH2D, {axisCentrality, axisMultiplicityPV});
265+
histos.add("hGenMultEta05VsCentralityFT0C", ";FT0C Centrality (%); Generated multiplicity (|#it{#eta}| < 0.5)", kTH2D, {axisCentrality, axisMultiplicityPV});
266+
histos.add("hGenMultEta05VsCentralityFT0CVar1", ";FT0CVar1 Centrality (%); Generated multiplicity (|#it{#eta}| < 0.5)", kTH2D, {axisCentrality, axisMultiplicityPV});
267+
histos.add("hGenMultEta05VsCentralityFT0CVar2", ";FT0CVar2 Centrality (%); Generated multiplicity (|#it{#eta}| < 0.5)", kTH2D, {axisCentrality, axisMultiplicityPV});
268+
histos.add("hGenMultEta05VsCentralityFDDM", ";FDDM Centrality (%); Generated multiplicity (|#it{#eta}| < 0.5)", kTH2D, {axisCentrality, axisMultiplicityPV});
269+
histos.add("hGenMultEta05VsCentralityNTPV", ";NTPV Centrality (%); Generated multiplicity (|#it{#eta}| < 0.5)", kTH2D, {axisCentrality, axisMultiplicityPV});
270+
histos.add("hGenMultEta05VsCentralityNGlobal", ";NGlobal Centrality (%); Generated multiplicity (|#it{#eta}| < 0.5)", kTH2D, {axisCentrality, axisMultiplicityPV});
271+
histos.add("hGenMultEta05VsCentralityMFT", ";MFT Centrality (%); Generated multiplicity (|#it{#eta}| < 0.5)", kTH2D, {axisCentrality, axisMultiplicityPV});
272+
273+
histos.add("hGenMultVsCentralityFV0A", ";FV0A Centrality (%); Generated multiplicity FV0A", kTH2D, {axisCentrality, axisMultiplicityPV});
274+
histos.add("hGenMultVsCentralityFT0M", ";FT0M Centrality (%); Generated multiplicity FT0M", kTH2D, {axisCentrality, axisMultiplicityPV});
275+
histos.add("hGenMultVsCentralityFT0A", ";FT0A Centrality (%); Generated multiplicity FT0A", kTH2D, {axisCentrality, axisMultiplicityPV});
276+
histos.add("hGenMultVsCentralityFT0C", ";FT0C Centrality (%); Generated multiplicity FT0C", kTH2D, {axisCentrality, axisMultiplicityPV});
277+
histos.add("hGenMultVsCentralityFT0CVar1", ";FT0CVar1 Centrality (%); Generated multiplicity FT0C", kTH2D, {axisCentrality, axisMultiplicityPV});
278+
histos.add("hGenMultVsCentralityFT0CVar2", ";FT0CVar2 Centrality (%); Generated multiplicity FT0C", kTH2D, {axisCentrality, axisMultiplicityPV});
279+
histos.add("hGenMultVsCentralityFDDM", ";FDDM Centrality (%); Generated multiplicity FDDM", kTH2D, {axisCentrality, axisMultiplicityPV});
280+
histos.add("hGenMultVsCentralityNTPV", ";NTPV Centrality (%); Generated multiplicity NTPV", kTH2D, {axisCentrality, axisMultiplicityPV});
281+
histos.add("hGenMultVsCentralityNGlobal", ";NGlobal Centrality (%); Generated multiplicity global tracks", kTH2D, {axisCentrality, axisMultiplicityPV});
282+
histos.add("hGenMultVsCentralityMFT", ";MFT Centrality (%); Generated multiplicity MFT", kTH2D, {axisCentrality, axisMultiplicityPV});
260283
}
261284
}
262285

@@ -861,6 +884,8 @@ struct CentralityQa {
861884
histos.fill(HIST("hCentProfileFV0A"), centFV0A, col.multNTracksPVetaHalf());
862885
histos.fill(HIST("hMultEta05VsCentFV0A"), centFV0A, col.multNTracksPVetaHalf());
863886
histos.fill(HIST("hMultEta05VsGenMultFV0A"), mcCol.multMCFV0A(), col.multNTracksPVetaHalf());
887+
histos.fill(HIST("hGenMultEta05VsCentralityFV0A"), col.centFV0A(), mcCol.multMCNParticlesEta05());
888+
histos.fill(HIST("hGenMultVsCentralityFV0A"), col.centFV0A(), mcCol.multMCFV0A());
864889
}
865890
PROCESS_SWITCH(CentralityQa, processMonteCarloRun3_FV0A, "Process with Run 3 FV0A estimator", false);
866891

@@ -881,6 +906,8 @@ struct CentralityQa {
881906
histos.fill(HIST("hCentProfileFT0M"), centFT0M, col.multNTracksPVetaHalf());
882907
histos.fill(HIST("hMultEta05VsCentFT0M"), centFT0M, col.multNTracksPVetaHalf());
883908
histos.fill(HIST("hMultEta05VsGenMultFT0M"), mcCol.multMCFT0A() + mcCol.multMCFT0C(), col.multNTracksPVetaHalf());
909+
histos.fill(HIST("hGenMultEta05VsCentralityFT0M"), col.centFT0M(), mcCol.multMCNParticlesEta05());
910+
histos.fill(HIST("hGenMultVsCentralityFT0M"), col.centFT0M(), mcCol.multMCFT0A() + mcCol.multMCFT0C());
884911
}
885912
PROCESS_SWITCH(CentralityQa, processMonteCarloRun3_FT0M, "Process with Run 3 FT0M estimator", false);
886913

@@ -901,6 +928,8 @@ struct CentralityQa {
901928
histos.fill(HIST("hCentProfileFT0A"), centFT0A, col.multNTracksPVetaHalf());
902929
histos.fill(HIST("hMultEta05VsCentFT0A"), centFT0A, col.multNTracksPVetaHalf());
903930
histos.fill(HIST("hMultEta05VsGenMultFT0A"), mcCol.multMCFT0A(), col.multNTracksPVetaHalf());
931+
histos.fill(HIST("hGenMultEta05VsCentralityFT0A"), col.centFT0A(), mcCol.multMCNParticlesEta05());
932+
histos.fill(HIST("hGenMultVsCentralityFT0A"), col.centFT0A(), mcCol.multMCFT0A());
904933
}
905934
PROCESS_SWITCH(CentralityQa, processMonteCarloRun3_FT0A, "Process with Run 3 FT0A estimator", false);
906935

@@ -921,6 +950,8 @@ struct CentralityQa {
921950
histos.fill(HIST("hCentProfileFT0C"), centFT0C, col.multNTracksPVetaHalf());
922951
histos.fill(HIST("hMultEta05VsCentFT0C"), centFT0C, col.multNTracksPVetaHalf());
923952
histos.fill(HIST("hMultEta05VsGenMultFT0C"), mcCol.multMCFT0C(), col.multNTracksPVetaHalf());
953+
histos.fill(HIST("hGenMultEta05VsCentralityFT0C"), col.centFT0C(), mcCol.multMCNParticlesEta05());
954+
histos.fill(HIST("hGenMultVsCentralityFT0C"), col.centFT0C(), mcCol.multMCFT0C());
924955
}
925956
PROCESS_SWITCH(CentralityQa, processMonteCarloRun3_FT0C, "Process with Run 3 FT0C estimator", false);
926957

@@ -941,6 +972,8 @@ struct CentralityQa {
941972
histos.fill(HIST("hCentProfileFT0CVar1"), centFT0Cvar1, col.multNTracksPVetaHalf());
942973
histos.fill(HIST("hMultEta05VsCentFT0CVar1"), centFT0Cvar1, col.multNTracksPVetaHalf());
943974
histos.fill(HIST("hMultEta05VsGenMultFT0CVar1"), mcCol.multMCFT0C(), col.multNTracksPVetaHalf());
975+
histos.fill(HIST("hGenMultEta05VsCentralityFT0CVar1"), col.centFT0CVariant1(), mcCol.multMCNParticlesEta05());
976+
histos.fill(HIST("hGenMultVsCentralityFT0CVar1"), col.centFT0CVariant1(), mcCol.multMCFT0C());
944977
}
945978
PROCESS_SWITCH(CentralityQa, processMonteCarloRun3_FT0CVar1, "Process with Run 3 FT0CVar1 estimator", false);
946979

@@ -961,6 +994,8 @@ struct CentralityQa {
961994
histos.fill(HIST("hCentProfileFT0CVar2"), centFT0Cvar2, col.multNTracksPVetaHalf());
962995
histos.fill(HIST("hMultEta05VsCentFT0CVar2"), centFT0Cvar2, col.multNTracksPVetaHalf());
963996
histos.fill(HIST("hMultEta05VsGenMultFT0CVar2"), mcCol.multMCFT0C(), col.multNTracksPVetaHalf());
997+
histos.fill(HIST("hGenMultEta05VsCentralityFT0CVar2"), col.centFT0CVariant2(), mcCol.multMCNParticlesEta05());
998+
histos.fill(HIST("hGenMultVsCentralityFT0CVar2"), col.centFT0CVariant2(), mcCol.multMCFT0C());
964999
}
9651000
PROCESS_SWITCH(CentralityQa, processMonteCarloRun3_FT0CVar2, "Process with Run 3 FT0CVar2 estimator", false);
9661001

@@ -981,6 +1016,8 @@ struct CentralityQa {
9811016
histos.fill(HIST("hCentProfileFDDM"), centFDDM, col.multNTracksPVetaHalf());
9821017
histos.fill(HIST("hMultEta05VsCentFDDM"), centFDDM, col.multNTracksPVetaHalf());
9831018
histos.fill(HIST("hMultEta05VsGenMultFDDM"), mcCol.multMCFDDA() + mcCol.multMCFDDC(), col.multNTracksPVetaHalf());
1019+
histos.fill(HIST("hGenMultEta05VsCentralityFDDM"), col.centFDDM(), mcCol.multMCNParticlesEta05());
1020+
histos.fill(HIST("hGenMultVsCentralityFDDM"), col.centFDDM(), mcCol.multMCFDDA() + mcCol.multMCFDDC());
9841021
}
9851022
PROCESS_SWITCH(CentralityQa, processMonteCarloRun3_FDDM, "Process with Run 3 FDDM estimator", false);
9861023

@@ -1000,6 +1037,8 @@ struct CentralityQa {
10001037
histos.fill(HIST("hCentProfileNTPV"), centNTPV, col.multNTracksPVetaHalf());
10011038
histos.fill(HIST("hMultEta05VsCentNTPV"), centNTPV, col.multNTracksPVetaHalf());
10021039
histos.fill(HIST("hMultEta05VsGenMultNTPV"), mcCol.multMCNParticlesEta08(), col.multNTracksPVetaHalf());
1040+
histos.fill(HIST("hGenMultEta05VsCentralityNTPV"), col.centNTPV(), mcCol.multMCNParticlesEta05());
1041+
histos.fill(HIST("hGenMultVsCentralityFDDM"), col.centNTPV(), mcCol.multMCNParticlesEta08());
10031042
}
10041043
PROCESS_SWITCH(CentralityQa, processMonteCarloRun3_NTPV, "Process with Run 3 NTPV estimator", false);
10051044

@@ -1020,6 +1059,8 @@ struct CentralityQa {
10201059
histos.fill(HIST("hCentProfileNGlobal"), centNGlo, col.multNTracksPVetaHalf());
10211060
histos.fill(HIST("hMultEta05VsCentNGlobal"), centNGlo, col.multNTracksPVetaHalf());
10221061
histos.fill(HIST("hMultEta05VsGenMultNGlobal"), mcCol.multMCNParticlesEta08(), col.multNTracksPVetaHalf());
1062+
histos.fill(HIST("hGenMultEta05VsCentralityNGlobal"), col.centNGlobal(), mcCol.multMCNParticlesEta05());
1063+
histos.fill(HIST("hGenMultVsCentralityNGlobal"), col.centNGlobal(), mcCol.multMCNParticlesEta08());
10231064
}
10241065
PROCESS_SWITCH(CentralityQa, processMonteCarloRun3_NGlobal, "Process with Run 3 NGlobal estimator", false);
10251066

@@ -1031,7 +1072,7 @@ struct CentralityQa {
10311072
return;
10321073
}
10331074

1034-
// const auto& mcCol = col.mcCollision_as<soa::Join<aod::McCollisions, aod::MultMCExtras>>(); // FIXME: uncomment when MC MFT mult is added in aod::MultMCExtras
1075+
const auto& mcCol = col.mcCollision_as<soa::Join<aod::McCollisions, aod::MultMCExtras>>(); // FIXME: uncomment when MC MFT mult is added in aod::MultMCExtras
10351076
Estimator mft = initEstimator(col, "MFT");
10361077
const float centMFT = mft.getCentrality(col.mftNtracks(), col.centMFT());
10371078

@@ -1040,6 +1081,8 @@ struct CentralityQa {
10401081
histos.fill(HIST("hCentProfileMFT"), centMFT, col.multNTracksPVetaHalf());
10411082
histos.fill(HIST("hMultEta05VsCentMFT"), centMFT, col.multNTracksPVetaHalf());
10421083
// histos.fill(HIST("hMultEta05VsGenMultMFT"), mcCol.multMCMFT(), col.multNTracksPVetaHalf()); // FIXME: uncomment when MC MFT mult is added in aod::MultMCExtras
1084+
histos.fill(HIST("hGenMultEta05VsCentralityMFT"), col.centMFT(), mcCol.multMCNParticlesEta05());
1085+
// histos.fill(HIST("hGenMultVsCentralityMFT"), col.centMFT(), mcCol.multMCMFT()); // FIXME: uncomment when MC MFT mult is added in aod::MultMCExtras
10431086
}
10441087
PROCESS_SWITCH(CentralityQa, processMonteCarloRun3_MFT, "Process with Run 3 MFT estimator", false);
10451088

0 commit comments

Comments
 (0)