Use pointer to event from evdata->GetRawDataBuffer() instead
  of making a local copy of the event before scanning it for
  scaler data

  There was already an AddEvtType in the base class.
  Should use SetEvtType to set the event type.  If the handler is to
  look for more than one event type, use AddEvtType for the extra
  event types.

**This class calculates kinematic variables related to the secondary (hadron)
detected particle

When making changes to support THcPrimaryKine added line
track->SetEnergy
which set the track energy to energy determined from
measurement of spectrometer delta.
But this energy should be the energy measured by the
Shower calorimeter.
Deleted the line in THcHallCSpectrometer since it is not
needed anyway.

Add class THcPrimaryKine to calculate physics kinematics quantities for single arm

Modified  THcHallCSpectrometer so that it would fill track energy and momentum 3 vector

   Gets rid of vars.push_back(RVarDef ... stuff

Also do some whitespace cleanup

  This is so that hodoscope paddles in Eve are illuminated

the best track.

-Changed the name of a Good tdc variables
-Added all good adc and tdc variables
modified:   src/THcScintillatorPlane.cxx
modified:   src/THcScintillatorPlane.h

example:
-fTotNumPosAdcHits (raw hodo+ multiplicity)
-fTotNumNegAdcHits (raw hodo- multiplicity)
-fTotNumAdcHits (raw hodo multiplicity)

Indentation cleanup

   Added to THcScintillatorPlane.h
   Initialize and clear in THcScintillatorPlane.cxx
Add raw multiplicities for ADCs and TDCs

  With this commit, db_run.dat is not needed in hcana
  If THcRunParms is used, hcana will get the beam momentum from the
  parameter "gpbeam" instead of db_run.dat.  To use, use the
  THaRun::SetRunParamClass method to set the run parameter class
  to "THcRunParameters".

  THaRun* run = new THaRun(RunFileName);
  run->SetRunParamClass("THcRunParameters");

autogeneration of src/HallC_LinkDef.h - all using SCons

Accidentally deleted line
btdcraw_pos = kTRUE;

in previouse commit

-Added fAdcTimeWindow Min/Max to THcHodoscope.cxx (.h)
-Added fAdcTimeWindow Min/Max to THcScintillatorPlane.cxx (.h)
-Initialized loop over number of hits inside the FillDynamic_Pedestal() method
-added some variables in the THcScintillatorPlane.cxx Dynamic Pedestal Method
-Added loop over adc_pos and adc_neg hits, and required the hits to fall
within fAdcTimeWindow Min/Max
-Added the requirement of at least a valid (adc && tdc) on either end of a bar
-re-defined boolean operator for adc_pos/adc_neg to be True if  they passed the errorflag and timeCUT requirements in the fADC_Dynamic Pedestal method.
-- deleted adc_pos/neg initial definitions, and set their initial values to -999.

  Add "Good" variables to calorimeter class
     occupancies and multiplicities
     required to pass time and threshold cut
  Change some variables from arrays to vectors
  Modify how thresholds are calculated
  Remove obsolete variables: posadchits, negadchits
  Add flash adc debug flag
  Correctly delete pointers in destructor

  Add "Good" variables to calorimeter class
     occupancies and multiplicities
     required to pass time and threshold cut
  Change some variables from arrays to vectors
  Modify how thresholds are calculated
  Remove obsolete variables: posadchits, negadchits
  Add flash adc debug flag
  Correctly delete pointers in destructor

Make matched cluster variables of the array part of calorimeter private.

Initialize fY variable in the THcShowerHit::THcShowerHit default
constructor.

Correct typo in THcShowerArray::Clear, assign -1000 to fMatchClY.

Edit debug output in THcShowerArray::CoarseProcess and remove obsolete
comments THcShowerArray::MatchCluster.

Clear fSizeClustArray and fNblockHighEnergy quantitites in
THcShower::Clear, for just in case.

So files don't split.
Seems to be working, hcana has not crashed yet with file sizes
upto 2.6 GB.

Main purpose is to calculate the DC efficiency per wire
using the track information of the golden track.


THcDC.h
-----------
1) add method EfficiciencyPerWire
2) Add variables fWire_hit_did and fWire_hit_should

THcDC.cxx
------------
1) fWire_hit_did and fWire_hit_should are fixed length arrays
   the size of NPlanes that in set in method Init
2) Add variables wireHitDid and wireHitShould to the tree
3) In method SetFocalPlaneBestTrack add call to EfficiencyPerWire method
4) added EfficiencyPerWire method which loops over the hits
    and compares the wire number for the hit to the the track
    position at the plane (convert back to wire number).
    Then fills fWire_hit_did with the wire number if
    equal to the wire position of the hit equals the track.
    Also does loop over all planes and fills fWire_hit_should
    with the wire number of the track in each plane.

THcDriftChamberPlane.h and cxx
---------------------
1) Added method CalcWireFromPos

The purpose is to add a call in THcHallCSpectrometer:TrackCalc
to fill variables in THcDC for the golden track.

THcHallCSpectrometer.h
----------------------
1) Add object THcDC* fDC
2) Eliminate object fShower since it was not being used.

THcHallCSpectrometer.cxx
----------------------
1) in ReadDatabase method cast fDC in detector "dc" is defined
2) in FindVertices method add call
    to fDC->SetFocalPlaneBestTrack(hit_gold_track)
    where hit_gold_track is the HaTracks array index
    of the golden track

THcDC.h
-------
1) Add method SetFocalPlaneBestTrack

THcDC.cxx
---------
1) Eliminate filling of best focal plane quantities
   in CoarseTrack which defined "best" by lowest chi2
2) Create method SetFocalPlaneBestTrack
   a)fills the best focal plane quantities using the golden track
   b) fills fResiduals using golden track

THcShowerHit.h and THcShowerHit.cxx
---------------
1) Added fY for horizontal position
2) Added fY to THcShowerHit

THcShower.h
-----------
1) Added variables
2) Added method addY, clY to get the cluster Y position

THcShower.cxx
--------------
1) Added variables to the tree for "layer" detector
    either SHMS preshower or HMS layers
2) Added variables if the SHMS "array" is used
3) Added y position to THcShowerHit in CoarseProcess
4) In FineProcess, if SHMS "array" used then fill X,Y info
     for the Cluster and the "best" track that matched the cluster

THcShowerArray.h and THcShowerArray.cxx
---------------------------------------
1) Added methods to Get the array cluster size, X and Y position
   and block number that has the highest energy in the cluster
2) Added array fBlockClusterID for all blocks. initialized to -1
   for each block for each event and then filled with
   the cluster number if the block is included in the cluster
3) Add variable fNhits and fNgoodhits

Main purpose was to have option for per wire tzero offsets

Modified THcDriftChamberPlane.h
------------------------------
1) added fUsingTzeroPerWire as parameter flag to switch on per wire tzero offsets when parameter equals 1
2) fTzeroWire is array of tzero offsets for each wire in plane

Modified THcDriftChamberPlane.cxx
------------------------------
1) THcDriftChamberPlane::ReadDatabase
    a)set a default value of fUsingTzeroPerWire=0
    b)added to DBRequest list option to read in fUsingTzeroPerWire
    c) if fUsingTzeroPerWire=1 then separate DBRequest to read in fTzeroWire
   d) if fUsingTzeroPerWire=0 then all fTzeroWire=0

2) THcDriftChamberPlane::ProcessHits
------------------------------------
   a) Modify calcualation of drift time to subtract fTzeroWire
Double_t time = -StartTime- rawtdc*fNSperChan + fPlaneTimeZero - fTzeroWire[wireNum-1];

  Then the # tubes of the Aerogel is known.
  Remove MaxNumXXXAeroPmt constants from THcAerogel.h

The line fNpeSum += fNegNpeSum
was basically double counting the fNegNpeSum
Same for fNpeSum += fPosNpeSum
since they were in a loop over Neg and Pos hits.

Just set
fNpeSum = fNegNpeSum + fPosNpeSum
after the two loops over Neg and Pos hits

THcShowerArray
1) in ProcessHits methodfFixed mistake in index of fThresh
   was using padnum instead of padnum-1
2) add parameter pcal_arr_AdcThreshold so that threshold
   can be set by a parameter. Default is zero which gives
   the best resolution.

THcShowerPlane
1) Added the parameter cal_AdcNegThreshold and cal_AdcPosThreshold
   so that the threshold can be set through a parameter
   Previously it was hardcoded.