#include <TofCalibration.h>

Public Member Functions
	TofCalibration (const int npar)

	~TofCalibration ()

virtual const string &	name () const

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

void	fillTxt (const char *file)

void	fillRoot (const char *file, bool isbarrel)

	TofCalibration (const int npar)

	~TofCalibration ()

virtual const string &	name () const

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

void	fillTxt (const char *file)

void	fillRoot (const char *file, bool isbarrel)

Protected Member Functions
virtual void	calculate_record (const Record *r, unsigned int icounter)=0

virtual void	calculate_result (unsigned int icounter)=0

virtual void	calculate_record (const Record *r, unsigned int icounter)=0

virtual void	calculate_result (unsigned int icounter)=0

Protected Attributes
string	m_name

int	Npar

unsigned int	Npar2

unsigned int	nHistPerCounter

unsigned int	nGraphPerCounter

unsigned int	nHistogram

unsigned int	nGraph

HepVector	X

std::vector< TH1F * >	histograms

std::vector< TGraph * >	graphs

std::vector< HepVector >	m_result

Detailed Description

Definition at line 13 of file InstallArea/include/tofcalgsec/tofcalgsec/TofCalibration.h.

Constructor & Destructor Documentation

◆ TofCalibration() [1/2]

TofCalibration::TofCalibration ( const int npar )

Definition at line 11 of file TofCalibration.cxx.

                                              :Npar(npar) {
  X = HepVector(Npar,0);
  m_name = string("calibration");
}

◆ ~TofCalibration() [1/2]

TofCalibration::~TofCalibration ( )

Definition at line 16 of file TofCalibration.cxx.

                                {
  std::vector<TH1F*>::iterator iter = histograms.begin();
  for( ; iter != histograms.end(); iter++ ) {
    delete (*iter);
  }
  histograms.clear();
  m_result.clear();
}

◆ TofCalibration() [2/2]

TofCalibration::TofCalibration ( const int npar )

◆ ~TofCalibration() [2/2]

TofCalibration::~TofCalibration ( )

Member Function Documentation

◆ calculate() [1/2]

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

virtual

Definition at line 26 of file TofCalibration.cxx.

                                                                        {
 
  std::cout << setiosflags(ios::left) << setw(10) << icounter << setw(8) << data->size() << setw(30) << name() << std::endl;
 
  if( data->size() <= 0 ) {
    for( int i=0; i<X.num_row(); i++ ) {
      X[i] = 1.0;
    }
  }
  else {
    std::vector<Record*>::iterator iter = data->begin();
    for( ; iter!=data->end(); iter++ ) {
      calculate_record( (*iter), icounter );
    }
 
    calculate_result( icounter );
  }
 
  return;
}

◆ calculate() [2/2]

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

virtual

◆ calculate_record() [1/2]

virtual void TofCalibration::calculate_record	(	const Record *	r,
		unsigned int	icounter
	)

protectedpure virtual

Referenced by calculate().

◆ calculate_record() [2/2]

virtual void TofCalibration::calculate_record	(	const Record *	r,
		unsigned int	icounter
	)

protectedpure virtual

◆ calculate_result() [1/2]

virtual void TofCalibration::calculate_result ( unsigned int icounter )

protectedpure virtual

Referenced by calculate().

◆ calculate_result() [2/2]

virtual void TofCalibration::calculate_result ( unsigned int icounter )

protectedpure virtual

◆ fillRoot() [1/2]

void TofCalibration::fillRoot	(	const char *	file,
		bool	isbarrel
	)

Definition at line 78 of file TofCalibration.cxx.

                                                               {
 
  unsigned int ncounter = 1;
  if( isbarrel ) {
    ncounter = NBarrel;
  }
  else {
    ncounter = NEndcap;
  }
 
  unsigned int nhist = histograms.size();
  if( nhist != (ncounter*nHistPerCounter + nHistogram) ) {
    std::cout<<" tofcalgsec::TofCalibration: the number of histograms is NOT integral number of times of counter number!"<<std::endl;
    exit(0);
  }
  unsigned int ngraph = graphs.size();
  if( ngraph != (ncounter*nGraphPerCounter + nGraph) ) {
    std::cout<<" tofcalgsec::TofCalibration: the number of histograms is NOT integral number of times of counter number!"<<std::endl;
    exit(0);
  }
 
  TFile f(file,"RECREATE");
 
  gStyle->SetOptFit(1111);
  gStyle->SetOptStat(1111);
  gStyle->SetLabelSize(0.03,"x");
  gStyle->SetLabelSize(0.03,"y");
 
  char dirname[256];
  if( nHistPerCounter>1 ) {
    std::vector<TH1F*>::iterator iter = histograms.begin();
    for( unsigned int i=0; i<ncounter; i++ ) {
      sprintf( dirname, "tofid%i", i );
      TDirectory* cdresult = f.mkdir( dirname );
      cdresult->cd();
      
      for( unsigned int i=0; i<nHistPerCounter; i++ ) {
    (*(iter+i))->Write();
      }
      iter = iter + nHistPerCounter;
    }
  }
  else {
    sprintf( dirname, "tofid" );
    TDirectory* cdresult = f.mkdir( dirname );
    cdresult->cd();
    
    std::vector<TH1F*>::iterator iter = histograms.begin();
    for( ; iter != histograms.end(); iter++ ) {
      (*iter)->Write();
    }
  }
  
  f.Close();
 
  return;
}

◆ fillRoot() [2/2]

void TofCalibration::fillRoot	(	const char *	file,
		bool	isbarrel
	)

◆ fillTxt() [1/2]

void TofCalibration::fillTxt ( const char * file )

Definition at line 48 of file TofCalibration.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( unsigned int i=0; i<(*it).num_row(); i++ ) {
      for( unsigned int i=0; i<Npar2; 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( unsigned int i=0; i<(*it).num_row(); i++ ) {
      for( unsigned int i=0; i<Npar2; i++ ) {
    cout << (*it)[i] << " ";
      }
      cout << std::endl;
    }
  }
 
  return;
}