WebSVN – f9daq – Blame – //cvi/instr/HISTO/H1D.c

Rev	Author	Line No.	Line
209	f9daq	1	#include <userint.h>
		2
		3
		4	#ifdef _CVI_
		5	# include <ansi_c.h>
		6	# else _CVI_
		7	# include <stdlib.h>
		8	# include <stdio.h>
		9	# include <string.h>
		10	#endif _CVI_
		11
		12	#include "H1D.h"
		13
		14	#define H1DMAX 500
		15	H1D *h1[H1DMAX];
		16	//int Printf(char *format, ...);
		17
		18	int _VI_FUNC H1D_Clear(int h1d) {
		19	if (!h1[h1d]) return -1;
		20	memset(h1[h1d]->data, 0,h1[h1d]->size );
		21	h1[h1d]->min=0;
		22	h1[h1d]->max=0;
		23	h1[h1d]->nentries=0;
		24	return 0;
		25
		26	}
		27
		28	int _VI_FUNC H1D_Print(int h1d) {
		29	if (!h1[h1d]) return -1;
		30	//Printf("PrintH1D_ nx=%d minx=%f stepx=%f ny=%d miny=%f stepy=%f size=%d\n", h1[h1d]->nx, h1[h1d]->minx, h1[h1d]->stepx, h1[h1d]->ny, h1[h1d]->miny, h1[h1d]->stepy, h1[h1d]->size ) ;
		31	return 0;
		32
		33	}
		34
		35	int _VI_FUNC H1D_Exist(int h) {
		36	if (h1[h]) return 1;
		37	else return 0;
		38	}
		39
		40
		41	int _VI_FUNC H1D_GetBin(int h, double value) {
		42	int nx,xmin,dx;
		43	int bin;
		44
		45	nx = H1D_GetNbinsX(h);
		46	xmin= H1D_GetMinX(h);
		47	dx = H1D_GetStepX(h);
		48
		49	if (dx==0) return -1;
		50	if (value<xmin) return -1;
		51	bin = (int)((value-xmin)/dx);
		52	if (bin>=nx) return -1;
		53	else return bin;
		54	}
		55
		56	int _VI_FUNC H1D_Fill(int h1d,double x, double val) {
		57
		58	int ix;
		59	if (!h1[h1d]) return -1;
		60
		61	ix = H1D_GetBin(h1d,x);
		62	if (ix<0) return ix;
		63
		64	h1[h1d]->data[ix]+=val;
		65	//Printf("%d %d data %f %f\n",x,y,val, h1[h1d]->data[idx]);
		66	if (h1[h1d]->data[ix]>h1[h1d]->max) h1[h1d]->max= h1[h1d]->data[ix];
		67	if (h1[h1d]->data[ix]<h1[h1d]->min) h1[h1d]->min= h1[h1d]->data[ix];
		68	h1[h1d]->nentries++;
		69	return 0;
		70	}
		71
		72
		73	int _VI_FUNC H1D_FillBin(int h1d,int x, double val) {
		74
		75	if (!h1[h1d]) {
		76	//Printf("FillH1D_ error h1d is not initialized\n");
		77	return -1;
		78	}
		79
		80	h1[h1d]->data[x]+=val;
		81	//Printf("%d %d data %f %f\n",x,y,val, h1[h1d]->data[idx]);
		82	if (h1[h1d]->data[x]>h1[h1d]->max) h1[h1d]->max= h1[h1d]->data[x];
		83	if (h1[h1d]->data[x]<h1[h1d]->min) h1[h1d]->min= h1[h1d]->data[x];
		84	h1[h1d]->nentries++;
		85	return 0;
		86	}
		87
		88	int _VI_FUNC H1D_SetBinContent(int h1d,int x, double val) {
		89
		90	if (!h1[h1d]) {
		91	//Printf("FillH1D_ error h1d is not initialized\n");
		92	return -1;
		93	}
		94
		95	h1[h1d]->data[x]=val;
		96	//Printf("%d %d data %f %f\n",x,y,val, h1[h1d]->data[idx]);
		97	if (h1[h1d]->data[x]>h1[h1d]->max) h1[h1d]->max= h1[h1d]->data[x];
		98	if (h1[h1d]->data[x]<h1[h1d]->min) h1[h1d]->min= h1[h1d]->data[x];
		99	h1[h1d]->nentries++;
		100	return 0;
		101	}
		102
		103
		104
		105	double _VI_FUNC H1D_GetBinContent(int h1d,int atx) {
		106
		107	int idx;
		108	if (!h1[h1d]) return 0;
		109	if (h1[h1d]->nx <= atx) return 0;
		110	if (atx<0) return 0;
		111	if (atx*sizeof(double) < h1[h1d]->size ) return h1[h1d]->data[atx];
		112
		113
		114	return 0;
		115	}
		116
		117
		118	int _VI_FUNC H1D_Init(int h1d,char name, char title,int nx, double minx, double maxx) {
		119
		120	if (h1[h1d]) {
		121
		122	free(h1[h1d]->data);
		123	free(h1[h1d]);
		124	h1[h1d] = NULL;
		125	}
		126	// if (h1d!=H1D_MAX-1) printf("InitH1D_ hID=%d\n",h1d);
		127	h1[h1d] = (H1D *) malloc(sizeof(H1D));
		128	//h2 =h1d;
		129
		130	H1D_SetTitle(h1d,title);
		131	H1D_SetName(h1d,name);
		132	H1D_SetTitleX(h1d,"x");
		133	;
		134	h1[h1d]->id=H1D_ID;
		135	h1[h1d]->nx = nx;
		136
		137
		138	h1[h1d]->minx = minx;
		139	h1[h1d]->stepx = (nx>1)?(maxx-minx)/(nx-1):0;
		140
		141
		142	h1[h1d]->size = h1[h1d]->nx*sizeof(double);
		143	h1[h1d]->data = (double *) malloc(h1[h1d]->size);
		144	h1[h1d]->len=sizeof(H1D)-sizeof(double *)+h1[h1d]->size;
		145	H1D_Clear(h1d);
		146	H1D_Print(h1d);
		147	//Printf("InitH1D 0x%x\n", h1d );
		148	return 0;
		149
		150	}
		151
		152
		153	double _VI_FUNC H1D_GetXBinCenter(int h1d,int xbin) {
		154	return h1[h1d]->minx+xbin*h1[h1d]->stepx;
		155	}
		156
		157
		158	int _VI_FUNC H1D_Write2File(int h1d,FILE *fp) {
		159
		160	if (!fp) return -1;
		161	//printf("H1D sizeof(H1D)=%lu len-datasize=%d len=%lu datasize=%d\t",sizeof(H1D)-sizeof(double *),h1[h1d]->len-h1[h1d]->size,h1[h1d]->len,h1[h1d]->size);
		162	//printf("H1D sz=%d %d\n",sizeof(double),sizeof(double *));
		163	if (!H1D_Exist(h1d)){
		164	printf("Histogram H1D=%d is not initialized\n",h1d);
		165	return -1;
		166	}
		167	fwrite (h1[h1d], 1, sizeof(H1D)-sizeof(double *), fp);
		168	fwrite (h1[h1d]->data, 1, h1[h1d]->size, fp);
		169	return 0;
		170	}
		171
		172	int _VI_FUNC H1D_Write(int h1d,const char fname,const char opt) {
		173	FILE *fp=fopen(fname,opt);
		174	H1D_Write2File(h1d,fp);
		175	fclose(fp);
		176	return 0;
		177	}
		178
		179
		180
		181
		182	int _VI_FUNC H1D_SetTitle(int h1d,char *title) {
		183	if (!h1[h1d]) {
		184	printf("h1d %d does not exist %s\n",h1d, title);
		185	return 0;
		186	}
		187	sprintf(h1[h1d]->title,"%s",title);
		188	return 0;
		189	}
		190
		191
		192	int _VI_FUNC H1D_SetTitleX(int h1d,char *title) {
		193	sprintf(h1[h1d]->titlex,"%s",title);
		194	return 0;
		195	}
		196
		197
		198	int _VI_FUNC H1D_SetTitleY(int h1d,char *title) {
		199	sprintf(h1[h1d]->titley,"%s",title);
		200	return 0;
		201	}
		202
		203
		204	char * _VI_FUNC H1D_GetTitleX(int h1d) {
		205	return h1[h1d]->titlex;
		206	}
		207
		208	char * _VI_FUNC H1D_GetTitleY(int h1d) {
		209	return h1[h1d]->titley;
		210	}
		211
		212	char * _VI_FUNC H1D_GetTitle(int h1d) {
		213	return h1[h1d]->title;
		214	}
		215
		216
		217
		218	int _VI_FUNC H1D_SetName(int h1d,char *title) {
		219	sprintf(h1[h1d]->name,"%s",title);
		220	return 0;
		221	}
		222
		223	int _VI_FUNC H1D_GetNbinsX(int h) {
		224	if (h1[h]) return h1[h]->nx;
		225	else return 0;
		226	}
		227
		228
		229
		230	double _VI_FUNC H1D_GetMinX(int h) {
		231	if (h1[h]) return h1[h]->minx;
		232	else return 0;
		233	}
		234
		235	double _VI_FUNC H1D_GetMaxX(int h){
		236	return h1[h]->minx+(h1[h]->nx-1)*h1[h]->stepx;
		237	}
		238
		239
		240	double _VI_FUNC H1D_GetStepX(int h) {
		241	if (h1[h]) return h1[h]->stepx;
		242	else return 0;
		243	}
		244
		245
		246	double _VI_FUNC H1D_GetMin(int h) {
		247	if (h1[h]) return h1[h]->min;
		248	else return 0;
		249	}
		250
		251	double _VI_FUNC H1D_GetMax(int h) {
		252	if (h1[h]) return h1[h]->max;
		253	else return 0;
		254	}
		255
		256	double * _VI_FUNC H1D_GetData(int h) {
		257	if (h1[h]) return h1[h]->data;
		258	else return NULL;
		259	}
		260
		261
		262	int _VI_FUNC H1D_Draw(int histogram,int panel, int control, int *plothandle) {
		263	const unsigned int hcolors[8]= {VAL_RED, VAL_GREEN, VAL_BLUE, VAL_CYAN, VAL_MAGENTA, VAL_YELLOW, VAL_GRAY, VAL_BLACK};
		264	if (!H1D_Exist(histogram)) {
		265	printf("1D Histogram %d does not exist!\n",histogram);
		266	return 0;
		267	}
		268
		269	#ifdef _CVI_
		270	if (plothandle> 0 ) DeleteGraphPlot (panel, control, plothandle, VAL_IMMEDIATE_DRAW);
		271
		272	*plothandle = PlotWaveform (panel, control, H1D_GetData(histogram), H1D_GetNbinsX(histogram), VAL_DOUBLE, 1, 0,
		273	H1D_GetMinX(histogram), H1D_GetStepX(histogram), VAL_FAT_LINE, VAL_SMALL_SOLID_SQUARE, VAL_SOLID, 1,
		274	hcolors[histogram%8]);
		275	RefreshGraph (panel, control);
		276	ProcessSystemEvents ();
		277	#endif
		278
		279	return *plothandle;
		280
		281	}
		282

Subversion Repositories f9daq

(root)//cvi/instr/HISTO/H1D.c – Rev 209