void getdatapointfromgraph(string graphname, double& max, int& indexofpoint)
 {
 	GraphPage pg;
 	pg = Project.GraphPages(graphname);
 	
 	
 	if( pg.IsValid() )
	{
			string setcmd;
 			GraphLayer		gl = pg.Layers(); // Get active layer
			printf("%s has %d layer(s)\n", pg.GetName(), pg.Layers.Count());
 			setcmd.Format("win -a %s", graphname);
 			LT_execute(setcmd); //Set graph active
 			setcmd.Format("win -z %s", graphname);
 			LT_execute(setcmd); //Maximize Graph
 			//DataPlot dp = gl.DataPlots(0);
 			//dp.SetActive(0,20);

 			
 			int nPts=1;
 			// Declare vectors to hold coordinates
			vector vecX, vecY, vecIndex;
			
			// Call function to get points
			int iRet = get_screen_reader_points(nPts, vecX, vecY, vecIndex, graphname);
			if( iRet < 0 )
				printf("Function call to get screen points returned error: %d\n", iRet);

			printf("User clicked on %d points\n", iRet);
		
			if( 0 != iRet )
			{
				printf("X, Y coordinates of points are:\n");
				for(int ii = 0; ii < iRet; ii++)
					{
					indexofpoint = vecIndex[ii];
					max=vecY[ii];
					printf(" %d: X = %f, Y = %f Index= %d\n", ii, vecX[ii], vecY[ii], indexofpoint); 
					}
			} 
			
			setcmd.Format("win -ia %s", graphname);
			LT_execute(setcmd); //Minimize Graph and activate previous window
 			
 				
 				
	}
	else 
	{
		printf("Invalid Graphwindow (%s) for selecting Maximum ", graphname)
	}
	

 } 

static int get_screen_reader_points(int nPts, vector& vecX, vector& vecY, vector& vecIndex, string graphname)
{
	// Return error if number of points not appropriate
	// Arbitrary upper limit of 99 - change as desired
	if( nPts < 1 || nPts > 99) return -1;

	// Run the LabTalk section in this file to give user control to click points
	string strThisFile = __FILE__;
	string strCMD;
	strCMD.Format("run.section(%s,GetScreenPts,%d,%s)", strThisFile, nPts, graphname);
	LT_execute(strCMD);

	// On return, the LabTalk variable tmp_scr_count has number of points user clicked + 1
	double dCount;
	LT_get_var("tmp_scr_count", &dCount);
	int nCount = (int) (dCount - 1);
	if( 0 == nCount ) return 0;

	// Get values of x, y coordinates from LabTalk variables and delete the vars
	vecX.SetSize(nCount);
	vecY.SetSize(nCount);
	vecIndex.SetSize(nCount);
	string strX, strY, strIndex;
	double dX, dY, dIndex;
	for(int ii = 1; ii <= nCount; ii++)
	{
		strX.Format("tmp_scr_pos_x%d", ii);
		strY.Format("tmp_scr_pos_y%d", ii);
		strIndex.Format("tmp_scr_pos_index%d", ii);
		LT_get_var(strX, &vecX[ii - 1]);
		LT_get_var(strY, &vecY[ii - 1]);
		LT_get_var(strIndex, &vecIndex[ii - 1]);
		strCMD.Format("del -v %s", strX);
		LT_execute(strCMD);
		strCMD.Format("del -v %s", strY);
		LT_execute(strCMD);
		strCMD.Format("del -v %s", strIndex);
		LT_execute(strCMD);
	}

	// Delete the count LabTalk variable
	LT_execute("del -v tmp_scr_count;");

	// Return the number of points clicked by user
	return nCount;
} 



// This section of LabTalk will be ignored by compiler since it is contained within an ifdef block.
// The section of script within this block can then be successfully called by a run.section LabTalk command.
#ifdef LABTALK 
[GetScreenPts]
// Number of points to get is passed as 1st parameter
// Select screen reader tool
dotool 3;  // 2 screen Reader, 3 data reader

// Main loop to accept points
tmp_scr_count = 1;
def quittoolbox
{
	// Check if user has clicked required number of points
	if( 1 == tmp_scr_count )
		done = 1;
	if( %1 != tmp_scr_count )
		return;
}


// Get points and store in temp string variables
def pointproc
{
	tmp_scr_pos_x$(tmp_scr_count)=x;
	tmp_scr_pos_y$(tmp_scr_count)=y;
	tmp_scr_pos_index$(tmp_scr_count)=index;
	if( tmp_scr_count >= %1 )
	{
		dotool 0;
		// Set the variable to break the infinite loop
		done = 1;
	}
	tmp_scr_count += 1;
}

// Set up an "infinite loop" so code waits for user to click points
for( done = 0; done == 0; )
{
// wait a while
sec -p 0.1;
} 
#endif 

void testdataget()
{
	//Create wks
	Worksheet wks1;
	wks1.Create();
	wks1.GetPage().Rename("wks1");
	
	
	//Fill col1
	Dataset ds1(wks1,0);
	wks1.Columns(0).SetName("col1");
	ds1.SetSize(20);
	
	for(int ii=0; ii<20; ii++)
	{
		ds1[ii] = ii;
	}

	//Fill col2
	Dataset ds2(wks1,1);
	wks1.Columns(1).SetName("col2");
	ds2.SetSize(20);
	for(ii=0; ii<20; ii++)
	{
		ds2[ii] = ii;
	}
	
	
	string strTemplate = "origin.otp" ;
	
	
	GraphPage grphData;
	bool thegrph = grphData.Create(strTemplate);
	GraphLayer grphLayer = grphData.Layers();
	Curve crvData(wks1, 0, 1);
	int	nPlot = grphLayer.AddPlot(crvData, IDM_PLOT_LINE);
	grphLayer.Rescale(); 
	grphData.Rename("Graph1");
	
	LT_execute("win -a Graph1");
	LT_execute("win -z Graph1");
	read_screen_points();
	

}







void read_screen_points(int nPts = 1)
{
// Declare vectors to hold coordinates
vector vecX, vecY;
// Call function to get points
int iRet = get_screen_reader_points(nPts, vecX, vecY);
if( iRet < 0 )
printf("Function call to get screen points returned error: %d\n", iRet);

printf("User clicked on %d points\n", iRet);
if( 0 != iRet )
{
printf("X, Y coordinates of points are:\n");
for(int ii = 0; ii < iRet; ii++)
printf(" %d: X = %f, Y = %f\n", ii, vecX[ii], vecY[ii]);
}
}




// This static function handles calling LabTalk to get user to click on the screen.
// Parameters:
// nPts: no. of points that user should click
// vecX: vector to hold x coordinates of points, passed by reference
// vecX: vector to hold y coordinates of points, passed by reference
// Return:
// -1: no. of points inappropriate
// >=0: no. of points user actually clicked
static int get_screen_reader_points(int nPts, vector& vecX, vector& vecY)
{
// Return error if number of points not appropriate
// Arbitrary upper limit of 99 - change as desired
if( nPts < 1 || nPts > 99) return -1;

// Run the LabTalk section in this file to give user control to click points
string strThisFile = __FILE__;
string strCMD;
strCMD.Format("run.section(%s,GetScreenPts,%d)", strThisFile, nPts);
LT_execute(strCMD);

// On return, the LabTalk variable tmp_scr_count has number of points user clicked + 1
double dCount;
LT_get_var("tmp_scr_count", &dCount);
int nCount = (int) (dCount - 1);
if( 0 == nCount ) return 0;

// Get values of x, y coordinates from LabTalk variables and delete the vars
vecX.SetSize(nCount);
vecY.SetSize(nCount);
string strX, strY;
double dX, dY;
for(int ii = 1; ii <= nCount; ii++)
{
strX.Format("tmp_scr_pos_x%d", ii);
strY.Format("tmp_scr_pos_y%d", ii);
LT_get_var(strX, &vecX[ii - 1]);
LT_get_var(strY, &vecY[ii - 1]);
strCMD.Format("del -v %s", strX);
LT_execute(strCMD);
strCMD.Format("del -v %s", strY);
LT_execute(strCMD);
}

// Delete the count LabTalk variable
LT_execute("del -v tmp_scr_count;");

// Return the number of points clicked by user
return nCount;
}

// This section of LabTalk will be ignored by compiler since it is contained within an ifdef block.
// The section of script within this block can then be successfully called by a run.section LabTalk command.
#ifdef LABTALK
[GetScreenPts]
// Number of points to get is passed as 1st parameter
// Select screen reader tool
dotool 3;

// Main loop to accept points
tmp_scr_count = 1;
def quittoolbox
{
// Check if user has clicked required number of points
if( 1 == tmp_scr_count )
done = 1;
if( %1 != tmp_scr_count )
return;
}

// Get points and store in temp string variables
def pointproc
{
tmp_scr_pos_x$(tmp_scr_count)=x;
tmp_scr_pos_y$(tmp_scr_count)=y;
if( tmp_scr_count >= %1 )
{
dotool 0;
// Set the variable to break the infinite loop
done = 1;
}
tmp_scr_count += 1;
}
// Set up an "infinite loop" so code waits for user to click points
for( done = 0; done == 0; )
{
// wait a while
sec -p 0.1;
}
#endif 

 

string strThisFile = __FILE__;
	string strCMD;
	strCMD.Format("run.section(%s,GetScreenPts,%d)", strThisFile, nPts); 
	LT_execute(strCMD);

 

 
if( MessageBox( GetWindow(), "Manual Selection of Maxima?", "Attention!" ,
	MB_YESNO | MB_ICONQUESTION | MB_DEFBUTTON1 ) == IDYES )
	{
		setmanual=1;
	}
	
	if (setmanual==1)
	{
		double maxN=0;
		int indexN=0;
		double maxS=0;
		int indexS=0;
		
		//Korrfkt Nord
		LT_execute("win -a grph_Ausl_AbsTheta_Nord");
		LT_execute("win -z grph_Ausl_AbsTheta_Nord");
		getdatapointfromgraph("grph_Ausl_AbsTheta_Nord", maxN, indexN);
		printf("Maximum N Y=%f Index=%d\n", maxN, indexN);
		
		
		dsXRMRKorrfktN_AbsTheta=dsXRMRnormAuslN_AbsTheta/maxN;
		
		int iRowN=indexN;
		int sizeN=dsXRMRKorrfktN_AbsTheta.GetSize();
		for (int ii = iRowN+1; ii < sizeN ; ii++)
		{
			dsXRMRKorrfktN_AbsTheta[ii] = 1.0;
		}
		
		// Korrfkt Sued
		LT_execute("win -a grph_Ausl_AbsTheta_Sued");
		LT_execute("win -z grph_Ausl_AbsTheta_Sued");
		getdatapointfromgraph("grph_Ausl_AbsTheta_Sued", maxS, indexS);
		printf("Maximum S Y=%f Index=%d\n", maxS, indexS);
		
		
		dsXRMRKorrfktS_AbsTheta=dsXRMRnormAuslS_AbsTheta/maxS;
		
		int iRowS=indexS;
		int sizeS=dsXRMRKorrfktS_AbsTheta.GetSize();
		for (ii = iRowS+1; ii < sizeS ; ii++)
		{
			dsXRMRKorrfktS_AbsTheta[ii] = 1.0;
		}
		
		//Doubleplot of the functions used to normalize
		colx = wks2.Columns("AbsTheta");                 
		coly1 = wks2.Columns("Korrfkt.N"); 
		coly2 = wks2.Columns("Korrfkt.S");
		strColNameX = colx.GetName();                  
		strColNameY1 = coly1.GetName();  
		strColNameY2 = coly2.GetName();  
		GraphTemplate = "Ausl";
		GraphName = "Ausl_SAbs_KorrFkt";
		xColNum = colx.GetIndex();
		y1ColNum = coly1.GetIndex();
		y2ColNum = coly2.GetIndex();
		CreateDoublePlot(wks2, xColNum, y1ColNum, y2ColNum, GraphName, GraphTemplate, strPathToTemplate);
	}
	else
	{
		dsXRMRKorrfktN_AbsTheta=dsXRMRnormAuslN_AbsTheta;
		dsXRMRKorrfktS_AbsTheta=dsXRMRnormAuslS_AbsTheta;
	}

 
if( MessageBox( GetWindow(), "Manual Selection of Maxima?", "Attention!" ,
	MB_YESNO | MB_ICONQUESTION | MB_DEFBUTTON1 ) == IDYES )

 
LT_execute("win -a grph_Ausl_AbsTheta_Nord");
		LT_execute("win -z grph_Ausl_AbsTheta_Nord");
		getdatapointfromgraph("grph_Ausl_AbsTheta_Nord", maxN, indexN);
The window appears maximized and accepts the selection. The next call 
LT_execute("win -a grph_Ausl_AbsTheta_Sued");
		LT_execute("win -z grph_Ausl_AbsTheta_Sued");
		getdatapointfromgraph("grph_Ausl_AbsTheta_Sued", maxS, indexS);