set good defaults

LDMX-Software · Nov 6, 2023 · 340c10d · 340c10d
1 parent c5ee6c4
commit 340c10d
Show file tree

Hide file tree

Showing 4 changed files with 3 additions and 5 deletions.
diff --git 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
@@ -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
@@ -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
@@ -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]);