示例#1
0
thStatus thRemoveHists(char *block_name){
  thHistSpecList *thishist,*nexthist;
  thHBlockList *thisblock,*nextblock,**lastblockp;

  if(block_name)
    if(*block_name=='\0')
      block_name = 0;
    else
      lastblockp = &thHBlockListP;

  if(thHBlockListP == 0){
    return(S_FAILURE);		/* No histograms defined */
  } else {
    thisblock = thHBlockListP;
    if(!block_name)
      HDELET(0);		/* Free all space */
    while(thisblock) {
      if(block_name)
	if(strcasecmp(block_name,thisblock->blockname)!=0){
	  lastblockp = &thisblock->next;
	  thisblock = thisblock->next;
	  continue;
	}
      thishist = thisblock->var->opaque;
      while(thishist){
	nexthist = thishist->next;
	if(block_name) {
	  HDELET(*(DAINT *) thishist->varname->varptr);
	}
	free(thishist);
	thishist = nexthist;
      }
      nextblock = thisblock->next;
      if(block_name){
	if(strcasecmp(block_name,thisblock->blockname)==0) {
	  *lastblockp = nextblock;
	  free(thisblock);
	}
	return(S_SUCCESS);
      } else {
	free(thisblock);
      }
      thisblock = nextblock;
    }
    if(block_name) {
      fprintf(stderr,"Hist block %s not found\n",block_name);
      return(S_FAILURE);
    }
    thHBlockListP = (thHBlockList *) NULL;
  }
  return(S_SUCCESS);
}
示例#2
0
文件: antidbc02.C 项目: krafczyk/AMS
//----
// function to display PMT pulse (flash-ADC array arr[]) :
//
void ANTI2DBc::displ_a(char comm[], int id, int mf, const geant arr[]){
  integer i;
  geant tm,a(0.);
  geant tb,tbi;
  char name[80], buf[10];
  sprintf(buf, "%4d\n",id);
  strcpy(name,comm);
  strcat(name,buf);
  tb=geant(mf)*ANTI2DBc::fadcbw();
  tbi=ANTI2DBc::fadcbw();
  HBOOK1(2638,name,100,0.,100*tb,0.);
  for(i=1;i<=ANTI2C::ANFADC;i++){
    if(i%mf==0){
      a+=arr[i-1];
      tm=i*tbi-0.5*tb;
      HF1(2638,tm,a/geant(mf));
      a=0.;
    }
    else{
      a+=arr[i-1];
    }
  }
  HPRINT(2638);
  HDELET(2638);
  return ;
}
示例#3
0
/*----------------------------- get_next_scan -------------------------------*/
int get_next_scan( int nscan)
  {

  if ( nscan != 1)
    {              // we don't need to fit and save at the begining 
    fit_thresh(nscan-1);
    save_hist(nscan-1);
    print_scan_result(nscan-1);
    HDELET(0);
    }
  init_hist (nscan);
  return 0;
  }
示例#4
0
/*----------------------------- save_hist_19 ----------------------------------*/
int save_hist_19()
{

  char file_name_element[] = "/test_results/temp/test_19_";
  char *file_name;
  file_name = malloc(strlen(datadir_translation) +
                     strlen(file_name_element) + 10);
  if (file_name == NULL)
    {
      printf("malloc() failed for file_name.\n");
      return -1;
    }
  sprintf(file_name,"%s%s%02d.his",datadir_translation,
          file_name_element,1);
  /* Save histograms */
  printf("Saving histograms in file %s\n", file_name);
  HRPUT(0, file_name, "N");
  HDELET(0);
  free(file_name);
  return 0;
}
示例#5
0
int test_24_finish(void)
  {
  int      j, ilayer, isegment, icfeb, n, nfail;
  logical  channel_ok;
  float    limits[7];
  char     timestring[40];
  time_t   t;
  char *file_his_name;

  /* Save histograms */
  file_his_name = malloc(strlen(datadir_translation) + 
			 strlen(test_24_file_his_el) + 1);
  if (file_his_name == NULL)
    {
    printf("malloc() failed for file_his_name.  Could not save histograms.\n");
    }
  else
    {
    sprintf(file_his_name,"%s%s",
	    datadir_translation, test_24_file_his_el);
    printf("Saving histograms\n");
    HRPUT(0, file_his_name, "N");
    free(file_his_name);    
    }

/* Look for values outside limits */
  fprintf(fp, "\nList of sections with peak outside limits %d-%d, or with "
	  "chisq/ndf > 10:\n", LOWER_LIMIT, UPPER_LIMIT);
  nfail = 0;

  for (ilayer = 0; ilayer < NLAYER; ilayer++)
    {
    for (isegment = 0; isegment < nsegments; isegment++)
      {
      for (icfeb = 0; icfeb < ncfebs; icfeb++)
	{
        channel_ok = _TRUE;
        j = ncfebs * isegment + icfeb;

        if (peak[ilayer][j] < LOWER_LIMIT || peak[ilayer][j] > UPPER_LIMIT)
	  {
          channel_ok = _FALSE;
          fprintf(fp, "Layer %d segment %d cfeb %d has peak= %4.0f\n", 
		  ilayer+1, isegment+1, icfeb+1, peak[ilayer][j]);
	  }

        if (chisq[ilayer][j] > 10)
	  {
          channel_ok = _FALSE;
          fprintf(fp, "Layer %d segment %d cfeb %d has chisq=%8.2f\n", 
		  ilayer+1, isegment+1, icfeb+1, chisq[ilayer][j]);
	  }

        if (!channel_ok) nfail++;
	}
      }
    }
  fprintf(fp, "\n");
  fprintf(fp, "Total bad channels = %d\n", nfail);

/* Print voltage Equalization table */
  fprintf(fp, "\nEqualization table rounded to nearest 10 volts:\n");
  if (nsegments == 5) fprintf(fp, "HV segment:     1     2     3     4     5    all\n");
  else                fprintf(fp, "HV segment:     1     2     3    all\n");
  fprintf(fp, "------------------------------------------------\n");
  for (ilayer = 0; ilayer < NLAYER; ilayer++)
    {
    fprintf(fp, "Layer %d:   ", ilayer+1);
    for (isegment = 0; isegment < nsegments; isegment++)
      {
      if (dv[ilayer][isegment]  >= 0) 
	n = 10 * (int)((dv[ilayer][isegment] / 10) + 0.5);
      else n = 10 * (int)((dv[ilayer][isegment] / 10) - 0.5);
      fprintf(fp, "  %4d", n);
      }
    if (dv_all[ilayer] >= 0) 
      n = 10 * (int)((dv_all[ilayer] / 10) + 0.5);
    else n = 10 * (int)((dv_all[ilayer] / 10) - 0.5);
    fprintf(fp, "  %4d\n", n);
    }
  fprintf(fp, "\n");


/* Plot results */
  limits[0] = 0;                   /* layer histo ymin */
  limits[1] = 1400;                /* layer histo ymax */
  limits[2] = 0;                   /* summary histo ymin */
  limits[3] = 1400;                /* summary histo ymax */
  limits[4] = LOWER_LIMIT;         /* lower limit line */
  limits[5] = UPPER_LIMIT;         /* upper limit line */
  limits[6] = nsegments * ncfebs;
  
  one_page_plot_v3(24, 1, GAIN_PLOT, (float *)ppeak, (float *)ppeak, limits, 
		   csc_name, "Landau peaks", "Landau peak, ADC "
		   "counts (5x5 sum)", 0);

/* Print mandatory summary statement */
  if (nfail > 0) fprintf(fp, "Test result = Failed\n\n");
  else fprintf(fp, "Test result = OK\n\n");

/* Add the date and time */
  time(&t);
  sprintf(timestring, "%s", ctime(&t));
  fprintf(fp, "%s\n\n", timestring);

  fclose(fp);

  HDELET(0);

  return 0;
  }