+// -*- mode: c++ -*-
/**************************************************************************
* This file is property of and copyright by the Experimental Nuclear *
* Physics Group, Dep. of Physics *
#include <iostream>
using namespace std;
+#include "AliCaloConstants.h"
ClassImp( AliCaloRawAnalyzerFastFit )
-AliCaloRawAnalyzerFastFit::AliCaloRawAnalyzerFastFit() : AliCaloRawAnalyzer("Fast Fit (Alexei)", "FF")
+
+AliCaloRawAnalyzerFastFit::AliCaloRawAnalyzerFastFit() : AliCaloRawAnalyzerFitter("Fast Fit (Alexei)", "FF")
{
// Comment
-
- for(int i=0; i < 1008; i++)
- {
- fXAxis[i] = i;
- }
-
+ fAlgo= Algo::kFastFit;
}
+
AliCaloRawAnalyzerFastFit::~AliCaloRawAnalyzerFastFit()
{
if( maxf < fAmpCut || ( maxamp - ped) > fOverflowCut ) // (maxamp - ped) > fOverflowCut = Close to saturation (use low gain then)
{
- return AliCaloFitResults( maxamp, ped, AliCaloFitResults::kCrude, maxf, timebinOffset);
+ return AliCaloFitResults( maxamp, ped, Algo::kCrude, maxf, timebinOffset);
}
else if ( maxf >= fAmpCut ) // no if statement needed really; keep for readability
{
int last = 0;
int maxrev = maxampindex - bunchvector.at(index).GetStartBin();
- SelectSubarray( fReversed, bunchvector.at(index).GetLength(), maxrev , &first, &last);
+ SelectSubarray( fReversed, bunchvector.at(index).GetLength(), maxrev , &first, &last, fFitArrayCut);
+
int nsamples = last - first + 1;
if( ( nsamples ) >= fNsampleCut )
Double_t chi2 = 0;
Double_t dTau = 2.35; // time-bin units
- AliCaloFastAltroFitv0::FastFit(fXAxis, ordered , nsamples,
+ AliCaloFastAltroFitv0::FastFit(fXaxis, ordered , nsamples,
eSignal, dTau, dAmp, eAmp, dTime0, eTime, chi2);
Double_t dTimeMax = dTime0 + timebinOffset - (maxrev - first) // abs. t0
+ dTau; // +tau, makes sum tmax
- return AliCaloFitResults(maxamp, ped, AliCaloFitResults::kFitPar, dAmp, dTimeMax, timebinOffset, chi2, AliCaloFitResults::kDummy,
- AliCaloFitResults::kDummy, AliCaloFitSubarray(index, maxrev, first, last) );
+ return AliCaloFitResults(maxamp, ped, Ret::kFitPar, dAmp, dTimeMax, timebinOffset, chi2, Ret::kDummy,
+ Ret::kDummy, AliCaloFitSubarray(index, maxrev, first, last) );
} // samplecut
else
{
Float_t chi2 = CalculateChi2(maxf, maxrev, first, last);
Int_t ndf = last - first - 1; // nsamples - 2
- return AliCaloFitResults( maxamp, ped, AliCaloFitResults::kCrude, maxf, timebinOffset,
- timebinOffset, chi2, ndf, AliCaloFitResults::kDummy, AliCaloFitSubarray(index, maxrev, first, last) );
+ return AliCaloFitResults( maxamp, ped, Ret::kCrude, maxf, timebinOffset,
+ timebinOffset, chi2, ndf, Ret::kDummy, AliCaloFitSubarray(index, maxrev, first, last) );
}
} // ampcut
} // bunch index
- return AliCaloFitResults(AliCaloFitResults::kInvalid , AliCaloFitResults::kInvalid);
+ return AliCaloFitResults( Ret::kInvalid , Ret::kInvalid );
}