HLT/PHOS/AliHLTPHOSRawAnalyzerLMS.h

   1 //-*- Mode: C++ -*-
   2 // $Id$
   3
   4
   5 #ifndef ALIHLTPHOSRAWANALYZERLMS_H
   6 #define ALIHLTPHOSRAWANALYZERLMS_H
   7 #include <Rtypes.h>
   8 #include "TObject.h"
   9 #include "AliHLTPHOSRawAnalyzer.h"
  10
  11 /* Copyright(c) 2006, ALICE Experiment at CERN, All rights reserved. *
  12  * See cxx source for full Copyright notice                          */
  13
  14 class AliHLTPHOSRawAnalyzerLMS : public AliHLTPHOSRawAnalyzer
  15 {
  16  public:
  17   AliHLTPHOSRawAnalyzerLMS();
  18
  19   /**
  20    * Main constructor
  21    * @param dataPtr Data array for wich a subarray will be taken to perform the fit
  22    * @param fs the sampling frequency in entities of MHz. Needed in order to calculate physical time
  23    **/
  24   AliHLTPHOSRawAnalyzerLMS(double *dataPtr, double fs);
  25   AliHLTPHOSRawAnalyzerLMS(const AliHLTPHOSRawAnalyzerLMS & );
  26   AliHLTPHOSRawAnalyzerLMS & operator = (const AliHLTPHOSRawAnalyzerLMS)
  27     {
  28       return *this;
  29     }
  30
  31   virtual ~AliHLTPHOSRawAnalyzerLMS();
  32
  33
  34   /**
  35    * Extraction of timing and energy using the LMS method.
  36    * The. The parameters "start" and "length" defines a sub array  of the data array
  37    * that will be used for the the fit. If start+length must not exeed the total length
  38    * of the Data array. "start" must be chosen as close as possible to t0.
  39    * The baseline must also be subtracted.
  40    * The length of "tVector" and "aVector" mus be equal to length.
  41    * "index + length" must not exeed the length of the data array set in the constructor.
  42    * @param start the start index of the subarray of the data array.
  43    * @param length the number of samples to use starting from index
  44    **/
  45
  46   virtual void Evaluate(int start = 0, int lenght = 100);
  47
  48  private:
  49   double   fMCovarPtrPtr[1008][1008]; /**<Covariance matrix of the measurements*/
  50   double   fPCovarPtrPtr[1008][1008];   /**<Covariance matrix of the estimated parameters*/
  51   ClassDef(AliHLTPHOSRawAnalyzerLMS, 2)
  52
  53     };
  54
  55 #endif