raytrace_inel.h

Go to the documentation of this file.

#pragma once
#include "raytrace.h"
#include "superarray.h"
 
namespace GOAT
{
    namespace raytracing
    {
        constexpr int INEL_CALCPHASE_EXCITATION=1;
        constexpr int INEL_CALCPHASE_RRT=2;
        constexpr int INEL_CALCPHASE_EXCITATION_ONLY=3;
 
        constexpr int INEL_MAX_NREFLEX=1;
        constexpr int INEL_SAVE_ABSE=1;
 
        constexpr int INEL_RADIATION_TYPE_FLOURESCENCE=1;
        constexpr int INEL_RADIATION_TYPE_RAMAN=2;
 
        constexpr int INEL_RADIATION_COHERENT=1;
        constexpr int INEL_RADIATION_INCOHERENT=2;
 
        constexpr int INEL_EXPORT_EXCITATION_FIELD_ABS=0;
        constexpr int INEL_EXPORT_EXCITATION_FIELD_VECTOR=1;
 

 
         struct RRTParms
        {
            maths::Vector<double> P; 
            maths::Vector<double> n; 
            maths::Vector<double> e1, e2; 
            double wvlinel=1.0; 
            int radiationType = INEL_RADIATION_TYPE_RAMAN; 
            int coherency = INEL_RADIATION_INCOHERENT; 
        } ;
 

        RRTParms calcDetDirParms (double theta, double phi, double wvlinel); 

        class Raytrace_Inel :
            public Raytrace
        {
        public:
            Raytrace_Inel();
            Raytrace_Inel(const Scene &S); 
            void trace(RRTParms D); 
            void traceEnterObject(); 
            void traceLeaveObject(); 
            void sceneChanged(const Scene &S); 
            void exportExcitation(std::string fname, int savetype = INEL_EXPORT_EXCITATION_FIELD_ABS); 
            double inel1, inel2; 
            void resetCalculation() { calcphase = INEL_CALCPHASE_EXCITATION; } 
            void setExcitationFieldOnly() { calcphase = INEL_CALCPHASE_EXCITATION_ONLY; } 
            void unsetExcitationFieldOnly() { calcphase = INEL_CALCPHASE_RRT;  } 
            SuperArray<maths::Vector<std::complex<double> > >* SGE=nullptr;  
             SuperArray<GOAT::maths::Vector<std::complex<double>>>*
                GetFieldArray() 
            {
                return this->SGE;
            }
 
             void setScene(const Scene& S);
        private:
            std::complex<double>  gewichte(maths::Vector<std::complex<double> > E, maths::Vector<std::complex<double> > p);
            void initExcitation();
            void initRRT();
 
            void traceExcitation();
            void saveExcitation();
            void saveRRT();
            void traceRRT();
            int calcphase;
            SuperArray<maths::Vector<std::complex<double> > >* SGRRT1 = nullptr;
            SuperArray<maths::Vector<std::complex<double> > >* SGRRT2 = nullptr;
            
            bool* active;
            int n;
            int iR;
            LightSrc** LSExcit;
            int nLSExcit;
            RRTParms parms;
        };
    }
}