calib_barrel_sigma Class Reference

#include <calib_barrel_sigma.h>

Inheritance diagram for calib_barrel_sigma:

List of all members.

Public Member Functions
	calib_barrel_sigma (const unsigned int nzbin)
	~calib_barrel_sigma ()
void	calculate (RecordSet *&data, unsigned int icounter)
const string &	name () const
void	fillTxt (const char *file)
void	fillRoot (const char *file)
HepVector	tcorrelation ()
void	setTCorrelation (HepVector tc)
Protected Attributes
int	m_npar
unsigned int	nKind
unsigned int	nBinPerCounter
unsigned int	nHistPerCounter
unsigned int	nCanvasPerCounter
std::vector< unsigned int >	nGraphPerCanvasPerCounter
unsigned int	nHistogram
unsigned int	nCanvas
std::vector< unsigned int >	nGraphPerCanvas
string	m_name
HepVector	X
std::vector< TH1F * >	m_histograms
std::vector< TH1F * >	m_graphs
std::vector< HepVector >	m_result
std::vector< string >	CanvasPerCounterName
std::vector< string >	CanvasName
HepVector	m_tcorrelation
Private Member Functions
void	fillRecord (const Record *r, unsigned int icounter)
void	fitHistogram (unsigned int icounter)
void	fillGraph (unsigned int icounter)
void	fitGraph (unsigned int icounter)
void	updateData (Record *r, unsigned int icounter)
void	fillRecordT0 (const Record *r, unsigned int icounter)
void	fitHistogramT0 (unsigned int icounter)
void	fillGraphT0 (unsigned int icounter)
void	fitGraphT0 (unsigned int icounter)
Private Attributes
std::vector< double >	zpos
std::vector< double >	zposerr
double	zstep
std::vector< HepVector >	m_fitresult

---

Detailed Description

Definition at line 10 of file calib_barrel_sigma.h.

---

Constructor & Destructor Documentation

calib_barrel_sigma::calib_barrel_sigma

(

const unsigned int

nzbin

)

Definition at line 22 of file calib_barrel_sigma.cxx.

References TofCalibFit::CanvasPerCounterName, genRecEmupikp::i, iter(), ganga-rec::j, m_fitresult, TofCalibFit::m_histograms, TofCalibFit::m_name, TofCalibFit::m_result, NBarrel, nBarrelSigma, TofCalibFit::nBinPerCounter, TofCalibFit::nCanvas, TofCalibFit::nCanvasPerCounter, TofCalibFit::nGraphPerCanvasPerCounter, nGraphTotalSigma, TofCalibFit::nHistogram, TofCalibFit::nHistPerCounter, TofCalibFit::nKind, nParSigma, deljobs::string, zbegin, zend, zpos, zposerr, and zstep.

                                                                :TofCalibFit( true, nBarrelSigma ) {

  nKind          = 5;    // 0:tleft, 1:tright, 2:t0, 3:tplus, 4:tminus
  nBinPerCounter = nzbin;

  nHistPerCounter   = nKind*nBinPerCounter;
  nCanvasPerCounter = 4;
  CanvasPerCounterName.push_back( static_cast<string>("Barrel-offset") );  
  CanvasPerCounterName.push_back( static_cast<string>("Offset-TimeCorrelation") );  
  CanvasPerCounterName.push_back( static_cast<string>("Barrel-sigma") );  
  CanvasPerCounterName.push_back( static_cast<string>("Sigma-TimeCorrelation") );  
  nGraphPerCanvasPerCounter.push_back(3);
  nGraphPerCanvasPerCounter.push_back(2);
  nGraphPerCanvasPerCounter.push_back(3);
  nGraphPerCanvasPerCounter.push_back(3);

  nHistogram      = 0;
  nCanvas         = 0;

  int numGraphs = 0;
  std::vector<unsigned int>::iterator iter = nGraphPerCanvasPerCounter.begin();
  for( ; iter!=nGraphPerCanvasPerCounter.end(); iter++ ) {
    numGraphs = numGraphs + (*iter);
  }
  if( numGraphs != nGraphTotalSigma ) {
    cout << "tofcalgsec::calib_barrel_sigma: the number of Graphs is NOT reasonable!!!" << endl;
    exit(0);
  }

  m_name = string("calib_barrel_sigma");

  const int tbin   = 150;
  const double tbegin = -1.5;
  const double tend   =  1.5;

  // histograms per counter
  char hname[256];
  for( unsigned int i=0; i<NBarrel; i++ ) {
    m_result.push_back( HepVector(nBarrelSigma,0) );
    for( unsigned int j=0; j<nKind; j++ ) {
      for( unsigned int k=0; k<nBinPerCounter; k++ ) {
        if( j==0 )      { sprintf( hname, "tleft-id%i-z%i",  i, k); }
        else if( j==1 ) { sprintf( hname, "tright-id%i-z%i", i, k); }
        else if( j==2 ) { sprintf( hname, "t0-id%i-z%i",     i, k); }
        else if( j==3 ) { sprintf( hname, "tplus-id%i-z%i",  i, k); }
        else if( j==4 ) { sprintf( hname, "tminus-id%i-z%i", i, k); }
        m_histograms.push_back( new TH1F( hname, hname, tbin, tbegin, tend ) );

        m_fitresult.push_back( HepVector(nParSigma,0) );
      }
    }
  }

  zpos.resize( nBinPerCounter );
  zposerr.resize( nBinPerCounter );
  zstep = ( zend - zbegin )/nBinPerCounter;
  for( unsigned int i=0; i<nBinPerCounter; i++ ) {
    zpos[i] = zbegin + ( i+0.5 )*zstep;
    zposerr[i] = 0.5*zstep;
  }

}

calib_barrel_sigma::~calib_barrel_sigma

(

)

Definition at line 86 of file calib_barrel_sigma.cxx.

References m_fitresult, zpos, and zposerr.

                                        {
  m_fitresult.clear();
  zpos.clear();
  zposerr.clear();
}

---

Member Function Documentation

void calib_barrel_sigma::calculate	(	RecordSet *&	data,
		unsigned int	icounter
	)			[virtual]

Implements TofCalibFit.

Definition at line 93 of file calib_barrel_sigma.cxx.

References fillGraph(), fillGraphT0(), fillRecord(), fillRecordT0(), fitGraph(), fitGraphT0(), fitHistogram(), fitHistogramT0(), iter(), TofCalibFit::name(), and updateData().

                                                                            {

  std::cout << setiosflags(ios::left) << setw(10) << icounter << setw(8) << data->size() << setw(30) << name() << std::endl;

  if( data->size() > 0 ) {
    std::vector<Record*>::iterator iter = data->begin();
    for( ; iter!=data->end(); iter++ ) {
      fillRecord( (*iter), icounter );
    }
  }
  fitHistogram( icounter );
  fillGraph( icounter );
  fitGraph( icounter );

  if( data->size() > 0 ) {
    std::vector<Record*>::iterator iter = data->begin();
    for( ; iter!=data->end(); iter++ ) {
      updateData( (*iter), icounter );
      fillRecordT0( (*iter), icounter );
    }
  }
  fitHistogramT0( icounter );
  fillGraphT0( icounter );
  fitGraphT0( icounter );

  return;
}

void calib_barrel_sigma::fillGraph

(

unsigned int

icounter

)

[private]

Definition at line 170 of file calib_barrel_sigma.cxx.

References iter(), ganga-rec::j, m_fitresult, TofCalibFit::m_graphs, TofCalibFit::nBinPerCounter, TofCalibFit::nKind, zbegin, and zend.

Referenced by calculate().

                                                          {

  char gname1[256], gname2[256];

  // fill graphs
  // 4 canvas per counter,
  //   1. offset of tleft, tright and t0 vs z 
  //   2. sigma of tleft, tright and t0 vs z
  //   3. offset of tplus and tminus vs z
  //   4. sigma of tplus, tminus and T_Correlation vs z
  std::vector<double> toffset, toffseterr;
  std::vector<double> tsigma, tsigmaerr;
  toffset.resize( nBinPerCounter );
  toffseterr.resize( nBinPerCounter );
  tsigma.resize( nBinPerCounter );
  tsigmaerr.resize( nBinPerCounter );

  unsigned int number = 0;
  std::vector<HepVector>::iterator iter = m_fitresult.begin() + icounter*nKind*nBinPerCounter;
  for( unsigned int j=0; j<nKind; j++ ) {
    if( j==0 )      { sprintf( gname1, "tleft-offset-tofid-%i",  icounter ); }
    else if( j==1 ) { sprintf( gname1, "tright-offset-tofid-%i", icounter ); }
    else if( j==2 ) { sprintf( gname1, "t0-offset-tofid-%i",     icounter ); }
    else if( j==3 ) { sprintf( gname1, "tplus-offset-tofid-%i",  icounter ); }
    else if( j==4 ) { sprintf( gname1, "tminus-offset-tofid-%i", icounter ); }

    m_graphs.push_back( new TH1F( gname1, gname1, nBinPerCounter, zbegin, zend ) );
    std::vector<TH1F*>::iterator itgraph = m_graphs.end() - 1;

    for( unsigned int k=0; k<nBinPerCounter; k++ ) {
      number = j*nBinPerCounter + k;
      toffset[k]    = (*(iter+number))[0];
      toffseterr[k] = (*(iter+number))[1];
      (*itgraph)->SetBinContent( k+1, toffset[k]    );
      (*itgraph)->SetBinError(   k+1, toffseterr[k] );
    }
    (*itgraph)->SetMarkerSize(1.5);
    if( j==0 || j==3) {
      (*itgraph)->SetMarkerStyle(20);
      (*itgraph)->SetMarkerColor(2);
      (*itgraph)->SetMaximum( 0.15 );
      (*itgraph)->SetMinimum(-0.15 );
    }
    else if( j==1 || j==4 ) {
      (*itgraph)->SetMarkerStyle(21);
      (*itgraph)->SetMarkerColor(4);
    }
    else if( j==2 ) {
      (*itgraph)->SetMarkerStyle(4);
      (*itgraph)->SetMarkerColor(2);
    }
  }
    
  for( unsigned int j=0; j<nKind; j++ ) {
    if( j==0 )      { sprintf( gname2, "tleft-sigma-tofid-%i",  icounter ); }
    else if( j==1 ) { sprintf( gname2, "tright-sigma-tofid-%i", icounter ); }
    else if( j==2 ) { sprintf( gname2, "t0-sigma-tofid-%i",     icounter ); }
    else if( j==3 ) { sprintf( gname2, "tplus-sigma-tofid-%i",  icounter ); }
    else if( j==4 ) { sprintf( gname2, "tminus-sigma-tofid-%i", icounter ); }
    m_graphs.push_back( new TH1F( gname2, gname2, nBinPerCounter, zbegin, zend ) );
    std::vector<TH1F*>::iterator itgraph = m_graphs.end() - 1;

    for( unsigned int k=0; k<nBinPerCounter; k++ ) {
      number = j*nBinPerCounter + k;
      tsigma[k]     = (*(iter+number))[2];
      tsigmaerr[k]  = (*(iter+number))[3];
      (*itgraph)->SetBinContent( k+1, tsigma[k]    );
      (*itgraph)->SetBinError(   k+1, tsigmaerr[k] );
    }
    (*itgraph)->SetMarkerSize(1.5);
    if( j==0 || j==3 ) {
      (*itgraph)->SetMarkerStyle(20);
      (*itgraph)->SetMarkerColor(2);
      (*itgraph)->SetMaximum( 0.3 );
      (*itgraph)->SetMinimum( 0.0 );
    }
    else if( j==1 || j==4 ) {
      (*itgraph)->SetMarkerStyle(21);
      (*itgraph)->SetMarkerColor(4);
    }
    else if( j==2 ) {
      (*itgraph)->SetMarkerStyle(4);
      (*itgraph)->SetMarkerColor(2);
    }
  }

  sprintf( gname2, "sigma-tofid-%i", icounter );
  m_graphs.push_back( new TH1F( gname2, gname2, nBinPerCounter, zbegin, zend ) );
  std::vector<TH1F*>::iterator itgraph = m_graphs.end() - 1;
  for( unsigned int k=0; k<nBinPerCounter; k++ ) {
    number = (nKind-1)*nBinPerCounter + k;
    double sigPlus  = (*(iter+number-nBinPerCounter))[2];
    double sigMinus = (*(iter+number))[2];
    double sigErrPlus  = (*(iter+number-nBinPerCounter))[3];
    double sigErrMinus = (*(iter+number))[3];
    tsigma[k]     = sqrt( sigPlus*sigPlus - sigMinus*sigMinus );
    tsigmaerr[k]  = sqrt( sigErrPlus*sigErrPlus + sigErrMinus*sigErrMinus );
    (*itgraph)->SetBinContent( k+1, tsigma[k]    );
    (*itgraph)->SetBinError(   k+1, tsigmaerr[k] );
  }
  (*itgraph)->SetMarkerSize(1.5);
  (*itgraph)->SetMarkerStyle(4);
  (*itgraph)->SetMarkerColor(2);

  return;
}

void calib_barrel_sigma::fillGraphT0

(

unsigned int

icounter

)

[private]

Definition at line 372 of file calib_barrel_sigma.cxx.

References genRecEmupikp::i, iter(), m_fitresult, TofCalibFit::m_graphs, TofCalibFit::nBinPerCounter, nGraphTotalSigma, and TofCalibFit::nKind.

Referenced by calculate().

                                                            {
  char gname1[256], gname2[256];

  // fill graphs
  // 4 canvas per counter,
  //   1. offset of tleft, tright and t0 vs z 
  //   2. sigma of tleft, tright and t0 vs z
  //   3. offset of tplus and tminus vs z
  //   4. sigma of tplus, tminus and T_Correlation vs z
  std::vector<double> toffset, toffseterr;
  std::vector<double> tsigma, tsigmaerr;
  toffset.resize( nBinPerCounter );
  toffseterr.resize( nBinPerCounter );
  tsigma.resize( nBinPerCounter );
  tsigmaerr.resize( nBinPerCounter );

  std::vector<HepVector>::iterator iter = m_fitresult.begin() + icounter*nKind*nBinPerCounter + 2*nBinPerCounter;
  for( unsigned int k=0; k<nBinPerCounter; k++ ) {
    toffset[k]    = (*(iter+k))[0];
    toffseterr[k] = (*(iter+k))[1];
    tsigma[k]     = (*(iter+k))[2];
    tsigmaerr[k]  = (*(iter+k))[3];
  }

  sprintf( gname1, "offset-tofid-%i", icounter );
  std::vector<TH1F*>::iterator itgraph = m_graphs.begin() + icounter*nGraphTotalSigma + 2;
  for( unsigned int i=0; i<nBinPerCounter; i++ ) {
    //    (*itgraph)->SetPoint( i, zpos[i], toffset[i] );
    //    (*itgraph)->SetPointError( i, zposerr[i], toffseterr[i] );
    (*itgraph)->SetBinContent( i+1, toffset[i]    );
    (*itgraph)->SetBinError(   i+1, toffseterr[i] );
  }

  sprintf( gname2, "sigma-tofid-%i", icounter );
  itgraph = m_graphs.begin() + icounter*nGraphTotalSigma + 7;
  for( unsigned int i=0; i<nBinPerCounter; i++ ) {
    //    (*itgraph)->SetPoint( i, zpos[i], tsigma[i] );
    //    (*itgraph)->SetPointError( i, zposerr[i], tsigmaerr[i] );
    (*itgraph)->SetBinContent( i+1, tsigma[i]    );
    (*itgraph)->SetBinError(   i+1, tsigmaerr[i] );
  }

  return;
}

void calib_barrel_sigma::fillRecord	(	const Record *	r,
		unsigned int	icounter
	)			[private]

Definition at line 122 of file calib_barrel_sigma.cxx.

References iter(), TofCalibFit::m_histograms, TofCalibFit::nBinPerCounter, TofCalibFit::nKind, Record::tleft(), Record::tright(), zbegin, zend, Record::zrhit(), and zstep.

Referenced by calculate().

                                                                            {

  double zhit = r->zrhit();
  if( zhit<zbegin || zhit>zend ) return;
  int zbin = static_cast<int>((zhit-zbegin)/zstep);
  if( zbin<0 ) { zbin = 0; }
  else if( zbin>static_cast<int>(nBinPerCounter-1) ) {
    cout << "tofcalgsec::calib_barrel_sigma:fillRecord: zhit is out of range, zhit=" << zhit << " zbin=" << zbin << endl;
    return;
  }

  std::vector<TH1F*>::iterator iter = m_histograms.begin();
  unsigned int number = icounter*nKind*nBinPerCounter + zbin;
  (*(iter+number))->Fill( r->tleft() );
  (*(iter+nBinPerCounter+number))->Fill( r->tright() );
  (*(iter+3*nBinPerCounter+number))->Fill( (r->tleft()+r->tright())/2.0 );
  (*(iter+4*nBinPerCounter+number))->Fill( (r->tleft()-r->tright())/2.0 );

  return;
}

void calib_barrel_sigma::fillRecordT0	(	const Record *	r,
		unsigned int	icounter
	)			[private]

Definition at line 335 of file calib_barrel_sigma.cxx.

References iter(), TofCalibFit::m_histograms, TofCalibFit::nBinPerCounter, TofCalibFit::nKind, Record::t0(), zbegin, zend, Record::zrhit(), and zstep.

Referenced by calculate().

                                                                              {
  double zhit = r->zrhit();
  if( zhit<zbegin || zhit>zend ) return;
  int zbin = static_cast<int>((zhit-zbegin)/zstep);
  if( zbin<0 ) { zbin = 0; }
  else if( zbin>static_cast<int>(nBinPerCounter-1) ) {
    cout << "tofcalgsec::calib_barrel_sigma:fillRecordT0: zhit is out of range, zhit=" << zhit << " zbin=" << zbin << endl;
    return;
  }

  std::vector<TH1F*>::iterator iter = m_histograms.begin();
  unsigned int number = icounter*nKind*nBinPerCounter + 2*nBinPerCounter + zbin;
  (*(iter+number))->Fill( r->t0() );
 
  return;
}

void TofCalibFit::fillRoot

(

const char *

file

)

[inherited]

Definition at line 80 of file TofCalibFit.cxx.

References TofCalibFit::CanvasName, TofCalibFit::CanvasPerCounterName, delete_small_size::dirname, genRecEmupikp::i, iter(), ganga-rec::j, TofCalibFit::m_graphs, TofCalibFit::m_histograms, TofCalibFit::m_name, TofCalibFit::nCanvas, TofCalibFit::nCanvasPerCounter, TofCalibFit::nCounter, TofCalibFit::nGraphPerCanvas, TofCalibFit::nGraphPerCanvasPerCounter, TofCalibFit::nHistogram, and TofCalibFit::nHistPerCounter.

                                             {

  unsigned int nhist = m_histograms.size();
  if( nhist != (nCounter*nHistPerCounter + nHistogram) ) {
    std::cout<<" tofcalgsec::TofCalibFit:" << m_name << ": the number of histograms is NOT same as the number of histograms saved!" << " nhist=" << nhist << " calculated=" << (nCounter*nHistPerCounter + nHistogram) << " nCounter=" << nCounter << " nHistPerCounter=" << nHistPerCounter << " nHistogram=" << nHistogram << std::endl;
    exit(0);
  }

  unsigned int numgraph1 = 0;
  unsigned int numgraph2 = 0;
  if( nCanvasPerCounter!=0 ) {
    std::vector<unsigned int>::iterator iter = nGraphPerCanvasPerCounter.begin();
    for( ; iter!=nGraphPerCanvasPerCounter.end(); iter++ ) {
      numgraph1 = numgraph1 + (*iter);
    }
  }
  if( nCanvas!=0 ) {
    std::vector<unsigned int>::iterator iter = nGraphPerCanvas.begin();
    for( ; iter!=nGraphPerCanvas.end(); iter++ ) {
      numgraph2 = numgraph2 + (*iter);
    }
  }
  unsigned int ngraph = m_graphs.size();
  if( ngraph != ( nCounter*numgraph1+numgraph2) ) {
    std::cout<<" tofcalgsec::TofCalibFit:"<< m_name << ": the number of graphs is NOT same as the number of graphs saved!"<<std::endl;
    exit(0);
  }

  TFile f(file,"RECREATE");

  gStyle->SetOptStat(2211);
  gStyle->SetOptFit(1111);
  gStyle->SetLabelSize(0.03,"x");
  gStyle->SetLabelSize(0.03,"y");

  char dirname[256];
  char canvasname[256];
  if( nHistPerCounter>0 || nCanvasPerCounter>0 ) {
    std::vector<TH1F*>::iterator iter1 = m_histograms.begin();
    std::vector<TH1F*>::iterator iter2 = m_graphs.begin();
    for( unsigned int i=0; i<nCounter; i++ ) {
      sprintf( dirname, "tofid%i", i );
      TDirectory* cdresult = f.mkdir( dirname );
      cdresult->cd();
      
      for( unsigned int j=0; j<nHistPerCounter; j++ ) {
        (*(iter1+j))->Write();
      }
      iter1 = iter1 + nHistPerCounter;

      for( unsigned int j=0; j<nCanvasPerCounter; j++ ) {
        std::vector<string>::iterator it1 = CanvasPerCounterName.begin() + j;
        std::vector<unsigned int>::iterator it2 = nGraphPerCanvasPerCounter.begin() + j;
        sprintf( canvasname, "%s-tofid-%i", (*it1).c_str(), i );
        TCanvas* c1 = new TCanvas( canvasname, canvasname, 1);
        c1->SetFillColor(10);
        for( unsigned int k=0; k<(*it2); k++ ) {
          if( k==0 ) {
            (*(iter2+k))->Draw("E");
          }
          else {
            (*(iter2+k))->Draw("Esame");
          }
        }
        iter2 = iter2 + (*it2);
        c1->Write();
      }
    }
  }

  if( nHistogram>0 || nCanvas>0 ) {
    sprintf( dirname, "summary" );
    TDirectory* cdresult = f.mkdir( dirname );
    cdresult->cd();

    std::vector<TH1F*>::iterator iter1 = m_histograms.begin() + nCounter*nHistPerCounter;
    std::vector<TH1F*>::iterator iter2 = m_graphs.begin() + nCounter*numgraph1;
    for( ; iter1 != m_histograms.end(); iter1++ ) {
      (*iter1)->Write();
    }

    for( unsigned int j=0; j<nCanvas; j++ ) {
      std::vector<string>::iterator it1 = CanvasName.begin() + j;
      std::vector<unsigned int>::iterator it2 = nGraphPerCanvas.begin() + j;
      sprintf( canvasname, (*it1).c_str() );
      TCanvas* c1 = new TCanvas( canvasname, canvasname, 1);
      c1->SetFillColor(10);
      for( unsigned int k=0; k<(*it2); k++ ) {
        if( k==0 ) {
          (*(iter2+k))->Draw("E");
        }
        else {
          (*(iter2+k))->Draw("Esame");
        }
      }
      iter2 = iter2 + (*it2);
      c1->Write();
    }

  }

  f.Close();

  return;
}

void TofCalibFit::fillTxt

(

const char *

file

)

[inherited]

Definition at line 52 of file TofCalibFit.cxx.

References genRecEmupikp::i, TofCalibFit::m_result, and TofCalibFit::name().

                                            {
  std::ofstream out(file,ios::out);
  if( out ) {
    std::vector<HepVector>::iterator it;
    for( it=m_result.begin(); it!=m_result.end(); it++ ) {
      for( int i=0; i<(*it).num_row(); i++ ) {
        out << (*it)[i] << " ";
      }
      out << std::endl;
    }
    out.close();
  }
  else{
    cerr << "error when open file " << file << " for write in " << name() << "::fillTxt()" << std::endl;
    cout << "print all parameters to srceen: in total " << m_result.size() << " items" << std::endl;
    std::vector<HepVector>::iterator it;
    for( it=m_result.begin(); it!=m_result.end(); it++ ) {
      for( int i=0; i<(*it).num_row(); i++ ) {
        cout << (*it)[i] << " ";
      }
      cout << std::endl;
    }
  }

  return;
}

void calib_barrel_sigma::fitGraph

(

unsigned int

icounter

)

[private]

Definition at line 278 of file calib_barrel_sigma.cxx.

References genRecEmupikp::i, iter(), TofCalibFit::m_graphs, TofCalibFit::m_npar, TofCalibFit::m_result, TofCalibFit::nGraphPerCanvasPerCounter, nGraphTotalSigma, TofCalibFit::X, zbegin, and zend.

Referenced by calculate().

                                                         {

  TF1* fsingle = new TF1("fsingle", "pol4");
  fsingle->SetLineColor(1);
  fsingle->SetLineWidth(1);

  std::vector<unsigned int>::iterator itnumber = nGraphPerCanvasPerCounter.begin();
  std::vector<TH1F*>::iterator        itgraph  = m_graphs.begin() + icounter*nGraphTotalSigma + (*itnumber) + (*(itnumber+1));

  fsingle->SetParameter( 0, 0.14    );
  fsingle->SetParameter( 1, -4.0e-4 );
  (*itgraph)->Fit( "fsingle", "QR", "", zbegin, zend );
  X = HepVector( m_npar, 0 );
  for( unsigned int i=0; i<5; i++ ) {
    X[i] = fsingle->GetParameter(i);
  }

  fsingle->SetParameter( 0, 0.14   );
  fsingle->SetParameter( 1, 4.0e-4 );
  (*(itgraph+1))->Fit( "fsingle", "QR", "", zbegin, zend );
  for( unsigned int i=0; i<5; i++ ) {
    X[i+5] = fsingle->GetParameter(i);
  }

  std::vector<HepVector>::iterator iter = m_result.begin() + icounter;
  (*iter) = X;

  return;
}

void calib_barrel_sigma::fitGraphT0

(

unsigned int

icounter

)

[private]

Definition at line 418 of file calib_barrel_sigma.cxx.

References iter(), TofCalibFit::m_graphs, TofCalibFit::m_result, TofCalibFit::nGraphPerCanvasPerCounter, nGraphTotalSigma, zbegin, and zend.

Referenced by calculate().

                                                           {

  //  TF1 *fdouble = new TF1( "fdouble", doubleEndFunc, zbegin, zend, 3 );
  TF1 *fdouble = new TF1( "fdouble", "pol4", zbegin, zend );
  fdouble->SetLineColor(1);
  fdouble->SetLineWidth(1);

  std::vector<unsigned int>::iterator itnumber = nGraphPerCanvasPerCounter.begin();
  std::vector<TH1F*>::iterator itgraph = m_graphs.begin() + icounter*nGraphTotalSigma + (*itnumber) + (*(itnumber+1)) + 2;
  (*itgraph)->Fit( "fdouble", "Q", "", zbegin, zend );

  std::vector<HepVector>::iterator iter = m_result.begin() + icounter;
  (*iter)[10] = fdouble->GetParameter(0);
  (*iter)[11] = fdouble->GetParameter(1);
  (*iter)[12] = fdouble->GetParameter(2);
  (*iter)[13] = fdouble->GetParameter(3);
  (*iter)[14] = fdouble->GetParameter(4);

  return;
}

void calib_barrel_sigma::fitHistogram

(

unsigned int

icounter

)

[private]

Definition at line 144 of file calib_barrel_sigma.cxx.

References genRecEmupikp::i, ganga-rec::j, m_fitresult, TofCalibFit::m_histograms, TofCalibFit::nBinPerCounter, and TofCalibFit::nKind.

Referenced by calculate().

                                                             {
  TF1* g = new TF1("g", "gaus");
  g->SetLineColor(2);
  g->SetLineWidth(1);

  std::vector<TH1F*>::iterator     iter1 = m_histograms.begin() + icounter*nKind*nBinPerCounter;
  std::vector<HepVector>::iterator iter2 = m_fitresult.begin() + icounter*nKind*nBinPerCounter;
  for( unsigned int i=0; i<nKind; i++ ) {
    for( unsigned int j=0; j<nBinPerCounter; j++ ) {
      if( i!=2 ) {
        (*iter1)->Fit( g, "Q");
        (*iter2)[0] = g->GetParameter(1);
        (*iter2)[1] = g->GetParError(1);
        (*iter2)[2] = g->GetParameter(2);
        (*iter2)[3] = g->GetParError(2);
      }
      iter1++;
      iter2++;
    }
  }

  return;

}

void calib_barrel_sigma::fitHistogramT0

(

unsigned int

icounter

)

[private]

Definition at line 353 of file calib_barrel_sigma.cxx.

References ganga-rec::j, m_fitresult, TofCalibFit::m_histograms, TofCalibFit::nBinPerCounter, and TofCalibFit::nKind.

Referenced by calculate().

                                                               {
  TF1* g = new TF1("g", "gaus");
  g->SetLineColor(2);
  g->SetLineWidth(1);

  std::vector<TH1F*>::iterator     iter1 = m_histograms.begin() + icounter*nKind*nBinPerCounter + 2*nBinPerCounter;
  std::vector<HepVector>::iterator iter2 = m_fitresult.begin() + icounter*nKind*nBinPerCounter + 2*nBinPerCounter;
  for( unsigned int j=0; j<nBinPerCounter; j++, iter1++, iter2++ ) {
    (*iter1)->Fit( g, "Q");
    (*iter2)[0] = g->GetParameter(1);
    (*iter2)[1] = g->GetParError(1);
    (*iter2)[2] = g->GetParameter(2);
    (*iter2)[3] = g->GetParError(2);
  }

  return;
}

const string& TofCalibFit::name

(

void

)

const [inline, inherited]

Definition at line 27 of file TofCalibFit.h.

References TofCalibFit::m_name.

Referenced by TofCalibManager::addCalib(), calib_etf_bunch::calculate(), calib_endcap_sigma::calculate(), calib_endcap_atten::calculate(), calculate(), calib_barrel_q0::calculate(), calib_barrel_common::calculate(), and TofCalibFit::fillTxt().

{ return m_name; }

void TofCalibFit::setTCorrelation

(

HepVector

tc

)

[inline, inherited]

Definition at line 35 of file TofCalibFit.h.

References TofCalibFit::m_tcorrelation.

{ m_tcorrelation = tc; }

HepVector TofCalibFit::tcorrelation

(

)

[inline, inherited]

Definition at line 34 of file TofCalibFit.h.

References TofCalibFit::X.

{ return X;            }

void calib_barrel_sigma::updateData	(	Record *	r,
		unsigned int	icounter
	)			[private]

Definition at line 309 of file calib_barrel_sigma.cxx.

References genRecEmupikp::i, iter(), TofCalibFit::m_result, TofCalibFit::m_tcorrelation, Record::setT0(), Record::tleft(), Record::tright(), and Record::zrhit().

Referenced by calculate().

                                                                      {
  double zhit = r->zrhit();
  double t1 = r->tleft();
  double t2 = r->tright();

  double par1[5], par2[5];
  std::vector<HepVector>::iterator iter = m_result.begin() + icounter;
  for( unsigned int i=0; i<5; i++ ) {
    par1[i] = (*iter)[i];
    par2[i] = (*iter)[i+5];
  }

  double tsigma1 = par1[0]+par1[1]*zhit+par1[2]*pow(zhit,2)+par1[3]*pow(zhit,3) + par1[4]*pow(zhit,4);
  double tsigma2 = par2[0]+par2[1]*zhit+par2[2]*pow(zhit,2)+par2[3]*pow(zhit,3) + par2[4]*pow(zhit,4);
  double tc      = m_tcorrelation[0];

  double weight1 = (tsigma2*tsigma2-tc*tc)/(tsigma1*tsigma1+tsigma2*tsigma2-2.0*tc*tc);
  double weight2 = (tsigma1*tsigma1-tc*tc)/(tsigma1*tsigma1+tsigma2*tsigma2-2.0*tc*tc);
  double t0 = weight1*t1 + weight2*t2;

  r->setT0( t0 );
  
  return;
}

---

Member Data Documentation

std::vector<string> TofCalibFit::CanvasName [protected, inherited]

Definition at line 58 of file TofCalibFit.h.

Referenced by calib_barrel_common::calib_barrel_common(), calib_barrel_q0::calib_barrel_q0(), calib_endcap_atten::calib_endcap_atten(), calib_etf_bunch::calib_etf_bunch(), calib_barrel_q0::fillGraph(), and TofCalibFit::fillRoot().

std::vector<string> TofCalibFit::CanvasPerCounterName [protected, inherited]

Definition at line 57 of file TofCalibFit.h.

Referenced by calib_barrel_sigma(), calib_endcap_atten::calib_endcap_atten(), calib_endcap_sigma::calib_endcap_sigma(), and TofCalibFit::fillRoot().

std::vector<HepVector> calib_barrel_sigma::m_fitresult [private]

Definition at line 34 of file calib_barrel_sigma.h.

Referenced by calib_barrel_sigma(), fillGraph(), fillGraphT0(), fitHistogram(), fitHistogramT0(), and ~calib_barrel_sigma().

std::vector<TH1F*> TofCalibFit::m_graphs [protected, inherited]

Definition at line 54 of file TofCalibFit.h.

Referenced by calib_etf_bunch::fillGraph(), calib_endcap_sigma::fillGraph(), calib_endcap_atten::fillGraph(), fillGraph(), calib_barrel_q0::fillGraph(), calib_barrel_common::fillGraph(), calib_endcap_atten::fillGraphQ0(), fillGraphT0(), TofCalibFit::fillRoot(), calib_etf_bunch::fitGraph(), calib_endcap_sigma::fitGraph(), calib_endcap_atten::fitGraph(), fitGraph(), calib_barrel_common::fitGraph(), fitGraphT0(), and TofCalibFit::~TofCalibFit().

std::vector<TH1F*> TofCalibFit::m_histograms [protected, inherited]

Definition at line 53 of file TofCalibFit.h.

Referenced by calib_barrel_common::calib_barrel_common(), calib_barrel_q0::calib_barrel_q0(), calib_barrel_sigma(), calib_endcap_atten::calib_endcap_atten(), calib_endcap_sigma::calib_endcap_sigma(), calib_etf_bunch::calib_etf_bunch(), calib_etf_bunch::fillRecord(), calib_endcap_sigma::fillRecord(), calib_endcap_atten::fillRecord(), fillRecord(), calib_barrel_q0::fillRecord(), calib_barrel_common::fillRecord(), calib_endcap_atten::fillRecordQ0(), fillRecordT0(), TofCalibFit::fillRoot(), calib_etf_bunch::fitHistogram(), calib_endcap_sigma::fitHistogram(), calib_endcap_atten::fitHistogram(), fitHistogram(), calib_barrel_q0::fitHistogram(), calib_barrel_common::fitHistogram(), calib_endcap_atten::fitHistogramQ0(), fitHistogramT0(), and TofCalibFit::~TofCalibFit().

string TofCalibFit::m_name [protected, inherited]

Definition at line 50 of file TofCalibFit.h.

Referenced by calib_barrel_common::calib_barrel_common(), calib_barrel_q0::calib_barrel_q0(), calib_barrel_sigma(), calib_endcap_atten::calib_endcap_atten(), calib_endcap_sigma::calib_endcap_sigma(), calib_etf_bunch::calib_etf_bunch(), TofCalibFit::fillRoot(), TofCalibFit::name(), and TofCalibFit::TofCalibFit().

int TofCalibFit::m_npar [protected, inherited]

Definition at line 38 of file TofCalibFit.h.

Referenced by calib_endcap_sigma::fitGraph(), calib_endcap_atten::fitGraph(), fitGraph(), and TofCalibFit::TofCalibFit().

std::vector<HepVector> TofCalibFit::m_result [protected, inherited]

Definition at line 55 of file TofCalibFit.h.

Referenced by calib_barrel_q0::calib_barrel_q0(), calib_barrel_sigma(), calib_endcap_atten::calib_endcap_atten(), calib_endcap_sigma::calib_endcap_sigma(), calib_endcap_atten::fillGraphQ0(), TofCalibFit::fillTxt(), calib_etf_bunch::fitGraph(), calib_endcap_sigma::fitGraph(), calib_endcap_atten::fitGraph(), fitGraph(), calib_barrel_q0::fitGraph(), calib_barrel_common::fitGraph(), fitGraphT0(), calib_endcap_atten::updateData(), updateData(), and TofCalibFit::~TofCalibFit().

HepVector TofCalibFit::m_tcorrelation [protected, inherited]

Definition at line 60 of file TofCalibFit.h.

Referenced by TofCalibFit::setTCorrelation(), TofCalibFit::TofCalibFit(), and updateData().

unsigned int TofCalibFit::nBinPerCounter [protected, inherited]

Definition at line 41 of file TofCalibFit.h.

Referenced by calib_barrel_common::calib_barrel_common(), calib_barrel_q0::calib_barrel_q0(), calib_barrel_sigma(), calib_endcap_atten::calib_endcap_atten(), calib_endcap_sigma::calib_endcap_sigma(), calib_etf_bunch::calib_etf_bunch(), calib_etf_bunch::fillGraph(), calib_endcap_sigma::fillGraph(), calib_endcap_atten::fillGraph(), fillGraph(), calib_barrel_common::fillGraph(), calib_endcap_atten::fillGraphQ0(), fillGraphT0(), calib_etf_bunch::fillRecord(), calib_endcap_sigma::fillRecord(), calib_endcap_atten::fillRecord(), fillRecord(), calib_barrel_common::fillRecord(), calib_endcap_atten::fillRecordQ0(), fillRecordT0(), calib_endcap_sigma::fitHistogram(), calib_endcap_atten::fitHistogram(), fitHistogram(), calib_endcap_atten::fitHistogramQ0(), fitHistogramT0(), and TofCalibFit::TofCalibFit().

unsigned int TofCalibFit::nCanvas [protected, inherited]

Definition at line 47 of file TofCalibFit.h.

Referenced by calib_barrel_common::calib_barrel_common(), calib_barrel_q0::calib_barrel_q0(), calib_barrel_sigma(), calib_endcap_atten::calib_endcap_atten(), calib_endcap_sigma::calib_endcap_sigma(), calib_etf_bunch::calib_etf_bunch(), TofCalibFit::fillRoot(), and TofCalibFit::TofCalibFit().

unsigned int TofCalibFit::nCanvasPerCounter [protected, inherited]

Definition at line 44 of file TofCalibFit.h.

Referenced by calib_barrel_common::calib_barrel_common(), calib_barrel_q0::calib_barrel_q0(), calib_barrel_sigma(), calib_endcap_atten::calib_endcap_atten(), calib_endcap_sigma::calib_endcap_sigma(), calib_etf_bunch::calib_etf_bunch(), TofCalibFit::fillRoot(), and TofCalibFit::TofCalibFit().

std::vector<unsigned int> TofCalibFit::nGraphPerCanvas [protected, inherited]

Definition at line 48 of file TofCalibFit.h.

Referenced by calib_barrel_common::calib_barrel_common(), calib_barrel_q0::calib_barrel_q0(), calib_endcap_atten::calib_endcap_atten(), calib_etf_bunch::calib_etf_bunch(), and TofCalibFit::fillRoot().

std::vector<unsigned int> TofCalibFit::nGraphPerCanvasPerCounter [protected, inherited]

Definition at line 45 of file TofCalibFit.h.

Referenced by calib_barrel_sigma(), calib_endcap_atten::calib_endcap_atten(), calib_endcap_sigma::calib_endcap_sigma(), TofCalibFit::fillRoot(), fitGraph(), and fitGraphT0().

unsigned int TofCalibFit::nHistogram [protected, inherited]

Definition at line 46 of file TofCalibFit.h.

Referenced by calib_barrel_common::calib_barrel_common(), calib_barrel_q0::calib_barrel_q0(), calib_barrel_sigma(), calib_endcap_atten::calib_endcap_atten(), calib_endcap_sigma::calib_endcap_sigma(), calib_etf_bunch::calib_etf_bunch(), TofCalibFit::fillRoot(), and TofCalibFit::TofCalibFit().

unsigned int TofCalibFit::nHistPerCounter [protected, inherited]

Definition at line 43 of file TofCalibFit.h.

Referenced by calib_barrel_common::calib_barrel_common(), calib_barrel_q0::calib_barrel_q0(), calib_barrel_sigma(), calib_endcap_atten::calib_endcap_atten(), calib_endcap_sigma::calib_endcap_sigma(), calib_etf_bunch::calib_etf_bunch(), TofCalibFit::fillRoot(), and TofCalibFit::TofCalibFit().

unsigned int TofCalibFit::nKind [protected, inherited]

Definition at line 40 of file TofCalibFit.h.

Referenced by calib_barrel_common::calib_barrel_common(), calib_barrel_q0::calib_barrel_q0(), calib_barrel_sigma(), calib_endcap_atten::calib_endcap_atten(), calib_endcap_sigma::calib_endcap_sigma(), calib_etf_bunch::calib_etf_bunch(), calib_etf_bunch::fillGraph(), calib_endcap_sigma::fillGraph(), calib_endcap_atten::fillGraph(), fillGraph(), calib_barrel_common::fillGraph(), fillGraphT0(), calib_endcap_sigma::fillRecord(), calib_endcap_atten::fillRecord(), fillRecord(), calib_barrel_common::fillRecord(), calib_endcap_atten::fillRecordQ0(), fillRecordT0(), calib_endcap_sigma::fitHistogram(), calib_endcap_atten::fitHistogram(), fitHistogram(), calib_endcap_atten::fitHistogramQ0(), fitHistogramT0(), and TofCalibFit::TofCalibFit().

HepVector TofCalibFit::X [protected, inherited]

Definition at line 51 of file TofCalibFit.h.

Referenced by calib_etf_bunch::fitGraph(), calib_endcap_sigma::fitGraph(), calib_endcap_atten::fitGraph(), fitGraph(), calib_barrel_common::fitGraph(), calib_barrel_common::fitHistogram(), TofCalibFit::tcorrelation(), and TofCalibFit::TofCalibFit().

std::vector<double> calib_barrel_sigma::zpos [private]

Definition at line 30 of file calib_barrel_sigma.h.

Referenced by calib_barrel_sigma(), and ~calib_barrel_sigma().

std::vector<double> calib_barrel_sigma::zposerr [private]

Definition at line 31 of file calib_barrel_sigma.h.

Referenced by calib_barrel_sigma(), and ~calib_barrel_sigma().

double calib_barrel_sigma::zstep [private]

Definition at line 32 of file calib_barrel_sigma.h.

Referenced by calib_barrel_sigma(), fillRecord(), and fillRecordT0().

---