Set 2d histogram limits on tracking results.

   Add track coordinates to track structure.
   Add code for single stubs
   Try to deal with plane number/index confusion

Since we have not yet overwritten the = operator for space points
had to add code to do this.

  Uses ROOT matrix math instead of solve_four_by_four
  Nothing is done with found tracks yet.  Will need to stuff these into
     podd style tracks.  (Or a class inherited from podd style tracks.)
  Residuals not saved anywhere
  Single stub mode not yet implemented
  No printing out of trackfit results (hdebugtrackprint)

Set fDebugDriftCh in DefineVariables
Put if fDebugDriftCh in front of all cout statements
In THcDriftChamber::DestroyPoorSpacePoints replace line (*fSpacePoints)[isp] = (*fSpacePoints)[osp];
In THcDriftChamber::SelectSpacePoints() replace line (*fSpacePoints)[sp_count] = (*fSpacePoints)[isp];

Set fDebugDC in THcDC::ReadDatabase

   for first hit. Added an else branch to increase finalnum
2)  In THcDriftChamber::LeftRight() when using the small angle
    approximation then assignment of sign to the y planes was reversed.

  Compiles and doesn't crash.

  THcDriftChamber::SpacePoint -> THcSpacePoint
  That way space points can be put in a TClonesArray and be made known
  to THcDC when stub information is added.

  1)   Moved hardcoding of Yplane plane number to Addplane since it is called
     before ReadDatabase
  2) In FindEasySpacePoint add variable max_dist to match the code in
     h_find_easy_space_point.f more and fix problem
  3) As in h_find_easy_space_point  use max_dist instead
     of FSpacePointCriterion in IF statement to pick matching "x" hits.
In THCDC.h
  1) Set GetPicth and GetCentralWire as Double_t
  2) Set  GetSpacePointCriterion to return fSpace_Point_Criterion2[chamber-1]
     instead of sqrt

The arrays are  frPosTDCHits, frNegTDCHits,frPosADCHits and frNegADCHits
There are filled in THcScintillatorPlane:ProcessHits
2) Also modified THcScintillatorPlane::ProcessHits explicitly set variable
padnum=hit->fCounter and index=padnum-1 so the the TClone arrays for
scintillator variable set padnum=1,16 while the index=0-15
Want pdnum=1-16 because it is used in a hitogram fo counts versus
paddle number
3) Because of change to filling PaddleNumber in ProcessHits had
to modify THcScintillatorPlane::PulseHeightCorrection so the
index variable "j" is set to etPaddleNumber()-1 .

In the ENGINE, this can happen if a hit belongs to more than one space point.

  New class THcDC is whole set of wire chambers.
  THcDriftChamber is one chamber (6 planes)
  Implement space point finding from h_pattern_recognition up to point
     before where space point lists from the two chambers are combined.
  Includes
      FindEasySpacePoint   (h_find_easy_space_point)
      FindHardSpacePoints  (find_space_points)
      DestroyPoorSpacePoints (h_sp_destroy)
      SpacePointMultiwire  (h_sp_multiwire)
      ChooseSingleHit 	   (h_choose_single_hit)
      SelectSpacePoint	   (select_space_points)

Space points are arrays of structures in the THcDriftChamber class.  Will
probably need to take space points out class so that they can be looked at
at the THcDC level.

New code is only tested to see that it compiles and that space points
are sometimes found.

properly with 4. Before it was filled with 1,2,3,4 for each plane 1,2,3,4
2) Only if fTofUsingInvAdc==1 the read database for hodo_pos_invadc_offset, hodo_pos_invadc_linear
,hodo_pos_invadc_adc, hodo_neg_invadc_offset, hodo_neg_invadc_linear
,hodo_neg_invadc_adc