/*****************************************************************
Function to read AC physics data v1.1
Author: Petter Hofverberg, petter@particle.kth.se

0408 v1.0 alive
0409 Added shifting operators to fix problem with odd data (BBAC|11DD)
                                                              -----
Errors: (returned to the main program as the variable "err")
err=
0 - Ok
1 - Data found, but corrupt
2 - Data dont found in buffer

ToDo: i)  do it endian-independent
      ii) check for more events in the same file?
******************************************************************/

#include <stdio.h>
#include "AC.h"

int i,err,found,j;
//unsigned short buffer[1000];
//unsigned short tmp1,tmp2;
unsigned short buffer[1000];
unsigned short tmp1;

struct physicsstruct physicsdata;
struct physicsstruct *physicspointer;

struct calibstruct calibdata;
struct calibstruct *calibpointer;


//int ACphysics(int length,unsigned short* datapointer,struct physicsstruct* physicspointer)
int ACphysics(int length,unsigned char* datapointer,struct physicsstruct* physicspointer)
{

 int iter=2;

  *physicspointer=physicsdata; 

  //look for header AC11, then cut out the event to a temp buffer
  while(found==0 && iter<length)
    {
      tmp1 = ((*datapointer) << 8) | *(datapointer+1);
      
      if(tmp1==0xAC11)
      //if(tmp1==0xAC11) //flipped
	{ 
	  found=1;
	  for(i=0;i<64;i++)
	    {
	      buffer[i]=((*datapointer) << 8) | *(datapointer+1);
	      datapointer=datapointer+2;
	    }
	  tmp1=0x0000;
	}
      else{
	datapointer++;
	iter=iter+1;}
    }
  
  //check errors
  if(found==0)
    err=2;
  else if(buffer[1]==0x22AC || buffer[1]==0xAC22) 
    err=0;
  else
    err=1;

  //Fill physicsdata from buffer
  for(i=0;i<2;i++) physicsdata.header[i]=buffer[i];
  physicsdata.status=buffer[2];
  physicsdata.hitmap=buffer[3];
  for(i=0;i<6;i++) physicsdata.regist[i]=buffer[4+i];
  for(i=0;i<16;i++) physicsdata.shift[i]=buffer[10+i];
  for(i=0;i<16;i++) physicsdata.counters[i]=buffer[26+i];
  for(i=0;i<8;i++) physicsdata.coinc[i]=buffer[42+i];
  physicsdata.trigg=buffer[50];
  for(i=0;i<2;i++) physicsdata.clock[i]=buffer[51+i];
  for(i=0;i<2;i++) physicsdata.temp[i]=buffer[53+i];
  for(i=0;i<8;i++) physicsdata.DAC[i]=buffer[55+i];
  physicsdata.CRC=buffer[63];

  return err;

}


int ACcalib(int length, unsigned char* datapointer, struct calibstruct* calibpointer)
{

  int iter=2;
	    
  //point struct-pointer to calibdata
  *calibpointer=calibdata;  

  //look for header 0xACCA
  while(found==0 && iter<length)
    {
      tmp1 = ((*datapointer) << 8) | *(datapointer+1);
      
      if(tmp1==0xACCA)
	{
	  found=1;
	  for(i=0;i<8256;i++)
	    {
	      buffer[i]=((*datapointer) << 8) | *(datapointer+1);
	      datapointer=datapointer+2;
	    }
	  tmp1=0x0000;
	}
      else{
      datapointer++;
      iter=iter+1;}
    }

  //check errors
  if(found==0)
    err=2;
  else if(buffer[8254]==0xAC22)
    err=0;
  else
    err=1;

  printf("length: %d\n iter: %d\n",length,iter);

  //Fill calibdata from buffer
  calibdata.header=buffer[0];
  for(i=0;i<5;i++) calibdata.status[i]=buffer[1+i];
  for(i=0;i<8;i++) calibdata.temp[i]=buffer[6+i];
  for(i=0;i<8;i++) calibdata.DAC1[i]=buffer[14+i];
  for(i=0;i<8;i++) calibdata.DAC2[i]=buffer[22+i];
  for(i=0;i<6;i++) calibdata.regist[i]=buffer[30+i];
  for(i=0;i<8;i++) calibdata.time[i]=buffer[36+i];
  calibdata.n_tr=buffer[44];
  for(i=0;i<16;i++) calibdata.hitmap_tr[i]=buffer[45+i];
  for(i=0;i<4096;i++) calibdata.curve1[i]=buffer[61+i];
  for(i=0;i<4096;i++) calibdata.curve2[i]=buffer[4157+i];
  calibdata.iCRC=buffer[8253];
  calibdata.tail=buffer[8254];
  calibdata.CRC=buffer[8255];

  /*  
  //debug
  if(err!=2)
    {
      printf("header.....   %hx\n",calibdata.header);
      printf("status.....   %hx\n",calibdata.status[0]);
      printf("temp.......   %hx\n",calibdata.temp[0]);
      printf("DAC1.......   %hx\n",calibdata.DAC1[0]);
      printf("DAC2.......   %hx\n",calibdata.DAC2[0]);
      printf("regist.....   %hx\n",calibdata.regist[0]);
      printf("time.......   %hx\n",calibdata.time[0]);
      printf("n_tr.......   %hx\n",calibdata.n_tr);
      printf("hitmap_tr..   %hx\n",calibdata.hitmap_tr[0]);
      printf("curve1.....   %hx\n",calibdata.curve1[0]);
      printf("curve2.....   %hx\n",calibdata.curve2[0]);
      printf("iCRC.......   %hx\n",calibdata.iCRC);
      printf("tail.......   %hx\n",calibdata.tail);
      printf("CRC........   %hx\n",calibdata.CRC);
    }
  printf("err: %d\n",err);
  */

  return err;
}



//transposing function
int flip(int i16)
{
  int temp;
  
  temp = 0;
  temp |= 0xFF00 & ((0x00FF & i16) << 8);
  temp |= 0x00FF & ((0xFF00 & i16) >> 8);
  return temp;
}