particle.h

/*
 * particle.h
 *
 *      Author: kappl
 */

#ifndef PARTICLE_H_
#define PARTICLE_H_

#include "Iparticle.h"

class Particle : public IParticle
{
        private:
        void step(Random&);
        double getDeltaT() const;
        long getDeltaTime() const;

        protected:
        double r;
        double theta;
        double phi;
        double time;
        Tproperties spec;
        TstepData data;
        double dt;
        double heliosphereBoundary;
        double perPart;
        double kappa0;
        double BfieldScaling;
        int polarity;
        double angle;

        virtual double getDeltaR(Random&) const = 0;
        double getDeltaTheta(Random&) const;
        double getDeltaPhi() const;

        void calculate();

        double getV() const;
        double getP2() const;
        double getRigidity() const;
        double getBeta() const;
        double getPolarity() const;
        double getAngle() const;
        double getTanPsi() const;
        double convertDeg(double) const;
        double getLarmorRadius() const;
        double getB0OverB() const;
        double getB() const;
        double getTheta() const;
        double getPhi() const;
        double getTime() const;
        double getMass() const;
        double getCharge() const;
        int getHeaviside() const;
        double getKappaPerpendicular() const;
        double getKappaThetaTheta() const;

        public:
        Particle(const TmodelInformation&, const Tproperties&);
        virtual ~Particle();

        long simulate(Random&);
        double getR() const;
        double getT() const;
        double getHeliosphereBoundary() const;
        double getISFlux(int) const;
        double getISFlux(bool) const;
        double getISFlux(bool, int) const;
};


#endif /* PARTICLE_H_ */