////////////////////////////////////////////////////////////////////////
/// \class RAT::Classifiers::Isotropy
///
/// \brief Calculates the Isotropy parameters for an event
///
/// \author name Gene Beier <geneb@hep.upenn.edu>
///
/// REVISION HISTORY:\n
/// 09/26/2014 : Gene Beier - New file /n
///
/// \detail Calculates theta_ij, beta_1, beta_2, beta_3, beta_4 and
/// snobeta_14 for an event using the prompt time cut
/// specified in CLASSIFIER_ISOTROPY.ratdb. The algorithm and
/// time cut for snobeta_14 reproduce the SNO beta_14
/// exactly.
///
////////////////////////////////////////////////////////////////////////

#ifndef __RAT_Classifiers_Isotropy__
#define __RAT_Classifiers_Isotropy__

#include <RAT/SeededClassifier.hh>
#include <RAT/SelectorClassifier.hh>
#include <RAT/DU/LightPathCalculator.hh>
#include <RAT/DS/FitResult.hh>

namespace RAT
{
  namespace DS
  {
    class FitResult;
  }

  namespace Classifiers
  {

    class Isotropy : public SeededClassifier, public SelectorClassifier
    {
    public:
      virtual std::string GetName() const { return Isotropy::Name(); }

      static std::string Name() { return std::string( "isotropy" ); }

      void Initialise(const std::string&){ };

      void BeginOfRun(DS::Run& ) { };

      void EndOfRun(DS::Run&) { };

      virtual DS::ClassifierResult GetClassification();
      /// Set the seed to the default
      void DefaultSeed();
      /// Set the seed from seed
      void SetSeed( const DS::FitResult& seed );

    private:
      DS::FitVertex fSeedVertex;
      DU::LightPathCalculator fLightPath;
      TVector3 fEventPos;
    };

  } // ::Classifiers
} // :: RAT

#endif // __RAT_Classifiers_Isotropy__