JShower3D.hh

Go to the documentation of this file.

#ifndef __JSHOWER3D__
#define __JSHOWER3D__
 
#include <istream>
#include <ostream>
 
#include "JIO/JSerialisable.hh"
 
#include "JGeometry3D/JAxis3D.hh"
#include "JGeometry3D/JTime.hh"
#include "JGeometry3D/JVertex3D.hh"
 
#include "JPhysics/JConstants.hh"
 
#include "JMath/JMathToolkit.hh"
 
 
/**
 * \author mdejong
 */
 
namespace JGEOMETRY3D {}
namespace JPP { using namespace JGEOMETRY3D; }
 
namespace JGEOMETRY3D {
 
  using JIO::JReader;
  using JIO::JWriter;
 
  /**
   * 3D shower.
   */
  class JShower3D :
    public JAxis3D,
    public JTime
  {
  public:
 
    using JTime::getT;
    using JTime::add;
    using JTime::sub;
    using JPosition3D::add;
    using JPosition3D::sub;
    using JDirection3D::getDirection;
 
 
    /**
     * Default constructor.
     */
    JShower3D() :
      JAxis3D(),
      JTime()
    {}
 
 
    /**
     * Constructor.
     *
     * \param  vertex        vertex
     * \param  dir           direction
     */
    JShower3D(const JVertex3D& vertex,
              const JVersor3D& dir) :
      JAxis3D(vertex.getPosition(), dir),
      JTime(vertex.getT())
    {}
 
 
    /**
     * Get vertex of this track.
     *
     * \return               vertex
     */
    JVertex3D getVertex() const
    {
      return JVertex3D(getPosition(), getT());
    }
 
 
    /**
     * Type conversion operator.
     *
     * \return               vertex
     */
    operator JVertex3D() const
    {
      return getVertex();
    }
 
 
    /**
     * Get arrival time of Cherenkov light at given position.
     *
     * \param  pos              position [m]
     * \return                  time     [ns]
     */
    inline double getT(const JVector3D& pos) const
    {
      using namespace JPHYSICS;
 
      return this->getT() + JMATH::getDistance(this->getPosition(), pos) * getIndexOfRefraction() * getInverseSpeedOfLight();
    }
 
 
    /**
     * Get photon direction of Cherenkov light on PMT.
     *
     * \param  pos              PMT position
     * \return                  direction
     */
    inline JVersor3D getDirection(const JVector3D& pos) const
    {
      JPosition3D D(pos);
 
      D.sub(this->getPosition());
 
      return JDirection3D(D);
    }
 
 
    /**
     * Get cosine angle of impact of Cherenkov light on PMT.
     *
     * \param  axis             PMT axis
     * \return                  cosine angle of impact
     */
    inline double getDot(const JAxis3D& axis) const
    {
      return getDirection(axis.getPosition()).getDot(axis.getDirection());
    }
 
 
    /**
     * Get cosine angle of emission of Cherenkov light on PMT.
     *
     * \param  pos              PMT position
     * \return                  cosine angle of impact
     */
    inline double getDot(const JVector3D& pos) const
    {
      return JDirection3D(pos - getPosition()).getDot(getDirection());
    }
 
 
    /**
     * Read shower from input.
     *
     * \param  in            input stream
     * \param  shower        shower
     * \return               input stream
     */
    friend inline std::istream& operator>>(std::istream& in, JShower3D& shower)
    {
      in >> static_cast<JAxis3D&>(shower);
      in >> static_cast<JTime&>  (shower);
 
      return in;
    }
 
 
    /**
     * Write shower to output.
     *
     * \param  out           output stream
     * \param  shower        shower
     * \return               output stream
     */
    friend inline std::ostream& operator<<(std::ostream& out, const JShower3D& shower)
    {
      out << static_cast<const JAxis3D&>(shower);
      out << ' ';
      out << static_cast<const JTime&>  (shower);
 
      return out;
    }
 
 
    /**
     * Read shower from input.
     *
     * \param  in            reader
     * \param  shower        shower
     * \return               reader
     */
    friend inline JReader& operator>>(JReader& in, JShower3D& shower)
    {
      in >> static_cast<JAxis3D&>(shower);
      in >> static_cast<JTime&>  (shower);
 
      return in;
    }
 
 
    /**
     * Write shower to output.
     *
     * \param  out           writer
     * \param  shower        shower
     * \return               writer
     */
    friend inline JWriter& operator<<(JWriter& out, const JShower3D& shower)
    {
      out << static_cast<const JAxis3D&>(shower);
      out << static_cast<const JTime&>  (shower);
 
      return out;
    }
  };
}
 
#endif