KVAlphaCalibration.h

/*
  $Id: KVAlphaCalibration.h, v 1.0 2019/05/14 17:00:00 lemarie Exp $
  $Revision: 1.0 $
  $Date: 2019/05/14 17:00:00 $
  $Author: lemarie $
*/
 
#ifndef _KVCALIBRATION_H_
#define _KVCALIBRATION_H_
 
#include "TH1.h"
#include "TF1.h"
#include "TSpectrum.h"
#include "TGraph.h"
#include <string>
#include <TMath.h>
#include "TCanvas.h"
#include "TSystem.h"
#include "TFile.h"
#include "TStyle.h"
#include "TROOT.h"
#include <vector>
#include <cstdlib>
#include <iostream>
#include <limits>
 
class KVAlphaCalibration {
 
protected :
    double SigmaOfTSpectrum;
    double SigmaOfGaussian;
    double ThresholdOfTSpectrum;
   bool IsOriginAtZero;
    int NumberOfPeak;
    TSpectrum* spec;
    double InitializationFitResults[2];
    std::vector<double> GaussianFitResults;
    std::vector<double> GaussianFitResultsError;
   double FunctionToFit(double* x, double* par); //the model we use to fit the peak
   std::vector<double> InitializationPeak;
 
 
private :
    TH1* Histo;
    TF1* InitializationFit;
    TF1* GaussianFit;
    TGraph* factorGraph;
    std::vector<double> MeanOfPeak;
    std::vector<double> IntensityOfPeak;
   void HistoInit(TH1* h); //-
   TF1 GetFunction(void);
   void SetFunction(std::string FunctionName = "x*[0]");
 
public :
 
   KVAlphaCalibration(int NumberOfPeak_); //-
   KVAlphaCalibration(int NumberOfPeak_, TH1* h); //-
   ~KVAlphaCalibration();
 
   void Init(int);
   void SetHisto(TH1* h);
   void SetParameters(double SigmaOfTSpectrum_ = 1., double SigmaOfGaussian_ = 1., double ThresholdOfTSpectrum_ = 0.5, double IsOriginAtZero_ = false);
   void AddPeak(double Energy_, double Intensity_); //-
   void SetHistRange(double xmin, double xmax);
   double GetGaussianFitParameter(int); //-
   double GetGaussianFitParError(int);
   double GetInitializationFitParameter(int);
 
 
 
   TGraph* FitInit(bool debug_ = false);
   void FitSpectrum(bool debug_ = false);
   void FitAll(bool debug_ = false);
 
   void DrawResult(bool WhatToDraw = true);
   void PrintResult(void);
 
};
 
#endif