1 /**************************************************************************
2 * Copyright(c) 2006, ALICE Experiment at CERN, All rights reserved. *
4 * Author: Per Thomas Hille for the ALICE HLT Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
16 #include "AliHLTPHOSRawAnalyzerPeakFinder.h"
22 ClassImp(AliHLTPHOSRawAnalyzerPeakFinder)
25 AliHLTPHOSRawAnalyzerPeakFinder::AliHLTPHOSRawAnalyzerPeakFinder(const AliHLTPHOSRawAnalyzerPeakFinder&):AliHLTPHOSRawAnalyzer(),tVector(0), aVector(0)
32 * The AliHLTPHOSPeakfinder class is the class for extracting the basic signal parameters
33 * "timing" and "energy" from the PHOS raw data. Physical data will for a given readout channel be
34 * a sequense of ADC digitized 10 bit integer values, however for performance reasons all values used in
35 * calculation is of type double.
37 AliHLTPHOSRawAnalyzerPeakFinder::AliHLTPHOSRawAnalyzerPeakFinder():AliHLTPHOSRawAnalyzer(),tVector(0), aVector(0)
39 // cout <<"PeakFinder:You cannot invoke the Fitter without arguments"<<endl;;
44 AliHLTPHOSRawAnalyzerPeakFinder::~AliHLTPHOSRawAnalyzerPeakFinder()
47 } //end AliHLTPHOSRawAnalyzerPeakFinder
50 AliHLTPHOSRawAnalyzerPeakFinder::SetTVector(double *tVec)
56 AliHLTPHOSRawAnalyzerPeakFinder::SetAVector(double *aVec)
63 * Extraction of timing and energy using the Peakfinde Algorithm.
64 * The. The parameters "start" and "length" defines a sub array of the data array
65 * that will be used for the the fit. If start+length must not exeed the total length
66 * of the Data array. "start" must be chosen as close as possible to t0.
67 * The baseline must also be subtracted.
68 * The length of "tVector" and "aVector" mus be equal to length.
69 * "index + length" must not exeed the length of the data array set in the constructor.
70 * @param start the start index of the subarray of the data array.
71 * @param length the number of samples to use starting from index
72 * @param tVector the peakfinder vector for timing
73 * @param aVector the peakfinder vector for amplitude (energy)
76 AliHLTPHOSRawAnalyzerPeakFinder::Evaluate(int start, int length)
78 printf("\n AliHLTPHOSRawAnalyzerPeakFinder::Evaluat from index %d to %d\n", start, start + length);
83 if(tVector == 0 || aVector == 0)
85 printf("\nError: the peakfinder vectors are not specified, aborting !!!\n");
90 printf("\nstart = %d, length = %d\n", start, length);
93 for(int i=0; i < length; i++)
95 fDAmpl += aVector[i]*fFloatDataPtr[i];
98 for(int i=0; i < length; i++)
100 tmpTime[i] = tVector[i]*fFloatDataPtr[i];
101 fDTof = fDTof + tmpTime[i];
104 fDTof = fDTof/fDAmpl;
109 } //end FitPeakFinder