#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <sstream>
#include <vector>
#include "TFile.h"
#include "TTree.h"
#include "include/MdcCalibConst.h"
#include "include/fun.h"

Functions
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)
vector< string >	getHistList ()
vector< string >	getHistList (string path)

Variables
vector< double >	XMEAS
vector< double >	TBINCEN
vector< double >	ERR
double	Tmax
double	Dmax
vector< double >	XMEASED
vector< double >	TBINCENED
vector< double >	ERRED
int	gNEntr [43]
double	gTimeShift = 0.0
double	gTesMin = 0.0
double	gTesMax = 9999.0
int	gFgIniCalConst = 2
bool	gPreT0SetTm = true
double	gInitT0 = 50.0
double	gT0Shift = 0.0
double	gTminFitChindf = 20.0
double	gTmaxFitChindf = 20.0
int	gResiType = 0
int	gCalSigma = 1
int	gFixXtC0 = 0
int	gFgCalib [NLAYER]
double	gTminFitRange [NLAYER][2]
double	gTmaxFitRange [NLAYER][2]
double	gInitTm [NLAYER]
double	gQmin [NLAYER]
double	gQmax [NLAYER]

Function Documentation

◆ fcnXT()

void fcnXT	(	Int_t &	npar,
		Double_t *	gin,
		Double_t &	f,
		Double_t *	par,
		Int_t	iflag )

Definition at line 65 of file MdcCalibAlg/share/distcalib/src/fun.cpp.

                                                                                  {
  unsigned int i;
  int          ord;
  Double_t     deta;
  Double_t     chisq = 0.0;
  Double_t     dfit;
 
  for ( i = 0; i < TBINCEN.size(); i++ )
  {
    dfit = 0;
    for ( ord = 0; ord <= 5; ord++ ) { dfit += par[ord] * pow( TBINCEN[i], ord ); }
    deta = ( dfit - XMEAS[i] ) / ERR[i];
    chisq += deta * deta;
  }
 
  f = chisq;
}

Referenced by XtCalib::calib().

◆ fcnXtEdge()

void fcnXtEdge	(	Int_t &	npar,
		Double_t *	gin,
		Double_t &	f,
		Double_t *	par,
		Int_t	iflag )

Definition at line 83 of file MdcCalibAlg/share/distcalib/src/fun.cpp.

                                                                                      {
  unsigned int i;
  Double_t     deta;
  Double_t     chisq = 0.0;
  Double_t     dfit;
 
  for ( i = 0; i < TBINCENED.size(); i++ )
  {
    dfit = par[0] * ( TBINCENED[i] - Tmax ) + Dmax;
    deta = ( dfit - XMEASED[i] ) / ERRED[i];
    chisq += deta * deta;
  }
 
  f = chisq;
}

Referenced by XtCalib::calib().

◆ getHistList() [1/2]

vector< string > getHistList ( )

Definition at line 187 of file MdcCalibAlg/share/distcalib/src/fun.cpp.

                             {
  vector<string> fnames;
 
  string command( "JobOutputDir=`/bin/ls -dt1 joboutput-* 2>/dev/null | head -1`\n"
                  "if [ -d \"${JobOutputDir}\" ]; then\n"
                  "  find ${JobOutputDir} -name hist.root\n"
                  "fi\n" );
 
  stringstream fnstream;
 
  char* fnbuf   = new char[1024];
  FILE* fstream = popen( command.c_str(), "r" );
 
  while ( fgets( fnbuf, 1024, fstream ) != NULL ) { fnstream << fnbuf; }
 
  string fname;
  while ( !( fnstream >> fname ).eof() ) { fnames.push_back( fname ); }
 
  pclose( fstream );
  delete[] fnbuf;
 
  if ( fnames.empty() )
  {
    cout << "WARNING: Failed to retrieve hist files in the current directory!" << endl;
    //       exit(1);
  }
  return fnames;
}

◆ getHistList() [2/2]

vector< string > getHistList ( string path )

Definition at line 216 of file MdcCalibAlg/share/distcalib/src/fun.cpp.

                                          {
  vector<string>    fnames;
  string            newpath = path;
  string::size_type strl    = newpath.length();
  if ( ( strl > 1 ) && ( '/' == newpath[strl - 1] ) ) newpath.erase( strl - 1 );
 
  char com1[500];
  sprintf( com1, "JobOutputDir=`/bin/ls -dt1 %s/joboutput-* 2>/dev/null | head -1`\n",
           newpath.c_str() );
  string       command1( com1 );
  string       command2( "if [ -d \"${JobOutputDir}\" ]; then\n"
                               "  find ${JobOutputDir} -name hist.root\n"
                               "fi\n" );
  string       command = command1 + command2;
  stringstream fnstream;
 
  char* fnbuf   = new char[1024];
  FILE* fstream = popen( command.c_str(), "r" );
 
  while ( fgets( fnbuf, 1024, fstream ) != NULL ) { fnstream << fnbuf; }
 
  string fname;
  while ( !( fnstream >> fname ).eof() ) { fnames.push_back( fname ); }
 
  pclose( fstream );
  delete[] fnbuf;
 
  if ( fnames.empty() )
  {
    cout << "ERROR: Failed to retrieve hist files!" << endl;
    exit( 1 );
  }
  return fnames;
}

◆ writeConst()

void writeConst	(	MdcCalibConst *	calconst,
		TObjArray *	newXtList,
		TObjArray *	r2tList )

Definition at line 110 of file MdcCalibAlg/share/distcalib/src/fun.cpp.

                                                                                     {
  TFile fout( "MdcCalibConst_new.root", "recreate" );
 
  int    key;
  double xtpar;
  TTree* xttree = new TTree( "XtTree", "XtTree" );
  xttree->Branch( "xtkey", &key, "key/I" );
  xttree->Branch( "xtpar", &xtpar, "xtpar/D" );
  for ( int lay = 0; lay < 43; lay++ )
  {
    for ( int entr = 0; entr < 18; entr++ )
    {
      for ( int lr = 0; lr < 3; lr++ )
      {
        for ( int ord = 0; ord < 8; ord++ )
        {
          key   = calconst->getXtKey( lay, entr, lr, ord );
          xtpar = calconst->getXtpar( lay, entr, lr, ord );
          xttree->Fill();
        }
      }
    }
  }
 
  double t0;
  double delt0;
  TTree* t0tree = new TTree( "T0Tree", "T0Tree" );
  t0tree->Branch( "t0", &t0, "t0/D" );
  t0tree->Branch( "delt0", &delt0, "delt0/D" );
  for ( int wid = 0; wid < 6796; wid++ )
  {
    t0    = calconst->getT0( wid );
    delt0 = calconst->getDelT0( wid );
    t0tree->Fill();
  }
 
  double qtval[2];
  TTree* qttree = new TTree( "QtTree", "QtTree" );
  qttree->Branch( "qtpar0", &( qtval[0] ), "qtpar0/D" );
  qttree->Branch( "qtpar1", &( qtval[1] ), "qtpar1/D" );
  for ( int lay = 0; lay < 43; lay++ )
  {
    qtval[0] = calconst->getQtpar0( lay );
    qtval[1] = calconst->getQtpar1( lay );
    qttree->Fill();
  }
 
  double sdpar;
  TTree* sdtree = new TTree( "SdTree", "SdTree" );
  sdtree->Branch( "sdkey", &key, "key/I" );
  sdtree->Branch( "sdpar", &sdpar, "sdpar/D" );
  for ( int lay = 0; lay < 43; lay++ )
  {
    for ( int entr = 0; entr < 6; entr++ )
    {
      for ( int lr = 0; lr < 2; lr++ )
      {
        for ( int bin = 0; bin < 24; bin++ )
        {
          key   = calconst->getSdKey( lay, entr, lr, bin );
          sdpar = calconst->getSdpar( lay, entr, lr, bin );
          sdtree->Fill();
        }
      }
    }
  }
 
  fout.cd();
  xttree->Write();
  t0tree->Write();
  qttree->Write();
  sdtree->Write();
  if ( ( newXtList->GetEntries() ) > 0 ) newXtList->Write();
  if ( ( r2tList->GetEntries() ) > 0 ) r2tList->Write();
  fout.Close();
}

Referenced by main().

◆ xtFitEdge()

Double_t xtFitEdge	(	Double_t *	x,
		Double_t	par[] )

Definition at line 105 of file MdcCalibAlg/share/distcalib/src/fun.cpp.

                                                  {
  double val = Dmax + ( x[0] - Tmax ) * par[0];
  return val;
}

Referenced by GrXtCalib::calib().

◆ xtFitFun()

Double_t xtFitFun	(	Double_t *	x,
		Double_t	par[] )

Definition at line 99 of file MdcCalibAlg/share/distcalib/src/fun.cpp.

                                                 {
  Double_t val = 0.0;
  for ( Int_t ord = 0; ord < 6; ord++ ) { val += par[ord] * pow( x[0], ord ); }
  return val;
}

Referenced by GrXtCalib::calib().

◆ xtFun()

double xtFun	(	double	t,
		double	xtpar[] )

Definition at line 47 of file MdcCalibAlg/share/distcalib/src/fun.cpp.

                                         {
  int    ord;
  double dist = 0.0;
  double tm   = xtpar[6];
 
  if ( t < tm )
  {
    for ( ord = 0; ord <= 5; ord++ ) { dist += xtpar[ord] * pow( t, ord ); }
  }
  else
  {
    for ( ord = 0; ord <= 5; ord++ ) { dist += xtpar[ord] * pow( tm, ord ); }
    dist += xtpar[7] * ( t - tm );
  }
 
  return dist;
}