[PWGCF] FemtoUniverse -- Fix of the delta phi wrapping (AliceO2Group#9364)

Co-authored-by: Zuzanna Chochulska <[email protected]>

Author: zchochul

PWGCF/FemtoUniverse/Core/FemtoUniverseContainer.h

@@ -175,7 +175,12 @@ class FemtoUniverseContainer
     deltaEta = part1.eta() - part2.eta();
     deltaPhi = part1.phi() - part2.phi();
 
-    deltaPhi = RecoDecay::constrainAngle(deltaPhi, 0);
+    while (deltaPhi < mPhiLow) {
+      deltaPhi += o2::constants::math::TwoPI;
+    }
+    while (deltaPhi > mPhiHigh) {
+      deltaPhi -= o2::constants::math::TwoPI;
+    }
 
     mHistogramRegistry->fill(HIST(FolderSuffix[EventType]) + HIST(o2::aod::femtouniverse_mc_particle::MCTypeName[mc]) + HIST("/relPairDist"), femtoObs, weight);
     mHistogramRegistry->fill(HIST(FolderSuffix[EventType]) + HIST(o2::aod::femtouniverse_mc_particle::MCTypeName[mc]) + HIST("/relPairkT"), kT, weight);

PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackPhi.cxx

@@ -454,14 +454,14 @@ struct FemtoUniversePairTaskTrackPhi {
     registryMCreco.add("MCrecoPpos", "MC reco proton;#it{p}_{T} (GeV/c); #eta", {HistType::kTH2F, {{500, 0, 5}, {400, -1.0, 1.0}}});
     registryMCreco.add("MCrecoPneg", "MC reco antiproton;#it{p}_{T} (GeV/c); #eta", {HistType::kTH2F, {{500, 0, 5}, {400, -1.0, 1.0}}});
 
-    trackHistoPartPhi.init(&qaRegistry, confBinsTempFitVarpT.value, confBinsTempFitVarInvMass.value, ConfBothTracks.confIsMC.value, 333);
-    trackHistoPartTrack.init(&qaRegistry, confBinsTempFitVarpT.value, confBinsTempFitVar.value, ConfBothTracks.confIsMC.value, ConfTrack.confTrackPDGCode.value);
+    trackHistoPartPhi.init(&qaRegistry, confBinsTempFitVarpT, confBinsTempFitVarInvMass, ConfBothTracks.confIsMC.value, 333);
+    trackHistoPartTrack.init(&qaRegistry, confBinsTempFitVarpT, confBinsTempFitVar, ConfBothTracks.confIsMC.value, ConfTrack.confTrackPDGCode.value);
 
     mixQaRegistry.add("MixingQA/hSECollisionBins", ";bin;Entries", kTH1F, {{120, -0.5, 119.5}});
     mixQaRegistry.add("MixingQA/hMECollisionBins", ";bin;Entries", kTH1F, {{120, -0.5, 119.5}});
 
-    sameEventAngularCont.init(&resultRegistry, confBinskstar.value, confBinsMult.value, confBinskT.value, confBinsmT.value, confBins3Dmult.value, confBins3DmT.value, ConfBothTracks.confBinsEta.value, ConfBothTracks.confBinsPhi.value, ConfBothTracks.confIsMC.value, ConfBothTracks.confUse3D.value);
-    mixedEventAngularCont.init(&resultRegistry, confBinskstar.value, confBinsMult.value, confBinskT.value, confBinsmT.value, confBins3Dmult.value, confBins3DmT.value, ConfBothTracks.confBinsEta.value, ConfBothTracks.confBinsPhi.value, ConfBothTracks.confIsMC.value, ConfBothTracks.confUse3D.value);
+    sameEventAngularCont.init(&resultRegistry, confBinskstar, confBinsMult, confBinskT, confBinsmT, confBins3Dmult, confBins3DmT, ConfBothTracks.confBinsEta, ConfBothTracks.confBinsPhi, ConfBothTracks.confIsMC, ConfBothTracks.confUse3D);
+    mixedEventAngularCont.init(&resultRegistry, confBinskstar, confBinsMult, confBinskT, confBinsmT, confBins3Dmult, confBins3DmT, ConfBothTracks.confBinsEta, ConfBothTracks.confBinsPhi, ConfBothTracks.confIsMC, ConfBothTracks.confUse3D);
 
     sameEventAngularCont.setPDGCodes(333, ConfTrack.confTrackPDGCode.value);
     mixedEventAngularCont.setPDGCodes(333, ConfTrack.confTrackPDGCode.value);