FitSpectra.h

Go to the documentation of this file.

#ifndef FITSPECTRA_H
#define FITSPECTRA_H

#include <vector>
#include <algorithm>

#include "TSpectrum.h"
#include "TString.h"
#include "TH1F.h"
#include "TCanvas.h"
#include "TStyle.h"
#include "TFile.h"
#include "TMath.h"
#include "TF1.h"
#include "TVirtualFitter.h"
#include "TROOT.h"
#include "TPad.h"
#include "TGraph.h"

class FitSpectra
{
public:
    FitSpectra(TString type = "Dino", TString source = "60Co");
    FitSpectra(TH1 *hist, TString type = "Dino", TString source = "60Co");
    ~FitSpectra();

    TVirtualPad *CurrentPad;

    double fWidth;
    void SetWidth(double w){fWidth = w;}
    double fThreshold;
    void SetThreshold(double t){fThreshold = t;}
    double fxmin;
    void SetXmin(double x){fxmin = x;}
    double fxmax;
    void SetXmax(double x){fxmax = x;}

    TH1 *fHistToFit;
    TH1 *GetHistogram(){return fHistToFit;}

    TString fHistName;
    TString GetHistName(){return fHistName;}

    TString fFitFctName;
    TString fBackGroundName;

    int fNPeaksToFit;
    std::vector< double > fEPeaks;
    std::vector < double > GetEPeaks(){return fEPeaks;}
    int fNFound;

    TH1F *BackGroundHist;
    TF1 *GausFct;
    TF1 *BGFct;

    TString fSourceName;
    TString fFitType;

    std::vector < float > fXPeaks;
    std::vector < float > fYPeaks;

    void InitSourceProperties();

    void SetHistogram(TH1 *h);
    void SetHistogram(TString HistName);
    void GetCurrentHistogram();

    void SetType(TString type){fFitType = type;}
    void SetSource(TString source);

    bool fDefinePersoBackground;
    int fNumberIterations;
    int fDirection;
    int fFilterOrder;
    bool fSmoothing;
    int fsmoothingWindow;
    bool fCompton;

    void SetPersoBackground(int NumberIterations = 100, int Direction = TSpectrum::kBackDecreasingWindow,int FilterOrder = TSpectrum::kBackOrder8,bool Smoothing = true,int smoothingWindow = TSpectrum::kBackSmoothing3,bool Compton = false);
    void UnsetPersoBackground(){fDefinePersoBackground = false;}

    TF1 *fFitFct;
    void Fit();

    int FitStatus;
    int GetFitStatus(){return FitStatus;}


    double BGPlusCBDS(double *xx, double *par);
    double BGPlusGaus(double *xx, double *par);
    double BGPlusDinoFct(double *xx, double *par);


    double MyGaus(double *xx, double *par);
    double fnc_dscb(double *xx, double *pp);
    double DinoFct(double*xx,double*pp);


    std::vector < double > fMean;
    std::vector < double > GetMean(){return fMean;}
    std::vector < double > fSigma;
    std::vector < double > GetSigma(){return fSigma;}
    std::vector < double > fFWHM;
    std::vector < double > GetFWHM(){return fFWHM;}
    std::vector < double > fFWHMCal;
    std::vector < double > GetFWHMCal(){return fFWHMCal;}
    std::vector < double > fFWHM_10;
    std::vector < double > GetFWHM_10(){return fFWHM_10;}
    std::vector < double > fFWHM_10Cal;
    std::vector < double > GetFWHM_10Cal(){return fFWHM_10Cal;}
    std::vector < double > fLeftTailProp;
    std::vector < double > GetLeftTailProp(){return fLeftTailProp;}
    std::vector < double > fRightTailProp;
    std::vector < double > GetRightTailProp(){return fRightTailProp;}
    std::vector < double > fArea;
    std::vector < double > GetArea(){return fArea;}
    std::vector < double > fEcal;
    std::vector < double > GetECal(){return fEcal;}

    double fCalibSlope;
    double GetCalibSlope(){return fCalibSlope;}
    double fCalibOffset;
    double GetCalibOffset(){return fCalibOffset;}


    void Calib();

    bool fUseOffset;
    void UseOffset(bool off){fUseOffset = off;}

    void ReduceListOfPeak();

//    ClassDef(FitSpectra,0)  // Watcher to work on Raw Apmlitude Signals spectra
};

#endif // FITSPECTRA_H