@@ -27,7 +27,7 @@ using namespace o2::framework::expressions;
2727
2828struct ThreePartCorr {
2929 Service<o2::ccdb::BasicCCDBManager> CCDB;
30-
30+
3131 // Histogram registry
3232 HistogramRegistry MECorrRegistry{" MECorrRegistry" , {}, OutputObjHandlingPolicy::AnalysisObject, false , true };
3333 HistogramRegistry SECorrRegistry{" SECorrRegistry" , {}, OutputObjHandlingPolicy::AnalysisObject, false , true };
@@ -86,7 +86,7 @@ struct ThreePartCorr {
8686 // Efficiency histograms
8787 TH1D** hEffPions = new TH1D*[2 ];
8888 TH1D** hEffKaons = new TH1D*[2 ];
89- TH1D** hEffProtons = new TH1D*[2 ];
89+ TH1D** hEffProtons = new TH1D*[2 ];
9090
9191 // Correlation variables
9292 Int_t T_Sign;
@@ -99,7 +99,7 @@ struct ThreePartCorr {
9999
100100 void init (InitContext const &)
101101 {
102-
102+
103103 const AxisSpec CentralityAxis{ConfCentBins};
104104 const AxisSpec ZvtxAxis{ConfZvtxBins};
105105 const AxisSpec PhiAxis{36 , (-1 . / 2 ) * M_PI, (3 . / 2 ) * M_PI};
@@ -164,7 +164,7 @@ struct ThreePartCorr {
164164 hEffKaons[0 ] = static_cast <TH1D*>(EfficiencyList->FindObject (" hEfficiencyKaonP" ));
165165 hEffKaons[1 ] = static_cast <TH1D*>(EfficiencyList->FindObject (" hEfficiencyKaonN" ));
166166 hEffProtons[0 ] = static_cast <TH1D*>(EfficiencyList->FindObject (" hEfficiencyProtonP" ));
167- hEffProtons[1 ] = static_cast <TH1D*>(EfficiencyList->FindObject (" hEfficiencyProtonN" ));
167+ hEffProtons[1 ] = static_cast <TH1D*>(EfficiencyList->FindObject (" hEfficiencyProtonN" ));
168168 }
169169
170170 // ================================================================================================================================================================================================================
@@ -224,19 +224,19 @@ struct ThreePartCorr {
224224
225225 if (CandMass >= massLambda - 4 * DGaussSigma && CandMass <= massLambda + 4 * DGaussSigma) {
226226 if (A_PID[0 ] == 0.0 ) { // Pions
227- SECorrRegistry.fill (HIST (" hSameLambdaPion_SGNL" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffPions, associate.sign (), associate.pt ()));
227+ SECorrRegistry.fill (HIST (" hSameLambdaPion_SGNL" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffPions, associate.sign (), associate.pt ()));
228228 } else if (A_PID[0 ] == 1.0 ) { // Kaons
229- SECorrRegistry.fill (HIST (" hSameLambdaKaon_SGNL" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffKaons, associate.sign (), associate.pt ()));
229+ SECorrRegistry.fill (HIST (" hSameLambdaKaon_SGNL" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffKaons, associate.sign (), associate.pt ()));
230230 } else if (A_PID[0 ] == 2.0 ) { // Protons
231- SECorrRegistry.fill (HIST (" hSameLambdaProton_SGNL" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffProtons, associate.sign (), associate.pt ()));
231+ SECorrRegistry.fill (HIST (" hSameLambdaProton_SGNL" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffProtons, associate.sign (), associate.pt ()));
232232 }
233233 } else if (CandMass >= massLambda - 8 * DGaussSigma && CandMass <= massLambda + 8 * DGaussSigma) {
234234 if (A_PID[0 ] == 0.0 ) { // Pions
235- SECorrRegistry.fill (HIST (" hSameLambdaPion_SB" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffPions, associate.sign (), associate.pt ()));
235+ SECorrRegistry.fill (HIST (" hSameLambdaPion_SB" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffPions, associate.sign (), associate.pt ()));
236236 } else if (A_PID[0 ] == 1.0 ) { // Kaons
237- SECorrRegistry.fill (HIST (" hSameLambdaKaon_SB" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffKaons, associate.sign (), associate.pt ()));
237+ SECorrRegistry.fill (HIST (" hSameLambdaKaon_SB" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffKaons, associate.sign (), associate.pt ()));
238238 } else if (A_PID[0 ] == 2.0 ) { // Protons
239- SECorrRegistry.fill (HIST (" hSameLambdaProton_SB" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffProtons, associate.sign (), associate.pt ()));
239+ SECorrRegistry.fill (HIST (" hSameLambdaProton_SB" ), DeltaPhi, DeltaEta, collision.centFT0C (), collision.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffProtons, associate.sign (), associate.pt ()));
240240 }
241241 }
242242 }
@@ -269,19 +269,19 @@ struct ThreePartCorr {
269269
270270 if (CandMass >= massLambda - 4 * DGaussSigma && CandMass <= massLambda + 4 * DGaussSigma) {
271271 if (A_PID[0 ] == 0.0 ) { // Pions
272- MECorrRegistry.fill (HIST (" hMixLambdaPion_SGNL" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffPions, associate.sign (), associate.pt ()));
272+ MECorrRegistry.fill (HIST (" hMixLambdaPion_SGNL" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffPions, associate.sign (), associate.pt ()));
273273 } else if (A_PID[0 ] == 1.0 ) { // Kaons
274- MECorrRegistry.fill (HIST (" hMixLambdaKaon_SGNL" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffKaons, associate.sign (), associate.pt ()));
274+ MECorrRegistry.fill (HIST (" hMixLambdaKaon_SGNL" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffKaons, associate.sign (), associate.pt ()));
275275 } else if (A_PID[0 ] == 2.0 ) { // Protons
276- MECorrRegistry.fill (HIST (" hMixLambdaProton_SGNL" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffProtons, associate.sign (), associate.pt ()));
276+ MECorrRegistry.fill (HIST (" hMixLambdaProton_SGNL" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffProtons, associate.sign (), associate.pt ()));
277277 }
278278 } else if (CandMass >= massLambda - 8 * DGaussSigma && CandMass <= massLambda + 8 * DGaussSigma) {
279279 if (A_PID[0 ] == 0.0 ) { // Pions
280- MECorrRegistry.fill (HIST (" hMixLambdaPion_SB" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffPions, associate.sign (), associate.pt ()));
280+ MECorrRegistry.fill (HIST (" hMixLambdaPion_SB" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffPions, associate.sign (), associate.pt ()));
281281 } else if (A_PID[0 ] == 1.0 ) { // Kaons
282- MECorrRegistry.fill (HIST (" hMixLambdaKaon_SB" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffKaons, associate.sign (), associate.pt ()));
282+ MECorrRegistry.fill (HIST (" hMixLambdaKaon_SB" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffKaons, associate.sign (), associate.pt ()));
283283 } else if (A_PID[0 ] == 2.0 ) { // Protons
284- MECorrRegistry.fill (HIST (" hMixLambdaProton_SB" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 ./ TrackEff (hEffProtons, associate.sign (), associate.pt ()));
284+ MECorrRegistry.fill (HIST (" hMixLambdaProton_SB" ), DeltaPhi, DeltaEta, coll_1.centFT0C (), coll_1.posZ (), T_Sign, associate.sign (), 1 . / TrackEff (hEffProtons, associate.sign (), associate.pt ()));
285285 }
286286 }
287287 }
@@ -396,7 +396,7 @@ struct ThreePartCorr {
396396
397397 return 0 ;
398398 }
399-
399+
400400 template <class TrackCand >
401401 Double_t* TrackPID (const TrackCand& Track)
402402 {
@@ -479,25 +479,25 @@ struct ThreePartCorr {
479479 std::array<float , 3 > DMomArray;
480480 std::array<float , 3 > AMomArray = Track.pVector ();
481481 if (TrackPID (Track)[0 ] == 0.0 ) {
482- MassArray = {o2::constants::physics::MassProton, o2::constants::physics::MassPionCharged};
483-
484- if (V0Sign (V0) == 1 && Track.sign () == -1 ) { // Lambda - Pi_min
485- const auto & dTrack = V0.template posTrack_as <MyFilteredTracks>();
486- DMomArray = dTrack.pVector ();
487- } else if (V0Sign (V0) == -1 && Track.sign () == 1 ) { // Antilambda - Pi_plus
488- const auto & dTrack = V0.template negTrack_as <MyFilteredTracks>();
489- DMomArray = dTrack.pVector ();
490- }
491- } else if (TrackPID (Track)[0 ] == 2.0 ) {
492- MassArray = {o2::constants::physics::MassPionCharged, o2::constants::physics::MassProton};
493-
494- if (V0Sign (V0) == 1 && Track.sign () == 1 ) { // Lambda - Proton
495- const auto & dTrack = V0.template negTrack_as <MyFilteredTracks>();
496- DMomArray = dTrack.pVector ();
497- } else if (V0Sign (V0) == -1 && Track.sign () == -1 ) { // Antilambda - Antiproton
498- const auto & dTrack = V0.template posTrack_as <MyFilteredTracks>();
499- DMomArray = dTrack.pVector ();
500- }
482+ MassArray = {o2::constants::physics::MassProton, o2::constants::physics::MassPionCharged};
483+
484+ if (V0Sign (V0) == 1 && Track.sign () == -1 ) { // Lambda - Pi_min
485+ const auto & dTrack = V0.template posTrack_as <MyFilteredTracks>();
486+ DMomArray = dTrack.pVector ();
487+ } else if (V0Sign (V0) == -1 && Track.sign () == 1 ) { // Antilambda - Pi_plus
488+ const auto & dTrack = V0.template negTrack_as <MyFilteredTracks>();
489+ DMomArray = dTrack.pVector ();
490+ }
491+ } else if (TrackPID (Track)[0 ] == 2.0 ) {
492+ MassArray = {o2::constants::physics::MassPionCharged, o2::constants::physics::MassProton};
493+
494+ if (V0Sign (V0) == 1 && Track.sign () == 1 ) { // Lambda - Proton
495+ const auto & dTrack = V0.template negTrack_as <MyFilteredTracks>();
496+ DMomArray = dTrack.pVector ();
497+ } else if (V0Sign (V0) == -1 && Track.sign () == -1 ) { // Antilambda - Antiproton
498+ const auto & dTrack = V0.template posTrack_as <MyFilteredTracks>();
499+ DMomArray = dTrack.pVector ();
500+ }
501501 }
502502
503503 Double_t M = RecoDecay::m (std::array{DMomArray, AMomArray}, MassArray);
0 commit comments