@@ -95,7 +95,7 @@ struct ThreeParticleCorrelations {
9595 TH1D** hEffPions = new TH1D*[2 ];
9696 TH1D** hEffKaons = new TH1D*[2 ];
9797 TH1D** hEffProtons = new TH1D*[2 ];
98-
98+
9999 // Correlation variables
100100 int triggSign;
101101 double candMass;
@@ -305,7 +305,7 @@ struct ThreeParticleCorrelations {
305305
306306 void processMixed (MyFilteredCollisions const &, MyFilteredV0s const &, MyFilteredTracks const &, aod::BCsWithTimestamps const &)
307307 {
308-
308+
309309 // Start of the Mixed-events Correlations
310310 for (const auto & [coll_1, v0_1, coll_2, track_2] : pair) {
311311
@@ -472,11 +472,10 @@ struct ThreeParticleCorrelations {
472472 }
473473 LOGF (info, " Retrieved GRP for timestamp %llu with magnetic field of %d kG" , timestamp, grpo->getNominalL3Field ());
474474 }
475-
475+
476476 return 0.1 * (grpo->getNominalL3Field ()); // 1 T = 10 kG
477477 }
478478
479-
480479 double trackEff (TH1D** efficiencies, int sign, double pT)
481480 {
482481
@@ -573,7 +572,7 @@ struct ThreeParticleCorrelations {
573572 return kFALSE ;
574573 }
575574 } else if (track.pt () > 1.5 && track.pt () < 2.3 ) {
576- if (track.tofNSigmaPi () < -4.0 ||track.tofNSigmaPi () > 0.0 ) {
575+ if (track.tofNSigmaPi () < -4.0 || track.tofNSigmaPi () > 0.0 ) {
577576 return kFALSE ;
578577 }
579578 } else if (track.pt () > 2.3 ) {
@@ -691,7 +690,7 @@ struct ThreeParticleCorrelations {
691690 dProton = v0.template negTrack_as <MyFilteredTracks>();
692691 dPion = v0.template posTrack_as <MyFilteredTracks>();
693692 }
694-
693+
695694 double dEtaProton = dProton.eta () - track.eta ();
696695 double dEtaPion = dPion.eta () - track.eta ();
697696 if (std::abs (dEtaProton) > 0.02 && std::abs (dEtaPion) > 0.02 ) {
@@ -701,41 +700,41 @@ struct ThreeParticleCorrelations {
701700 double phiStarProton, phiStarPion;
702701 double dPhiProton = dProton.phi () - track.phi ();
703702 double dPhiPion = dPion.phi () - track.phi ();
704- double phaseProton = (-0.3 *B* dProton.sign ()) / (2 * dProton.pt ());
705- double phasePion = (-0.3 *B* dPion.sign ()) / (2 * dPion.pt ());
706- double phaseTrack = (-0.3 *B* track.sign ()) / (2 * track.pt ());
703+ double phaseProton = (-0.3 * B * dProton.sign ()) / (2 * dProton.pt ());
704+ double phasePion = (-0.3 * B * dPion.sign ()) / (2 * dPion.pt ());
705+ double phaseTrack = (-0.3 * B * track.sign ()) / (2 * track.pt ());
707706
708707 for (double r = 0.8 ; r <= 2.5 ; r += 0.01 ) {
709- phiStarProton = RecoDecay::constrainAngle (dPhiProton + std::asin (phaseProton* r) - std::asin (phaseTrack* r), -constants::math::PIHalf);
710- phiStarPion = RecoDecay::constrainAngle (dPhiPion + std::asin (phasePion* r) - std::asin (phaseTrack* r), -constants::math::PIHalf);
711-
708+ phiStarProton = RecoDecay::constrainAngle (dPhiProton + std::asin (phaseProton * r) - std::asin (phaseTrack * r), -constants::math::PIHalf);
709+ phiStarPion = RecoDecay::constrainAngle (dPhiPion + std::asin (phasePion * r) - std::asin (phaseTrack * r), -constants::math::PIHalf);
710+
712711 if (r == 0.8 ) {
713- if (!Mix) {
714- rPhiStarRegistry.fill (HIST (" hSEProtonPreCut" ), phiStarProton, dEtaProton);
715- rPhiStarRegistry.fill (HIST (" hSEPionPreCut" ), phiStarPion, dEtaPion);
716- } else {
717- rPhiStarRegistry.fill (HIST (" hMEProtonPreCut" ), phiStarProton, dEtaProton);
718- rPhiStarRegistry.fill (HIST (" hMEPionPreCut" ), phiStarPion, dEtaPion);
719- }
712+ if (!Mix) {
713+ rPhiStarRegistry.fill (HIST (" hSEProtonPreCut" ), phiStarProton, dEtaProton);
714+ rPhiStarRegistry.fill (HIST (" hSEPionPreCut" ), phiStarPion, dEtaPion);
715+ } else {
716+ rPhiStarRegistry.fill (HIST (" hMEProtonPreCut" ), phiStarProton, dEtaProton);
717+ rPhiStarRegistry.fill (HIST (" hMEPionPreCut" ), phiStarPion, dEtaPion);
718+ }
720719 }
721-
720+
722721 if ((std::abs (dEtaProton) <= 0.02 && std::abs (phiStarProton) <= 0.045 ) || (std::abs (dEtaPion) <= 0.02 && std::abs (phiStarPion) <= 0.045 )) {
723- return kFALSE ;
722+ return kFALSE ;
724723 }
725-
724+
726725 if (r == 0.8 ) {
727- if (!Mix) {
728- rPhiStarRegistry.fill (HIST (" hSEProtonPostCut" ), phiStarProton, dEtaProton);
729- rPhiStarRegistry.fill (HIST (" hSEPionPostCut" ), phiStarPion, dEtaPion);
730- } else {
731- rPhiStarRegistry.fill (HIST (" hMEProtonPostCut" ), phiStarProton, dEtaProton);
732- rPhiStarRegistry.fill (HIST (" hMEPionPostCut" ), phiStarPion, dEtaPion);
733- }
726+ if (!Mix) {
727+ rPhiStarRegistry.fill (HIST (" hSEProtonPostCut" ), phiStarProton, dEtaProton);
728+ rPhiStarRegistry.fill (HIST (" hSEPionPostCut" ), phiStarPion, dEtaPion);
729+ } else {
730+ rPhiStarRegistry.fill (HIST (" hMEProtonPostCut" ), phiStarProton, dEtaProton);
731+ rPhiStarRegistry.fill (HIST (" hMEPionPostCut" ), phiStarPion, dEtaPion);
732+ }
734733 }
735734 }
736735
737736 return kTRUE ;
738- }
737+ }
739738};
740739
741740// ==================================================================================================================================================================================================================
0 commit comments