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, bool top, double c=0): slsReceiverData<uint32_t>(xpixels, ypixels, np, bsize),
	xtalk(c), bufferSize(bsize), dataSize(dsize), dynamicRange(dr), numberOfPackets(np){


		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;


		if(top){
			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;
						}
					}
				}
			}
		}

/*
		else{
			iPacket1 = (totalNumberOfBytes/2) - 1040 - 8;
			iPacket2 = totalNumberOfBytes - 1040 - 8;
			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);//increment now but decrement before by 1024 before and after
						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 frame number for the given dataset.
	     \param buff pointer to the dataset
	     \returns frame number
	 */
	int getFrameNumber(char *buff){
		return(*(unsigned int*)(((eiger_packet_header *)((char*)buff))->num1));
	};


	/** 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){
		/*other way roud if its a bottom*/

#ifdef VERY_DEBUG
		cprintf(RED, "\n0x%x - %d - %d",
			(*(uint8_t*)(((eiger_packet_header *)((char*)(buff)))->num3)),
			(*(uint8_t*)(((eiger_packet_header *)((char*)(buff)))->num4)),
			(*(uint16_t*)(((eiger_packet_header *)((char*)(buff)))->num2)));
#endif
		//32 bit packet number written in num2
		if(dynamicRange == 32){
			//both ports have same packet numbers, so reconstruct
			if((*(uint8_t*)(((eiger_packet_header *)((char*)buff))->num3))){
				return ((*(uint16_t*)(((eiger_packet_header *)((char*)buff))->num2))+(numberOfPackets/2) +1);
			}else{
				return ((*(uint16_t*)(((eiger_packet_header *)((char*)buff))->num2))+1);
			}
		}
		else{
			//both ports have same packet numbers, so reconstruct
			if((*(uint8_t*)(((eiger_packet_header *)((char*)buff))->num3))){
				return ((*(uint8_t*)(((eiger_packet_header *)((char*)buff))->num4))+(numberOfPackets/2) +1);
			}else{
				return ((*(uint8_t*)(((eiger_packet_header *)((char*)buff))->num4))+1);
			}
		}
	};


	/**
    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;
	const int numberOfPackets;


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

};



#endif