#include <stdlib.h>
#include <stdio.h>
#include <math.h>

int main(){
  double x,newet,oldet,et,newdnpdy,olddnpdy,dnpdy;
  double sigmay,oldsigmay,newsigmay;
  int i,j,column[5],oldcolumn[5],nevents,newnevents,oldnevents;
  FILE *oldfptr,*newfptr,*tempfptr;
  oldfptr=fopen("temp_data/oldscatvstime.dat","r");
  newfptr=fopen("temp_data/scatvstime.dat","r");
  tempfptr=fopen("temp_data/tempscatvstime.dat","w");
  fscanf(newfptr,"%lf %lf %lf %d",&newet,&newdnpdy,&newsigmay,&newnevents);
  fscanf(oldfptr,"%lf %lf %lf %d",&oldet,&olddnpdy,&oldsigmay,&oldnevents);
  et=newet+oldet;
  sigmay=oldsigmay+newsigmay;
  dnpdy=olddnpdy+newdnpdy;
  nevents=oldnevents+1;
  //printf("oldnevents=%d   newnevents=%d\n",oldnevents,newnevents);
  printf("average et=%g  average dnpdy=%g average sigmay=%g\n",
	 et/nevents,dnpdy/nevents,sigmay/nevents);
  fprintf(tempfptr,"%g %g %g %d\n",et,dnpdy,sigmay,nevents);
  for(i=0;i<50;i++){
    for(j=0;j<5;j++){
      fscanf(newfptr,"%d",&column[j]);
      fscanf(oldfptr,"%d",&oldcolumn[j]);
      if(j>0) column[j]=column[j]+oldcolumn[j];
    }
    fprintf(tempfptr,"%d %d %d %d %d\n",
	    column[0],column[1],column[2],column[3],column[4]);
  }
  fclose(newfptr);
  fclose(oldfptr);
  fclose(tempfptr);
  return 0;
}