From 1f6c5c26b4d40be7d423a72f5e350d8ff018946d Mon Sep 17 00:00:00 2001
From: Dhanya Maliakal <dhanya.maliakal@psi.ch>
Date: Tue, 14 Apr 2015 17:13:09 +0200
Subject: [PATCH] eiger mapping for 32 bit done

---
 slsDetectorCalibration/eigerHalfModuleData.h | 112 ++++++++++++++-----
 slsDetectorCalibration/slsDetectorData.h     |  10 +-
 slsDetectorCalibration/slsReceiverData.h     |   2 +-
 3 files changed, 87 insertions(+), 37 deletions(-)

diff --git a/slsDetectorCalibration/eigerHalfModuleData.h b/slsDetectorCalibration/eigerHalfModuleData.h
index 1c878af8c..553e9dbd3 100644
--- a/slsDetectorCalibration/eigerHalfModuleData.h
+++ b/slsDetectorCalibration/eigerHalfModuleData.h
@@ -19,8 +19,8 @@ public:
 	 */
 
 
-	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){
+	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;
@@ -42,29 +42,60 @@ public:
 		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;
+
+		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;
+						}
+					}
+				}
+			}
+		}
+*/
 
 
 
@@ -87,7 +118,7 @@ public:
 	     \returns frame number
 	 */
 	int getFrameNumber(char *buff){
-		return htonl(*(unsigned int*)((eiger_image_header *)((char*)(buff)))->fnum);
+		return(*(unsigned int*)(((eiger_packet_header *)((char*)buff))->num1));
 	};
 
 
@@ -96,7 +127,31 @@ public:
      	 \returns packet number
 	 */
 	int getPacketNumber(char *buff){
-		return htonl(*(unsigned int*)((eiger_packet_header *)((char*)(buff)))->num2);
+		/*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);
+			}
+		}
 	};
 
 
@@ -146,21 +201,16 @@ private:
 	const int bufferSize;
 	const int dataSize;
 	const int dynamicRange;
+	const int numberOfPackets;
 
 
-	/** 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];
+		unsigned char num2[2];
+		unsigned char num3[1];
+		unsigned char num4[1];
 	} eiger_packet_header;
 
 };
diff --git a/slsDetectorCalibration/slsDetectorData.h b/slsDetectorCalibration/slsDetectorData.h
index 04ddc43d5..cb364bc04 100644
--- a/slsDetectorCalibration/slsDetectorData.h
+++ b/slsDetectorCalibration/slsDetectorData.h
@@ -211,14 +211,14 @@ class slsDetectorData {
 
   */
   virtual dataType getChannel(char *data, int ix, int iy, int dr) {
-	  dataType m=0, d=0;
+	  dataType m=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];
-
+		  //cout <<"ix:"<<dec<<ix<<" iy:"<<iy<<endl;
 		  numBytes = (nx * iy + ix);
 		  divFactor=2;
 		  if(dr == 4) divFactor = 16;
@@ -228,10 +228,10 @@ class slsDetectorData {
 		  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;
+		  //cout <<"pixelval:"<<dec <<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;
+		 // cout<<"t:"<<hex<<t<<endl;
 
 	  }else
 		  cprintf(RED,"outside limits\n");
diff --git a/slsDetectorCalibration/slsReceiverData.h b/slsDetectorCalibration/slsReceiverData.h
index 4781cb36e..15d324ef6 100644
--- a/slsDetectorCalibration/slsReceiverData.h
+++ b/slsDetectorCalibration/slsReceiverData.h
@@ -67,7 +67,7 @@ public:
 		while (dd<=(dsize-packetSize)) {
 			pnum=getPacketNumber(p);
 			fn=getFrameNumber(p);
-
+			//cout <<"pnum:"<<pnum<<" fn:"<<fn<<"\t";
 
 			if (pnum<1 || pnum>nPackets) {
 				cout << "Bad packet number " << pnum << " frame "<< fn << endl;