histplot.c

Go to the documentation of this file.

 
/*****************************************************************************/
#include "tpcclibConfig.h"
/*****************************************************************************/
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <string.h>
/*****************************************************************************/
#include "tpcextensions.h"
#include "tpctac.h"
#include "libtpcsvg.h"
/*****************************************************************************/
#include "tpctacmod.h"
/*****************************************************************************/
 
/*****************************************************************************/
extern int SVG_INLINE;
/*****************************************************************************/
int tacPlotHistogramSVG(
  TAC *d,
  const char *main_title,
  const double x1,
  const double x2,
  const double y1,
  const double y2,
  const char *fname,
  TPCSTATUS *status
) {
  if(status==NULL) return TPCERROR_FAIL;
  if(status->verbose>1)
    printf("%s(tac, '%s', %g, %g, %g, %g, '%s')\n", __func__, main_title, x1, x2, y1, y2, fname);
  if(d==NULL || d->sampleNr<1 || d->tacNr<1) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_NO_DATA); return(status->error);}
  if(fname==NULL || strnlen(fname, 1)<1) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_INVALID_FILENAME); return(status->error);}
 
  /* Determine the plot min and max x and y values */
  double minx, maxx, miny, maxy, tx1, tx2, ty1, ty2;
  if(tacSampleXRange(d, &tx1, &tx2)) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_INVALID_XRANGE); return(status->error);}
  minx=tx1; maxx=tx2;
  if(minx>0.0) {double f=maxx-minx; minx-=0.05*f; if(minx<0.0) minx=0.0;}
  if(isfinite(x1)) minx=x1; 
  if(isfinite(x2)) maxx=x2;
  if(status->verbose>1) {printf("  minx := %g\n  maxx := %g\n", minx, maxx); fflush(stdout);}
  /* Determine the plot min and max y values inside the x range */
  if(tacYRangeInXRange(d, 0, minx, maxx, &ty1, &ty2, NULL, NULL, NULL, NULL)) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_INVALID_XRANGE); return(status->error);}
//  if(ty1>0.0) ty1=0.0;
//  if(ty2<0.0) ty2=0.0;
  miny=ty1; maxy=ty2;
  if(miny>0.0) {double f=maxy-miny; miny-=0.05*f; if(miny<0.0) miny=0.0;}
  else if(maxy<0.0) maxy=0.0;
  if(isfinite(y1)) miny=y1; 
  if(isfinite(y2)) maxy=y2;
  if(status->verbose>1) {printf("  miny := %g\n  maxy := %g\n", miny, maxy); fflush(stdout);}
 
  /* Calculate the axis ticks */
  struct svg_viewports viewports; svg_init_viewports(&viewports);
  viewports.x.min=minx; viewports.x.max=maxx;
  viewports.y.min=miny; viewports.y.max=maxy;
  viewports.label_area_viewport.is=0; // no plot legends
  if(svg_calculate_axes(&viewports, status->verbose-1)) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_INVALID_VALUE); return(status->error);}
 
  /* Set the plot window and window area sizes */
  if(svg_define_viewports(0, 0, strlen(main_title), 0, 0, 0, &viewports, status->verbose-1)) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_INVALID_VALUE); return(status->error);}
 
  /* Initiate graphics file */
  FILE *fp=svg_initiate(fname, 0, 0, &viewports, NULL, status->verbose-1);
  if(fp==NULL) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_CANNOT_OPEN); return(status->error);}
 
  /* Put the graph titles into their own viewports */
  if(svg_create_main_title(fp, main_title, "", &viewports, NULL, status->verbose-2)) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_CANNOT_WRITE);
    fclose(fp); return(status->error);
  }
  /* Put the plot into its own viewport */
  if(svg_start_plot_viewport(fp, &viewports, NULL, status->verbose-2)) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_CANNOT_WRITE);
    fclose(fp); return(status->error);
  }
  /*  Start coordinate area viewport */
  if(svg_start_coordinate_viewport(fp, &viewports, NULL, status->verbose-3)) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_CANNOT_WRITE);
    fclose(fp); return(status->error);
  }
  /* Write plot axes */
  if(svg_write_axes(fp, &viewports, NULL, status->verbose-3)) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_CANNOT_WRITE);
    fclose(fp); return(status->error);
  }
 
  /*
   *  Draw the data
   */
  char ilc[9], tmp[1024];
  if(SVG_INLINE) strcpy(ilc, "svg:"); else strcpy(ilc, "");
  /* Initiate the data object group */
  sprintf(tmp, "\n<%sg stroke=\"black\" stroke-width=\"25\" fill=\"black\">\n", ilc);
  if(svg_write(fp, tmp, NULL, status->verbose-5)!=0) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_CANNOT_WRITE);
    fclose(fp); return(status->error);
  }
  if(d->isframe==0) { // vertical lines
    /* Open the path */
    sprintf(tmp, "<%spath fill=\"none\" d=\"", ilc);
    svg_write(fp, tmp, NULL, status->verbose-5);
    for(int i=0; i<d->sampleNr; i++) {
      if(!isfinite(d->x[i]) || !isfinite(d->c[0].y[i])) continue;
      double nx1, ny1, nx2, ny2;
      nx1=nx2=viewports.x.origo+d->x[i]*viewports.x.scale;
      ny1=viewports.coordinate_area_viewport.h-(viewports.y.origo+viewports.y.scale*0.0);
      ny2=viewports.coordinate_area_viewport.h-(viewports.y.origo+viewports.y.scale*d->c[0].y[i]);
      /* Do not plot if x is outside viewport */
      if(nx1<0 && nx2<=0) continue;
      if(nx1>viewports.coordinate_area_viewport.w+1 || nx1>viewports.coordinate_area_viewport.w+1)
        continue;
      /* Do not plot if both y's are outside viewport */
      if((ny1<0 || ny2>viewports.coordinate_area_viewport.h+1) &&
         (ny2<0 || ny1>viewports.coordinate_area_viewport.h+1)) continue;
      /* Set y limit to the viewport border */
      if(ny1<0) ny1=0.0; else if(ny2<0) ny2=0.0;
      if(ny1>viewports.coordinate_area_viewport.h+1) ny1=viewports.coordinate_area_viewport.h+1;
      else if(ny2>viewports.coordinate_area_viewport.h+1) ny2=viewports.coordinate_area_viewport.h+1;
      /* Draw */
      sprintf(tmp, " M%.0f %.0f L%.0f %.0f", nx1, ny1, nx2, ny2);
      svg_write(fp, tmp, NULL, status->verbose-5);
    }
    /* Close the path */
    strcpy(tmp, "\" />\n");
    svg_write(fp, tmp, NULL, status->verbose-5);
  } else { // rectangles
    for(int i=0; i<d->sampleNr; i++) {
      if(!isfinite(d->x[i]) || !isfinite(d->c[0].y[i])) continue;
      double nx1, ny1, nx2, ny2;
      nx1=viewports.x.origo+d->x1[i]*viewports.x.scale;
      nx2=viewports.x.origo+d->x2[i]*viewports.x.scale;
      ny1=viewports.coordinate_area_viewport.h-(viewports.y.origo+viewports.y.scale*0.0);
      ny2=viewports.coordinate_area_viewport.h-(viewports.y.origo+viewports.y.scale*d->c[0].y[i]);
      /* Do not plot if x is outside viewport */
      if(nx1<0 && nx2<=0) continue;
      if(nx1>viewports.coordinate_area_viewport.w+1 || nx1>viewports.coordinate_area_viewport.w+1)
        continue;
      /* Do not plot if both y's are outside viewport */
      if((ny1<0 || ny2>viewports.coordinate_area_viewport.h+1) &&
         (ny2<0 || ny1>viewports.coordinate_area_viewport.h+1)) continue;
      /* Set y limit to the viewport border */
      if(ny1<0) ny1=0.0; else if(ny2<0) ny2=0.0;
      if(ny1>viewports.coordinate_area_viewport.h+1) ny1=viewports.coordinate_area_viewport.h+1;
      else if(ny2>viewports.coordinate_area_viewport.h+1) ny2=viewports.coordinate_area_viewport.h+1;
      if(ny1>ny2) {double f=ny1; ny1=ny2; ny2=f;}
      /* Draw */
      sprintf(tmp, " <rect x=\"%.0f\" y=\"%.0f\" width=\"%.0f\" height=\"%.0f\" />\n",
              nx1, ny1, nx2-nx1, fabs(ny2-ny1));
      svg_write(fp, tmp, NULL, status->verbose-5);
    }
  }
  /* Close the data object group */
  sprintf(tmp, "</%sg>\n", ilc);
  if(svg_write(fp, tmp, NULL, status->verbose-5)!=0) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_CANNOT_WRITE);
    fclose(fp); return(status->error);
  }
 
  /* Close the coordinate viewport */
  if(svg_end_coordinate_viewport(fp, NULL, status->verbose-1)) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_CANNOT_WRITE);
    fclose(fp); return(status->error);
  }
 
  /* Write the axis ticks */
  if(svg_write_xticks(fp, &viewports, NULL, status->verbose-3)!=0 ||
     svg_write_yticks(fp, &viewports, NULL, status->verbose-3)!=0) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_CANNOT_WRITE);
    fclose(fp); return(status->error);
  }
 
  /* Close the plot viewport */
  if(svg_end_plot_viewport(fp, NULL, status->verbose-2)) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_CANNOT_WRITE);
    fclose(fp); return(status->error);
  }
 
  /* Close the SVG file */
  if(svg_close(fp, NULL, status->verbose-1)) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_CANNOT_WRITE);
    fclose(fp); return(status->error);
  }
 
  if(status->verbose>0) fprintf(stdout, "  written file %s\n", fname);
 
  statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_OK);
  return(TPCERROR_OK);
}
/*****************************************************************************/
 
/*****************************************************************************/