Add number of calorimeter tracks (track associated clusters)

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

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) {

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