Index: branches/Mars_MC/macros/fact/fake_drs_calibration.C
===================================================================
--- branches/Mars_MC/macros/fact/fake_drs_calibration.C	(revision 17051)
+++ branches/Mars_MC/macros/fact/fake_drs_calibration.C	(revision 17051)
@@ -0,0 +1,58 @@
+#include "DrsCalib.h"
+
+
+using namespace std;
+
+// --------------------------------------------------------------------------
+
+int fake_drs_calibration()
+{
+    DrsCalibration drscalib;
+    const string filename = "/fhgfs/groups/app/testingdrsfile/test150samples_2.drs.fits";
+    uint16_t fRoi = 150;
+    uint16_t fNumTm = 0;
+    float baselineMean = 3280 * 2000.0/4096.0;
+    float baselineRms = 0;
+    float gainMean = 1907.35 * 16;
+    float gainRms = 0;
+    float triggerOffsetMean = 0;
+    float triggerOffsetRms = 0;
+    float triggerOffsetTmMean = 0;
+    float triggerOffsetTmRms = 0;
+
+    const size_t n = 1440*1024*4 + 1440*fRoi*2 + fNumTm*fRoi*2 + 3;
+    vector<float> vec;
+    vec.resize(n);
+
+    reinterpret_cast<uint32_t*>(vec.data())[0] = fRoi;
+    reinterpret_cast<uint32_t*>(vec.data())[1] = 1;
+    reinterpret_cast<uint32_t*>(vec.data())[2] = 2;
+    reinterpret_cast<uint32_t*>(vec.data())[3] = 3;
+
+    for (int i = 0 ; i < 1024*1440 ; i++)
+    {
+        vec.data()[i+4]             = baselineMean;
+        vec.data()[i+4+1*1024*1440] = baselineRms;
+        vec.data()[i+4+2*1024*1440] = gainMean;
+        vec.data()[i+4+3*1024*1440] = gainRms;
+    }
+    for (int i = 0 ; i < fRoi*1440 ; i++)
+    {
+        vec.data()[i+4+4*1024*1440]             = triggerOffsetMean;
+        vec.data()[i+4+4*1024*1440+1*fRoi*1440] = triggerOffsetRms;
+    }
+    for (int i = 0 ; i < fRoi*fNumTm ; i++)
+    {
+        vec.data()[i+4+2*1024*1440+2*fRoi*1440]             = triggerOffsetTmMean;
+        vec.data()[i+4+3*1024*1440+2*fRoi*1440+fRoi*fNumTm] = triggerOffsetTmRms;
+    }
+
+    drscalib.fRoi   = fRoi;
+    drscalib.fNumTm = fNumTm;
+    drscalib.fStep  = 2;
+
+    drscalib.WriteFitsImp(filename,vec);
+
+    return 0;
+}
+