**************************************************************************/
#include "AliHLTPHOSRawAnalyzer.h"
-#include <iostream>
+//#include <iostream>
-using std::cout;
-using std::endl;
+//using std::cout;
+//using std::endl;
-
-
-AliHLTPHOSRawAnalyzer:: AliHLTPHOSRawAnalyzer():fFloatDataPtr(0), fSampleFrequency(10), fTau(2), fDTof(99999), fDAmpl(99999), n(99999)
+AliHLTPHOSRawAnalyzer:: AliHLTPHOSRawAnalyzer(): AliHLTPHOSBase(), fIntDataPtr(0), fSampleFrequency(10), fTau(2), fDTof(99999), fDAmpl(99999)
{
- fFloatDataPtr = new double[1008];
+ fIntDataPtr = new UInt_t[1008];
}
AliHLTPHOSRawAnalyzer::~AliHLTPHOSRawAnalyzer()
{
-
+ delete[] fIntDataPtr;
}
-AliHLTPHOSRawAnalyzer::AliHLTPHOSRawAnalyzer(AliHLTPHOSRawAnalyzer const&):fFloatDataPtr(0), fSampleFrequency(10), fTau(2), fDTof(99999), fDAmpl(99999), n(99999)
-{
-}
/**
* Main constructor
-* @param dataPtr Data array for wich a subarray will be taken to perform the fit
+* @param dtaPtr Data array for wich a subarray will be taken to perform the fit
* @param fs the sampling frequency in entities of MHz. Needed in order to calculate physical time
**/
-AliHLTPHOSRawAnalyzer::AliHLTPHOSRawAnalyzer(double *dtaPtr, double fs):fFloatDataPtr(0), fSampleFrequency(10), fTau(2), fDTof(99999), fDAmpl(99999), n(99999)
+AliHLTPHOSRawAnalyzer::AliHLTPHOSRawAnalyzer(double *dtaPtr, double fs): AliHLTPHOSBase(), fIntDataPtr(0), fSampleFrequency(10), fTau(2), fDTof(99999), fDAmpl(99999), fStartIndex(0)
{
- fFloatDataPtr = dtaPtr;
fSampleFrequency = fs;
} //end
void
AliHLTPHOSRawAnalyzer::BaselineCorrection(double *dataPtr, int N)
{
- fFloatDataPtr = dataPtr;
- n=N;
cout << "Baseline correction not yet implemeted" << endl;
} //end BaselineCorrection
/**
-* Shifts the basline with the amount given by baselineValue
+* Shifts the baseline with the amount given by baselineValue
* If pretrigger samples are not present then the basline correction will be incorrect.
* @param dataPtr array for wich to correct the basline
-* @param BaslineValue the basline value to subtract..
+* @param baselineValue the basline value to subtract..
**/
void
AliHLTPHOSRawAnalyzer::BaselineCorrection(double *dataPtr, double baselineValue)
{
- fFloatDataPtr = dataPtr;
printf("\nbaselineValue = %f\n", baselineValue);
cout << "Baseline correction not yet implemeted" << endl;
} //end BaslineCorrection
* Physical time is found by multiplying with the sampling intervall (Ts).
**/
float
-AliHLTPHOSRawAnalyzer::GetTiming()
+AliHLTPHOSRawAnalyzer::GetTiming() const
{
return fDTof;
} //end GetTiming
* Absolute enrgy is found by multiplying with offline calibration constants.
**/
float
-AliHLTPHOSRawAnalyzer::GetEnergy()
+AliHLTPHOSRawAnalyzer::GetEnergy() const
{
return fDAmpl;
} //end GetEnergy
+/**
+ * Set data array. Overrides data data array set in the constructor.
+ **/
+void
+AliHLTPHOSRawAnalyzer::SetData(UInt_t *data)
+{
+ fIntDataPtr = data;
+}
+
/**
* Set data array. Overrides data data array set in the constructor.
**/
void
AliHLTPHOSRawAnalyzer::SetData(double *data)
{
- // double *dta;
- // dta = data;
fFloatDataPtr = data;
- // cout << "Set data not yet implemented" << endl;
}
+
+
void
AliHLTPHOSRawAnalyzer::SetSampleFreq(double freq)
{
int
AliHLTPHOSRawAnalyzer::FindStartIndex(double treshold)
{
- printf("\ntreshold = %f \n", treshold);
cout << "Find Start index not yet implemented" << endl;
return 0;
} //end FindStartIndex
void
AliHLTPHOSRawAnalyzer::MakeInitialGuess(int treshold)
{
- printf("\ntreshold = %d\n", treshold);
cout << "Make initial guess not yet implemeted" << endl;
}
+
+
+void
+AliHLTPHOSRawAnalyzer::SetStartIndex(int index)
+{
+ fStartIndex = index;
+}
+
+
+
+void
+AliHLTPHOSRawAnalyzer::SetTVector(Double_t *tVector, Int_t size)
+{
+ cout <<"ERROR: AliHLTPHOSRawAnalyzer::SetTVector: You cannot set the peakfindervector here, must be set in derived class peakfinder"<<endl;
+}
+
+
+
+void
+AliHLTPHOSRawAnalyzer::SetAVector(Double_t *aVector, Int_t size)
+{
+ cout <<"ERROR: AliHLTPHOSRawAnalyzer::SetAVector: You cannot set the peakfindervector here, must be set in derived class peakfinder"<<endl;
+}
+
+/*
+UInt_t
+AliHLTPHOSRawAnalyzer::GetMaxValue(UInt_t *dta, Int_t size) const
+{
+
+ Double_t tmpMax = 0;
+
+ for(int i = 0; i < size; i++)
+ {
+ if(dta[i] > tmpMax)
+ {
+ tmpMax = dta[i];
+ }
+ }
+
+ return tmpMax;
+
+}
+*/