2828
2929#include " TPDGCode.h"
3030
31- #include < string>
3231#include < algorithm>
32+ #include < string>
3333#include < vector>
3434
3535using namespace o2 ;
@@ -443,7 +443,7 @@ struct ThreeParticleCorrelations {
443443 deltaPhi = RecoDecay::constrainAngle (trigger.phi () - associate.phi (), -constants::math::PIHalf);
444444 deltaEta = trigger.eta () - associate.eta ();
445445
446- if (candMass >= MassLambda0 - invMassNSigma* dGaussSigma && candMass <= MassLambda0 + invMassNSigma* dGaussSigma) {
446+ if (candMass >= MassLambda0 - invMassNSigma * dGaussSigma && candMass <= MassLambda0 + invMassNSigma * dGaussSigma) {
447447 if (assocPID[0 ] == pionID) { // Pions
448448 rSECorrRegistry.fill (HIST (" hSameLambdaPion_SGNL" ), deltaPhi, deltaEta, collision.centFT0C (), collision.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionPions, associate, collision.centFT0C ()) / (trackEff (hEffPions, associate, collision.centFT0C ()) * v0Efficiency));
449449 } else if (assocPID[0 ] == kaonID) { // Kaons
@@ -452,28 +452,28 @@ struct ThreeParticleCorrelations {
452452 rSECorrRegistry.fill (HIST (" hSameLambdaProton_SGNL" ), deltaPhi, deltaEta, collision.centFT0C (), collision.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionProtons, associate, collision.centFT0C ()) / (trackEff (hEffProtons, associate, collision.centFT0C ()) * v0Efficiency));
453453 }
454454
455- } else if (candMass >= MassLambda0 - 2 * invMassNSigma* dGaussSigma && candMass <= MassLambda0 + 2 * invMassNSigma* dGaussSigma) {
455+ } else if (candMass >= MassLambda0 - 2 * invMassNSigma * dGaussSigma && candMass <= MassLambda0 + 2 * invMassNSigma * dGaussSigma) {
456456 if (assocPID[0 ] == pionID) { // Pions
457457 rSECorrRegistry.fill (HIST (" hSameLambdaPion_SB" ), deltaPhi, deltaEta, collision.centFT0C (), collision.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionPions, associate, collision.centFT0C ()) / (trackEff (hEffPions, associate, collision.centFT0C ()) * v0Efficiency));
458- if (candMass >= MassLambda0 - 2 * invMassNSigma* dGaussSigma && candMass < MassLambda0 - invMassNSigma* dGaussSigma) {
458+ if (candMass >= MassLambda0 - 2 * invMassNSigma * dGaussSigma && candMass < MassLambda0 - invMassNSigma * dGaussSigma) {
459459 rSECorrRegistry.fill (HIST (" hSameLambdaPion_leftSB" ), deltaPhi, deltaEta, collision.centFT0C (), collision.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionPions, associate, collision.centFT0C ()) / (trackEff (hEffPions, associate, collision.centFT0C ()) * v0Efficiency));
460- } else if (candMass > MassLambda0 + invMassNSigma* dGaussSigma && candMass <= MassLambda0 + 2 * invMassNSigma* dGaussSigma) {
460+ } else if (candMass > MassLambda0 + invMassNSigma * dGaussSigma && candMass <= MassLambda0 + 2 * invMassNSigma * dGaussSigma) {
461461 rSECorrRegistry.fill (HIST (" hSameLambdaPion_rightSB" ), deltaPhi, deltaEta, collision.centFT0C (), collision.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionPions, associate, collision.centFT0C ()) / (trackEff (hEffPions, associate, collision.centFT0C ()) * v0Efficiency));
462462 }
463463
464464 } else if (assocPID[0 ] == kaonID) { // Kaons
465465 rSECorrRegistry.fill (HIST (" hSameLambdaKaon_SB" ), deltaPhi, deltaEta, collision.centFT0C (), collision.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionKaons, associate, collision.centFT0C ()) / (trackEff (hEffKaons, associate, collision.centFT0C ()) * v0Efficiency));
466- if (candMass >= MassLambda0 - 2 * invMassNSigma* dGaussSigma && candMass < MassLambda0 - invMassNSigma* dGaussSigma) {
466+ if (candMass >= MassLambda0 - 2 * invMassNSigma * dGaussSigma && candMass < MassLambda0 - invMassNSigma * dGaussSigma) {
467467 rSECorrRegistry.fill (HIST (" hSameLambdaKaon_leftSB" ), deltaPhi, deltaEta, collision.centFT0C (), collision.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionKaons, associate, collision.centFT0C ()) / (trackEff (hEffKaons, associate, collision.centFT0C ()) * v0Efficiency));
468- } else if (candMass > MassLambda0 + invMassNSigma* dGaussSigma && candMass <= MassLambda0 + 2 * invMassNSigma* dGaussSigma) {
468+ } else if (candMass > MassLambda0 + invMassNSigma * dGaussSigma && candMass <= MassLambda0 + 2 * invMassNSigma * dGaussSigma) {
469469 rSECorrRegistry.fill (HIST (" hSameLambdaKaon_rightSB" ), deltaPhi, deltaEta, collision.centFT0C (), collision.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionKaons, associate, collision.centFT0C ()) / (trackEff (hEffKaons, associate, collision.centFT0C ()) * v0Efficiency));
470470 }
471471
472472 } else if (assocPID[0 ] == protonID) { // Protons
473473 rSECorrRegistry.fill (HIST (" hSameLambdaProton_SB" ), deltaPhi, deltaEta, collision.centFT0C (), collision.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionProtons, associate, collision.centFT0C ()) / (trackEff (hEffProtons, associate, collision.centFT0C ()) * v0Efficiency));
474- if (candMass >= MassLambda0 - 2 * invMassNSigma* dGaussSigma && candMass < MassLambda0 - invMassNSigma* dGaussSigma) {
474+ if (candMass >= MassLambda0 - 2 * invMassNSigma * dGaussSigma && candMass < MassLambda0 - invMassNSigma * dGaussSigma) {
475475 rSECorrRegistry.fill (HIST (" hSameLambdaProton_leftSB" ), deltaPhi, deltaEta, collision.centFT0C (), collision.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionProtons, associate, collision.centFT0C ()) / (trackEff (hEffProtons, associate, collision.centFT0C ()) * v0Efficiency));
476- } else if (candMass > MassLambda0 + invMassNSigma* dGaussSigma && candMass <= MassLambda0 + 2 * invMassNSigma* dGaussSigma) {
476+ } else if (candMass > MassLambda0 + invMassNSigma * dGaussSigma && candMass <= MassLambda0 + 2 * invMassNSigma * dGaussSigma) {
477477 rSECorrRegistry.fill (HIST (" hSameLambdaProton_rightSB" ), deltaPhi, deltaEta, collision.centFT0C (), collision.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionProtons, associate, collision.centFT0C ()) / (trackEff (hEffProtons, associate, collision.centFT0C ()) * v0Efficiency));
478478 }
479479 }
@@ -520,7 +520,7 @@ struct ThreeParticleCorrelations {
520520 deltaPhi = RecoDecay::constrainAngle (trigger.phi () - associate.phi (), -constants::math::PIHalf);
521521 deltaEta = trigger.eta () - associate.eta ();
522522
523- if (candMass >= MassLambda0 - invMassNSigma* dGaussSigma && candMass <= MassLambda0 + invMassNSigma* dGaussSigma) {
523+ if (candMass >= MassLambda0 - invMassNSigma * dGaussSigma && candMass <= MassLambda0 + invMassNSigma * dGaussSigma) {
524524 if (assocPID[0 ] == pionID) { // Pions
525525 rMECorrRegistry.fill (HIST (" hMixLambdaPion_SGNL" ), deltaPhi, deltaEta, coll_1.centFT0C (), coll_1.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionPions, associate, coll_1.centFT0C ()) / (trackEff (hEffPions, associate, coll_1.centFT0C ()) * v0Efficiency));
526526 } else if (assocPID[0 ] == kaonID) { // Kaons
@@ -529,28 +529,28 @@ struct ThreeParticleCorrelations {
529529 rMECorrRegistry.fill (HIST (" hMixLambdaProton_SGNL" ), deltaPhi, deltaEta, coll_1.centFT0C (), coll_1.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionProtons, associate, coll_1.centFT0C ()) / (trackEff (hEffProtons, associate, coll_1.centFT0C ()) * v0Efficiency));
530530 }
531531
532- } else if (candMass >= MassLambda0 - 2 * invMassNSigma* dGaussSigma && candMass <= MassLambda0 + 2 * invMassNSigma* dGaussSigma) {
532+ } else if (candMass >= MassLambda0 - 2 * invMassNSigma * dGaussSigma && candMass <= MassLambda0 + 2 * invMassNSigma * dGaussSigma) {
533533 if (assocPID[0 ] == pionID) { // Pions
534534 rMECorrRegistry.fill (HIST (" hMixLambdaPion_SB" ), deltaPhi, deltaEta, coll_1.centFT0C (), coll_1.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionPions, associate, coll_1.centFT0C ()) / (trackEff (hEffPions, associate, coll_1.centFT0C ()) * v0Efficiency));
535- if (candMass >= MassLambda0 - 2 * invMassNSigma* dGaussSigma && candMass < MassLambda0 - invMassNSigma* dGaussSigma) {
535+ if (candMass >= MassLambda0 - 2 * invMassNSigma * dGaussSigma && candMass < MassLambda0 - invMassNSigma * dGaussSigma) {
536536 rMECorrRegistry.fill (HIST (" hMixLambdaPion_leftSB" ), deltaPhi, deltaEta, coll_1.centFT0C (), coll_1.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionPions, associate, coll_1.centFT0C ()) / (trackEff (hEffPions, associate, coll_1.centFT0C ()) * v0Efficiency));
537- } else if (candMass > MassLambda0 + invMassNSigma* dGaussSigma && candMass <= MassLambda0 + 2 * invMassNSigma* dGaussSigma) {
537+ } else if (candMass > MassLambda0 + invMassNSigma * dGaussSigma && candMass <= MassLambda0 + 2 * invMassNSigma * dGaussSigma) {
538538 rMECorrRegistry.fill (HIST (" hMixLambdaPion_rightSB" ), deltaPhi, deltaEta, coll_1.centFT0C (), coll_1.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionPions, associate, coll_1.centFT0C ()) / (trackEff (hEffPions, associate, coll_1.centFT0C ()) * v0Efficiency));
539539 }
540540
541541 } else if (assocPID[0 ] == kaonID) { // Kaons
542542 rMECorrRegistry.fill (HIST (" hMixLambdaKaon_SB" ), deltaPhi, deltaEta, coll_1.centFT0C (), coll_1.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionKaons, associate, coll_1.centFT0C ()) / (trackEff (hEffKaons, associate, coll_1.centFT0C ()) * v0Efficiency));
543- if (candMass >= MassLambda0 - 2 * invMassNSigma* dGaussSigma && candMass < MassLambda0 - invMassNSigma* dGaussSigma) {
543+ if (candMass >= MassLambda0 - 2 * invMassNSigma * dGaussSigma && candMass < MassLambda0 - invMassNSigma * dGaussSigma) {
544544 rMECorrRegistry.fill (HIST (" hMixLambdaKaon_leftSB" ), deltaPhi, deltaEta, coll_1.centFT0C (), coll_1.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionKaons, associate, coll_1.centFT0C ()) / (trackEff (hEffKaons, associate, coll_1.centFT0C ()) * v0Efficiency));
545- } else if (candMass > MassLambda0 + invMassNSigma* dGaussSigma && candMass <= MassLambda0 + 2 * invMassNSigma* dGaussSigma) {
545+ } else if (candMass > MassLambda0 + invMassNSigma * dGaussSigma && candMass <= MassLambda0 + 2 * invMassNSigma * dGaussSigma) {
546546 rMECorrRegistry.fill (HIST (" hMixLambdaKaon_rightSB" ), deltaPhi, deltaEta, coll_1.centFT0C (), coll_1.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionKaons, associate, coll_1.centFT0C ()) / (trackEff (hEffKaons, associate, coll_1.centFT0C ()) * v0Efficiency));
547547 }
548548
549549 } else if (assocPID[0 ] == protonID) { // Protons
550550 rMECorrRegistry.fill (HIST (" hMixLambdaProton_SB" ), deltaPhi, deltaEta, coll_1.centFT0C (), coll_1.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionProtons, associate, coll_1.centFT0C ()) / (trackEff (hEffProtons, associate, coll_1.centFT0C ()) * v0Efficiency));
551- if (candMass >= MassLambda0 - 2 * invMassNSigma* dGaussSigma && candMass < MassLambda0 - invMassNSigma* dGaussSigma) {
551+ if (candMass >= MassLambda0 - 2 * invMassNSigma * dGaussSigma && candMass < MassLambda0 - invMassNSigma * dGaussSigma) {
552552 rMECorrRegistry.fill (HIST (" hMixLambdaProton_leftSB" ), deltaPhi, deltaEta, coll_1.centFT0C (), coll_1.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionProtons, associate, coll_1.centFT0C ()) / (trackEff (hEffProtons, associate, coll_1.centFT0C ()) * v0Efficiency));
553- } else if (candMass > MassLambda0 + invMassNSigma* dGaussSigma && candMass <= MassLambda0 + 2 * invMassNSigma* dGaussSigma) {
553+ } else if (candMass > MassLambda0 + invMassNSigma * dGaussSigma && candMass <= MassLambda0 + 2 * invMassNSigma * dGaussSigma) {
554554 rMECorrRegistry.fill (HIST (" hMixLambdaProton_rightSB" ), deltaPhi, deltaEta, coll_1.centFT0C (), coll_1.posZ (), triggSign, associate.sign (), ratioCorrection (hCorrectionProtons, associate, coll_1.centFT0C ()) / (trackEff (hEffProtons, associate, coll_1.centFT0C ()) * v0Efficiency));
555555 }
556556 }
@@ -1175,7 +1175,7 @@ struct ThreeParticleCorrelations {
11751175 }
11761176
11771177 double invMass = RecoDecay::m (std::array{dMomArray, aMomArray}, massArray);
1178- if (invMass >= MassLambda0 - invMassNSigma* dGaussSigma && invMass <= MassLambda0 + invMassNSigma* dGaussSigma) {
1178+ if (invMass >= MassLambda0 - invMassNSigma * dGaussSigma && invMass <= MassLambda0 + invMassNSigma * dGaussSigma) {
11791179 return false ;
11801180 }
11811181 }
0 commit comments