#include <TofCalibFit.h>

Inheritance diagram for TofCalibFit:

Public Member Functions
	TofCalibFit (bool isbarrel, const int npar)

	~TofCalibFit ()

const string &	name () const

virtual void	calculate (RecordSet *&data, unsigned int icounter)=0

void	fillTxt (const char *file)

void	fillRoot (const char *file)

HepVector	tcorrelation ()

void	setTCorrelation (HepVector tc)

	TofCalibFit (bool isbarrel, const int npar)

	~TofCalibFit ()

const string &	name () const

virtual void	calculate (RecordSet *&data, unsigned int icounter)=0

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

Detailed Description

Definition at line 21 of file InstallArea/include/tofcalgsec/tofcalgsec/TofCalibFit.h.

Constructor & Destructor Documentation

◆ TofCalibFit() [1/2]

TofCalibFit::TofCalibFit	(	bool	isbarrel,
		const int	npar
	)

Definition at line 17 of file TofCalibFit.cxx.

                                                        {
  if( isbarrel ) { nCounter = NBarrel; }
  else { nCounter = NEndcap; }
 
  nKind          = 0;
  nBinPerCounter = 0;
 
  nHistPerCounter   = 0;
  nCanvasPerCounter = 0;
  nHistogram        = 0;
  nCanvas           = 0;
 
  m_npar = npar;
  X = HepVector(m_npar,0);
  m_name = string("calibfit");
 
  m_tcorrelation = HepVector( nBarrelCommon,0);
}

◆ ~TofCalibFit() [1/2]

TofCalibFit::~TofCalibFit ( )

Definition at line 37 of file TofCalibFit.cxx.

                          {
  std::vector<TH1F*>::iterator iter1 = m_histograms.begin();
  for( ; iter1 != m_histograms.end(); iter1++ ) {
    delete (*iter1);
  }
  m_histograms.clear();
  std::vector<TH1F*>::iterator iter2 = m_graphs.begin();
  for( ; iter2 != m_graphs.end(); iter2++ ) {
    delete (*iter2);
  }
  m_graphs.clear();
  m_result.clear();
}

◆ TofCalibFit() [2/2]

TofCalibFit::TofCalibFit	(	bool	isbarrel,
		const int	npar
	)

◆ ~TofCalibFit() [2/2]

TofCalibFit::~TofCalibFit ( )

Member Function Documentation

◆ calculate() [1/2]

virtual void TofCalibFit::calculate	(	RecordSet *&	data,
		unsigned int	icounter
	)

pure virtual

Implemented in calib_etf_bunch, calib_etf_bunch, calib_barrel_common, calib_barrel_q0, calib_barrel_sigma, calib_endcap_atten, calib_endcap_sigma, calib_barrel_common, calib_barrel_q0, calib_barrel_sigma, calib_endcap_atten, and calib_endcap_sigma.

◆ calculate() [2/2]

virtual void TofCalibFit::calculate	(	RecordSet *&	data,
		unsigned int	icounter
	)

pure virtual

Implemented in calib_etf_bunch, calib_etf_bunch, calib_barrel_common, calib_barrel_q0, calib_barrel_sigma, calib_endcap_atten, calib_endcap_sigma, calib_barrel_common, calib_barrel_q0, calib_barrel_sigma, calib_endcap_atten, and calib_endcap_sigma.

◆ fillRoot() [1/2]

void TofCalibFit::fillRoot ( const char * file )

Definition at line 80 of file TofCalibFit.cxx.

                                             {
 
  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;
}

◆ fillRoot() [2/2]

void TofCalibFit::fillRoot ( const char * file )

◆ fillTxt() [1/2]

void TofCalibFit::fillTxt ( const char * file )

Definition at line 52 of file TofCalibFit.cxx.

                                            {
  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;
}