Index: trunk/MagicSoft/Mars/Changelog
===================================================================
--- trunk/MagicSoft/Mars/Changelog	(revision 4809)
+++ trunk/MagicSoft/Mars/Changelog	(revision 4810)
@@ -23,4 +23,5 @@
    * mcalib/MCalibrationChargeCalc.cc
    * mcalib/MCalibrationTestCalc.cc
+   * mcalib/MCalibrationRelTimeCalc.cc
    * mcalib/MHCalibrationTestTimeCam.cc
    * mcalib/MHCalibrationTestCam.cc
Index: trunk/MagicSoft/Mars/mcalib/MCalibrationTestCalc.cc
===================================================================
--- trunk/MagicSoft/Mars/mcalib/MCalibrationTestCalc.cc	(revision 4809)
+++ trunk/MagicSoft/Mars/mcalib/MCalibrationTestCalc.cc	(revision 4810)
@@ -235,22 +235,15 @@
   const UInt_t nsectors = fGeom->GetNumSectors();
 
-  Double_t lowlim      [nareas];
-  Double_t upplim      [nareas];
-  Double_t areaphotons [nareas], sectorphotons [nsectors];
-  Double_t areavars    [nareas], sectorvars    [nsectors];
-  Double_t fittedmean  [nareas], fittedsigma   [nareas];
-  Int_t   numareavalid[nareas], numsectorvalid[nsectors];
-
-  memset(lowlim        ,0, nareas   * sizeof(Double_t));
-  memset(upplim        ,0, nareas   * sizeof(Double_t));
-  memset(fittedmean    ,0, nareas   * sizeof(Double_t));
-  memset(fittedsigma   ,0, nareas   * sizeof(Double_t));
-  memset(areaphotons   ,0, nareas   * sizeof(Double_t));
-  memset(areavars      ,0, nareas   * sizeof(Double_t));
-  memset(numareavalid  ,0, nareas   * sizeof(Int_t  ));
-  memset(sectorphotons ,0, nsectors * sizeof(Double_t));
-  memset(sectorvars    ,0, nsectors * sizeof(Double_t));
-  memset(numsectorvalid,0, nsectors * sizeof(Int_t  ));
-  
+  TArrayD lowlim      (nareas);
+  TArrayD upplim      (nareas);
+  TArrayD areaphotons (nareas); 
+  TArrayD sectorphotons(nsectors);
+  TArrayD areavars    (nareas); 
+  TArrayD sectorvars  (nsectors);
+  TArrayD fittedmean  (nareas); 
+  TArrayD fittedsigma (nareas);
+  TArrayI numareavalid(nareas); 
+  TArrayI numsectorvalid(nsectors);
+
   //
   // First loop: Get mean number of calibrated photons and the RMS
@@ -380,7 +373,7 @@
   *fLog << endl;
 
-  memset(numareavalid,0,nareas*sizeof(Int_t));
-  memset(areaphotons ,0,nareas*sizeof(Double_t));
-  memset(areavars    ,0,nareas*sizeof(Double_t));
+  numareavalid.Reset();
+  areaphotons .Reset();
+  areavars    .Reset();
 
   //