@@ -68,7 +68,7 @@ struct ThreeParticleCorrelations {
6868
6969 // Table aliases - MC
7070 using MyFilteredMCGenCollisions = soa::Filtered<soa::Join<aod::McCollisions, aod::CentFT0Ms>>;
71- using MyFilteredMCGenCollision = MyFilteredMCGenCollisions::iterator;
71+ using MyFilteredMCGenCollision = MyFilteredMCGenCollisions::iterator;
7272 using MyFilteredMCParticles = soa::Filtered<aod::McParticles>;
7373 using MyFilteredMCRecCollision = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels>>::iterator;
7474 using MyFilteredMCTracks = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TrackSelection, aod::McTrackLabels,
@@ -357,46 +357,46 @@ struct ThreeParticleCorrelations {
357357 void processMCSame (MyFilteredMCGenCollision const & collision, MyFilteredMCParticles const & particles)
358358 {
359359
360- 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 ;
360+ 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 ;
361361 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 ) ||
362- ((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 ));
362+ ((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 ));
364364 mcTriggers.bindTable (particles);
365365 mcAssociates.bindTable (particles);
366-
366+
367367 // Start of the MC Same-Event Correlations
368368 for (const auto & trigger : mcTriggers) {
369369 if (trigger.isPhysicalPrimary ()) {
370-
371- if (trigger.pdgCode () > 0 ) {
372- triggSign = 1 ;
373- rQARegistry.fill (HIST (" hNLambdas_MC" ), 1 , trigger.pt (), collision.centFT0M ());
374- } else if (trigger.pdgCode () < 0 ) {
375- triggSign = -1 ;
376- rQARegistry.fill (HIST (" hNLambdas_MC" ), -1 , trigger.pt (), collision.centFT0M ());
377- }
378-
379- for (const auto & associate : mcAssociates) {
380- if (associate.isPhysicalPrimary ()) {
381-
382- if (associate.pdgCode () > 0 ) {
383- assocSign = 1 ;
384- } else if (associate.pdgCode () < 0 ) {
385- assocSign = -1 ;
386- }
387-
388- deltaPhi = RecoDecay::constrainAngle (trigger.phi () - associate.phi (), -constants::math::PIHalf);
389- deltaEta = trigger.eta () - associate.eta ();
390-
391- if (std::abs (associate.pdgCode ()) == kPiPlus ) {
392- rSECorrRegistry.fill (HIST (" hSameLambdaPion_MC" ), deltaPhi, deltaEta, collision.centFT0M (), collision.posZ (), triggSign, assocSign);
393- } else if (std::abs (associate.pdgCode ()) == kKPlus ) {
394- rSECorrRegistry.fill (HIST (" hSameLambdaKaon_MC" ), deltaPhi, deltaEta, collision.centFT0M (), collision.posZ (), triggSign, assocSign);
395- } else if (std::abs (associate.pdgCode ()) == kProton ) {
396- rSECorrRegistry.fill (HIST (" hSameLambdaProton_MC" ), deltaPhi, deltaEta, collision.centFT0M (), collision.posZ (), triggSign, assocSign);
397- }
398- }
399- }
370+
371+ if (trigger.pdgCode () > 0 ) {
372+ triggSign = 1 ;
373+ rQARegistry.fill (HIST (" hNLambdas_MC" ), 1 , trigger.pt (), collision.centFT0M ());
374+ } else if (trigger.pdgCode () < 0 ) {
375+ triggSign = -1 ;
376+ rQARegistry.fill (HIST (" hNLambdas_MC" ), -1 , trigger.pt (), collision.centFT0M ());
377+ }
378+
379+ for (const auto & associate : mcAssociates) {
380+ if (associate.isPhysicalPrimary ()) {
381+
382+ if (associate.pdgCode () > 0 ) {
383+ assocSign = 1 ;
384+ } else if (associate.pdgCode () < 0 ) {
385+ assocSign = -1 ;
386+ }
387+
388+ deltaPhi = RecoDecay::constrainAngle (trigger.phi () - associate.phi (), -constants::math::PIHalf);
389+ deltaEta = trigger.eta () - associate.eta ();
390+
391+ if (std::abs (associate.pdgCode ()) == kPiPlus ) {
392+ rSECorrRegistry.fill (HIST (" hSameLambdaPion_MC" ), deltaPhi, deltaEta, collision.centFT0M (), collision.posZ (), triggSign, assocSign);
393+ } else if (std::abs (associate.pdgCode ()) == kKPlus ) {
394+ rSECorrRegistry.fill (HIST (" hSameLambdaKaon_MC" ), deltaPhi, deltaEta, collision.centFT0M (), collision.posZ (), triggSign, assocSign);
395+ } else if (std::abs (associate.pdgCode ()) == kProton ) {
396+ rSECorrRegistry.fill (HIST (" hSameLambdaProton_MC" ), deltaPhi, deltaEta, collision.centFT0M (), collision.posZ (), triggSign, assocSign);
397+ }
398+ }
399+ }
400400 }
401401 }
402402 // End of the MC Same-Event Correlations
@@ -407,49 +407,49 @@ struct ThreeParticleCorrelations {
407407
408408 // Start of the MC Mixed-events Correlations
409409 for (const auto & [coll_1, v0_1, coll_2, track_2] : pairMC) {
410- 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 ;
410+ 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 ;
411411 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 ) ||
412- ((aod::mcparticle::pdgCode == static_cast <int >(kKPlus ) || aod::mcparticle::pdgCode == static_cast <int >(kKMinus )) && aod::mcparticle::pt > 0 .5f && aod::mcparticle::pt < 2 .5f ) ||
413- ((aod::mcparticle::pdgCode == static_cast <int >(kProton ) || aod::mcparticle::pdgCode == static_cast <int >(kProtonBar )) && aod::mcparticle::pt > 0 .5f ));
412+ ((aod::mcparticle::pdgCode == static_cast <int >(kKPlus ) || aod::mcparticle::pdgCode == static_cast <int >(kKMinus )) && aod::mcparticle::pt > 0 .5f && aod::mcparticle::pt < 2 .5f ) ||
413+ ((aod::mcparticle::pdgCode == static_cast <int >(kProton ) || aod::mcparticle::pdgCode == static_cast <int >(kProtonBar )) && aod::mcparticle::pt > 0 .5f ));
414414 mcTriggers.bindTable (v0_1);
415415 mcAssociates.bindTable (track_2);
416-
416+
417417 for (const auto & [trigger, associate] : soa::combinations (soa::CombinationsFullIndexPolicy (mcTriggers, mcAssociates))) {
418- if (trigger.isPhysicalPrimary () && associate.isPhysicalPrimary ()) {
419-
420- if (trigger.pdgCode () > 0 ) {
421- triggSign = 1 ;
422- } else if (trigger.pdgCode () < 0 ) {
423- triggSign = -1 ;
424- }
425- if (associate.pdgCode () > 0 ) {
426- assocSign = 1 ;
427- } else if (associate.pdgCode () < 0 ) {
428- assocSign = -1 ;
429- }
430-
431- deltaPhi = RecoDecay::constrainAngle (trigger.phi () - associate.phi (), -constants::math::PIHalf);
432- deltaEta = trigger.eta () - associate.eta ();
433-
434- if (std::abs (associate.pdgCode ()) == kPiPlus ) {
435- rMECorrRegistry.fill (HIST (" hMixLambdaPion_MC" ), deltaPhi, deltaEta, coll_1.centFT0M (), coll_1.posZ (), triggSign, assocSign);
436- } else if (std::abs (associate.pdgCode ()) == kKPlus ) {
437- rMECorrRegistry.fill (HIST (" hMixLambdaKaon_MC" ), deltaPhi, deltaEta, coll_1.centFT0M (), coll_1.posZ (), triggSign, assocSign);
438- } else if (std::abs (associate.pdgCode ()) == kProton ) {
439- rMECorrRegistry.fill (HIST (" hMixLambdaProton_MC" ), deltaPhi, deltaEta, coll_1.centFT0M (), coll_1.posZ (), triggSign, assocSign);
440- }
441- }
418+ if (trigger.isPhysicalPrimary () && associate.isPhysicalPrimary ()) {
419+
420+ if (trigger.pdgCode () > 0 ) {
421+ triggSign = 1 ;
422+ } else if (trigger.pdgCode () < 0 ) {
423+ triggSign = -1 ;
424+ }
425+ if (associate.pdgCode () > 0 ) {
426+ assocSign = 1 ;
427+ } else if (associate.pdgCode () < 0 ) {
428+ assocSign = -1 ;
429+ }
430+
431+ deltaPhi = RecoDecay::constrainAngle (trigger.phi () - associate.phi (), -constants::math::PIHalf);
432+ deltaEta = trigger.eta () - associate.eta ();
433+
434+ if (std::abs (associate.pdgCode ()) == kPiPlus ) {
435+ rMECorrRegistry.fill (HIST (" hMixLambdaPion_MC" ), deltaPhi, deltaEta, coll_1.centFT0M (), coll_1.posZ (), triggSign, assocSign);
436+ } else if (std::abs (associate.pdgCode ()) == kKPlus ) {
437+ rMECorrRegistry.fill (HIST (" hMixLambdaKaon_MC" ), deltaPhi, deltaEta, coll_1.centFT0M (), coll_1.posZ (), triggSign, assocSign);
438+ } else if (std::abs (associate.pdgCode ()) == kProton ) {
439+ rMECorrRegistry.fill (HIST (" hMixLambdaProton_MC" ), deltaPhi, deltaEta, coll_1.centFT0M (), coll_1.posZ (), triggSign, assocSign);
440+ }
441+ }
442442 }
443443 }
444444 // End of the MC Mixed-events Correlations
445- }
445+ }
446446
447447 void processMCGen (MyFilteredMCGenCollision const &, MyFilteredMCParticles const & particles)
448448 {
449449
450450 Partition<MyFilteredMCParticles> mcParticles = aod::mcparticle::pt > 0 .2f && aod::mcparticle::pt < 3 .0f ;
451451 mcParticles.bindTable (particles);
452-
452+
453453 // Start of the Monte-Carlo generated QA
454454 for (const auto & particle : mcParticles) {
455455 if (particle.isPhysicalPrimary ()) {
0 commit comments