OpenPaperOpt/Component.h

Go to the documentation of this file.

#ifndef COMPONENT_H
#define COMPONENT_H
#include "StructureObject.h"
#include "Basesheet.h"
#include "WavePacket.h"
#include "Distribution.h"
#include "Fluorescence.h"
#include <vector>
#include "Transform.h"
#include "utilities.h"
#include <complex>
using namespace std;
class Component
{
public:
Fluorescence* mpFluorescence;
void SetFluorescent(Fluorescence* pFluor) { mpFluorescence = pFluor; }

protected:
        Component();

        Component(const StructureObject* const pStructureObject, vector<complex_LC>* pRefrIndex, 
                vector<double>* pQScat, vector<double>* pQAbs, vector<double>* pQAbsFluo,
                Distributions::Distribution* pDistScat, double ContactReduction, V2<float>* pRoughness);

public:
        virtual ~Component(void);

        virtual int Process(WavePacket& wp) const=0;
        virtual int Process(WavePacket& wp, int ComponentNumber, bool EntryScattering = false, complex_LC* pPrevRefrIndex = 0) const = 0;

protected:
        enum Occurrence
        {
                NOTHING,
                SCATTERING,
                ABSORPTION
        };
        
        const StructureObject *mpStructureObject;

        vector<complex_LC> *mpRefrIndex;
                
        V2<float>* mRoughness;


    vector<double>* mpQScat;
        vector<double>* mpQAbs;
        vector<double>* mpQAbsFluo;
        vector<double> mpQExt;   //No point in having this as a pointer since it is not passed into constructor.
        Distributions::Distribution *mpDistScat;
        Occurrence CheckFreePath(const V3<double>& rPrevPos, V3<double>& rCurPos, int NbrLambda, double* pT) const;

        int ScatterEnvelope(V3<double>& rDirection, V3<double>& rPolarization, double Lambda, complex_LC N1,  complex_LC N2, const V3<double>& rNormal) const;

        /*
         * Scatter direction.
         *
         * Estimates new direction after a scattering according to the distribution in m_dist_scat.
         *
         * @param rDir Current direction vector.
         * @param rPol Current polarization vector.
         * @param NbrLambda Wavelength number.
         *
         * @return void
         */
        void ScatterDirection(V3<double>& rDir, V3<double>& rPol, int NbrLambda = 0) const;
        virtual int WhatWillHappenNext(double Depth, int ComponentNumber) const;

//      Distributions::Distribution* mpDistScat;
        double mContactReduction;
        
        long mIterationLimit;

        bool mAbsFluoOrNot;
        

}; // Component

#endif // COMPONENT_H