* provided "as is" without express or implied warranty. *
**************************************************************************/
-/* $Id: $ */
+/* $Id$ */
// This class extracts the signal parameters (energy, time, quality)
// from ALTRO samples. Energy is in ADC counts, time is in time bin units.
// of bunches in the signal.
//
// AliPHOSRawFitterv0 *fitterv0=new AliPHOSRawFitterv0();
-// fitterv0->SetSamples(sig,sigStart,sigLength);
-// fitterv0->SetNBunches(nBunches);
// fitterv0->SetChannelGeo(module,cellX,cellZ,caloFlag);
// fitterv0->SetCalibData(fgCalibData) ;
-// fitterv0->Eval();
+// fitterv0->Eval(sig,sigStart,sigLength);
// Double_t amplitude = fitterv0.GetEnergy();
// Double_t time = fitterv0.GetTime();
// Bool_t isLowGain = fitterv0.GetCaloFlag()==0;
//-----------------------------------------------------------------------------
AliPHOSRawFitterv0::AliPHOSRawFitterv0():
TObject(),
- fSignal(0),
fModule(0),
fCellX(0),
fCellZ(0),
fCaloFlag(0),
- fStart(0),
- fLength(0),
fNBunches(0),
fPedSubtract(kFALSE),
fEnergy(-111),
AliPHOSRawFitterv0::~AliPHOSRawFitterv0()
{
//Destructor
- delete [] fSignal;
}
//-----------------------------------------------------------------------------
AliPHOSRawFitterv0::AliPHOSRawFitterv0(const AliPHOSRawFitterv0 &phosFitter ):
TObject(),
- fSignal (phosFitter.fSignal),
fModule (phosFitter.fModule),
fCellX (phosFitter.fCellX),
fCellZ (phosFitter.fCellZ),
fCaloFlag (phosFitter.fCaloFlag),
- fStart (phosFitter.fStart),
- fLength (phosFitter.fLength),
fNBunches (phosFitter.fNBunches),
fPedSubtract (phosFitter.fPedSubtract),
fEnergy (phosFitter.fEnergy),
//Assignment operator.
if(this != &phosFitter) {
- fSignal = phosFitter.fSignal;
fModule = phosFitter.fModule;
fCellX = phosFitter.fCellX;
fCellZ = phosFitter.fCellZ;
fCaloFlag = phosFitter.fCaloFlag;
- fStart = phosFitter.fStart;
- fLength = phosFitter.fLength;
fNBunches = phosFitter.fNBunches;
fPedSubtract = phosFitter.fPedSubtract;
fEnergy = phosFitter.fEnergy;
//-----------------------------------------------------------------------------
-void AliPHOSRawFitterv0::SetSamples(const UShort_t *sig, Int_t sigStart, Int_t sigLength)
-{
- // Set the sample array, its start and length in time bin units
-
- fStart = sigStart;
- fLength = sigLength;
- fSignal = new UShort_t[fLength];
- for (Int_t i=0; i<fLength; i++) {
- fSignal[i] = sig[i];
- }
-}
-//-----------------------------------------------------------------------------
-
void AliPHOSRawFitterv0::SetChannelGeo(const Int_t module, const Int_t cellX,
const Int_t cellZ, const Int_t caloFlag)
{
}
//-----------------------------------------------------------------------------
-Bool_t AliPHOSRawFitterv0::Eval()
+Bool_t AliPHOSRawFitterv0::Eval(const UShort_t *signal, Int_t sigStart, Int_t sigLength)
{
// Calculate signal parameters (energy, time, quality) from array of samples
// Energy is a maximum sample minus pedestal 9
// Time is the first time bin
// Signal overflows is there are at least 3 samples of the same amplitude above 900
+ fOverflow= kFALSE ;
fEnergy = 0;
if (fNBunches > 1) {
fQuality = 1000;
UShort_t maxSample = 0;
Int_t nMax = 0;
- for (Int_t i=0; i<fLength; i++) {
- if (i<kPreSamples) {
+ for (Int_t i=0; i<sigLength; i++) {
+ if (i>sigLength-kPreSamples) { //inverse signal time order
nPed++;
- pedMean += fSignal[i];
- pedRMS += fSignal[i]*fSignal[i] ;
+ pedMean += signal[i];
+ pedRMS += signal[i]*signal[i] ;
}
-
- if(fSignal[i] > maxSample) maxSample = fSignal[i];
- if(fSignal[i] == maxSample) nMax++;
+ if(signal[i] > maxSample){ maxSample = signal[i]; nMax=0;}
+ if(signal[i] == maxSample) nMax++;
- if(fPedSubtract) {
- if( (fSignal[i]-(Float_t)(pedMean/nPed)) >kBaseLine ) fTime = (Double_t)i;
- }
- else //ZS
- if( (fSignal[i]-(Float_t)fAmpOffset) >kBaseLine ) fTime = (Double_t)i;
}
+
fEnergy = (Double_t)maxSample;
if (maxSample > 900 && nMax > 2) fOverflow = kTRUE;
+ Double_t pedestal = 0 ;
if (fPedSubtract) {
if (nPed > 0) {
fPedestalRMS=(pedRMS - pedMean*pedMean/nPed)/nPed ;
if(fPedestalRMS > 0.)
fPedestalRMS = TMath::Sqrt(fPedestalRMS) ;
- fEnergy -= (Double_t)(pedMean/nPed); // pedestal subtraction
+ pedestal = (Double_t)(pedMean/nPed);
}
else
return kFALSE;
}
else {
//take pedestals from DB
- Double_t pedestal = (Double_t) fAmpOffset ;
+ pedestal = (Double_t) fAmpOffset ;
if (fCalibData) {
Float_t truePed = fCalibData->GetADCpedestalEmc(fModule, fCellZ, fCellX) ;
Int_t altroSettings = fCalibData->GetAltroOffsetEmc(fModule, fCellZ, fCellX) ;
pedestal += truePed - altroSettings ;
}
else{
- AliWarning(Form("Can not read data from OCDB")) ;
+ AliDebug(2,Form("Pedestal and offset are not read from OCDB. Use 0 for their values.")) ;
}
- fEnergy-=pedestal ;
}
+ fEnergy-=pedestal ;
if (fEnergy < kBaseLine) fEnergy = 0;
-
+
+ //Evaluate time
+ fTime = sigStart-sigLength-3;
+ const Int_t nLine= 6 ; //Parameters of fitting
+ const Float_t eMinTOF = 10. ; //Choosed from beam-test and cosmic analyis
+ const Float_t kAmp=0.35 ; //Result slightly depends on them, so no getters
+ // Avoid too low peak:
+ if(fEnergy < eMinTOF){
+ return kTRUE;
+ }
+
+ // Find index posK (kLevel is a level of "timestamp" point Tk):
+ Int_t posK =sigLength-1 ; //last point before crossing k-level
+ Double_t levelK = pedestal + kAmp*fEnergy;
+ while(signal[posK] <= levelK && posK>=0){
+ posK-- ;
+ }
+ posK++ ;
+
+ if(posK == 0 || posK==sigLength-1){
+ return kTRUE;
+ }
+
+ // Find crosing point by solving linear equation (least squares method)
+ Int_t np = 0;
+ Int_t iup=posK-1 ;
+ Int_t idn=posK ;
+ Double_t sx = 0., sy = 0., sxx = 0., sxy = 0.;
+ Double_t x,y ;
+
+ while(np<nLine){
+ //point above crossing point
+ if(iup>=0){
+ x = sigLength-iup-1;
+ y = signal[iup];
+ sx += x;
+ sy += y;
+ sxx += (x*x);
+ sxy += (x*y);
+ np++ ;
+ iup-- ;
+ }
+ //Point below crossing point
+ if(idn<sigLength){
+ if(signal[idn]<pedestal){
+ idn=sigLength-1 ; //do not scan further
+ idn++ ;
+ continue ;
+ }
+ x = sigLength-idn-1;
+ y = signal[idn];
+ sx += x;
+ sy += y;
+ sxx += (x*x);
+ sxy += (x*y);
+ np++;
+ idn++ ;
+ }
+ if(idn>=sigLength && iup<0){
+ break ; //can not fit futher
+ }
+ }
+
+ Double_t det = np*sxx - sx*sx;
+ if(det == 0){
+ return kTRUE;
+ }
+ if(np == 0){
+ return kFALSE;
+ }
+ Double_t c1 = (np*sxy - sx*sy)/det; //slope
+ Double_t c0 = (sy-c1*sx)/np; //offset
+ if(c1 == 0){
+ return kTRUE;
+ }
+
+ // Find where the line cross kLevel:
+ fTime += (levelK - c0)/c1-5. ; //5: mean offset between k-Level and start times
return kTRUE;
}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff