HLT/PHOS/AliHLTPHOSRawAnalyzerCrude.cxx

   1 /**************************************************************************
   2  * Copyright(c) 2006, ALICE Experiment at CERN, All rights reserved.      *
   3  *                                                                        *
   4  * Author: Per Thomas Hille for the ALICE HLT Project.                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16 #include "AliHLTPHOSRawAnalyzerCrude.h"
  17
  18 //ClassImp(AliHLTPHOSRawAnalyzerCrude)
  19
  20
  21 /**
  22  * The AliHLTPHOSPeakfinder class is the class for extracting the basic signal parameters
  23  * "timing" and "energy" from the PHOS raw data. Physical data will for a given readout channel be
  24  * a sequense of ADC digitized 10 bit integer values, however for performance reasons all values used in
  25  * calculation is of type double.
  26  **/
  27 //____________________________________________________________________________
  28 AliHLTPHOSRawAnalyzerCrude::AliHLTPHOSRawAnalyzerCrude():AliHLTPHOSRawAnalyzer()
  29 {
  30
  31 }
  32
  33
  34 //____________________________________________________________________________
  35 AliHLTPHOSRawAnalyzerCrude::~AliHLTPHOSRawAnalyzerCrude()
  36 {
  37
  38 } //end AliHLTPHOSRawAnalyzerCrude
  39
  40
  41 /**
  42 * Extraction of timing and energy using Crude estimate.
  43 * The. The parameters "start" and "length" defines a sub array  of the data array
  44 * that will be used for the the fit. If start+length must not exeed the total length
  45 * of the Data array. "start" must be chosen as close as possible to t0.
  46 * The baseline must also be subtracted.                                                              .
  47 * "index + length" must not exeed the length of the data array set in the constructor.
  48 * @param start the start index of the subarray of the data array.
  49 * @param length the number of samples to use starting from index
  50 **/
  51 // //____________________________________________________________________________
  52 // void
  53 // AliHLTPHOSRawAnalyzerCrude::Evaluate(int start, int length)
  54 // {
  55 //   double tmpAmplitudeMax =0;
  56 //   double tmpTime = 0;
  57
  58 //   for(int i=start; i<length; i++)
  59 //     {
  60 //       if(fIntDataPtr[i] >  tmpAmplitudeMax && i > 5)
  61 //      {
  62 //        tmpAmplitudeMax = fIntDataPtr[i];
  63 //        tmpTime = i;
  64 //      }
  65 //     }
  66
  67 //   fDAmpl = tmpAmplitudeMax;
  68 //   fDTof =  tmpTime;
  69
  70 //   //thats all
  71 // } //end Crude
  72
  73
  74 void
  75 AliHLTPHOSRawAnalyzerCrude::Evaluate(int start, int length)
  76 {
  77   //  cout << "AliHLTPHOSRawAnalyzerCrude::Evaluate TP0"  << endl;
  78
  79   //DumpData(T
  80
  81   //  DumpData(fDoubleDataPtr,50, 25);
  82
  83   double tmpAmplitudeMax =0;
  84   double tmpTime = 0;
  85
  86   for(int i=start; i<length; i++)
  87     {
  88         //   if(fDoubleDataPtr[i] >  tmpAmplitudeMax && i > 5)
  89       if(fIntDataPtr[i] >  tmpAmplitudeMax && i > 5)
  90         {
  91           //      tmpAmplitudeMax = fDoubleDataPtr[i];
  92           tmpAmplitudeMax = fIntDataPtr[i];
  93           tmpTime = i;
  94         }
  95     }
  96
  97
  98   fDAmpl = tmpAmplitudeMax;
  99   fDTof =  tmpTime;
 100
 101   //thats all
 102 } //end Crude