diff --git a/slsDetectorCalibration/IntMap.h b/slsDetectorCalibration/IntMap.h
new file mode 100644
index 000000000..253ee95e5
--- /dev/null
+++ b/slsDetectorCalibration/IntMap.h
@@ -0,0 +1,113 @@
+#ifndef INTMAP_H
+#define INTMAP_H
+#define N 60
+#include <stdio.h>
+#include <math.h>
+
+
+
+
+class IntMap{
+
+
+ public:
+
+  IntMap(){};
+
+  ~IntMap(){};
+
+  void Init(){
+    //lookup table of the output intensity for IR laser
+    //measurements performed by Dominic april 2014
+    //intensity[59]=intensity with 0 Optical Density
+    //intensity[0]=intensity with 5.9 Optical Density
+    intensity[0]=29;//5.9
+    intensity[1]=21;//5.8
+    intensity[2]=31;//5.7
+    intensity[3]=43;//5.6
+    intensity[4]=60;//5.5
+    intensity[5]=91;//5.4
+    intensity[6]=69;//5.3
+    intensity[7]=102;//5.2
+    intensity[8]=136;//5.1
+    intensity[9]=196;//5.0
+    intensity[10]=425;//4.9
+    intensity[11]=311;//4.8
+    intensity[12]=462;//4.7
+    intensity[13]=653;//4.6
+    intensity[14]=926;//4.5
+    intensity[15]=1423;//4.4
+    intensity[16]=1072;//4.3
+    intensity[17]=1592;//4.2
+    intensity[18]=2142;//4.1
+    intensity[19]=3085;//4.0
+    intensity[20]=729;//3.9
+    intensity[21]=533;//3.8
+    intensity[22]=793;//3.7
+    intensity[23]=1121;//3.6
+    intensity[24]=1588;//3.5
+    intensity[25]=2439;//3.4
+    intensity[26]=1842;//3.3
+    intensity[27]=2730;//3.2
+    intensity[28]=3663;//3.1
+    intensity[29]=5271;//3.0
+    intensity[30]=8102;//2.9
+    intensity[31]=5933;//2.8
+    intensity[32]=8789;//2.7
+    intensity[33]=12350;//2.6
+    intensity[34]=17358;//2.5
+    intensity[35]=26300;//2.4
+    intensity[36]=20029;//2.3
+    intensity[37]=29414;//2.2
+    intensity[38]=39202;//2.1
+    intensity[39]=55724;//2.0
+    intensity[40]=15697;//1.9
+    intensity[41]=11541;//1.8
+    intensity[42]=16976;//1.7
+    intensity[43]=23866;//1.6
+    intensity[44]=33478;//1.5
+    intensity[45]=50567;//1.4
+    intensity[46]=38552;//1.3
+    intensity[47]=56394;//1.2
+    intensity[48]=74897;//1.1
+    intensity[49]=106023;//1.0
+    intensity[50]=157384;//0.9
+    intensity[51]=117677;//0.8
+    intensity[52]=171101;//0.7
+    intensity[53]=236386;//0.6
+    intensity[54]=327248;//0.5
+    intensity[55]=492781;//0.4
+    intensity[56]=379641;//0.3
+    intensity[57]=546927;//0.2
+    intensity[58]=717203;//0.1
+    intensity[59]=1000000;//0.
+    return;
+  };
+    //_od is the total Optical Density
+  int getIntensity(float _od){
+    int _int(-1);
+
+    //these lines are to take into account rounding errors with floats
+    float hun_od = 100.*_od;
+    int Ihun_od = (int)round(hun_od); 
+    float R_od =(float) Ihun_od/10.;
+    int I_od = (int)R_od;
+
+    if(I_od >-1 && I_od <60){
+      int index=59-I_od;
+      cerr<<index<<endl;
+      _int=intensity[index];
+    }else{
+      cerr<<"Optical density out of range!"<<endl;
+    }
+    return _int;
+  };
+
+
+ private:
+
+  int intensity[N];
+
+};
+
+#endif