pulsecalculation_field.h

Go to the documentation of this file.

#pragma once
 
#include "superarray.h"
#include "fft.h"
#include "raytrace_field_usp.h"
#include <vector>
#include "superarray.h"
#include <sstream>
#include <locale>
#include <iomanip>
 
 
namespace GOAT
{
    namespace raytracing
    {
        class pulseCalculation_Field
        {
            public:
                pulseCalculation_Field(Scene S);
                double findHitTime(int ObjNo); 
                void fieldCalculation(); 
                void fieldCalculation(double omega); 
                void setPulseWidth(double dt); 
                void setSpatialResolution(double dx); 
                void setRefractiveIndexFunctions(std::vector<std::function<std::complex<double>(double) > > nList); 
                void setNumReflex(int numReflex);
                void setSpectralRanges(int n); 
                void setNumWavelengthsPerRange(int nS);
                void setCenterWavelength(double wvl); 
                void setBandwidth(double dWvl); 
                void setRepetitionRate(double rep);
                void setPeriod (double time);
 
                void field(double t); 
                void reset(); 
                void setReferenceTime(double tref);
                Trafo trafo;
                SuperArray<GOAT::maths::Vector<std::complex<double> > > SAres;
                    std::vector<SuperArray<std::vector<gridEntry> > >  SA; 
                    Raytrace_Field_usp rt;
                    void addBoxDetector(Box* box) { BoxDetector.push_back(box); }
                    void saveIntensity(std::string FName, int i);
                    TrafoParms getTrafoParms() { return trafoparms; }
                    double getReferenceTime() { return tref; }
                    INDEX_TYPE getNumCellsPerDirection() { return nn; }
 
 
            private:    
                /* In this function the default values (trafoparms) for the calculations are set as follows:
                * dt  : 1E-14s
                * wvl : 1.0�m
                * nI  : 1
                * 
                */
                void setDefaults();                              
                void calcTrafoParms();
                // std::vector< std::vector<SuperArray<std::vector<gridEntry> > > > SA;
                double domega; 
                double dWvl=0.02;  
                double dRWvl;      
                INDEX_TYPE  nn;    // number of cells over the whole width of the calculation space (i.e. 2*r0)
                Scene S;
                
                bool raytracingDone = false; 
                std::vector<Box*> BoxDetector;
                TrafoParms trafoparms;
                double tref = 0.0; 
                int numReflex = INEL_MAX_NREFLEX;                              
        };
    }
}