From b12e1dbed59d1492f697ddf42b4d3d710d1661f1 Mon Sep 17 00:00:00 2001 From: Mario Ciacco Date: Mon, 23 Jun 2025 14:33:21 +0200 Subject: [PATCH 1/2] fix histogram filling using eta --- PWGLF/TableProducer/Strangeness/strangeTreeCreator.cxx | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/PWGLF/TableProducer/Strangeness/strangeTreeCreator.cxx b/PWGLF/TableProducer/Strangeness/strangeTreeCreator.cxx index a18397f214e..cd1a40a9379 100644 --- a/PWGLF/TableProducer/Strangeness/strangeTreeCreator.cxx +++ b/PWGLF/TableProducer/Strangeness/strangeTreeCreator.cxx @@ -457,14 +457,16 @@ struct StrangeTreeCreator { } float ptPos = std::hypot(momPos[0], momPos[1]); + float momPos = std::hypot(momPos[0], momPos[1], momPos[2]); + float etaPos = 0.5 * std::log((momPos + momPos[2]) / (momPos - momPos[2])); float deltaMass = mK0Short - o2::constants::physics::MassK0; float massKineBias = deltaMass / kineFactor(momV0, momPos, momNeg, o2::constants::physics::MassPiMinus, o2::constants::physics::MassPiMinus, false); float massKineReso = deltaMass / kineFactor(momV0, momPos, momNeg, o2::constants::physics::MassPiMinus, o2::constants::physics::MassPiMinus, true); histos.fill(HIST("QA/massLambda"), centrality, ptV0, mLambda); histos.fill(HIST("QA/massK0s"), ptV0, mK0Short); - histos.fill(HIST("massKineBias"), ptPos, massKineBias); - histos.fill(HIST("massKineReso"), ptPos, massKineReso); + histos.fill(HIST("massKineBias"), etaPos, ptPos, massKineBias); + histos.fill(HIST("massKineReso"), etaPos, ptPos, massKineReso); CandidateV0 candV0; candV0.pt = matter > 0. ? ptV0 : -ptV0; From eabeccd78dc5beb57bc16a9d063d7b2c877f7e5f Mon Sep 17 00:00:00 2001 From: Mario Ciacco Date: Mon, 23 Jun 2025 14:43:21 +0200 Subject: [PATCH 2/2] change variable name to avoid clash with existing one --- PWGLF/TableProducer/Strangeness/strangeTreeCreator.cxx | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/PWGLF/TableProducer/Strangeness/strangeTreeCreator.cxx b/PWGLF/TableProducer/Strangeness/strangeTreeCreator.cxx index cd1a40a9379..455da9e0719 100644 --- a/PWGLF/TableProducer/Strangeness/strangeTreeCreator.cxx +++ b/PWGLF/TableProducer/Strangeness/strangeTreeCreator.cxx @@ -457,8 +457,8 @@ struct StrangeTreeCreator { } float ptPos = std::hypot(momPos[0], momPos[1]); - float momPos = std::hypot(momPos[0], momPos[1], momPos[2]); - float etaPos = 0.5 * std::log((momPos + momPos[2]) / (momPos - momPos[2])); + float pPos = std::hypot(momPos[0], momPos[1], momPos[2]); + float etaPos = 0.5 * std::log((pPos + momPos[2]) / (pPos - momPos[2])); float deltaMass = mK0Short - o2::constants::physics::MassK0; float massKineBias = deltaMass / kineFactor(momV0, momPos, momNeg, o2::constants::physics::MassPiMinus, o2::constants::physics::MassPiMinus, false); float massKineReso = deltaMass / kineFactor(momV0, momPos, momNeg, o2::constants::physics::MassPiMinus, o2::constants::physics::MassPiMinus, true);