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

/* Event structure:

       4 bytes event number
      26 bytes DUT temperatures
  5x 658 bytes motherboard data block:

       6 bytes MAC header
       2 bytes Status register
       4 bytes Event number
       4 bytes TDC value
       320 bytes ADC values ASIC 0
       320 bytes ADC values ASIC 1
       2 bytes NTC reading 
*/

#define ADlength  13

#define Nstations  5
#define MBlength 658


#define IPlength   6
#define SRlength   2
#define ENlength   4
#define TMlength   4
#define A0length 160
#define A1length 160
#define TTlength   2

#define BLlength (4 + 2*ADlength + Nstations * MBlength)

  unsigned Offset[] = {
127.94,
147.20,
153.55,
106.16,
139.55,
150.16,
141.81,
158.51,
122.72,
131.20,
160.81,
163.63,
0.
};


  unsigned char *buffer;
  unsigned char *Value;

int main(int argc, char **argv)
{
  FILE *fp;

  int ind;
  char *Fname;
  long lSize;

  size_t result;

  int Nblocks;

  if (argc!=2) {fputs ("enter file name\n", stderr); exit(0);}

  Fname = argv[1];

  fp = fopen (Fname, "rb");
  if (fp==NULL) {fputs ("File error", stderr); exit (1);}

  /* obtain file size */

  fseek (fp , 0 , SEEK_END);
  lSize = ftell (fp);
  rewind (fp);

  /* allocate memory to contain one event */

  buffer = (unsigned char*) malloc(sizeof(unsigned char)*BLlength);

  if (buffer==NULL) {fputs ("Memory error", stderr); exit (2);}

/* Loop over the entire data sample */

  Nblocks = lSize/BLlength;  // integer division to cut off corrupted tail !!!

  for (ind=0; ind<Nblocks; ind++)
    {
      int i, MBA;

      int Gevt[4], Gevent;
      int Tmp[ADlength], Temp[ADlength];

/* copy the data portion into the buffer */

     result = fread (buffer, 1, BLlength, fp);
     if (result != BLlength) {fputs ("Reading error", stderr); exit (3);}

/* set pointer to the beginning of data field */

     Value = buffer;

/* Global event number - needed when the local event number
is not available (motherboard readout causes time out) */

     Gevt[0] = (int)*Value++;
     Gevt[1] = (int)*Value++;     
     Gevt[2] = (int)*Value++;     
     Gevt[3] = (int)*Value++;

     Gevent = Gevt[3]<<24 | Gevt[2]<<16 | Gevt[1]<<8 | Gevt[0];

/* DUT temperatures */

     for (i=0; i<ADlength; i++)
        {
          int ADRESH;
          int ADRESL;
          int ADRESM;

          ADRESH = (int)*Value++;
          ADRESL = (int)*Value++;

          ADRESM = (ADRESH & 0x1C)<<3;
          ADRESH = (ADRESH & 0x03)<<8;

          Tmp[i] = ADRESH | ADRESM | ADRESL;
        }

/* Mapping temperature channels and pedestal subtraction */

     Temp[0]  = Tmp[7]  - Tmp[2];
     Temp[1]  = Tmp[6]  - Tmp[2];
     Temp[2]  = Tmp[5]  - Tmp[2];
     Temp[3]  = Tmp[4]  - Tmp[2];
     Temp[4]  = Tmp[3]  - Tmp[2];
     Temp[5]  = Tmp[1]  - Tmp[2];
     Temp[6]  = Tmp[0]  - Tmp[2];
     Temp[7]  = Tmp[11] - Tmp[2];
     Temp[8]  = Tmp[9]  - Tmp[2];
     Temp[9]  = Tmp[8]  - Tmp[2];
     Temp[10] = Tmp[10] - Tmp[2];
     Temp[11] = Tmp[12] - Tmp[2];
     Temp[12] = Tmp[2];

// for (i=0; i<ADlength-1; i++) printf("%8.3f ", ((float)Temp[i]-Offset[i])/(-5.94));
//                              printf("\n");

/* Retrieve frontend data */

     for (MBA=0; MBA<Nstations; MBA++)
        {
          int i;

          int MAC[IPlength];

          int Status, Phase;
          int Levt[4], Levent;

          int TDC[4], Time, Tcor;

          int ASIC0[A0length];
          int ASIC1[A1length];

          int NTC;


/* MAC address with the motherboard number and timeout flags */

          for (i=0; i<IPlength; i++) MAC[i] = *Value++;

/* Status register of each motherboard */

          Status = (int)*Value++;
          Phase  = (int)*Value++;

          Tcor = (Status>>1)&0x1;

//             printf(":%02X ", Status);
//             printf(":%02X ", Phase);
// if (MBA==4) printf("\n");

/* Local event number */

          Levt[3] = (int)*Value++;
          Levt[2] = (int)*Value++;     
          Levt[1] = (int)*Value++;     
          Levt[0] = (int)*Value++;

          Levent = Levt[3]<<24 | Levt[2]<<16 | Levt[1]<<8 | Levt[0];

 if (MBA==0) printf("%d ", Gevent+1);
             printf("%d ", Levent);
 if (MBA==4) printf("\n");


/* TDC values */

          TDC[3] = (int)*Value++;
          TDC[2] = (int)*Value++;
          TDC[1] = (int)*Value++;
          TDC[0] = (int)*Value++;

          Time = 2*(TDC[0] + TDC[1] + TDC[2] + TDC[3]) + Tcor;

//             printf("%d ", Time);
// if (MBA==4) printf("\n");


/* Amplitudes */

          for (i=0; i<A0length; i++)
             {
	       ASIC0[i] = (int)*Value++;
	       ASIC0[i] = ASIC0[i]<<8 | (int)*Value++;
             }

          for (i=0; i<A1length; i++)
             {
	       ASIC1[i] = (int)*Value++;
	       ASIC1[i] = ASIC1[i]<<8 | (int)*Value++;
             }

// if (Levent==1 && MBA==2) for(i=0; i<A0length; i++) printf("%d %d\n", i, ASIC0[i]);
// if (Levent==1 && MBA==2) for(i=0; i<A1length; i++) printf("%d %d\n", i, ASIC1[i]);

// if (MBA==2) printf("%d %d\n", Time, ASIC0[35]);

/* NTC temperatures */

          NTC = (int)*Value++;
          NTC = NTC<<8 | (int)*Value++;

//             printf("%d ", NTC);
// if (MBA==4) printf("\n");


        } /* end for(MBA) */
    } /* end for(Nblocks) */

  fclose (fp);
  free (buffer);
  return 0;
}