Skip to content

Commit a118602

Browse files
JStaaJStaa
committed
Fixed bugs in processMCSame and processMCMixed functions
1 parent 10cef7b commit a118602

File tree

1 file changed

+57
-48
lines changed

1 file changed

+57
-48
lines changed

‎PWGCF/MultiparticleCorrelations/Tasks/threeParticleCorrelations.cxx‎

Lines changed: 57 additions & 48 deletions
Original file line numberDiff line numberDiff line change
@@ -26,6 +26,8 @@
2626
#include "RecoDecay.h"
2727
#include "TPDGCode.h"
2828

29+
#include "Common/DataModel/McCollisionExtra.h"
30+
2931
using namespace o2;
3032
using namespace o2::framework;
3133
using namespace o2::framework::expressions;
@@ -41,7 +43,7 @@ struct ThreeParticleCorrelations {
4143
HistogramRegistry rQARegistry{"QARegistry", {}, OutputObjHandlingPolicy::AnalysisObject, false, true};
4244

4345
// Collision & Event filters
44-
Filter collCent = aod::cent::centFT0M > 0.0f && aod::cent::centFT0M < 90.0f;
46+
Filter collCent = aod::cent::centFT0C > 0.0f && aod::cent::centFT0C < 90.0f;
4547
Filter collZvtx = nabs(aod::collision::posZ) < 7.0f;
4648
Filter mcCollZvtx = nabs(aod::mccollision::posZ) < 7.0f;
4749
Filter evSelect = aod::evsel::sel8 == true;
@@ -59,15 +61,15 @@ struct ThreeParticleCorrelations {
5961
Filter particleEta = nabs(aod::mcparticle::eta) < 0.8f;
6062

6163
// Table aliases - Data
62-
using MyFilteredCollisions = soa::Filtered<soa::Join<aod::Collisions, aod::CentFT0Ms, aod::EvSels>>;
64+
using MyFilteredCollisions = soa::Filtered<soa::Join<aod::Collisions, aod::CentFT0Cs, aod::EvSels>>;
6365
using MyFilteredCollision = MyFilteredCollisions::iterator;
6466
using MyFilteredV0s = soa::Filtered<aod::V0Datas>;
6567
using MyFilteredTracks = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TrackSelection,
6668
aod::pidTPCPi, aod::pidTPCKa, aod::pidTPCPr,
6769
aod::pidTOFFullPi, aod::pidTOFFullKa, aod::pidTOFFullPr, aod::pidTOFbeta>>;
6870

6971
// Table aliases - MC
70-
using MyFilteredMCGenCollisions = soa::Filtered<soa::Join<aod::McCollisions, aod::CentFT0Ms>>;
72+
using MyFilteredMCGenCollisions = soa::Filtered<soa::Join<aod::McCollisions, aod::McCollsExtra>>;
7173
using MyFilteredMCGenCollision = MyFilteredMCGenCollisions::iterator;
7274
using MyFilteredMCParticles = soa::Filtered<aod::McParticles>;
7375
using MyFilteredMCRecCollision = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels>>::iterator;
@@ -79,11 +81,11 @@ struct ThreeParticleCorrelations {
7981
SliceCache cache;
8082
ConfigurableAxis confCentBins{"confCentBins", {VARIABLE_WIDTH, 0.0f, 10.0f, 20.0f, 30.0f, 40.0f, 50.0f, 60.0f, 70.0f, 80.0f, 90.0f}, "ME Centrality binning"};
8183
ConfigurableAxis confZvtxBins{"confZvtxBins", {VARIABLE_WIDTH, -7.0f, -5.0f, -3.0f, -1.0f, 0.0f, 1.0f, 3.0f, 5.0f, 7.0f}, "ME Zvtx binning"};
82-
using BinningType = ColumnBinningPolicy<aod::cent::CentFT0M, aod::collision::PosZ>;
83-
84+
85+
using BinningType = ColumnBinningPolicy<aod::cent::CentFT0C, aod::collision::PosZ>;
86+
using BinningTypeMC = ColumnBinningPolicy<aod::mccollisionprop::BestCollisionCentFT0C, aod::mccollision::PosZ>;
8487
BinningType collBinning{{confCentBins, confZvtxBins}, true};
85-
Pair<MyFilteredCollisions, MyFilteredV0s, MyFilteredTracks, BinningType> pairData{collBinning, 5, -1, &cache};
86-
SameKindPair<MyFilteredMCGenCollisions, MyFilteredMCParticles, BinningType> pairMC{collBinning, 5, -1, &cache};
88+
BinningTypeMC collBinningMC{{confCentBins, confZvtxBins}, true};
8789

8890
// Process configurables
8991
Configurable<bool> confFilterSwitch{"confFilterSwitch", false, "Switch for the fakeV0Filter function"};
@@ -123,6 +125,7 @@ struct ThreeParticleCorrelations {
123125
rQARegistry.add("hTrackEta", "hTrackEta", {HistType::kTH1D, {{100, -1, 1}}});
124126
rQARegistry.add("hTrackPhi", "hTrackPhi", {HistType::kTH1D, {{100, (-1. / 2) * constants::math::PI, (5. / 2) * constants::math::PI}}});
125127
rQARegistry.add("hEventCentrality", "hEventCentrality", {HistType::kTH1D, {{centralityAxis}}});
128+
rQARegistry.add("hEventCentrality_MC", "hEventCentrality_MC", {HistType::kTH1D, {{centralityAxis}}});
126129
rQARegistry.add("hEventZvtx", "hEventZvtx", {HistType::kTH1D, {{zvtxAxis}}});
127130

128131
rQARegistry.add("hdEdx", "hdEdx", {HistType::kTH2D, {{56, 0.2, 3.0}, {180, 20, 200}}});
@@ -146,7 +149,7 @@ struct ThreeParticleCorrelations {
146149

147150
rQARegistry.add("hInvMassLambda", "hInvMassLambda", {HistType::kTH3D, {{lambdaInvMassAxis}, {v0PtAxis}, {centralityAxis}}});
148151
rQARegistry.add("hInvMassAntiLambda", "hInvMassAntiLambda", {HistType::kTH3D, {{lambdaInvMassAxis}, {v0PtAxis}, {centralityAxis}}});
149-
rQARegistry.add("hNLambdas_MC", "hNLambdas_MC", {HistType::kTH3D, {{2, -2, 2}, {v0PtAxis}, {centralityAxis}}});
152+
rQARegistry.add("hNLambdas", "hNLambdas", {HistType::kTH3D, {{2, -2, 2}, {v0PtAxis}, {centralityAxis}}});
150153

151154
// PhiStar
152155
rPhiStarRegistry.add("hSEProtonPreCut", "hSEProtonPreCut", {HistType::kTH2D, {{80, -0.2, 0.2}, {40, -0.1, 0.1}}});
@@ -223,7 +226,7 @@ struct ThreeParticleCorrelations {
223226

224227
auto bc = collision.bc_as<aod::BCsWithTimestamps>();
225228
auto bField = getMagneticField(bc.timestamp());
226-
rQARegistry.fill(HIST("hEventCentrality"), collision.centFT0M());
229+
rQARegistry.fill(HIST("hEventCentrality"), collision.centFT0C());
227230
rQARegistry.fill(HIST("hEventZvtx"), collision.posZ());
228231

229232
// Start of the Track QA
@@ -261,17 +264,17 @@ struct ThreeParticleCorrelations {
261264
}
262265
// End of the Track QA
263266

264-
// Start of the Same-Event Correlations
267+
// Start of the Same-Event correlations
265268
for (const auto& trigger : v0s) {
266269
if (v0Filters(trigger)) {
267270

268271
triggSign = v0Sign(trigger);
269272
if (triggSign == 1) {
270273
candMass = trigger.mLambda();
271-
rQARegistry.fill(HIST("hInvMassLambda"), trigger.mLambda(), trigger.pt(), collision.centFT0M());
274+
rQARegistry.fill(HIST("hInvMassLambda"), trigger.mLambda(), trigger.pt(), collision.centFT0C());
272275
} else if (triggSign == -1) {
273276
candMass = trigger.mAntiLambda();
274-
rQARegistry.fill(HIST("hInvMassAntiLambda"), trigger.mAntiLambda(), trigger.pt(), collision.centFT0M());
277+
rQARegistry.fill(HIST("hInvMassAntiLambda"), trigger.mAntiLambda(), trigger.pt(), collision.centFT0C());
275278
}
276279

277280
for (const auto& associate : tracks) {
@@ -284,33 +287,36 @@ struct ThreeParticleCorrelations {
284287

285288
if (candMass >= massLambda - 4 * dGaussSigma && candMass <= massLambda + 4 * dGaussSigma) {
286289
if (assocPID[0] == 0.0) { // Pions
287-
rSECorrRegistry.fill(HIST("hSameLambdaPion_SGNL"), deltaPhi, deltaEta, collision.centFT0M(), collision.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffPions, associate.sign(), associate.pt()));
290+
rSECorrRegistry.fill(HIST("hSameLambdaPion_SGNL"), deltaPhi, deltaEta, collision.centFT0C(), collision.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffPions, associate.sign(), associate.pt()));
288291
} else if (assocPID[0] == 1.0) { // Kaons
289-
rSECorrRegistry.fill(HIST("hSameLambdaKaon_SGNL"), deltaPhi, deltaEta, collision.centFT0M(), collision.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffKaons, associate.sign(), associate.pt()));
292+
rSECorrRegistry.fill(HIST("hSameLambdaKaon_SGNL"), deltaPhi, deltaEta, collision.centFT0C(), collision.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffKaons, associate.sign(), associate.pt()));
290293
} else if (assocPID[0] == 2.0) { // Protons
291-
rSECorrRegistry.fill(HIST("hSameLambdaProton_SGNL"), deltaPhi, deltaEta, collision.centFT0M(), collision.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffProtons, associate.sign(), associate.pt()));
294+
rSECorrRegistry.fill(HIST("hSameLambdaProton_SGNL"), deltaPhi, deltaEta, collision.centFT0C(), collision.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffProtons, associate.sign(), associate.pt()));
292295
}
293296
} else if (candMass >= massLambda - 8 * dGaussSigma && candMass <= massLambda + 8 * dGaussSigma) {
294297
if (assocPID[0] == 0.0) { // Pions
295-
rSECorrRegistry.fill(HIST("hSameLambdaPion_SB"), deltaPhi, deltaEta, collision.centFT0M(), collision.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffPions, associate.sign(), associate.pt()));
298+
rSECorrRegistry.fill(HIST("hSameLambdaPion_SB"), deltaPhi, deltaEta, collision.centFT0C(), collision.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffPions, associate.sign(), associate.pt()));
296299
} else if (assocPID[0] == 1.0) { // Kaons
297-
rSECorrRegistry.fill(HIST("hSameLambdaKaon_SB"), deltaPhi, deltaEta, collision.centFT0M(), collision.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffKaons, associate.sign(), associate.pt()));
300+
rSECorrRegistry.fill(HIST("hSameLambdaKaon_SB"), deltaPhi, deltaEta, collision.centFT0C(), collision.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffKaons, associate.sign(), associate.pt()));
298301
} else if (assocPID[0] == 2.0) { // Protons
299-
rSECorrRegistry.fill(HIST("hSameLambdaProton_SB"), deltaPhi, deltaEta, collision.centFT0M(), collision.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffProtons, associate.sign(), associate.pt()));
302+
rSECorrRegistry.fill(HIST("hSameLambdaProton_SB"), deltaPhi, deltaEta, collision.centFT0C(), collision.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffProtons, associate.sign(), associate.pt()));
300303
}
301304
}
302305
}
303306
}
304307
}
305308
}
306309
}
307-
// End of the Same-Event Correlations
310+
// End of the Same-Event correlations
308311
}
309312

310-
void processMixed(MyFilteredCollisions const&, MyFilteredV0s const&, MyFilteredTracks const&, aod::BCsWithTimestamps const&)
313+
void processMixed(MyFilteredCollisions const& collisions, MyFilteredV0s const& v0s, MyFilteredTracks const& tracks, aod::BCsWithTimestamps const&)
311314
{
312315

313-
// Start of the Mixed-events Correlations
316+
auto tracksTuple = std::make_tuple(v0s, tracks);
317+
Pair<MyFilteredCollisions, MyFilteredV0s, MyFilteredTracks, BinningType> pairData{collBinning, 5, -1, collisions, tracksTuple, &cache};
318+
319+
// Start of the Mixed-Event correlations
314320
for (const auto& [coll_1, v0_1, coll_2, track_2] : pairData) {
315321

316322
auto bc = coll_1.bc_as<aod::BCsWithTimestamps>();
@@ -332,51 +338,51 @@ struct ThreeParticleCorrelations {
332338

333339
if (candMass >= massLambda - 4 * dGaussSigma && candMass <= massLambda + 4 * dGaussSigma) {
334340
if (assocPID[0] == 0.0) { // Pions
335-
rMECorrRegistry.fill(HIST("hMixLambdaPion_SGNL"), deltaPhi, deltaEta, coll_1.centFT0M(), coll_1.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffPions, associate.sign(), associate.pt()));
341+
rMECorrRegistry.fill(HIST("hMixLambdaPion_SGNL"), deltaPhi, deltaEta, coll_1.centFT0C(), coll_1.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffPions, associate.sign(), associate.pt()));
336342
} else if (assocPID[0] == 1.0) { // Kaons
337-
rMECorrRegistry.fill(HIST("hMixLambdaKaon_SGNL"), deltaPhi, deltaEta, coll_1.centFT0M(), coll_1.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffKaons, associate.sign(), associate.pt()));
343+
rMECorrRegistry.fill(HIST("hMixLambdaKaon_SGNL"), deltaPhi, deltaEta, coll_1.centFT0C(), coll_1.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffKaons, associate.sign(), associate.pt()));
338344
} else if (assocPID[0] == 2.0) { // Protons
339-
rMECorrRegistry.fill(HIST("hMixLambdaProton_SGNL"), deltaPhi, deltaEta, coll_1.centFT0M(), coll_1.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffProtons, associate.sign(), associate.pt()));
345+
rMECorrRegistry.fill(HIST("hMixLambdaProton_SGNL"), deltaPhi, deltaEta, coll_1.centFT0C(), coll_1.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffProtons, associate.sign(), associate.pt()));
340346
}
341347
} else if (candMass >= massLambda - 8 * dGaussSigma && candMass <= massLambda + 8 * dGaussSigma) {
342348
if (assocPID[0] == 0.0) { // Pions
343-
rMECorrRegistry.fill(HIST("hMixLambdaPion_SB"), deltaPhi, deltaEta, coll_1.centFT0M(), coll_1.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffPions, associate.sign(), associate.pt()));
349+
rMECorrRegistry.fill(HIST("hMixLambdaPion_SB"), deltaPhi, deltaEta, coll_1.centFT0C(), coll_1.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffPions, associate.sign(), associate.pt()));
344350
} else if (assocPID[0] == 1.0) { // Kaons
345-
rMECorrRegistry.fill(HIST("hMixLambdaKaon_SB"), deltaPhi, deltaEta, coll_1.centFT0M(), coll_1.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffKaons, associate.sign(), associate.pt()));
351+
rMECorrRegistry.fill(HIST("hMixLambdaKaon_SB"), deltaPhi, deltaEta, coll_1.centFT0C(), coll_1.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffKaons, associate.sign(), associate.pt()));
346352
} else if (assocPID[0] == 2.0) { // Protons
347-
rMECorrRegistry.fill(HIST("hMixLambdaProton_SB"), deltaPhi, deltaEta, coll_1.centFT0M(), coll_1.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffProtons, associate.sign(), associate.pt()));
353+
rMECorrRegistry.fill(HIST("hMixLambdaProton_SB"), deltaPhi, deltaEta, coll_1.centFT0C(), coll_1.posZ(), triggSign, associate.sign(), 1. / trackEff(hEffProtons, associate.sign(), associate.pt()));
348354
}
349355
}
350356
}
351357
}
352358
}
353359
}
354-
// End of the Mixed-events Correlations
360+
// End of the Mixed-Event Correlations
355361
}
356362

357363
void processMCSame(MyFilteredMCGenCollision const& collision, MyFilteredMCParticles const& particles)
358364
{
359365

366+
rQARegistry.fill(HIST("hEventCentrality_MC"), collision.bestCollisionCentFT0C());
360367
Partition<MyFilteredMCParticles> mcTriggers = (aod::mcparticle::pdgCode == static_cast<int>(kLambda0) || aod::mcparticle::pdgCode == static_cast<int>(kLambda0Bar)) && aod::mcparticle::pt > 0.6f && aod::mcparticle::pt < 12.0f && nabs(aod::mcparticle::eta) < 0.72f;
361368
Partition<MyFilteredMCParticles> mcAssociates = (((aod::mcparticle::pdgCode == static_cast<int>(kPiPlus) || aod::mcparticle::pdgCode == static_cast<int>(kPiMinus)) && aod::mcparticle::pt > 0.3f && aod::mcparticle::pt < 2.3f) ||
362369
((aod::mcparticle::pdgCode == static_cast<int>(kKPlus) || aod::mcparticle::pdgCode == static_cast<int>(kKMinus)) && aod::mcparticle::pt > 0.5f && aod::mcparticle::pt < 2.5f) ||
363-
((aod::mcparticle::pdgCode == static_cast<int>(kProton) || aod::mcparticle::pdgCode == static_cast<int>(kProtonBar)) && aod::mcparticle::pt > 0.5f));
370+
((aod::mcparticle::pdgCode == static_cast<int>(kProton) || aod::mcparticle::pdgCode == static_cast<int>(kProtonBar)) && aod::mcparticle::pt > 0.5f));
364371
mcTriggers.bindTable(particles);
365372
mcAssociates.bindTable(particles);
366-
367-
// Start of the MC Same-Event Correlations
368-
for (const auto& trigger : mcTriggers) {
373+
374+
// Start of the MC Same-Event correlations
375+
for (const auto& trigger : mcTriggers) {
369376
if (trigger.isPhysicalPrimary()) {
370-
377+
371378
if (trigger.pdgCode() > 0) {
372379
triggSign = 1;
373-
rQARegistry.fill(HIST("hNLambdas_MC"), 1, trigger.pt(), collision.centFT0M());
374380
} else if (trigger.pdgCode() < 0) {
375381
triggSign = -1;
376-
rQARegistry.fill(HIST("hNLambdas_MC"), -1, trigger.pt(), collision.centFT0M());
377382
}
383+
rQARegistry.fill(HIST("hNLambdas"), triggSign, trigger.pt(), collision.bestCollisionCentFT0C());
378384

379-
for (const auto& associate : mcAssociates) {
385+
for (const auto& associate : mcAssociates) {
380386
if (associate.isPhysicalPrimary()) {
381387

382388
if (associate.pdgCode() > 0) {
@@ -389,22 +395,25 @@ struct ThreeParticleCorrelations {
389395
deltaEta = trigger.eta() - associate.eta();
390396

391397
if (std::abs(associate.pdgCode()) == kPiPlus) {
392-
rSECorrRegistry.fill(HIST("hSameLambdaPion_MC"), deltaPhi, deltaEta, collision.centFT0M(), collision.posZ(), triggSign, assocSign);
398+
rSECorrRegistry.fill(HIST("hSameLambdaPion_MC"), deltaPhi, deltaEta, collision.bestCollisionCentFT0C(), collision.posZ(), triggSign, assocSign);
393399
} else if (std::abs(associate.pdgCode()) == kKPlus) {
394-
rSECorrRegistry.fill(HIST("hSameLambdaKaon_MC"), deltaPhi, deltaEta, collision.centFT0M(), collision.posZ(), triggSign, assocSign);
400+
rSECorrRegistry.fill(HIST("hSameLambdaKaon_MC"), deltaPhi, deltaEta, collision.bestCollisionCentFT0C(), collision.posZ(), triggSign, assocSign);
395401
} else if (std::abs(associate.pdgCode()) == kProton) {
396-
rSECorrRegistry.fill(HIST("hSameLambdaProton_MC"), deltaPhi, deltaEta, collision.centFT0M(), collision.posZ(), triggSign, assocSign);
402+
rSECorrRegistry.fill(HIST("hSameLambdaProton_MC"), deltaPhi, deltaEta, collision.bestCollisionCentFT0C(), collision.posZ(), triggSign, assocSign);
397403
}
398-
}
399-
}
404+
}
405+
}
400406
}
401407
}
402408
// End of the MC Same-Event Correlations
403409
}
404410

405-
void processMCMixed(MyFilteredMCGenCollisions const&, MyFilteredMCParticles const&)
411+
void processMCMixed(MyFilteredMCGenCollisions const& collisions, MyFilteredMCParticles const& particles)
406412
{
407413

414+
auto particlesTuple = std::make_tuple(particles);
415+
SameKindPair<MyFilteredMCGenCollisions, MyFilteredMCParticles, BinningTypeMC> pairMC{collBinningMC, 5, -1, collisions, particlesTuple, &cache};
416+
408417
// Start of the MC Mixed-events Correlations
409418
for (const auto& [coll_1, v0_1, coll_2, track_2] : pairMC) {
410419
Partition<MyFilteredMCParticles> mcTriggers = (aod::mcparticle::pdgCode == static_cast<int>(kLambda0) || aod::mcparticle::pdgCode == static_cast<int>(kLambda0Bar)) && aod::mcparticle::pt > 0.6f && aod::mcparticle::pt < 12.0f && nabs(aod::mcparticle::eta) < 0.72f;
@@ -432,11 +441,11 @@ struct ThreeParticleCorrelations {
432441
deltaEta = trigger.eta() - associate.eta();
433442

434443
if (std::abs(associate.pdgCode()) == kPiPlus) {
435-
rMECorrRegistry.fill(HIST("hMixLambdaPion_MC"), deltaPhi, deltaEta, coll_1.centFT0M(), coll_1.posZ(), triggSign, assocSign);
444+
rMECorrRegistry.fill(HIST("hMixLambdaPion_MC"), deltaPhi, deltaEta, coll_1.bestCollisionCentFT0C(), coll_1.posZ(), triggSign, assocSign);
436445
} else if (std::abs(associate.pdgCode()) == kKPlus) {
437-
rMECorrRegistry.fill(HIST("hMixLambdaKaon_MC"), deltaPhi, deltaEta, coll_1.centFT0M(), coll_1.posZ(), triggSign, assocSign);
446+
rMECorrRegistry.fill(HIST("hMixLambdaKaon_MC"), deltaPhi, deltaEta, coll_1.bestCollisionCentFT0C(), coll_1.posZ(), triggSign, assocSign);
438447
} else if (std::abs(associate.pdgCode()) == kProton) {
439-
rMECorrRegistry.fill(HIST("hMixLambdaProton_MC"), deltaPhi, deltaEta, coll_1.centFT0M(), coll_1.posZ(), triggSign, assocSign);
448+
rMECorrRegistry.fill(HIST("hMixLambdaProton_MC"), deltaPhi, deltaEta, coll_1.bestCollisionCentFT0C(), coll_1.posZ(), triggSign, assocSign);
440449
}
441450
}
442451
}
@@ -449,7 +458,7 @@ struct ThreeParticleCorrelations {
449458

450459
Partition<MyFilteredMCParticles> mcParticles = aod::mcparticle::pt > 0.2f && aod::mcparticle::pt < 3.0f;
451460
mcParticles.bindTable(particles);
452-
461+
453462
// Start of the Monte-Carlo generated QA
454463
for (const auto& particle : mcParticles) {
455464
if (particle.isPhysicalPrimary()) {
@@ -553,8 +562,8 @@ struct ThreeParticleCorrelations {
553562

554563
PROCESS_SWITCH(ThreeParticleCorrelations, processSame, "Process same-event correlations", true);
555564
PROCESS_SWITCH(ThreeParticleCorrelations, processMixed, "Process mixed-event correlations", true);
556-
PROCESS_SWITCH(ThreeParticleCorrelations, processMCSame, "Process MC same-event correlations", true);
557-
PROCESS_SWITCH(ThreeParticleCorrelations, processMCMixed, "Process MC mixed-event correlations", true);
565+
PROCESS_SWITCH(ThreeParticleCorrelations, processMCSame, "Process MC same-event correlations", false);
566+
PROCESS_SWITCH(ThreeParticleCorrelations, processMCMixed, "Process MC mixed-event correlations", false);
558567
PROCESS_SWITCH(ThreeParticleCorrelations, processMCGen, "Process Monte-Carlo, generator level", false);
559568
PROCESS_SWITCH(ThreeParticleCorrelations, processMCRec, "Process Monte-Carlo, reconstructed level", false);
560569

0 commit comments

Comments
 (0)