1) Add flag fDebugDC if set=1 in THcDC::Setup will write out info

src/THcDC.cxx

@@ -126,6 +126,7 @@ void THcDC::Setup(const char* name, const char* description)
     }
     fPlanes.push_back(newplane);
     newplane->SetDebug(fDebug);
+    fDebugDC=0;
     cout << "Created Drift Chamber Plane " << fPlaneNames[i] << ", " << desc << endl;
 
   }
@@ -399,7 +400,8 @@ Int_t THcDC::Decode( const THaEvData& evdata )
 {
 
   ClearEvent();
-
+  Int_t num_event = evdata.GetEvNum();
+  if (fDebugDC) cout << " event num = " << num_event << endl;
   // Get the Hall C style hitlist (fRawHitList) for this event
   fNhits = THcHitList::DecodeToHitList(evdata);
 
@@ -447,7 +449,6 @@ Int_t THcDC::CoarseTrack( TClonesArray& /* tracks */ )
   //
   //  static const Double_t sqrt2 = TMath::Sqrt(2.);
   for(Int_t i=0;i<fNChambers;i++) {
-
     fChambers[i]->FindSpacePoints();
     fChambers[i]->CorrectHitTimes();
     fChambers[i]->LeftRight();
@@ -500,18 +501,20 @@ void THcDC::LinkStubs()
     Int_t nchamber=fChambers[ich]->GetChamberNum();
     TClonesArray* spacepointarray = fChambers[ich]->GetSpacePointsP();
     for(Int_t isp=0;isp<fChambers[ich]->GetNSpacePoints();isp++) {
+      if (fDebugDC) cout << " Chamber = " << nchamber << " number of space pts = " << fNSp << endl;
       fSp.push_back(static_cast<THcSpacePoint*>(spacepointarray->At(isp)));
       fSp[fNSp]->fNChamber = nchamber;
       fNSp++;
     }
   }
+  //  fDebugDC=0;
   Int_t ntracks_fp=0;		// Number of Focal Plane tracks found
   Double_t stubminx = 999999;
   Double_t stubminy = 999999;
   Double_t stubminxp = 999999;
   Double_t stubminyp = 999999;
   Int_t stub_tracks[MAXTRACKS];
-
+  if (fDebugDC) cout << " single stub , fsnp" << fSingleStub << fNSp << endl;
   if(!fSingleStub) {
     for(Int_t isp1=0;isp1<fNSp-1;isp1++) {
       Int_t sptracks=0;
@@ -552,6 +555,7 @@ void THcDC::LinkStubs()
 	       && (TMath::Abs(dposxp) < fXptTrCriterion)
 	       && (TMath::Abs(dposyp) < fYptTrCriterion)) {
 	      if(newtrack) {
+  if (fDebugDC) cout << " new track" << endl;
 		assert(sptracks==0);
 		//stubtest=1;  Used in h_track_tests.f
 		// Make a new track if there are not to many
@@ -575,6 +579,7 @@ void THcDC::LinkStubs()
 		  return;
 		}
 	      } else {
+  if (fDebugDC) cout << " check if another space point in same chamber" << endl;
 		// Check if there is another space point in the same chamber
 		for(Int_t itrack=0;itrack<sptracks;itrack++) {
 		  Int_t track=stub_tracks[itrack];