Make shower exclusion fcl-toggled but defaulted on

Author: henrylay97

sbncode/Calibration/TrackCaloSkimmer.h

@@ -223,6 +223,7 @@ class sbn::TrackCaloSkimmer : public art::EDAnalyzer {
   double fTopCRTDistanceCutPassing;
   double fSideCRTDistanceCutStopping;
   double fSideCRTDistanceCutPassing;  
+  bool fAllowShowerLikePFPs;
 
   // tools
   std::vector<std::unique_ptr<sbn::ITCSSelectionTool>> fSelectionTools;

sbncode/Calibration/TrackCaloSkimmer_module.cc

@@ -92,7 +92,8 @@ sbn::TrackCaloSkimmer::TrackCaloSkimmer(fhicl::ParameterSet const& p)
   fTopCRTDistanceCutPassing = p.get<double>("TopCRTDistanceCut_throughgoing", 100.);
   fSideCRTDistanceCutStopping = p.get<double>("SideCRTDistanceCut_stopping", 100.);
   fSideCRTDistanceCutPassing = p.get<double>("SideCRTDistanceCut_throughgoing", 100.);
-  
+  fAllowShowerLikePFPs = p.get<bool>("AllowShowerLikePFPs", false);
+
   if (fTailFitResidualRange > 10.) {
     std::cout << "sbn::TrackCaloSkimmer: Bad tail fit residual range config :(" << fTailFitResidualRange << "). Fits will not be meaningful.\n";
   }
@@ -315,7 +316,7 @@ void sbn::TrackCaloSkimmer::analyze(art::Event const& e)
   for (art::Ptr<recob::PFParticle> p_pfp: PFParticleList) {
     const recob::PFParticle &pfp = *p_pfp;
 
-    if(p_pfp->PdgCode() == 11)
+    if(p_pfp->PdgCode() == 11 && !fAllowShowerLikePFPs)
       continue;
 
     const std::vector<art::Ptr<recob::Track>> thisTrack = fmTracks.at(p_pfp.key());