@@ -25,7 +25,7 @@ using namespace o2::framework::expressions;
2525
2626struct ThreePartCorr {
2727 Service<o2::ccdb::BasicCCDBManager> CCDB;
28-
28+
2929 // Histogram registry
3030 HistogramRegistry MECorrRegistry{" MECorrRegistry" , {}, OutputObjHandlingPolicy::AnalysisObject, false , true };
3131 HistogramRegistry SECorrRegistry{" SECorrRegistry" , {}, OutputObjHandlingPolicy::AnalysisObject, false , true };
@@ -84,7 +84,7 @@ struct ThreePartCorr {
8484 // Efficiency histograms
8585 TH1D** hEffPions = new TH1D*[2 ];
8686 TH1D** hEffKaons = new TH1D*[2 ];
87- TH1D** hEffProtons = new TH1D*[2 ];
87+ TH1D** hEffProtons = new TH1D*[2 ];
8888
8989 // Correlation variables
9090 Int_t T_Sign;
@@ -97,7 +97,7 @@ struct ThreePartCorr {
9797
9898 void init (InitContext const &)
9999 {
100-
100+
101101 const AxisSpec CentralityAxis{ConfCentBins};
102102 const AxisSpec ZvtxAxis{ConfZvtxBins};
103103 const AxisSpec PhiAxis{36 , (-1 . / 2 ) * M_PI, (3 . / 2 ) * M_PI};
@@ -162,7 +162,7 @@ struct ThreePartCorr {
162162 hEffKaons[0 ] = static_cast <TH1D*>(EfficiencyList->FindObject (" hEfficiencyKaonP" ));
163163 hEffKaons[1 ] = static_cast <TH1D*>(EfficiencyList->FindObject (" hEfficiencyKaonN" ));
164164 hEffProtons[0 ] = static_cast <TH1D*>(EfficiencyList->FindObject (" hEfficiencyProtonP" ));
165- hEffProtons[1 ] = static_cast <TH1D*>(EfficiencyList->FindObject (" hEfficiencyProtonN" ));
165+ hEffProtons[1 ] = static_cast <TH1D*>(EfficiencyList->FindObject (" hEfficiencyProtonN" ));
166166 }
167167
168168 // ================================================================================================================================================================================================================
@@ -222,19 +222,19 @@ struct ThreePartCorr {
222222
223223 if (CandMass >= massLambda - 4 * DGaussSigma && CandMass <= massLambda + 4 * DGaussSigma) {
224224 if (A_PID[0 ] == 0.0 ) { // Pions
225- SECorrRegistry.fill (HIST (" hSameLambdaPion_SGNL" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffPions, associate.sign (), associate.pt ()));
225+ SECorrRegistry.fill (HIST (" hSameLambdaPion_SGNL" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffPions, associate.sign (), associate.pt ()));
226226 } else if (A_PID[0 ] == 1.0 ) { // Kaons
227- SECorrRegistry.fill (HIST (" hSameLambdaKaon_SGNL" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffKaons, associate.sign (), associate.pt ()));
227+ SECorrRegistry.fill (HIST (" hSameLambdaKaon_SGNL" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffKaons, associate.sign (), associate.pt ()));
228228 } else if (A_PID[0 ] == 2.0 ) { // Protons
229- SECorrRegistry.fill (HIST (" hSameLambdaProton_SGNL" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffProtons, associate.sign (), associate.pt ()));
229+ SECorrRegistry.fill (HIST (" hSameLambdaProton_SGNL" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffProtons, associate.sign (), associate.pt ()));
230230 }
231231 } else if (CandMass >= massLambda - 8 * DGaussSigma && CandMass <= massLambda + 8 * DGaussSigma) {
232232 if (A_PID[0 ] == 0.0 ) { // Pions
233- SECorrRegistry.fill (HIST (" hSameLambdaPion_SB" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffPions, associate.sign (), associate.pt ()));
233+ SECorrRegistry.fill (HIST (" hSameLambdaPion_SB" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffPions, associate.sign (), associate.pt ()));
234234 } else if (A_PID[0 ] == 1.0 ) { // Kaons
235- SECorrRegistry.fill (HIST (" hSameLambdaKaon_SB" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffKaons, associate.sign (), associate.pt ()));
235+ SECorrRegistry.fill (HIST (" hSameLambdaKaon_SB" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffKaons, associate.sign (), associate.pt ()));
236236 } else if (A_PID[0 ] == 2.0 ) { // Protons
237- SECorrRegistry.fill (HIST (" hSameLambdaProton_SB" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffProtons, associate.sign (), associate.pt ()));
237+ SECorrRegistry.fill (HIST (" hSameLambdaProton_SB" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffProtons, associate.sign (), associate.pt ()));
238238 }
239239 }
240240 }
@@ -267,19 +267,19 @@ struct ThreePartCorr {
267267
268268 if (CandMass >= massLambda - 4 * DGaussSigma && CandMass <= massLambda + 4 * DGaussSigma) {
269269 if (A_PID[0 ] == 0.0 ) { // Pions
270- MECorrRegistry.fill (HIST (" hMixLambdaPion_SGNL" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffPions, associate.sign (), associate.pt ()));
270+ MECorrRegistry.fill (HIST (" hMixLambdaPion_SGNL" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffPions, associate.sign (), associate.pt ()));
271271 } else if (A_PID[0 ] == 1.0 ) { // Kaons
272- MECorrRegistry.fill (HIST (" hMixLambdaKaon_SGNL" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffKaons, associate.sign (), associate.pt ()));
272+ MECorrRegistry.fill (HIST (" hMixLambdaKaon_SGNL" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffKaons, associate.sign (), associate.pt ()));
273273 } else if (A_PID[0 ] == 2.0 ) { // Protons
274- MECorrRegistry.fill (HIST (" hMixLambdaProton_SGNL" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffProtons, associate.sign (), associate.pt ()));
274+ MECorrRegistry.fill (HIST (" hMixLambdaProton_SGNL" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffProtons, associate.sign (), associate.pt ()));
275275 }
276276 } else if (CandMass >= massLambda - 8 * DGaussSigma && CandMass <= massLambda + 8 * DGaussSigma) {
277277 if (A_PID[0 ] == 0.0 ) { // Pions
278- MECorrRegistry.fill (HIST (" hMixLambdaPion_SB" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffPions, associate.sign (), associate.pt ()));
278+ MECorrRegistry.fill (HIST (" hMixLambdaPion_SB" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffPions, associate.sign (), associate.pt ()));
279279 } else if (A_PID[0 ] == 1.0 ) { // Kaons
280- MECorrRegistry.fill (HIST (" hMixLambdaKaon_SB" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffKaons, associate.sign (), associate.pt ()));
280+ MECorrRegistry.fill (HIST (" hMixLambdaKaon_SB" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffKaons, associate.sign (), associate.pt ()));
281281 } else if (A_PID[0 ] == 2.0 ) { // Protons
282- MECorrRegistry.fill (HIST (" hMixLambdaProton_SB" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffProtons, associate.sign (), associate.pt ()));
282+ MECorrRegistry.fill (HIST (" hMixLambdaProton_SB" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffProtons, associate.sign (), associate.pt ()));
283283 }
284284 }
285285 }
@@ -394,7 +394,7 @@ struct ThreePartCorr {
394394
395395 return 0 ;
396396 }
397-
397+
398398 template <class TrackCand >
399399 Double_t* TrackPID (const TrackCand& Track)
400400 {
@@ -477,25 +477,25 @@ struct ThreePartCorr {
477477 std::array<float , 3 > DMomArray;
478478 std::array<float , 3 > AMomArray = Track.pVector ();
479479 if (TrackPID (Track)[0 ] == 0.0 ) {
480- MassArray = {o2::constants::physics::MassProton, o2::constants::physics::MassPionCharged};
481-
482- if (V0Sign (V0) == 1 && Track.sign () == -1 ) { // Lambda - Pi_min
483- const auto & dTrack = V0.template posTrack_as <MyFilteredTracks>();
484- DMomArray = dTrack.pVector ();
485- } else if (V0Sign (V0) == -1 && Track.sign () == 1 ) { // Antilambda - Pi_plus
486- const auto & dTrack = V0.template negTrack_as <MyFilteredTracks>();
487- DMomArray = dTrack.pVector ();
488- }
489- } else if (TrackPID (Track)[0 ] == 2.0 ) {
490- MassArray = {o2::constants::physics::MassPionCharged, o2::constants::physics::MassProton};
491-
492- if (V0Sign (V0) == 1 && Track.sign () == 1 ) { // Lambda - Proton
493- const auto & dTrack = V0.template negTrack_as <MyFilteredTracks>();
494- DMomArray = dTrack.pVector ();
495- } else if (V0Sign (V0) == -1 && Track.sign () == -1 ) { // Antilambda - Antiproton
496- const auto & dTrack = V0.template posTrack_as <MyFilteredTracks>();
497- DMomArray = dTrack.pVector ();
498- }
480+ MassArray = {o2::constants::physics::MassProton, o2::constants::physics::MassPionCharged};
481+
482+ if (V0Sign (V0) == 1 && Track.sign () == -1 ) { // Lambda - Pi_min
483+ const auto & dTrack = V0.template posTrack_as <MyFilteredTracks>();
484+ DMomArray = dTrack.pVector ();
485+ } else if (V0Sign (V0) == -1 && Track.sign () == 1 ) { // Antilambda - Pi_plus
486+ const auto & dTrack = V0.template negTrack_as <MyFilteredTracks>();
487+ DMomArray = dTrack.pVector ();
488+ }
489+ } else if (TrackPID (Track)[0 ] == 2.0 ) {
490+ MassArray = {o2::constants::physics::MassPionCharged, o2::constants::physics::MassProton};
491+
492+ if (V0Sign (V0) == 1 && Track.sign () == 1 ) { // Lambda - Proton
493+ const auto & dTrack = V0.template negTrack_as <MyFilteredTracks>();
494+ DMomArray = dTrack.pVector ();
495+ } else if (V0Sign (V0) == -1 && Track.sign () == -1 ) { // Antilambda - Antiproton
496+ const auto & dTrack = V0.template posTrack_as <MyFilteredTracks>();
497+ DMomArray = dTrack.pVector ();
498+ }
499499 }
500500
501501 Double_t M = RecoDecay::m (std::array{DMomArray, AMomArray}, MassArray);
0 commit comments