////////////////////////////////////////////////////////////////////////
/// \class RAT::PDFs::PositionDirectionPDF
///
/// \brief  PDF that combines position (GV1D) and direction (FIT_DIR) pdfs.
///
/// \author Matt Mottram <m.mottram@qmul.ac.uk> -- contact person
///
/// REVISION HISTORY:\n
///     2014/03/20: Matt Mottram - new file \n
///
////////////////////////////////////////////////////////////////////////

#ifndef __RAT_PDF_PositionDirectionPDF_
#define __RAT_PDF_PositionDirectionPDF_

#include <RAT/PDF.hh>
#include <RAT/DU/LightPathCalculator.hh>

#include <string>

namespace RAT
{

class FitterPMT;

namespace DS
{
  class FitVertex;
}

namespace PDFs
{

class PositionDirectionPDF : public PDF
{
public:
  virtual double
  GetProbability( const FitterPMT& pmt,
                  const DS::FitVertex& vertex );

  virtual std::string GetName() const { return Name(); }
  static std::string Name() { return std::string( "positionDirectionPDF" ); }

  void BeginOfRun( DS::Run& run );

  void EndOfRun( DS::Run& run );

  /// Initialise the fAngularProbability and fTimingProbability PDFs
  virtual void
  Initialise( const std::string& param );
protected:
  G4PhysicsOrderedFreeVector* fAngularProbability; /// <PDF of hit angles relative to primary direction
  G4PhysicsOrderedFreeVector* fTimingProbability; /// <PDF of hit time residuals
  double fConcentratorRadius; /// <PMT concentrator radius
  std::string fIndex; ///< Optional database index
  DU::LightPathCalculator fLightPath; ///< Light path calculator for a given run
};

} //::PDFs

} //::RAT

#endif