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

Fix mistake in THcShowerArray and make the threshold a parameter

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.

Fix bug for corner case for gcc 4.4 and ROOT5

New parameters are added.  hallc_replay of 3/31/2017 or later is needed.

As of now, the size of the TClonesArray objects are hard-coded in the
respective header files. This is not okay in the long run. Ideally,
the size of the TClonesArrays should be defined by the parameter which
dictates how many readout detectors exist for the respective
detector. I tried doing this in the Init method however, podd
complained regarding the size of the TClonesArrays in the DEF-file
histograms. I was unsuccessful in my first attempt. This should be
addressed in the future even though things work for now.

  Process slot with clock first so that time is available to other slots
  Add option to use data in first scaler event
     (handler->SetUseFirstEvent(Bool_t)) instead of subtracting first event
     from subsequent events

Modify THcHodoscope.cxx

The change was in part of code which calculates the fiducial
 area of the scintillators where the drift chamber should
 have had a track. This sets fGoodScintHits=1 or 0
 Used to calculate the DC track efficiency

1) Eliminate variables fScinShould and fScindid. These will
   be defined in the cuts file.
2) Eliminate code which determined bestXpScin and bestYpScin since
   since these variables were never used.
3) Initialized  fHitSweet1X,fHitSweet1Y,fHitSweet2X and
   fHitSweet2Y to zero.
4) Change so that loops over number of paddles do not
    assume that X (Y) planes for Ch1 and Ch2 have same number
    of paddles.

Modify THcShower, THcShowerPlane and THcShowerArray classes

Main purposes:
1) Add the possibility of multiple ADC hits per channel and
   select best ADC by a time window.
2) Change ProcessHits method in ShowerPlane and ShowerArray to only
    fill "raw" TClonesArrays.
3) Add CoarseProcessHits method to ShowerPlane and ShowerArray
    to fill the "good" data arrays from call in THcShower:CoarseProcess
4) Moved code to fill Track object with energy from FineProcess
   to CoarseProcess so that the Track->Get(Energy) can be used by other
    detectors in their FineProcess.
5) FineProcess loops over Track object add gets Track energy
6) Add AdcErrorFlag TClonesArray to be filled to check for
   problems with FADC




THcShower.h
1) Add methods GetADCMode(),GetAdcTimeWindowMin(),GetAdcTimeWindowMax()
    so that the ShowerPlane can access these parameters
2) Add parameters ADCMode,fAdcTimeWindowMin,fAdcTimeWindowMax
3) Add variables
   fNclustTrack : NUmber of cluster that matches best track
   fXclustTrack : X pos of cluster that matches best track
   fXTrack      : X pos of best track that matches cluster
   fEtrack  : Cluster energy associated with best track
   fEtrackNorm  : Cluster energy/Track momentum associated with best track
THcShower.cxx
1) Add calls to THcShowerPlane and THcShowerArray CoarseProcessHits
    in CoarseProcess
2) Moved filling of Energy in Track object to CoarseProcess
3) In FineProcess get the "best" Track object and match
   Track energy and trajectory to cluster.

THcShowerArray.h and cxx
1)add parameters fADCMode, fAdcTimeWindowMin and fAdcTimeWindowMax
2) Add method CoarseProcessHits
3) In ProcessHits fill fPosThresh and fNegThresh with a fix
    value of 250 integrated channels above pedestal.
     Need to make this a parameter.

THcShowerPlane.h and cxx
1)Gets fADCMode, fAdcTimeWindowMin and fAdcTimeWindowMax from THcShower
2) Add method CoarseProcessHits
3) In ProcessHits fill fPosThresh and fNegThresh with a fix
    value of 250 integrated channels above pedestal.
     Need to make this a parameter.

Add method to THcRawAdcHit

Add method
GetF250_PeakPedestalRatio() {return fPeakPedestalRatio;};

so that a detector class can calculate the pedestal that was used
in determining the pedestal subtracted FADC Pulse Integral.

This is useful if the detector class wants to set a threshold
above the pedestal.

Modified THcHallCSpectrometer.cxx and THcHodoscope.cxx

The change was made so that NonTrackingDetector could have information
on the "best" track in their FineProcess method.

In THcHallCSpectrometer::FindVertices the "best" track will be selected
and the Track object will its SetIndex method to set the index to 0 for
the best track. The other track will have index=1 or greater.

In the NonTrackingDetector FineProcess, the track that is selected to
be the "best" can be found by looping through the tracks checking
track index using Track.GetIndex() for find the track with Index=0.

In  the NonTrackingDetector's CoarseProcess they need to set quantities
for the Track object needed for selection of "best" track for the
"prune" and "usingscin" track selection methods.

1) In THcHallCSpectrometer:
  a) Moved selection methods of "best" track from TrackCalc to FindVertices.
  b) Removed setting of Golden Track in the selection methods.
  c) Add setting of Track Index=0 for "best" track in each selection method.
  d) Modified TrackCalc method to loop through Track object and find
     the track with Index=0 and set the Golden Track equal to this Track.

2) In THcHodoscope.cxx
  a) Moved all calculations in FineProcess into CoarseProcess
  b) CoarseProcess calculates the scintillator Beta for each Track
     and does Track->SetBeta.
  c) In determination "zcor" for corrected hits times
      use fBetaNominal instead fBetap from track. Still use Track
      to get distance.
  c) Eliminate fBetap variable which was using the Track momentum
  d) Add variable fBeta which is just from scintillator times.
  e) In FineProcess fill fBeta with the the beta associated
     with the "best" track determined by selecting Track with index=0

Modified HcTrigDet::Init , Decode and ReadDatabase

In HcTrigDet::Init moved InitHitList after gHcDetectorMap->FillMap so that the refindex can be filled for the trigger TDC modules
In ReadDatabase add parameters fTdcOffset and fTdcChanperNS
In Decode set fTdcTime[cnt] = rawTdcHit.GetTime()*fTdcChanperNS-fTdcOffset
Ths makes the time into a positive number wiht units of ns
Need to add the parameter to either thms.param or tshms.param or
it will use defaults of 0.1 and 300.
Need to eventually get the fTdcChanperNS from event 125.

	In THcAerogel.cxx, mean pedestals are calculated if
	fMinPeds>0. Set fMinPeds to default zero in
	THcAerogel::InitializePedestals. Reinitialize it with
	paero_min_peds parameter in paero.param file.

Modify THcDC::CoarseTrack and THcDC::TrackFit

Add variables fX_fp_best,fY_fp_best,fXp_fp_best,fYp_fp_best to THcDC.h
Modified THcDC::CoarseTrack so that the variables are filled with the
quantities from the track with the smallest chi-squared.
Modified THcDC::TrackFit so that the fResidual per plane is filled with
for the track with the lowest chi-squared.
Modified THcDC::DefineVariables so that access the variables as
dc.x_fp,dc.y_fp,dc.xp_fp and dc.yp_fp

Modified THcCherenkov.cxx and THcCherenkov.h

Add array fADC_hit that is a flag if there is a hit in the ADC
fADC_hit is set in Decode method
The FADC has a threshold on accepting ADC hit
Modified CoarseProcess method to eliminate old FASTBUS bias in the analysis

Modify ShowerPlane and ShowerArray ProcessHits Methods

Remove the hc_hodo_calib directory

The hc_hodo_calib directory is now under hallc_replay/CALIBRATION

Modified hodoscope and scintillatorplane code for hodo calibration

Temporarily set THcHodoscope.cxx so that only the first three planes are used to calculate the focal
plane time and beta.
In THcHodoscope , Modified format for writing out the dump file used in hodo calibration.
Add cosmicflag to THcScintillatorPLane so the corrected time is done right for cosmics.
Modified hc_hodo_calib/tofcal.f so that it uses detector 10 as the reference
and that it writes out better defualts for parameters in there are not enough events for a PMT.

Modify THcHodoHit.h and THcScintillatorPlane.cxx

THcHodoHit.h add methods to set and get the ADC peak amp and time
THcScintillatorPlane::ProcessHits changed to set the ADC peak amp and time in the Hodohit
DefineVariable added variables to tree

   Used by hodtest.C example (for old HMS data)

Modified THcHodoscope::FineProcess for dumping data for hodo calibrations

Moved the writing out of the "0" to the fDumpOut for hodo calibrations
to immediatiately after the TOF with track info calculations.
Before it had been after then lines

  if ( !fChern || !fShower ) {
    return 0;
  }

so if Cerenkov and shower were not used in replay the code would not get
to write out the "0" which is used to denote the end of an event data for
the calibration script.

Added capability to set F250 parameters from PSE125.

PSE125 holds information about setting of F250 modules. Now this
information is used to set `fNPedestalSamples` and `fNPeakSamples`
parameters of `THcRawAdcHit`.

This is done through `THcHitList`, since it has information about
which crate and module each channel is in. Currently, the information
is set for each event, because of the hitlist and analyzer design.

I added the neccessary support for these changes to all raw hit
classes.