NCerNPE and fNormETot be optional parameters
When there is not the  !fChern || !fShower and instead just
not increment the counter.

If calorimeter and cherenkov is not added to SOS or HMS then it prints
a warning message and does not calcualte the efficiency for that
spectrometer.

Minimm value of normalized energy and minum value of number of
cherenkov photo electrons is added to parameter files for the
calculation of tracking efficiency.

  Variable found in a report template is first tested to see if it is a
  Hall C style parameter (in gHcParms).  If not, then tested to see if it
  is a global variable (in gHaVars).

  Move some stuff outside of # tracks > 0 test

  Clean up THcShower.h and THcShower.cxx.
  typedef vector<THcShowerHit*> THcShowerHitList by
  typedef set<THcShowerHit*> THcShowerHitList.
  Remove inheritance of class THcShowerCluster from THcShowerHitList.
  Define THcShowerCluster as THcShowerHitList, aka set<THcShowerHit*>.
  Remove inheritance of class THcShowerClusterList from
        THcShClusterList, a.k.a. vector<THcShowerCluster*> container.
  Define THcShowerClusterList as vector<THcShowerCluster*>.
  Rename THcShowerHitSet by THcShowerHitList.
  Rename HitList (related to the HMS calorimeter) to HitSet.
  Add operator< in THcShowerHit class, in order to have THcShowerHitSet
    objects to be properly sorted.

Modified THcShower and THcHallCSpectrometer

    THcHallCSpectrometer inherits from THaSpectrometer which has method
    GetPcentral and a member fPcentral.
     THcHallCSpectrometer had been using fPCentral but decided to
    use fPcentral.
    Delete fPCentral in THcHallCSpectrometer.h
    Replace all fPCentral  by fPcentral in ReadDatabase and FindVertices methods.

    Add variable fEtot= the total energy of the calorimeter which
      is calculated in Decode method
    Add variable fEtotnorm= the total energy of the calorimeter
      divided by central spectrometer which is calculated in Decode method
    Add both to DefineVariables list as etot and etotnorm

      Declare loop variables in the for statement
      Go back to UInt_t for some variables
      Don't prefix non-member variables with f

  Adds tests to THcHodoscope fiducial tracking efficiency
  Adds tests to THcCherenkov for Cerenkov efficiency

       State neighbours hits in a row separated by single block in
       THcShower::IsNeighbour method. This allows for inclusion of a
       stand alone hit of low energy deposition (presumably coming
       from noise in a channel) into a big cluster, hence reducing
       differences between engine and hcana.

in hcana, in part of the HMS calorimeter.

   Consistently use variables of the track proper
   in the HMS calibartion code.

  When it was UInt_t, variable length global variables didn't get the
  length correctly and there were run time errors.

  For the most part this is making sure that the loop variable in
  a for loop is of the same type as the variable that it is compared
  against (Often a UInt_t).

This reverts commit e9ce1040.

Conflicts:
	src/THcHallCSpectrometer.cxx

  For the most part this is making sure that the loop variable in
  a for loop is of the same type as the variable that it is compared
  against (Often a UInt_t).

This is a temporary fix for problem that arise when compiling on the
ifarm machines. In the DefinedValue method when using
kCutScaler and kCutNCalled the switch statement goes to the default
case.

On ifarm gcc version 4.4.6 20110731 (Red Hat 4.4.6-3) (GCC)
There are compiler warnings
 src/THcFormula.cxx: In member function 'virtual Double_t THcFormula::DefinedValue(Int_t)':
src/THcFormula.cxx:131: warning: case label value exceeds maximum value for type
src/THcFormula.cxx:134: warning: case label value exceeds maximum value for type
One permenant fix would be to add kCutScaler and kCutNCalled to the
enum EVariableType in THaFormula.h.

Add THaSpectrometer::DefineVariables( mode );
to define those tree variables

First, in h_pattern_recognition.f we have:

https://github.com/zahmed9/engine/blob/mkj/HTRACKING/h_pattern_recognition.f#L139
(hdc_wire_center(yy)-hdc_wire_center(yyprime))**2 .lt. (hspace_point_criterion(ich))

So I changed
https://github.com/zahmed9/hcana/blob/develop/src/THcDriftChamber.cxx#L292
TMath::Abs(fHits[plane_hitind]->GetPos() - fHits[planep_hitind]->GetPos()) < fSpacePointCriterion

to

pow((fHits[plane_hitind]->GetPos() - fHits[planep_hitind]->GetPos()),2) < fSpacePointCriterion

Second, the value of MAX_SPACE_POINTS in ENGINE is 100 so in hcana I changed it from 50 to 100.

  This is because the classes in these header files are only
  used by THcShower.

     Edit debug outputs in THcSHower and THcShowerPlane classes. Group
     outputs preferably by the end of methods, for smooth code
     flow. Structure tests of outputs for easy reading.

the HMS calorimeter.

    Add src/THcShowerCluster.h.float, a version of THcShowerCluster
    class with float precision calculation in clX method.

first track from the HMS shower counter code.

    Remove variables and code related to the maximum cluster of hits,
    and association of a cluster to the 1-st spectrometer track in
    THcShower.h and THcShower.cxx. Those were used for testing of the
    code, not needed any more.

track in Engine and the golden track in hcana.

      Add directory examples/gtrk_e holding the code.

consistent with Engine.

	   Reverse order of cycling through clusters in THcShower::MatchCluster.

Access the energy of the golden track from THaGoldenTrack.

    Include THaTrack.h in THaGoldenTrack.h.
    Add GetEnergy method in THaGoldenTrack.h.
    Define variable name "e" for the energy of golden track in THaGoldenTrack.C.