@@ -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
@@ -578,7 +577,7 @@ struct ThreeParticleCorrelations {
578577 return kFALSE ;
579578 }
580579 } else if (track.pt () > 1.5 && track.pt () < 2.3 ) {
581- if (track.tofNSigmaPi () < -4.0 ||track.tofNSigmaPi () > 0.0 ) {
580+ if (track.tofNSigmaPi () < -4.0 || track.tofNSigmaPi () > 0.0 ) {
582581 return kFALSE ;
583582 }
584583 } else if (track.pt () > 2.3 ) {
@@ -696,7 +695,7 @@ struct ThreeParticleCorrelations {
696695 const auto & dProton = v0.template negTrack_as <MyFilteredTracks>();
697696 const auto & dPion = v0.template posTrack_as <MyFilteredTracks>();
698697 }
699-
698+
700699 double dEtaProton = dProton.eta () - track.eta ();
701700 double dEtaPion = dPion.eta () - track.eta ();
702701 if (std::abs (dEtaProton) > 0.02 && std::abs (dEtaPion) > 0.02 ) {
@@ -706,41 +705,41 @@ struct ThreeParticleCorrelations {
706705 double phiStarProton, phiStarPion;
707706 double dPhiProton = dProton.phi () - track.phi ();
708707 double dPhiPion = dPion.phi () - track.phi ();
709- double phaseProton = (-0.3 *B* dProton.sign ()) / (2 * dProton.pt ());
710- double phasePion = (-0.3 *B* dPion.sign ()) / (2 * dPion.pt ());
711- double phaseTrack = (-0.3 *B* track.sign ()) / (2 * track.pt ());
708+ double phaseProton = (-0.3 * B * dProton.sign ()) / (2 * dProton.pt ());
709+ double phasePion = (-0.3 * B * dPion.sign ()) / (2 * dPion.pt ());
710+ double phaseTrack = (-0.3 * B * track.sign ()) / (2 * track.pt ());
712711
713712 for (double r = 0.8 ; r <= 2.5 ; r += 0.01 ) {
714- phiStarProton = RecoDecay::constrainAngle (dPhiProton + std::asin (phaseProton* r) - std::asin (phaseTrack* r), -constants::math::PIHalf);
715- phiStarPion = RecoDecay::constrainAngle (dPhiPion + std::asin (phasePion* r) - std::asin (phaseTrack* r), -constants::math::PIHalf);
716-
713+ phiStarProton = RecoDecay::constrainAngle (dPhiProton + std::asin (phaseProton * r) - std::asin (phaseTrack * r), -constants::math::PIHalf);
714+ phiStarPion = RecoDecay::constrainAngle (dPhiPion + std::asin (phasePion * r) - std::asin (phaseTrack * r), -constants::math::PIHalf);
715+
717716 if (r == 0.8 ) {
718- if (!Mix) {
719- rPhiStarRegistry.fill (HIST (" hSEProtonPreCut" ), phiStarProton, dEtaProton);
720- rPhiStarRegistry.fill (HIST (" hSEPionPreCut" ), phiStarPion, dEtaPion);
721- } else {
722- rPhiStarRegistry.fill (HIST (" hMEProtonPreCut" ), phiStarProton, dEtaProton);
723- rPhiStarRegistry.fill (HIST (" hMEPionPreCut" ), phiStarPion, dEtaPion);
724- }
717+ if (!Mix) {
718+ rPhiStarRegistry.fill (HIST (" hSEProtonPreCut" ), phiStarProton, dEtaProton);
719+ rPhiStarRegistry.fill (HIST (" hSEPionPreCut" ), phiStarPion, dEtaPion);
720+ } else {
721+ rPhiStarRegistry.fill (HIST (" hMEProtonPreCut" ), phiStarProton, dEtaProton);
722+ rPhiStarRegistry.fill (HIST (" hMEPionPreCut" ), phiStarPion, dEtaPion);
723+ }
725724 }
726-
725+
727726 if ((std::abs (dEtaProton) <= 0.02 && std::abs (phiStarProton) <= 0.045 ) || (std::abs (dEtaPion) <= 0.02 && std::abs (phiStarPion) <= 0.045 )) {
728- return kFALSE ;
727+ return kFALSE ;
729728 }
730-
729+
731730 if (r == 0.8 ) {
732- if (!Mix) {
733- rPhiStarRegistry.fill (HIST (" hSEProtonPostCut" ), phiStarProton, dEtaProton);
734- rPhiStarRegistry.fill (HIST (" hSEPionPostCut" ), phiStarPion, dEtaPion);
735- } else {
736- rPhiStarRegistry.fill (HIST (" hMEProtonPostCut" ), phiStarProton, dEtaProton);
737- rPhiStarRegistry.fill (HIST (" hMEPionPostCut" ), phiStarPion, dEtaPion);
738- }
731+ if (!Mix) {
732+ rPhiStarRegistry.fill (HIST (" hSEProtonPostCut" ), phiStarProton, dEtaProton);
733+ rPhiStarRegistry.fill (HIST (" hSEPionPostCut" ), phiStarPion, dEtaPion);
734+ } else {
735+ rPhiStarRegistry.fill (HIST (" hMEProtonPostCut" ), phiStarProton, dEtaProton);
736+ rPhiStarRegistry.fill (HIST (" hMEPionPostCut" ), phiStarPion, dEtaPion);
737+ }
739738 }
740739 }
741740
742741 return kTRUE ;
743- }
742+ }
744743};
745744
746745// ==================================================================================================================================================================================================================
0 commit comments