detector.h

Go to the documentation of this file.

#ifndef DETECTOR_H
#define DETECTOR_H
 
#include "vector.h"
#include "box.h"
 
#include <iostream>
#include <vector>
 
namespace GOAT
{
    namespace raytracing
    {
        #define DETECTOR_PLANE 20000
        #define DETECTOR_ANGLE 20001
 

    class Detector
    {
    public:
        Detector(void);
        Detector(int n1, int n2); 
        Detector(const Detector& Det); 
        Detector& operator= (const Detector& D);
        maths::Vector<std::complex<double> >& operator () (int i1, int i2) { return D[i1][i2]; } 
        ~Detector(void);
 
        void clean(); 
        void clear(); 
        virtual bool cross(maths::Vector<double> P, maths::Vector<double> k, int& i1, int& i2, double& l) = 0;
        int N1(); 
        int N2(); 
        void setN1(int n1); 
        void setN2(int n2); 
        void setN(int n1, int n2); 
        void setPosition(maths::Vector<double> Pos) {this->P=Pos;} 
        double D1(); 
        double D2(); 
 
        void setD(double d1, double d2);
 
        void setD1(double d1);
 
        void setD2(double d2);
 
        int Type() { return type; } 
        bool load(const char* fn); 
        void save(const char* fn); 
        void saveabs(const char* fn);  
        void savePhase(const char* fn, int coord);  
        void savereal(const char* fn, int coord); 
        void saveimag(const char* fn, int coord); 
        std::vector <std::vector<maths::Vector<std::complex<double>>>> D; 
        //  maths::Vector<std::complex<double> >** D=nullptr; ///< Here, the data will be stored
            maths::Vector<double> position() { return P; } 
            maths::Vector<double> norm() { return n; } 
        friend std::ostream& operator << (std::ostream& os, Detector& D);
        std::string fname;
        void mult(double fac); 
        maths::Vector<double> gete1() { return e1; } 
        maths::Vector<double> gete2() { return e2; } 
    protected:
        maths::Vector<double> e1; 
        maths::Vector<double> e2; 
        maths::Vector<double> P; 
        maths::Vector<double> n;  
        void init(int n1, int n2); 
        double d1=0, d2=0;
        int n1=0, n2=0;
        int type=-1;
        friend class DetectorPlane;
    };
 

    class DetectorPlane : public Detector
    {
    public:
        DetectorPlane(void);
        DetectorPlane(maths::Vector<double> P, maths::Vector<double> n, double d, int N);
        DetectorPlane(maths::Vector<double> P, maths::Vector<double> e1, maths::Vector<double> e2, int n1, int n2);
        void setNorm(maths::Vector<double> n); 
        bool cross(maths::Vector<double> P, maths::Vector<double> k, int& i1, int& i2, double& l);  
    };
 
    /*class DetectorBox : public Detector
    {
    public:
        DetectorBox(maths::Vector<double> P, maths::Vector<double> d, maths::Vector<int> n);
        bool cross(maths::Vector<double> P, maths::Vector<double> k, int& i1, int& i2, double& l);
         
        std::vector<std::vector<std::vector<std::complex<double> > > > data;
        Box box;
    };*/
#define SAVE_X 0
#define SAVE_Y 1
#define SAVE_Z 2
#define SAVE_PHASE_X 3
#define SAVE_PHASE_Y 4
#define SAVE_PHASE_Z 5
#define SAVE_ABS 6
}
}
#endif