Goldent track is moved to THcHallCSpectrometer::TrackCalc()
  if  fSelUsingScin == 0 then golden track is calculated just
  like podd. i.e. it is the first track with minimum chi2/ndf
  with sorting ON

  if fSelUsingScin == 1 then golden track is calculetd just like
  engine/HTRACKING/h_select_best_track_using_scin.h. This method
  gives the best track with minimum value of chi2/ndf but with
  additional cuts on the tracks. These cuts are on dedx, beta
  and on energy.

  This commit also corrects and modifies the beta code written
  in the THcHodoscope::FineProcess

  At this point we have godlend track and in additoion to
  previous vlues the following values are set to a track:
  Energy, beta, chi2 of beta and dedx[track][first hit].

NOTE: For some tracks we have large differences in value of chi
      square between ENGINE and hcana and this is already
      discussed in the last group meeting. Following emails are
      exchanged regarding this difference of chisquare.

email to Mark:

Hi Mark,
I tried to trace the origin of difference of the chi square value between
ENGINE and hcana. Please see the attached file, it has long table which I
cannot present in this email. Or you can see the file at the following
link, thanks.

https://userweb.jlab.org/~zahmed/chi2diff.txt

Ahmed

Mark's reply:
Ahmed,
        Thanks. I agree with what you see, but there is also another reason
for the difference. For the space point in ch2 the engine and hcana
choose different left/right combinations. That is why you see the large
jump in the summed chi2 from planes 9 and 10. The two codes pick
different left/right because there is a cut on the difference between the
measured stub xp and an expected xp = sp->GetX()/875  and one code has a combination
slightly below the cut value and the other slightly above. I do not see a way
to make the codes agree more exactly since the stub fitting involves a matrix
inversion.

                Cheers,
                       Mark

Goldent track is moved to THcHallCSpectrometer::TrackCalc()
  if  fSelUsingScin == 0 then golden track is calculated just
  like podd. i.e. it is the first track with minimum chi2/ndf
  with sorting ON

  if fSelUsingScin == 1 then golden track is calculetd just like
  engine/HTRACKING/h_select_best_track_using_scin.h. This method
  gives the best track with minimum value of chi2/ndf but with
  additional cuts on the tracks. These cuts are on dedx, beta
  and on energy.

  This commit also corrects and modifies the beta code written
  in the THcHodoscope::FineProcess

  At this point we have godlend track and in additoion to
  previous vlues the following values are set to a track:
  Energy, beta, chi2 of beta and dedx[track][first hit].

NOTE: For some tracks we have large differences in value of chi
      square between ENGINE and hcana and this is already
      discussed in the last group meeting. Following emails are
      exchanged regarding this difference of chisquare.

email to Mark:

Hi Mark,
I tried to trace the origin of difference of the chi square value between
ENGINE and hcana. Please see the attached file, it has long table which I
cannot present in this email. Or you can see the file at the following
link, thanks.

https://userweb.jlab.org/~zahmed/chi2diff.txt

Ahmed

Mark's reply:
Ahmed,
        Thanks. I agree with what you see, but there is also another reason
for the difference. For the space point in ch2 the engine and hcana
choose different left/right combinations. That is why you see the large
jump in the summed chi2 from planes 9 and 10. The two codes pick
different left/right because there is a cut on the difference between the
measured stub xp and an expected xp = sp->GetX()/875  and one code has a combination
slightly below the cut value and the other slightly above. I do not see a way
to make the codes agree more exactly since the stub fitting involves a matrix
inversion.

                Cheers,
                       Mark

Bugs fixed.

  Goldent track is moved to THcHallCSpectrometer::TrackCalc()
  if  fSelUsingScin == 0 then golden track is calculated just
  like podd. i.e. it is the first track with minimum chi2/ndf
  with sorting ON

  if fSelUsingScin == 1 then golden track is calculetd just like
  engine/HTRACKING/h_select_best_track_using_scin.h. This method
  gives the best track with minimum value of chi2/ndf but with
  additional cuts on the tracks. These cuts are on dedx, beta
  and on energy.

  This commit also corrects and modifies the beta code written
  in the THcHodoscope::FineProcess

  At this point we have godlend track and in additoion to
  previous vlues the following values are set to a track:
  Energy, beta, chi2 of beta and dedx[track][first hit].

NOTE: For some tracks we have large differences in value of chi
      square between ENGINE and hcana and this is already
      discussed in the last group meeting. Following emails are
      exchanged regarding this difference of chisquare.

email to Mark:

Hi Mark,
I tried to trace the origin of difference of the chi square value between
ENGINE and hcana. Please see the attached file, it has long table which I
cannot present in this email. Or you can see the file at the following
link, thanks.

https://userweb.jlab.org/~zahmed/chi2diff.txt

Ahmed

Mark's reply:
Ahmed,
        Thanks. I agree with what you see, but there is also another reason
for the difference. For the space point in ch2 the engine and hcana
choose different left/right combinations. That is why you see the large
jump in the summed chi2 from planes 9 and 10. The two codes pick
different left/right because there is a cut on the difference between the
measured stub xp and an expected xp = sp->GetX()/875  and one code has a combination
slightly below the cut value and the other slightly above. I do not see a way
to make the codes agree more exactly since the stub fitting involves a matrix
inversion.

                Cheers,
                       Mark

---------------------------------------------------------------------------------------

This commmit is only for Beta. THcHodoscope is modified and few
bugs are fixed. I will not make pull request for this commit.
Pull request will be made after the completion of Golden track
in THcHallCSpectrometer.

shower counter's FineProcess.

       Remove THcHallCSpectrometer::TrackEnergies method.
       Add GetShEnergy call for each tack in THcShower::fineProcess.

memory consumption with event number.

       Add THcShowerClusterList::purge() method to purge and clean
       THcShower::fClusterList by the end of each event, in
       THcShower::Clear().

crashes.

    Add TrackEnergies method in THcHallCSpectrometer class. It is
    called from TrackCalc metod.

calorimeter.

Add list of clusters fClusterList in THcShower class.
Add method THcShower::GetShEnergy(THaTrack*).
Add method THcShowerClusterList::clear() in THcShowerCluster.h.

Increment fNtracks++ each time a cluster matches a spectrometer
track. Only best matches were counted before. Move fNtracks++ from
THcShower::CoarceProcess to THcShower:MatchCluster.

In tree they are called tnhit and trawhit.
Created two new variables for THcDriftChamber, fNthits and fN_True_RawHits.
In tree they are called tnhit and trawhit.
THcDC defined subdetector THcDriftChamber to have name Ch1 and Ch2 to do used in the tree.
THcDriftChamberPlane defined new method GetNRawhits to return fNRawhits for each plane.
Added member fNRawhits to THcDriftChamberPlane

    Time difference variables now defined in output.def

Remove track related variables from the HMS Shower class, use ntuples
from the track (H.tr) block instead in the HMS calorimeter calibration
scripts.

- Remove a block of duplicated track related variables from the
THcShower class which have been added for calorimeter calibration
purposes.

- In hcal_calib directory, replace those variables with correspondent
variables from H.tr block in the calorimeter calibration class
(THcShowerCalib.h), in the cut (hcal_replay_cuts.def) and output
(output_hcal_replay.def) definitions.

      Compute time differences between planes
      Calculate beta and Chisquard for each track in the hodoscope.
      New variables: fBeta, fBetaChisq
      New histograms for time differences

Conversion of h_tof.f to THcHodoscope::CoarseProcess method.

 Quick summary of third version at  branch zafar1:

 New functions in THcScintillatorPlane.h
 ---------------------------------------
    TClonesArray* GetPosADC() { return fPosADCHits;}; // Ahmed
    TClonesArray* GetNegADC() { return fNegADCHits;}; // Ahmed
    TClonesArray* GetPosTDC() { return fPosTDCHits;}; // Ahmed
    TClonesArray* GetNegTDC() { return fNegTDCHits;}; // Ahmed

 New data members in THcHodoscope.h
 ----------------------------------
    TClonesArray* scinPosADC; // Ahmed
    TClonesArray* scinNegADC; // Ahmed
    TClonesArray* scinPosTDC; // Ahmed
    TClonesArray* scinNegTDC; // Ahmed

 Action in THcHodoscope::CoarseProcess
 -------------------------------------
 scinPosADC = fPlanes[ip]-GetPosADC();
 hitPaddle  = (THcSignalHit*)scinPosADC->At(ihit))->GetPaddleNumber()-1;

Added the raw raster ADC values to the histos and trees. Modified the ranges and binning of the
histogram definitions in output_bpw_.def accordingly.

In raster_jun04.map, reassigned the raster channels into a single plane instead of different planes.
This makes it easier for comparision with the engine.
Modified raster_test.C to analyze run 52947 in the cache directory.
Some minor clean up in the THcRaster.cxx class- removing comments and unused variables.

  Don't put [1] at end of variable name as is done with
  numerical parameters with Define.

  Use new CCDB methods to get table row and column counts

Read gpbeam in THcRaster::ReadDatabase
Add jan05.kinematics
Modified DBASE/raster_test.database to fill g_ctp_kinematics_filename
with DBASE/jan05.kinematics
Modified raster_test.C to load g_ctp_kinematics_filename

  Title now copied to THcParmList variable
  String variables copied.  (Only first element of string array copied)
  Parameter created as int or double depending on CCDB table type
  Multiple column tables ignored.

  CCDB must be installed and CCDB_HOME defined before compiling.
  If CCDB_HOME is not defined, hcana is compiled without CCDB support
  scons has not been updated to compile with CCDB support
  To use, instead of gHcParms->Load(filename), do
    gHcParms->OpenCCDB(runnum, ccdb_connection_string)
    gHcParms->LoadCCDBDirectory("hms","h")
    gHcParms->LoadCCDBDirectory("sos","s")
    gHcParms->LoadCCDBDirectory("gen","g")
  Parmaters will be loaded into the same internal gHcParms list as with
  ThcParmList::Load, so no changes to detector classes are needed.
  Scripts to convert existing parameters into CCDB are under development.
  Currently all ccdb constants are loaded as doubles as there is not
  a way to easily query ccdb to get the data type.
  Also strings are not loaded.

  Fixes not necessarily correct, just avoid the errors.

  Files with the prefix raster_ are coppies of original *.database,
  and *.param files which now includes changes nessacery for
  decoding raster data. I created seperate files for these because I
  need to merge this code upstream before I can compare the raster
  decoding results with that of the analyzer. For e.g. the map files are
  different from the original map files. So Steve needs to create a perl
  script that can propgate these changes in the map files to the map
  file of the analyzer inorder for me to do the analyzer hcana
  comparision. Once the comparision is done and we are satisifies with
  the results, we will merge the changes with the original files and
  remove the files with the prefix raster_*.

  gbeam.param now included by general.param.  No need to load
  a separate parameter set to test raster/BPM code.

  examples/MAPS/raster_jun04.map adds new detector IDs for the raster
  and BPMs as follows:

Add the following IDs to the misc elements,
BEAM
 RASTER_ID= 18          ADC
 BPM_ID= 19             ADC

In the channel section did the following modification:
 sub , plane, counter, sig
  14, 1, 13, 1 !ADC ROC#1 Sl#15 Ch#14 Fast Raster X-sync
  15, 1, 14, 2 !ADC ROC#1 Sl#15 Ch#15 Fast Raster X-signal
  12, 1, 15, 3 !ADC ROC#1 Sl#15 Ch#12 Fast Raster Y-sync
  13, 1, 16, 4 !ADC ROC#1 Sl#15 Ch#13 Fast Raster Y-signal

  THcRasterRawHit class to helps with the decoding of the raster
     signals via THcHitList class.
  THcRasteredBeam class derives from THaApparatus to handle rastered beam.
  Modify Makefile and HallC_LinkDef.h to compile new classes

Added coding for read parameters from the hcana.param file.

In the hcana.param file added the four raster channels
graster_channels = "xsync xsig ysync ysig"

Created raster_test.C script to run the analyzer for raster code testing purpose. This was coppied from the hodotest.C

add flags debugflagpr and debugflagstubs
add method THcDC::PrintStubs
modify THcDC::PrintSpacePoints
Set THcDC::LinkStubs() to only set fNDCTracks=0 and retrun when
  fNDCTracks>10 for HMS style chambers to match engine.

In CalcFptime changes for loop to use  fNScinGoodHits
In calcFptime if not a good hit ( both PMT within time window) set to -10000.
In Clear set fpTime=-10000.