From 5d3d0ef277a584409f6207fccf18ebc804e894c9 Mon Sep 17 00:00:00 2001
From: "Stephen A. Wood" <saw@jlab.org>
Date: Wed, 21 Aug 2013 16:12:45 -0400
Subject: [PATCH] Tweak copying of tracks fitting results into THaTrack
Results are focal plane coordinate (transport) system, not detector system
Turn on golden track selection.
---
src/THcDC.cxx | 7 ++++---
src/THcHallCSpectrometer.cxx | 2 +-
2 files changed, 5 insertions(+), 4 deletions(-)
diff --git a/src/THcDC.cxx b/src/THcDC.cxx
index d596df0..4dc3af6 100644
--- a/src/THcDC.cxx
+++ b/src/THcDC.cxx
@@ -481,14 +481,15 @@ Int_t THcDC::CoarseTrack( TClonesArray& tracks )
for(Int_t itrack=0;itrack<fNDCTracks;itrack++) {
THaTrack* theTrack = NULL;
theTrack = AddTrack(tracks, 0.0, 0.0, 0.0, 0.0); // Leaving off trackID
- // Should we add stubs with AddCluster?
+ // Should we add stubs with AddCluster? Could we do this
+ // by having stubs inherit from cluster
THcDCTrack *tr = static_cast<THcDCTrack*>( fDCTracks->At(itrack));
- theTrack->SetD(tr->GetX(), tr->GetY(), tr->GetXP(), tr->GetYP());
+ theTrack->Set(tr->GetX(), tr->GetY(), tr->GetXP(), tr->GetYP());
theTrack->SetFlag((UInt_t) 0);
Int_t nhits=tr->GetNHits();
// Need to look at how engine does chi2 and track selection. Reduced?
- theTrack->SetChi2(tr->GetChisq(),nhits-4); // Nconstraints - Nparameters
+ theTrack->SetChi2(tr->GetChisq(),tr->GetNFree());
// CalcFocalPlaneCoords. Aren't our tracks already in focal plane coords
// We should have some kind of track ID so that the THaTrack can be
// associate back with the DC track
diff --git a/src/THcHallCSpectrometer.cxx b/src/THcHallCSpectrometer.cxx
index 53c141c..a23020b 100644
--- a/src/THcHallCSpectrometer.cxx
+++ b/src/THcHallCSpectrometer.cxx
@@ -68,7 +68,7 @@ THcHallCSpectrometer::THcHallCSpectrometer( const char* name, const char* descri
//sc_ref = static_cast<THaScintillator*>(GetDetector("s1"));
- SetTrSorting(kFALSE);
+ SetTrSorting(kTRUE);
}
//_____________________________________________________________________________
--
GitLab