MdcCalibAlg/MdcCalibAlg-01-00-00/share/distcalib/src/include/fun.h

Go to the documentation of this file.

#ifndef FUN_H
#define FUN_H
 
#include <string>
#include <cstring>
#include <vector>
#include "TROOT.h"
#include "TObjArray.h"
 
#include "MdcCalibConst.h"
 
/* using namespace std; */
const int NLAYER = 43;
const int NWIRE = 6796;
const int NLR = 3; 
const int NXTBIN = 50;
const int NXTPAR = 8;
const int NENTRXT = 18;
const int NENTRSD = 6;
const int NSDBIN = 24;
 
const int gfixXtC0 = 1;
const double gTbinw = 10.0;
 
extern std::vector<double> XMEAS;
extern std::vector<double> TBINCEN;
extern std::vector<double> ERR;
extern double Tmax;
extern double Dmax;
extern std::vector<double> XMEASED;
extern std::vector<double> TBINCENED;
extern std::vector<double> ERRED;
 
extern int gNEntr[43];
 
/* from calib config file */
extern double gTimeShift;   /* if T<0 after subtracting Tes, use this */
extern double gTesMin;      /* minimun Tes for calibration */
extern double gTesMax;      /* maximun Tes for calibration */
extern int    gFgIniCalConst;   /* effective for IniMdcCalib */
extern bool   gPreT0SetTm;  /* flag for updating Tm in PreT0Calib */
extern double gInitT0;      /* initial value of T0 fit */
extern double gT0Shift;     /* t0 shift based on leading edge fitting */
extern double gTminFitChindf;   /* chisquare cut for Tmin fit */
extern double gTmaxFitChindf;   /* chisquare cut for Tmax fit */
extern int    gResiType;    /* 0: including measurement point; 1: excluding */
extern int    gCalSigma;    /* flag for update sigma */
extern int    gFixXtC0;     /* 1: fix c0 at 0 */
extern int    gFgCalib[NLAYER];
extern double gTminFitRange[NLAYER][2];
extern double gTmaxFitRange[NLAYER][2];
extern double gInitTm[NLAYER];      /* initial value of Tm fit */
extern double gQmin[NLAYER];        /* QT fit range */
extern double gQmax[NLAYER];        /* QT fit range */
 
double xtFun(double t, double xtpar[]);
void fcnXT(Int_t &npar, Double_t *gin, Double_t &f, Double_t *par, Int_t iflag);
void fcnXtEdge(Int_t &npar, Double_t *gin, Double_t &f, Double_t *par, Int_t iflag);
Double_t xtFitFun(Double_t *x, Double_t par[]);
Double_t xtFitEdge(Double_t *x, Double_t par[]);
void writeConst(MdcCalibConst* calconst, TObjArray* newXtList, TObjArray* r2tList);
std::vector<std::string> getHistList();
std::vector<std::string> getHistList(std::string path);
 
#endif