EvtVector3C.hh

Go to the documentation of this file.

 
/***********************************************************************
* Copyright 1998-2020 CERN for the benefit of the EvtGen authors       *
*                                                                      *
* This file is part of EvtGen.                                         *
*                                                                      *
* EvtGen is free software: you can redistribute it and/or modify       *
* it under the terms of the GNU General Public License as published by *
* the Free Software Foundation, either version 3 of the License, or    *
* (at your option) any later version.                                  *
*                                                                      *
* EvtGen is distributed in the hope that it will be useful,            *
* but WITHOUT ANY WARRANTY; without even the implied warranty of       *
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the        *
* GNU General Public License for more details.                         *
*                                                                      *
* You should have received a copy of the GNU General Public License    *
* along with EvtGen.  If not, see <https://www.gnu.org/licenses/>.     *
***********************************************************************/
 
#ifndef EVTVECTOR3C_N
#define EVTVECTOR3C_N
 
#include "EvtGenBase/EvtComplex.hh"
#include "EvtGenBase/EvtVector3R.hh"
 
#include <iosfwd>
 
class EvtVector3C final {
    friend EvtVector3C rotateEuler( const EvtVector3C& v, double phi,
                                    double theta, double ksi );
 
    inline friend EvtVector3C operator*( const EvtComplex& c,
                                         const EvtVector3C& v2 );
    inline friend EvtVector3C operator*( const EvtComplex& c,
                                         const EvtVector3R& v2 );
    inline friend EvtComplex operator*( const EvtVector3R& v1,
                                        const EvtVector3C& v2 );
    inline friend EvtComplex operator*( const EvtVector3C& v1,
                                        const EvtVector3R& v2 );
    inline friend EvtComplex operator*( const EvtVector3C& v1,
                                        const EvtVector3C& v2 );
    inline friend EvtVector3C operator+( const EvtVector3C& v1,
                                         const EvtVector3C& v2 );
    inline friend EvtVector3C operator-( const EvtVector3C& v1,
                                         const EvtVector3C& v2 );
    inline friend EvtVector3C operator*( const EvtVector3C& v1,
                                         const EvtComplex& c );
 
  public:
    EvtVector3C();
    EvtVector3C( const EvtComplex&, const EvtComplex&, const EvtComplex& );
    inline void set( const int, const EvtComplex& );
    inline void set( const EvtComplex&, const EvtComplex&, const EvtComplex& );
    inline void set( double, double, double );
    inline EvtVector3C& operator*=( const EvtComplex& c );
    inline EvtVector3C& operator/=( const EvtComplex& c );
    inline EvtVector3C& operator+=( const EvtVector3C& v2 );
    inline EvtVector3C& operator-=( const EvtVector3C& v2 );
    inline EvtVector3C( const EvtVector3R& v1 );
    void applyRotateEuler( double phi, double theta, double ksi );
    inline const EvtComplex& get( int ) const;
    inline EvtVector3C conj() const;
    EvtVector3C cross( const EvtVector3C& v2 );
    friend std::ostream& operator<<( std::ostream& c, const EvtVector3C& v );
    double dot( const EvtVector3C& p2 );
 
  private:
    EvtComplex m_v[3];
};
 
inline EvtVector3C::EvtVector3C( const EvtVector3R& v1 )
{
    m_v[0] = EvtComplex( v1.get( 0 ), 0.0 );
    m_v[1] = EvtComplex( v1.get( 1 ), 0.0 );
    m_v[2] = EvtComplex( v1.get( 2 ), 0.0 );
}
 
inline void EvtVector3C::set( const int i, const EvtComplex& c )
{
    m_v[i] = c;
}
 
inline void EvtVector3C::set( const EvtComplex& x, const EvtComplex& y,
                              const EvtComplex& z )
{
    m_v[0] = x;
    m_v[1] = y;
    m_v[2] = z;
}
 
inline void EvtVector3C::set( double x, double y, double z )
{
    m_v[0] = EvtComplex( x );
    m_v[1] = EvtComplex( y );
    m_v[2] = EvtComplex( z );
}
 
inline const EvtComplex& EvtVector3C::get( int i ) const
{
    return m_v[i];
}
 
inline EvtVector3C& EvtVector3C::operator*=( const EvtComplex& c )
{
    m_v[0] *= c;
    m_v[1] *= c;
    m_v[2] *= c;
    return *this;
}
 
inline EvtVector3C& EvtVector3C::operator/=( const EvtComplex& c )
{
    m_v[0] /= c;
    m_v[1] /= c;
    m_v[2] /= c;
    return *this;
}
 
inline EvtVector3C& EvtVector3C::operator+=( const EvtVector3C& v2 )
{
    m_v[0] += v2.m_v[0];
    m_v[1] += v2.m_v[1];
    m_v[2] += v2.m_v[2];
    return *this;
}
 
inline EvtVector3C& EvtVector3C::operator-=( const EvtVector3C& v2 )
{
    m_v[0] -= v2.m_v[0];
    m_v[1] -= v2.m_v[1];
    m_v[2] -= v2.m_v[2];
    return *this;
}
 
inline EvtVector3C operator+( const EvtVector3C& v1, const EvtVector3C& v2 )
{
    return EvtVector3C( v1 ) += v2;
}
 
inline EvtVector3C operator-( const EvtVector3C& v1, const EvtVector3C& v2 )
{
    return EvtVector3C( v1 ) -= v2;
}
 
inline EvtVector3C operator*( const EvtVector3C& v1, const EvtComplex& c )
{
    return EvtVector3C( v1 ) *= c;
}
 
inline EvtVector3C operator*( const EvtComplex& c, const EvtVector3C& v2 )
{
    return EvtVector3C( v2 ) *= c;
}
 
inline EvtVector3C operator*( const EvtComplex& c, const EvtVector3R& v2 )
{
    return EvtVector3C( v2 ) *= c;
}
 
inline EvtComplex operator*( const EvtVector3R& v1, const EvtVector3C& v2 )
{
    return v1.get( 0 ) * v2.m_v[0] + v1.get( 1 ) * v2.m_v[1] +
           v1.get( 2 ) * v2.m_v[2];
}
 
inline EvtComplex operator*( const EvtVector3C& v1, const EvtVector3R& v2 )
{
    return v1.m_v[0] * v2.get( 0 ) + v1.m_v[1] * v2.get( 1 ) +
           v1.m_v[2] * v2.get( 2 );
}
 
inline EvtComplex operator*( const EvtVector3C& v1, const EvtVector3C& v2 )
{
    return v1.m_v[0] * v2.m_v[0] + v1.m_v[1] * v2.m_v[1] + v1.m_v[2] * v2.m_v[2];
}
 
inline EvtVector3C EvtVector3C::conj() const
{
    return EvtVector3C( ::conj( m_v[0] ), ::conj( m_v[1] ), ::conj( m_v[2] ) );
}
 
#endif