+// $Id$
+
/**************************************************************************
* This file is property of and copyright by the Experimental Nuclear *
* Physics Group, Dep. of Physics *
* provided "as is" without express or implied warranty. *
**************************************************************************/
-
#include "AliHLTPHOSRawAnalyzerPeakFinder.h"
-//#include <iostream>
#include <cmath>
+#include "AliHLTCaloUtilities.h"
-
-using std::cout;
-using std::endl;
+//using std::cout;
+//using std::endl;
ClassImp(AliHLTPHOSRawAnalyzerPeakFinder)
-
/**
* The AliHLTPHOSPeakfinder class is the class for extracting the basic signal parameters
* "timing" and "energy" from the PHOS raw data. Physical data will for a given readout channel be
* a sequense of ADC digitized 10 bit integer values, however for performance reasons all values used in
* calculation is of type double.
**/
-AliHLTPHOSRawAnalyzerPeakFinder::AliHLTPHOSRawAnalyzerPeakFinder():AliHLTPHOSRawAnalyzer(), fTVectorPtr(0), fAVectorPtr(0), fTVectorSize(0), fAVectorSize(0)
+AliHLTPHOSRawAnalyzerPeakFinder::AliHLTPHOSRawAnalyzerPeakFinder():AliHLTPHOSRawAnalyzer(),
+ fTVectorPtr(0),
+ fAVectorPtr(0),
+ fTVectorSize(0),
+ fAVectorSize(0)
{
- // cout <<"PeakFinder:You cannot invoke the Fitter without arguments"<<endl;;
+
}
} //end AliHLTPHOSRawAnalyzerPeakFinder
-
void
AliHLTPHOSRawAnalyzerPeakFinder::SetTVector(Double_t *tVec, Int_t size)
{
}
-//___________________________________________________________________
void
AliHLTPHOSRawAnalyzerPeakFinder::SetAVector(Double_t *aVec, Int_t size)
{
-
+ //comment
fAVectorSize = size;
if(fAVectorPtr != 0)
}
}
-
-
-//___________________________________________________________________
void
AliHLTPHOSRawAnalyzerPeakFinder::Evaluate(Int_t /*start*/, Int_t length)
{
+ //comment
fDTof = 0;
fDAmpl = 0;
Int_t tmpLength;
for(int i=0; i < tmpLength; i++)
{
- //fDAmpl += fAVectorPtr[i]*fIntDataPtr[i]; removed 18 april 2008
fDAmpl += fAVectorPtr[i]*fDoubleDataPtr[i];
}
for(int i=0; i < tmpLength; i++)
{
- // fDTof += fTVectorPtr[i]*fIntDataPtr[i]; removed 18 april 2008
fDTof += fTVectorPtr[i]*fDoubleDataPtr[i];
}
if(fDAmpl > 900)
{
- // Double_t tmpMax = MaxValue(const_cast<unsigned int*>(fIntDataPtr), tmpLength); removed 18 april 2008
- double tmpMax = MaxValue(fDoubleDataPtr, tmpLength);
+ double tmpMax = AliHLTCaloUtilities::MaxValue(fDoubleDataPtr, tmpLength);
if(tmpMax == 1023)
{
fDAmpl = tmpMax;
}
}
-
fDTof = fDTof/fDAmpl;
-
}
} //end Evaluate