From 340c10d6c4585cfcd6eb266f332292826204da58 Mon Sep 17 00:00:00 2001
From: Christian Herwig <therwig@users.noreply.github.com>
Date: Mon, 6 Nov 2023 08:13:42 -0800
Subject: [PATCH] set good defaults

---
 Recon/src/Recon/PFEcalClusterProducer.cxx | 2 +-
 Recon/src/Recon/PFHcalClusterProducer.cxx | 2 +-
 Recon/src/Recon/PFTrackProducer.cxx       | 2 +-
 Recon/src/Recon/ParticleFlow.cxx          | 2 --
 4 files changed, 3 insertions(+), 5 deletions(-)

diff --git a/Recon/src/Recon/PFEcalClusterProducer.cxx b/Recon/src/Recon/PFEcalClusterProducer.cxx
index c163c3f63..c631e21b9 100644
--- a/Recon/src/Recon/PFEcalClusterProducer.cxx
+++ b/Recon/src/Recon/PFEcalClusterProducer.cxx
@@ -18,7 +18,7 @@ void PFEcalClusterProducer::configure(framework::config::Parameters& ps) {
   //DBScan parameters
   minClusterHitMult_ = ps.getParameter<int>("minClusterHitMult");
   clusterHitDist_ = ps.getParameter<double>("clusterHitDist");
-  clusterZBias_ = ps.getParameter<double>("clusterZBias");
+  clusterZBias_ = ps.getParameter<double>("clusterZBias",1);
   minHitEnergy_ = ps.getParameter<double>("minHitEnergy");
 }
 
diff --git a/Recon/src/Recon/PFHcalClusterProducer.cxx b/Recon/src/Recon/PFHcalClusterProducer.cxx
index 87b7b1dfe..230c4a8b4 100644
--- a/Recon/src/Recon/PFHcalClusterProducer.cxx
+++ b/Recon/src/Recon/PFHcalClusterProducer.cxx
@@ -22,7 +22,7 @@ void PFHcalClusterProducer::configure(framework::config::Parameters& ps) {
   //DBScan parameters
   minClusterHitMult_ = ps.getParameter<int>("minClusterHitMult");
   clusterHitDist_ = ps.getParameter<double>("clusterHitDist");
-  clusterZBias_ = ps.getParameter<double>("clusterZBias");
+  clusterZBias_ = ps.getParameter<double>("clusterZBias",1);
   minHitEnergy_ = ps.getParameter<double>("minHitEnergy");
 }
 
diff --git a/Recon/src/Recon/PFTrackProducer.cxx b/Recon/src/Recon/PFTrackProducer.cxx
index d3e9a381a..a0306e1e7 100644
--- a/Recon/src/Recon/PFTrackProducer.cxx
+++ b/Recon/src/Recon/PFTrackProducer.cxx
@@ -25,7 +25,7 @@ void PFTrackProducer::produce(framework::Event& event) {
   if(truthTracking_){ 
     for(const auto &spHit : ecalSpHits){
       if ( spHit.getTrackID()!=1 || fabs(240-spHit.getPosition()[2])>0.1  || spHit.getMomentum()[2] <= 0 ) continue;
-      if ( spHit.getPdgID() == 22 || spHit.getPdgID() == 2112 ) continue
+      if ( spHit.getPdgID() == 22 || spHit.getPdgID() == 2112 ) continue;
       pfTracks.push_back(spHit);
       break;
     }
diff --git a/Recon/src/Recon/ParticleFlow.cxx b/Recon/src/Recon/ParticleFlow.cxx
index 4b79f1a4a..62abea4bc 100644
--- a/Recon/src/Recon/ParticleFlow.cxx
+++ b/Recon/src/Recon/ParticleFlow.cxx
@@ -238,7 +238,6 @@ void ParticleFlow::produce(framework::Event& event) {
     // std::vector<ldmx::PFCandidate> chargedMatch;
     // std::vector<ldmx::PFCandidate> chargedUnmatch;
     for(int i=0; i<tracks.size(); i++){
-      //if (tracks[i].getPdgID()==22 || tracks[i].getPdgID()==2112) continue;  // abandon build track if photon or neutron
       ldmx::PFCandidate cand;
       FillCandTrack(cand, tracks[i]);  // append track info to candidate
       
@@ -259,7 +258,6 @@ void ParticleFlow::produce(framework::Event& event) {
     for(int i=0; i<ecalClusters.size(); i++){
       // already linked with ECal in the previous step
       if(EMIsTkLinked[i]) continue;
-      //if((EMIsTkLinked[i]) && (!(tracks[i].getPdgID()==22 || tracks[i].getPdgID()==2112))) continue;  // must also be charged for EMIsTkLinked to make sense
 
       ldmx::PFCandidate cand;
       FillCandEMCalo(cand, ecalClusters[i]);