a trial version of eiger mapping data

slsDetectorCalibration/eigerHalfModuleData.h

@@ -0,0 +1,172 @@
+#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 np, int dr, int t, int bsize, int dsize, double c=0): slsReceiverData<uint32_t>(xpixels, ypixels, nPackets, bsize),
+	xtalk(c), nPackets(np), dynamicRange(dr), bufferSize(bsize), dataSize(dsize){
+
+
+		int **dMap;
+		uint32_t **dMask;
+		int ix, iy;
+
+		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(ix=0; ix<ypixels; ++ix)
+			for(iy=0; iy<xpixels; ++iy)
+				dMask[ix][iy] = 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 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 nPackets;
+	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

slsDetectorCalibration/slsDetectorData.h

@@ -197,6 +197,32 @@ class slsDetectorData {
     return d^m;
   };
 
+
+  /**
+
+     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, int dr) {
+    dataType m=0, d=0;
+    uint64_t* s;
+    if (ix>=0 && ix<nx && iy>=0 && iy<ny && dataMap[iy][ix]>=0 && dataMap[iy][ix]<dataSize) {
+      m=dataMask[iy][ix];
+
+      int k = (((ix * iy + iy) * (dr / 8)) % 8);//which byte (1 to 8)
+      int index = (8 - k) * 8;
+      s = be64toh(((uint64_t)(*(((uint64_t*)data)+(dataMap[iy][ix]-index)))));
+
+      d=*((dataType*)(s+(k*8)));
+    }
+    return d^m;
+  };
+
   /**
 
      Returns the value of the selected channel for the given dataset as double.
@@ -208,6 +234,18 @@ class slsDetectorData {
   */
   virtual double getValue(char *data, int ix, int iy=0) {return (double)getChannel(data, ix, iy);};
 
+  /**
+
+     Returns the value of the selected channel for the given dataset as double.
+     \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 as double
+
+  */
+  virtual double getValue(char *data, int ix, int iy, int dr) {return (double)getChannel(data, ix, iy, dr);};
+
    /**
 
      Returns the frame number for the given dataset. Purely virtual func.