slsDetectorPackage/slsDetectorCalibration/eigerHalfModuleData.h

#ifndef EIGERMODULEDATA_H
#define  EIGERMODULEDATA_H
#include "slsReceiverData.h"



class eigerHalfModuleData : public slsReceiverData<uint32_t> {
public:




	/**
     Implements the slsReceiverData structure for the eiger prototype read out by a half module i.e. using the slsReceiver
     (256*256 pixels, 512 packets for 16 bit mode, 256 for 8, 128 for 4, 1024 for 32, 1040  etc.)
     \param d dynamic range
     \param c crosstalk parameter for the output buffer

	 */


	eigerHalfModuleData(int dr, int np, int bsize, int dsize, double c=0): slsReceiverData<uint32_t>(xpixels, ypixels, np, bsize),
	xtalk(c), dynamicRange(dr), bufferSize(bsize), dataSize(dsize){


		int **dMap;
		uint32_t **dMask;

		dMap=new int*[ypixels];
		dMask=new uint32_t*[ypixels];


		for (int i = 0; i < ypixels; i++) {
			dMap[i] = new int[xpixels];
			dMask[i] = new uint32_t[xpixels];
		}

		//Map
		int totalNumberOfBytes = 1040 * dynamicRange * 16 *2; //for both 1g and 10g
		int iPacket1 = 8;
		int iPacket2 = (totalNumberOfBytes/2) + 8;
		int iData1 = 0, iData2 = 0;
		int iPort;

		for (int ir=0; ir<ypixels; ir++) {
			for (int ic=0; ic<xpixels; ic++) {
				iPort = ic / (xpixels/2);
				if(!iPort){
					dMap[ir][ic]  = iPacket1;
					iPacket1 += (dynamicRange / 8);
					iData1 +=(dynamicRange / 8);
					if(iData1 >= dataSize){
						iPacket1 += 16;
						iData1 = 0;
					}
				}else{
					dMap[ir][ic]  = iPacket2;
					iPacket2 += (dynamicRange / 8);
					iData2 +=(dynamicRange / 8);
					if(iData2 >= dataSize){
						iPacket2 += 16;
						iData2 = 0;
					}
				}
			}
		}




		//Mask
		for(int ir=0; ir<ypixels; ++ir)
			for(int ic=0; ic<xpixels; ++ic)
				dMask[ir][ic] = 0x0;

		setDataMap(dMap);
		setDataMask(dMask);




	};

  /**

     Returns the value of the selected channel for the given dataset. Virtual function, can be overloaded.
     \param data pointer to the dataset (including headers etc)
     \param ix pixel number in the x direction
     \param iy pixel number in the y direction
     \param dr dynamic range
     \returns data for the selected channel, with inversion if required

  */
  virtual dataType getChannel(char *data, int ix, int iy=0, int dr=0) {
	  dataType m=0, d=0;
	  uint64_t t;
	  int numBytes,divFactor,newix,pixelval;


	  if (ix>=0 && ix<nx && iy>=0 && iy<ny && dataMap[iy][ix]>=0 && dataMap[iy][ix]<dataSize) {
		  m=dataMask[iy][ix];

		  numBytes = (nx * iy + ix);
		  divFactor=2;
		  if(dr == 4) divFactor = 16;
		  else if (dr == 8) divFactor = 8;
		  else if (dr == 16) divFactor = 4;

		  pixelval = numBytes % divFactor;
		  newix = ix - pixelval;

		  //cout <<"pixelval:"<<pixelval<<" newix:"<<newix<<endl;
		  //cout <<"64:"<< hex<<((uint64_t)(*((uint64_t*)(((char*)data)+(dataMap[iy][newix])))))<<endl;
		  t = (be64toh((uint64_t)(*((uint64_t*)(((char*)data)+(dataMap[iy][newix]))))));
		  //cout<<"t:"<<t<<endl;

	  }else
		  cprintf(RED,"outside limits\n");

	  if(dr == 4)
		  //uint8_t value =  t >> (pixelval*4); cout <<"value:"<< value << endl;
		  return ((t >> (pixelval*4)) & 0xf)^m;
	  else if(dr == 8)
		  //uint8_t value =  t >> (pixelval*8); cout <<"value:"<< value << endl;
		  return ((t >> (pixelval*8)) & 0xff)^m;
	  else if(dr == 16){
		  //uint16_t value =  t >> (pixelval*16); cout <<"value:"<< value << endl;
		  return ((t >> (pixelval*16)) & 0xffff)^m;
	  }else{
		  //uint32_t value =  t >> (pixelval*32); cout <<"value:"<< value << endl;
		  return ((t >> (pixelval*32)) & 0xffffffff)^m;
	  }
  };


	/** Returns the frame number for the given dataset.
	     \param buff pointer to the dataset
	     \returns frame number
	 */
	int getFrameNumber(char *buff){
		return htonl(*(unsigned int*)((eiger_image_header *)((char*)(buff)))->fnum);
	};


	





	/** gets the packets number (last packet is labelled with 0 and is replaced with 40)
     	 \param buff pointer to the memory
     	 \returns packet number
	 */
	int getPacketNumber(char *buff){
		return htonl(*(unsigned int*)((eiger_packet_header *)((char*)(buff)))->num2);
	};


	/**
    returns the pixel value as double correcting for the output buffer crosstalk
     \param data pointer to the memory
     \param ix coordinate in the x direction
     \param iy coordinate in the y direction
     \returns channel value as double

	 */
	double getValue(char *data, int ix, int iy=0) {
		//  cout << "##" << (void*)data << " " << ix << " " <<iy << endl;
		if (xtalk==0)
			  return slsDetectorData<uint32_t>::getValue(data, ix, iy);
		else
			return slsDetectorData<uint32_t>::getValue(data, ix, iy)-xtalk*slsDetectorData<uint32_t>::getValue(data, ix-1, iy);
	};



	/** sets the output buffer crosstalk correction parameter
      \param c output buffer crosstalk correction parameter to be set
      \returns current value for the output buffer crosstalk correction parameter

	 */
	double setXTalk(double c) {xtalk=c; return xtalk;}


	/** gets the output buffer crosstalk parameter
      \returns current value for the output buffer crosstalk correction parameter
	 */
	double getXTalk() {return xtalk;}







private:

	double xtalk; /**<output buffer crosstalk correction parameter */

	const static int xpixels = 1024;
	const static int ypixels = 256;
	const int bufferSize;
	const int dataSize;
	const int dynamicRange;


	/** structure of an eiger image header*/
	typedef struct
	{
		unsigned char header_before[20];
		unsigned char  fnum[4];
		unsigned char  header_after[24];
	} eiger_image_header;

	/** structure of an eiger image header*/
	typedef struct
	{
		unsigned char num1[4];
		unsigned char num2[4];
	} eiger_packet_header;

};



#endif