WebSVN – f9daq – Blame – /cvi/instr/DRS/drs4.c

Rev	Author	Line No.	Line
217	f9daq	1	#include "H1D.h"
		2	#include "H2D.h"
203	f9daq	3	#include <formatio.h>
195	f9daq	4	#include <utility.h>
		5	#include <ansi_c.h>
203	f9daq	6	#include <cvirte.h>
195	f9daq	7	#include <userint.h>
		8	#include "drs4.h"
203	f9daq	9	#include "drsread.h"
195	f9daq	10
203	f9daq	11	#define MIKRO
		12
		13	#ifdef MIKRO
		14	#include "MIKRO.h"
		15	#endif
		16
195	f9daq	17	static int daq_on;
203	f9daq	18	static int ph, p2,plothandle[4]= {0,0,0,0};
217	f9daq	19	static int scanplothandle;
195	f9daq	20	static int tfID;
		21	static int controlID;
		22
203	f9daq	23	#define MAX_THREADS 10
195	f9daq	24
203	f9daq	25	static CmtThreadPoolHandle poolHandle = 0;
195	f9daq	26
217	f9daq	27	static float gSum[4]={0,0,0,0};
195	f9daq	28	int main (int argc, char *argv[]) {
203	f9daq	29	short port;
		30	if (InitCVIRTE (0, argv, 0) == 0)
		31	return -1; /* out of memory */
		32	if ((ph = LoadPanel (0, "drs4.uir", PANEL)) < 0)
		33	return -1;
		34	if ((p2 = LoadPanel (0, "drs4.uir", SCAN)) < 0)
		35	return -1;
		36	SetStdioPort (CVI_STDIO_WINDOW);
		37	SetSleepPolicy(VAL_SLEEP_MORE);
		38	CmtNewThreadPool (MAX_THREADS, &poolHandle);
		39
		40	DisplayPanel (ph);
		41	DisplayPanel (p2);
		42
		43	#ifdef MIKRO
		44	GetCtrlVal(p2, SCAN_PORT, &port);
		45	if (MIKRO_Open (port)) MessagePopup ("Error", "Mikro Port Not found !\n Change in the GUI") ;
		46	MIKRO_Init(1,0);
		47	MIKRO_Init(2,0);
		48	#endif
		49
		50	RunUserInterface ();
		51	DiscardPanel (ph);
		52	DiscardPanel (p2);
		53	CmtDiscardThreadPool (poolHandle);
		54	MIKRO_Close ();
		55	return 0;
195	f9daq	56	}
		57
		58
210	f9daq	59	char strbuf[0xFF];
195	f9daq	60
210	f9daq	61	int gLog=0;
195	f9daq	62
210	f9daq	63	int printf(const char *format, ...) {
		64	va_list aptr;
		65	int ret;
		66	FILE *flog;
195	f9daq	67
210	f9daq	68	va_start(aptr, format);
		69	ret = vsprintf(strbuf, format, aptr);
		70	va_end(aptr);
		71	SetCtrlVal(ph,PANEL_STDIO,strbuf);
		72
		73	if (gLog) {
		74	flog = fopen ("stdio.log", "a");
		75	fprintf (flog, "%s", strbuf);
		76	fclose (flog);
		77	}
		78	return(ret);
		79	}
		80
		81
		82
195	f9daq	83	static void start_timer (double tout) {
203	f9daq	84	SetCtrlAttribute (ph, PANEL_TIMER, ATTR_INTERVAL, tout);
		85	SetCtrlAttribute (ph, PANEL_TIMER, ATTR_ENABLED, 1);
195	f9daq	86	}
		87
		88	static void stop_timer ( void ) {
203	f9daq	89	SetCtrlAttribute (ph, PANEL_TIMER, ATTR_ENABLED, 0);
		90	DRSSetTimeout();
195	f9daq	91	}
		92
		93
		94
		95
		96	void CVICALLBACK EndOfThread ( CmtThreadPoolHandle poolhandle,
203	f9daq	97	CmtThreadFunctionID functionID, unsigned int event,
		98	int value, void *callbackData ) {
195	f9daq	99
203	f9daq	100	daq_on=0;
		101	//SetDimming(0);
		102	printf("End of Thread \n");
		103	return ;
195	f9daq	104
		105	}
		106
		107
217	f9daq	108	static float xs[4][1024];
203	f9daq	109	int CVICALLBACK daq(void *functionData) {
195	f9daq	110
		111
203	f9daq	112	int neve;
		113	char filename[0xff];
		114	int imask[4];
		115	unsigned long mask;
		116	int frequency;
		117	double trgdelay;
		118	double trglevel;
		119	int trgtype;
		120	int trgchannel;
		121	int trgpolarity;
		122	int verbose;
		123	double range;
		124	int pfreq;
		125	int enabledoutput;
		126	int neveold = 0;
		127	double rate;
		128	int args = (int ) functionData;
217	f9daq	129	double twin[2];
201	f9daq	130
217	f9daq	131
203	f9daq	132	GetCtrlVal(ph, PANEL_CH0, &imask[0] );
		133	GetCtrlVal(ph, PANEL_CH1, &imask[1] );
		134	GetCtrlVal(ph, PANEL_CH2, &imask[2] );
		135	GetCtrlVal(ph, PANEL_CH3, &imask[3] );
		136	mask = 0;
		137	for (int i=0; i<4; i++) {
		138	if (imask[i]) mask \|= (1<<i);
201	f9daq	139	}
203	f9daq	140
		141	GetCtrlVal(ph,PANEL_NEVE, &neve);
		142	GetCtrlVal(ph,PANEL_DEBUG, &verbose);
		143	GetCtrlVal(ph,PANEL_PFREQ, &pfreq);
		144	GetCtrlVal(ph,PANEL_ENABLEDOUTPUT, &enabledoutput);
		145
		146	GetCtrlVal(ph,PANEL_FREQUENCY, &frequency);
		147	GetCtrlVal(ph,PANEL_TRGDELAY, &trgdelay);
		148	GetCtrlVal(ph,PANEL_TRGCHANNEL, &trgchannel);
		149
		150	GetCtrlVal(ph,PANEL_TRGTYPE, &trgtype);
		151	GetCtrlVal(ph,PANEL_TRGLEVEL, &trglevel);
		152	GetCtrlVal(ph,PANEL_TRGPOLARITY, &trgpolarity);
		153	GetCtrlVal(ph,PANEL_RANGE, &range);
		154
217	f9daq	155	GetCtrlVal(ph,PANEL_TWIN0, &twin[0]);
		156	GetCtrlVal(ph,PANEL_TWIN1, &twin[1]);
203	f9daq	157	//printf("mask=0x%x\n",mask);
		158
		159	DRSSetMask( (unsigned char)( mask & 0xF ) );
		160
		161	DRSSetFrequency( frequency );
		162	DRSSetTriggerDelay(trgdelay );
		163	DRSSetTriggerChannel(trgchannel );
		164	DRSSetTriggerType( trgtype );
		165	DRSSetTriggerLevel(trglevel);
		166	DRSSetTriggerPolarity(trgpolarity);
		167
		168	DRSSetRange ( range );
		169
		170
		171
		172	FILE *fp= NULL;
		173
		174	if (enabledoutput) {
217	f9daq	175	if (args[0]) {
		176	GetCtrlVal(ph, PANEL_FILENAME, filename );
		177	fp = fopen(filename,"wb");
		178	} else {
		179	GetCtrlVal(p2, SCAN_FILENAME, filename );
		180	fp = fopen(filename,"ab");
		181	}
203	f9daq	182	}
		183
		184	static unsigned char *buffer;
		185
		186	int buffer_size = 0;
		187	const int nBoards=1;
		188	const int waveDepth=1024;
207	f9daq	189
203	f9daq	190	if (buffer_size == 0) {
207	f9daq	191	buffer_size = 8; // file header + time header
		192	buffer_size += nBoards * (4 + 4(4+waveDepth4)); // bin widths
		193	buffer_size += 24 + nBoards * (8 + 4(8+waveDepth2));
		194	buffer = (unsigned char *)malloc(buffer_size);
		195	}
		196
203	f9daq	197	time_t t=0,told=0, tstart=0;
		198
		199	if (!DRSInit()) {
		200	time(&tstart);
		201	told=tstart;
217	f9daq	202	int nev=0;
		203	for (int k = 0;k<4;k++){
		204	gSum[k]=0;
		205	}
		206	for (int i=0; i<neve; i++) {
203	f9daq	207	start_timer(1);// 1 s timeout
		208	int retval = DRSRead(0);
		209	stop_timer();
		210	int nb = ( retval == 0 && fp ) ? DRSToBuffer( buffer , i ) : 0;
		211	SetCtrlVal(ph,PANEL_CEVE,i);
		212	if (retval) i--;
217	f9daq	213	nev++;
203	f9daq	214	if (!daq_on) break;
		215	time(&t);
		216	if (t!=told ) {
		217	rate = (i-neveold);
		218	printf("%d events in %2.2f min (%d s) Rate %f Hz %s ",i+1, (double)(t-tstart)/60.,(t-tstart), rate , ctime(&t));
		219	GetCtrlVal(ph,PANEL_PFREQ, &pfreq);
		220	neveold = i;
		221	}
		222	told=t;
		223	// Save data
		224	if (nb>0 && fp) fwrite(buffer, 1,nb ,fp);
		225	// Plot Data
217	f9daq	226	for (int k=0; k<4; k++) {
		227	if ( (mask & ( 0x1<<k )) ){
		228	float *t=DRSGetTime(k);
		229	float *x=DRSGetWave(k);
		230	;
		231	for (int j=0 ; j<1024 ; j++) {
		232	xs[k][j]= x[j]*1e-3;
		233
		234	if (t[j]> twin[0] && t[j] < twin[1]) gSum[k]+= fabs(xs[k][j]);
		235	if (verbose) printf("[%d] %d. x= %3.2f y=%3.2f\n", k, i, t[j], x[j] );
		236	//h[k]->Fill( t[i], x[i]*1e-3);
		237	}
203	f9daq	238
217	f9daq	239	if (i % pfreq == 0) {
		240
		241	const int col[4]= {VAL_WHITE,VAL_RED,VAL_GREEN,VAL_BLUE};
		242	if (plothandle[k]) DeleteGraphPlot (ph, PANEL_GRAPH, plothandle[k], VAL_IMMEDIATE_DRAW);
		243	plothandle[k] = PlotXY (ph, PANEL_GRAPH, t, xs[k], 1024, VAL_FLOAT, VAL_FLOAT, VAL_THIN_LINE, VAL_NO_POINT, VAL_SOLID, 1, col[k]);
		244	}
203	f9daq	245	}
217	f9daq	246
203	f9daq	247	}
		248
		249
		250	}
		251	time(&t);
217	f9daq	252	printf("%d events in %2.2f min (%d s) %s",nev, (double)(t-tstart)/60.,t-tstart, ctime(&t));
203	f9daq	253	DRSEnd();
		254	}
		255
		256	if (fp) fclose(fp);
		257
		258	free(buffer);
		259
		260	return 0;
		261
195	f9daq	262	}
		263
		264
		265
203	f9daq	266	int CVICALLBACK scan(void *functionData) {
195	f9daq	267
207	f9daq	268	int dx[3]={0,0,0};
		269	int nx[3]={0,0,0};
		270	int x0[3]={0,0,0};
		271	int ix[3]={0,0,0};
		272	int idx[3]={0,0,0};
203	f9daq	273	int size;
207	f9daq	274	char posrec[4]="POSR";
		275	char runbuf[4]="PRUN";
203	f9daq	276
207	f9daq	277	int n[3];
203	f9daq	278	char filename[0xFF];
		279	int enabledoutput;
		280
		281	FILE *fp;
217	f9daq	282	GetCtrlVal(p2, SCAN_FILENAME, filename );
203	f9daq	283
		284	GetCtrlVal(ph,PANEL_ENABLEDOUTPUT, &enabledoutput);
		285
		286	if ( GetFileInfo(filename,&size) ) {
217	f9daq	287	MessagePopup ("Warning","File exist. Remove it first or choose another file");
		288	return 0;
203	f9daq	289	}
		290	GetCtrlVal(p2, SCAN_STEPX, &dx[0]);
		291	GetCtrlVal(p2, SCAN_STEPY, &dx[1]);
		292	GetCtrlVal(p2, SCAN_NSTEPSX, &nx[0]);
		293	GetCtrlVal(p2, SCAN_NSTEPSY, &nx[1]);
		294	GetCtrlVal(p2, SCAN_STARTX, &x0[0]);
		295	GetCtrlVal(p2, SCAN_STARTY, &x0[1]);
217	f9daq	296	for (int k=0;k<4;k++) {
		297	H2D_Init(k, "charge","Induced charge",
		298	nx[0], x0[0] - dx[0]0.5 ,x0[0] + dx[0] ( nx[0] - 0.5) ,
		299	nx[1], x0[1] - dx[1]0.5 ,x0[1] + dx[1] ( nx[1] - 0.5));
		300	H1D_Init(k, "charge","Induced charge projection x",
		301	nx[0], x0[0] - dx[0]0.5 ,x0[0] + dx[0] ( nx[0] - 0.5) );
		302	H1D_Init(100 + k, "charge","Induced charge projection y",
		303	nx[1], x0[1] - dx[1]0.5 ,x0[1] + dx[1] ( nx[1] - 0.5) );
		304	}
207	f9daq	305	if (enabledoutput) {
		306	fp = fopen(filename,"ab");
		307	if (fp) {
		308	size=36;
		309	fwrite(runbuf, 1,4 ,fp);
		310	fwrite(&size , 1,4 ,fp);
		311	fwrite(x0 , 1,4*3 ,fp);
		312	fwrite(dx , 1,4*3 ,fp);
		313	fwrite(nx , 1,4*3 ,fp);
		314	fclose(fp);
		315	}
		316	}
		317
203	f9daq	318	for (int i=0; i<nx[0]; i++) {
		319
		320	ix[0]= x0[0]+i*dx[0];
		321	#ifdef MIKRO
		322	MIKRO_MoveTo(1,ix[0]);
		323	#endif
		324	SetCtrlVal (p2, SCAN_IX, i);
		325	for (int j=0; j<nx[1]; j++) {
		326
		327	SetCtrlVal (p2, SCAN_IY, j);
		328
		329	ix[1]= x0[1]+j*dx[1];
		330	#ifdef MIKRO
		331	MIKRO_MoveTo(2,ix[1]);
		332
		333	MIKRO_GetPosition(1,&n[0]);
		334	SetCtrlVal (p2, SCAN_XP, n[0]);
		335	MIKRO_GetPosition(2,&n[1]);
		336	SetCtrlVal (p2, SCAN_YP, n[1]);
207	f9daq	337
		338	if (enabledoutput) {
203	f9daq	339	fp = fopen(filename,"ab");
		340	if (fp) {
207	f9daq	341	idx[0]=i;
		342	idx[1]=j;
		343	size=24;
		344	fwrite(posrec, 1,4 ,fp);
		345	fwrite(&size , 1,4 ,fp);
		346	fwrite(n , 1,4*3 ,fp);
		347	fwrite(idx , 1,4*3 ,fp);
203	f9daq	348	fclose(fp);
		349	}
		350	}
		351
		352	#endif
		353	daq(functionData);
217	f9daq	354
		355	for (int k=0;k<4;k++) {
		356	H2D_Fill(k, ix[0], ix[1] , gSum[k] );
		357	H1D_Fill(k, ix[0] , gSum[k] );
		358	H1D_Fill(100+ k, ix[1] , gSum[k] );
		359	}
		360	PlotScanHistogramCB(0,0,EVENT_COMMIT,NULL, 0,0);
203	f9daq	361	if (!daq_on) break;
		362	}
		363	if (!daq_on) break;
		364	}
		365
		366	return 0;
195	f9daq	367	}
		368
203	f9daq	369
195	f9daq	370	int CVICALLBACK StartCB (int panel, int control, int event,
203	f9daq	371	void *callbackData, int eventData1, int eventData2) {
		372	ThreadFunctionPtr mythread = NULL;
		373	switch (event) {
		374
		375	case EVENT_COMMIT:
		376
		377	controlID=0;
		378	if (panel == ph && control == PANEL_START) {
		379	mythread = daq;
		380	controlID= control;
		381	}
		382	if (panel == p2 && control == SCAN_SCAN) mythread = scan;
		383	if (mythread!=NULL) {
		384	printf("New Thread panel=%d button=%d\n", panel, control);
		385
		386	// SetDimming(1);
		387
		388	daq_on=1;
		389	CmtScheduleThreadPoolFunctionAdv (poolHandle, mythread, &controlID,
		390	DEFAULT_THREAD_PRIORITY,
		391	EndOfThread,
		392	EVENT_TP_THREAD_FUNCTION_END,
		393	NULL, RUN_IN_SCHEDULED_THREAD,
		394	&tfID);
		395	}
		396	break;
		397	}
		398	return 0;
195	f9daq	399	}
		400
		401	int CVICALLBACK StopCB (int panel, int control, int event,
203	f9daq	402	void *callbackData, int eventData1, int eventData2) {
		403	switch (event) {
		404	case EVENT_COMMIT:
		405	daq_on=0;
		406	break;
		407	}
		408	return 0;
195	f9daq	409	}
		410
		411	int CVICALLBACK ExitCB (int panel, int control, int event,
203	f9daq	412	void *callbackData, int eventData1, int eventData2) {
		413	switch (event) {
		414	case EVENT_COMMIT:
		415	QuitUserInterface (0);
		416	break;
		417	}
		418	return 0;
195	f9daq	419	}
203	f9daq	420
		421	int CVICALLBACK MoveStageCB (int panel, int control, int event,
		422	void *callbackData, int eventData1, int eventData2) {
		423	int axis=0, step=1000, direction=1, n;
		424	switch (event) {
		425	case EVENT_COMMIT:
		426
		427	if (panel == p2) {
		428	switch (control) {
		429	case SCAN_BR :
		430	axis = 1;
		431	direction = 1;
		432	GetCtrlVal(p2, SCAN_STEPX, &step);
		433	break;
		434	case SCAN_BL :
		435	axis = 1;
		436	direction = -1;
		437	GetCtrlVal(p2, SCAN_STEPX, &step);
		438	break;
		439	case SCAN_BU :
		440	axis = 2;
		441	direction = 1;
		442	GetCtrlVal(p2, SCAN_STEPY, &step);
		443	break;
		444	case SCAN_BD :
		445	axis = 2;
		446	direction = -1;
		447	GetCtrlVal(p2, SCAN_STEPY, &step);
		448	break;
		449	}
		450	#ifdef MIKRO
		451	MIKRO_MoveFor(axis, direction*step );
		452	MIKRO_GetPosition(axis,&n);
		453	if (axis == 1) SetCtrlVal (p2, SCAN_XP, n);
		454	if (axis == 2) SetCtrlVal (p2, SCAN_YP, n);
		455	#endif // MIKRO
		456	}
		457
		458	break;
		459	}
		460	return 0;
		461	}
		462
		463
		464
		465
		466
		467	int CVICALLBACK GoXCB (int panel, int control, int event,
		468	void *callbackData, int eventData1, int eventData2) {
		469	int n2;
		470	switch (event) {
		471	case EVENT_COMMIT:
		472	GetCtrlVal (p2, SCAN_XG, &n2);
		473	#ifdef MIKRO
		474	MIKRO_MoveTo(1,n2);
		475	MIKRO_GetPosition(1,&n2);
		476
		477	#endif
		478	SetCtrlVal (p2, SCAN_XP, n2);
		479	break;
		480	}
		481	return 0;
		482	}
		483
		484	int CVICALLBACK GoYCB (int panel, int control, int event,
		485	void *callbackData, int eventData1, int eventData2) {
		486	int n2;
		487	switch (event) {
		488	case EVENT_COMMIT:
		489	GetCtrlVal (p2, SCAN_YG, &n2);
		490	#ifdef MIKRO
		491	MIKRO_MoveTo(2,n2);
		492	MIKRO_GetPosition(2,&n2);
		493
		494	#endif
		495	SetCtrlVal (p2, SCAN_YP, n2);
		496	break;
		497	}
		498	return 0;
		499	}
		500
		501	int CVICALLBACK GetCurrentPositionCB (int panel, int control, int event,
		502	void *callbackData, int eventData1, int eventData2) {
		503
		504	int n[2];
		505	switch (event) {
		506	case EVENT_COMMIT:
		507	#ifdef MIKRO
		508	MIKRO_GetPosition(1,&n[0]);
		509	SetCtrlVal (p2, SCAN_XP, n[0]);
		510	MIKRO_GetPosition(2,&n[1]);
		511	SetCtrlVal (p2, SCAN_YP, n[1]);
		512	#endif
		513	break;
		514	}
		515	return 0;
		516	}
		517
		518	int CVICALLBACK HomeCB (int panel, int control, int event,
		519	void *callbackData, int eventData1, int eventData2) {
		520	switch (event) {
		521	case EVENT_COMMIT:
		522	#ifdef MIKRO
		523	MIKRO_ReferenceMove(1);
		524	MIKRO_ReferenceMove(2);
		525	GetCurrentPositionCB(panel, control, event, NULL, 0, 0);
		526	#endif
		527	break;
		528	}
		529	return 0;
		530	}
217	f9daq	531
		532	int CVICALLBACK PlotScanHistogramCB (int panel, int control, int event,
		533	void *callbackData, int eventData1, int eventData2) {
		534	int hid=0;
		535	int nx;
		536	int ny;
		537	switch (event) {
		538	case EVENT_COMMIT:
		539	GetCtrlVal (p2, SCAN_CHANNEL, &hid);
		540	GetCtrlVal(p2, SCAN_NSTEPSX, &nx);
		541	GetCtrlVal(p2, SCAN_NSTEPSY, &ny);
		542	if (nx>1 && ny>1) {
		543	H2D_Draw(hid,p2,SCAN_GRAPH,&scanplothandle);
		544	printf("redraw 2d\n");
		545	} else {
		546	if (nx>1) {
		547	printf("redraw 1d x\n");
		548	H1D_Draw(hid,p2,SCAN_GRAPH,&scanplothandle);
		549	}
		550	if (ny>1) {
		551	printf("redraw 1d x\n");
		552	H1D_Draw(100+hid,p2,SCAN_GRAPH,&scanplothandle);
		553	}
		554	}
		555	break;
		556	}
		557	return 0;
		558	}
		559
		560	int CVICALLBACK OpenGuiCB (int panel, int control, int event,
		561	void *callbackData, int eventData1, int eventData2) {
		562	switch (event) {
		563	case EVENT_COMMIT:
		564	DisplayPanel (ph);
		565	DisplayPanel (p2);
		566	break;
		567	}
		568	return 0;
		569	}

Subversion Repositories f9daq

(root)/cvi/instr/DRS/drs4.c – Rev 217