THcDriftChamberPlane

  1) Remove subtraction of hodoscope starttime from ProcessHits
    which is called by THcDC::Decode.
  The hodoscope start time is determined in hodoscope Decode
   which is called after THcDC::Decode. So drift chamber
   was using the starttime from the previous event.

   2) Create method SubtractStartTime which is called by THcDC:CoarseTrack

THcDriftChamber
    1) Move fChambers[ic]->ProcessHits() to CoarseTrack
    2) Before looping through fChambers[ic]->ProcessHits() in CoarseTrack
        loop through all planes and call fPlanes[ip]->SubtractStartTime()

Removed cuts in calculation of fXAtCer and fYAtCer

Add time windows for pos and neg ADC selection.

Added fGoodPosAdcMult and fGoodNegAdcMult variables

Remove cuts for  calculation of fXAtAero and fYAtAero

  Correction is done in the hitlist by comparing the trigger time
  from the TI to the trigger time in the FADC.  These generally have
  a fixed offset.  A FUDGE offset is included to make the difference
  between the FADC and TI times zero or close to zero.

  The FADC delivers times in units of 0.0625 ns.  The trigger time is in
  units of 4ns.  So the FADCs will be shifted by 64 units for each unit
  of trigger time slippage.

fFPTimeAll is the average focal plane time from all
  scintillator times. In EstimatedFocalPlaneTime fFPTimeAll is set to
  starttime which is the average FP time without using the DC track.

If there are tracks, then a focal plane time is calculated
for each track using the track information to improve the focal
plane time calculation and save for each track. Previously
 fFPTimeAll was being filled for each track, so in effect
fFPTimeAll was being filled with the FP time of the last track.
ALso there was no check that nFPTimeSum was not zero when
calculating average FP time, so one could get NaN values for
fFPTimeAll

To fix the problem , added  and modified line

Double_t fptime=-1000;
if (nFPTimeSum>0) fptime = FPTimeSum/nFPTimeSum;

And moved setting fFPTimeAll to THcHodoscope::FineProcess
Now fFPTimeAll is set the the average focal plane time of
the golden track.

  Newer root version
  Do make and scons
  Do ldd and version check after building
  Print python and scons versions

  Usual trigger time offsets for the 3 FADC crates are set
  in slipstudy.C

  This event type handler can be used to gather statistics on the slippage
  of trigger times in FADC250 headers.
  It looks at the relative trigger times between the ADCs and the TI modules
  and counts the number of times the relative times shift.

  The ADC trigger times will have a small offset relative to the TI trigger
  times.  The method AddExpectedOffset(roc, offset) allows the setting
  of this offset on a crate by crate basis.  If expected offsets are not
  set, the first event is used to determine the offsets.

  (Otherwise breaks with root 5)

  We'll start with version 0.90.0 and increment the second number up to 99.
  Goal is to declare a verison 1.0.0 sometime over the summer.

  To change the version, edit both SConscript.py and Makefile.

Modify THcDCTrack::TrackFit to calculate residuals for each plane
   without the plane involved in the track.

Add tree variable residualExclPlane

   THcFormula::DefinedValue is removed as
     THaFormula::DefinedValue now handles the kCutScaler and kCutNCalled types

      Includes THaFormula patches needed for THcFormula simplification

  Fix annoyance of spurious extra prompt on exit
  Add -l, -v and --version command line options

1) Add to tree variable arrays stub_x,stub_xp,stub_y and stub_yp
    and ncombos. These are arrays of the number of spacepoints in
    the chamber.

2) In ReadDatabase
   a) new variable  fRatio_xpfp_to_xfp which is set differently
       for HMS and SHMS. Used in method LeftRight
   b) Set default value of optional parameter fStubMaxXPDiff = 999.

3)  Modified method LeftRight
   a) Previously only for the old fHMSStyleChambers would
      the code only select LR combinations where
      the difference between stub_xp to space_point_X*ratio is
        with in the fStubMaxXPDiff .
      The ratio came from the HMS optics.
   b) Changed LeftRight so that if fStubMaxXPDiff < 100
         then the code will only select LR combinations where
      the difference between stub_xp to space_point_X*ratio is
        with in the fStubMaxXPDiff . The ratio is
      set in ReadDataBase according to the spectrometer.


5) In FindSpacePoints add some comments.

Modified THcSpacePoint.h to have public member that keeps track
   of the spacepoint ID number for each chamber.
   Used in THcDC::LinkStubs

Modified THcDC.cxx and h
1) Add tree variables sp1_id and sp2_id which are the integer
    ID in the spacepoints in each chamber used in the golden track.
     The ID refers to index in array of spacepoints created in
     THcDriftChamber.cxx

2) Modify LinkStub to fill sp1_id and sp2_id for each track

3) For golden track fill sp1_ID_best and sp2_ID_best

Add methods for setting the spacepoint integer identifier
  for two spacepoints used in the track.

Add calls:
+  Double_t GetStubX() {return fStub[0];};
+  Double_t GetStubXP() {return fStub[2];};
+  Double_t GetStubY() {return fStub[1];};
+  Double_t GetStubYP() {return fStub[3];};

  Just replaces Process method.

Add AdcTimeWindowMin and Max for each pmt

Add tree variables GoodPosAdcMult and fGoodPosAdcHitUsed

Fix sign mistake in correction for vertical pointing offset
   in THcExtTarCor.cxx

Write Hall C version of THcReactionPoint

Add new ntuple variables
goodAdcMult == the number of hits in raw ADc for that channel
goodAdcHitUsed == which hit number of the raw ADC hits passed the
                      the ADC time window cut

Modified fAdcTimeWindowMin  and fAdcTimeWindowMax to be arrays
  that are the size of the number of Cerenkov pmts.

  Statistics for missing ADC reference times are still accumulated
  and still printed at the end of the run and stored in parameters.

Moved the old track efficiency tests to OriginalTrackEffTest method and
this is presently commented out. Eventually can be eliminated.

The old track efficiency assumed that both planes had matching
  paddlle spacing which is not true for SHMS.

Added new method TrackEffTest which calculates the number of clusters
  in each plane and the position for the cluster for scintillator
  paddle that had good times in both ends. Selects events with
  one cluster with size of 1 or 2 in either 4 of 4 planes
  or 3 of 4 planes. Also checks that the events are inside
  a given range of scintillator paddles.

  So that the counts can be reported in report files.

   First hit above this cut is taken as the reference time.
   Cuts are per detector and read as an optional parameter in
   each detector class.  TDCs and Flash ADC times have separate cuts.

Implement a second method of selecting best reference time.
  If InitHitList is passed a possitive reference time cut,
    use the first reference time above that cut.  If none of the
    hits pass the cut, then no reference time is set
  If InitHitList is passed a negative cut, the cut is taken to be
    the absolute value of the cut, but a reference time is garuanteed
    to be set as long as there is at least one reference time hit.
    If none of the reference time hits pass the cut, then the last
    hit is taken to be the reference time.  (Should be the largest.)

Allow reference time cuts for all the spectrometer detector classes
  Hodoscopes, drift chambers, Aerogels, Gas Cherenkovs, Calorimeters

  Add End method to each detector that calls MissReport which
  prints the number of events with missing TDC and ADC reference times.

  in THcDummySpectrometer and THcTrigApp
  by removing calls to DefineVarsFromList that are not needed because
  those classes define no variables.