pem.h

#ifndef PEM_H
#define PEM_H

#include "collideableobject.h"
#include <QTimer>
#include <random>
#include <utility>

typedef Eigen::Matrix<std::complex<double>, 2, 2> Matrix22cd;

/*!
 * \brief The PEM class
 * This is the class which defines the methods that a PEM instance needs to
 * implement.
 */
class PEM : public CollideableObject {
  Q_OBJECT
private:
  Matrix22cd m_polarizationMatrix;
  std::complex<double> m_phaseAmplitude;
  std::complex<double> m_angularFrequency;
  double time, m_radius;
  std::complex<double> m_i = {0, 1};
  int noiseState = 0;
  std::normal_distribution<> dist;
  std::mt19937 noise_gen;
  Eigen::Vector2cd Er;

  void calculatePolarisationMatrix();

public:
  PEM(Eigen::Vector3d location, int side, Eigen::Vector3d normal, double radius,
      std::complex<double> phaseAmplitude,
      std::complex<double> angularFrequency);

  ~PEM() override = default;
  ObjectType getType() override;

  void collide(Ray &ray, Eigen::Vector3d &pointOfInterception) override;
  bool intersect(Ray &ray, Eigen::Vector3d &pointOfInterception) override;
  void incrementTime();

  void setNoiseState(int state);
  void newNoise(std::normal_distribution<> d, std::mt19937 gen);
  Matrix22cd getPolarizationMatrix();

signals:
  void outputPolarisationUpdated(Matrix4cd polarisation);

};

#endif // PEM_H