LaserCalibratorFit Class Reference

A class do laser calibration fit provide different fitting method (under development) More...

#include <LaserCalibratorFit.h>

Public Member Functions
	LaserCalibratorFit (unsigned moduleID)
	Constructor.
	~LaserCalibratorFit ()
	Destructor.
void	setHist (const std::vector< TH1F * > &hist)
	set time hist of all channels in one moduleID
void	setFitMethod (const std::string &method)
	set time fit function
void	setFitRange (double xmin=-200, double xmax=200)
	set x range in the fit
int	fitChannel (unsigned channel)
	fit for a specific channel
void	writeFile (const std::string &outfile)
	write fit result to a root file
int	getModuleID () const
	get the module ID in the fit
double	getFitChisq (unsigned channel)
	get chi^2 in the fit
std::string	getFitMethod () const
	get the fit method used in recent module
std::vector< double >	getMaxPos ()
	get the center positon of hist max bin
double	getFitT ()
	get mean positon after fit
double	getFitTErr ()
	returns the error mean positon after fit

Private Member Functions
TF1 *	makeGFit (unsigned channel)
	Fit process using single gaussian function.
TF1 *	makeCBFit (unsigned channel)
	Fit process using single Crystal Ball fuction.
TF1 *	makeCB2Fit (unsigned channel, bool minOut)
	Fit process using double Crystal Ball fuction.

Private Attributes
unsigned	m_moduleID = 0
	one moduleID/slot
std::vector< double >	m_maxpos
	center positon of hist max bin
std::vector< double >	m_maxpos_error
	error on the center positon of hist max bin
double	m_xmin = 0
	fitting low-edge
double	m_xmax = 0
	fitting upper-edge
double	m_fitT = 0
	mean position after fit
double	m_fitTErr = 0.
	error on the mean position estimated by the fit
std::vector< TH1F * >	m_hist
	time hist of 512 channels
std::vector< TF1 * >	m_func
	fitting function
std::string	m_fitMethod
	fitting method

Detailed Description

A class do laser calibration fit provide different fitting method (under development)

Definition at line 44 of file LaserCalibratorFit.h.

Constructor & Destructor Documentation

LaserCalibratorFit()

LaserCalibratorFit ( unsigned  moduleID )

explicit

Constructor.

Parameters

moduleID

module ID

Definition at line 54 of file LaserCalibratorFit.cc.

                                                           :
      m_moduleID(moduleID)
    {
    }

~LaserCalibratorFit()

~LaserCalibratorFit

(

)

Destructor.

Definition at line 59 of file LaserCalibratorFit.cc.

    {
      for (auto f : m_func) {
        delete f;
      }
      m_func.clear();
 
      for (auto h : m_hist) {
        delete h;
      }
      m_hist.clear();
    }

Member Function Documentation

fitChannel()

int fitChannel

(

unsigned

channel

)

fit for a specific channel

Parameters

channel

a specific channel (0 - 511)

Definition at line 92 of file LaserCalibratorFit.cc.

    {
      if (channel > m_hist.size()) {
        B2WARNING("Wrong channel!");
        return 0;
      }
      if (!m_hist[channel]) {
        B2WARNING("No Evts. in channel " << channel);
        return 0;
      }
      //require at least 10 evts for each fit
      if (m_hist[channel]->GetEntries() < 10) {
        B2WARNING("Too little Evts. for fitting (Evts. < 10)!");
        return 0;
      }
      m_hist[channel]->SetAxisRange(m_xmin, m_xmax);
      m_maxpos[channel] = m_hist[channel]->GetXaxis()->GetBinCenter(m_hist[channel]->GetMaximumBin());
      m_hist[channel]->SetAxisRange(m_maxpos[channel] - 1.5, m_maxpos[channel] + 2);
      m_maxpos_error[channel] = m_hist[channel]->GetXaxis()->GetBinWidth(m_hist[channel]->GetMaximumBin()); // Error = bin width ?
      if (m_fitMethod == "gauss" || m_fitMethod == "gaus") {
        m_func[channel] = makeGFit(channel);
      } else if (m_fitMethod == "cb") {
        m_func[channel] = makeCBFit(channel);
      } else if (m_fitMethod == "cb2") {
        m_func[channel] = makeCB2Fit(channel, 0);
      } else {
        B2WARNING("No matched fit method!");
        return 0;
      }
 
      return 1;
    }

getFitChisq()

double getFitChisq

(

unsigned

channel

)

get chi^2 in the fit

Parameters

channel

chi^2 for a channel

Definition at line 83 of file LaserCalibratorFit.cc.

    {
      if (m_func[channel]) {
        return m_func[channel]->GetChisquare();
      } else {
        return -9;
      };
    }

getFitMethod()

std::string getFitMethod ( ) const

inline

get the fit method used in recent module

Definition at line 102 of file LaserCalibratorFit.h.

{return m_fitMethod;}

getFitT()

double getFitT ( )

inline

get mean positon after fit

Definition at line 112 of file LaserCalibratorFit.h.

{return m_fitT;}

getFitTErr()

double getFitTErr ( )

inline

returns the error mean positon after fit

Definition at line 117 of file LaserCalibratorFit.h.

{return m_fitTErr;}

getMaxPos()

std::vector< double > getMaxPos ( )

inline

get the center positon of hist max bin

Definition at line 107 of file LaserCalibratorFit.h.

{return m_maxpos;}

getModuleID()

int getModuleID ( ) const

inline

get the module ID in the fit

Definition at line 91 of file LaserCalibratorFit.h.

{return m_moduleID;}

makeCB2Fit()

TF1 * makeCB2Fit ( unsigned  channel, bool  minOut  )

private

Fit process using double Crystal Ball fuction.

Parameters

channel	refers to a typical time separation of two main peaks
minOut	minimum printing need more studies

Definition at line 291 of file LaserCalibratorFit.cc.

    {
      TH1F* h = m_hist[channel];
      auto func = new TF1("fcnCB2", fcnCB2, m_xmin, m_xmax, 8);
 
      const double vdt[8] = {0.272, 0.242, 0.208, 0.178, 0.113, 0.082, 0.0485, 0.017}; //input para. from MC study, need further studies
      double parms[8];
      //the input para. need further studies
      parms[0] = 1;
      parms[1] = 0.6;
      parms[2] = 0.1;
      parms[3] = -0.5;
      parms[4] = 4;
      parms[5] = 1;
      parms[6] = vdt[(channel) / 64]; //refers to a typical time separation of two main peaks
      //parms[6] = 0.3;
      parms[7] = 0.1;
 
      func->SetParameters(parms);
      func->SetNpx(1000);
 
      func->SetParName(0, "norm1");
      func->SetParName(1, "time1");
      func->SetParName(2, "sigma1");
      func->SetParName(3, "alpha1_CB");
      func->SetParName(4, "n1_CB");
      func->SetParName(5, "norm2");
      func->SetParName(6, "dt12");
      func->SetParName(7, "sigma2");
      func->SetParName(8, "alpha2_CB");
      func->SetParName(9, "n2_CB");
 
      func->SetParLimits(1, parms[1] - 0.25, parms[1] + 0.15);
      //func->SetParLimits(2,par[2]-0.02,par[2]+0.08);
      func->SetParLimits(3, -5, 0);
      func->SetParLimits(4, 0, 10);
      func->SetParLimits(6, parms[6] * 0.7, parms[6] * 1.5);
      func->SetParLimits(8, -5, 0);
      func->SetParLimits(9, 0, 10);
 
      ROOT::Math::MinimizerOptions::SetDefaultMaxFunctionCalls(50000); //set to a large number for a test
      if (minOut) {
        h->Fit(func, "Q", "");
      } else {
        h->Fit(func);
      }
      m_fitT = parms[1];
      m_fitTErr = func->GetParError(1);
      return func;
    }

makeCBFit()

TF1 * makeCBFit ( unsigned  channel )

private

Fit process using single Crystal Ball fuction.

Parameters

channel

channel number

Definition at line 253 of file LaserCalibratorFit.cc.

    {
      TH1F* h = m_hist[channel];
      auto func = new TF1("fcnCB", fcnCB, m_xmin, m_xmax, 5);
 
      double parms[5];
      parms[0] = 1;
      parms[1] = h->GetMean();
      parms[2] = h->GetRMS();
      parms[3] = -0.5;
      parms[4] = 4;
 
      func->SetParameters(parms);
      func->SetNpx(1000);
 
      func->SetParName(0, "norm");
      func->SetParName(1, "time");
      func->SetParName(2, "sigma");
      func->SetParName(3, "alpha_CB");
      func->SetParName(4, "n_CB");
 
      //func->SetParLimits(1,par[1]-0.2,par[1]+0.2);
      //func->SetParLimits(2,par[2]-0.02,par[2]+0.08);
      func->SetParLimits(3, -5, 0);
      func->SetParLimits(4, 0, 10);
 
      h->Fit(func, "Q", "");
      if (fabs(m_maxpos[channel] - parms[1]) < parms[2]) {
        m_fitT = parms[1];
        m_fitTErr = func->GetParError(1);
      } else {
        m_fitT = m_maxpos[channel];
        m_fitTErr = m_maxpos_error[channel];
      }
      return func;
    }

makeGFit()

TF1 * makeGFit ( unsigned  channel )

private

Fit process using single gaussian function.

Parameters

channel

channel number

Definition at line 235 of file LaserCalibratorFit.cc.

    {
      TH1F* h = m_hist[channel];
      double m = h->GetMean();
      double w = h->GetRMS();
      h->Fit("gaus", "Q", "", m - 2.*w, m + 4.*w);
      double parms[3];
      TF1* func = h->GetFunction("gaus");
      func->GetParameters(parms);
      func->SetNpx(1000);
      h->Fit("gaus", "Q", "", m - w, m + w);
      func = h->GetFunction("gaus");
      m_fitT = parms[1];
      m_fitTErr = func->GetParError(1);
      return func;
    }

setFitMethod()

void setFitMethod ( const std::string &  method )

inline

set time fit function

Parameters

method

gauss: single gaussian; cb: single Crystal Ball; cb2: double Crystal Ball

Definition at line 67 of file LaserCalibratorFit.h.

{m_fitMethod = method;}

setFitRange()

void setFitRange ( double  xmin = -200, double  xmax = 200  )

inline

set x range in the fit

Parameters

xmin	lower x
xmax	upper x

Definition at line 74 of file LaserCalibratorFit.h.

{m_xmin = xmin; m_xmax = xmax;}

setHist()

void setHist

(

const std::vector< TH1F * > &

hist

)

set time hist of all channels in one moduleID

Parameters

hist	time hist vector

Definition at line 73 of file LaserCalibratorFit.cc.

    {
      m_hist = hist;
      for (unsigned i = 0; i < hist.size(); i++) {
        m_func.push_back(0);
        m_maxpos.push_back(0);
        m_maxpos_error.push_back(0);
      }
    }

writeFile()

void writeFile

(

const std::string &

outfile

)

write fit result to a root file

Parameters

outfile

output file name

Definition at line 125 of file LaserCalibratorFit.cc.

    {
      auto file = new TFile(outfile.c_str(), "RECREATE");
      auto otree = new TTree("fits", "fitted times");
 
      unsigned channel = 0;
      double maxpos = 0;
 
      otree->Branch("maxpos", &maxpos, "maxpos/D");
 
      if (m_fitMethod == "gauss") {
        double parms[3];
        double parmErrs[3];
        otree->Branch("channel", &channel, "channel/I");
        otree->Branch("norm", &(parms[0]), "norm/D");
        otree->Branch("time", &(parms[1]), "time/D");
        otree->Branch("reso", &(parms[2]), "reso/D");
        otree->Branch("normErr", &(parmErrs[0]), "normErr/D");
        otree->Branch("timeErr", &(parmErrs[1]), "timeErr/D");
        otree->Branch("resoErr", &(parmErrs[2]), "resoErr/D");
 
        for (auto& f : m_func) {
          maxpos = m_maxpos[channel];
          if (f) {
            f->GetParameters(parms);
            for (int iPar = 0; iPar < 3; iPar++)
              parmErrs[iPar] = f->GetParError(iPar); // cppcheck-suppress unreadVariable
            otree->Fill();
            channel++;
          }
        }
      } else if (m_fitMethod == "cb2") {
        double parms[8];
        double parmErrs[8];
 
        double time2 = 0;
        otree->Branch("channel", &channel, "channel/I");
        otree->Branch("norm1", &(parms[0]), "norm1/D");
        otree->Branch("time1", &(parms[1]), "time1/D");
        otree->Branch("reso1", &(parms[2]), "reso1/D");
        otree->Branch("alpha_CB", &(parms[3]), "alpha_CB/D");
        otree->Branch("n_CB", &(parms[4]), "n_CB/D");
        otree->Branch("norm2", &(parms[5]), "norm2/D");
        otree->Branch("dt12", &(parms[6]), "dt12/D");
        otree->Branch("reso2", &(parms[7]), "reso2/D");
        otree->Branch("time2", &time2, "time2/D");
 
        otree->Branch("norm1Err", &(parmErrs[0]), "norm1Err/D");
        otree->Branch("time1Err", &(parmErrs[1]), "time1Err/D");
        otree->Branch("reso1Err", &(parmErrs[2]), "reso1Err/D");
        otree->Branch("alpha_CBErr", &(parmErrs[3]), "alpha_CBErr/D");
        otree->Branch("n_CBErr", &(parmErrs[4]), "n_CBErr/D");
        otree->Branch("norm2Err", &(parmErrs[5]), "norm2Err/D");
        otree->Branch("dt12Err", &(parmErrs[6]), "dt12Err/D");
        otree->Branch("reso2Err", &(parmErrs[7]), "reso2Err/D");
 
 
        for (auto& f : m_func) {
          maxpos = m_maxpos[channel];
          if (f) {
            f->GetParameters(parms);
            for (int iPar = 0; iPar < 8; iPar++)
              parmErrs[iPar] = f->GetParError(iPar); // cppcheck-suppress unreadVariable
            otree->Fill();
            channel++;
          }
        }
      } else if (m_fitMethod == "cb") {
        double parms[5];
        double parmErrs[5];
 
        otree->Branch("channel", &channel, "channel/I");
        otree->Branch("norm", &(parms[0]), "norm/D");
        otree->Branch("time", &(parms[1]), "time/D");
        otree->Branch("reso", &(parms[2]), "reso/D");
        otree->Branch("alpha_CB", &(parms[3]), "alpha_CB/D");
        otree->Branch("n_CB", &(parms[4]), "n_CB/D");
 
        otree->Branch("normErr", &(parmErrs[0]), "normErr/D");
        otree->Branch("timeErr", &(parmErrs[1]), "timeErr/D");
        otree->Branch("resoErr", &(parmErrs[2]), "resoErr/D");
        otree->Branch("alpha_CBErr", &(parmErrs[3]), "alpha_CBErr/D");
        otree->Branch("n_CBErr", &(parmErrs[4]), "n_CBErr/D");
 
        for (auto& f : m_func) {
          maxpos = m_maxpos[channel];
          if (f) {
            f->GetParameters(parms);
            for (int iPar = 0; iPar < 5; iPar++)
              parmErrs[iPar] = f->GetParError(iPar); // cppcheck-suppress unreadVariable
            otree->Fill();
            channel++;
          }
        }
      } else {
        B2WARNING("Nothing be written to files!");
        return;
      }
      for (auto& h : m_hist) {
        if (h && h->GetEntries() >= 10) {
          h->Write();
        }
      }
      otree->Write();
      delete otree;
      file->Close();
      delete file;
    }

Member Data Documentation

m_fitMethod

std::string m_fitMethod

private

fitting method

Definition at line 151 of file LaserCalibratorFit.h.

m_fitT

double m_fitT = 0

private

mean position after fit

Definition at line 147 of file LaserCalibratorFit.h.

m_fitTErr

double m_fitTErr = 0.

private

error on the mean position estimated by the fit

Definition at line 148 of file LaserCalibratorFit.h.

m_func

std::vector<TF1*> m_func

private

fitting function

Definition at line 150 of file LaserCalibratorFit.h.

m_hist

std::vector<TH1F*> m_hist

private

time hist of 512 channels

Definition at line 149 of file LaserCalibratorFit.h.

m_maxpos

std::vector<double> m_maxpos

private

center positon of hist max bin

Definition at line 143 of file LaserCalibratorFit.h.

m_maxpos_error

std::vector<double> m_maxpos_error

private

error on the center positon of hist max bin

Definition at line 144 of file LaserCalibratorFit.h.

m_moduleID

unsigned m_moduleID = 0

private

one moduleID/slot

Definition at line 142 of file LaserCalibratorFit.h.

m_xmax

double m_xmax = 0

private

fitting upper-edge

Definition at line 146 of file LaserCalibratorFit.h.

m_xmin

double m_xmin = 0

private

fitting low-edge

Definition at line 145 of file LaserCalibratorFit.h.

---

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

• top/modules/TOPLaserCalibrator/include/LaserCalibratorFit.h
• top/modules/TOPLaserCalibrator/src/LaserCalibratorFit.cc