HLT/PHOS/AliHLTPHOSRawAnalyzerCrude.cxx

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