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

#ifndef spectrometersH
#define spectrometersH

#include "filters.h"
#include "procObj.h"
#include "objIcon.h"
#include "input.h"
#include "classDefs.h"
#include "tasksAndRoutines.h"
#include <valarray.h>
#include <Classes.hpp>

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



class spectrometer;

class spectrum;

class spectrometerIcon : public objIcon
{
        public:
        spectrometerIcon(spectrometer *owner,int h,int w);
        spectrometer *owner;
        void __fastcall initConns(bool b);
        virtual void __fastcall draw();
};

class initSpectrometer : public routine
{
        public:
        initSpectrometer(spectrometer *Owner);
        spectrometer *owner;
        virtual void __fastcall Execute();
};



class spectrometer : public input
{
        protected:

        public:
                spectroscopySuite *spectroscopyTool;        
                virtual int __fastcall connect()=0;
                spectrometer(process *owner,String name,int addr,int samplingInterval,int h,int w,input *tempSensor);
                input *tempSensor;
                void __fastcall (__closure *onScan)();
                virtual void __fastcall initSpectra();
                void __fastcall makeIcon(int h,int w);
                String storedSampleRef;
                String storedDarkRef;
                String storedLightRef;

                virtual void __fastcall loadSample(String filename)
                {
                        if(loadSpectrum(filename,sampleRef))
                        {
                                storedSampleRef = filename;
                        }
                }

                bool INITED;
                virtual void __fastcall getPixelConv(){};
                bool __fastcall loadSpectrum(String filename,spectrum *s);
                bool __fastcall loadSpectrumFromString(String input,spectrum *s);                
                spectrum *__fastcall getRelAbs();
                spectrum *relTrans;
                spectrum *relAbs;

                spectrum *darkRef;
                spectrum *sample;
                spectrum *sampleRef;
                const unsigned int ARRSIZE;

                unsigned short refAvg;
                unsigned short sampleAvg;
                virtual void __fastcall setArrSize(const unsigned int size);
                __property unsigned int arrSize={read=ARRSIZE,write=setArrSize};
                unsigned short iTime;
                virtual void __fastcall getDarkRef()=0;
                virtual void __fastcall getSampleRef()=0;
                virtual int __fastcall avgSpec(spectrum *s, int points)=0;
                virtual double __fastcall pixelToWavelength(int pixel)=0;
                virtual double __fastcall wavelengthToPixel(double wl)=0;

                double saturationThreshold;
                DynamicArray<spectrum *>spectra;
                unsigned int nextReferenceTime;
                unsigned int referenceInterval;
                virtual void __fastcall rereference()
                {
                        getDarkRef();
                }

                String refFileName;
                void __fastcall logReference(spectrum *ref);
                double maxY;
                virtual void __fastcall setConnected(bool state)=0;
                virtual bool __fastcall getConnected()=0;

                void __fastcall (__closure*setInited)(bool);
                __property bool initialized={read=getInited};
                virtual void __fastcall initialize()=0;
                routine *initSpectrometer;

                virtual int __fastcall getSpectrum()=0;
//
//                virtual void __fastcall getState();
                bool CONNECTED;
                __property bool connected={read=getConnected,write=setConnected};
                TprogressForm *progForm;
                virtual void __fastcall makeNodes(TTreeNode *parent){}

                bool *initCriteria;
                unsigned int numInitCriteria;

                virtual bool __fastcall getInited()
                {
                        if(!INITED)
                        {
                                bool result=true;
                                int i=-1;
                                while(++i<numInitCriteria && result)
                                {
                                        result &=initCriteria[i];
                                }
                                INITED = result;
                                if(setInited!=NULL)
                                {
                                        setInited(result);
                                }
                        }
                        return INITED;
                }
};



#endif

//---------------------------------------------------------------------------
#include "objConn.h"
#include "utils.h"
#include "dir.h"
#include "procMenu.h"
#include "objTasks.h"
#include "paramClasses.h"
#include "processClasses.h"
#include "spectroscopyTools.h"
#include "progForm.h"
#include "appMain.h"
#include "spectrum.h"

#include "flowPath.h"
#pragma hdrstop

#pragma package(smart_init)
#include "spectrometers.h"

initSpectrometer::initSpectrometer(spectrometer *Owner) : routine("Initialize "+Owner->name,Owner->owner)
{
        owner = Owner;
}

void __fastcall initSpectrometer::Execute()
{
        while(!Terminated)
        {
                owner->getDarkRef();  
                owner->getSampleRef();
        }
}


bool __fastcall spectrometer::loadSpectrum(String fname,spectrum *s)
{
        double test;
        char l[65536];
        FILE *f = fopen(fname.c_str(),"r");
        DynamicArray<String>line;
        bool ok = true;
        while(fgets(l,65536,f)&&ok)
        {
                line = explode((String)l,",");
                for(int i =0;i<line.Length;i++)
                {
                        String str =Trim(line[i]);
                        if(!TryStrToFloat(str,test))
                        {
                                MessageDlg("Invalid Data Format",mtError,TMsgDlgButtons()<<mbOK,0);
                                ok=false;
                                i=line.Length;
                        }
                }
        }
        fclose(f);
        if(ok)
        {
                f = fopen(fname.c_str(),"r");
                DynamicArray<String>line;
                int d=0;
                while(fgets(l,999999,f)&&ok&&d<arrSize)
                {
                        line = explode((String)l,",");
                        for(int i =0;i<line.Length;i++)
                        {
                                double val;
                                String str =Trim(line[i]);
                                if(!TryStrToFloat(Trim(line[i]),val))
                                {
                                        MessageDlg("Invalid Data Format",mtError,TMsgDlgButtons()<<mbOK,0);
                                        ok=false;
                                        i=line.Length;
                                }
                                else
                                {
                                        s->data[d] = val;
                                }
                                d++;
                        }
                }
        }
        return ok;
}

bool __fastcall spectrometer::loadSpectrumFromString(String input,spectrum *s)
{
        double test;
        bool ok;
        int d =0;
        DynamicArray<String>line;
        line = explode((String)input,",");        
        if(line.Length!=arrSize)
        {
                MessageDlg("Array length mismatch - cannot load spectrum for "+s->title,mtError,TMsgDlgButtons()<<mbOK,0);
                return false;
        }
        for(int i =0;i<line.Length;i++)
        {
                double val;
                String str =Trim(line[i]);
                if(!TryStrToFloat(Trim(line[i]),val))
                {
                        MessageDlg("Invalid Data Format",mtError,TMsgDlgButtons()<<mbOK,0);
                        ok=false;
                        i=line.Length;
                }
                else
                {
                        s->data[d] = val;
                }
                d++;
        }
        return ok;
}


spectrometer::spectrometer(process *owner,String name,int addr,int samplingInterval,int h,int w,input *TempSensor): input(owner,name,addr,samplingInterval),
INITED(false),
onScan(NULL),
numInitCriteria(0),
tempSensor(TempSensor)
{
        arrSize=1;
        relTrans= new spectrum(this,arrSize,"Relative Transmission");
        relAbs= new spectrum(this,arrSize,"Relative Absorbance");
        sample= new spectrum(this,arrSize,"Current Spectrum");
        darkRef= new spectrum(this,arrSize,"Dark Reference");
        sampleRef= new spectrum(this,arrSize,"Sample Reference");
        CONNECTED=false;
        refFileName = name+"\ References.csv";
        initSpectra();
        if((h+w)>0)
        {
                makeIcon(h,w);
                owner->syncUpdateDisplays(this);
        }
}


void __fastcall spectrometer::makeIcon(int h,int w)
{
        icon=new spectrometerIcon(this,h,w);
        loadPopup();
        draw();        
}

void __fastcall spectrometer::logReference(spectrum *ref)
{
        String timeString;
        DateTimeToString(timeString,"yyyy-mm-dd hh:nn:ss.zz",Now());
        writeToFile(refFileName,ref->title+" "+timeString+"\n","a+");
        writeToFile(refFileName,ref->asString()+"\n","a+");
}

/*
void __fastcall spectrometer::getState()
{
        if(CONNECTED)
        {
                getSpectrum();
        }
}
*/



void __fastcall spectrometer::setArrSize(const unsigned int size)
{
        if(size!=ARRSIZE)
        {
                *(int *)&ARRSIZE=size;
                initSpectra();
        }
}

void __fastcall spectrometer::initSpectra()
{

        for(int i=0;i<spectra.Length;++i)
        {
                spectra[i]->size=arrSize;
                delete[] spectra[i]->data;
                spectra[i]->data=new double[arrSize];
                spectra[i]->fill(0);
        }
}


spectrometerIcon::spectrometerIcon(spectrometer *Owner,int h,int w) : objIcon(Owner,h,w)
{
        owner = Owner;
        initConns(NULL);
}

void __fastcall spectrometerIcon::initConns(bool b)
{
        conns.Length = 4;
        conns[0]=new objConn(owner,Width/2,0);
        conns[1]=new objConn(owner,Width,Height/2);
        conns[2]=new objConn(owner,Width/2,Height);
        conns[3]=new objConn(owner,0,Height/2);
}

void __fastcall spectrometerIcon::draw()
{
        clear();
        int th = Canvas->TextHeight(owner->name);
        int tw = Canvas->TextWidth(owner->name);

        int titleheight=0,specheight=0,border=2,specWidth=Width-2*border;
        int tpos = (Width-tw)/border;
        if(Height/4>2*border+th)
        {
                titleheight=2*border+th;
                specheight=Height-titleheight-2*border;
        }


        Canvas->Brush->Color = clSilver;
        Canvas->FillRect(TRect(0,0,Width,Height));
        Canvas->Pen->Color = 0x00FFFFFF;

        Canvas->MoveTo(0,Height);
        Canvas->LineTo(0,0);
        Canvas->LineTo(Width-1,0);

        Canvas->Pen->Color = clDkGray;
        Canvas->LineTo(Width-1,Height-1);
        Canvas->LineTo(0,Height-1);
        Canvas->Pen->Color = clMedGray;
        Canvas->MoveTo(3,Height-border);
        Canvas->LineTo(Width-border,Height-border);
        Canvas->LineTo(Width-border,1);
        Canvas->Brush->Color = 0x00FFFFFF;

        Canvas->RoundRect(tpos-8,border,tpos+tw+8,titleheight-border,8,8);
        Canvas->Brush->Style = bsClear;
        Canvas->TextOut(tpos,(titleheight-th)/border,owner->name);
        Canvas->Brush->Style = bsSolid;
        int xo = 2*border;
        Canvas->FillRect(TRect(xo,titleheight,Width-xo,Height-4));

        Canvas->Pen->Color = clMedGray;
        Canvas->MoveTo(xo,Height-3);
        Canvas->LineTo(xo,titleheight);
        Canvas->LineTo(Width-xo,titleheight);
        if(owner->CONNECTED)
        {
                int r = owner->arrSize/specWidth;
                int h= specheight-4;
                Canvas->MoveTo(4,Height-3);
                Canvas->Pen->Color = 0x00FF0000;
                for(int i=0;i<specWidth;i++)
                {
                        int v1 = owner->sample->data[(int)r*i]*h;
                        Canvas->LineTo(i+xo+1,Height-4-v1);
                }
        }
        else
        {
                TColor color = Canvas->Font->Color;
                Canvas->Font->Color=0x000000DF;
                Canvas->Font->Style = TFontStyles()<< fsBold;
                String msg = "NOT CONNECTED";
                int tw = Canvas->TextWidth(msg);
                int th = Canvas->TextHeight(msg);
                Canvas->TextOut(xo+(Width-border*xo-tw)/2,(specheight-1+th)/2,msg);
                Canvas->Font->Color=color;
                Canvas->Font->Style = TFontStyles();
        }
}

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

#ifndef BWTEKUSBH
#define BWTEKUSBH
#include "spectrometers.h"

class bwTekUSB : public spectrometer
{
        protected:
        String dllPath;
        public:
        virtual int __fastcall bwtest();
        virtual double __fastcall pixelToWavelength(int pixel);
        virtual double __fastcall wavelengthToPixel(double wl);
        virtual int __fastcall avgSpec(spectrum *s, int points);

        typedef int _stdcall (*bwtekSetupChannel)(int channelNumber,int unknown);
        virtual void __fastcall getPixelConv();
        DynamicArray<double>pixelToWavelengthCoeffs;
        DynamicArray<double>wavelengthToPixelCoeffs;
        typedef int _stdcall (*bwtekTestUSB)
        (
                int nUSBTiming, // USB Interface timing option
                int nPixelNo, // number of pixels of a detector to be readout
                int nInputMode, // signal conditioning stage gain value
                int nChannel, // channel to get data from
                int pParam // extra setting parameters only for RS232
        );

        typedef long _stdcall (*bwtekSetTimeUSB)
        (
                long lTime, // integration time setting
                int nChannel // channel to get data from
        );
        typedef int _stdcall (*bwtekDataReadUSB)
        (
                int nTriggerMode, // Base address for plug-in data acquisition board
                unsigned short *pArray, // data value array from the read operation stored
                int nChannel // channel to get data from
        );

        typedef int _stdcall (*bwtekReadResultUSB)
        (
                int nTriggerMode, // setting for trigger mode
                int nAverage, // number of scans to be averaged
                int nTypeSmoothing, // setting for smoothing type
                int nValueSmoothing, // setting for smoothing parameter
                unsigned short *pArray, // an array memory space for scan data
                int nChannel=0 // setting for channel number
        );

        typedef int _stdcall (*bwtekCloseUSB)
        (
                int nchannel // channel to get data from
        );


        typedef int _stdcall (*bwtekReadEEPROMUSB)
        (
                char *OutFileName, // The filename in which data from EEPROM would be saved
                int nChannel=0 // channel number to get data from
        );

        typedef int _stdcall (*bwtekSetTimingsUSB)
        (
                long lTriggerExit, // Setting external trigger timeout
                int nMultiple, // A multifly factor for long integration time
                int nChannel=0 //channel number to get data from
        );

        typedef int __stdcall (* bwtekPolyFit)
        (
        	double *x, // Array of independent variables
        	double *y, // Array of dependent variables
        	int const numPts, // Number of points in independent and dependent arrays
        	double *coefs, // Pointer to array containing calculated coeffi cients [index from 0 to order]
        	int const order // Desired order of polynomial fi t
        );

        typedef void _stdcall (*bwtekPolyCalc)
        (
        	double *coefs,
        	int const order,
        	int const x,
        	double *y
        );

        typedef int _stdcall (*bwtekGetExtStatus)
        (
        	int nchannel //channel to get data from
        );

        typedef int _stdcall (*bwtekSetExtLaser)
        (
        	int nOnOff, //switch for On/Off
        	int nchannel //channel to get data from
        );

        typedef int _stdcall (*bwtekSetExtShutter)
        (
        	int nOnOff, // switch for On/Off
        	int nchannel // channel to get data from
        );

        typedef int _stdcall (*bwtekSetExtSync)
        (
        	int nOnOff, //switch for On/Off
        	int nchannel //channel to get data from
        );

        typedef int _stdcall (*bwtekGatedMode)
        (
        	int nGatedTime, // Gated Time between an External trigger and a CCD trigger
        	int nchannel // channel to get data from
        );

        typedef int _stdcall (*bwtekSetExtPulse)
        (
        	int nOnOff, // switch to turn on or off pulse train
        	int nDelayTime, // delay time between the CCD trigger and the fi rst pulse
        	int nHigh, // time duration for high pulse width
        	int nLow, // time duration for low pulse width
        	int nPulse, // number of pulses to be generated
        	int nInverse, // pulse output will be inverted
        	int nChannel // channel to get data from
        );
        
        virtual void __fastcall rereference();
        virtual void __fastcall maintMode(bool on);
        virtual void __fastcall loadPopup();        
        typedef int _stdcall (*bwtekShutterOpen)(int);
//        enum initializationCriteria{iTimeOK=1,darkOK=2,sampleOK=4,pixelConv=8,OK=15}inited;
        const unsigned int iTimeOK;
        const unsigned int darkOK;
        const unsigned int sampleOK;
        const unsigned int pixelConvOK;
        void __fastcall processSpectrum();
        void __fastcall checkRef();
        typedef int _stdcall (*bwtekShutterClose)(int);
        unsigned int iTimeOffset;
        virtual void __fastcall initialize();
        virtual void __fastcall autoSetITime();
        virtual void __fastcall getRelTrans();

        virtual int __fastcall getSpectrum();
        int __fastcall connect();
        void __fastcall setITime(int ITime);
        int __fastcall setExtShutter(int state);

        virtual void __fastcall getSampleRef();
        virtual void __fastcall getDarkRef();
        virtual int setShutter(unsigned short state);
        void __fastcall loadDLL(String path);
        void __fastcall initPix(const int size);
        

        bwTekUSB(String name,process *Owner,int Channel,int arrSize ,int ADCTYPE,int TimingMode,bool connectNow,input *tempSensor,int h,int w);
        int _stdcall (*bwTest)
        (
                int nUSBTiming, // USB Interface timing option
                int nPixelNo, // number of pixels of a detector to be readout
                int nInputMode, // signal conditioning stage gain value
                int nChannel, // channel to get data from
                int pParam // extra setting parameters only for RS232
        );

        int _stdcall (*readEeprom)
        (
                char *OutFileName, // The filename in which data from EEPROM would be saved
                int nChannel=0 // channel number to get data from
        );

        int _stdcall (*setupChannel)(int channel,int unknown);
        int _stdcall (*openShutter)(int);

        int _stdcall (*closeShutter)(int);

        long _stdcall (*bwSetTime)
        (
                long lTime, // integration time setting
                int nChannel // channel to get data from
        );

        int _stdcall (*bwReadData)
        (
                int nTriggerMode, // Base address for plug-in data acquisition board
                unsigned short *pArray, // data value array from the read operation stored
                int nChannel // channel to get data from
        );

        int _stdcall (*bwClose)
        (
                int nchannel // channel to get data from
        );

        int _stdcall (*bwReadAvgResult)
        (
                int nTriggerMode, // setting for trigger mode
                int nAverage, // number of scans to be averaged
                int nTypeSmoothing, // setting for smoothing type
                int nValueSmoothing, // setting for smoothing parameter
                unsigned short *pArray, // an array memory space for scan data
                int nChannel=0 // setting for channel number
        );

        int _stdcall (*bwReadEEPROM)
        (
                char *OutFileName, // The filename in which data from EEPROM would be saved
                int nChannel=0 // channel number to get data from
        );

        int _stdcall (*bwSetTiming)
        (
                long lTriggerExit, // Setting external trigger timeout
                int nMultiple, // A multifly factor for long integration time
                int nChannel=0 //channel number to get data from
        );

        int _stdcall (*bwPolyFit)
        (
        	double *x, // Array of independent variables
        	double *y, // Array of dependent variables
        	int const numPts, // Number of points in independent and dependent arrays
        	double *coefs, // Pointer to array containing calculated coeffi cients [index from 0 to order]
        	int const order // Desired order of polynomial fi t
        );

        void _stdcall (*bwPolyCalc)
        (
        	double *coefs,
        	int const order,
        	int const x,
        	double *y
        );

         int _stdcall (*bwGetStatus)
        (
        	int nchannel //channel to get data from
        );

         int _stdcall (*bwSetExtLaser)
        (
        	int nOnOff, //switch for On/Off
        	int nchannel //channel to get data from
        );

         int _stdcall (*bwSetExtShutter)
        (
        	int nOnOff, // switch for On/Off
        	int nchannel // channel to get data from
        );

         int _stdcall (*bwSetExtSync)
        (
        	int nOnOff, //switch for On/Off
        	int nchannel //channel to get data from
        );

         int _stdcall (*bwGatedMode)
        (
        	int nGatedTime, // Gated Time between an External trigger and a CCD trigger
        	int nchannel // channel to get data from
        );

         int _stdcall (*bwSetPulse)
        (
        	int nOnOff, // switch to turn on or off pulse train
        	int nDelayTime, // delay time between the CCD trigger and the fi rst pulse
        	int nHigh, // time duration for high pulse width
        	int nLow, // time duration for low pulse width
        	int nPulse, // number of pulses to be generated
                int nInverse, // pulse output will be inverted
	        int nChannel // channel to get data from
        );
        void __fastcall scanOnTimer(TObject *Sender);

        unsigned short shutterState;
        unsigned short timingMode;
        unsigned short adcType;

        double *tmp;
        unsigned short *rawData;
        double maxY;

        void __fastcall setChannel(int chan);
        virtual bool __fastcall getInited()
        {
                if(!INITED)
                {
                        bool result=true;
                        int i=0;
                        while(i<numInitCriteria && result)
                        {
                                result &=initCriteria[i];
                                i++;
                        }
                        INITED = result;
                        if(setInited!=NULL)
                        {
                                setInited(result);
                        }

                }
                return INITED;
        }
        void __fastcall (__closure*setInited)(bool);
        __property bool initialized={read=getInited};

        virtual void __fastcall setConnected(bool state);
        virtual bool __fastcall getConnected();
        __property bool connected={read=getConnected,write=setConnected};
        unsigned short shutterWait;
        void __fastcall refreshIcon();
};
//---------------------------------------------------------------------------
class bwTek261E_512 : public bwTekUSB
{
        public:
        bwTek261E_512(String name,process *Owner,int chan,bool connectNow,input *tempSensor,int h, int w):bwTekUSB(name,Owner,chan,512,13,3,connectNow,tempSensor,h,w)
        {
                HINSTANCE HInst=LoadLibrary("BWTEKUSB.dll");
                if( HInst>0 )
                {
                        getMode=(bwtekGetInGaAsMode)GetProcAddress(HInst,"bwtekGetInGaAsMode");
                        setMode=(bwtekSetInGaAsMode) GetProcAddress(HInst,"bwtekSetInGaAsMode");
                }
        }

        typedef int _stdcall (*bwtekGetInGaAsMode)
        (
                int *nMode, // select detector mode
                int nChannel // channel to get data from
        );
        typedef int _stdcall (*bwtekSetInGaAsMode)
        (
                int nMode, // select feedback capacitor mode
                int nChannel // channel to get data from
        );
        int _stdcall (*getMode)
        (
                int *nMode, // select detector mode
                int nChannel // channel to get data from
        );
        int _stdcall (*setMode)
        (
                int nMode, // select feedback capacitor mode
                int nChannel // channel to get data from
        );
};

class initBwTekS255_256 : public initSpectrometer
{
        public:
        initBwTekS255_256(bwTekS255_256 *Owner);
        bwTekS255_256 *owner;
        virtual void __fastcall Execute();
};

class bwTekS255_256 : public bwTekUSB
{
        public:
        bwTekS255_256(process *Owner,String name, int chan,bool connectNow,input *tempSensor,int h,int w);
        typedef int _stdcall (*bwtekShutterControl)
        (
                int sampleShutterState,
                int referenceShutterState,
                int channel
        );

        void __fastcall loadDark(String filename)
        {
                if(loadSpectrum(filename,dark0))
                {
                        storedDarkRef = filename;
                        firstDark = false;
                }
        }

        void __fastcall loadLight(String filename)
        {
                if(loadSpectrum(filename,ref0))
                {
                        storedLightRef = filename;
                        firstLight = false;                        
                }
        }
        virtual int __fastcall bwtest();
        int _stdcall (*bwtekSetShutters)(int,int,int);
        virtual int setShutter(unsigned short state);
        unsigned short refShutterState;
        virtual void __fastcall loadPopup();
        void __fastcall setRefShutter(unsigned short state);
        virtual void __fastcall autoSetITime();
        virtual void __fastcall getDarkRef();
        virtual void __fastcall getLightRef();
        virtual void __fastcall getSampleRef();
        virtual void __fastcall getRelTrans();
        virtual void __fastcall initialize();
        bool firstDark;
        bool firstLight;
        spectrum *ref0;
        spectrum *dark0;        
        spectrum *lastRef;
        virtual void __fastcall rereference();
};

//---------------------------------------------------------------------------
#endif

//---------------------------------------------------------------------------
#include "vclUtils.h"
#include "permissionSet.h"
#include "utils.h"
#include "dir.h"
#include "procMenu.h"
#include "objTasks.h"
#include "paramClasses.h"
#include "processClasses.h"
#include "spectroscopyTools.h"
#include "progForm.h"
#include "spectrum.h"
#include "appMain.h"
#include "flowPath.h"
#include "tasksAndRoutines.h"
#pragma hdrstop

#include "BWTEKUSB.h"
#pragma package(smart_init)

bwTekUSB::bwTekUSB(String name,process *Owner,int Channel,int arrSize ,int ADCTYPE,int TimingMode,bool connectNow,input *tempSensor,int h, int w) :
        spectrometer(Owner,name,Channel,iTime+10,h,w,tempSensor),
        iTimeOK(0),
        darkOK(1),
        sampleOK(2),
        pixelConvOK(3)
{

        get->addTask(new vvPtrTask(name+" Check reference time",&checkRef));
        int *returnCode = new int();
        get->addTask(new avPtrTask<int>(name+" Get Spectrum",&getSpectrum,returnCode));
        get->addTask(
                new compare<int>(
                        name+" compare return code",
                        returnCode,">=",new int(0),
                        new vvPtrTask(name+" processSpectrum",&processSpectrum),
                        NULL
                )
        );


        numInitCriteria = 4;
        initCriteria = new bool[numInitCriteria];
        for(int i=0;i<numInitCriteria;i++)
        {
                initCriteria[i]=false;
        }
        maxY=65535;
        progForm = Owner->progForm;

        connected=false;
        int i=0;

        referenceInterval = 3600*1000*4;
        nextReferenceTime = GetTickCount()+referenceInterval;

        shutterWait = 1000;

        saturationThreshold = 0.91552734375;
//        pathnode->state = NULL;

        refAvg = 10;
        iTime = 30;
        iTimeOffset = 500;
        sampleInterval = iTime+iTimeOffset;
        shutterState =0;
        adcType = ADCTYPE;
        timingMode = TimingMode;
        initPix(arrSize);
        dllPath = "BWTEKUSB.dll";
        loadDLL(dllPath);
        if(connectNow)
        {
                connect();
        }

        intParam *ip = new intParam((int *)&refAvg,"Spectra to average in measurement",NULL,NULL,NULL);
        params->add(ip);

        ip = new intParam((int *)&referenceInterval,"Time between spectral references",NULL,NULL,NULL);
        params->add(ip);

        ip = new intParam((int *)&iTime,"Integration Time",new String("ms"),NULL,NULL);
        ip->setFunc = &setITime;
        params->add(ip);
        params->add(new dblParam(&saturationThreshold,"Saturation Threshold",new String("%"),NULL,NULL));
        ip= new intParam(&iAddr,"Channel Number",NULL,NULL,NULL);
        ip->setFunc = &setChannel;
        params->add(ip);

//        draw();

}

void __fastcall bwTekUSB::setConnected(bool state)
{
        CONNECTED = state;
}

bool __fastcall bwTekUSB::getConnected()
{
        try
        {
                int i=setupChannel(iAddr,NULL);
                if(i<=0)
                {
                        connect();
                }
        }
        catch(...)
        {
        }
        return CONNECTED;
}

void __fastcall bwTekUSB::setChannel(int chan)
{
        try
        {
                bwClose(iAddr);
                iAddr = chan;
                connect();
        }
        catch(...)
        {
        }
}




double __fastcall bwTekUSB::pixelToWavelength(int pixel)
{
        double wl=0;
        for(int j=0;j<pixelToWavelengthCoeffs.Length;j++)
        {
                double x=pixelToWavelengthCoeffs[j];
                wl+=x*pow((double)pixel,j);
        }
        return wl;
}

double __fastcall bwTekUSB::wavelengthToPixel(double wl)
{
        double pixel=0;
        for(int j=0;j<wavelengthToPixelCoeffs.Length;j++)
        {
                double x=wavelengthToPixelCoeffs[j];
                pixel+=x*pow(wl,j);
        }
        return pixel;
}


void __fastcall bwTekUSB::getPixelConv()
{
        char cd[MAXPATH];
        getcwd(cd,MAXPATH);
        String fname = (String)cd+"\\bwtekCurrentSettings.ini";
        int result = readEeprom(fname.c_str(),iAddr);
        if(result>0)
        {
                TIniFile *f = new TIniFile(fname);
                String cname = "coefs_";
                String sname = "COMMON";
                if(f->SectionExists(sname))
                {
                        pixelToWavelengthCoeffs.Length=0;
                        bool ok = true;
                        for(int i=0;i<4;i++)
                        {
                                pixelToWavelengthCoeffs.Length++;
                                pixelToWavelengthCoeffs[i]=f->ReadFloat(sname,cname+"a"+(String)i,NULL);
                                wavelengthToPixelCoeffs.Length++;
                                wavelengthToPixelCoeffs[i]=f->ReadFloat(sname,cname+"b"+(String)i,NULL);
                                ok &= (pixelToWavelengthCoeffs[i]!=NULL && wavelengthToPixelCoeffs[i]!=NULL);
                        }
                        if(!ok)
                        {
                                MessageDlg("Pixel to wavelength coefficients not found in "+fname,mtError,TMsgDlgButtons()<<mbOK,0);
                        }
                        else
                        {
                                initCriteria[pixelConvOK]=true;
                        }
                }
                else
                {
                        MessageDlg("Section "+sname+" does not exist in "+fname,mtError,TMsgDlgButtons()<<mbOK,0);
                }        
                delete f;
        }
        else
        {
                if(result == -1)
                {
                        MessageDlg("Communication error with "+name,mtError,TMsgDlgButtons()<<mbOK,0);
                }
                if(result ==-2)
                {
                        MessageDlg("No data in "+name+" EEPROM",mtError,TMsgDlgButtons()<<mbOK,0);                
                }
        }
        bool b = initialized;
}

void __fastcall bwTekUSB::maintMode(bool on)
{

}

void __fastcall bwTekUSB::loadPopup()
{
        procObj::loadPopup();
}




void __fastcall bwTekUSB::initPix(const int size)
{
        arrSize = size;
        tmp = new double[size];
        rawData = new unsigned short[size];
        for(int i =0;i<arrSize;i++)
        {
                tmp[i]=0;
                rawData[i]=0;
        }
}

void __fastcall bwTekUSB::loadDLL(String path)
{
        HINSTANCE HInst=LoadLibrary(path.c_str());
        if( HInst>0 )
        {
                readEeprom = (bwtekReadEEPROMUSB)GetProcAddress(HInst,"bwtekReadEEPROMUSB");
                setupChannel=(bwtekSetupChannel)GetProcAddress(HInst,"bwtekSetupChannel");
                bwTest=(bwtekTestUSB)GetProcAddress(HInst,"bwtekTestUSB");
                bwSetTime=(bwtekSetTimeUSB) GetProcAddress(HInst,"bwtekSetTimeUSB");
                bwReadData=(bwtekDataReadUSB)GetProcAddress(HInst,"bwtekDataReadUSB");
                bwClose=(bwtekCloseUSB)GetProcAddress(HInst,"bwtekCloseUSB");
                bwReadAvgResult=(bwtekReadResultUSB)GetProcAddress(HInst,"bwtekReadResultUSB");
                bwPolyFit = (bwtekPolyFit)GetProcAddress(HInst,"bwtekPolyFit");
                bwPolyCalc = (bwtekPolyCalc)GetProcAddress(HInst,"bwtekPolyCalc");
                bwGetStatus = (bwtekGetExtStatus)GetProcAddress(HInst,"bwtekGetExtStatus");
                bwSetExtLaser = (bwtekSetExtLaser)GetProcAddress(HInst,"bwtekSetExtLaser");
                bwSetExtShutter = (bwtekSetExtShutter)GetProcAddress(HInst,"bwtekSetExtShutter");
                bwSetExtSync = (bwtekSetExtSync)GetProcAddress(HInst,"bwtekSetExtSync");
                bwGatedMode = (bwtekGatedMode)GetProcAddress(HInst,"bwtekGatedMode");
                bwSetPulse = (bwtekSetExtPulse)GetProcAddress(HInst,"bwtekSetExtPulse");
                closeShutter = (bwtekShutterOpen)GetProcAddress(HInst,"bwtekShutterOpen");
                openShutter = (bwtekShutterClose)GetProcAddress(HInst,"bwtekShutterClose");
                dllPath = path;
        }
        else
        {
                TFileOpen *o = new TFileOpen(owner);
                o->Caption = "Please Locate File bwtekusb.dll";
                o->Dialog->Filter ="Dynamically Linked Library (bwtekusb.dll)|BWTEKUSB.DLL";
                o->Dialog->Title = "Please Locate File bwtekusb.dll";
                open:
                if(o->Dialog->Execute())
                {
                        loadDLL(o->Dialog->FileName);
                }
                else
                {
                        goto open;
                }
                delete o;
        }
}


int bwTekUSB::setShutter(unsigned short state)
{
        int i=0;
        while(shutterState!=state && i<100)
        {
                if(state>0)
                {
                        shutterState=openShutter(iAddr);
                }
                else
                {
                        shutterState=closeShutter(iAddr);
                        shutterState^=1;
                }
        }
        if(i==100&&shutterState!=state)
        {
                shutterState=-1;
        }
        return shutterState;
}

int __fastcall bwTekUSB::avgSpec(spectrum *s, int points)
{
        if(connected)
        {
                for(int i =0;i<arrSize;i++)
                {
                        tmp[i]=0;
                }
                for(int i= 0;i<points;i++)
                {
                        bwReadData(0,rawData,iAddr);
                        for(int j=0;j<arrSize;j++)
                        {
                                tmp[j] = (double)(tmp[j]*i+rawData[j])/(double)(i+1);
                        }
                }
                unsigned short ti;
                for(int i =0;i<arrSize;i++)
                {
                        s->data[i] = (double)tmp[i]/maxY;
                }
        }
}

void __fastcall bwTekUSB::getDarkRef()
{
        if(connected)
        {
                setShutter(0);
                int then=GetTickCount()+shutterWait;
                while(GetTickCount()<then){};
                avgSpec(darkRef,refAvg);
                logReference(darkRef);
                setShutter(1);
        }
}

void __fastcall bwTekUSB::getSampleRef()
{
        if(connected)
        {
                setShutter(1);
                avgSpec(sampleRef,refAvg);
                logReference(sampleRef);
        }
}

spectrum *__fastcall spectrometer::getRelAbs()
{
        for(int i=0;i<arrSize;i++)
        {
                if(sample->data[i]==0)
                {
                        relAbs->data[i]=0;
                }
                else
                {
                        if(relTrans->data[i]==0)
                        {
                                relAbs->data[i]=65535;
                        }
                        else
                        {
                                try
                                {
                                        relAbs->data[i]=-log10(relTrans->data[i]);
                                }
                                catch(...)
                                {
                                        relAbs->data[i]=0;
                                }
                        }
                }
        }
        return relAbs;
}

int __fastcall bwTekUSB::setExtShutter(int state)
{
        if(connected)
        {
                return bwSetExtShutter(state,iAddr);
        }
}

void __fastcall bwTekUSB::setITime(int ITime)
{
        if(connected)
        {
                int result =bwSetTime(ITime,iAddr);
                if(result==ITime)
                {
                        iTime = ITime;
                }
        }
}

int __fastcall bwTekUSB::bwtest()
{
        int returnCode = bwTest(timingMode,arrSize,adcType,iAddr,NULL);
        if(returnCode>=0)
        {
                setShutter(shutterState);
        }
        return returnCode;
}

int __fastcall bwTekUSB::connect()
{
        int returnCode = bwtest();
//        int returnCode = -1;
        if(returnCode>=0)
        {
                connected = true;
        }
        else
        {
                connected = false;
        }
        return CONNECTED;
}

int __fastcall bwTekUSB::getSpectrum()
{
        if(connected)
        {
                try
                {
                        if(this->shutterState!=1)
                        {
                                setShutter(1);
                        }
                        int returnCode = bwReadData(0,rawData,iAddr);
                        if(returnCode>=0)
                        {
                                for(int i =0;i<arrSize;i++)
                                {
                                        sample->data[i]=(double)rawData[i]/65535;
                                }
                        }
                        return returnCode;
                }
                catch(...){}
        }
}


void __fastcall bwTekUSB::checkRef()
{
        if(GetTickCount()>=nextReferenceTime)
        {
                rereference();
                nextReferenceTime = GetTickCount()+referenceInterval;
        }
}
void __fastcall bwTekUSB::processSpectrum()
{
        (*measNo)++;
        getRelTrans();
        getRelAbs();
        if(onScan!=NULL)
        {
                onScan();
        }
        nextSampleTime=GetTickCount()+sampleInterval;
}

void __fastcall bwTekUSB::scanOnTimer(TObject *Sender)
{
        getState();
}

void __fastcall bwTekUSB::getRelTrans()
{
        for(int i =0;i<arrSize;i++)
        {
                if(darkRef->data[i]>=sample->data[i]||darkRef->data[i]>=sampleRef->data[i])
                {
                        relTrans->data[i]=0;
                }
                else
                {
                        float f =((sample->data[i]-darkRef->data[i]))/((sampleRef->data[i]-darkRef->data[i]));
                        if(f<1)
                        {
                                relTrans->data[i]=f;
                        }
                        else
                        {
                                relTrans->data[i]=1;
                        }
                }
        }
}

void __fastcall bwTekUSB::autoSetITime()
{
        if(connected)
        {
        if(shutterState<1)
        {
                setShutter(1);
        }
        double tmax = saturationThreshold;
        double imax =0;
        int tm=tmax*100;
        int im=imax*100;
        int di = 0;

        while(im!=tm)
        {
                di=0;
                setITime(iTime);
                avgSpec(sample,refAvg);
                for(int i =0;i<arrSize-1;i++)
                {
                        if(imax<sample->data[i])
                        {
                                imax = sample->data[i];
                        }
                }
                if(imax>tmax)
                {
                        di =-1;
                }
                if(imax<tmax)
                {
                        di = 1;
                }
                iTime+=di;
                im=imax*100;

        }
        setITime(iTime);
        sampleInterval = iTime+iTimeOffset;        
        initCriteria[iTimeOK]=true;
        bool b = initialized;
        }
}


void __fastcall bwTekUSB::initialize()
{
        if(connected)
        {
                getPixelConv();
                autoSetITime();
                getDarkRef();
                getSampleRef();
                getState();
                nextReferenceTime = GetTickCount()+referenceInterval;
                bool b = initialized;
        }
}

void __fastcall bwTekUSB::rereference()
{
        getDarkRef();
}


void __fastcall bwTekUSB::refreshIcon()
{
        int xo = (Width-128)/2;
        icon->Canvas->FillRect(TRect(xo,Height/4,Width-xo,Height-4));
        icon->Canvas->Pen->Color = clMedGray;
        icon->Canvas->MoveTo(xo,Height-3);
        icon->Canvas->LineTo(xo,Height/4);
        icon->Canvas->LineTo(Width-xo,Height/4);

        icon->Canvas->MoveTo(4,Height-3);
        icon->Canvas->Pen->Color = 0x00FF0000;
        int r = arrSize/128;
        int h= 3*Height/4;
        for(int i=0;i<128;i++)
        {
                int v1 = this->rawData[r*i]*h/maxY;
                icon->Canvas->LineTo(i+xo+1,Height-4-v1);
        }
        pathnode->draw();
}


bwTekS255_256::bwTekS255_256(process *Owner,String name, int chan,bool connectNow,input *tempSensor,int h,int w) : bwTekUSB(name,Owner,chan,256,14,3,connectNow,tempSensor,h,w)
{
        params->editable = false;
        firstLight = true;
        firstDark = true;
        fileParam *fp = new fileParam(&storedSampleRef,"Sample Reference File");
        fp->setFunc = &loadSample;
        params->add(fp);

        internalPerms->onLogin->addTask(new setVar<bool>("Set Sample Reference File Parameter Editable",&params->editable,new bool(true)));
        internalPerms->onLogout->addTask(new setVar<bool>("Set Sample Reference File Parameter Editable",&params->editable,new bool(false)));


        fp = new fileParam(&storedDarkRef,"Dark Reference File");
        fp->setFunc = &loadDark;
        params->add(fp);

        fp = new fileParam(&storedLightRef,"Light Reference File");
        fp->setFunc = &loadLight;
        params->add(fp);

        HINSTANCE HInst=LoadLibrary(dllPath.c_str());
        if( HInst>0 )
        {
                bwtekSetShutters=(bwtekShutterControl)GetProcAddress(HInst,"bwtekShutterControl");
        }
        refShutterState = 0;
        ref0=new spectrum(this,arrSize,"Original Reference Channel Spectrum");
        dark0=new spectrum(this,arrSize,"Original Dark Reference Spectrum");
        lastRef = new spectrum(this,arrSize,"Current Reference Channel Spectrum");

        spectroscopyTool = new spectroscopySuite(owner->owner,name+" Tools",this);
        MainForm->addProcess(spectroscopyTool);

        internalPerms->onLogin->addTask(new twoArgPtrTask<void,TTabSheet *,bool>("Set Spectroscopy Tab Visible",&setVisible,spectroscopyTool->tab,true));
        internalPerms->onLogout->addTask(new twoArgPtrTask<void,TTabSheet *,bool>("Set Spectroscopy Tab Visible",&setVisible,spectroscopyTool->tab,false));
        internalPerms->onLogin->addTask(new twoArgPtrTask<void,TMenuItem *,bool>("Set Spectroscopy Tab Visible",&setVisible,spectroscopyTool->procMainMenu,true));
        internalPerms->onLogout->addTask(new twoArgPtrTask<void,TMenuItem *,bool>("Set Spectroscopy Tab Visible",&setVisible,spectroscopyTool->procMainMenu,false));
        internalPerms->onLogout->addTask(new twoArgPtrTask<void,TTimer *,bool>("Set Spectroscopy Tab Visible",&setEnabled,((TspectroscopyTool *)spectroscopyTool->PnID)->timer,false));

        spectroscopyTool->tab->TabVisible = false;
        spectroscopyTool->procMainMenu->Visible = false;

        initSpectrometer = new initBwTekS255_256(this);
}



int bwTekS255_256::setShutter(unsigned short state)
{
        if(connected)
        {
                try
                {
                        if(bwtekSetShutters(state,-1,iAddr))
                        {
                                shutterState = state;
                        }
                }
                catch(...)
                {
                        
                }
        }
}

void __fastcall bwTekS255_256::setRefShutter(unsigned short state)
{
        if(connected)
        {
                try
                {
                        if(bwtekSetShutters(-1,state,iAddr))
                        {
                                refShutterState =state;
                        }
                }
                catch(...)
                {
                }
        }
}

void __fastcall bwTekS255_256::autoSetITime()
{
        if(connected)
        {
                setShutter(0);
                setRefShutter(1);
                double tmax = saturationThreshold;
                double imax =0;
                int di = 0;
                double lastErr =tmax+1;
                double err=tmax;
                while(lastErr>=err&&imax!=tmax)
                {
                        imax=0;
                        di=0;
                        lastErr = err;
                        setITime(iTime);
                        avgSpec(sample,refAvg);
                        for(int i =0;i<arrSize-1;i++)
                        {
                                if(imax<sample->data[i])
                                {
                                        imax = sample->data[i];
                                }
                        }
                        if(imax>tmax)
                        {
                                di =-1;
                        }
                        if(imax<tmax)
                        {
                                di = 1;
                        }
                        err =fabs(tmax-imax);
                        if(err>lastErr)
                        {
                                di-=2*di;
                        }
                        if(imax==tmax)
                        {
                                di=0;
                        }
                        iTime+=di;
                }
                setITime(iTime);
                sampleInterval = iTime+iTimeOffset;
                initCriteria[iTimeOK]=true;
                bool b = initialized;        
        }
}

void __fastcall bwTekS255_256::getLightRef()
{
        if(connected)
        {
                setShutter(0);
                setRefShutter(1);
                int then=GetTickCount()+shutterWait;
                while(GetTickCount()<then){}
                avgSpec(lastRef,refAvg);
                logReference(lastRef);
                setRefShutter(0);
                setShutter(1);
                bool b = initialized;
                if(firstLight)
                {
                        firstLight = false;
                        for(int i=0;i<arrSize;i++)
                        {
                                ref0->data[i]=lastRef->data[i];
                        }
                }
        }
}
void __fastcall bwTekS255_256::getDarkRef()
{
        if(connected)
        {
                setShutter(0);
                setRefShutter(0);
                int then=GetTickCount()+shutterWait;
                while(GetTickCount()<then){}
                avgSpec(darkRef,refAvg);
                logReference(darkRef);
                initCriteria[darkOK]=true;
                setRefShutter(0);
                setShutter(1);
                bool b = initialized;
                if(firstDark)
                {
                        firstDark = false;
                        for(int i=0;i<arrSize;i++)
                        {
                                dark0->data[i]=darkRef->data[i];
                        }
                }
        }
}

void __fastcall bwTekS255_256::getSampleRef()
{
        if(connected)
        {
                setShutter(1);
                setRefShutter(0);
                int then=GetTickCount()+shutterWait;
                while(GetTickCount()<then){};
                avgSpec(sampleRef,refAvg);
                logReference(sampleRef);
                setRefShutter(0);
                setShutter(1);
                initCriteria[sampleOK]=true;
                bool b = initialized;

        }
}

void __fastcall bwTekS255_256::initialize()
{
        if(connected)
        {
                initSpectrometer->Resume();
        }
}

void __fastcall bwTekS255_256::getRelTrans()
{
        for(int i =0;i<arrSize;i++)
        {
                if(darkRef->data[i]>=sample->data[i]||darkRef->data[i]>=sampleRef->data[i]||darkRef->data[i]>=lastRef->data[i]||dark0->data[i]>=ref0->data[i])
                {
                        relTrans->data[i]=0;
                }
                else
                {
                        float f =((sample->data[i]-darkRef->data[i]))/((sampleRef->data[i]-dark0->data[i]))*((lastRef->data[i]-darkRef->data[i])/(ref0->data[i]-dark0->data[i]));
                        relTrans->data[i]=f;
                }
        }
}

void __fastcall bwTekS255_256::loadPopup()
{
        bwTekUSB::loadPopup();
}

void __fastcall bwTekS255_256::rereference()
{
        if(connected)
        {
                getDarkRef();
                getLightRef();
        }
}

int __fastcall bwTekS255_256::bwtest()
{
        try
        {
                int returnCode = bwTest(timingMode,arrSize,adcType,iAddr,NULL);
                if(returnCode>=0)
                {
                        setShutter(shutterState);
                        setRefShutter(refShutterState);
                }
                return returnCode;
        }
        catch(...)
        {
        }
}

initBwTekS255_256::initBwTekS255_256(bwTekS255_256 *Owner) : initSpectrometer(Owner)
{
        owner = Owner;
}

void __fastcall initBwTekS255_256::Execute()
{
        while(!Terminated)
        {

                owner->progForm->progBar->POS=0;
                owner->progForm->progBar->min=0;
                owner->progForm->progBar->max = 5;

                owner->progForm->labelCaption="Reading EEPROM";
                Synchronize(owner->progForm->updateLabelCaption);
                Synchronize(owner->getPixelConv);
                owner->progForm->progBar->POS++;
                Synchronize(owner->progForm->progBar->Paint);
                owner->progForm->labelCaption="Detecting integration time...";
                Synchronize(owner->progForm->updateLabelCaption);

                Synchronize(owner->autoSetITime);
                owner->progForm->progBar->POS++;

                owner->progForm->labelCaption="Getting dark reference...";
                Synchronize(owner->progForm->progBar->Paint);
                Synchronize(owner->progForm->updateLabelCaption);
                Synchronize(owner->getDarkRef);


                owner->progForm->labelCaption="Getting light reference...";
                Synchronize(owner->progForm->updateLabelCaption);
                owner->progForm->progBar->POS++;
                Synchronize(owner->progForm->progBar->Paint);

                Synchronize(owner->getLightRef);

                owner->progForm->labelCaption="Getting water background...";
                Synchronize(owner->progForm->updateLabelCaption);

                owner->progForm->progBar->POS++;
                Synchronize(owner->progForm->progBar->Paint);

                Synchronize(owner->getSampleRef);
                owner->progForm->progBar->POS++;
                Synchronize(owner->progForm->progBar->Paint);

                owner->progForm->labelCaption="Finishing initialization";
                Synchronize(owner->progForm->updateLabelCaption);

                Synchronize(owner->getState);
                for(int i=0;i<owner->arrSize;i++)
                {
                        owner->ref0->data[i]=owner->lastRef->data[i];
                        owner->dark0->data[i]=owner->darkRef->data[i];
                }
                owner->nextReferenceTime = GetTickCount()+owner->referenceInterval;

                owner->progForm->close();
                Suspend();
        }
}

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

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

#ifndef spectrometersH
#define spectrometersH

#include "filters.h"
#include "procObj.h"
#include "objIcon.h"
#include "input.h"
#include "classDefs.h"
#include "tasksAndRoutines.h"
#include <valarray.h>
#include <Classes.hpp>

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



class spectrometer;

class spectrum;

class spectrometerIcon : public objIcon
{
        public:
        spectrometerIcon(spectrometer *owner,int h,int w);
        spectrometer *owner;
        void __fastcall initConns(bool b);
        virtual void __fastcall draw();
};

class initSpectrometer : public routine
{
        public:
        initSpectrometer(spectrometer *Owner);
        spectrometer *owner;
        virtual void __fastcall Execute();
};



class spectrometer : public input
{
        protected:

        public:
                spectroscopySuite *spectroscopyTool;        
                virtual int __fastcall connect()=0;
                spectrometer(process *owner,String name,int addr,int samplingInterval,int h,int w,input *tempSensor);
                input *tempSensor;
                void __fastcall (__closure *onScan)();
                virtual void __fastcall initSpectra();
                void __fastcall makeIcon(int h,int w);
                String storedSampleRef;
                String storedDarkRef;
                String storedLightRef;

                virtual void __fastcall loadSample(String filename)
                {
                        if(loadSpectrum(filename,sampleRef))
                        {
                                storedSampleRef = filename;
                        }
                }

                bool INITED;
                virtual void __fastcall getPixelConv(){};
                bool __fastcall loadSpectrum(String filename,spectrum *s);
                bool __fastcall loadSpectrumFromString(String input,spectrum *s);                
                spectrum *__fastcall getRelAbs();
                spectrum *relTrans;
                spectrum *relAbs;

                spectrum *darkRef;
                spectrum *sample;
                spectrum *sampleRef;
                const unsigned int ARRSIZE;

                unsigned short refAvg;
                unsigned short sampleAvg;
                virtual void __fastcall setArrSize(const unsigned int size);
                __property unsigned int arrSize={read=ARRSIZE,write=setArrSize};
                unsigned short iTime;
                virtual void __fastcall getDarkRef()=0;
                virtual void __fastcall getSampleRef()=0;
                virtual int __fastcall avgSpec(spectrum *s, int points)=0;
                virtual double __fastcall pixelToWavelength(int pixel)=0;
                virtual double __fastcall wavelengthToPixel(double wl)=0;

                double saturationThreshold;
                DynamicArray<spectrum *>spectra;
                unsigned int nextReferenceTime;
                unsigned int referenceInterval;
                virtual void __fastcall rereference()
                {
                        getDarkRef();
                }

                String refFileName;
                void __fastcall logReference(spectrum *ref);
                double maxY;
                virtual void __fastcall setConnected(bool state)=0;
                virtual bool __fastcall getConnected()=0;

                void __fastcall (__closure*setInited)(bool);
                __property bool initialized={read=getInited};
                virtual void __fastcall initialize()=0;
                routine *initSpectrometer;

                virtual int __fastcall getSpectrum()=0;
//
//                virtual void __fastcall getState();
                bool CONNECTED;
                __property bool connected={read=getConnected,write=setConnected};
                TprogressForm *progForm;
                virtual void __fastcall makeNodes(TTreeNode *parent){}

                bool *initCriteria;
                unsigned int numInitCriteria;

                virtual bool __fastcall getInited()
                {
                        if(!INITED)
                        {
                                bool result=true;
                                int i=-1;
                                while(++i<numInitCriteria && result)
                                {
                                        result &=initCriteria[i];
                                }
                                INITED = result;
                                if(setInited!=NULL)
                                {
                                        setInited(result);
                                }
                        }
                        return INITED;
                }
};



#endif

//---------------------------------------------------------------------------
#include "objConn.h"
#include "utils.h"
#include "dir.h"
#include "procMenu.h"
#include "objTasks.h"
#include "paramClasses.h"
#include "processClasses.h"
#include "spectroscopyTools.h"
#include "progForm.h"
#include "appMain.h"
#include "spectrum.h"

#include "flowPath.h"
#pragma hdrstop

#pragma package(smart_init)
#include "spectrometers.h"

initSpectrometer::initSpectrometer(spectrometer *Owner) : routine("Initialize "+Owner->name,Owner->owner)
{
        owner = Owner;
}

void __fastcall initSpectrometer::Execute()
{
        while(!Terminated)
        {
                owner->getDarkRef();  
                owner->getSampleRef();
        }
}


bool __fastcall spectrometer::loadSpectrum(String fname,spectrum *s)
{
        double test;
        char l[65536];
        FILE *f = fopen(fname.c_str(),"r");
        DynamicArray<String>line;
        bool ok = true;
        while(fgets(l,65536,f)&&ok)
        {
                line = explode((String)l,",");
                for(int i =0;i<line.Length;i++)
                {
                        String str =Trim(line[i]);
                        if(!TryStrToFloat(str,test))
                        {
                                MessageDlg("Invalid Data Format",mtError,TMsgDlgButtons()<<mbOK,0);
                                ok=false;
                                i=line.Length;
                        }
                }
        }
        fclose(f);
        if(ok)
        {
                f = fopen(fname.c_str(),"r");
                DynamicArray<String>line;
                int d=0;
                while(fgets(l,999999,f)&&ok&&d<arrSize)
                {
                        line = explode((String)l,",");
                        for(int i =0;i<line.Length;i++)
                        {
                                double val;
                                String str =Trim(line[i]);
                                if(!TryStrToFloat(Trim(line[i]),val))
                                {
                                        MessageDlg("Invalid Data Format",mtError,TMsgDlgButtons()<<mbOK,0);
                                        ok=false;
                                        i=line.Length;
                                }
                                else
                                {
                                        s->data[d] = val;
                                }
                                d++;
                        }
                }
        }
        return ok;
}

bool __fastcall spectrometer::loadSpectrumFromString(String input,spectrum *s)
{
        double test;
        bool ok;
        int d =0;
        DynamicArray<String>line;
        line = explode((String)input,",");        
        if(line.Length!=arrSize)
        {
                MessageDlg("Array length mismatch - cannot load spectrum for "+s->title,mtError,TMsgDlgButtons()<<mbOK,0);
                return false;
        }
        for(int i =0;i<line.Length;i++)
        {
                double val;
                String str =Trim(line[i]);
                if(!TryStrToFloat(Trim(line[i]),val))
                {
                        MessageDlg("Invalid Data Format",mtError,TMsgDlgButtons()<<mbOK,0);
                        ok=false;
                        i=line.Length;
                }
                else
                {
                        s->data[d] = val;
                }
                d++;
        }
        return ok;
}


spectrometer::spectrometer(process *owner,String name,int addr,int samplingInterval,int h,int w,input *TempSensor): input(owner,name,addr,samplingInterval),
INITED(false),
onScan(NULL),
numInitCriteria(0),
tempSensor(TempSensor)
{
        arrSize=1;
        relTrans= new spectrum(this,arrSize,"Relative Transmission");
        relAbs= new spectrum(this,arrSize,"Relative Absorbance");
        sample= new spectrum(this,arrSize,"Current Spectrum");
        darkRef= new spectrum(this,arrSize,"Dark Reference");
        sampleRef= new spectrum(this,arrSize,"Sample Reference");
        CONNECTED=false;
        refFileName = name+"\ References.csv";
        initSpectra();
        if((h+w)>0)
        {
                makeIcon(h,w);
                owner->syncUpdateDisplays(this);
        }
}


void __fastcall spectrometer::makeIcon(int h,int w)
{
        icon=new spectrometerIcon(this,h,w);
        loadPopup();
        draw();        
}

void __fastcall spectrometer::logReference(spectrum *ref)
{
        String timeString;
        DateTimeToString(timeString,"yyyy-mm-dd hh:nn:ss.zz",Now());
        writeToFile(refFileName,ref->title+" "+timeString+"\n","a+");
        writeToFile(refFileName,ref->asString()+"\n","a+");
}

/*
void __fastcall spectrometer::getState()
{
        if(CONNECTED)
        {
                getSpectrum();
        }
}
*/



void __fastcall spectrometer::setArrSize(const unsigned int size)
{
        if(size!=ARRSIZE)
        {
                *(int *)&ARRSIZE=size;
                initSpectra();
        }
}

void __fastcall spectrometer::initSpectra()
{

        for(int i=0;i<spectra.Length;++i)
        {
                spectra[i]->size=arrSize;
                delete[] spectra[i]->data;
                spectra[i]->data=new double[arrSize];
                spectra[i]->fill(0);
        }
}


spectrometerIcon::spectrometerIcon(spectrometer *Owner,int h,int w) : objIcon(Owner,h,w)
{
        owner = Owner;
        initConns(NULL);
}

void __fastcall spectrometerIcon::initConns(bool b)
{
        conns.Length = 4;
        conns[0]=new objConn(owner,Width/2,0);
        conns[1]=new objConn(owner,Width,Height/2);
        conns[2]=new objConn(owner,Width/2,Height);
        conns[3]=new objConn(owner,0,Height/2);
}

void __fastcall spectrometerIcon::draw()
{
        clear();
        int th = Canvas->TextHeight(owner->name);
        int tw = Canvas->TextWidth(owner->name);

        int titleheight=0,specheight=0,border=2,specWidth=Width-2*border;
        int tpos = (Width-tw)/border;
        if(Height/4>2*border+th)
        {
                titleheight=2*border+th;
                specheight=Height-titleheight-2*border;
        }


        Canvas->Brush->Color = clSilver;
        Canvas->FillRect(TRect(0,0,Width,Height));
        Canvas->Pen->Color = 0x00FFFFFF;

        Canvas->MoveTo(0,Height);
        Canvas->LineTo(0,0);
        Canvas->LineTo(Width-1,0);

        Canvas->Pen->Color = clDkGray;
        Canvas->LineTo(Width-1,Height-1);
        Canvas->LineTo(0,Height-1);
        Canvas->Pen->Color = clMedGray;
        Canvas->MoveTo(3,Height-border);
        Canvas->LineTo(Width-border,Height-border);
        Canvas->LineTo(Width-border,1);
        Canvas->Brush->Color = 0x00FFFFFF;

        Canvas->RoundRect(tpos-8,border,tpos+tw+8,titleheight-border,8,8);
        Canvas->Brush->Style = bsClear;
        Canvas->TextOut(tpos,(titleheight-th)/border,owner->name);
        Canvas->Brush->Style = bsSolid;
        int xo = 2*border;
        Canvas->FillRect(TRect(xo,titleheight,Width-xo,Height-4));

        Canvas->Pen->Color = clMedGray;
        Canvas->MoveTo(xo,Height-3);
        Canvas->LineTo(xo,titleheight);
        Canvas->LineTo(Width-xo,titleheight);
        if(owner->CONNECTED)
        {
                int r = owner->arrSize/specWidth;
                int h= specheight-4;
                Canvas->MoveTo(4,Height-3);
                Canvas->Pen->Color = 0x00FF0000;
                for(int i=0;i<specWidth;i++)
                {
                        int v1 = owner->sample->data[(int)r*i]*h;
                        Canvas->LineTo(i+xo+1,Height-4-v1);
                }
        }
        else
        {
                TColor color = Canvas->Font->Color;
                Canvas->Font->Color=0x000000DF;
                Canvas->Font->Style = TFontStyles()<< fsBold;
                String msg = "NOT CONNECTED";
                int tw = Canvas->TextWidth(msg);
                int th = Canvas->TextHeight(msg);
                Canvas->TextOut(xo+(Width-border*xo-tw)/2,(specheight-1+th)/2,msg);
                Canvas->Font->Color=color;
                Canvas->Font->Style = TFontStyles();
        }
}

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

#ifndef BWTEKUSBH
#define BWTEKUSBH
#include "spectrometers.h"

class bwTekUSB : public spectrometer
{
        protected:
        String dllPath;
        public:
        virtual int __fastcall bwtest();
        virtual double __fastcall pixelToWavelength(int pixel);
        virtual double __fastcall wavelengthToPixel(double wl);
        virtual int __fastcall avgSpec(spectrum *s, int points);

        typedef int _stdcall (*bwtekSetupChannel)(int channelNumber,int unknown);
        virtual void __fastcall getPixelConv();
        DynamicArray<double>pixelToWavelengthCoeffs;
        DynamicArray<double>wavelengthToPixelCoeffs;
        typedef int _stdcall (*bwtekTestUSB)
        (
                int nUSBTiming, // USB Interface timing option
                int nPixelNo, // number of pixels of a detector to be readout
                int nInputMode, // signal conditioning stage gain value
                int nChannel, // channel to get data from
                int pParam // extra setting parameters only for RS232
        );

        typedef long _stdcall (*bwtekSetTimeUSB)
        (
                long lTime, // integration time setting
                int nChannel // channel to get data from
        );
        typedef int _stdcall (*bwtekDataReadUSB)
        (
                int nTriggerMode, // Base address for plug-in data acquisition board
                unsigned short *pArray, // data value array from the read operation stored
                int nChannel // channel to get data from
        );

        typedef int _stdcall (*bwtekReadResultUSB)
        (
                int nTriggerMode, // setting for trigger mode
                int nAverage, // number of scans to be averaged
                int nTypeSmoothing, // setting for smoothing type
                int nValueSmoothing, // setting for smoothing parameter
                unsigned short *pArray, // an array memory space for scan data
                int nChannel=0 // setting for channel number
        );

        typedef int _stdcall (*bwtekCloseUSB)
        (
                int nchannel // channel to get data from
        );


        typedef int _stdcall (*bwtekReadEEPROMUSB)
        (
                char *OutFileName, // The filename in which data from EEPROM would be saved
                int nChannel=0 // channel number to get data from
        );

        typedef int _stdcall (*bwtekSetTimingsUSB)
        (
                long lTriggerExit, // Setting external trigger timeout
                int nMultiple, // A multifly factor for long integration time
                int nChannel=0 //channel number to get data from
        );

        typedef int __stdcall (* bwtekPolyFit)
        (
        	double *x, // Array of independent variables
        	double *y, // Array of dependent variables
        	int const numPts, // Number of points in independent and dependent arrays
        	double *coefs, // Pointer to array containing calculated coeffi cients [index from 0 to order]
        	int const order // Desired order of polynomial fi t
        );

        typedef void _stdcall (*bwtekPolyCalc)
        (
        	double *coefs,
        	int const order,
        	int const x,
        	double *y
        );

        typedef int _stdcall (*bwtekGetExtStatus)
        (
        	int nchannel //channel to get data from
        );

        typedef int _stdcall (*bwtekSetExtLaser)
        (
        	int nOnOff, //switch for On/Off
        	int nchannel //channel to get data from
        );

        typedef int _stdcall (*bwtekSetExtShutter)
        (
        	int nOnOff, // switch for On/Off
        	int nchannel // channel to get data from
        );

        typedef int _stdcall (*bwtekSetExtSync)
        (
        	int nOnOff, //switch for On/Off
        	int nchannel //channel to get data from
        );

        typedef int _stdcall (*bwtekGatedMode)
        (
        	int nGatedTime, // Gated Time between an External trigger and a CCD trigger
        	int nchannel // channel to get data from
        );

        typedef int _stdcall (*bwtekSetExtPulse)
        (
        	int nOnOff, // switch to turn on or off pulse train
        	int nDelayTime, // delay time between the CCD trigger and the fi rst pulse
        	int nHigh, // time duration for high pulse width
        	int nLow, // time duration for low pulse width
        	int nPulse, // number of pulses to be generated
        	int nInverse, // pulse output will be inverted
        	int nChannel // channel to get data from
        );
        
        virtual void __fastcall rereference();
        virtual void __fastcall maintMode(bool on);
        virtual void __fastcall loadPopup();        
        typedef int _stdcall (*bwtekShutterOpen)(int);
//        enum initializationCriteria{iTimeOK=1,darkOK=2,sampleOK=4,pixelConv=8,OK=15}inited;
        const unsigned int iTimeOK;
        const unsigned int darkOK;
        const unsigned int sampleOK;
        const unsigned int pixelConvOK;
        void __fastcall processSpectrum();
        void __fastcall checkRef();
        typedef int _stdcall (*bwtekShutterClose)(int);
        unsigned int iTimeOffset;
        virtual void __fastcall initialize();
        virtual void __fastcall autoSetITime();
        virtual void __fastcall getRelTrans();

        virtual int __fastcall getSpectrum();
        int __fastcall connect();
        void __fastcall setITime(int ITime);
        int __fastcall setExtShutter(int state);

        virtual void __fastcall getSampleRef();
        virtual void __fastcall getDarkRef();
        virtual int setShutter(unsigned short state);
        void __fastcall loadDLL(String path);
        void __fastcall initPix(const int size);
        

        bwTekUSB(String name,process *Owner,int Channel,int arrSize ,int ADCTYPE,int TimingMode,bool connectNow,input *tempSensor,int h,int w);
        int _stdcall (*bwTest)
        (
                int nUSBTiming, // USB Interface timing option
                int nPixelNo, // number of pixels of a detector to be readout
                int nInputMode, // signal conditioning stage gain value
                int nChannel, // channel to get data from
                int pParam // extra setting parameters only for RS232
        );

        int _stdcall (*readEeprom)
        (
                char *OutFileName, // The filename in which data from EEPROM would be saved
                int nChannel=0 // channel number to get data from
        );

        int _stdcall (*setupChannel)(int channel,int unknown);
        int _stdcall (*openShutter)(int);

        int _stdcall (*closeShutter)(int);

        long _stdcall (*bwSetTime)
        (
                long lTime, // integration time setting
                int nChannel // channel to get data from
        );

        int _stdcall (*bwReadData)
        (
                int nTriggerMode, // Base address for plug-in data acquisition board
                unsigned short *pArray, // data value array from the read operation stored
                int nChannel // channel to get data from
        );

        int _stdcall (*bwClose)
        (
                int nchannel // channel to get data from
        );

        int _stdcall (*bwReadAvgResult)
        (
                int nTriggerMode, // setting for trigger mode
                int nAverage, // number of scans to be averaged
                int nTypeSmoothing, // setting for smoothing type
                int nValueSmoothing, // setting for smoothing parameter
                unsigned short *pArray, // an array memory space for scan data
                int nChannel=0 // setting for channel number
        );

        int _stdcall (*bwReadEEPROM)
        (
                char *OutFileName, // The filename in which data from EEPROM would be saved
                int nChannel=0 // channel number to get data from
        );

        int _stdcall (*bwSetTiming)
        (
                long lTriggerExit, // Setting external trigger timeout
                int nMultiple, // A multifly factor for long integration time
                int nChannel=0 //channel number to get data from
        );

        int _stdcall (*bwPolyFit)
        (
        	double *x, // Array of independent variables
        	double *y, // Array of dependent variables
        	int const numPts, // Number of points in independent and dependent arrays
        	double *coefs, // Pointer to array containing calculated coeffi cients [index from 0 to order]
        	int const order // Desired order of polynomial fi t
        );

        void _stdcall (*bwPolyCalc)
        (
        	double *coefs,
        	int const order,
        	int const x,
        	double *y
        );

         int _stdcall (*bwGetStatus)
        (
        	int nchannel //channel to get data from
        );

         int _stdcall (*bwSetExtLaser)
        (
        	int nOnOff, //switch for On/Off
        	int nchannel //channel to get data from
        );

         int _stdcall (*bwSetExtShutter)
        (
        	int nOnOff, // switch for On/Off
        	int nchannel // channel to get data from
        );

         int _stdcall (*bwSetExtSync)
        (
        	int nOnOff, //switch for On/Off
        	int nchannel //channel to get data from
        );

         int _stdcall (*bwGatedMode)
        (
        	int nGatedTime, // Gated Time between an External trigger and a CCD trigger
        	int nchannel // channel to get data from
        );

         int _stdcall (*bwSetPulse)
        (
        	int nOnOff, // switch to turn on or off pulse train
        	int nDelayTime, // delay time between the CCD trigger and the fi rst pulse
        	int nHigh, // time duration for high pulse width
        	int nLow, // time duration for low pulse width
        	int nPulse, // number of pulses to be generated
                int nInverse, // pulse output will be inverted
	        int nChannel // channel to get data from
        );
        void __fastcall scanOnTimer(TObject *Sender);

        unsigned short shutterState;
        unsigned short timingMode;
        unsigned short adcType;

        double *tmp;
        unsigned short *rawData;
        double maxY;

        void __fastcall setChannel(int chan);
        virtual bool __fastcall getInited()
        {
                if(!INITED)
                {
                        bool result=true;
                        int i=0;
                        while(i<numInitCriteria && result)
                        {
                                result &=initCriteria[i];
                                i++;
                        }
                        INITED = result;
                        if(setInited!=NULL)
                        {
                                setInited(result);
                        }

                }
                return INITED;
        }
        void __fastcall (__closure*setInited)(bool);
        __property bool initialized={read=getInited};

        virtual void __fastcall setConnected(bool state);
        virtual bool __fastcall getConnected();
        __property bool connected={read=getConnected,write=setConnected};
        unsigned short shutterWait;
        void __fastcall refreshIcon();
};
//---------------------------------------------------------------------------
class bwTek261E_512 : public bwTekUSB
{
        public:
        bwTek261E_512(String name,process *Owner,int chan,bool connectNow,input *tempSensor,int h, int w):bwTekUSB(name,Owner,chan,512,13,3,connectNow,tempSensor,h,w)
        {
                HINSTANCE HInst=LoadLibrary("BWTEKUSB.dll");
                if( HInst>0 )
                {
                        getMode=(bwtekGetInGaAsMode)GetProcAddress(HInst,"bwtekGetInGaAsMode");
                        setMode=(bwtekSetInGaAsMode) GetProcAddress(HInst,"bwtekSetInGaAsMode");
                }
        }

        typedef int _stdcall (*bwtekGetInGaAsMode)
        (
                int *nMode, // select detector mode
                int nChannel // channel to get data from
        );
        typedef int _stdcall (*bwtekSetInGaAsMode)
        (
                int nMode, // select feedback capacitor mode
                int nChannel // channel to get data from
        );
        int _stdcall (*getMode)
        (
                int *nMode, // select detector mode
                int nChannel // channel to get data from
        );
        int _stdcall (*setMode)
        (
                int nMode, // select feedback capacitor mode
                int nChannel // channel to get data from
        );
};

class initBwTekS255_256 : public initSpectrometer
{
        public:
        initBwTekS255_256(bwTekS255_256 *Owner);
        bwTekS255_256 *owner;
        virtual void __fastcall Execute();
};

class bwTekS255_256 : public bwTekUSB
{
        public:
        bwTekS255_256(process *Owner,String name, int chan,bool connectNow,input *tempSensor,int h,int w);
        typedef int _stdcall (*bwtekShutterControl)
        (
                int sampleShutterState,
                int referenceShutterState,
                int channel
        );

        void __fastcall loadDark(String filename)
        {
                if(loadSpectrum(filename,dark0))
                {
                        storedDarkRef = filename;
                        firstDark = false;
                }
        }

        void __fastcall loadLight(String filename)
        {
                if(loadSpectrum(filename,ref0))
                {
                        storedLightRef = filename;
                        firstLight = false;                        
                }
        }
        virtual int __fastcall bwtest();
        int _stdcall (*bwtekSetShutters)(int,int,int);
        virtual int setShutter(unsigned short state);
        unsigned short refShutterState;
        virtual void __fastcall loadPopup();
        void __fastcall setRefShutter(unsigned short state);
        virtual void __fastcall autoSetITime();
        virtual void __fastcall getDarkRef();
        virtual void __fastcall getLightRef();
        virtual void __fastcall getSampleRef();
        virtual void __fastcall getRelTrans();
        virtual void __fastcall initialize();
        bool firstDark;
        bool firstLight;
        spectrum *ref0;
        spectrum *dark0;        
        spectrum *lastRef;
        virtual void __fastcall rereference();
};

//---------------------------------------------------------------------------
#endif

//---------------------------------------------------------------------------
#include "vclUtils.h"
#include "permissionSet.h"
#include "utils.h"
#include "dir.h"
#include "procMenu.h"
#include "objTasks.h"
#include "paramClasses.h"
#include "processClasses.h"
#include "spectroscopyTools.h"
#include "progForm.h"
#include "spectrum.h"
#include "appMain.h"
#include "flowPath.h"
#include "tasksAndRoutines.h"
#pragma hdrstop

#include "BWTEKUSB.h"
#pragma package(smart_init)

bwTekUSB::bwTekUSB(String name,process *Owner,int Channel,int arrSize ,int ADCTYPE,int TimingMode,bool connectNow,input *tempSensor,int h, int w) :
        spectrometer(Owner,name,Channel,iTime+10,h,w,tempSensor),
        iTimeOK(0),
        darkOK(1),
        sampleOK(2),
        pixelConvOK(3)
{

        get->addTask(new vvPtrTask(name+" Check reference time",&checkRef));
        int *returnCode = new int();
        get->addTask(new avPtrTask<int>(name+" Get Spectrum",&getSpectrum,returnCode));
        get->addTask(
                new compare<int>(
                        name+" compare return code",
                        returnCode,">=",new int(0),
                        new vvPtrTask(name+" processSpectrum",&processSpectrum),
                        NULL
                )
        );


        numInitCriteria = 4;
        initCriteria = new bool[numInitCriteria];
        for(int i=0;i<numInitCriteria;i++)
        {
                initCriteria[i]=false;
        }
        maxY=65535;
        progForm = Owner->progForm;

        connected=false;
        int i=0;

        referenceInterval = 3600*1000*4;
        nextReferenceTime = GetTickCount()+referenceInterval;

        shutterWait = 1000;

        saturationThreshold = 0.91552734375;
//        pathnode->state = NULL;

        refAvg = 10;
        iTime = 30;
        iTimeOffset = 500;
        sampleInterval = iTime+iTimeOffset;
        shutterState =0;
        adcType = ADCTYPE;
        timingMode = TimingMode;
        initPix(arrSize);
        dllPath = "BWTEKUSB.dll";
        loadDLL(dllPath);
        if(connectNow)
        {
                connect();
        }

        intParam *ip = new intParam((int *)&refAvg,"Spectra to average in measurement",NULL,NULL,NULL);
        params->add(ip);

        ip = new intParam((int *)&referenceInterval,"Time between spectral references",NULL,NULL,NULL);
        params->add(ip);

        ip = new intParam((int *)&iTime,"Integration Time",new String("ms"),NULL,NULL);
        ip->setFunc = &setITime;
        params->add(ip);
        params->add(new dblParam(&saturationThreshold,"Saturation Threshold",new String("%"),NULL,NULL));
        ip= new intParam(&iAddr,"Channel Number",NULL,NULL,NULL);
        ip->setFunc = &setChannel;
        params->add(ip);

//        draw();

}

void __fastcall bwTekUSB::setConnected(bool state)
{
        CONNECTED = state;
}

bool __fastcall bwTekUSB::getConnected()
{
        try
        {
                int i=setupChannel(iAddr,NULL);
                if(i<=0)
                {
                        connect();
                }
        }
        catch(...)
        {
        }
        return CONNECTED;
}

void __fastcall bwTekUSB::setChannel(int chan)
{
        try
        {
                bwClose(iAddr);
                iAddr = chan;
                connect();
        }
        catch(...)
        {
        }
}




double __fastcall bwTekUSB::pixelToWavelength(int pixel)
{
        double wl=0;
        for(int j=0;j<pixelToWavelengthCoeffs.Length;j++)
        {
                double x=pixelToWavelengthCoeffs[j];
                wl+=x*pow((double)pixel,j);
        }
        return wl;
}

double __fastcall bwTekUSB::wavelengthToPixel(double wl)
{
        double pixel=0;
        for(int j=0;j<wavelengthToPixelCoeffs.Length;j++)
        {
                double x=wavelengthToPixelCoeffs[j];
                pixel+=x*pow(wl,j);
        }
        return pixel;
}


void __fastcall bwTekUSB::getPixelConv()
{
        char cd[MAXPATH];
        getcwd(cd,MAXPATH);
        String fname = (String)cd+"\\bwtekCurrentSettings.ini";
        int result = readEeprom(fname.c_str(),iAddr);
        if(result>0)
        {
                TIniFile *f = new TIniFile(fname);
                String cname = "coefs_";
                String sname = "COMMON";
                if(f->SectionExists(sname))
                {
                        pixelToWavelengthCoeffs.Length=0;
                        bool ok = true;
                        for(int i=0;i<4;i++)
                        {
                                pixelToWavelengthCoeffs.Length++;
                                pixelToWavelengthCoeffs[i]=f->ReadFloat(sname,cname+"a"+(String)i,NULL);
                                wavelengthToPixelCoeffs.Length++;
                                wavelengthToPixelCoeffs[i]=f->ReadFloat(sname,cname+"b"+(String)i,NULL);
                                ok &= (pixelToWavelengthCoeffs[i]!=NULL && wavelengthToPixelCoeffs[i]!=NULL);
                        }
                        if(!ok)
                        {
                                MessageDlg("Pixel to wavelength coefficients not found in "+fname,mtError,TMsgDlgButtons()<<mbOK,0);
                        }
                        else
                        {
                                initCriteria[pixelConvOK]=true;
                        }
                }
                else
                {
                        MessageDlg("Section "+sname+" does not exist in "+fname,mtError,TMsgDlgButtons()<<mbOK,0);
                }        
                delete f;
        }
        else
        {
                if(result == -1)
                {
                        MessageDlg("Communication error with "+name,mtError,TMsgDlgButtons()<<mbOK,0);
                }
                if(result ==-2)
                {
                        MessageDlg("No data in "+name+" EEPROM",mtError,TMsgDlgButtons()<<mbOK,0);                
                }
        }
        bool b = initialized;
}

void __fastcall bwTekUSB::maintMode(bool on)
{

}

void __fastcall bwTekUSB::loadPopup()
{
        procObj::loadPopup();
}




void __fastcall bwTekUSB::initPix(const int size)
{
        arrSize = size;
        tmp = new double[size];
        rawData = new unsigned short[size];
        for(int i =0;i<arrSize;i++)
        {
                tmp[i]=0;
                rawData[i]=0;
        }
}

void __fastcall bwTekUSB::loadDLL(String path)
{
        HINSTANCE HInst=LoadLibrary(path.c_str());
        if( HInst>0 )
        {
                readEeprom = (bwtekReadEEPROMUSB)GetProcAddress(HInst,"bwtekReadEEPROMUSB");
                setupChannel=(bwtekSetupChannel)GetProcAddress(HInst,"bwtekSetupChannel");
                bwTest=(bwtekTestUSB)GetProcAddress(HInst,"bwtekTestUSB");
                bwSetTime=(bwtekSetTimeUSB) GetProcAddress(HInst,"bwtekSetTimeUSB");
                bwReadData=(bwtekDataReadUSB)GetProcAddress(HInst,"bwtekDataReadUSB");
                bwClose=(bwtekCloseUSB)GetProcAddress(HInst,"bwtekCloseUSB");
                bwReadAvgResult=(bwtekReadResultUSB)GetProcAddress(HInst,"bwtekReadResultUSB");
                bwPolyFit = (bwtekPolyFit)GetProcAddress(HInst,"bwtekPolyFit");
                bwPolyCalc = (bwtekPolyCalc)GetProcAddress(HInst,"bwtekPolyCalc");
                bwGetStatus = (bwtekGetExtStatus)GetProcAddress(HInst,"bwtekGetExtStatus");
                bwSetExtLaser = (bwtekSetExtLaser)GetProcAddress(HInst,"bwtekSetExtLaser");
                bwSetExtShutter = (bwtekSetExtShutter)GetProcAddress(HInst,"bwtekSetExtShutter");
                bwSetExtSync = (bwtekSetExtSync)GetProcAddress(HInst,"bwtekSetExtSync");
                bwGatedMode = (bwtekGatedMode)GetProcAddress(HInst,"bwtekGatedMode");
                bwSetPulse = (bwtekSetExtPulse)GetProcAddress(HInst,"bwtekSetExtPulse");
                closeShutter = (bwtekShutterOpen)GetProcAddress(HInst,"bwtekShutterOpen");
                openShutter = (bwtekShutterClose)GetProcAddress(HInst,"bwtekShutterClose");
                dllPath = path;
        }
        else
        {
                TFileOpen *o = new TFileOpen(owner);
                o->Caption = "Please Locate File bwtekusb.dll";
                o->Dialog->Filter ="Dynamically Linked Library (bwtekusb.dll)|BWTEKUSB.DLL";
                o->Dialog->Title = "Please Locate File bwtekusb.dll";
                open:
                if(o->Dialog->Execute())
                {
                        loadDLL(o->Dialog->FileName);
                }
                else
                {
                        goto open;
                }
                delete o;
        }
}


int bwTekUSB::setShutter(unsigned short state)
{
        int i=0;
        while(shutterState!=state && i<100)
        {
                if(state>0)
                {
                        shutterState=openShutter(iAddr);
                }
                else
                {
                        shutterState=closeShutter(iAddr);
                        shutterState^=1;
                }
        }
        if(i==100&&shutterState!=state)
        {
                shutterState=-1;
        }
        return shutterState;
}

int __fastcall bwTekUSB::avgSpec(spectrum *s, int points)
{
        if(connected)
        {
                for(int i =0;i<arrSize;i++)
                {
                        tmp[i]=0;
                }
                for(int i= 0;i<points;i++)
                {
                        bwReadData(0,rawData,iAddr);
                        for(int j=0;j<arrSize;j++)
                        {
                                tmp[j] = (double)(tmp[j]*i+rawData[j])/(double)(i+1);
                        }
                }
                unsigned short ti;
                for(int i =0;i<arrSize;i++)
                {
                        s->data[i] = (double)tmp[i]/maxY;
                }
        }
}

void __fastcall bwTekUSB::getDarkRef()
{
        if(connected)
        {
                setShutter(0);
                int then=GetTickCount()+shutterWait;
                while(GetTickCount()<then){};
                avgSpec(darkRef,refAvg);
                logReference(darkRef);
                setShutter(1);
        }
}

void __fastcall bwTekUSB::getSampleRef()
{
        if(connected)
        {
                setShutter(1);
                avgSpec(sampleRef,refAvg);
                logReference(sampleRef);
        }
}

spectrum *__fastcall spectrometer::getRelAbs()
{
        for(int i=0;i<arrSize;i++)
        {
                if(sample->data[i]==0)
                {
                        relAbs->data[i]=0;
                }
                else
                {
                        if(relTrans->data[i]==0)
                        {
                                relAbs->data[i]=65535;
                        }
                        else
                        {
                                try
                                {
                                        relAbs->data[i]=-log10(relTrans->data[i]);
                                }
                                catch(...)
                                {
                                        relAbs->data[i]=0;
                                }
                        }
                }
        }
        return relAbs;
}

int __fastcall bwTekUSB::setExtShutter(int state)
{
        if(connected)
        {
                return bwSetExtShutter(state,iAddr);
        }
}

void __fastcall bwTekUSB::setITime(int ITime)
{
        if(connected)
        {
                int result =bwSetTime(ITime,iAddr);
                if(result==ITime)
                {
                        iTime = ITime;
                }
        }
}

int __fastcall bwTekUSB::bwtest()
{
        int returnCode = bwTest(timingMode,arrSize,adcType,iAddr,NULL);
        if(returnCode>=0)
        {
                setShutter(shutterState);
        }
        return returnCode;
}

int __fastcall bwTekUSB::connect()
{
        int returnCode = bwtest();
//        int returnCode = -1;
        if(returnCode>=0)
        {
                connected = true;
        }
        else
        {
                connected = false;
        }
        return CONNECTED;
}

int __fastcall bwTekUSB::getSpectrum()
{
        if(connected)
        {
                try
                {
                        if(this->shutterState!=1)
                        {
                                setShutter(1);
                        }
                        int returnCode = bwReadData(0,rawData,iAddr);
                        if(returnCode>=0)
                        {
                                for(int i =0;i<arrSize;i++)
                                {
                                        sample->data[i]=(double)rawData[i]/65535;
                                }
                        }
                        return returnCode;
                }
                catch(...){}
        }
}


void __fastcall bwTekUSB::checkRef()
{
        if(GetTickCount()>=nextReferenceTime)
        {
                rereference();
                nextReferenceTime = GetTickCount()+referenceInterval;
        }
}
void __fastcall bwTekUSB::processSpectrum()
{
        (*measNo)++;
        getRelTrans();
        getRelAbs();
        if(onScan!=NULL)
        {
                onScan();
        }
        nextSampleTime=GetTickCount()+sampleInterval;
}

void __fastcall bwTekUSB::scanOnTimer(TObject *Sender)
{
        getState();
}

void __fastcall bwTekUSB::getRelTrans()
{
        for(int i =0;i<arrSize;i++)
        {
                if(darkRef->data[i]>=sample->data[i]||darkRef->data[i]>=sampleRef->data[i])
                {
                        relTrans->data[i]=0;
                }
                else
                {
                        float f =((sample->data[i]-darkRef->data[i]))/((sampleRef->data[i]-darkRef->data[i]));
                        if(f<1)
                        {
                                relTrans->data[i]=f;
                        }
                        else
                        {
                                relTrans->data[i]=1;
                        }
                }
        }
}

void __fastcall bwTekUSB::autoSetITime()
{
        if(connected)
        {
        if(shutterState<1)
        {
                setShutter(1);
        }
        double tmax = saturationThreshold;
        double imax =0;
        int tm=tmax*100;
        int im=imax*100;
        int di = 0;

        while(im!=tm)
        {
                di=0;
                setITime(iTime);
                avgSpec(sample,refAvg);
                for(int i =0;i<arrSize-1;i++)
                {
                        if(imax<sample->data[i])
                        {
                                imax = sample->data[i];
                        }
                }
                if(imax>tmax)
                {
                        di =-1;
                }
                if(imax<tmax)
                {
                        di = 1;
                }
                iTime+=di;
                im=imax*100;

        }
        setITime(iTime);
        sampleInterval = iTime+iTimeOffset;        
        initCriteria[iTimeOK]=true;
        bool b = initialized;
        }
}


void __fastcall bwTekUSB::initialize()
{
        if(connected)
        {
                getPixelConv();
                autoSetITime();
                getDarkRef();
                getSampleRef();
                getState();
                nextReferenceTime = GetTickCount()+referenceInterval;
                bool b = initialized;
        }
}

void __fastcall bwTekUSB::rereference()
{
        getDarkRef();
}


void __fastcall bwTekUSB::refreshIcon()
{
        int xo = (Width-128)/2;
        icon->Canvas->FillRect(TRect(xo,Height/4,Width-xo,Height-4));
        icon->Canvas->Pen->Color = clMedGray;
        icon->Canvas->MoveTo(xo,Height-3);
        icon->Canvas->LineTo(xo,Height/4);
        icon->Canvas->LineTo(Width-xo,Height/4);

        icon->Canvas->MoveTo(4,Height-3);
        icon->Canvas->Pen->Color = 0x00FF0000;
        int r = arrSize/128;
        int h= 3*Height/4;
        for(int i=0;i<128;i++)
        {
                int v1 = this->rawData[r*i]*h/maxY;
                icon->Canvas->LineTo(i+xo+1,Height-4-v1);
        }
        pathnode->draw();
}


bwTekS255_256::bwTekS255_256(process *Owner,String name, int chan,bool connectNow,input *tempSensor,int h,int w) : bwTekUSB(name,Owner,chan,256,14,3,connectNow,tempSensor,h,w)
{
        params->editable = false;
        firstLight = true;
        firstDark = true;
        fileParam *fp = new fileParam(&storedSampleRef,"Sample Reference File");
        fp->setFunc = &loadSample;
        params->add(fp);

        internalPerms->onLogin->addTask(new setVar<bool>("Set Sample Reference File Parameter Editable",&params->editable,new bool(true)));
        internalPerms->onLogout->addTask(new setVar<bool>("Set Sample Reference File Parameter Editable",&params->editable,new bool(false)));


        fp = new fileParam(&storedDarkRef,"Dark Reference File");
        fp->setFunc = &loadDark;
        params->add(fp);

        fp = new fileParam(&storedLightRef,"Light Reference File");
        fp->setFunc = &loadLight;
        params->add(fp);

        HINSTANCE HInst=LoadLibrary(dllPath.c_str());
        if( HInst>0 )
        {
                bwtekSetShutters=(bwtekShutterControl)GetProcAddress(HInst,"bwtekShutterControl");
        }
        refShutterState = 0;
        ref0=new spectrum(this,arrSize,"Original Reference Channel Spectrum");
        dark0=new spectrum(this,arrSize,"Original Dark Reference Spectrum");
        lastRef = new spectrum(this,arrSize,"Current Reference Channel Spectrum");

        spectroscopyTool = new spectroscopySuite(owner->owner,name+" Tools",this);
        MainForm->addProcess(spectroscopyTool);

        internalPerms->onLogin->addTask(new twoArgPtrTask<void,TTabSheet *,bool>("Set Spectroscopy Tab Visible",&setVisible,spectroscopyTool->tab,true));
        internalPerms->onLogout->addTask(new twoArgPtrTask<void,TTabSheet *,bool>("Set Spectroscopy Tab Visible",&setVisible,spectroscopyTool->tab,false));
        internalPerms->onLogin->addTask(new twoArgPtrTask<void,TMenuItem *,bool>("Set Spectroscopy Tab Visible",&setVisible,spectroscopyTool->procMainMenu,true));
        internalPerms->onLogout->addTask(new twoArgPtrTask<void,TMenuItem *,bool>("Set Spectroscopy Tab Visible",&setVisible,spectroscopyTool->procMainMenu,false));
        internalPerms->onLogout->addTask(new twoArgPtrTask<void,TTimer *,bool>("Set Spectroscopy Tab Visible",&setEnabled,((TspectroscopyTool *)spectroscopyTool->PnID)->timer,false));

        spectroscopyTool->tab->TabVisible = false;
        spectroscopyTool->procMainMenu->Visible = false;

        initSpectrometer = new initBwTekS255_256(this);
}



int bwTekS255_256::setShutter(unsigned short state)
{
        if(connected)
        {
                try
                {
                        if(bwtekSetShutters(state,-1,iAddr))
                        {
                                shutterState = state;
                        }
                }
                catch(...)
                {
                        
                }
        }
}

void __fastcall bwTekS255_256::setRefShutter(unsigned short state)
{
        if(connected)
        {
                try
                {
                        if(bwtekSetShutters(-1,state,iAddr))
                        {
                                refShutterState =state;
                        }
                }
                catch(...)
                {
                }
        }
}

void __fastcall bwTekS255_256::autoSetITime()
{
        if(connected)
        {
                setShutter(0);
                setRefShutter(1);
                double tmax = saturationThreshold;
                double imax =0;
                int di = 0;
                double lastErr =tmax+1;
                double err=tmax;
                while(lastErr>=err&&imax!=tmax)
                {
                        imax=0;
                        di=0;
                        lastErr = err;
                        setITime(iTime);
                        avgSpec(sample,refAvg);
                        for(int i =0;i<arrSize-1;i++)
                        {
                                if(imax<sample->data[i])
                                {
                                        imax = sample->data[i];
                                }
                        }
                        if(imax>tmax)
                        {
                                di =-1;
                        }
                        if(imax<tmax)
                        {
                                di = 1;
                        }
                        err =fabs(tmax-imax);
                        if(err>lastErr)
                        {
                                di-=2*di;
                        }
                        if(imax==tmax)
                        {
                                di=0;
                        }
                        iTime+=di;
                }
                setITime(iTime);
                sampleInterval = iTime+iTimeOffset;
                initCriteria[iTimeOK]=true;
                bool b = initialized;        
        }
}

void __fastcall bwTekS255_256::getLightRef()
{
        if(connected)
        {
                setShutter(0);
                setRefShutter(1);
                int then=GetTickCount()+shutterWait;
                while(GetTickCount()<then){}
                avgSpec(lastRef,refAvg);
                logReference(lastRef);
                setRefShutter(0);
                setShutter(1);
                bool b = initialized;
                if(firstLight)
                {
                        firstLight = false;
                        for(int i=0;i<arrSize;i++)
                        {
                                ref0->data[i]=lastRef->data[i];
                        }
                }
        }
}
void __fastcall bwTekS255_256::getDarkRef()
{
        if(connected)
        {
                setShutter(0);
                setRefShutter(0);
                int then=GetTickCount()+shutterWait;
                while(GetTickCount()<then){}
                avgSpec(darkRef,refAvg);
                logReference(darkRef);
                initCriteria[darkOK]=true;
                setRefShutter(0);
                setShutter(1);
                bool b = initialized;
                if(firstDark)
                {
                        firstDark = false;
                        for(int i=0;i<arrSize;i++)
                        {
                                dark0->data[i]=darkRef->data[i];
                        }
                }
        }
}

void __fastcall bwTekS255_256::getSampleRef()
{
        if(connected)
        {
                setShutter(1);
                setRefShutter(0);
                int then=GetTickCount()+shutterWait;
                while(GetTickCount()<then){};
                avgSpec(sampleRef,refAvg);
                logReference(sampleRef);
                setRefShutter(0);
                setShutter(1);
                initCriteria[sampleOK]=true;
                bool b = initialized;

        }
}

void __fastcall bwTekS255_256::initialize()
{
        if(connected)
        {
                initSpectrometer->Resume();
        }
}

void __fastcall bwTekS255_256::getRelTrans()
{
        for(int i =0;i<arrSize;i++)
        {
                if(darkRef->data[i]>=sample->data[i]||darkRef->data[i]>=sampleRef->data[i]||darkRef->data[i]>=lastRef->data[i]||dark0->data[i]>=ref0->data[i])
                {
                        relTrans->data[i]=0;
                }
                else
                {
                        float f =((sample->data[i]-darkRef->data[i]))/((sampleRef->data[i]-dark0->data[i]))*((lastRef->data[i]-darkRef->data[i])/(ref0->data[i]-dark0->data[i]));
                        relTrans->data[i]=f;
                }
        }
}

void __fastcall bwTekS255_256::loadPopup()
{
        bwTekUSB::loadPopup();
}

void __fastcall bwTekS255_256::rereference()
{
        if(connected)
        {
                getDarkRef();
                getLightRef();
        }
}

int __fastcall bwTekS255_256::bwtest()
{
        try
        {
                int returnCode = bwTest(timingMode,arrSize,adcType,iAddr,NULL);
                if(returnCode>=0)
                {
                        setShutter(shutterState);
                        setRefShutter(refShutterState);
                }
                return returnCode;
        }
        catch(...)
        {
        }
}

initBwTekS255_256::initBwTekS255_256(bwTekS255_256 *Owner) : initSpectrometer(Owner)
{
        owner = Owner;
}

void __fastcall initBwTekS255_256::Execute()
{
        while(!Terminated)
        {

                owner->progForm->progBar->POS=0;
                owner->progForm->progBar->min=0;
                owner->progForm->progBar->max = 5;

                owner->progForm->labelCaption="Reading EEPROM";
                Synchronize(owner->progForm->updateLabelCaption);
                Synchronize(owner->getPixelConv);
                owner->progForm->progBar->POS++;
                Synchronize(owner->progForm->progBar->Paint);
                owner->progForm->labelCaption="Detecting integration time...";
                Synchronize(owner->progForm->updateLabelCaption);

                Synchronize(owner->autoSetITime);
                owner->progForm->progBar->POS++;

                owner->progForm->labelCaption="Getting dark reference...";
                Synchronize(owner->progForm->progBar->Paint);
                Synchronize(owner->progForm->updateLabelCaption);
                Synchronize(owner->getDarkRef);


                owner->progForm->labelCaption="Getting light reference...";
                Synchronize(owner->progForm->updateLabelCaption);
                owner->progForm->progBar->POS++;
                Synchronize(owner->progForm->progBar->Paint);

                Synchronize(owner->getLightRef);

                owner->progForm->labelCaption="Getting water background...";
                Synchronize(owner->progForm->updateLabelCaption);

                owner->progForm->progBar->POS++;
                Synchronize(owner->progForm->progBar->Paint);

                Synchronize(owner->getSampleRef);
                owner->progForm->progBar->POS++;
                Synchronize(owner->progForm->progBar->Paint);

                owner->progForm->labelCaption="Finishing initialization";
                Synchronize(owner->progForm->updateLabelCaption);

                Synchronize(owner->getState);
                for(int i=0;i<owner->arrSize;i++)
                {
                        owner->ref0->data[i]=owner->lastRef->data[i];
                        owner->dark0->data[i]=owner->darkRef->data[i];
                }
                owner->nextReferenceTime = GetTickCount()+owner->referenceInterval;

                owner->progForm->close();
                Suspend();
        }
}

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