diff --git a/setup.sh b/setup.sh
index 8dea0627324e17d80dd117bec5b16a307feed42c..40260e59864bc4a2a6f7e30fcf733e12b02fdf94 100644
--- a/setup.sh
+++ b/setup.sh
@@ -7,4 +7,5 @@ export HCANALYZER=${MYDIR}
 if [ ! ${LD_LIBRARY_PATH} ]; then
   export LD_LIBRARY_PATH=""
 fi
-export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:$ANALYZER:$HCANALYZER
+#export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:$ANALYZER:$HCANALYZER
+export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:$ANALYZER:$HCANALYZER:$PATH
diff --git a/src/THcShower.cxx b/src/THcShower.cxx
index 558ce075c3ce51e2ce8070e04295c4644cb5703f..82a35561354a5710bfff0b899c112675c47730fb 100644
--- a/src/THcShower.cxx
+++ b/src/THcShower.cxx
@@ -425,6 +425,7 @@ Int_t THcShower::DefineVariables( EMode mode )
  //   { "asum_p", "Sum of ped-subtracted ADCs",         "fAsum_p" },
  //   { "asum_c", "Sum of calibrated ADCs",             "fAsum_c" },
     { "nclust", "Number of clusters",                            "fNclust" },
+    { "ntracks", "Number of shower tracks",                      "fNtracks" },
     { "emax",   "Energy of largest cluster",                     "fE" },
     { "eprmax",   "Preshower Energy of largest cluster",         "fEpr" },
     { "xmax",      "x-position (cm) of largest cluster",         "fX" },
@@ -501,6 +502,7 @@ void THcShower::Clear(Option_t* opt)
 
   fNhits = 0;
   fNclust = 0;
+  fNtracks = 0;
   fMult = 0;
   fE = 0.;
   fEpr = 0.;
@@ -757,6 +759,8 @@ Int_t THcShower::CoarseProcess( TClonesArray& tracks)
 
     Int_t mclust = MatchCluster(theTrack, ClusterList, Xtr, Ytr);
 
+    if (mclust >= 0) fNtracks++;  // number of shower tracks
+
     // Do this for the 1-st track only for now.
     //
     if (itrk==0) {
diff --git a/src/THcShower.h b/src/THcShower.h
index 3ce02c7857d78b49149c864469105039bdafd695..fa13670a9f516fc296ed0986f14284c5bc13e912 100644
--- a/src/THcShower.h
+++ b/src/THcShower.h
@@ -127,6 +127,8 @@ protected:
 
   Int_t fNhits;              // Total number of hits
   Int_t fNclust;             // Number of clusters
+  Int_t fNtracks;            // Number of shower tracks, i.e. tracks with
+                             // associated clusters
   Double_t fE;               // Energy of the largest cluster
   Double_t fEpr;             // Preshower Energy of the largest cluster
   Double_t fX;               // x-position (cm) of the largest cluster