diff --git a/PWGCF/Flow/Tasks/CMakeLists.txt b/PWGCF/Flow/Tasks/CMakeLists.txt index 0f5f97f64fb..64f235e7337 100644 --- a/PWGCF/Flow/Tasks/CMakeLists.txt +++ b/PWGCF/Flow/Tasks/CMakeLists.txt @@ -114,7 +114,7 @@ o2physics_add_dpl_workflow(pid-flow-pt-corr PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore O2Physics::AnalysisCCDB O2Physics::GFWCore COMPONENT_NAME Analysis) -o2physics_add_dpl_workflow(flow-gpol-ese +o2physics_add_dpl_workflow(flow-directed-flow-task SOURCES flowDirectedFlowTask.cxx PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore COMPONENT_NAME Analysis) diff --git a/PWGCF/Flow/Tasks/flowDirectedFlowTask.cxx b/PWGCF/Flow/Tasks/flowDirectedFlowTask.cxx index 28f32417010..253174ef5bb 100644 --- a/PWGCF/Flow/Tasks/flowDirectedFlowTask.cxx +++ b/PWGCF/Flow/Tasks/flowDirectedFlowTask.cxx @@ -146,7 +146,6 @@ struct flowDirectedFlowTask { using EventCandidates = soa::Filtered>; using AllTrackCandidates = soa::Filtered>; using ResoV0s = aod::V0Datas; - using BCsRun3 = soa::Join; void init(InitContext&) { @@ -222,17 +221,6 @@ struct flowDirectedFlowTask { histos.add("hSparseAntiLambdaAvgUyQ1", "AntiLambda vs q1", HistType::kTHnSparseF, axesPolSPQ1, true); } - float getPhiInRange(float phi) - { - while (phi <= -o2::constants::math::PI) { - phi += o2::constants::math::TwoPI; - } - while (phi > o2::constants::math::PI) { - phi -= o2::constants::math::TwoPI; - } - return phi; - } - template bool eventSelected(TCollision collision) { @@ -388,9 +376,6 @@ struct flowDirectedFlowTask { if (std::abs(v0.yLambda()) > cfgV0Rap) { return false; } - if (std::abs(v0.eta()) > 0.8) { - return false; - } return true; } @@ -439,12 +424,12 @@ struct flowDirectedFlowTask { const float cosThetaStar = daughterStar.Pz() / daughterStar.P(); const float sinThetaStar = std::sqrt(std::max(0.0, 1.0 - cosThetaStar * cosThetaStar)); - const float sinPhiStar = std::sin(getPhiInRange(phiStar)); - const float cosPhiStar = std::cos(getPhiInRange(phiStar)); + const float sinPhiStar = std::sin(phiStar); + const float cosPhiStar = std::cos(phiStar); - float polEP_A = std::sin(getPhiInRange(phiStar - psiA)); - float polEP_C = std::sin(getPhiInRange(phiStar - psiC)); - float polEP = std::sin(getPhiInRange(phiStar - psiFull)); + float polEP_A = std::sin(phiStar - psiA); + float polEP_C = std::sin(phiStar - psiC); + float polEP = std::sin(phiStar - psiFull); const float qxFull = qxC - qxA; const float qyFull = qyC - qyA; @@ -464,8 +449,8 @@ struct flowDirectedFlowTask { polEP_A /= accDen; polEP_C /= accDen; - const float cosPsi = std::cos(getPhiInRange(psiFull)); - const float sinPsi = std::sin(getPhiInRange(psiFull)); + const float cosPsi = std::cos(psiFull); + const float sinPsi = std::sin(psiFull); if (isLambda) { histos.fill(HIST("hSparseLambdaPolSPQ1"), mass, v0.pt(), polSP, centrality, q1, wgt); @@ -502,7 +487,7 @@ struct flowDirectedFlowTask { } } - void processData(EventCandidates::iterator const& collision, AllTrackCandidates const& tracks, ResoV0s const& v0s, BCsRun3 const&) + void processData(EventCandidates::iterator const& collision, AllTrackCandidates const& tracks, ResoV0s const& v0s, aod::BCsWithTimestamps const&) { if (!eventSelected(collision)) { return; @@ -539,7 +524,7 @@ struct flowDirectedFlowTask { float dotAC = qxA * qxC + qyA * qyC; float resDot = dotAC / (magA * magC); - histos.fill(HIST("hpResCosAC"), centrality, std::cos(getPhiInRange(psiA - psiC))); + histos.fill(HIST("hpResCosAC"), centrality, std::cos(psiA - psiC)); histos.fill(HIST("hpResDotAC"), centrality, resDot); histos.fill(HIST("hpQxAQxC"), centrality, qxA * qxC); histos.fill(HIST("hpQyAQyC"), centrality, qyA * qyC); @@ -563,9 +548,9 @@ struct flowDirectedFlowTask { histos.fill(HIST("hV1SPAQ1"), centrality, track.pt(), track.eta(), q1, v1SPA); histos.fill(HIST("hV1SPCQ1"), centrality, track.pt(), track.eta(), q1, v1SPC); - histos.fill(HIST("hV1EPFullQ1"), centrality, track.pt(), track.eta(), q1, std::cos(getPhiInRange(phi - psiFull))); - histos.fill(HIST("hV1EPAQ1"), centrality, track.pt(), track.eta(), q1, std::cos(getPhiInRange(phi - psiA))); - histos.fill(HIST("hV1EPCQ1"), centrality, track.pt(), track.eta(), q1, std::cos(getPhiInRange(phi - psiC))); + histos.fill(HIST("hV1EPFullQ1"), centrality, track.pt(), track.eta(), q1, std::cos(phi - psiFull)); + histos.fill(HIST("hV1EPAQ1"), centrality, track.pt(), track.eta(), q1, std::cos(phi - psiA)); + histos.fill(HIST("hV1EPCQ1"), centrality, track.pt(), track.eta(), q1, std::cos(phi - psiC)); } for (const auto& v0 : v0s) {