#include <calib_endcap_atten.h>

Inheritance diagram for calib_endcap_atten:

Public Member Functions
	calib_endcap_atten (const unsigned int nrbin)

	~calib_endcap_atten ()

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

	calib_endcap_atten (const unsigned int nrbin)

	~calib_endcap_atten ()

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

Public Member Functions inherited from TofCalibFit
	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)

Additional Inherited Members
Protected Attributes inherited from TofCalibFit
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 10 of file InstallArea/include/tofcalgsec/tofcalgsec/calib_endcap_atten.h.

Constructor & Destructor Documentation

◆ calib_endcap_atten() [1/2]

calib_endcap_atten::calib_endcap_atten ( const unsigned int nrbin )

Definition at line 13 of file calib_endcap_atten.cxx.

                                                                :TofCalibFit( false, nEndcapAtten ) {
 
  nKind          = 1;    // pulse height
  nBinPerCounter = nrbin + 1;
 
  nHistPerCounter   = nKind*nBinPerCounter;
  nCanvasPerCounter = 2;
  CanvasPerCounterName.push_back( static_cast<string>("Pulse Height Most Probable Value") );  
  CanvasPerCounterName.push_back( static_cast<string>("Pulse Height Sigma") );  
  nGraphPerCanvasPerCounter.push_back(1);
  nGraphPerCanvasPerCounter.push_back(1);
 
  nHistogram      = 0;
  nCanvas         = 2;
  CanvasName.push_back( static_cast<string>("Pulse Height Most Probable Value vs TOF Counter Number") );  
  CanvasName.push_back( static_cast<string>("Pulse Height Sigma vs TOF Counter Number") );  
  nGraphPerCanvas.push_back(1);
  nGraphPerCanvas.push_back(1);
 
  int numGraphs1 = 0;
  std::vector<unsigned int>::iterator iter = nGraphPerCanvasPerCounter.begin();
  for( ; iter!=nGraphPerCanvasPerCounter.end(); iter++ ) {
    numGraphs1 = numGraphs1 + (*iter);
  }
  if( numGraphs1 != nGraphEcAtten ) {
    cout << "tofcalgsec::calib_endcap_atten: the number of Graphs is NOT reasonable!!!" << endl;
    exit(0);
  }
  int numGraphs2 = 0;
  iter = nGraphPerCanvas.begin();
  for( ; iter!=nGraphPerCanvas.end(); iter++ ) {
    numGraphs2 = numGraphs2 + (*iter);
  }
  if( numGraphs2 != nGraphEcAtten ) {
    cout << "tofcalgsec::calib_endcap_atten: the number of Graphs is NOT reasonable!!!" << endl;
    exit(0);
  }
 
  m_name = string("calib_endcap_atten");
 
  const int qbin      = 100;
  const double qbegin = 0.0;
  const double qend   = 5000.0;
 
  // histograms per counter
  char hname[256];
  for( unsigned int i=0; i<NEndcap; i++ ) {
    m_result.push_back( HepVector(nEndcapAtten,0) );
    for( unsigned int k=0; k<nrbin; k++ ) {
      sprintf( hname, "Q-id%i-r%i", i, k);
      m_histograms.push_back( new TH1F( hname, hname, qbin, qbegin, qend ) );
      m_fitresult.push_back( HepVector(nParEcAtten,0) );
    }
    sprintf( hname, "Q0-id%i", i );
    m_histograms.push_back( new TH1F( hname, hname, qbin, qbegin, qend ) );
    m_fitresult.push_back( HepVector(nParEcAtten,0) );
  }
 
  rpos.resize( nrbin );
  rposerr.resize( nrbin );
  rstep = ( rend - rbegin )/nrbin;
  for( unsigned int i=0; i<nrbin; i++ ) {
    rpos[i] = rbegin + ( i+0.5 )*rstep;
    rposerr[i] = 0.5*rstep;
  }
  itofid.resize( NEndcap );
  itofiderr.resize( NEndcap );
  itofidstep = 1.0;
  for( unsigned int i=0; i<NEndcap; i++ ) {
    itofid[i] = i*1.0;
    itofiderr[i] = 0.5;
  }
 
}

◆ ~calib_endcap_atten() [1/2]

calib_endcap_atten::~calib_endcap_atten ( )

Definition at line 89 of file calib_endcap_atten.cxx.

                                        {
  m_fitresult.clear();
  rpos.clear();
  rposerr.clear();
  itofid.clear();
  itofiderr.clear();
}

◆ calib_endcap_atten() [2/2]

calib_endcap_atten::calib_endcap_atten ( const unsigned int nrbin )

◆ ~calib_endcap_atten() [2/2]

calib_endcap_atten::~calib_endcap_atten ( )

Member Function Documentation

◆ calculate() [1/2]

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

virtual

Implements TofCalibFit.

Definition at line 98 of file calib_endcap_atten.cxx.

                                                                            {
 
  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 );
  }
  else {
    fillGraph( icounter );  // keep the m_graphs is not empty()
  }
 
  if( data->size() > 0 ) {
    std::vector<Record*>::iterator iter = data->begin();
    for( ; iter!=data->end(); iter++ ) {
      updateData( (*iter), icounter );
      fillRecordQ0( (*iter), icounter );
    }
    fitHistogramQ0( icounter );
  }
 
  if( icounter==(NEndcap-1) ) {
    fillGraphQ0();
  }
 
  return;
}

◆ calculate() [2/2]

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