#include <utility.h>
#include <ansi_c.h>
#include <cvirte.h>		
#include <userint.h>
#include <formatio.h>
#include <string.h>
#include <future.h>

#include "YLCStuff.h"
#include "FXAnalyse.h" 
#include "Plot.h"
#include "Allan.h"
#include "DDSBes.h"
#include "DDS4xAD9912.h"
#include "DDS_Fox.h" 
#include "muParserDLL.h"

#define FXLINELENGTH 123
#define LOGFILEPATH "C:\\Femto\\Software\\FXQE80"

#define DEDRIFT_DDS_FREQUENCY 70000000

//==============================================================================
// Static global variables

/* multi threading stuff */  // to be used when I will replace timer by thread for counder readout
//static int ghPool;
//static volatile int vgExiting = 0;

/* panel handling stuff */
static PanelHandle MainPanel ;
static PanelHandle CalcN1Panel ;
static PanelHandle CalcN2Panel ;
static PanelHandle CalcN3Panel ;
static PanelHandle EstimateN3Panel;

//==============================================================================
// Static functions

//static int CVICALLBACK ReadCounterThread(void *functionData);	   // to be used when I will replace timer by thread for counder readout  

//==============================================================================
// Global variables

char LogFileName[MAX_PATHNAME_LEN];
char ExtraMathFileName[MAX_PATHNAME_LEN];

double utc = 0;
double mjd = 0;

double Ch1, Ch2, Ch3, Ch4 ; // le type long double est identique au double (8 octets) sous CVI !!! 
double Math1, Math2, Math3, Math4, Math5; // N est en fait un entier en pratique...
double N1, N2, N3, N4 ;

Plot_Data PlotCh1, PlotCh2, PlotCh3, PlotCh4, PlotMath1, PlotMath2, PlotMath3, PlotMath4, PlotMath5 ;  
Allan_Data AllanCh1, AllanCh2, AllanCh3, AllanCh4, AllanMath1, AllanMath2, AllanMath3, AllanMath4, AllanMath5 ; 

// 1xAD9956 DDS box
DDSParameter DDS1xAD9956;
// 4xAD9912 DDS box
DDS4xAD9912_Data DDS4xAD9912;

muParserHandle_t MathParser1, MathParser2, MathParser3, MathParser4, MathParser5 ;

volatile bool Acquiring = FALSE ;

long OldLogFilePtr = 0 ;
double Ndiv=8.0;

///////////////////
//Pour les calculs des N

enum {
	N_MEASUREMENT_STEP_0,
	N_MEASUREMENT_STEP_1,
	N_MEASUREMENT_STEP_2,
	N_MEASUREMENT_STEP_3,
	N_MEASUREMENT_STEP_4,
};

enum {
	N_MEASUREMENT_NONE,
	N_MEASUREMENT_INIT,
	N_MEASUREMENT_SLOPE,
	N_MEASUREMENT_ADJUST_FREQ_PLUS,
	N_MEASUREMENT_FREP_PLUS,
	N_MEASUREMENT_ADJUST_FREQ_MINUS,
	N_MEASUREMENT_FREP_MINUS,
};

int settling = 0;

volatile bool Measuring_1 = FALSE,Step1_1 = FALSE,Step2_1 = FALSE,Step3_1 = FALSE, Measuring_2 = FALSE,Step1_2 = FALSE,Step2_2 = FALSE,Step3_2 = FALSE;
volatile bool Measuring_3 = FALSE,Step1_3 = FALSE,Step2_3 = FALSE,Step3_3 = FALSE;
volatile bool TimetoBegin=FALSE;

double FrequDDS1=110000000.0, FrequDDS2=15300000.0, FrequDDS3=150400000.0, FrequDDS4=110000000.0;
double Slope_1=0.0,Slope_2=0.0,Slope_3=0.0,Beatslope_2=0.0; 
double SlopeTime1=40.0, SlopeTime2=40.0; SlopeTime3=40.0;

double N_1=0.0, N_2=0.0, N_3=0.0;
double DeltaT_1=20.0, DeltakHz_1=500.0, t1_1=0.0, t2_1=0.0, t3_1=0.0, Frepplus_1=0.0, Frepminus_1=0.0;
double DeltaT_2=20.0, DeltakHz_2=500.0, t1_2=0.0, t2_2=0.0, t3_2=0.0, Frepplus_2=0.0, Frepminus_2=0.0;
double DeltaT_3=20.0, DeltakHz_3=500.0, t1_3=0.0, t2_3=0.0, t3_3=0.0, Frepplus_3=0.0, Frepminus_3=0.0;

int n_1=0, n_2=0, n_3=0;
double Frequ_slope_1=0.0,Moy_slope_1=0.0,Slope_slope_1=0.0,Frequ_slope_2=0.0,Moy_slope_2=0.0,Slope_slope_2=0.0,Frequ_slope_3=0.0,Moy_slope_3=0.0,Slope_slope_3=0.0;  
int N_slope_1=0,N_slope_2=0,N_slope_3=0;
double  Beat_slope_2=0.0 ,Moy_Beatslope_2=0.0,Slope_Beatslope_2=0.0; 
double Ch4_slope=0.0,Moy_Ch4slope_1=0.0,Slope_Ch4slope_1=0.0,Ch4Slope=0.0;

double FrequencyDDSBes=0.0,FrequencyDDSBesInit=0.0;
double FrequencyDDS3=0.0,FrequencyDDS3Init=0.0;
double DDSBesChanged1=FALSE,DDSBesChanged2=FALSE;

double DeltaDDS3=0.0,Delta10K_Plus=0.0,Delta10K_Minus=0.0;
double Nu1=0.0, Nu2= 200000-147000+282143746.557455e6;  

double Step1=800000.0,Step2=800000.0;

double Ch4Plus=0.0,Ch4Minus=0.0;

int nDDSChange=0, nstabilization=0;

//////
//Pour les signes de fb

double Frequencystep1=10000.0, tbegin1=0.0, Frepbefore1=0.0, Frequency1=0.0;
double Frequencystep2=10.0, tbegin2=0.0, Frepbefore2=0.0, Ch2before=0.0, Frequency2=0.0;
double Frequencystep3=100000.0, tbegin3=0.0, Frepbefore3=0.0, Frequency3=0.0;

volatile bool Getsign1=FALSE,Getsign2=FALSE,Getsign3=FALSE;
double Signe1=1.0,Signe2=1.0,Signe3=0.0;





/////////////////////////
// Slope Cancelling Variables

volatile bool SlopeMeasuring=FALSE;
volatile bool AutoStopSlopeCancellingIfDelock=TRUE;
volatile bool OnSlopeCancelling =FALSE;  
double SlopeMath2=0.0;
double TimetoSlope=60.0;
double SlopeMeasuringTimeBegin=0.0;

double Slope_Math2slope=0.0,Math2_slope=0.0,Moy_Math2slope=0.0;
int N_Math2slope=0.0,nstabilisationSlopeMeasuring=0;

double LimitToDelock=5.0;
double SlopeCorrection=0.0;
double limitotakoff=70.0;

int ratio=10; //Recentre la frequence tous les ratios

volatile bool FrequCorrec=FALSE,KeepFrequ=TRUE, KeepSlope=TRUE ;
int Nratio=-1;
double MoyMath2 = 0.0, CenterFrequencyCh2 = 0.0;
volatile bool CenterFrequencyCh2ToDetermine=FALSE;



enum {
	SLOPE_REFERENCE_MICROWAVE = 0,
	SLOPE_REFERENCE_HG_CAVITY = 1,
};

int slopeReference = SLOPE_REFERENCE_MICROWAVE;

//////////////////////////////////
//  Keep Centered Variables

volatile bool KeepCentering=FALSE; 
double Timetorecenter275K=3600.0*10; 
double Timetorecenter10K=3600.0*3;  
double CenteringTimeBegin275K=0.0;
double CenteringTimeBegin10K=0.0; 


int DDSFox_Set(DDSParameter *DDS, double Frequency, double Sweeprate)
{
	DDS->Frequency = Frequency; 
	DDS->SweepRate = Sweeprate;
	
	DDSFox_SetFreq(DDS);
	DDSFox_SetDf(DDS);
	
	return 0;
}


void initMathParser(muParserHandle_t * pMathParser) 
{
	(* pMathParser) = mupCreate() ;		// Math parser number 1, normally used for frep
	mupDefineOprtChars( * pMathParser, "abcdefghijklmnopqrstuvwxyzµ"
                              		   "ABCDEFGHIJKLMNOPQRSTUVWXYZ"
                            		   "+-*^/?<>=#!$%&|~'_");
	mupDefineVar( * pMathParser, "UTC", &utc);
	mupDefineVar( * pMathParser, "MJD", &mjd);
	mupDefineVar( * pMathParser, "Ch1", &Ch1);
	mupDefineVar( * pMathParser, "Ch2", &Ch2);
	mupDefineVar( * pMathParser, "Ch3", &Ch3);
	mupDefineVar( * pMathParser, "Ch4", &Ch4);
	mupDefineVar( * pMathParser, "DDS1", &(DDS4xAD9912.Frequency1));   
	mupDefineVar( * pMathParser, "DDS2", &(DDS4xAD9912.Frequency2));
	mupDefineVar( * pMathParser, "DDS3", &(DDS4xAD9912.Frequency3));  
	mupDefineVar( * pMathParser, "DDS4", &(DDS4xAD9912.Frequency4));  
	mupDefineVar( * pMathParser, "N1", &N1);
	mupDefineVar( * pMathParser, "N2", &N2);
	mupDefineVar( * pMathParser, "N3", &N3);
	mupDefineVar( * pMathParser, "Nu1", &Nu1);
	mupDefineVar( * pMathParser, "Nu2", &Nu2);
	mupDefineVar( * pMathParser, "DeltaDDS3", &DeltaDDS3);
	mupDefineVar( * pMathParser, "Signe1", &Signe1);
	mupDefineVar( * pMathParser, "Signe2", &Signe2); 
	mupDefineVar( * pMathParser, "Ndiv", &Ndiv); 
	mupDefinePostfixOprt( * pMathParser, "P", &Peta , 1);
	mupDefinePostfixOprt( * pMathParser, "T", &Tera , 1);
	mupDefinePostfixOprt( * pMathParser, "G", &Giga , 1);
	mupDefinePostfixOprt( * pMathParser, "M", &Mega , 1);
	mupDefinePostfixOprt( * pMathParser, "k", &kilo , 1);
	mupDefinePostfixOprt( * pMathParser, "m", &milli, 1);
	mupDefinePostfixOprt( * pMathParser, "u", &micro, 1);
	mupDefinePostfixOprt( * pMathParser, "µ", &micro, 1);
	mupDefinePostfixOprt( * pMathParser, "n_1", &nano,  1);
	mupDefinePostfixOprt( * pMathParser, "p", &pico,  1);
	mupDefinePostfixOprt( * pMathParser, "f", &femto, 1);
}

int main (int argc, char *argv[])
{
	//int ThreadId ;   // to be used when I will replace timer by thread for counder readout
	
	double initfreqDDS ;
	char * MathString1 = "                                                             " ;
	char * MathString2 = "                                                             " ;
	char * MathString3 = "                                                             " ;
	char * MathString4 = "                                                             " ;  
	char * MathString5 = "                                                             " ;
	
	if (InitCVIRTE (0, argv, 0) == 0)
		return -1;	/* out of memory */
	if ((MainPanel = LoadPanel (0, "FXAnalyse.uir", PANEL)) < 0)
		return -1;
	if ((CalcN1Panel = LoadPanel (0, "FXAnalyse.uir", CALCN1)) < 0)
		return -1;
	if ((CalcN2Panel = LoadPanel (0, "FXAnalyse.uir", CALCN2)) < 0)
		return -1;
	if ((CalcN3Panel = LoadPanel (0, "FXAnalyse.uir", CALCN3)) < 0)
		return -1;
	if ((EstimateN3Panel = LoadPanel (MainPanel, "FXAnalyse.uir", ESTIMATEN3)) < 0)
		return -1;
	
	DisplayPanel (MainPanel);
	
	// initialize 4x AD9959 DDS box
	DDS4xAD9912_Reset(&DDS4xAD9912) ;
	GetCtrlVal(MainPanel, PANEL_DDS1, &initfreqDDS) ;
	DDS4xAD9912_SetFrequency(&DDS4xAD9912, 1, initfreqDDS) ;
	GetCtrlVal(MainPanel, PANEL_DDS2, &initfreqDDS) ;
	DDS4xAD9912_SetFrequency(&DDS4xAD9912, 2, initfreqDDS) ;
	GetCtrlVal(MainPanel, PANEL_DDS3, &initfreqDDS) ;
	DDS4xAD9912_SetFrequency(&DDS4xAD9912, 3, initfreqDDS) ;
	GetCtrlVal(MainPanel, PANEL_DDS4, &initfreqDDS) ;
	DDS4xAD9912_SetFrequency(&DDS4xAD9912, 4, initfreqDDS) ;
	
	
	// initialyze 1x AD9956 DDS box
	strcpy(DDS1xAD9956.ip, "145.238.205.58"); 
	DDS1xAD9956.Port = 6665;
	DDS1xAD9956.Profil = 7; // configuration profile: single freq=0, sweep=7
	DDS1xAD9956.Clock = 200000000;
	DDS1xAD9956.Frequency = DEDRIFT_DDS_FREQUENCY; 
	DDS1xAD9956.Delta_T = 0.01;
	DDS1xAD9956.SweepRate = 0.0;
	DDS_Initialize(&DDS1xAD9956);
	
	// Initialize the 5 (separate) math parsers for frep, Nu_Hg, Nu_Sr, N3 and ExtraMath

	Ch1=0.0 ; Ch2=0.0 ; Ch3=0.0 ; Ch4=0.0 ;
	GetCtrlVal(MainPanel, PANEL_N1CHOICE, &N1) ;
	GetCtrlVal(MainPanel, PANEL_N2CHOICE, &N2) ;

	GetCtrlVal(MainPanel, PANEL_N3CHOICE, &N3) ;
	
	initMathParser(&MathParser1);
	GetCtrlVal(MainPanel, PANEL_MATHSTRING1, MathString1);
	mupSetExpr(MathParser1, MathString1);
	
	initMathParser(&MathParser2) ;
	mupDefineVar( MathParser2, "Math1", &Math1);
	mupDefineVar( MathParser2, "DDS", &(DDS4xAD9912.Frequency1));
	GetCtrlVal(MainPanel, PANEL_MATHSTRING2, MathString2) ;
	mupSetExpr(MathParser2, MathString2);
	
	initMathParser(&MathParser3) ;
	mupDefineVar( MathParser3, "Math1", &Math1);
	mupDefineVar( MathParser3, "Math2", &Math2);
	mupDefineVar( MathParser3, "DDS", &(DDS4xAD9912.Frequency2));
	GetCtrlVal(MainPanel, PANEL_MATHSTRING3, MathString3) ;
	mupSetExpr(MathParser3, MathString3);
	
	initMathParser(&MathParser4) ;
	mupDefineVar( MathParser4, "Math1", &Math1);
	mupDefineVar( MathParser4, "Math2", &Math2);
	mupDefineVar( MathParser4, "Math3", &Math3);
	GetCtrlVal(MainPanel, PANEL_MATHSTRING4, MathString4) ;
	mupSetExpr(MathParser4, MathString4);
	
	initMathParser(&MathParser5) ;
	mupDefineVar( MathParser5, "Math1", &Math1);
	mupDefineVar( MathParser5, "Math2", &Math2);
	mupDefineVar( MathParser5, "Math3", &Math3);
	mupDefineVar( MathParser5, "Math4", &Math4);
	GetCtrlVal(MainPanel, PANEL_MATHSTRING5, MathString5) ;
	mupSetExpr(MathParser5, MathString5);
	
	
	RunUserInterface ();
	
	// Stops everything at the end
	
	DiscardPanel (MainPanel);
	return 0;
}


void OnCloseViewPanel(int panel){
	
	if (panel==PlotCh1.PlotPanel) 		{ SetCtrlVal(MainPanel,PANEL_CHECKBOX_FREQ1PLOT, FALSE) ; } ; 
	if (panel==PlotCh2.PlotPanel) 		{ SetCtrlVal(MainPanel,PANEL_CHECKBOX_FREQ2PLOT, FALSE) ; } ; 
	if (panel==PlotCh3.PlotPanel) 		{ SetCtrlVal(MainPanel,PANEL_CHECKBOX_FREQ3PLOT, FALSE) ; } ; 
	if (panel==PlotCh4.PlotPanel) 		{ SetCtrlVal(MainPanel,PANEL_CHECKBOX_FREQ4PLOT, FALSE) ; } ;
	if (panel==PlotMath1.PlotPanel) 	{ SetCtrlVal(MainPanel,PANEL_CHECKBOX_MATH1PLOT, FALSE) ; } ; 
	if (panel==PlotMath2.PlotPanel) 	{ SetCtrlVal(MainPanel,PANEL_CHECKBOX_MATH2PLOT, FALSE) ; } ;
	if (panel==PlotMath3.PlotPanel) 	{ SetCtrlVal(MainPanel,PANEL_CHECKBOX_MATH3PLOT, FALSE) ; } ; 
	if (panel==PlotMath4.PlotPanel) 	{ SetCtrlVal(MainPanel,PANEL_CHECKBOX_MATH4PLOT, FALSE) ; } ;
	if (panel==PlotMath5.PlotPanel) 	{ SetCtrlVal(MainPanel,PANEL_CHECKBOX_MATH5PLOT, FALSE) ; } ;
	
	if (panel==AllanCh1.AllanPanel) 	{ SetCtrlVal(MainPanel,PANEL_CHECKBOX_FREQ1ALLAN, FALSE) ; } ; 
	if (panel==AllanCh2.AllanPanel) 	{ SetCtrlVal(MainPanel,PANEL_CHECKBOX_FREQ2ALLAN, FALSE) ; } ; 
	if (panel==AllanCh3.AllanPanel) 	{ SetCtrlVal(MainPanel,PANEL_CHECKBOX_FREQ3ALLAN, FALSE) ; } ; 
	if (panel==AllanCh4.AllanPanel) 	{ SetCtrlVal(MainPanel,PANEL_CHECKBOX_FREQ4ALLAN, FALSE) ; } ; 
	if (panel==AllanMath1.AllanPanel) 	{ SetCtrlVal(MainPanel,PANEL_CHECKBOX_MATH1ALLAN, FALSE) ; } ; 
	if (panel==AllanMath2.AllanPanel) 	{ SetCtrlVal(MainPanel,PANEL_CHECKBOX_MATH2ALLAN, FALSE) ; } ; 
	if (panel==AllanMath3.AllanPanel) 	{ SetCtrlVal(MainPanel,PANEL_CHECKBOX_MATH3ALLAN, FALSE) ; } ; 
	if (panel==AllanMath4.AllanPanel) 	{ SetCtrlVal(MainPanel,PANEL_CHECKBOX_MATH4ALLAN, FALSE) ; } ; 
	if (panel==AllanMath5.AllanPanel) 	{ SetCtrlVal(MainPanel,PANEL_CHECKBOX_MATH5ALLAN, FALSE) ; } ;  
	
	return ;
}





int CVICALLBACK QuitCallback (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	switch (event)
		{
		case EVENT_COMMIT:
			QuitUserInterface (0);
			mupRelease(MathParser1);
			mupRelease(MathParser2); 
			mupRelease(MathParser3);
			mupRelease(MathParser4);
			mupRelease(MathParser5);
			break;
		case EVENT_RIGHT_CLICK:

			break;
		}
	return 0;
}

int  CVICALLBACK CB_OnEventMain(int panel, int event, void *callbackData, int eventData1, int eventData2)
{
	int ActiveControl ;
	int StepIndex ;
	double Step ;
	
	switch (event) {
		case EVENT_CLOSE:
			QuitUserInterface (0);
			mupRelease(MathParser1);
			mupRelease(MathParser2); 
			mupRelease(MathParser3);
			break;
		case EVENT_KEYPRESS:
			switch (eventData1)	// ie the code of the key which has been stroke 
				{
				case VAL_RIGHT_ARROW_VKEY :
					ActiveControl = GetActiveCtrl(panel);
					if (ActiveControl==PANEL_DDS1 || ActiveControl==PANEL_DDS1STEP) {
						GetCtrlIndex(MainPanel, PANEL_DDS1STEP, &StepIndex);
						if (StepIndex<14){
							SetCtrlIndex(MainPanel, PANEL_DDS1STEP, ++StepIndex) ;
							GetCtrlVal(MainPanel, PANEL_DDS1STEP, &Step);
							SetCtrlAttribute(MainPanel, PANEL_DDS1, ATTR_INCR_VALUE, Step) ;
						};
					};
					if (ActiveControl==PANEL_DDS2 || ActiveControl==PANEL_DDS2STEP) {
						GetCtrlIndex(MainPanel, PANEL_DDS2STEP, &StepIndex);
						if (StepIndex<14){
							SetCtrlIndex(MainPanel, PANEL_DDS2STEP, ++StepIndex) ;
							GetCtrlVal(MainPanel, PANEL_DDS2STEP, &Step);
							SetCtrlAttribute(MainPanel, PANEL_DDS2, ATTR_INCR_VALUE, Step) ;
						};
					};
					if (ActiveControl==PANEL_DDS3|| ActiveControl==PANEL_DDS3STEP) {
						GetCtrlIndex(MainPanel, PANEL_DDS3STEP, &StepIndex);
						if (StepIndex<14){
							SetCtrlIndex(MainPanel, PANEL_DDS3STEP, ++StepIndex) ;
							GetCtrlVal(MainPanel, PANEL_DDS3STEP, &Step);
							SetCtrlAttribute(MainPanel, PANEL_DDS3, ATTR_INCR_VALUE, Step) ;
						};
					};
					if (ActiveControl==PANEL_DDS4|| ActiveControl==PANEL_DDS4STEP) {
						GetCtrlIndex(MainPanel, PANEL_DDS4STEP, &StepIndex);
						if (StepIndex<14){
							SetCtrlIndex(MainPanel, PANEL_DDS4STEP, ++StepIndex) ;
							GetCtrlVal(MainPanel, PANEL_DDS4STEP, &Step);
							SetCtrlAttribute(MainPanel, PANEL_DDS4, ATTR_INCR_VALUE, Step) ;
						};
					};
					break;
				case VAL_LEFT_ARROW_VKEY :
					ActiveControl = GetActiveCtrl(panel);
					if (ActiveControl==PANEL_DDS1 || ActiveControl==PANEL_DDS1STEP) {
						GetCtrlIndex(MainPanel, PANEL_DDS1STEP, &StepIndex);
						if (StepIndex>0){
							SetCtrlIndex(MainPanel, PANEL_DDS1STEP, --StepIndex) ;
							GetCtrlVal(MainPanel, PANEL_DDS1STEP, &Step);
							SetCtrlAttribute(MainPanel, PANEL_DDS1, ATTR_INCR_VALUE, Step) ;
						};
					};
					if (ActiveControl==PANEL_DDS2 || ActiveControl==PANEL_DDS2STEP) {
						GetCtrlIndex(MainPanel, PANEL_DDS2STEP, &StepIndex);
						if (StepIndex>0){
							SetCtrlIndex(MainPanel, PANEL_DDS2STEP, --StepIndex) ;
							GetCtrlVal(MainPanel, PANEL_DDS2STEP, &Step);
							SetCtrlAttribute(MainPanel, PANEL_DDS2, ATTR_INCR_VALUE, Step) ;
						};
					};
					if (ActiveControl==PANEL_DDS3 || ActiveControl==PANEL_DDS3STEP) {
						GetCtrlIndex(MainPanel, PANEL_DDS3STEP, &StepIndex);
						if (StepIndex>0){
							SetCtrlIndex(MainPanel, PANEL_DDS3STEP, --StepIndex) ;
							GetCtrlVal(MainPanel, PANEL_DDS3STEP, &Step);
							SetCtrlAttribute(MainPanel, PANEL_DDS3, ATTR_INCR_VALUE, Step) ;
						};
					};
					if (ActiveControl==PANEL_DDS4 || ActiveControl==PANEL_DDS4STEP) {
						GetCtrlIndex(MainPanel, PANEL_DDS4STEP, &StepIndex);
						if (StepIndex>0){
							SetCtrlIndex(MainPanel, PANEL_DDS4STEP, --StepIndex) ;
							GetCtrlVal(MainPanel, PANEL_DDS4STEP, &Step);
							SetCtrlAttribute(MainPanel, PANEL_DDS4, ATTR_INCR_VALUE, Step) ;
						};
					};
					break;
				case VAL_F2_VKEY :
					SetActiveCtrl(MainPanel, PANEL_DDS1);
					break;
				case VAL_F3_VKEY :
					SetActiveCtrl(MainPanel, PANEL_DDS2);
					break;
				case VAL_F4_VKEY :
					SetActiveCtrl(MainPanel, PANEL_DDS3);
					break;
				case VAL_F5_VKEY :
					SetActiveCtrl(MainPanel, PANEL_DDS4);
					break;
				};
			break;
			
	}
	return 0;
}


void CurrentFileName(char *fname)
{
	char day[3], month[3], year[3];
	char *date = DateStr();
	Scan(date, "%s>%s[w2]-%s[w2]-20%s[w2]", month, day, year);
	Fmt(fname, "%s<%s\\%s%s%s_Frequ.txt", LOGFILEPATH, year, month, day);
}


int CVICALLBACK CB_OnStart (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	switch (event)
	{
		case EVENT_COMMIT:
			if (Acquiring) {
				PlotCh1.IndexPoint = 0;
				PlotCh2.IndexPoint = 0;
				PlotCh3.IndexPoint = 0;
				PlotCh4.IndexPoint = 0;
				PlotMath1.IndexPoint = 0;
				PlotMath2.IndexPoint = 0; 
				PlotMath3.IndexPoint = 0; 
				PlotMath4.IndexPoint = 0; 
				PlotMath5.IndexPoint = 0;
				Allan_Reset(&AllanCh1);
				Allan_Reset(&AllanCh2);
				Allan_Reset(&AllanCh3);
				Allan_Reset(&AllanCh4);
				Allan_Reset(&AllanMath1);
				Allan_Reset(&AllanMath2);
				Allan_Reset(&AllanMath3);
				Allan_Reset(&AllanMath4);
				Allan_Reset(&AllanMath5);
			}
			Acquiring = TRUE;
			SetCtrlAttribute(MainPanel, PANEL_STARTBUTTON, ATTR_LABEL_TEXT, "__RESET"); 
			
			CurrentFileName(LogFileName);
			GetFileInfo(LogFileName, &OldLogFilePtr);
			OldLogFilePtr -= OldLogFilePtr%FXLINELENGTH + FXLINELENGTH - 2;
			break;
		}
	return 0;
}

int CVICALLBACK CB_OnStop (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	switch (event)
		{
		case EVENT_COMMIT:
			Acquiring = FALSE ;
			SetCtrlAttribute(MainPanel, PANEL_STARTBUTTON, ATTR_LABEL_TEXT, "__START");
			break;
		case EVENT_RIGHT_CLICK:

			break;
		}
	return 0;
}

int CVICALLBACK CB_OnTimer (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	int LogFile, FileOpt, FileHg, FileExtraMath;
	long LogFileSize;
	char LineBuffer[FXLINELENGTH+10] = "\r\n_1 ";
	
	char TimeTag[] = "100103 000000.000";   // K+K time tag meaning here 2010 january the 3rd at 00:00:00.000
	char Date[] = "03/01/2010" ;
	char Time[] = "00:00:00.000" ;
	char Year[] = "2010";
	char ShortYear[] = "10";  // the last 2 digits of calendar year only
	char Month[] = "01";
	char Day[] = "03";
	char Hour[] = "00";
	char Min[] = "00" ;
	char Sec[] = "00.000";
	struct tm LocalTime ;
	time_t utcTime;
	char * ReportString = "dd.mm.yy\t00:00:00\t3481610838.000\t0000000000000.00000000";

	int BoxChecked = FALSE; 

	double FoxFrequ=0.0;
	
	double FrequencyToChange=0.0;
	double DeltaCh2=0.0;
	double DeltaCh4=0.0;  
	
	double ActualFrequ=0.0;
	double ActualSlope=0.0;
	
	switch (event)
	{
		case EVENT_TIMER_TICK:
			switch(Acquiring) 
				{
				case TRUE:
					
					GetFileInfo(LogFileName, &LogFileSize) ;
					
					if (LogFileSize > OldLogFilePtr+2*FXLINELENGTH-2) {  // if a complete newline has been written
						
						SuspendTimerCallbacks() ; 
						
						// Open Log file and get to the beginning of newly completed line
						LogFile = OpenFile(LogFileName, VAL_READ_ONLY, VAL_OPEN_AS_IS, VAL_ASCII) ;
						OldLogFilePtr += FXLINELENGTH;
						SetFilePtr(LogFile, OldLogFilePtr, 0) ;
					
						// return the last complete string from the log file and scan it for date and time information 
						
						// first, the time tag, and store it in various formats
						ReadFile(LogFile, TimeTag, 17);
						
						CopyBytes(Date,0,TimeTag,4,2);
						CopyBytes(Date,3,TimeTag,2,2);   
						CopyBytes(Date,8,TimeTag,0,2);
						CopyBytes(Time,0,TimeTag,7,2);
						CopyBytes(Time,3,TimeTag,9,2);
						CopyBytes(Time,6,TimeTag,11,6);
						SetCtrlVal(MainPanel, PANEL_DATE, Date) ;
						SetCtrlVal(MainPanel, PANEL_TIME, Time) ; 
						
						CopyBytes(Year,2,TimeTag,0,2); // first 2 bytes of year string remains "20"
						CopyBytes(ShortYear,0,TimeTag,0,2);
						CopyBytes(Month,0,TimeTag,2,2);
						CopyBytes(Day,0,TimeTag,4,2);
						CopyBytes(Hour,0,TimeTag,7,2);
						CopyBytes(Min,0,TimeTag,9,2);
						CopyBytes(Sec,0,TimeTag,11,6);
						Fmt(&LocalTime.tm_year, "%d<%s", Year);
						Fmt(&LocalTime.tm_mon,  "%d<%s", Month);
						Fmt(&LocalTime.tm_mday, "%d<%s", Day);
						Fmt(&LocalTime.tm_hour, "%d<%s", Hour);
						Fmt(&LocalTime.tm_min,  "%d<%s", Min);
						Fmt(&LocalTime.tm_sec,  "%d<%s", "00");    // special case to handle non integer number of UTC seconds
						LocalTime.tm_hour += 0;
						LocalTime.tm_min  -= 0;  
						LocalTime.tm_sec  -= 0;  
						LocalTime.tm_mon  -= 1 ;  	// january is month 0 for tm struct
						LocalTime.tm_year -= 1900 ; // year is number of years since 1900 for tm struct
						LocalTime.tm_isdst = -1;  // daylight saving flag MUST be set to -1 (unallocated is bugging and +1 is making 1 hour error in summer)
						utcTime = mktime (&LocalTime);
						utc = (double) utcTime + strtod(Sec,NULL) ;
						mjd=utc/86400.;		//nb de jours depuis l'origine d'UTC (01/01/1900 à 00h00 GMT)
						mjd+=15020;			//date MJD de la date origine d'UTC
						SetCtrlVal(MainPanel, PANEL_UTC, utc);
						SetCtrlVal(MainPanel, PANEL_MJD, mjd);
				
						// scan the line for counters's channels information
						
						ReadLine(LogFile, LineBuffer, FXLINELENGTH+9) ;
						CloseFile(LogFile);
						
						Scan(LineBuffer, "%f%f%f%f", &Ch1, &Ch2, &Ch3, &Ch4) ;
						Ch1 = 1000*Ch1;
						Ch2 = 1000*Ch2;
						Ch3 = 1000*Ch3;
						Ch4 = 1000*Ch4;
						
						SetCtrlVal(MainPanel, PANEL_FREQ1, Ch1) ;
						SetCtrlVal(MainPanel, PANEL_FREQ2, Ch2) ;
						SetCtrlVal(MainPanel, PANEL_FREQ3, Ch3) ;  
						SetCtrlVal(MainPanel, PANEL_FREQ4, Ch4) ;
						SetCtrlVal(MainPanel, PANEL_SIGN1, Signe1);
						SetCtrlVal(MainPanel, PANEL_SIGN2, Signe2); 
						SetCtrlVal(MainPanel, PANEL_SIGN3, Signe3); 
						SetCtrlVal(CalcN1Panel, CALCN1_N, N_1) ;
						SetCtrlVal(CalcN1Panel, CALCN1_SLOPE, Slope_1) ;
						SetCtrlVal(CalcN2Panel, CALCN2_N, N_2) ;
						SetCtrlVal(CalcN2Panel, CALCN2_SLOPE, Beatslope_2) ;
						SetCtrlVal(MainPanel, PANEL_SLOPETOCANCEL, SlopeMath2);
						SetCtrlVal(MainPanel, PANEL_CENTERFREQUENCY, CenterFrequencyCh2);  
						SetCtrlVal(MainPanel, PANEL_LED1, SlopeMeasuring); 	   
						SetCtrlVal(MainPanel, PANEL_LED2, OnSlopeCancelling);
						// Treat data
						
						Math1 = mupEval(MathParser1) ;
						SetCtrlVal(MainPanel,PANEL_MATH1, Math1) ;

						Math2 = mupEval(MathParser2) ;
						SetCtrlVal(MainPanel,PANEL_MATH2, Math2) ;
						
						Math3 = mupEval(MathParser3) ;
						SetCtrlVal(MainPanel,PANEL_MATH3, Math3) ;
						
						Math4 = mupEval(MathParser4) ;
						SetCtrlVal(MainPanel,PANEL_MATH4, Math4) ;
						
						Math5 = mupEval(MathParser5) ;
						SetCtrlVal(MainPanel,PANEL_MATH5, Math5) ;
						
						// Plot Data and calculus if required
						
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_FREQ1PLOT, &BoxChecked) ;
						if (BoxChecked) {
							Plot_AddFrequency(&PlotCh1, Ch1) ;
							}
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_FREQ1ALLAN, &BoxChecked) ;
						if (BoxChecked) {
							Allan_AddFrequency(&AllanCh1, Ch1) ;
							}

						GetCtrlVal(MainPanel, PANEL_CHECKBOX_FREQ2PLOT, &BoxChecked) ;
						if (BoxChecked) {
							Plot_AddFrequency(&PlotCh2, Ch2) ;
							}
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_FREQ2ALLAN, &BoxChecked) ;
						if (BoxChecked) {
							Allan_AddFrequency(&AllanCh2, Ch2) ;
							}
							
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_FREQ3PLOT, &BoxChecked) ;
						if (BoxChecked) {
							Plot_AddFrequency(&PlotCh3, Ch3) ;
							}
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_FREQ3ALLAN, &BoxChecked) ;
						if (BoxChecked) {
							Allan_AddFrequency(&AllanCh3, Ch3) ;
							}
							
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_FREQ4PLOT, &BoxChecked) ;
						if (BoxChecked) {
							Plot_AddFrequency(&PlotCh4, Ch4) ;
							}
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_FREQ4ALLAN, &BoxChecked) ;
						if (BoxChecked) {
							Allan_AddFrequency(&AllanCh4, Ch4) ;
							}
						
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_MATH1PLOT, &BoxChecked) ;
						if (BoxChecked) {
							Plot_AddFrequency(&PlotMath1, Math1) ;
							}
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_MATH1ALLAN, &BoxChecked) ;
						if (BoxChecked) {
							Allan_AddFrequency(&AllanMath1, Math1) ;
							}
	
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_MATH2PLOT, &BoxChecked) ;
						if (BoxChecked) {
							Plot_AddFrequency(&PlotMath2, Math2) ;
							}
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_MATH2ALLAN, &BoxChecked) ;
						if (BoxChecked) {
							Allan_AddFrequency(&AllanMath2, Math2) ;
							}
						
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_MATH3PLOT, &BoxChecked) ;
						if (BoxChecked) {
							Plot_AddFrequency(&PlotMath3, Math3) ;
							}
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_MATH3ALLAN, &BoxChecked) ;
						if (BoxChecked) {
							Allan_AddFrequency(&AllanMath3, Math3) ;
							}
						
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_MATH4PLOT, &BoxChecked) ;
						if (BoxChecked) {
							Plot_AddFrequency(&PlotMath4, Math4) ;
							}
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_MATH4ALLAN, &BoxChecked) ;
						if (BoxChecked) {
							Allan_AddFrequency(&AllanMath4, Math4) ;
							}
						
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_MATH5PLOT, &BoxChecked) ;
						if (BoxChecked) {
							Plot_AddFrequency(&PlotMath5, Math5) ;
							}
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_MATH5ALLAN, &BoxChecked) ;
						if (BoxChecked) {
							Allan_AddFrequency(&AllanMath5, Math5) ;
							}
						
						// Calcul de N
						
						switch (Measuring_1) {

							case N_MEASUREMENT_NONE:
								// not measuring
								break;
							
							case N_MEASUREMENT_INIT:
								// initialization step
								
								// set DDS1 to nominal frequency
								SetCtrlVal(MainPanel, PANEL_DDS1, FrequDDS1);
								DDS4xAD9912_SetFrequency(&DDS4xAD9912, 1, FrequDDS1);
									
								GetCtrlVal(MainPanel, PANEL_DDS2, &FrequencyDDSBesInit);
								t2_1 = t3_1 = 0.0;
								t1_1 = utc;
								Frequ_slope_1 = Math1;
								Moy_slope_1 = Frequ_slope_1;
								Ch4_slope = Ch4;
								Moy_Ch4slope_1 = Ch4_slope;
								N_slope_1 = 1;
								
								// next step
								Measuring_1 += 1;
								break;
								
							case N_MEASUREMENT_SLOPE:
								// slope measurement
								
								N_slope_1 = N_slope_1 + 1; 
								Frequ_slope_1 = Math1;
								Ch4_slope = Ch4; 
								Moy_slope_1 = ((N_slope_1-1)*Moy_slope_1 + Frequ_slope_1)/N_slope_1;
								Moy_Ch4slope_1 = ((N_slope_1-1)*Moy_Ch4slope_1 + Ch4_slope)/N_slope_1;
								Slope_slope_1 = (Slope_slope_1*(N_slope_1-2) + 6*(Frequ_slope_1-Moy_slope_1)/N_slope_1)/(N_slope_1+1);
								Slope_Ch4slope_1 = (Slope_Ch4slope_1*(N_slope_1-2) + 6*(Ch4_slope-Moy_Ch4slope_1)/N_slope_1)/(N_slope_1+1);
								
								if ((utc - t1_1) > SlopeTime1) {
									Slope_1 = Slope_slope_1;
									Ch4Slope =  Slope_Ch4slope_1;
									
									N_slope_1 = 0;
									Frequ_slope_1 = 0.0;
									Moy_slope_1 = 0.0;
									Slope_slope_1 = 0.0;
									Ch4_slope = 0.0;
									Moy_Ch4slope_1 = 0.0;
									Slope_Ch4slope_1 = 0.0;

									// frep positive step
									DDS4xAD9912_FrequencyRampe(&DDS4xAD9912,1, FrequDDS1,(FrequDDS1+DeltakHz_1*1000), Step1/Ndiv);
									SetCtrlVal(MainPanel, PANEL_DDS1, (FrequDDS1+DeltakHz_1*1000));
									DDS4xAD9912_SetFrequency(&DDS4xAD9912, 1, (FrequDDS1+DeltakHz_1*1000));
								
									// allow counter to settle
									settling = 3;
								
									// next step
									Measuring_1 += 1;
								}
								break;
								
							case N_MEASUREMENT_ADJUST_FREQ_PLUS:
							case N_MEASUREMENT_ADJUST_FREQ_MINUS:
								// adjust DDS frequency to keep beatnote within the bandpass filter

								if (settling > 0) {
									settling--;
									break;
								}
								
								double fDDS2 = FrequencyDDSBesInit + 275000 - Ch4;
								SetCtrlVal(MainPanel, PANEL_DDS2, fDDS2) ;
								DDS4xAD9912_SetFrequency(&DDS4xAD9912, 2, fDDS2);
								
								// allow counter to settle
								settling = 3;
								
								// next step
								Measuring_1 += 1;
								break;								
								

							case N_MEASUREMENT_FREP_PLUS:
								// frep positive step
								
								if (settling > 0) {
									settling--;
									break;
								}
								
								if (t2_1 == 0.0)
									t2_1 = utc;
								
								Frepplus_1 = Frepplus_1 + Math1 - Slope_1 * (utc - t2_1);
								Ch4Plus = Ch4Plus + Ch4 - Ch4Slope * (utc - t2_1);
								n_1 += 1;
									
								if ((utc - t2_1) > DeltaT_1) {
									Frepplus_1 = Frepplus_1 / n_1;
									Ch4Plus = Ch4Plus / n_1; 
									n_1 = 0;
									
									// frep negative step
									DDS4xAD9912_FrequencyRampe(&DDS4xAD9912,1, (FrequDDS1+DeltakHz_1*1000),(FrequDDS1-DeltakHz_1*1000), Step1/Ndiv);
									SetCtrlVal(MainPanel, PANEL_DDS1, (FrequDDS1-DeltakHz_1*1000));  
									DDS4xAD9912_SetFrequency(&DDS4xAD9912, 1, (FrequDDS1-DeltakHz_1*1000));
									
									// allow counter to settle
									settling = 3;
								
									// next step
									Measuring_1 += 1;
								}
								break;
												
							case N_MEASUREMENT_FREP_MINUS:
								// frep negative step
								
								if (settling > 0) {
									settling--;
									break;
								}

								if (t3_1 == 0.0)
									t3_1 = utc;
								
								Frepminus_1 = Frepminus_1 + Math1 - Slope_1 * (utc - t3_1);
								Ch4Minus = Ch4Minus + Ch4 - Ch4Slope * (utc - t3_1);
								n_1 += 1;
								
								if ((utc - t3_1) > DeltaT_1) {
									Frepminus_1 = Frepminus_1 / n_1;
									Ch4Minus = Ch4Minus / n_1;
									n_1 = 0;
									
									t1_1 = 0.0;
									t2_1 = 0.0;
									t3_1 = 0.0;
									Frepminus_1 = 0.0;
									Frepplus_1 = 0.0;
									
									// compute N1
									N_1 = Signe1 * (2*Ndiv * DeltakHz_1 * 1000)/(Frepminus_1 - Frepplus_1 - Slope_1 * (t3_1 - t2_1));
									
									DDS4xAD9912_FrequencyRampe(&DDS4xAD9912, 1, FrequDDS1-DeltakHz_1*1000,FrequDDS1, Step1/Ndiv );
									SetCtrlVal(MainPanel, PANEL_DDS1, FrequDDS1);
									DDS4xAD9912_SetFrequency(&DDS4xAD9912, 1, FrequDDS1);
									SetCtrlVal(MainPanel, PANEL_DDS2, FrequencyDDSBesInit);
									DDS4xAD9912_SetFrequency(&DDS4xAD9912, 2, FrequencyDDSBesInit);
									
									// done
									Measuring_1 = N_MEASUREMENT_NONE;
								}
								break;
						}
						
						if (Measuring_2==TRUE)
							{
								if (Step1_2==FALSE) {
									SetCtrlVal(MainPanel, PANEL_DDS1, FrequDDS1) ;
									DDS4xAD9912_SetFrequency(&DDS4xAD9912, 1, FrequDDS1);
									
								    GetCtrlVal(MainPanel, PANEL_DDS2, &FrequencyDDSBesInit) ;
									GetCtrlVal(MainPanel, PANEL_DDS3, &FrequencyDDS3Init) ;
									t1_2=utc;
									Frequ_slope_2=Math1;
									Beat_slope_2=Ch2;
									Moy_slope_2= Frequ_slope_2;
									Moy_Beatslope_2= Beat_slope_2;
									N_slope_2=1;
									Step1_2=TRUE ;
									Nu1=(250000000+Math1)*N1;
									}
								else {
									if(Step2_2==FALSE){
										N_slope_2=N_slope_2+1; 
										Frequ_slope_2=Math1;
										Beat_slope_2=Ch2;  
										Moy_slope_2=((N_slope_2-1)*Moy_slope_2 + Frequ_slope_2)/N_slope_2;
										Moy_Beatslope_2=((N_slope_2-1)*Moy_Beatslope_2 + Beat_slope_2)/N_slope_2;
										Slope_slope_2 = (Slope_slope_2*(N_slope_2-2) + 6*(Frequ_slope_2-Moy_slope_2)/N_slope_2)/(N_slope_2+1) ;
										Slope_Beatslope_2 = (Slope_Beatslope_2*(N_slope_2-2) + 6*(Beat_slope_2-Moy_Beatslope_2)/N_slope_2)/(N_slope_2+1) ;
										if (utc-t1_2>SlopeTime2) {
											Slope_2 = Slope_slope_2;
											Beatslope_2 = Slope_Beatslope_2;
											Step2_2=TRUE ;
											N_slope_2=0;
											Frequ_slope_2=0.0;
											Moy_slope_2=0.0;
											Slope_slope_2 =0.0;
											Moy_Beatslope_2=0.0;
											Slope_Beatslope_2 =0.0;
											Beat_slope_2=0.0;
											
											DDS4xAD9912_FrequencyRampe (&DDS4xAD9912,1, FrequDDS1,(FrequDDS1+DeltakHz_2*1000), Step2/Ndiv ) ;
											SetCtrlVal(MainPanel, PANEL_DDS1, (FrequDDS1+DeltakHz_2*1000)) ;   
											DDS4xAD9912_SetFrequency(&DDS4xAD9912, 1, (FrequDDS1+DeltakHz_2*1000));
											Delay(0.1); 
											DeltaDDS3 = -DeltakHz_2*1000*(-Signe1/Signe2)*Ndiv*(Nu2)/(Nu1) - Beatslope_2*(utc-t1_2);
											SetCtrlVal(MainPanel, PANEL_DDS3, (FrequencyDDS3Init+DeltaDDS3)) ; 
											DDS4xAD9912_SetFrequency (&DDS4xAD9912,3,(FrequencyDDS3Init+DeltaDDS3)) ;
											
											}
										}
									else {
									  if (DDSBesChanged1==FALSE){
										  
										  		if (nDDSChange<3)
												{ nDDSChange=nDDSChange+1;}
												
												else
												{
												nDDSChange=0;
													
												double DeltaFrep275 = 275000-Ch4;
												double DeltaFrep10 = 10000-Ch2;
												DeltaDDS3= DeltaDDS3 +(DeltaFrep10);
												FrequencyDDSBes=FrequencyDDSBesInit+(DeltaFrep275);
												FrequencyDDS3=FrequencyDDS3Init+DeltaDDS3;
												
												SetCtrlVal(MainPanel, PANEL_DDS2, FrequencyDDSBes) ;
												DDS4xAD9912_SetFrequency(&DDS4xAD9912,2,FrequencyDDSBes);
												
												Delay(0.1); 
												SetCtrlVal(MainPanel, PANEL_DDS3, FrequencyDDS3) ;    
												DDS4xAD9912_SetFrequency(&DDS4xAD9912,3,FrequencyDDS3);
												 
												DDSBesChanged1=TRUE;
												t2_2=utc;
							
												}
										  	
									  }
									   else{
										  
										if(Step3_2==FALSE){
										  if (nstabilization<3)
										  {nstabilization= nstabilization+1;}
										  else
										  {
											if (utc-t2_2<DeltaT_2) {
													Frepplus_2=Frepplus_2 +Math1+250000000-Slope_2*(utc-t2_2);
													Delta10K_Plus= Delta10K_Plus + 10000 - (Ch2 -Beatslope_2*(utc-t2_2));
													n_2=n_2+1;
													}
											else{
												Frepplus_2=Frepplus_2/n_2;
												Delta10K_Plus=Delta10K_Plus/n_2;
												n_2=0;
												Step3_2=TRUE ;
												nstabilization=0;
												DDS4xAD9912_FrequencyRampe ( &DDS4xAD9912,1, (FrequDDS1+DeltakHz_2*1000),(FrequDDS1-DeltakHz_2*1000), Step2/Ndiv ) ;
												SetCtrlVal(MainPanel, PANEL_DDS1, (FrequDDS1-DeltakHz_2*1000)) ; 
												DDS4xAD9912_SetFrequency (&DDS4xAD9912,1, (FrequDDS1-DeltakHz_2*1000) ) ;
												
												Delay(0.1); 
												DeltaDDS3 = (FrequencyDDS3Init+DeltakHz_2*1000*(-Signe1/Signe2)*Ndiv*(Nu2)/(Nu1)) - FrequencyDDS3;
												SetCtrlVal(MainPanel, PANEL_DDS3, FrequencyDDS3+DeltaDDS3) ; 
												DDS4xAD9912_SetFrequency (&DDS4xAD9912,3, FrequencyDDS3+DeltaDDS3 ) ; 
												
												}
										  }
										}
										
										else {
											if (DDSBesChanged2==FALSE){
												
												if (nDDSChange<3)
												{ nDDSChange=nDDSChange+1;}
												
												else
												{
												nDDSChange=0;
										
												
												double DeltaFrep275=275000-Ch4;
												double DeltaFrep10=10000-Ch2;
												DeltaDDS3 = DeltaDDS3+(DeltaFrep10) ;
												FrequencyDDSBes=FrequencyDDSBes+(DeltaFrep275) ;
												FrequencyDDS3=FrequencyDDS3+DeltaDDS3 ; 
												DDS4xAD9912_SetFrequency(&DDS4xAD9912,2,FrequencyDDSBes);
												SetCtrlVal(MainPanel, PANEL_DDS2, FrequencyDDSBes);
												
												Delay(0.1); 
												SetCtrlVal(MainPanel, PANEL_DDS3, FrequencyDDS3) ; 
												DDS4xAD9912_SetFrequency(&DDS4xAD9912,3,FrequencyDDS3);
												
												DDSBesChanged2=TRUE;
												t3_2=utc;

						
												}
											}
											
											else{
											  if (nstabilization<3)
											  { nstabilization=nstabilization+1;}
											  else
											  {
												if (utc-t3_2<DeltaT_2) {
													Frepminus_2=Frepminus_2 +Math1+250000000-Slope_2*(utc-t3_2);
													Delta10K_Minus=  Delta10K_Minus +10000 - ( Ch2 -Beatslope_2*(utc-t3_2));
													n_2=n_2+1;
													}
												else{
													Frepminus_2=Frepminus_2/(n_2);
													Delta10K_Minus= Delta10K_Minus/n_2;
													N_2 = (Signe2)*(-DeltaDDS3+Delta10K_Plus-Delta10K_Minus-Beatslope_2*(t3_2-t2_2) )/(Frepminus_2-Frepplus_2-Slope_2*(t3_2-t2_2));
													n_2=0;
													Frepminus_2=0.0;
													Frepplus_2=0.0;
													Delta10K_Minus=0.0;
													Delta10K_Plus=0.0;  
													DDS4xAD9912_FrequencyRampe (&DDS4xAD9912, 1, FrequDDS1-DeltakHz_2*1000,FrequDDS1, Step2/Ndiv ) ;
													SetCtrlVal(MainPanel, PANEL_DDS1, FrequDDS1) ;
													DDS4xAD9912_SetFrequency(&DDS4xAD9912,1,FrequDDS1);
													
													Delay(0.1); 
													
													SetCtrlVal(MainPanel, PANEL_DDS2, FrequencyDDSBesInit) ;  
													DDS4xAD9912_SetFrequency (&DDS4xAD9912, 2, FrequencyDDSBesInit ) ;
													
													Delay(0.1); 
													
													SetCtrlVal(MainPanel, PANEL_DDS3, FrequencyDDS3Init-Beatslope_2*(utc-t1_2)) ;
													DDS4xAD9912_SetFrequency (&DDS4xAD9912, 3, FrequencyDDS3Init-Beatslope_2*(utc-t1_2) ) ;
													
													Measuring_2=FALSE ;
													Step1_2=FALSE ;
													Step2_2=FALSE ;
													Step3_2=FALSE ;
													t1_2=0.0;
													t2_2=0.0;
													t3_2=0.0;
													DDSBesChanged1=FALSE;
													DDSBesChanged2=FALSE;
													FrequencyDDSBes=0.0;
													nstabilization=0;
													
													}
											  }
											}
											}
									  }
										}
									}
							}
						
						switch (Measuring_3) {

							case N_MEASUREMENT_STEP_0:
								// not measuring N3
								break;
							
							case N_MEASUREMENT_STEP_1: // init
								
								SetCtrlVal(MainPanel, PANEL_DDS4, FrequDDS4);
								DDS4xAD9912_SetFrequency(&DDS4xAD9912, 4, FrequDDS4);
								settling = 3;
								
								t1_3 = utc;
								N_slope_3 = 0;
								// record current DDS3 frequency
								GetCtrlVal(MainPanel, PANEL_DDS3, &FrequencyDDS3Init);

								 // step 1 done
								Measuring_3 = N_MEASUREMENT_STEP_2;
								break;
								
							case N_MEASUREMENT_STEP_2: // slope measurement
								
								if (settling > 0) {
									settling--;
									break;
								}
								
								N_slope_3++;
								Frequ_slope_3 = Ch2;
								Moy_slope_3 = ((N_slope_3-1)*Moy_slope_3 + Frequ_slope_3)/N_slope_3;
								Slope_slope_3 = (Slope_slope_3*(N_slope_3-2) + 6*(Frequ_slope_3-Moy_slope_3)/N_slope_3)/(N_slope_3+1) ;
								
								if (utc - t1_3 > SlopeTime3) {
									// slope measurement
									Slope_3 = Slope_slope_3;
									
									t2_3 = utc;
									N_slope_3 = 0;
									Frequ_slope_3 = 0.0;
									Moy_slope_3 = 0.0;
									Slope_slope_3 = 0.0;
									
									// step 2 done
									Measuring_3 = N_MEASUREMENT_STEP_3;
									
									// frep positive step
									SetCtrlVal(MainPanel, PANEL_DDS4, FrequDDS4 + DeltakHz_3 * 1000);
									DDS4xAD9912_SetFrequency(&DDS4xAD9912, 4, FrequDDS4 + DeltakHz_3 * 1000);
									// compensate with DDS3 to keep measured beatnote in counter box range
									double fDDS3 = FrequencyDDS3Init + N3/N1 * Ndiv * DeltakHz_3 * 1000;
									SetCtrlVal(MainPanel, PANEL_DDS3, fDDS3);
									DDS4xAD9912_SetFrequency(&DDS4xAD9912, 3, fDDS3);
									// allow counter to settle
									settling = 3;
								}
								break;
								
							case N_MEASUREMENT_STEP_3: // frep positive step
								
								if (settling > 0) {
									settling--;
									break;
								}
								
								n_3++;
								Frepplus_3 += Ch2 - Slope_3 * (utc - t2_3);
									
								if (utc - t2_3 > DeltaT_3) {
									// positive step measurement
									Frepplus_3 = Frepplus_3 / n_3;
									
									n_3 = 0;
									t3_3 = utc;
									
									// step 3 done
									Measuring_3 = N_MEASUREMENT_STEP_4;
									
									// frep negative step
									SetCtrlVal(MainPanel, PANEL_DDS4, FrequDDS4 - DeltakHz_3 * 1000);
									DDS4xAD9912_SetFrequency(&DDS4xAD9912, 4, FrequDDS4 - DeltakHz_3 * 1000);
									// compensate with DDS3 to keep measured beatnote in counter box range
									double fDDS3 = FrequencyDDS3Init - N3/N1 * Ndiv * DeltakHz_3 * 1000;
									SetCtrlVal(MainPanel, PANEL_DDS3, fDDS3);
									DDS4xAD9912_SetFrequency(&DDS4xAD9912, 3, fDDS3);
									// allow counter to settle
									settling = 3;
								}
								break;
								
								
							case N_MEASUREMENT_STEP_4: // frep negative step
								
								if (settling > 0) {
									settling--;
									break;
								}
								
								n_3++;
								Frepminus_3 += Ch2 - Slope_3 * (utc - t3_3);
												
								if (utc - t3_3 > DeltaT_3) {
									// positive step measurement
									Frepminus_3 = Frepminus_3 / n_3;
									
									// compute N3
									N_3 = 1000.0 * DeltakHz_3 / (Frepplus_3 - Frepminus_3 + (2 * N3/N1 * Ndiv * DeltakHz_3 * 1000));
									SetCtrlVal(CalcN3Panel, CALCN3_N, N_3);
									
									t1_3=0.0;
									t2_3=0.0;
									t3_3=0.0;
									n_3 = 0;
									Frepminus_3 = 0.0;
									Frepplus_3 = 0.0;
										
									// step 4 done
									Measuring_3 = N_MEASUREMENT_STEP_0;

									// back to nominal frep
									SetCtrlVal(MainPanel, PANEL_DDS4, FrequDDS4);
									DDS4xAD9912_SetFrequency(&DDS4xAD9912, 4, FrequDDS4);
									// back to initial DDS3 frequency
									SetCtrlVal(MainPanel, PANEL_DDS3, FrequencyDDS3Init);
									DDS4xAD9912_SetFrequency(&DDS4xAD9912, 3, FrequencyDDS3Init);
								}
								break;
						}
						
						// Calcul du signe de fb
						
						if (Getsign1 == TRUE) {
							if (utc > tbegin1+2) {
								if (Math1 > Frepbefore1)
									Signe1 = -1.0;
								else
									Signe1 = +1.0;
								SetCtrlVal(MainPanel, PANEL_DDS1, Frequency1) ; 
								DDS4xAD9912_SetFrequency(&DDS4xAD9912, 1, Frequency1);
								Getsign1 = FALSE;
							}
						}
						if (Getsign2 == TRUE) {
							if (utc > tbegin2+2){
								if (Math1 > Frepbefore2) { 
									if (Ch2 > Ch2before)
										Signe2 = +1.0;
									else
										Signe2 = -1.0;
								} else {
									if (Ch2 > Ch2before)
										Signe2 = -1.0;
									else 
										Signe2 = +1.0;
								}
								SetCtrlVal(MainPanel, PANEL_DDS1, Frequency2) ; 
								DDS4xAD9912_SetFrequency(&DDS4xAD9912, 1, Frequency2);
								Getsign2 = FALSE;
							}
						}
						if (Getsign3 == TRUE) {
							if (utc > tbegin3+2){
								if (Ch3 > Frepbefore3)
									Signe3 = -1.0;
								else
									Signe3 = +1.0;
								SetCtrlVal(MainPanel, PANEL_DDS3, Frequency3) ; 
								DDS4xAD9912_SetFrequency(&DDS4xAD9912, 4, Frequency3);
								Getsign3 = FALSE;
							}
						}
								
						
						
						
						/*
						 * Slope Cancelling 
						 */
						
						if (SlopeMeasuring == TRUE)
						{
							double currentFreq = 0.0;
							
							/* select reference */
							switch (slopeReference) {
								case SLOPE_REFERENCE_MICROWAVE:
									currentFreq = Math2;
									break;
								case SLOPE_REFERENCE_HG_CAVITY:
									currentFreq = Ch2 * 1062.5 / 1542.2;
									break;
							}
							
							if  (utc-SlopeMeasuringTimeBegin > TimetoSlope)
							{
								Nratio = Nratio + 1;
								
								if (Nratio >= 1)
								{
									MoyMath2 = MoyMath2 + Moy_Math2slope;
								}
								SlopeMath2 = SlopeMath2 + Slope_Math2slope;
								N_Math2slope = 0.0;
								Math2_slope = 0.0;
								Moy_Math2slope = 0.0;
								Slope_Math2slope = 0.0;
								
								if (Nratio == 1 && CenterFrequencyCh2ToDetermine == TRUE)
								{
									CenterFrequencyCh2 = MoyMath2;
									CenterFrequencyCh2ToDetermine = FALSE;
								}

								OnSlopeCancelling = TRUE;
								
								FoxFrequ = DDSFox_ReadFreq(&DDS1xAD9956);
								
								if (Nratio == ratio)
								{
									if (FrequCorrec == TRUE)
									{
										SlopeCorrection = (MoyMath2/ratio-CenterFrequencyCh2)/TimetoSlope;
										SlopeMath2 = SlopeMath2 + SlopeCorrection;
									}
									Nratio = 0;
									MoyMath2 = 0.0;
								}
								
								DDSFox_Set(&DDS1xAD9956, FoxFrequ, SlopeMath2);
								
								nstabilisationSlopeMeasuring = 0;
								SlopeMeasuringTimeBegin = utc;
								
							}
							else
							{
								if (nstabilisationSlopeMeasuring < 5)
								{
									nstabilisationSlopeMeasuring = nstabilisationSlopeMeasuring + 1;
									Math2_slope = currentFreq; 
								}
								else
								{
									if ((currentFreq-Math2_slope) < limitotakoff && (currentFreq-Math2_slope) > -limitotakoff)
									{
										N_Math2slope = N_Math2slope + 1;
										Math2_slope = currentFreq;
										Moy_Math2slope = ((N_Math2slope-1)*Moy_Math2slope + Math2_slope)/N_Math2slope;
										Slope_Math2slope = (Slope_Math2slope*(N_Math2slope-2) + 6*(Math2_slope-Moy_Math2slope)/N_Math2slope)/(N_Math2slope+1) ;
									}
									else
									{
										if (AutoStopSlopeCancellingIfDelock)
										{
											
											/* 
											 * stop slope cancelling if the comb is not locked
											 */
											
											if (KeepFrequ == TRUE)
											{
												ActualFrequ = DDSFox_ReadFreq(&DDS1xAD9956);
											}
											else
											{
												ActualFrequ = DEDRIFT_DDS_FREQUENCY;
											}  
												
											if (KeepSlope == TRUE)
											{
												ActualSlope = SlopeMath2;
											}
											else
											{
												ActualSlope = 0.0;
											    OnSlopeCancelling = FALSE;
											} 
												
											DDSFox_Set(&DDS1xAD9956, ActualFrequ, ActualSlope);
											
											SlopeMeasuring = FALSE;
											N_Math2slope = 0.0;
											Math2_slope = 0.0;
											MoyMath2 = 0.0; 
											Moy_Math2slope = 0.0;
											Slope_Math2slope = 0.0;
											CenterFrequencyCh2 = 0.0;
											CenterFrequencyCh2ToDetermine = TRUE;
											Nratio = -1;
											nstabilisationSlopeMeasuring = 0;
											
											SetCtrlVal(MainPanel, PANEL_STARTCANCEL, 0);
										}
									}
								}
							}
						}
						
						
						////   For re-centering
						
						if (KeepCentering) 
						{
							
							DeltaCh4=275000-Ch4;
							DeltaCh2=10000-Ch2;
							
							if (utc- CenteringTimeBegin275K > Timetorecenter275K && CenteringTimeBegin275K>10)
							{
								 GetCtrlVal(MainPanel, PANEL_DDS2, &FrequencyToChange) ;
								 SetCtrlVal(MainPanel, PANEL_DDS2, FrequencyToChange+DeltaCh4) ;
					 		   	 DDS4xAD9912_SetFrequency(&DDS4xAD9912, 2, FrequencyToChange+DeltaCh4);
								 CenteringTimeBegin275K=utc;
								
							
							}
							
							if (utc- CenteringTimeBegin10K > Timetorecenter10K && CenteringTimeBegin10K>10)
							{
								  GetCtrlVal(MainPanel, PANEL_DDS3, &FrequencyToChange) ;
								  SetCtrlVal(MainPanel, PANEL_DDS3, FrequencyToChange+DeltaCh2) ; 
					 			  DDS4xAD9912_SetFrequency(&DDS4xAD9912, 3, FrequencyToChange+DeltaCh2);
								  CenteringTimeBegin10K=utc; 
								
								
							}
							
							
						}
							
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_MATH2AUTOSAV, &BoxChecked) ;		   // AutoSave OL
						if (BoxChecked) {
							SetCtrlVal(MainPanel, PANEL_CHECKBOX_MATH2SAVE, TRUE) ;   // so that it will try to write it (at next block) if it seems reasonnable, even though it was off before
						}
						
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_MATH2SAVE, &BoxChecked) ;			   // Save OL (Math2)
						if (BoxChecked) {
							FileOpt = OpenFile("z:\\MeasuresFifi1\\OptCavity.txt", VAL_WRITE_ONLY, VAL_APPEND, VAL_ASCII) ;
							Fmt(ReportString, "%s\t%s\t%f[p3]\t%f[p3]", Date, Time, utc, Math2);
							WriteLine(FileOpt, ReportString, -1) ;
							CloseFile(FileOpt) ;
							FileOpt = OpenFile("C:\\Femto\\Results\\OptCavity.txt", VAL_WRITE_ONLY, VAL_APPEND, VAL_ASCII) ;   // a local backup for debugging
							Fmt(ReportString, "%s\t%s\t%f[p3]\t%f[p3]", Date, Time, utc, Math2);
							WriteLine(FileOpt, ReportString, -1) ;
							CloseFile(FileOpt) ;
						}
						
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_MATH3AUTOSAV, &BoxChecked) ;		   // AutoSave Hg (Math3)
						if (BoxChecked) {
							SetCtrlVal(MainPanel, PANEL_CHECKBOX_MATH3SAVE, TRUE) ;   // so that it will try to write it (at next block) if it seems reasonnable, even though it was off before
						}
						
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_MATH3SAVE, &BoxChecked) ;			   // Save Hg
						
						if (BoxChecked) {
							FileHg = OpenFile("z:\\MeasuresFifi1\\HgCavity.txt", VAL_WRITE_ONLY, VAL_APPEND, VAL_ASCII) ;
							Fmt(ReportString, "%s\t%s\t%f[p3]\t%f[p3]", Date, Time, utc, Math3);
							WriteLine(FileHg, ReportString, -1) ;
							CloseFile(FileHg) ;
							FileHg = OpenFile("C:\\Femto\\Results\\HgCavity.txt", VAL_WRITE_ONLY, VAL_APPEND, VAL_ASCII) ;
							Fmt(ReportString, "%s\t%s\t%f[p3]\t%f[p3]", Date, Time, utc, Math3);
							WriteLine(FileHg, ReportString, -1) ;
							CloseFile(FileHg) ;
						}
						
						GetCtrlVal(MainPanel, PANEL_CHECKBOX_MATH5SAVE, &BoxChecked) ;				// Save ExtraMath (Math5)
						if (BoxChecked) {
							FileExtraMath = OpenFile(ExtraMathFileName, VAL_WRITE_ONLY, VAL_APPEND, VAL_ASCII) ;
							Fmt(ReportString, "%s\t%s\t%f[p3]\t%f[p7]", Date, Time, utc, Math5);
							WriteLine(FileExtraMath, ReportString, -1) ;
							CloseFile(FileExtraMath) ;
						}
						
						// Special case to handle change of day at next second
						if ( LocalTime.tm_hour==23 && LocalTime.tm_min==59 && strtod(Sec,NULL)>=58 ) { 
							printf("Entering change of date substructure\r\n"); 
							printf("old file : %s %s %s \r\n",Date, Time, LogFileName);
							Acquiring = FALSE ;
							do {
								Delay(5.1);
								CurrentFileName(LogFileName) ;
							} while (!GetFileInfo(LogFileName, &OldLogFilePtr));
							printf("new file : %s \r\n\r\n",Date, Time, LogFileName);  
							Acquiring = TRUE ;
							OldLogFilePtr = 2;
						}
						
						ResumeTimerCallbacks() ; 
						
						} ;
					
					break;
				case FALSE:
					break;
				}
			break;
		}
	return 0;
}




int CVICALLBACK CB_OnFreqPlot (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	int BoxChecked ;
	Plot_Data * pPlot = NULL;
	char PlotTitle[] = "Ch? Frequency Plot" ;
	double PlotMin = 10e6 ;
	double PlotMax = 65e6 ;
	
	switch (event)
		{
		case EVENT_COMMIT:
		
			GetCtrlVal(MainPanel, control, &BoxChecked);
			switch (control) {
				case PANEL_CHECKBOX_FREQ1PLOT :
					pPlot = &PlotCh1 ;
					Fmt(PlotTitle, "Ch1 Frequency Plot") ;
					PlotMin = 54.999e6 ; PlotMax = 55.001e6 ;
					break ;
				case PANEL_CHECKBOX_FREQ2PLOT :
					pPlot = &PlotCh2 ;
					Fmt(PlotTitle, "Ch2 Frequency Plot") ;
					PlotMin = 8.0e3 ; PlotMax = 12.0e3 ;
					break ;
				case PANEL_CHECKBOX_FREQ3PLOT :
					pPlot = &PlotCh3 ;
					Fmt(PlotTitle, "Ch3 Frequency Plot") ;
					PlotMin = 8.0e3 ; PlotMax = 12.0e3 ;
					break ;
				case PANEL_CHECKBOX_FREQ4PLOT :
					pPlot = &PlotCh4 ;
					Fmt(PlotTitle, "Ch4 Frequency Plot") ;
					break ;
				case PANEL_CHECKBOX_MATH1PLOT :
					pPlot = &PlotMath1 ;
					Fmt(PlotTitle, "Math1 Plot") ;
					PlotMin = 765.0e6 ; PlotMax = 775.0e6 ;
					break ;
				case PANEL_CHECKBOX_MATH2PLOT :
					pPlot = &PlotMath2 ;
					Fmt(PlotTitle, "Math2 Plot") ;
					PlotMin = -1.0e9 ; PlotMax = 1.0e9 ;
					break ;	
				case PANEL_CHECKBOX_MATH3PLOT :
					pPlot = &PlotMath3 ;
					Fmt(PlotTitle, "Math3 Plot") ;
					PlotMin = -1.0e9 ; PlotMax = 1.0e9 ;
					break ;	
				case PANEL_CHECKBOX_MATH4PLOT :
					pPlot = &PlotMath4 ;
					Fmt(PlotTitle, "Math4 Plot") ;
					PlotMin = -1.0e9 ; PlotMax = 1.0e9 ;
					break ;	
				case PANEL_CHECKBOX_MATH5PLOT :
					pPlot = &PlotMath5 ;
					Fmt(PlotTitle, "Math5 Plot") ;
					PlotMin = -1.0e9 ; PlotMax = 1.0e9 ;
					break ;
			}
			
			if (BoxChecked) {
				Plot_InitPanel(pPlot, PlotTitle, PlotMin, PlotMax, &OnCloseViewPanel) ;	
				}
			else {
				Plot_ClosePanel(pPlot) ;
				} ;
			break;
			
		case EVENT_RIGHT_CLICK:

			break;
			
		}
	return 0;
}


int CVICALLBACK CB_OnAllanPlot (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	int BoxChecked ;
	Allan_Data * pAllan = NULL;
	char AllanTitle[] = "Ch? Allan Deviation  " ;
	double Normalizer = 300e12 ;
	
	switch (event)
		{
		case EVENT_COMMIT:
			
			GetCtrlVal(MainPanel, control, &BoxChecked);
			switch (control) {
				case PANEL_CHECKBOX_FREQ1ALLAN :
					pAllan = &AllanCh1 ;
					Fmt(AllanTitle, "Ch1 Allan Deviation") ;
					Normalizer = 1.84e12 ;
					break ;
				case PANEL_CHECKBOX_FREQ2ALLAN :
					pAllan = &AllanCh2 ;
					Fmt(AllanTitle, "Ch2 Allan Deviation") ;
					Normalizer = 10.e3 ; 
					break ;
				case PANEL_CHECKBOX_FREQ3ALLAN :
					pAllan = &AllanCh3 ;
					Fmt(AllanTitle, "Ch3 Allan Deviation") ;
					Normalizer = 429.228e12 ; 
					break ;
				case PANEL_CHECKBOX_FREQ4ALLAN :
					pAllan = &AllanCh4 ;
					Fmt(AllanTitle, "Ch4 Allan Deviation") ;
					Normalizer = 275.0e3 ; 
					break ;
				case PANEL_CHECKBOX_MATH1ALLAN :
					pAllan = &AllanMath1 ;
					Fmt(AllanTitle, "Math1 Allan Deviation") ;
					Normalizer = 250.0e6 ;
					break ;
				case PANEL_CHECKBOX_MATH2ALLAN :
					pAllan = &AllanMath2 ;
					Fmt(AllanTitle, "Math2 Allan Deviation") ;
					Normalizer = 194.395e12 ;
					break ;
				case PANEL_CHECKBOX_MATH3ALLAN :
					pAllan = &AllanMath3 ;
					Fmt(AllanTitle, "Math3 Allan Deviation") ;
					Normalizer = 282.143e12 ;
					break ;
				case PANEL_CHECKBOX_MATH4ALLAN :
					pAllan = &AllanMath4 ;
					Fmt(AllanTitle, "Math4 Allan Deviation") ;
					Normalizer = 429.228e12 ;
					break ;
				case PANEL_CHECKBOX_MATH5ALLAN :
					pAllan = &AllanMath5 ;
					Fmt(AllanTitle, "Math5 Allan Deviation") ;
					Normalizer = 429.228e12 ;
					break ;
			}			
			
			if (BoxChecked) {
				Allan_InitPanel(pAllan, AllanTitle, Normalizer, &OnCloseViewPanel) ;	
				}
			else {
				Allan_ClosePanel(pAllan) ;
				} ;
			break;
			
		case EVENT_RIGHT_CLICK:

			break;
			
		}
	return 0;
}

int  CVICALLBACK CB_ChangeDDSOut (int panel, int control, int event, 
		 void *callbackData, int eventData1, int eventData2)
{
	double frequency ;
	
	switch (event) {
		case EVENT_COMMIT:
			GetCtrlVal(MainPanel, control, &frequency);
			
			switch (control) {
				case PANEL_DDS1:
					DDS4xAD9912_SetFrequency(&DDS4xAD9912, 1, frequency);
					break;
				case PANEL_DDS2:
					DDS4xAD9912_SetFrequency(&DDS4xAD9912, 2, frequency);
					break;
				case PANEL_DDS3:
					DDS4xAD9912_SetFrequency(&DDS4xAD9912, 3, frequency);
					break;
				case PANEL_DDS4:
					DDS4xAD9912_SetFrequency(&DDS4xAD9912, 4, frequency);
					break;
			}
	}
	return 0;
}

int  CVICALLBACK CB_ChangeDDSStep (int panel, int control, int event, 
		 void *callbackData, int eventData1, int eventData2)
{
	double Step ;
	
	switch (event)
		{
		case EVENT_COMMIT:
			GetCtrlVal(MainPanel, control, &Step);
			if (control==PANEL_DDS1STEP)   { SetCtrlAttribute(panel, PANEL_DDS1,   ATTR_INCR_VALUE, Step); }   
			if (control==PANEL_DDS2STEP)   { SetCtrlAttribute(panel, PANEL_DDS2,   ATTR_INCR_VALUE, Step); }
			if (control==PANEL_DDS3STEP)   { SetCtrlAttribute(panel, PANEL_DDS3,   ATTR_INCR_VALUE, Step); }   
			if (control==PANEL_DDS4STEP)   { SetCtrlAttribute(panel, PANEL_DDS4,   ATTR_INCR_VALUE, Step); } 
			break;
		case EVENT_RIGHT_CLICK:

			break;
			
		}
	return 0;
}

int CVICALLBACK CB_ChangeMath (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	int len;
	char *string;
		
	switch (event)
	{
		case EVENT_COMMIT:
			GetCtrlAttribute(panel, control, ATTR_STRING_TEXT_LENGTH, &len);
			string = (char *)malloc(sizeof(char) * (len + 1));
			GetCtrlVal(panel, control, string);
			switch (control) {
				case PANEL_MATHSTRING1:
					mupSetExpr(MathParser1, string);
					break;
				case PANEL_MATHSTRING2:
					mupSetExpr(MathParser2, string);
					break;
				case PANEL_MATHSTRING3:
					mupSetExpr(MathParser3, string);
					break;
				case PANEL_MATHSTRING4:
					mupSetExpr(MathParser4, string);
					break;
				case PANEL_MATHSTRING5:
					mupSetExpr(MathParser5, string);
					break;
			}
			free(string);
			break;
	}
	return 0;
}

int CVICALLBACK CB_ChangeN (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{

	switch (event)
		{
		case EVENT_COMMIT:
			if (control==PANEL_N1CHOICE) GetCtrlVal(MainPanel, control, &N1) ;
			if (control==PANEL_N2CHOICE) GetCtrlVal(MainPanel, control, &N2) ;
			if (control==PANEL_N3CHOICE) GetCtrlVal(MainPanel, control, &N3) ;
			break;
		}
	return 0;
}

int CVICALLBACK CB_OnPlus10k (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	double Frequency ;
	
	switch (event)
	{
		case EVENT_COMMIT:
			switch (control)
			{
				case PANEL_PLUS10KDDS1:
					GetCtrlVal(MainPanel, PANEL_DDS1, &Frequency);
					Frequency += 10000.0;
					DDS4xAD9912_SetFrequency(&DDS4xAD9912, 1, Frequency);	
					SetCtrlVal(MainPanel, PANEL_DDS1, Frequency);
					break;
				case PANEL_PLUS10KDDS2:
					GetCtrlVal(MainPanel, PANEL_DDS2, &Frequency);
					Frequency += 10000.0;
					DDS4xAD9912_SetFrequency(&DDS4xAD9912, 2, Frequency);	
					SetCtrlVal(MainPanel, PANEL_DDS2, Frequency);
					break;
				case PANEL_PLUS10KDDS3:
					GetCtrlVal(MainPanel, PANEL_DDS3, &Frequency);
					Frequency += 10000.0;
					DDS4xAD9912_SetFrequency(&DDS4xAD9912, 3, Frequency);	
					SetCtrlVal(MainPanel, PANEL_DDS3, Frequency);
					break;
				case PANEL_PLUS10KDDS4:
					GetCtrlVal(MainPanel, PANEL_DDS4, &Frequency);
					Frequency += 10000.0;
					DDS4xAD9912_SetFrequency(&DDS4xAD9912, 4, Frequency);	
					SetCtrlVal(MainPanel, PANEL_DDS4, Frequency);
					break;
			}
			break;
	}
	return 0;
}


int CVICALLBACK CB_OnMinus10k (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	double Frequency;
	
	switch (event)
	{
		case EVENT_COMMIT:
			switch (control)
			{
				case PANEL_MINUS10KDDS1:
					GetCtrlVal(MainPanel, PANEL_DDS1, &Frequency);
					Frequency -= 10000.0;
					DDS4xAD9912_SetFrequency(&DDS4xAD9912, 1, Frequency);
					SetCtrlVal(MainPanel, PANEL_DDS1, Frequency);
					break;
				case PANEL_MINUS10KDDS2:
					GetCtrlVal(MainPanel, PANEL_DDS2, &Frequency);
					Frequency -= 10000.0;
					DDS4xAD9912_SetFrequency(&DDS4xAD9912, 2, Frequency);	
					SetCtrlVal(MainPanel, PANEL_DDS2, Frequency);
					break;
				case PANEL_MINUS10KDDS3:
					GetCtrlVal(MainPanel, PANEL_DDS3, &Frequency);
					Frequency -= 10000.0;
					DDS4xAD9912_SetFrequency(&DDS4xAD9912, 3, Frequency);	
					SetCtrlVal(MainPanel, PANEL_DDS3, Frequency);
					break;
				case PANEL_MINUS10KDDS4:
					GetCtrlVal(MainPanel, PANEL_DDS4, &Frequency);
					Frequency -= 10000.0;
					DDS4xAD9912_SetFrequency(&DDS4xAD9912, 4, Frequency);
					SetCtrlVal(MainPanel, PANEL_DDS4, Frequency);
					break;
			}
			break;
	}
	return 0;
}



int CVICALLBACK CB_OnSaveOpt (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	bool BoxChecked = FALSE ;
	
	switch (event)
		{
		case EVENT_COMMIT:
			
			GetCtrlVal(MainPanel, control, &BoxChecked);
			
			if (BoxChecked) {
				SetCtrlAttribute(MainPanel, PANEL_CHECKBOX_MATH2AUTOSAV, ATTR_DIMMED, FALSE) ;
				}
			
			if (!BoxChecked) {
				SetCtrlAttribute(MainPanel, PANEL_CHECKBOX_MATH2AUTOSAV, ATTR_DIMMED, TRUE) ;
				}
			
			break;
		}
	return 0;
}



 
int CVICALLBACK CB_OnSaveHg (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	bool BoxChecked = FALSE ;
	
	switch (event)
		{
		case EVENT_COMMIT:
			
			GetCtrlVal(MainPanel, control, &BoxChecked);
			
			if (BoxChecked) {
				SetCtrlAttribute(MainPanel, PANEL_CHECKBOX_MATH3AUTOSAV, ATTR_DIMMED, FALSE) ;
				}
			
			if (!BoxChecked) {
				SetCtrlAttribute(MainPanel, PANEL_CHECKBOX_MATH3AUTOSAV, ATTR_DIMMED, TRUE) ;
				}
			
			break;
		}
	return 0;
}

int CVICALLBACK CB_OnAutoSaveHg (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	bool BoxChecked = FALSE ;
	
	switch (event)
		{
		case EVENT_COMMIT:
			
			GetCtrlVal(MainPanel, control, &BoxChecked);
			
			if (BoxChecked) {
				SetCtrlAttribute(MainPanel, PANEL_CHECKBOX_MATH2SAVE, ATTR_DIMMED, TRUE) ;
				}
			
			if (!BoxChecked) {
				SetCtrlAttribute(MainPanel, PANEL_CHECKBOX_MATH2SAVE, ATTR_DIMMED, FALSE) ;
				}

			break;
		}
	return 0;
}

int CVICALLBACK CB_OnSaveSr (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	bool BoxChecked = FALSE ;
	
	switch (event)
		{
		case EVENT_COMMIT:
			
			GetCtrlVal(MainPanel, control, &BoxChecked);
			
			if (BoxChecked) {
				SetCtrlAttribute(MainPanel, PANEL_CHECKBOX_MATH3AUTOSAV, ATTR_DIMMED, FALSE) ;
				}
			
			if (!BoxChecked) {
				SetCtrlAttribute(MainPanel, PANEL_CHECKBOX_MATH3AUTOSAV, ATTR_DIMMED, TRUE) ;        
				}
			break;
		}
	return 0;
}

int CVICALLBACK CB_OnAutoSaveSr (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	bool BoxChecked = FALSE ; 
	
	switch (event)
		{
		case EVENT_COMMIT:

			GetCtrlVal(MainPanel, control, &BoxChecked);
			
			if (BoxChecked) {
				SetCtrlAttribute(MainPanel, PANEL_CHECKBOX_MATH3SAVE, ATTR_DIMMED, TRUE) ;
				}
			
			if (!BoxChecked) {
				SetCtrlAttribute(MainPanel, PANEL_CHECKBOX_MATH3SAVE, ATTR_DIMMED, FALSE) ;
				}
			
			break;
		}
	return 0;
}


int CVICALLBACK CB_OnAcceptN (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	switch (event)
	{
		case EVENT_COMMIT:
			if (panel == CalcN1Panel) {
				N1 = round(N_1);
				SetCtrlVal(MainPanel, PANEL_N1CHOICE, N1);
			}
			if (panel == CalcN2Panel) {
				N2 = round(N_2);
				SetCtrlVal(MainPanel, PANEL_N2CHOICE, N2);
			}
			if (panel == CalcN3Panel) {
				N3 = round(N_3);
				SetCtrlVal(MainPanel, PANEL_N3CHOICE, N3);
			} 
			break;
	}
	return 0;
}


int CVICALLBACK OnChooseSaveFile (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	switch (event)
		{
		case EVENT_COMMIT:

			FileSelectPopup("D:\\Manip Femto", "*.txt", "*.txt", "Selection of the ExtraMath saving file", VAL_SAVE_BUTTON, 0, 0, 0, FALSE, ExtraMathFileName) ;
			SetCtrlAttribute(MainPanel, PANEL_CHECKBOX_MATH5SAVE, ATTR_DIMMED, FALSE) ;
			
			break;
		}
	return 0;
}


int CVICALLBACK CB_OnNCalculus (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{

	int PanelIsVisible, PanelTop, PanelLeft, PanelWidth ;
	
	switch (event)
		{
		case EVENT_COMMIT:
			
			if (control==PANEL_N1CALCULUS){
				GetPanelAttribute (CalcN1Panel, ATTR_VISIBLE , &PanelIsVisible);
				if (PanelIsVisible) {
				}
				else {
					GetPanelAttribute (MainPanel, ATTR_TOP,   &PanelTop) ;
					GetPanelAttribute (MainPanel, ATTR_LEFT,  &PanelLeft) ;
					GetPanelAttribute (MainPanel, ATTR_WIDTH, &PanelWidth) ;
					SetPanelAttribute (CalcN1Panel, ATTR_TOP, PanelTop+170 ) ;
					SetPanelAttribute (CalcN1Panel, ATTR_LEFT, PanelLeft+PanelWidth+4 ) ;
					DisplayPanel(CalcN1Panel) ;
				} 
				}
			if (control==PANEL_N2CALCULUS){
				GetPanelAttribute (CalcN2Panel, ATTR_VISIBLE , &PanelIsVisible);
				if (PanelIsVisible) {
				}
				else {
					GetPanelAttribute (MainPanel, ATTR_TOP,   &PanelTop) ;
					GetPanelAttribute (MainPanel, ATTR_LEFT,  &PanelLeft) ;
					GetPanelAttribute (MainPanel, ATTR_WIDTH, &PanelWidth) ;
					SetPanelAttribute (CalcN2Panel, ATTR_TOP, PanelTop+310 ) ;
					SetPanelAttribute (CalcN2Panel, ATTR_LEFT, PanelLeft+PanelWidth+4 ) ;
					DisplayPanel(CalcN2Panel) ;
				} 
				}
			if (control==PANEL_N3CALCULUS){
				GetPanelAttribute (CalcN3Panel, ATTR_VISIBLE , &PanelIsVisible);
				if (PanelIsVisible) {
				}
				else {
					GetPanelAttribute (MainPanel, ATTR_TOP,   &PanelTop) ;
					GetPanelAttribute (MainPanel, ATTR_LEFT,  &PanelLeft) ;
					GetPanelAttribute (MainPanel, ATTR_WIDTH, &PanelWidth) ;
					SetPanelAttribute (CalcN3Panel, ATTR_TOP, PanelTop+480 ) ;
					SetPanelAttribute (CalcN3Panel, ATTR_LEFT, PanelLeft+PanelWidth+4 ) ;
					DisplayPanel(CalcN3Panel) ;
				} 
				}
			break;
		}
	
	return 0;
}


int CVICALLBACK CB_OnStartNCalculus (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	
	switch (event)
	{
		case EVENT_COMMIT:
			if (panel == CalcN1Panel)
				Measuring_1 = TRUE;
			if (panel == CalcN2Panel)
				Measuring_2 = TRUE;
			if (panel == CalcN3Panel)
				Measuring_3 = TRUE;
		break;
	}
	
	return 0;
}



int CVICALLBACK CB_OnNStop (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	int PanelIsVisible;
	
	switch (event)
	{
		case EVENT_COMMIT:
			if (panel == CalcN1Panel) {
				GetPanelAttribute (CalcN1Panel, ATTR_VISIBLE, &PanelIsVisible);
				if (PanelIsVisible)
					HidePanel(CalcN1Panel);
				
				Measuring_1=FALSE;
				Step1_1=FALSE ;
				Step2_1=FALSE ;
				Step3_1=FALSE ;
				Frepminus_1=0.0;
				Frepplus_1=0.0;
			  	t1_1=0.0;
			 	t2_1=0.0;
			 	t3_1=0.0; 
				N_slope_1=0;
				Frequ_slope_1=0.0;
				Moy_slope_1=0.0;
				Slope_slope_1 =0.0;
				Ch4_slope=0.0;
				Moy_Ch4slope_1=0.0;
				Slope_Ch4slope_1=0.0;
				DDSBesChanged1=FALSE;
				DDSBesChanged2=FALSE; 
				nDDSChange=0.0;
					
				SetCtrlVal(MainPanel, PANEL_DDS1, FrequDDS1) ;
				DDS4xAD9912_SetFrequency(&DDS4xAD9912, 1, FrequDDS1);
				SetCtrlVal(MainPanel, PANEL_DDS2, FrequencyDDSBesInit) ;
				DDS4xAD9912_SetFrequency(&DDS4xAD9912, 2, FrequencyDDSBesInit); 
				nstabilization=0;
			}	
			if (panel == CalcN2Panel) {
				GetPanelAttribute (CalcN2Panel, ATTR_VISIBLE, &PanelIsVisible);
				if (PanelIsVisible)
					HidePanel(CalcN2Panel);
				
				Measuring_2=FALSE;
				Step1_2=FALSE ;
				Step2_2=FALSE ;
				Step3_2=FALSE ;
				Frepminus_2=0.0;
				Delta10K_Minus=0.0;
				Frepplus_2=0.0;
				Delta10K_Plus=0.0; 
				DeltaDDS3=0.0; 
			 	t1_2=0.0;
			  	t2_2=0.0;
			  	t3_2=0.0; 
				N_slope_2=0;
				Frequ_slope_2=0.0;
				Moy_slope_2=0.0;
				Beat_slope_2=0.0;
				Moy_Beatslope_2=0.0;
				Slope_Beatslope_2 =0.0;
				Slope_slope_2 =0.0;
				DDSBesChanged1=FALSE;
				DDSBesChanged2=FALSE; 
				nDDSChange=0.0;
				
				SetCtrlVal(MainPanel, PANEL_DDS1, FrequDDS1) ;
				DDS4xAD9912_SetFrequency(&DDS4xAD9912, 1, FrequDDS1);
				SetCtrlVal(MainPanel, PANEL_DDS2, FrequencyDDSBesInit) ;
				DDS4xAD9912_SetFrequency(&DDS4xAD9912, 2, FrequencyDDSBesInit);
				SetCtrlVal(MainPanel, PANEL_DDS3, FrequencyDDS3Init) ;
				DDS4xAD9912_SetFrequency(&DDS4xAD9912, 3, FrequencyDDS3Init); 
				nstabilization=0;
			}
			if (panel == CalcN3Panel) {		
				GetPanelAttribute (CalcN3Panel, ATTR_VISIBLE, &PanelIsVisible);
				if (PanelIsVisible)										  
					HidePanel(CalcN3Panel);
				
				Measuring_3 = FALSE;
				Frepminus_3 = 0.0;
				Frepplus_3 = 0.0;
			   	t1_3 = 0.0;
			  	t2_3 = 0.0;
			   	t3_3 = 0.0; 
				N_slope_3 = 0;
				Frequ_slope_3 = 0.0;
				Moy_slope_3 = 0.0;
				Slope_slope_3 = 0.0;
			}
		break;
	}
	return 0;
}


int  CVICALLBACK CB_OnIntegrationTime (int panel, int control, int event, 
		 void *callbackData, int eventData1, int eventData2)
{
	double Time ;
	
	switch (event) {
		case EVENT_COMMIT:
			if (panel == CalcN1Panel) {
				GetCtrlVal(CalcN1Panel, CALCN1_INTEGRATIONTIME, &Time) ;
				DeltaT_1=Time;
			}
			if (panel == CalcN2Panel) {
				GetCtrlVal(CalcN2Panel, CALCN2_INTEGRATIONTIME, &Time) ;
				DeltaT_2=Time;
			}
			if (panel == CalcN3Panel) {
				GetCtrlVal(CalcN3Panel, CALCN3_INTEGRATIONTIME, &Time) ;
				DeltaT_3=Time;
			}
			break;
	}
	return 0;
}


int  CVICALLBACK CB_OnDeltaFreq (int panel, int control, int event, 
		 void *callbackData, int eventData1, int eventData2)
{
	double Frequ ;
	
	switch (event)
	{
		case EVENT_COMMIT:
			if (panel == CalcN1Panel) {
				GetCtrlVal(CalcN1Panel, CALCN1_DELTAFREQ, &Frequ);
				DeltakHz_1 = Frequ;
			}
			if (panel == CalcN2Panel) {
				GetCtrlVal(CalcN2Panel, CALCN2_DELTAFREQ, &Frequ) ;
				DeltakHz_2=Frequ;
			}
			if (panel == CalcN3Panel) {
				GetCtrlVal(CalcN3Panel, CALCN3_DELTAFREQ, &Frequ) ;
				DeltakHz_3=Frequ;
			}
			break;
	}
	return 0;
}


int  CVICALLBACK CB_OnFindSign (int panel, int control, int event, 
		void *callbackData, int eventData1, int eventData2)
{
	
	switch (event) {
		case EVENT_COMMIT:
			switch (control)
			{
				case PANEL_FINDSIGN1:
					tbegin1 = utc;
					Frepbefore1 = Math1;
					GetCtrlVal(MainPanel, PANEL_DDS1, &Frequency1) ;
					SetCtrlVal(MainPanel, PANEL_DDS1, Frequency1+Frequencystep1) ;  
					DDS4xAD9912_SetFrequency(&DDS4xAD9912, 1, Frequency1+Frequencystep1);
					Getsign1 = TRUE;
					break;
				case PANEL_FINDSIGN2:
					tbegin2 = utc;
					Frepbefore2 = Math1;
					Ch2before = Ch2;
					GetCtrlVal(MainPanel, PANEL_DDS1, &Frequency2) ;
					SetCtrlVal(MainPanel, PANEL_DDS1, Frequency2+Frequencystep2) ;
					DDS4xAD9912_SetFrequency(&DDS4xAD9912, 1, Frequency2+Frequencystep2);
					Getsign2 = TRUE;
					break;
				case PANEL_FINDSIGN3:
					tbegin3 = utc;
					Frepbefore3 = Math1;
					GetCtrlVal(MainPanel, PANEL_DDS4, &Frequency3) ;
					SetCtrlVal(MainPanel, PANEL_DDS4, Frequency3+Frequencystep3) ; 
					DDS4xAD9912_SetFrequency(&DDS4xAD9912, 4, Frequency3+Frequencystep3);
					Getsign3 = TRUE;
					break;
			}
			break;
	}
	return 0;
}


int  CVICALLBACK CB_OnFind275K (int panel, int control, int event, 
		 void *callbackData, int eventData1, int eventData2)
 {
	 
	 
	double FrequencyToChange,DeltaCh4,DeltaCh2;
	
	DeltaCh4=275000-Ch4;
	DeltaCh2=10000-Ch2;
	
	switch (event)
		{
		case EVENT_COMMIT:
			 if (control==PANEL_FIND275K2){
					 GetCtrlVal(MainPanel, PANEL_DDS2, &FrequencyToChange) ;
					 SetCtrlVal(MainPanel, PANEL_DDS2, FrequencyToChange+DeltaCh4) ;  
					 DDS4xAD9912_SetFrequency(&DDS4xAD9912, 2, FrequencyToChange+DeltaCh4);
					 
			 }
			 
			 if (control==PANEL_FIND10K3){
					 GetCtrlVal(MainPanel, PANEL_DDS3, &FrequencyToChange) ;
					 SetCtrlVal(MainPanel, PANEL_DDS3, FrequencyToChange+DeltaCh2) ; 
					 DDS4xAD9912_SetFrequency(&DDS4xAD9912, 3, FrequencyToChange+DeltaCh2);
					 
			 }
			
		break;
			
		}
	return 0;
} 


int  CVICALLBACK CB_OnSlopeTime (int panel, int control, int event, 
		 void *callbackData, int eventData1, int eventData2)
{
	double Time ;
	
	switch (event)
	{
		case EVENT_COMMIT:
			if (panel == CalcN1Panel) {
				GetCtrlVal(CalcN1Panel, CALCN1_SLOPETIME, &Time);
				SlopeTime1=Time;
			}
			if (panel == CalcN2Panel) {
				GetCtrlVal(CalcN2Panel, CALCN2_SLOPETIME, &Time);
				SlopeTime2=Time;
			}
			if (panel == CalcN3Panel) {
				GetCtrlVal(CalcN3Panel, CALCN3_SLOPETIME, &Time);
				SlopeTime3=Time;
			}
			break;
	}
	return 0;
}


 int  CVICALLBACK CB_OnChangeNdiv (int panel, int control, int event, 
		 void *callbackData, int eventData1, int eventData2)
{
	double NewN ;
	
	switch (event)
		{
		case EVENT_COMMIT:
			   if (control==PANEL_CHANGENDIV){ 
				GetCtrlVal(MainPanel, PANEL_CHANGENDIV, &NewN) ;
				Ndiv=NewN;
				FrequDDS1=880000000.0/Ndiv;
				DeltakHz_1=1800.0*8/Ndiv;
				DeltakHz_2=1500.0*8/Ndiv;
				SetCtrlVal(MainPanel, PANEL_DDS1, FrequDDS1) ;
				DDS4xAD9912_SetFrequency(&DDS4xAD9912,1,FrequDDS1);
				SetCtrlVal(CalcN1Panel, CALCN1_DELTAFREQ, DeltakHz_1) ;
				SetCtrlVal(CalcN2Panel, CALCN2_DELTAFREQ, DeltakHz_2) ; 
			   }

			break;
			
		}
	return 0;
}


 int  CVICALLBACK CB_OnStartSlopeCancelling (int panel, int control, int event, 
		 void *callbackData, int eventData1, int eventData2)
{

	double   ActualFrequ=0.0;
	double   ActualSlope=0.0; 
	int buttonstatu=0;
	
	switch (event)
		{
		case EVENT_COMMIT:
			   if (control==PANEL_STARTCANCEL){ 
				   
				   GetCtrlVal(MainPanel, PANEL_STARTCANCEL, &buttonstatu) ; 
				   
				    if (buttonstatu==1)
					{
						SlopeMeasuringTimeBegin=utc; 
						SlopeMeasuring=TRUE;
						ActualFrequ=DDSFox_ReadFreq(&DDS1xAD9956);
						ActualSlope= SlopeMath2;
						DDSFox_Set(&DDS1xAD9956, ActualFrequ, ActualSlope);
						Nratio=-1;
					 	CenterFrequencyCh2ToDetermine=TRUE ;
						CenterFrequencyCh2=0.0; 

					}
					else
					{
						if (KeepFrequ==TRUE) {ActualFrequ=DDSFox_ReadFreq(&DDS1xAD9956);}
						else {  ActualFrequ= DEDRIFT_DDS_FREQUENCY;}
						if (KeepSlope==TRUE) {ActualSlope= SlopeMath2;
											  OnSlopeCancelling =TRUE; }
						else {  ActualSlope= 0.0;
								OnSlopeCancelling =FALSE;     }
						 DDSFox_Set(&DDS1xAD9956, ActualFrequ, ActualSlope);
						 SlopeMeasuring=FALSE; 
						 N_Math2slope=0.0;
						 Math2_slope=0.0;
						 Moy_Math2slope=0.0;
						 Slope_Math2slope=0.0;
						 MoyMath2=0.0; 
						 Slope_Math2slope=0.0;
						 CenterFrequencyCh2= 0.0;
						 Nratio=-1;
						 nstabilisationSlopeMeasuring=0;
					}
			   }

			break;
			
		}
	return 0;
}

 
int  CVICALLBACK CB_OnResetSlope (int panel, int control, int event, 
		void *callbackData, int eventData1, int eventData2)
{

	switch (event) {
		case EVENT_COMMIT:
			SlopeMath2 = 0.0; 
			CenterFrequencyCh2 = 0.0;
			MoyMath2 = 0.0;
			Nratio = -1;
			CenterFrequencyCh2ToDetermine = TRUE ;
			OnSlopeCancelling = FALSE; 
			DDSFox_Set(&DDS1xAD9956, DEDRIFT_DDS_FREQUENCY, 0.0);
			break;
	}
	return 0;
}
  
  
int  CVICALLBACK CB_ChangeSlopeTime (int panel, int control, int event, 
		void *callbackData, int eventData1, int eventData2)
{

	switch (event) {
		case EVENT_COMMIT:
			GetCtrlVal(MainPanel, PANEL_SLOPETIME, &TimetoSlope); 
			break;
	}
	return 0;
}


int  CVICALLBACK CB_OnCROX (int panel, int control, int event, 
		 void *callbackData, int eventData1, int eventData2)
{
		bool BoxChecked = FALSE; 
		double ActualFrequ=0.0;

	switch (event)
		{
		case EVENT_COMMIT:
			
			if (control== PANEL_CHECKBOX_CORRFREQU)  
			{
				   	GetCtrlVal(MainPanel, PANEL_CHECKBOX_CORRFREQU, &BoxChecked) ;
						if (BoxChecked) {
							FrequCorrec=TRUE ;
							}
						else {FrequCorrec=FALSE ;}
							
			}
			
			if (control== PANEL_CHECKBOX_KEEP)  
			{
				   	GetCtrlVal(MainPanel, PANEL_CHECKBOX_KEEP, &BoxChecked) ;
						if (BoxChecked) {
							KeepFrequ=TRUE ;
							}
						else {KeepFrequ=FALSE ;}
							
			}
			
			
			if (control== PANEL_CHECKBOX_KEEPSLOPE)  
			{
				   	GetCtrlVal(MainPanel, PANEL_CHECKBOX_KEEPSLOPE, &BoxChecked) ;
						if (BoxChecked) {
							KeepSlope=TRUE ;
							if (SlopeMeasuring==FALSE && OnSlopeCancelling==FALSE)
							  {
								    ActualFrequ = DDSFox_ReadFreq(&DDS1xAD9956);
									DDSFox_Set(&DDS1xAD9956, ActualFrequ, SlopeMath2);
							 		OnSlopeCancelling=TRUE;
								  
								  
							  }
							}
						else {KeepSlope=FALSE ;
							  if (SlopeMeasuring==FALSE && OnSlopeCancelling==TRUE)
							  {
								  	  ActualFrequ=DDSFox_ReadFreq(&DDS1xAD9956);
									  DDSFox_Set(&DDS1xAD9956, ActualFrequ, 0.0);
							 		  OnSlopeCancelling=FALSE;
							  }
						
						}
							
			}
			
 			break;
			
		}
	return 0;
}



int  CVICALLBACK CB_OnReCentering (int panel, int control, int event, 
		 void *callbackData, int eventData1, int eventData2)
{
		bool BoxChecked = FALSE; 

	switch (event)
		{
		case EVENT_COMMIT:
			
			if (control== PANEL_CHECKBOX_RECENTER)  
			{
				   	GetCtrlVal(MainPanel, PANEL_CHECKBOX_RECENTER, &BoxChecked) ;
						if (BoxChecked) {
							KeepCentering=TRUE ;
							CenteringTimeBegin275K=utc;
							CenteringTimeBegin10K=utc;
							}
						else {KeepCentering=FALSE ;}
							
			}
			
			
 			break;
			
		}
	return 0;
}


int  CVICALLBACK CB_OnAutoStopSlopeCancellingIfDelock (int panel, int control, int event, 
		 void *callbackData, int eventData1, int eventData2)
{
	switch (event)
	{
		case EVENT_COMMIT:
			GetCtrlVal(MainPanel, PANEL_CHECKBOX_STOPIFAUTODE, &AutoStopSlopeCancellingIfDelock);
			break;
	}
	return 0;
}

int CVICALLBACK CB_OnSlopeReference (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	switch (event)
	{
		case EVENT_COMMIT:
			GetCtrlVal(MainPanel, PANEL_SLOPE_REFERENCE, &slopeReference);
			break;
	}
	return 0;
}


int CVICALLBACK CB_OnEstimateN (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	int visible;
	double wl;
	
	switch (event)
	{
		case EVENT_COMMIT:
			/* be prepaid to support more N estimates for different beat notes */
			switch (control)
			{
				case PANEL_ESTIMATE_N3:
					GetPanelAttribute(EstimateN3Panel, ATTR_VISIBLE , &visible);
					if (! visible) {
						DisplayPanel(EstimateN3Panel);
					}
					/* set current frep */
					SetCtrlVal(EstimateN3Panel, ESTIMATEN3_FREP, 250e6 + Math1);
					/* default wavelenght for Sr cavity */
					GetCtrlVal(EstimateN3Panel, ESTIMATEN3_WAVELENGTH, &wl);
					if (wl == 0.0)
						SetCtrlVal(EstimateN3Panel, ESTIMATEN3_WAVELENGTH, 698.50);
					/* reset N3 estimate */
					SetCtrlVal(EstimateN3Panel, ESTIMATEN3_N, 0.0);
					break;
			}
			break;
	}
	
	return 0;
}			

int CVICALLBACK CB_OnNEstimateCancel (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	int visible;
	
	switch (event)
	{
		case EVENT_COMMIT:
			GetPanelAttribute(panel, ATTR_VISIBLE, &visible);
			if (visible)
				HidePanel(panel);
			break;
	}
	return 0;
}

int CVICALLBACK CB_OnNEstimateSet (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	switch (event)
	{
		case EVENT_COMMIT:
			GetCtrlVal(panel, ESTIMATEN3_N, &N3);
			SetCtrlVal(MainPanel, PANEL_N3CHOICE, N3);
			HidePanel(panel);
			break;
	}
	return 0;
}

int CVICALLBACK CB_OnNEstimate (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	double frep, wl, N;
	
	switch (event)
	{
		case EVENT_COMMIT:
			GetCtrlVal(panel, ESTIMATEN3_FREP, &frep);
			GetCtrlVal(panel, ESTIMATEN3_WAVELENGTH, &wl);
			N = round(299792458.0 / wl / 1e-9 / frep);
			SetCtrlVal(panel, ESTIMATEN3_N, N);
			break;
	}
	return 0;
}

int CVICALLBACK CB_SetSlope (int panel, int control, int event,
		void *callbackData, int eventData1, int eventData2)
{
	double frequency, slope;
	
	switch (event)
	{
		case EVENT_COMMIT:
			GetCtrlVal(panel, control, &slope);
			frequency = DDSFox_ReadFreq(&DDS1xAD9956);
			DDSFox_Set(&DDS1xAD9956, frequency, slope);
			break;
	}
	return 0;
}