Subversion Repositories f9daq

Rev

Rev 221 | Go to most recent revision | Blame | Compare with Previous | Last modification | View Log | RSS feed

#ifdef _CVI_
#  include <ansi_c.h>
# else _CVI_
#  include <stdlib.h>
#  include <stdio.h>
#  include <string.h>
#endif _CVI_

#include "H3D.h"
#include "H2D.h"
#include "H1D.h"

#define H3DMAX 500
H3D *h3[H3DMAX];
//int Printf(char *format, ...);

int _VI_FUNC  H3D_Clear(int h3d) {
  if (!h3[h3d]) return -1;
  memset(h3[h3d]->data, 0,h3[h3d]->size );
  h3[h3d]->min=0;
  h3[h3d]->max=0;
  h3[h3d]->nentries=0;
  return  0;

}

int _VI_FUNC  H3D_Print(int h3d) {
  if (!h3[h3d]) return -1;
//Printf("PrintH3D nx=%d minx=%f stepx=%f ny=%d miny=%f stepy=%f size=%d\n", h3[h3d]->nx, h3[h3d]->minx, h3[h3d]->stepx, h3[h3d]->ny, h3[h3d]->miny, h3[h3d]->stepy, h3[h3d]->size ) ;
  return  0;

}

int _VI_FUNC  H3D_Exist(int h) {
  if (h3[h]) return 1;
  else return 0;
}


int _VI_FUNC  H3D_CalculateBin(int h, int axis, double value) {
  int nx=1,xmin=0,dx=0;
  int bin;
  switch (axis) {
    case 0:
      nx = H3D_GetNbinsX(h);
      xmin= H3D_GetMinX(h);
      dx = H3D_GetStepX(h);
      break;
    case 1:
      nx = H3D_GetNbinsY(h);
      xmin= H3D_GetMinY(h);
      dx = H3D_GetStepY(h);
      break;
    case 2:
      nx = H3D_GetNbinsZ(h);
      xmin= H3D_GetMinZ(h);
      dx = H3D_GetStepZ(h);
      break;

  }
  if (dx==0) return -1;
  if (value<xmin) return -1;
  bin = (int)((value-xmin)/dx);
  if (bin>=nx) return -1;
  else return bin;
}

int _VI_FUNC  H3D_Fill(int h3d,double x, double y, double z, double val) {

  int ix,iy, iz;
  if (!h3[h3d]) return -1;

  ix = H3D_CalculateBin(h3d,0,x);
  if (ix<0) return ix;

  iy = H3D_CalculateBin(h3d,1,y);
  if (iy<0) return iy;

  iz = H3D_CalculateBin(h3d,2,z);
  if (iz<0) return iz;

 
  H3D_FillBin(h3d, ix, iy, iz, val);
  return 0;
}

int _VI_FUNC  H3D_GetBin(int h3d,int x, int y, int z) {
  return x+h3[h3d]->nx *(y+z*h3[h3d]->ny);
}

int _VI_FUNC  H3D_FillBin(int h3d,int x, int y, int z, double val) {

  int idx;
  if (!h3[h3d]) {
//Printf("FillH3D error h3d is not initialized\n");
    return -1;
  }

  idx = H3D_GetBin(h3d,x,y,z);
  h3[h3d]->data[idx]+=val;
//Printf("%d %d data %f %f\n",x,y,val, h3[h3d]->data[idx]);
  if (h3[h3d]->data[idx]>h3[h3d]->max) h3[h3d]->max= h3[h3d]->data[idx];
  if (h3[h3d]->data[idx]<h3[h3d]->min) h3[h3d]->min= h3[h3d]->data[idx];
  h3[h3d]->nentries++;
  return 0;
}

double _VI_FUNC  H3D_GetBinContent(int h3d,int atx, int aty, int atz) {

  int idx;
  if (!h3[h3d]) return 0;
  if (h3[h3d]->nx <= atx)  return 0;
  if (h3[h3d]->ny <= aty)  return 0;
  if (h3[h3d]->nz <= atz)  return 0;
  idx = H3D_GetBin(h3d,atx,aty,atz);
  if (idx*sizeof(double) < h3[h3d]->size ) return h3[h3d]->data[idx];


  return 0;
}


int _VI_FUNC  H3D_Init(int h3d,char *name, char *title,
                       int nx, double minx, double maxx,
                       int ny, double miny, double maxy,
                       int nz, double minz, double maxz) {

  if (h3[h3d]) {

    free(h3[h3d]->data);
    free(h3[h3d]);
    h3[h3d] = NULL;
  }
  //printf("InitH3D hID=%d\n",h3d);
  h3[h3d] = (H3D *) malloc(sizeof(H3D));
//h2  =h3d;

  H3D_SetTitle(h3d,title);
  H3D_SetName(h3d,name);
  h3[h3d]->id=H3D_ID;
  h3[h3d]->nx = nx;
  h3[h3d]->ny = ny;
  h3[h3d]->nz = nz;

  h3[h3d]->minx = minx;
  h3[h3d]->miny = miny;
  h3[h3d]->minz = minz;

  h3[h3d]->stepx = (nx>1)?(maxx-minx)/(nx):0;
  h3[h3d]->stepy = (nx>1)?(maxy-miny)/(ny):0;
  h3[h3d]->stepz = (nx>1)?(maxz-minz)/(nz):0;

  h3[h3d]->size = h3[h3d]->nx*h3[h3d]->ny*h3[h3d]->nz*sizeof(double);
  h3[h3d]->data = (double *) malloc(h3[h3d]->size);
  h3[h3d]->len=sizeof(H3D)-sizeof(double *)+h3[h3d]->size;
  H3D_Clear(h3d);
  H3D_Print(h3d);
//Printf("InitH3D 0x%x\n", h3d );
  return  0;

}


int _VI_FUNC  H3D_SliXY(int histogram,int slice, int direction) {      
 
if (!H3D_Exist(histogram)) {
    printf("3D Histogram %d does not exist!\n",histogram);
    return -1;
  }

  int hid=499;
  switch (direction){
    case 0:// xy
      H2D_Init(hid,"projectionXY","projectionXY",
               H3D_GetNbinsX(histogram), H3D_GetMinX(histogram),H3D_GetMaxX(histogram),
               H3D_GetNbinsY(histogram), H3D_GetMinY(histogram),H3D_GetMaxY(histogram)
              );
      for (int i=0; i < H3D_GetNbinsX(histogram); i++ )
        for (int j=0; j < H3D_GetNbinsY(histogram); j++ )
          H2D_SetBinContent(hid,i,j,H3D_GetBinContent(histogram,i,j,slice));
      break;
    case 1:// xz
      H2D_Init(hid,"projectionXZ","projectionXZ",
               H3D_GetNbinsX(histogram), H3D_GetMinX(histogram),H3D_GetMaxX(histogram),
               H3D_GetNbinsZ(histogram), H3D_GetMinZ(histogram),H3D_GetMaxZ(histogram)
              );
      for (int i=0; i < H3D_GetNbinsX(histogram); i++ )
        for (int j=0; j < H3D_GetNbinsZ(histogram); j++ )
          H2D_SetBinContent(hid,i,j,H3D_GetBinContent(histogram,i,slice,j));
      break;
    case 2:// yz
    default:
      H2D_Init(hid,"projectionYZ","projectionYZ",
               H3D_GetNbinsY(histogram), H3D_GetMinY(histogram),H3D_GetMaxY(histogram),
               H3D_GetNbinsZ(histogram), H3D_GetMinZ(histogram),H3D_GetMaxZ(histogram)
              );
      for (int i=0; i < H3D_GetNbinsY(histogram); i++ )
        for (int j=0; j < H3D_GetNbinsZ(histogram); j++ )
          H2D_SetBinContent(hid,i,j,H3D_GetBinContent(histogram,slice,i,j));
     
      break;
  }
  return hid;
}

double  _VI_FUNC H3D_GetMaxX(int h){
  return h3[h]->minx+(h3[h]->nx)*h3[h]->stepx;
}

double  _VI_FUNC H3D_GetMaxY(int h){
  return h3[h]->miny+(h3[h]->ny)*h3[h]->stepy;
}

double  _VI_FUNC H3D_GetMaxZ(int h){
  return h3[h]->minz+(h3[h]->nz)*h3[h]->stepz;
}

double _VI_FUNC  H3D_GetXBinCenter(int h3d,int xbin) {
  return h3[h3d]->minx+(xbin+0.5)*h3[h3d]->stepx;
}

double _VI_FUNC  H3D_GetYBinCenter(int h3d,int ybin) {
  return h3[h3d]->miny+(ybin+0.5)*h3[h3d]->stepy;
}

double _VI_FUNC  H3D_GetZBinCenter(int h3d,int zbin) {
  return h3[h3d]->minz+(zbin+0.5)*h3[h3d]->stepz;
}


int _VI_FUNC  H3D_Write2File(int h3d,FILE *fp) {

  if (!fp) return -1;
//printf("H3D sizeof(H3D)=%lu len-datasize=%d len=%lu datasize=%d\t",sizeof(H3D)-sizeof(double *),h3[h3d]->len-h3[h3d]->size,h3[h3d]->len,h3[h3d]->size);
//printf("H3D sz=%d %d\n",sizeof(double),sizeof(double *));
  fwrite (h3[h3d], 1, sizeof(H3D)-sizeof(double *), fp);
  fwrite (h3[h3d]->data, 1, h3[h3d]->size, fp);
  return  0;
}

int _VI_FUNC  H3D_Write(int h3d,const char *fname,const char *opt) {
  FILE *fp=fopen(fname,opt);
  H3D_Write2File(h3d,fp);
  fclose(fp);
  return  0;
}




int _VI_FUNC  H3D_SetTitle(int h3d,char *title) {
  sprintf(h3[h3d]->title,"%s",title);
  return 0;
}


int _VI_FUNC  H3D_SetTitleX(int h3d,char *title) {
  sprintf(h3[h3d]->titlex,"%s",title);
  return 0;
}


int _VI_FUNC  H3D_SetTitleY(int h3d,char *title) {
  sprintf(h3[h3d]->titley,"%s",title);
  return 0;
}

int _VI_FUNC  H3D_SetTitleZ(int h3d,char *title) {
  sprintf(h3[h3d]->titlez,"%s",title);
  return 0;
}


int _VI_FUNC  H3D_SetName(int h3d,char *title) {
  sprintf(h3[h3d]->name,"%s",title);
  return 0;
}

int _VI_FUNC  H3D_GetNbinsY(int h) {
  if (h3[h]) return h3[h]->ny;
  else return 0;
}

int _VI_FUNC  H3D_GetNbinsX(int h) {
  if (h3[h]) return h3[h]->nx;
  else return 0;
}

int _VI_FUNC  H3D_GetNbinsZ(int h) {
  if (h3[h]) return h3[h]->nz;
  else return 0;
}



double _VI_FUNC  H3D_GetMinX(int h) {
  if (h3[h]) return h3[h]->minx;
  else return 0;
}

double _VI_FUNC  H3D_GetMinY(int h) {
  if (h3[h]) return h3[h]->miny;
  else return 0;
}

double _VI_FUNC  H3D_GetMinZ(int h) {
  if (h3[h]) return h3[h]->minz;
  else return 0;
}




double _VI_FUNC  H3D_GetStepX(int h) {
  if (h3[h]) return h3[h]->stepx;
  else return 0;
}

double _VI_FUNC  H3D_GetStepY(int h) {
  if (h3[h]) return h3[h]->stepy;
  else return 0;
}

double _VI_FUNC  H3D_GetStepZ(int h) {
  if (h3[h]) return h3[h]->stepz;
  else return 0;
}




double _VI_FUNC  H3D_GetMin(int h) {
  if (h3[h]) return h3[h]->min;
  else return 0;
}

double _VI_FUNC  H3D_GetMax(int h) {
  if (h3[h]) return h3[h]->max;
  else return 0;
}

double *  _VI_FUNC  H3D_GetSliceXYData(int h, int slice) {
  if (h3[h]) return (h3[h]->data+slice*h3[h]->nx*h3[h]->ny);
  else return NULL;
}

double *  _VI_FUNC  H3D_GetData(int h) {
  if (h3[h]) return h3[h]->data;
  else return NULL;
}


#ifdef _CVI_
// defined only in CVI
static HColorMap *colormap = NULL;


HColorMap * _VI_FUNC  H3D_GetColorMap(void) {
  return colormap;
}


int _VI_FUNC  H3D_SetRangeColors( double min, double max) {
  int i;


  if (colormap == NULL) {
    colormap = malloc(sizeof(HColorMap));

    colormap->numberofColors = 5;
    colormap->array = malloc(colormap->numberofColors*sizeof(ColorMapEntry));

    colormap->array[0].color = 0x0000ff; //BLUE
    colormap->array[1].color = 0x00ff00; //GREEN
    colormap->array[2].color = 0xffff00; //YELLOW
    colormap->array[3].color = 0xff8000; //ORANGE
    colormap->array[4].color = 0xff0000; //RED

    colormap->HiColor =colormap->array[colormap->numberofColors-1].color ;
  }
  if (colormap->numberofColors<2) return -1;
  double fx = (max-min)/(colormap->numberofColors-1);
  for  (i=0; i<colormap->numberofColors; i++) {
    colormap->array[i].dataValue.valDouble=i*fx+min;
  }
  return 0;
}

int  _VI_FUNC  H3D_DrawSliceXY(int histogram,int slice,int direction, int panel, int control, int *plothandle) {
 
 
 
  if (H3D_GetNbinsY(histogram)==1|| H3D_GetNbinsX(histogram)==1) {
    if (H3D_GetNbinsY(histogram)==1) {
      H1D_Init(499,"projection","projection", H3D_GetNbinsX(histogram), H3D_GetMinX(histogram),H3D_GetMaxX(histogram));
      for (int i=0; i < H3D_GetNbinsX(histogram); i++ ) H1D_SetBinContent(499,i,H3D_GetBinContent(histogram,i,0,slice));
    } else {
      H1D_Init(499,"projection","projection", H3D_GetNbinsY(histogram), H3D_GetMinY(histogram),H3D_GetMaxY(histogram));
      for (int i=0; i < H3D_GetNbinsY(histogram); i++ ) H1D_SetBinContent(499,i,H3D_GetBinContent(histogram,0,i,slice));
    }
    H1D_Draw(499, panel, control, plothandle);
  } else {

    if (*plothandle> 0 ) DeleteGraphPlot (panel, control, *plothandle, VAL_IMMEDIATE_DRAW);
    H3D_SetRangeColors(H3D_GetMin(histogram),H3D_GetMax(histogram));
    *plothandle = PlotScaledIntensity (panel, control,
                                       H3D_GetSliceXYData(histogram, slice),
                                       H3D_GetNbinsX(histogram),
                                       H3D_GetNbinsY(histogram),
                                       VAL_DOUBLE,
                                       H3D_GetStepY(histogram),
                                       H3D_GetYBinCenter(histogram,0),
                                       H3D_GetStepX(histogram),
                                       H3D_GetXBinCenter(histogram,0),
                                       colormap->array,
                                       colormap->HiColor,
                                       colormap->numberofColors, 1, 0);
  }
  ProcessSystemEvents ();
  return *plothandle;

}
#endif