JVectorNZ.hh

Go to the documentation of this file.

#ifndef __JFIT__JVECTORNZ__
#define __JFIT__JVECTORNZ__

#include <vector>
#include <ostream>

#include "Jeep/JPrint.hh"
#include "JFit/JLine1Z.hh"
#include "JMath/JVectorND.hh"
#include "JGeometry3D/JVector3D.hh"


/**
 * \author mdejong
 */

namespace JFIT {}
namespace JPP { using namespace JFIT; }

namespace JFIT {

  using JMATH::JVectorND;


  /**
   * Determination of the time residual vector of hits for a track along z-axis (JFIT::JLine1Z).
   */
  class JVectorNZ : 
    public JVectorND
  {
  public:
    /**
     * Default contructor.
     */
    JVectorNZ()
    {}
              

    /**
     * Constructor.
     *
     * The template argument <tt>T</tt> refers to an iterator of a data structure which should have the following member methods:
     *   - double getX();       // [m]
     *   - double getY();       // [m]
     *   - double getZ();       // [m]
     *   - double getT();       // [ns]
     *
     * \param  track            track
     * \param  __begin          begin of data
     * \param  __end            end   of data
     */
    template<class T>
    JVectorNZ(const JLine1Z& track,
              T              __begin,
              T              __end)
    {
      set(track, __begin, __end);
    }


    /**
     * Set time residual vector.
     *
     * The template argument <tt>T</tt> refers to an iterator of a data structure which should have the following member methods:
     *   - double getX();       // [m]
     *   - double getY();       // [m]
     *   - double getZ();       // [m]
     *   - double getT();       // [ns]
     *
     * \param  track            track
     * \param  __begin          begin of data
     * \param  __end            end   of data
     */
    template<class T>
    void set(const JLine1Z& track,
             T              __begin,
             T              __end)
    {
      using JGEOMETRY3D::JVector3D;

      this->clear();

      for (T hit = __begin; hit != __end; ++hit) {
        this->push_back(hit->getT() - track.getT(JVector3D(hit->getX(), 
                                                           hit->getY(), 
                                                           hit->getZ())));
      }
    }


    /**
     * Print ASCII formatted output.
     *
     * \param  out      output stream
     * \param  Y        Y vector
     * \return          output stream
     */
    friend inline std::ostream& operator<<(std::ostream& out, const JVectorNZ& Y)
    {
      using namespace std;

      JFlags flags(out);

      for (JVectorNZ::const_iterator i = Y.begin(); i != Y.end(); ++i) {
        out << showpos << SCIENTIFIC(10,2) << *i << ' ';
      }

      out << endl;

      return out;
    }
  };
}

#endif