TGo4FitPeakFinder.h

Go to the documentation of this file.

// $Id: TGo4FitPeakFinder.h 478 2009-10-29 12:26:09Z linev $
//-----------------------------------------------------------------------
//       The GSI Online Offline Object Oriented (Go4) Project
//         Experiment Data Processing at EE department, GSI
//-----------------------------------------------------------------------
// Copyright (C) 2000- GSI Helmholtzzentrum für Schwerionenforschung GmbH
//                     Planckstr. 1, 64291 Darmstadt, Germany
// Contact:            http://go4.gsi.de
//-----------------------------------------------------------------------
// This software can be used under the license agreements as stated
// in Go4License.txt file which is part of the distribution.
//-----------------------------------------------------------------------

#ifndef TGO4FITPEAKFINDER_H
#define TGO4FITPEAKFINDER_H

#include "TGo4FitterAction.h"

#include "TArrayD.h"

class TGo4Fitter;
class TGo4FitData;

class TGo4FitPeakFinder : public TGo4FitterAction {
   public:

      TGo4FitPeakFinder();

      TGo4FitPeakFinder(const char* Name, const char* DataName = 0, Bool_t ClearModels = kFALSE, Int_t PolOrder = -1);

      virtual ~TGo4FitPeakFinder();

      void SetDataName(const char* name) { fxDataName = name; }
      const char* GetDataName() { return fxDataName.Data(); }
      void SetClearModels(Bool_t clear) { fbClearModels = clear; }
      Bool_t GetClearModels() { return fbClearModels; }


      void SetupPolynomialBackground(Int_t PolynomOrder);

      void SetupForFirst(Double_t MaxAmplFactor, Double_t MinWidth, Double_t MaxWidth);

      void SetupForSecond(Double_t LineWidth);

      void SetupForThird(Double_t NoiseFactor, Double_t NoiseMinimum, Int_t ChannelSum);

      virtual void DoAction(TGo4FitterAbstract* Fitter);

      virtual Bool_t CanChangeFitter() const { return kTRUE; }

      // Following functions directly changes single parameters of peakfinders.
      // Normally you should not use them, because some or all of them can be changed later on
      // You should use Setup... routines to configure specific peakfinder

      void SetUsePolynom(Bool_t use) { fbUsePolynom = use; }
      Bool_t GetUsePolynom() { return fbUsePolynom; }
      void SetPolynomOrder(Int_t order) { fiPolynomOrder = order; }
      Int_t GetPolynomOrder() { return fiPolynomOrder; }

      void SetPeakFinderType(Int_t typ) { fiPeakFinderType = typ; }
      Int_t GetPeakFinderType() { return fiPeakFinderType; }

      void Set0MaxAmplFactor(Double_t factor) { fd0MaxAmplFactor = factor; }
      Double_t Get0MaxAmplFactor() { return fd0MaxAmplFactor; }
      void Set0MinWidth(Double_t min) { fd0MinWidth = min; }
      Double_t Get0MinWidth() { return fd0MinWidth; }
      void Set0MaxWidth(Double_t max) { fd0MaxWidth = max; }
      Double_t Get0MaxWidth() { return fd0MaxWidth; }

      void Set1LineWidth(Double_t width) { fd1LineWidth = width; }
      Double_t Get1LineWidth() { return fd1LineWidth; }

      void Set2NoiseFactor(Double_t factor) { fd2NoiseFactor = factor; }
      Double_t Get2NoiseFactor() { return fd2NoiseFactor; }
      void Set2NoiseMinimum(Double_t min) { fd2NoiseMinimum = min; }
      Double_t Get2NoiseMinimum() { return fd2NoiseMinimum; }
      void Set2ChannelSum(Int_t sum) { fi2ChannelSum = sum; }
      Int_t Get2ChannelSum() { return fi2ChannelSum; }



      virtual void Print(Option_t* option) const;

    private:

      void ROOTPeakFinder(TGo4Fitter* fitter,
                          TGo4FitData* data,
                          Int_t PolynomOrder,
                          Double_t Sigma);


      void SergeyLinevPeakFinder(TGo4Fitter* fitter,
                                 TGo4FitData* data,
                                 Int_t PolOrder,
                                 Double_t AmplThreshold,
                                 Double_t MinWidth,
                                 Double_t MaxWidth);

      void HansEsselPeakFinder(TGo4Fitter* fitter,
                               TGo4FitData* data,
                               Int_t MaxNumPeaks = 50,
                               Int_t ChannelSum = 1,
                               Double_t NoiseFactor = 2.,
                               Double_t NoiseMinimum = 10.,
                               Int_t MinimasOrder = -1);

      static void DefinePolynom(Int_t size,              // number of experim points
                                Double_t* bins,          // arrray of bins
                                Double_t* scales,        // array of scale values
                                TArrayD& Coef,           // array of polynomial coefficients
                                Double_t* weight = 0,    // statistical weight of each bin
                                Double_t* backgr = 0,    // background values
                                Char_t* use = 0);        // usage of each points

      static void DefinePolynomEx(Int_t size, Double_t* bins, Double_t* scales, Double_t* weight, Double_t* backgr,
                                  Int_t lbound, Int_t rbound, TArrayD& Coef);


      static Double_t CalcPolynom(const TArrayD& Coef, Double_t x);


      Int_t fiPeakFinderType;   // 0 - by S.Linev,  1 - ROOT TSpectrum , 2 - by H.Essel

      TString fxDataName;

      Bool_t fbClearModels;

      Bool_t fbUsePolynom;
      Int_t  fiPolynomOrder;

      Double_t fd0MinWidth;
      Double_t fd0MaxWidth;
      Double_t fd0MaxAmplFactor;

      Double_t fd1LineWidth;

      Double_t fd2NoiseFactor;
      Double_t fd2NoiseMinimum;
      Int_t  fi2ChannelSum;

   ClassDef(TGo4FitPeakFinder,1)
};

#endif // TGO4FITPEAKFINDER_H