energyCalibration.h

Go to the documentation of this file.

#ifndef ENERGYCALIBRATION_H
#define ENERGYCALIBRATION_H

#include <TROOT.h>
#include <TF1.h>

using namespace std;

class TH1F;
class TGraphErrors;





const float conven=1000./3.6; 
const float el=1.67E-4; 
class energyCalibrationFunctions {

 public:
  
  energyCalibrationFunctions(int s=-1) {setScanSign(s);};
  
  int setScanSign(int s=0) {if (s==1 || s==-1) sign=s; return sign;};;
  

  Double_t gaussChargeSharing(Double_t *x, Double_t *par);

Double_t erfFunction(Double_t *x, Double_t *par) ;

Double_t erfFunctionChargeSharing(Double_t *x, Double_t *par);
  
Double_t erfFuncFluo(Double_t *x, Double_t *par);


  static float median(float *x, int n);
  static int quick_select(int arr[], int n);
  static int kth_smallest(int *a, int n, int k);

  
  Double_t spectrum(Double_t *x, Double_t *par);


  Double_t scurve(Double_t *x, Double_t *par);



  Double_t scurveFluo(Double_t *x, Double_t *par);


 private:
  int sign;
    

};

class energyCalibration  {


 public:
  energyCalibration();
    
  ~energyCalibration();
    
  int setPlotFlag(int p=-1) {if (p>=0) plot_flag=p; return plot_flag;};

  int setScanSign(int s=0) {return funcs->setScanSign(s);};
  
  int setChargeSharing(int p=-1);
  
  void setFitRange(Double_t mi, Double_t ma){fit_min=mi; fit_max=ma;};

  void getFitRange(Double_t &mi, Double_t &ma){mi=fit_min; ma=fit_max;};


  void setStartParameters(Double_t *par);
  
  void getStartParameters(Double_t *par);

  TF1 *fitSCurve(TH1 *h1, Double_t *mypar, Double_t *emypar);


  TF1 *fitSpectrum(TH1 *h1, Double_t *mypar, Double_t *emypar);


  TGraphErrors* linearCalibration(int nscan, Double_t *en, Double_t *een, Double_t *fl, Double_t *efl, Double_t &gain, Double_t &off, Double_t &egain, Double_t &eoff);  

  TGraphErrors* calibrateScurves(int nscan, Double_t *en, Double_t *een, TH1F **h1, Double_t &gain, Double_t &off, Double_t &egain, Double_t &eoff){return calibrate(nscan, en, een, h1, gain, off, egain, eoff, 1);};    

  TGraphErrors* calibrateSpectra(int nscan, Double_t *en, Double_t *een, TH1F **h1, Double_t &gain, Double_t &off, Double_t &egain, Double_t &eoff){return calibrate(nscan, en, een, h1, gain, off, egain, eoff, 0);};  



 private:
  TGraphErrors* calibrate(int nscan, Double_t *en, Double_t *een, TH1F **h1, Double_t &gain, Double_t &off, Double_t &egain, Double_t &eoff, int integral=1);  


  void initFitFunction(TF1 *fun, TH1 *h1);


  TF1 *fitFunction(TF1 *fun, TH1 *h1, Double_t *mypar, Double_t *emypar);  


  int plot_flag; 
  int cs_flag; 
  Double_t fit_min; 
  Double_t fit_max; 
  Double_t bg_offset; 
  Double_t bg_slope; 
  Double_t flex; 
  Double_t noise; 
  Double_t ampl; 
  Double_t cs_slope; 
  energyCalibrationFunctions *funcs;

  TF1 *fscurve; 
  TF1 *fspectrum; 
};

#endif

---