diff --git a/examples/tracking.C b/examples/tracking.C
new file mode 100644
index 0000000000000000000000000000000000000000..74d40e286aaa0eba6fc3cb8047e3ce45ed7dd5d6
--- /dev/null
+++ b/examples/tracking.C
@@ -0,0 +1,19 @@
+{
+  TFile* f = new TFile("hodtest.root");
+  TTree *T=(TTree*)f->Get("T");
+  TCanvas *c1 = new TCanvas("c1", "Drift Tracking", 800, 1000); 
+  c1->Divide(2,3);
+
+  c1->cd(1);
+  T->Draw("H.dc.nhit");
+  c1->cd(2);
+  T->Draw("H.dc.ntrack");
+  c1->cd(3);
+  T->Draw("H.dc.x:H.dc.y");
+  c1->cd(4);
+  T->Draw("H.dc.x:H.dc.xp");
+  c1->cd(5);
+  T->Draw("H.dc.y:H.dc.yp");
+  c1->cd(6);
+  T->Draw("H.dc.xp:H.dc.yp");
+}
diff --git a/src/THcDC.cxx b/src/THcDC.cxx
index 812bdcf855a06325f37800b75f5c1076cace1831..e2804ff306e2e7f6e7bb5037b74da2052d55f479 100644
--- a/src/THcDC.cxx
+++ b/src/THcDC.cxx
@@ -328,6 +328,11 @@ Int_t THcDC::DefineVariables( EMode mode )
 
   RVarDef vars[] = {
     { "nhit", "Number of DC hits",  "fNhits" },
+    { "ntrack", "Number of Tracks", "fNDCTracks" },
+    { "x", "X at focal plane", "fDCTracks.THcDCTrack.GetX()"},
+    { "y", "Y at focal plane", "fDCTracks.THcDCTrack.GetY()"},
+    { "xp", "YP at focal plane", "fDCTracks.THcDCTrack.GetXP()"},
+    { "yp", "YP at focal plane", "fDCTracks.THcDCTrack.GetYP()"},
     { 0 }
   };
   return DefineVarsFromList( vars, mode );
diff --git a/src/THcDCTrack.h b/src/THcDCTrack.h
index 2b67faf1aacf7a94090d1ed975911859068dd6c7..a44b9494e01e9bec43148c764facb09c152e975e 100644
--- a/src/THcDCTrack.h
+++ b/src/THcDCTrack.h
@@ -35,6 +35,10 @@ public:
   Double_t GetCoord(Int_t ip)     const {return fCoords[ip];}
   Double_t GetResidual(Int_t ip)     const {return fResiduals[ip];}
   void GetRay(Double_t *ray) const {ray[0]=fX_fp; ray[1]=fY_fp; ray[2]=fXp_fp; ray[3]=fYp_fp;}
+  Double_t GetX()                 const {return fX_fp;}
+  Double_t GetY()                 const {return fY_fp;}
+  Double_t GetXP()                 const {return fXp_fp;}
+  Double_t GetYP()                 const {return fYp_fp;}
   void SetNFree(Int_t nfree)           {fNfree = nfree;}
   void SetCoord(Int_t ip, Double_t coord) {fCoords[ip] = coord;}
   void SetResidual(Int_t ip, Double_t coord) {fResiduals[ip] = coord;}