MdcCalibAlg/share/distcalib/src/include/fun.h

Go to the documentation of this file.

#ifndef FUN_H
#define FUN_H
 
#include "TObjArray.h"
#include "TROOT.h"
#include <cstring>
#include <string>
#include <vector>
 
#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