Properly zero out traceback coefficients and exponents.
  This insures that 5th exponent is zero if not specified on recon file
Golden track using prune:
      Selection of golden track using prune method is added.
      A bug is also fixed in THcHodoscope class
      Number of pmts hits, focal plane time, good time for plane 4
      and good time for plane 3 are set to the tracks in
      THcHodoscope class.
Supply SOS particle mass for run 50017 in test.database
Remove old source file from examples
Add DefineVariables to put beta and betachi2 into tree
Replace fixed length arrays with vectors and vectors of structures
    Structure TOFPInfo holds information for hits within a plane
    Structure TOFCalc holds hit information for all the hits.
    fdEdX changed from 2d array to vector of vectors
Remove obsoleted/commented out code in THcHodoscope
Convert some arrays to vectors and move from class to method
  Convert some arrays that were fixed length to vectors.
  Since these arrays are only used local to the method
  ThcHallCSpectrometer::TrackCalc(), they were moved out of the class
  into the method.
  The affected arrays are fX2D, fY2D, f2XHits and f2YHits
Convert Reconstruction coefs/exps to vector of structures
Add access to the golden track proper in THaGolden Track class.
Access the energy of the golden track from THaGoldenTrack.

    Include THaTrack.h in THaGoldenTrack.h.
    Add GetEnergy method in THaGoldenTrack.h.
    Define variable name "e" for the energy of golden track in THaGoldenTrack.C.
Make association of hit clusters to tracks in the HMS more
consistent with Engine.

	   Reverse order of cycling through clusters in THcShower::MatchCluster.
Add a script for event-by-event comparison of energies of the best
track in Engine and the golden track in hcana.

      Add directory examples/gtrk_e holding the code.
Remove analysis of the maximum cluster and cluster association to the
first track from the HMS shower counter code.

    Remove variables and code related to the maximum cluster of hits,
    and association of a cluster to the 1-st spectrometer track in
    THcShower.h and THcShower.cxx. Those were used for testing of the
    code, not needed any more.
Add float precision calculation of the X coordinate of hit clusters in
the HMS calorimeter.

    Add src/THcShowerCluster.h.float, a version of THcShowerCluster
    class with float precision calculation in clX method.
Remove softlinks to root files in examples/gtrk_e
Edit debug outputs in the Shower code.

     Edit debug outputs in THcSHower and THcShowerPlane classes. Group
     outputs preferably by the end of methods, for smooth code
     flow. Structure tests of outputs for easy reading.
Remove example/gtrk_e as scripts have been moved to hcana_replay
Move THcShowerCluster.h and THcShowerHit.h into THcShower.h
  This is because the classes in these header files are only
  used by THcShower.
Tracking corrections to agree with ENGINE.

First, in h_pattern_recognition.f we have:

https://github.com/zahmed9/engine/blob/mkj/HTRACKING/h_pattern_recognition.f#L139
(hdc_wire_center(yy)-hdc_wire_center(yyprime))**2 .lt. (hspace_point_criterion(ich))

So I changed
https://github.com/zahmed9/hcana/blob/develop/src/THcDriftChamber.cxx#L292
TMath::Abs(fHits[plane_hitind]->GetPos() - fHits[planep_hitind]->GetPos()) < fSpacePointCriterion

to

pow((fHits[plane_hitind]->GetPos() - fHits[planep_hitind]->GetPos()),2) < fSpacePointCriterion

Second, the value of MAX_SPACE_POINTS in ENGINE is 100 so in hcana I changed it from 50 to 100.
Fixed problems with DefineVariables in THcHallCSpectrometer