OpenPaperOpt/WavePacket.h

Go to the documentation of this file.

// WavePacket.h: interface for the WavePacket class.

#if !defined(AFX_WAVEPACKET_H__CE436FB5_8CBD_11D5_803C_00609705562A__INCLUDED_)
#define AFX_WAVEPACKET_H__CE436FB5_8CBD_11D5_803C_00609705562A__INCLUDED_

#if _MSC_VER >= 1000
#pragma once
#pragma warning(disable:4005)

#endif // _MSC_VER >= 1000
#include <math.h>
#include <list>
#include <iostream>
#include <fstream>
#include <string>

#include "V3.h"
#include "V2.h"
using namespace std;

//----------------------------------------------------------------------------------------

#define INCOMING -9999
#define ABSORBED 9999
#define REF_SURF 10
#define REF  11
#define TRA 12 
#define ERROR -10   

//#define REF_SURF 11

class WavePacket  
{
public:
        
        WavePacket() {mLen=0;mNumRefr=0;mNumScat=0;obj_num=INCOMING;mRefrInd=1;}
        //TODO check the refr index of wp if coming from other medium
        
        WavePacket& operator=(const WavePacket& wp){
                //mState = wp.mState;
                mLen =  wp.mLen;
                mNumRefr = wp.mNumRefr;
                mNumScat = wp.mNumScat;
                obj_num = wp.obj_num;
                mPos = wp.mPos;
                mDir = wp.mDir;
                mPol = wp.mPol;
                mRefrInd = wp.mRefrInd;
                return *this;
        }
        
        void WpError(std::string message){

                this->obj_num =ERROR;
                this->mErrorType = message;
        }
        
        // newPos is assumed to be in the positive travelling direction
        // -- no checks are made.
        void mTravelToPosition(V3 <double>& newPos) {
                
                V3 <double> temp = newPos-mPos;
                mLen += temp.LengthVector() * mRefrInd;
                mLastZpos = mPos.Z();
                mPos = newPos;
        }

        //TODO should this be a distribution?
        void mIsotropicScattering();

        void mLambertianReflectance();

        // Sposition is assumed to be in the positive travelling
        // direction -- no checks are made.
        void mTravelToSposition() {
                
                V3 <double> temp = (mSposition-mPos);
                mLen += temp.LengthVector() * mRefrInd;
                mPos = V3 <double> (mSposition);
        }
        
        // One can use t to increment mLen. However, don't
        // forget to update the other travelling functions as well!
        int mTravel(double t){
                if(t<0){
                        WpError("t negative in mTravel");
                        return(0);
                }
                mLastZpos = mPos.Z();
                mPos.X(mPos.X()+t*mDir.X());
                mPos.Y(mPos.Y()+t*mDir.Y());
                mPos.Z(mPos.Z()+t*mDir.Z());
                mLen += t * mRefrInd;
                return(1);
        }
        void mCalcPolSp();
        V3 <double> mPos;
        V3 <double> mDir;
        V3 <double> mPol;
        V3 <double> mSnormal;
        V3 <double> mSposition;
        string mErrorType;
        
        double mLen;            //traveled path 
        double mEnergy;         //relative energy of the wave packet
        double mLastZpos;       //last z-positio before leaving structure
        int   mNumRefr;         //Number of refractions and reflections
        int mNumScat;       //Number of scattering events in the bulks
        int   obj_num;      //object where the wp is
        float mLambda;

        int nbr_wp;                             //the processing number of the wp in the main loop
        int nbr_lambda;
        double mRefrInd;                //index of the medium where the wp is 

        //TODO: make sure mRefrInd is updated in all objects and components. Check also that mLen is always updated
        double mPpol,mSpol;     // polarisation

        bool hasfluored; //set to 1 when wp has been fluoresced
        int nbrFluo; // LC temp tbr 06/01/10
        int nbrAbs;
};      
#endif // !defined(AFX_WAVEPACKET_H__CE436FB5_8CBD_11D5_803C_00609705562A__INCLUDED_)