Implemented the clustering debug flag (Mark).

examples/hodtest_mkj.C

@@ -7,7 +7,8 @@ void hodtest_mkj(Int_t RunNumber=50017, Int_t MaxEventToReplay=5000) {
   if (RunNumber == 50017) {
     char* RunFileNamePattern="daq04_%d.log.0";
   } else {
-    char* RunFileNamePattern="/cache/mss/hallc/daq04/raw/daq04_%d.log.0";
+    //    char* RunFileNamePattern="/cache/mss/hallc/daq04/raw/daq04_%d.log.0";
+    char* RunFileNamePattern="/scratch/tadevosn/daq04_%d.log.0";
   }
   gHcParms->Define("gen_run_number", "Run Number", RunNumber);
   gHcParms->AddString("g_ctp_database_filename", "DBASE/test.database");

src/THcShower.cxx

@@ -161,6 +161,7 @@ Int_t THcShower::ReadDatabase( const TDatime& date )
       {"cal_num_neg_columns", &fNegCols, kInt},
       {"cal_slop", &fSlop, kDouble},
       {"cal_fv_test", &fvTest, kDouble},
+      {"dbg_clusters_cal", &fdbg_clusters_cal, kInt},
       {0}
     };
     gHcParms->LoadParmValues((DBRequest*)&list, prefix);
@@ -169,6 +170,7 @@ Int_t THcShower::ReadDatabase( const TDatime& date )
   cout << "Number of neg. columns   = " << fNegCols << endl;
   cout << "Slop parameter           = " << fSlop << endl;
   cout << "Fiducial volum test flag = " << fvTest << endl;
+  cout << "Cluster debug flag       = " << fdbg_clusters_cal  << endl;
 
   BlockThick = new Double_t [fNLayers];
   fNBlocks = new Int_t [fNLayers];
@@ -524,6 +526,7 @@ Int_t THcShower::CoarseProcess( TClonesArray&  ) //tracks
   //
   //  static const Double_t sqrt2 = TMath::Sqrt(2.);
   
+  if (fdbg_clusters_cal)
   cout << "THcShower::CoarseProcess called ---------------------------" <<endl;
 
   //  ApplyCorrections();
@@ -573,37 +576,43 @@ Int_t THcShower::CoarseProcess( TClonesArray&  ) //tracks
   //Print out hits before clustering.
   //
   fNhits = HitList.size();
-  cout << "Total hits:     " << fNhits << endl;
-  for (UInt_t i=0; i!=fNhits; i++) {
-    cout << "unclustered hit " << i << ": ";
-    (*(HitList.begin()+i))->show();
+
+  if (fdbg_clusters_cal) {
+    cout << "Total hits:     " << fNhits << endl;
+    for (UInt_t i=0; i!=fNhits; i++) {
+      cout << "unclustered hit " << i << ": ";
+      (*(HitList.begin()+i))->show();
+    }
   }
 
   THcShowerClusterList* ClusterList = new THcShowerClusterList;
   ClusterList->ClusterHits(HitList);
 
+  fNclust = (*ClusterList).NbClusters();
+
   //Print out the cluster list.
   //
-  fNclust = (*ClusterList).NbClusters();
-  cout << "Cluster_list size: " << fNclust << endl;
+  if (fdbg_clusters_cal) {
+    cout << "Cluster_list size: " << fNclust << endl;
 
-  for (UInt_t i=0; i!=fNclust; i++) {
+    for (UInt_t i=0; i!=fNclust; i++) {
 
-    THcShowerCluster* cluster = (*ClusterList).ListedCluster(i);
+      THcShowerCluster* cluster = (*ClusterList).ListedCluster(i);
 
-    cout << "Cluster #" << i 
-         <<":  E=" << (*cluster).clE() 
-         << "  Epr=" << (*cluster).clEpr()
-         << "  X=" << (*cluster).clX()
-         << "  Z=" << (*cluster).clZ()
-         << "  size=" << (*cluster).clSize()
-         << endl;
+      cout << "Cluster #" << i 
+	   <<":  E=" << (*cluster).clE() 
+	   << "  Epr=" << (*cluster).clEpr()
+	   << "  X=" << (*cluster).clX()
+	   << "  Z=" << (*cluster).clZ()
+	   << "  size=" << (*cluster).clSize()
+	   << endl;
+
+      for (UInt_t j=0; j!=(*cluster).clSize(); j++) {
+	THcShowerHit* hit = (*cluster).ClusteredHit(j);
+	cout << "  hit #" << j << ":  "; (*hit).show();
+      }
 
-    for (UInt_t j=0; j!=(*cluster).clSize(); j++) {
-      THcShowerHit* hit = (*cluster).ClusteredHit(j);
-      cout << "  hit #" << j << ":  "; (*hit).show();
     }
-
   }
 
   
@@ -612,9 +621,6 @@ Int_t THcShower::CoarseProcess( TClonesArray&  ) //tracks
 
   if (fNclust != 0) {
 
-    //    THcShowerCluster* MaxCluster;
-    //    for (UInt_t i=0; i!=fNclust; i++) {
-
     THcShowerCluster* MaxCluster =  (*ClusterList).ListedCluster(0);
     fE = (*MaxCluster).clE();
 
@@ -636,19 +642,12 @@ Int_t THcShower::CoarseProcess( TClonesArray&  ) //tracks
     fX = (*MaxCluster).clX();
   }
 
-  cout << fEpr << " " << fE << " " << fX << " PrSh" << endl;
-
-  /*
-    cout << "Cluster #" << i 
-         <<":  E=" << (*cluster).clE() 
-         << "  Epr=" << (*cluster).clEpr()
-         << "  X=" << (*cluster).clX()
-         << "  Z=" << (*cluster).clZ()
-         << "  size=" << (*cluster).clSize()
-         << endl;
-  */
-   
+  if (fdbg_clusters_cal)
+    cout << fEpr << " " << fE << " " << fX << " PrSh" << endl;
+
+  if (fdbg_clusters_cal)
   cout << "THcShower::CoarseProcess return ---------------------------" <<endl;
+
   return 0;
 }
 

src/THcShower.h

@@ -110,6 +110,8 @@ public:
   Double_t fSlop;               // Track to cluster vertical slop distance.
   Int_t fvTest;                 // fiducial volume test flag for tracking
 
+  Int_t fdbg_clusters_cal;      // Shower debug flag
+
   THcShowerPlane** fPlanes;     // [fNLayers] Shower Plane objects
 
   TClonesArray*  fTrackProj;    // projection of track onto plane