OpenPaperOpt/Fiber.h

Go to the documentation of this file.

#pragma once
#ifndef FIBER_H
#define FIBER_H

// INCLUDES
//
#include "Distribution.h"
#include "WavePacket.h"
#include "Transform.h"
#include "Basesheet.h"
#include "Component.h"
#include <string>
#include "FiberGenerator.h"
#include "Fluorescence.h"

class FiberGenerator;

class Fiber : public Component
{
public:
        Fiber(vector<complex_LC>& rRefrIndex, vector<complex_LC>& rLumenRefrIndex, vector<double>& rQScat, 
                vector<double>& rQAbs,vector<double>& rQAbsFluo, double ContactReduction, V2<float>* pRoughness,
                Distributions::Distribution* pDistScat, double Theta, double Phi,
                double PolAngle, double Ellipticity, double Length, 
                double Width, double WallThickness, const Basesheet* const pBasesheet);

        ~Fiber(void);

        int     Process(WavePacket& wp, int ComponentNumber, bool EntryScattering = false, complex_LC *pPrevRefrIndex = 0) const;
        int     Process(WavePacket& wp) const; //ML: Temp fix, this function must be defined because it is virtual
                                            //    in Component.

// ACCESS
        const V3<double>& GetDirectionVector() const { return this->mDirectionVector; }
        const double GetPolAngle() const { return this->mPolAngle; }
        const double GetEllipticity() const { return this->mEllipticity; }
        const double GetWidth() const { return this->mWidth; }
        void SetIterationLimit(long N) { mIterationLimit = N; }

private:
        enum Position
        {
                OUTER_ENVELOPE,
                INNER_ENVELOPE,
                LUMEN,
                FIBERWALL,
                OUTSIDE
        };


//      void ReadParamFile(std::string FilenameFiberwall, std::string FilenameLumen);
        Position InterceptFiberwall(const V3<double>& rDirection, V3<double>& rPosition, V3<double>& rNormal) const;
        void CreateFiber(WavePacket& wp, V3<double>& rLocalPos, V3<double>& rLocalDir, V3<double>& rLocalPol, V3<double>& rNormal) const;
        Position InterceptLumen(const V3<double>& rDirection, V3<double>& rPosition, V3<double>& rNormal) const;
        int CheckFiberends(const V3<double>& rPos1, V3<double>& rPos2) const;

        /* 
         * @param mpLumenRefrIndex      Lumen complex refractive index.
         * @param mB                            B Fiber b-axis size                             
         * @param mGlobalToLocal        
         * @param mDirectionVector      Local normal wave packet relative direction vector
         * @param mTheta                        Polar angle  of fiber.
         * @param mPhi                          Azimuthal angle  of fiber.
         * @param mPolAngle                     Tilt of fiber.
         * @param mWidth                        Fiber width.
         * @param mWallThickness        Fiberwall thickness in proportion to fiber width. Should be in the range [0, 1].
         * @param mIterationLimit       Iteration limit of fiber
*/
         vector<complex_LC> *mpLumenRefrIndex;
        double mTheta;
        double mPhi;
        double mPolAngle;
        double mEllipticity;
        double mLength;
        double mWidth;
        double mWallThickness;
        double mB;
        V3<double> mDirectionVector;

        //Used to change vector bases
        Transform mGlobalToLocal;
        
}; // Fiber

#endif // FIBER_H