From 3d1835187a068237ca4e27d4c322db25003b7355 Mon Sep 17 00:00:00 2001 From: blacw Date: Sun, 14 Sep 2025 18:43:27 +0800 Subject: [PATCH 1/4] add DCA fitter switch --- PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx | 34 ++++++++++--------- 1 file changed, 18 insertions(+), 16 deletions(-) diff --git a/PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx b/PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx index f50262bd998..f4649b049f2 100644 --- a/PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx +++ b/PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx @@ -609,22 +609,24 @@ struct he3HadronFemto { he3Hadcand.signHe3 = trackHe3.sign(); he3Hadcand.signHad = trackHad.sign(); - - // he3Hadcand.dcaxyHe3 = trackHe3.dcaXY(); - // he3Hadcand.dcaxyHad = trackHad.dcaXY(); - // he3Hadcand.dcazHe3 = trackHe3.dcaZ(); - // he3Hadcand.dcazHad = trackHad.dcaZ(); - auto trackCovHe3 = getTrackParCov(trackHe3); - auto trackCovHad = getTrackParCov(trackHad); - std::array dcaInfo; - o2::base::Propagator::Instance()->propagateToDCABxByBz({collisionVertex[0], collisionVertex[1], collisionVertex[2]}, trackCovHe3, 2.f, mFitter.getMatCorrType(), &dcaInfo); - he3Hadcand.dcaxyHe3 = dcaInfo[0]; - he3Hadcand.dcazHe3 = dcaInfo[1]; - o2::base::Propagator::Instance()->propagateToDCABxByBz({collisionVertex[0], collisionVertex[1], collisionVertex[2]}, trackCovHad, 2.f, mFitter.getMatCorrType(), &dcaInfo); - he3Hadcand.dcaxyHad = dcaInfo[0]; - he3Hadcand.dcazHad = dcaInfo[1]; - he3Hadcand.dcaPair = std::sqrt(std::abs(mFitter.getChi2AtPCACandidate())); - + if(!settingEnableDCAfitter){ + he3Hadcand.dcaxyHe3 = trackHe3.dcaXY(); + he3Hadcand.dcaxyHad = trackHad.dcaXY(); + he3Hadcand.dcazHe3 = trackHe3.dcaZ(); + he3Hadcand.dcazHad = trackHad.dcaZ(); + } else { + auto trackCovHe3 = getTrackParCov(trackHe3); + auto trackCovHad = getTrackParCov(trackHad); + std::array dcaInfo; + o2::base::Propagator::Instance()->propagateToDCABxByBz({collisionVertex[0], collisionVertex[1], collisionVertex[2]}, trackCovHe3, 2.f, mFitter.getMatCorrType(), &dcaInfo); + he3Hadcand.dcaxyHe3 = dcaInfo[0]; + he3Hadcand.dcazHe3 = dcaInfo[1]; + o2::base::Propagator::Instance()->propagateToDCABxByBz({collisionVertex[0], collisionVertex[1], collisionVertex[2]}, trackCovHad, 2.f, mFitter.getMatCorrType(), &dcaInfo); + he3Hadcand.dcaxyHad = dcaInfo[0]; + he3Hadcand.dcazHad = dcaInfo[1]; + he3Hadcand.dcaPair = std::sqrt(std::abs(mFitter.getChi2AtPCACandidate())); + } + he3Hadcand.tpcSignalHe3 = trackHe3.tpcSignal(); bool heliumPID = trackHe3.pidForTracking() == o2::track::PID::Helium3 || trackHe3.pidForTracking() == o2::track::PID::Alpha; float correctedTPCinnerParamHe3 = (heliumPID && settingCompensatePIDinTracking) ? trackHe3.tpcInnerParam() / 2.f : trackHe3.tpcInnerParam(); From f753af616c5032837a672ffb881d5c5d51a6d8ac Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Sun, 14 Sep 2025 10:52:51 +0000 Subject: [PATCH 2/4] Please consider the following formatting changes --- PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx b/PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx index f4649b049f2..3889fdcfc9f 100644 --- a/PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx +++ b/PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx @@ -609,11 +609,11 @@ struct he3HadronFemto { he3Hadcand.signHe3 = trackHe3.sign(); he3Hadcand.signHad = trackHad.sign(); - if(!settingEnableDCAfitter){ - he3Hadcand.dcaxyHe3 = trackHe3.dcaXY(); - he3Hadcand.dcaxyHad = trackHad.dcaXY(); - he3Hadcand.dcazHe3 = trackHe3.dcaZ(); - he3Hadcand.dcazHad = trackHad.dcaZ(); + if (!settingEnableDCAfitter) { + he3Hadcand.dcaxyHe3 = trackHe3.dcaXY(); + he3Hadcand.dcaxyHad = trackHad.dcaXY(); + he3Hadcand.dcazHe3 = trackHe3.dcaZ(); + he3Hadcand.dcazHad = trackHad.dcaZ(); } else { auto trackCovHe3 = getTrackParCov(trackHe3); auto trackCovHad = getTrackParCov(trackHad); @@ -626,7 +626,7 @@ struct he3HadronFemto { he3Hadcand.dcazHad = dcaInfo[1]; he3Hadcand.dcaPair = std::sqrt(std::abs(mFitter.getChi2AtPCACandidate())); } - + he3Hadcand.tpcSignalHe3 = trackHe3.tpcSignal(); bool heliumPID = trackHe3.pidForTracking() == o2::track::PID::Helium3 || trackHe3.pidForTracking() == o2::track::PID::Alpha; float correctedTPCinnerParamHe3 = (heliumPID && settingCompensatePIDinTracking) ? trackHe3.tpcInnerParam() / 2.f : trackHe3.tpcInnerParam(); From 6664838a499377f9e2b2696016b2daf201c8fa0b Mon Sep 17 00:00:00 2001 From: blacwovie Date: Sun, 14 Sep 2025 19:14:14 +0800 Subject: [PATCH 3/4] Update he3HadronFemto.cxx --- PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx | 1 + 1 file changed, 1 insertion(+) diff --git a/PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx b/PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx index 3889fdcfc9f..2fdb955b98c 100644 --- a/PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx +++ b/PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx @@ -65,6 +65,7 @@ #include // std::prev #include #include +#include using namespace o2; using namespace o2::framework; From 2ebe6657ec3114b5f9893226961bddc26ab0c56c Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Sun, 14 Sep 2025 11:14:45 +0000 Subject: [PATCH 4/4] Please consider the following formatting changes --- PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx b/PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx index 2fdb955b98c..975572e927e 100644 --- a/PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx +++ b/PWGLF/TableProducer/Nuspex/he3HadronFemto.cxx @@ -64,8 +64,8 @@ #include #include // std::prev #include -#include #include +#include using namespace o2; using namespace o2::framework;