#include "TMath.h"
#include "Peak.h"
#include "Database.h"
#include "opt_parameter.h"
#include <cstdlib>
#include <cstdio>
#include <cstring>

extern database db;

using namespace std;

void write_test_kumac(double beamx, double beamy) 
{
  FILE *fd;
  char fname[17];
  //  y_tg
  for(int y_peak=1;y_peak<=db.nysievpeak; y_peak++) {
    fname[0]=y_peak + '0';
    strcpy(fname+1,"_y_peaks.kumac");
    printf("%s\n",fname);
    fd=fopen(fname,"w");     
    for(int z_peak=1;z_peak<=db.nzpeak; z_peak++) {
      Peak pk(z_peak, y_peak);
      pk.set(beamx, beamy, 0.0);

      fprintf(fd,"line  %f -100000.0 %f 100000.0 \n",
	      pk.ypeak(),pk.ypeak());
    }
    fclose(fd);
  }
  // angles
  for(int z_peak=1;z_peak<=db.nzpeak; z_peak++) {
    fname[0]=z_peak + '0';
    strcpy(fname+1,"_p_peaks.kumac");
    printf("%s\n",fname);
    fd=fopen(fname,"w");
    for(int y_peak=1;y_peak<=db.nysievpeak; y_peak++) {
      Peak pk(z_peak, y_peak);
      pk.set(beamx, beamy, 0.0);
      fprintf(fd,"line  %f -0.5 %f 0.5 \n",
	      pk.phipeak(),pk.phipeak());
    }
    for(int x_peak=1;x_peak<=db.nxsievpeak; x_peak++){
      fprintf(fd,"line -0.5   %f 0.5  %f \n",
	      db.xsievpeak[x_peak],db.xsievpeak[x_peak]);
    }
    fclose(fd);
  } 
}