////////////////////////////////////////////////////////////////////////
/// \class RAT::Classifiers::IsotropyRegions
///
/// \brief  Classifies the isotropy of an event
///
/// \author name Phil G Jones <p.g.jones@qmul.ac.uk>
/// \author Evelina Arushanova <e.arushanova@qmul.ac.uk> -- contact person
///
/// REVISION HISTORY:\n
///  - 26/06/2011 : P G Jones - New file
///  - 2015/04/02 : M Mottram - Fix for divide by zero
///
/// \details Returns the chi squared difference in the observed hits in
///         regions to the expected (total hits/num regions). For all
///         the regions and the average for equal latitude regions.
///
////////////////////////////////////////////////////////////////////////

#ifndef __RAT_Classifiers_IsotropyRegions__
#define __RAT_Classifiers_IsotropyRegions__

#include <RAT/SeededClassifier.hh>

#include <TVector3.h>

namespace RAT
{
namespace DS
{
  class FitResult;
}

namespace Classifiers
{

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

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

  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:
  TVector3 fEventPos;
  double fEventTime;
};

} //::Classifier

} //::RAT

#endif