From eb5a590ef37e45d5f3a5c0d61798a7752ad09aaf Mon Sep 17 00:00:00 2001
From: Fangyi Guo <guofangyi@ihep.ac.cn>
Date: Fri, 7 Feb 2025 09:58:46 +0800
Subject: [PATCH 1/2] Update reclustering parameters
---
.../RecPFACyber/src/Algorithm/PFOReclusteringAlg.cpp | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/Reconstruction/RecPFACyber/src/Algorithm/PFOReclusteringAlg.cpp b/Reconstruction/RecPFACyber/src/Algorithm/PFOReclusteringAlg.cpp
index 996c56e3..73e77e32 100644
--- a/Reconstruction/RecPFACyber/src/Algorithm/PFOReclusteringAlg.cpp
+++ b/Reconstruction/RecPFACyber/src/Algorithm/PFOReclusteringAlg.cpp
@@ -17,8 +17,8 @@ StatusCode PFOReclusteringAlg::ReadSettings(Settings& m_settings){
if(settings.map_floatPars.find("SplitSigma")==settings.map_floatPars.end()) settings.map_floatPars["SplitSigma"] = 0.2;
if(settings.map_floatPars.find("NeutralMergeSigma")==settings.map_floatPars.end()) settings.map_floatPars["NeutralMergeSigma"] = 0.;
if(settings.map_floatPars.find("VirtualMergeSigma")==settings.map_floatPars.end()) settings.map_floatPars["VirtualMergeSigma"] = 0.4;
- if(settings.map_floatPars.find("MinAngleForNeuMerge")==settings.map_floatPars.end()) settings.map_floatPars["MinAngleForNeuMerge"] = 0.20;
- if(settings.map_floatPars.find("MinAngleForVirMerge")==settings.map_floatPars.end()) settings.map_floatPars["MinAngleForVirMerge"] = 0.24;
+ if(settings.map_floatPars.find("MinAngleForNeuMerge")==settings.map_floatPars.end()) settings.map_floatPars["MinAngleForNeuMerge"] = 0.18;
+ if(settings.map_floatPars.find("MinAngleForVirMerge")==settings.map_floatPars.end()) settings.map_floatPars["MinAngleForVirMerge"] = 0.20;
return StatusCode::SUCCESS;
--
GitLab
From d18198dab39024b577b3f55ec2d0820dba0f16b0 Mon Sep 17 00:00:00 2001
From: Fangyi Guo <guofangyi@ihep.ac.cn>
Date: Fri, 7 Feb 2025 09:59:24 +0800
Subject: [PATCH 2/2] Update charged PFO energy with mass
---
Reconstruction/RecPFACyber/src/Tools/OutputCreator.cpp | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/Reconstruction/RecPFACyber/src/Tools/OutputCreator.cpp b/Reconstruction/RecPFACyber/src/Tools/OutputCreator.cpp
index 43ed4967..3016a262 100644
--- a/Reconstruction/RecPFACyber/src/Tools/OutputCreator.cpp
+++ b/Reconstruction/RecPFACyber/src/Tools/OutputCreator.cpp
@@ -160,8 +160,8 @@ namespace Cyber{
TVector3 p3vec = vec_trks[trkIndex]->getP3();
edm4hep::Vector3f p3(p3vec.x(), p3vec.y(), p3vec.z());
m_pfo.setMomentum(p3);
- m_pfo.setEnergy( vec_trks[trkIndex]->getMomentum() );
- m_pfo.setMass(139.57039); //TODO: all charged particles are set to pion mass.
+ m_pfo.setMass(0.139570); //TODO: all charged particles are set to pion mass.
+ m_pfo.setEnergy( sqrt(vec_trks[trkIndex]->getMomentum()*vec_trks[trkIndex]->getMomentum() + 0.139570*0.139570) );
}
else{
TVector3 p3vec = vec_Pos*( (EcalClusE+HcalClusE)/vec_Pos.Mag() );
--
GitLab