From 232d3c5522d58cae9e64618798478d7953cafdd7 Mon Sep 17 00:00:00 2001
From: Mark Jones <jones@jlab.org>
Date: Thu, 20 Jun 2013 13:30:04 -0400
Subject: [PATCH] 1) Add flag fDebugDC if set=1 in THcDC::Setup will write out
info
---
src/THcDC.cxx | 11 ++++++++---
1 file changed, 8 insertions(+), 3 deletions(-)
diff --git a/src/THcDC.cxx b/src/THcDC.cxx
index f2c5f3e..7786201 100644
--- a/src/THcDC.cxx
+++ b/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];
--
GitLab