Class to perform fitting for each xt function. More...

#include <XT.h>

Public Member Functions
	XT (TProfile *h1)
	Initialized with TProfile histogram.

	XT (TProfile *h1, int mode)
	Initialized with TProfile histogram and mode.

	XT (TH1D *h1, int mode)
	Initialized with TH1D histogram and mode.

void	setP6 (double p6)
	Set Parameter 6 for polynomial fit.

void	setMode (int mode)
	Set XT mode.

void	BField (bool bfield)
	set to use BField

void	setXTParams (const double p[8])
	Set Parameters for fit.

void	setXTParams (double p0, double p1, double p2, double p3, double p4, double p5, double p6, double p7)
	Set Initial parameters for fitting.

void	setFitRange (double tmin, double tmax)
	Set Fit range.

void	setSmallestEntryRequired (int min)
	Set minimum number of entry required for fit.

void	setDebug (bool debug)
	Set Debug.

int	getFitStatus ()
	get fitted flag.

double	getProb ()
	Get the chi2 probability.

void	getFittedXTParams (double pa[8])
	get fit parameters.

TF1 *	getXTFunction (int mode)
	Get XT function.

TF1 *	getXTFunction ()
	Get XT function.

TProfile *	getFittedHisto ()
	Get histogram.

void	FitXT ()
	Do fitting.

void	FitXT (int mode)
	Do fitting.

void	FitPol5 ()
	Fit xt histogram incase 5th order polynomial is used.

void	FitChebyshev ()
	Fit xt histogram incase 5th order Chebeshev polynomial is used.

Private Attributes
TProfile *	m_h1
	Histogram of xt relation.

TF1 *	xtpol5 = new TF1("xtpol5", pol5pol1, 0.0, 700, 8)
	5th order polynomial function

TF1 *	xtCheb5 = new TF1("xtCheb5", Cheb5pol1, 0.0, 700, 8)
	5th order Cheb.

int	m_mode = 1
	XT mode, 0 is for 5th order polynomial, 1 is Chebshev polynomial.

bool	m_debug = true
	Print debug durring fitting or not.

bool	m_draw = false
	Draw and store png plot of each histo or not.

bool	m_BField = true
	With magnetic field or not.

int	m_minRequiredEntry = 800
	Minimum entry required for each histo.

double	m_XTParam [8] = {}
	Parameter fo xt.

double	m_FittedXTParams [8] = {}
	Fitted parameters.

int	m_fitflag = 0
	Fit Flag =-1: low statitic =1: good =0: Fit failure =2: Error Outer =3: Error Inner part;.

double	m_Prob = 0
	Chi2 prob of fitting.

double	m_tmin = 20
	lower boundary of fit range

double	m_tmax = m_XTParam[6] + 250
	upper boundary of fit range

Detailed Description

Class to perform fitting for each xt function.

Definition at line 60 of file XT.h.

Constructor & Destructor Documentation

◆ XT() [1/3]

XT ( TProfile * h1 )

inlineexplicit

Initialized with TProfile histogram.

Definition at line 65 of file XT.h.

  {
    m_h1 = (TProfile*)h1->Clone();
    m_h1->SetDirectory(0);
  }

◆ XT() [2/3]

XT	(	TProfile *	h1,
		int	mode
	)

inline

Initialized with TProfile histogram and mode.

Definition at line 73 of file XT.h.

  {
    m_h1 = (TProfile*)h1->Clone();
    m_h1->SetDirectory(0);
    m_mode = mode;
  }

◆ XT() [3/3]

XT	(	TH1D *	h1,
		int	mode
	)

inline

Initialized with TH1D histogram and mode.

Definition at line 83 of file XT.h.

  {
    m_h1 = (TProfile*)h1->Clone();
    m_h1->SetDirectory(0);
    m_mode = mode;
  }

Member Function Documentation

◆ BField()

void BField ( bool bfield )

inline

set to use BField

Definition at line 110 of file XT.h.

110{m_BField = bfield;}

XT::m_BField

bool m_BField

With magnetic field or not.

Definition: XT.h:231

◆ FitChebyshev()

void FitChebyshev ( )

Fit xt histogram incase 5th order Chebeshev polynomial is used.

Definition at line 332 of file XT.h.

{
  if (m_h1->GetEntries() < m_minRequiredEntry) {
    m_fitflag = -1;
    return;
  }
  //  m_tmax = m_XTParam[6] + 100;
  //xtCheb5->SetParameters(0.0, 0.005, 0., 0., 0., 0., m_XTParam[6], 0.001);
  double p[6];
  double par[8];
  xtCheb5->SetParLimits(7, 0., 0.001);
  m_h1->Fit("chebyshev5", "QME", "", m_tmin, m_XTParam[6]);
  m_h1->GetFunction("chebyshev5")->GetParameters(p);
  xtCheb5->SetParameters(p[0], p[1], p[2], p[3], p[4], p[5], m_XTParam[6], 0.000);
 
  double stat;
  for (int i = 0; i < 10; ++i) {
    stat = m_h1->Fit("xtCheb5", "MQ", "0", m_tmin, m_tmax);
    if (stat == 0) {
      xtCheb5->SetParameters(p[0], p[1], p[2], p[3], p[4], p[5], m_XTParam[6] - 20, 0.000);
      m_tmax -= 10;
      continue;
    }
    xtCheb5->GetParameters(par);
    /*Eval outer region,*/
    double fp6 = xtCheb5->Eval(par[6]);
    double fbehindp6 = xtCheb5->Eval(par[6] - 10) - 0.005;
    if (fp6 < fbehindp6 || fp6 > 1) { /*may be change to good value*/
      m_fitflag = 2;
      //      out += 1;
      xtCheb5->SetParameters(p[0], p[1], p[2], p[3], p[4], p[5], par[6] - 20, 0.000);
      xtCheb5->SetParLimits(6, par[6] - 50,  par[6] - 10);
      m_tmax -= 10;
      //      if (m_tmax < p6default + 30) {
      //  m_tmax = p6default + 30;
      // }
    } else {
      break;
    }
    //    m_tmax +=10;
    //if (stat != 0) break;
  }
  xtCheb5->GetParameters(m_FittedXTParams);
  m_Prob = xtCheb5->GetProb();
  if (stat == 0)
    m_fitflag = 0;
  else
    m_fitflag = 1;
  //  xtCheb5->DrawF1(0,400,"same");
  if (m_draw) {
    TString hname = m_h1->GetName();
    TString name = hname + ".pdf";
    TCanvas* c1 = new TCanvas("c1", "", 800, 600);
    m_h1->Draw();
    xtCheb5->SetLineColor(kBlack);
    xtCheb5->DrawF1(0, 400, "same");
    c1->SaveAs(name);
  }
}

◆ FitPol5()

void FitPol5 ( )

Fit xt histogram incase 5th order polynomial is used.

Definition at line 250 of file XT.h.

{
  double max_dif = 0.12;
  double max_dif2 = 0.05;
  double par[8];
  m_h1->Fit("pol1", "MQ", "", m_tmin, 50);
  TF1* f1 = (TF1*)m_h1->GetFunction("pol1");
  double p0 = f1->GetParameter(0);
  double p1 = f1->GetParameter(1);
  double f10 = f1->Eval(10);
  /****************************/
  int out = 0; /*how many time outer part change fit limit*/
  xtpol5->SetParameters(p0, p1, 0, 0, 0, 0, m_XTParam[6], 0);
  double p6default = m_XTParam[6];
  if (m_BField) {
    xtpol5->SetParLimits(7, 0.0, 0.001);
    xtpol5->SetParLimits(1, 0.0, 0.01);
  } else {
    xtpol5->SetParLimits(7, 0.0, 0.01);
    xtpol5->SetParLimits(1, 0.0, 0.01);
  }
  //  xtpol5->SetParLimits(6, p6default - 30,  p6default + 30);
  for (int i = 0; i < 10; ++i) {
    m_fitflag = 1;
    std::cout << "Fitting" << std::endl;
    double stat = m_h1->Fit(xtpol5, "M", "0", m_tmin, m_tmax);
    xtpol5->GetParameters(par);
 
    /*Eval outer region,*/
    double fp6 = xtpol5->Eval(par[6]);
    double fbehindp6 = xtpol5->Eval(par[6] - 12) - 0.01;
    if (fp6 < fbehindp6 || fp6 > 1) { /*may be change to good value*/
      m_fitflag = 2;
      out += 1;
      xtpol5->SetParameters(p0, p1, 0, 0, 0, 0, p6default, 0);
      xtpol5->SetParLimits(6, p6default - 10,  p6default + 10 - out / 2);
      m_tmax += 2;
      if (m_tmax < p6default + 30) {
        m_tmax = p6default + 30;
      }
    }
    /* EVal inner region*/
    /* compare p0 of linear fit vs p0 of xtfun fit and eval point t=10)*/
    if (fabs(par[0] - p0) > max_dif || fabs(f10 - xtpol5->Eval(10)) > max_dif2) {
      m_fitflag = 3;
      if (i == 9) std::cout << "ERROR XT FIT inner part" << std::endl;
      xtpol5->SetParameters(p0, p1, 0, 0, 0, 0, p6default, 0);
      xtpol5->SetParLimits(1, 0, 0.08);
      m_tmin -= 0.5;
 
      if (m_tmin < 14) {
        m_tmin = 14; std::cout << "ERROR: tmin so small, fit iter: " << std::endl;
      }
    }
 
    if (stat == 0) {
      m_fitflag = 0;
      if (m_debug) {std::cout << "FIT Failure; Layer: " << std::endl;}
    }
 
    if (m_debug) {printf("P6default= %3.2f, p6 fit = %3.2f", p6default, par[6]);}
 
    if (m_fitflag == 1) {
      if (m_debug) std::cout << "Fit success" << std::endl;
      xtpol5->GetParameters(m_FittedXTParams);
      m_Prob = xtpol5->GetProb();
      break;
    }
 
  } //end loop of fitting
  if (m_draw) {
    TString hname = m_h1->GetName();
    TString name = hname + ".pdf";
    TCanvas* c1 = new TCanvas("c1", "", 800, 600);
    m_h1->Draw();
    xtpol5->SetLineColor(kBlack);
    xtpol5->DrawF1(0, 400, "same");
    c1->SaveAs(name);
 
  }
}

◆ FitXT() [1/2]

void FitXT ( )

inline

Do fitting.

Definition at line 200 of file XT.h.

  {
    if (m_mode == 0) FitPol5();
    else FitChebyshev();
  }

◆ FitXT() [2/2]

void FitXT ( int mode )

inline

Do fitting.

Definition at line 209 of file XT.h.

  {
    if (mode == 0) FitPol5();
    else FitChebyshev();
  }

◆ getFitStatus()

int getFitStatus ( )

inline

get fitted flag.

Definition at line 157 of file XT.h.

  {
    return m_fitflag;
  }

◆ getFittedHisto()

TProfile * getFittedHisto ( )

inline

Get histogram.

Definition at line 195 of file XT.h.

195{return m_h1;}

◆ getFittedXTParams()

void getFittedXTParams ( double pa[8] )

inline

get fit parameters.

Definition at line 172 of file XT.h.

  {
    for (int i = 0; i < 8; ++i) {pa[i] = m_FittedXTParams[i];}
  }

◆ getProb()

double getProb ( )

inline

Get the chi2 probability.

Definition at line 165 of file XT.h.

  {
    return m_Prob;
  }

◆ getXTFunction() [1/2]

TF1 * getXTFunction ( )

inline

Get XT function.

Definition at line 187 of file XT.h.

  {
    if (m_mode == 0) return xtpol5;
    else return xtCheb5;
  };

◆ getXTFunction() [2/2]

TF1 * getXTFunction ( int mode )

inline

Get XT function.

Definition at line 179 of file XT.h.

  {
    if (mode == 0) return xtpol5;
    else return xtCheb5;
  }

◆ setDebug()

void setDebug ( bool debug )

inline

Set Debug.

Definition at line 150 of file XT.h.

  {
    m_debug = debug;
  }

◆ setFitRange()

void setFitRange	(	double	tmin,
		double	tmax
	)

inline

Set Fit range.

Definition at line 135 of file XT.h.

  {
    m_tmin = tmin;
    m_tmax = tmax;
  }

◆ setMode()

void setMode ( int mode )

inline

Set XT mode.

1 is 5th order Chebyshev polynomial. 0 is 5th order polynomial.

Definition at line 103 of file XT.h.

  {
    m_mode = mode;
  }

◆ setP6()

void setP6 ( double p6 )

inline

Set Parameter 6 for polynomial fit.

Definition at line 93 of file XT.h.

  {
    m_XTParam[6] = p6;
  }

◆ setSmallestEntryRequired()

void setSmallestEntryRequired ( int min )

inline

Set minimum number of entry required for fit.

Definition at line 143 of file XT.h.

  {
    m_minRequiredEntry = min;
  }

◆ setXTParams() [1/2]

void setXTParams ( const double p[8] )

inline

Set Parameters for fit.

Definition at line 115 of file XT.h.

  {
    for (int i = 0; i < 8; ++i) {m_XTParam[i] = p[i];}
    m_tmax = p[6] + 250;
  }

◆ setXTParams() [2/2]

void setXTParams	(	double	p0,
		double	p1,
		double	p2,
		double	p3,
		double	p4,
		double	p5,
		double	p6,
		double	p7
	)

inline

Set Initial parameters for fitting.

Definition at line 123 of file XT.h.

  {
    m_XTParam[0] = p0;     m_XTParam[1] = p1;     m_XTParam[2] = p2;
    m_XTParam[3] = p3;     m_XTParam[4] = p4;     m_XTParam[5] = p5;
    m_XTParam[6] = p6;     m_XTParam[7] = p7;
    m_tmax = p6 + 250;
  }

Member Data Documentation

◆ m_BField

bool m_BField = true

private

With magnetic field or not.

Definition at line 231 of file XT.h.

◆ m_debug

bool m_debug = true

private

Print debug durring fitting or not.

Definition at line 229 of file XT.h.

◆ m_draw

bool m_draw = false

private

Draw and store png plot of each histo or not.

Definition at line 230 of file XT.h.

◆ m_fitflag

int m_fitflag = 0

private

Fit Flag =-1: low statitic =1: good =0: Fit failure =2: Error Outer =3: Error Inner part;.

Definition at line 244 of file XT.h.

◆ m_FittedXTParams

double m_FittedXTParams[8] = {}

private

Fitted parameters.

Definition at line 234 of file XT.h.

◆ m_h1

TProfile* m_h1

private

Histogram of xt relation.

Definition at line 224 of file XT.h.

◆ m_minRequiredEntry

int m_minRequiredEntry = 800

private

Minimum entry required for each histo.

Definition at line 232 of file XT.h.

◆ m_mode

int m_mode = 1

private

XT mode, 0 is for 5th order polynomial, 1 is Chebshev polynomial.

Definition at line 228 of file XT.h.

◆ m_Prob

double m_Prob = 0

private

Chi2 prob of fitting.

Definition at line 245 of file XT.h.

◆ m_tmax

double m_tmax = m_XTParam[6] + 250

private

upper boundary of fit range

Definition at line 247 of file XT.h.

◆ m_tmin

double m_tmin = 20

private

lower boundary of fit range

Definition at line 246 of file XT.h.

◆ m_XTParam

double m_XTParam[8] = {}

private

Parameter fo xt.

Definition at line 233 of file XT.h.

◆ xtCheb5

TF1* xtCheb5 = new TF1("xtCheb5", Cheb5pol1, 0.0, 700, 8)

private

5th order Cheb.

polynomial function

Definition at line 226 of file XT.h.

◆ xtpol5

TF1* xtpol5 = new TF1("xtpol5", pol5pol1, 0.0, 700, 8)

private

5th order polynomial function

Definition at line 225 of file XT.h.

The documentation for this class was generated from the following file:

cdc/calibration/include/XT.h

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ XT() [1/3]

◆ XT() [2/3]

◆ XT() [3/3]

Member Function Documentation

◆ BField()

◆ FitChebyshev()

◆ FitPol5()

◆ FitXT() [1/2]

◆ FitXT() [2/2]

◆ getFitStatus()

◆ getFittedHisto()

◆ getFittedXTParams()

◆ getProb()

◆ getXTFunction() [1/2]

◆ getXTFunction() [2/2]

◆ setDebug()

◆ setFitRange()

◆ setMode()

◆ setP6()

◆ setSmallestEntryRequired()

◆ setXTParams() [1/2]

◆ setXTParams() [2/2]

Member Data Documentation

◆ m_BField

◆ m_debug

◆ m_draw

◆ m_fitflag

◆ m_FittedXTParams

◆ m_h1

◆ m_minRequiredEntry

◆ m_mode

◆ m_Prob

◆ m_tmax

◆ m_tmin

◆ m_XTParam

◆ xtCheb5

◆ xtpol5