////////////////////////////////////////////////////////////////////////
/// \class RAT::Classifiers::AlphaClassifierSeeded
///
/// \brief  Classifies post-reconstruction events as either Alpha or Electron
///
/// \author name Krish Majumdar <Krishanu.Majumdar@physics.ox.ac.uk>
///              based on MeanTime by Phil G Jones <p.g.jones@qmul.ac.uk>
///
/// REVISION HISTORY:\n
///     07/05/2012 : Krish Majumdar - New file \n
///
/// \details Looks at the distribution of Time Residuals, as a smaller number
///         of prompt PMT hits indicates an alpha.  Will return the ratio of
///         prompt PMT hits to the total number of PMT hits in the event.
///
////////////////////////////////////////////////////////////////////////

#ifndef __RAT_Classifiers_AlphaClassifierSeeded_
#define __RAT_Classifiers_AlphaClassifierSeeded_

#include <string>
#include <RAT/SeededClassifier.hh>
#include <RAT/DU/LightPathCalculator.hh>
#include <TVector3.h>

namespace RAT
{
  namespace DS
  {
    class FitResult;
  }

  namespace Classifiers
  {
    class AlphaClassifierSeeded : public SeededClassifier
    {

    public:
      virtual std::string GetName() const { return AlphaClassifierSeeded::Name(); }
      static std::string Name() { return std::string("SeededAlpha"); }
      void Initialise(const std::string&) { }
      void BeginOfRun(DS::Run& run);
      void EndOfRun(DS::Run&) { }
      virtual DS::ClassifierResult GetClassification();
      void DefaultSeed();
      void SetSeed(const DS::FitResult& seed);

    private:
      double ft1;  // Lower limit of the "prompt" PMT Hit range
      double ft2;  // Upper limit of the "prompt" PMT Hit range
      double ftmax;
      TVector3 fEventPos;
      double fEventTime;
      DU::LightPathCalculator fLightPath; ///< Light path calculator for a given run
    };
  } //::Classifier
} //::RAT

#endif