OpalParticle.h

Go to the documentation of this file.

//
// Class OpalParticle
//   This class represents the canonical coordinates of a particle.
//
// Copyright (c) 2008 - 2020, Paul Scherrer Institut, Villigen PSI, Switzerland
// All rights reserved
//
// This file is part of OPAL.
//
// OPAL 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.
//
// You should have received a copy of the GNU General Public License
// along with OPAL. If not, see <https://www.gnu.org/licenses/>.
//
#ifndef CLASSIC_OpalParticle_HH
#define CLASSIC_OpalParticle_HH
 
#include "Vektor.h"
 
class OpalParticle {
 
public:
 
    // Particle coordinate numbers.
    enum { X, Y, L, INVALID };
 
    //  Construct particle with the given coordinates.
    OpalParticle(int64_t id,
                 double x, double px,
                 double y, double py,
                 double z, double pz,
                 double time,
                 double q, double m);
 
    OpalParticle(int64_t id,
                 Vector_t const& R, Vector_t const& P,
                 double time, double q, double m);
 
    OpalParticle();
 
    void setX(double) ;
 
    void setPx(double);
 
    void setY(double) ;
 
    void setPy(double);
 
    void setZ(double) ;
 
    void setPz(double);
 
    void setR(Vector_t const&);
 
    void setP(Vector_t const&);
 
    void setTime(double t);
 
    int64_t getId() const;
 
    //  Access coordinate by index for constant particle.
    double operator[](unsigned int) const;
 
    double getX() const;
 
    double getPx() const;
 
    double getY() const;
 
    double getPy() const;
 
    double getZ() const;
 
    double getPz() const;
 
    const Vector_t& getR() const;
 
    const Vector_t& getP() const;
 
    double getTime() const;
 
    double getCharge() const;
 
    double getMass() const;
 
private:
    int64_t id_m;
    Vector_t R_m;
    Vector_t P_m;
    double time_m;
    double charge_m;
    double mass_m;
};
 
inline
void OpalParticle::setX(double val)
{
    R_m[X] = val;
}
 
inline
void OpalParticle::setY(double val)
{
    R_m[Y] = val;
}
 
inline
void OpalParticle::setZ(double val)
{
    R_m[L] = val;
}
 
inline
void OpalParticle::setPx(double val)
{
    P_m[X] = val;
}
 
inline
void OpalParticle::setPy(double val)
{
    P_m[Y] = val;
}
 
inline
void OpalParticle::setPz(double val)
{
    P_m[L]  = val;
}
 
inline
void OpalParticle::setR(Vector_t const& R)
{
    R_m = R;
}
 
inline
void OpalParticle::setP(Vector_t const& P)
{
    P_m = P;
}
 
inline
void OpalParticle::setTime(double t)
{
    time_m = t;
}
 
inline
int64_t OpalParticle::getId() const
{
    return id_m;
}
 
inline
double OpalParticle::operator[](unsigned int i) const
{
    PAssert_LT(i, 6u);
    return i % 2 == 0? R_m[i / 2]: P_m[i / 2];
}
 
inline
double OpalParticle::getX() const
{
    return R_m[X];
}
 
inline
double OpalParticle::getY() const
{
    return R_m[Y];
}
 
inline
double OpalParticle::getZ() const
{
    return R_m[L];
}
 
inline
double OpalParticle::getPx() const
{
    return P_m[X];
}
 
inline
double OpalParticle::getPy() const
{
    return P_m[Y];
}
 
inline
double OpalParticle::getPz() const
{
    return P_m[L];
}
 
inline
const Vector_t& OpalParticle::getR() const
{
    return R_m;
}
 
inline
const Vector_t& OpalParticle::getP() const
{
    return P_m;
}
 
inline
double OpalParticle::getTime() const
{
    return time_m;
}
 
inline
double OpalParticle::getCharge() const
{
    return charge_m;
}
 
inline
double OpalParticle::getMass() const
{
    return mass_m;
}
 
#endif // CLASSIC_OpalParticle_HH