EmcBhaCalibData Class Reference

#include <EmcBhaCalibData.h>

Public Member Functions
	EmcBhaCalibData (int numxtals)
	EmcBhaCalibData (int numxtals, int MsgFlag)
	~EmcBhaCalibData ()
double &	matrixMEle (int row, int column)
EmcLSSMatrix *	getMatrixM ()
double &	vectorR (int ind)
double *	getVectorR ()
int	xtalInd (int ind)
int &	nXtals ()
int &	nXtalsHit ()
int &	xtalHits (int ind)
int &	xtalHitsDir (int ind)
bool	reduce ()
void	reset ()
void	writeOut (ostream &OutM, ostream &outV)
void	readIn (istream &InM, istream &InV)
void	printVec (int number)

Detailed Description

This class holds the information that is needed to do the calibration of the EMC with Bhabha events. It stores the accumulated matrix, the accumulated vector, the number of hits percrystal and other stuff.

The data are filled in the method EmcMCBhaEvent::CollectBhabha().

This software was developed for the BesIII collaboration. If you use all or part of it, please give an appropriate acknowledgement.

Copyright (C) 2005 IHEP

See also

Author

(Chunxiu Liu) (originator/contributor etc.);

Definition at line 63 of file EmcBhaCalibData.h.

Constructor & Destructor Documentation

EmcBhaCalibData() [1/2]

EmcBhaCalibData::EmcBhaCalibData

(

int

numxtals

)

Definition at line 48 of file EmcBhaCalibData.cxx.

  : _nXtals(numxtals),
    _nXtalsHit(0)
{
  //  _matrixM = new EmcLSSMatrix(_nXtals,600);
  _matrixM = new EmcLSSMatrix(_nXtals,_nXtals);
  _vectorR = new double[_nXtals];
  _xtalInd = new int[_nXtals];
  _xtalHits = new int[_nXtals];
  _xtalHitsMax = new int[_nXtals];
  for (int i=0;i<_nXtals;i++) {
    _vectorR[i] = 0.;
    _xtalInd[i] = i;  
    _xtalHits[i] = 0;
    _xtalHitsMax[i] = 0;
  }
  _MsgFlag =0;
}

EmcBhaCalibData() [2/2]

EmcBhaCalibData::EmcBhaCalibData	(	int	numxtals,
		int	MsgFlag
	)

Definition at line 70 of file EmcBhaCalibData.cxx.

  : _nXtals(numxtals),
    _nXtalsHit(0)
{
  //  _matrixM = new EmcLSSMatrix(_nXtals,600);
  _matrixM = new EmcLSSMatrix(_nXtals,_nXtals, MsgFlag);
  _vectorR = new double[_nXtals];
  _xtalInd = new int[_nXtals];
  _xtalHits = new int[_nXtals];
  _xtalHitsMax = new int[_nXtals];
  for (int i=0;i<_nXtals;i++) {
    _vectorR[i] = 0.;
    _xtalInd[i] = i;
    _xtalHits[i] = 0;
    _xtalHitsMax[i] = 0;
  }
 
  _MsgFlag = MsgFlag;
}

~EmcBhaCalibData()

EmcBhaCalibData::~EmcBhaCalibData

(

)

Definition at line 93 of file EmcBhaCalibData.cxx.

                                  {
  if ( 0 != _matrixM) {
    delete _matrixM; 
    _matrixM = 0;
  }
  if ( 0 != _vectorR) {
    delete[] _vectorR;
    _vectorR = 0;
  }
  if ( 0 != _xtalInd) {
    delete[] _xtalInd;
    _xtalInd = 0;
  }
  if ( 0 != _xtalHits) {
    delete[] _xtalHits;
    _xtalHits = 0;
  }
  if ( 0 != _xtalHitsMax) {
    delete[] _xtalHitsMax;
    _xtalHitsMax = 0;
  }
}

Member Function Documentation

getMatrixM()

EmcLSSMatrix * EmcBhaCalibData::getMatrixM ( )

inline

Definition at line 81 of file EmcBhaCalibData.h.

{ return _matrixM; }

getVectorR()

double * EmcBhaCalibData::getVectorR ( )

inline

Definition at line 87 of file EmcBhaCalibData.h.

{ return _vectorR; }

matrixMEle()

double & EmcBhaCalibData::matrixMEle ( int  row, int  column  )

inline

Definition at line 77 of file EmcBhaCalibData.h.

{ return ((*_matrixM)(row,column)); }

nXtals()

int & EmcBhaCalibData::nXtals ( )

inline

Definition at line 93 of file EmcBhaCalibData.h.

{ return _nXtals; }

Referenced by EmcBhaCalib::execute().

nXtalsHit()

int & EmcBhaCalibData::nXtalsHit ( )

inline

Definition at line 95 of file EmcBhaCalibData.h.

{ return _nXtalsHit; }

Referenced by EmcBhaCalib::execute().

printVec()

void EmcBhaCalibData::printVec

(

int

number

)

Definition at line 265 of file EmcBhaCalibData.cxx.

                                    {
 
  if (_MsgFlag <= 3) {
    std::cout << "EmcBhaCalibData:: INFO "<< "R ( " << number << " ): " << _vectorR[number] 
    << "  hits: " << _xtalHits[number] 
    << "  direct hits: " << _xtalHitsMax[number]
    << endl;
  }
}

Referenced by EmcBhaCalib::execute().

readIn()

void EmcBhaCalibData::readIn	(	istream &	InM,
		istream &	InV
	)

Definition at line 174 of file EmcBhaCalibData.cxx.

                                                   {
 
  _matrixM->readIn( InM );
  
  _matrixM->print_NonZeros();
 
  //cout<<"num_filled_rows(48)="<<_matrixM->num_filled_rows(48)<<endl;
  //cout<<"num_filled_cols(48)="<<_matrixM->num_filled_cols(48)<<endl;
 
  //cout<<"matrix(47,47)"<<_matrixM->operator()(47,47)<<endl;
 
  double vect_ele;
  int nr_ele;
  int xtal_ind;
  int hits;
  int hits_max;
 
  InV >> nr_ele;
 
  for (int i=0; i<nr_ele; i++) {
 
    InV >> xtal_ind >> hits >> hits_max
    >> vect_ele ;
 
    if (_vectorR[xtal_ind] == 0. ) {
      _nXtalsHit++;
    }
    
    _vectorR[xtal_ind] += vect_ele ;
    _xtalInd[xtal_ind] = xtal_ind;
    _xtalHits[xtal_ind] += hits;
    _xtalHitsMax[xtal_ind] += hits_max;
 
 }
 
}

reduce()

bool EmcBhaCalibData::reduce

(

)

Definition at line 212 of file EmcBhaCalibData.cxx.

                        {
 
  bool successful = true;
 
  //delete all zero elements in matrix and vector R
  //save only non zero elements (this format is needed by slap),
  //keep the original xtal index in array xtalInd
 
  //array to store the new row/column reference for later matrix reduction 
  int* xRef = new int[_nXtals];
  memset(xRef, -1, (_nXtals*sizeof(int)));
 
  _nXtalsHit = 0;
  for ( int _arrayIndx = 0; 
    _arrayIndx < _nXtals; _arrayIndx++) {
 
    if ( _vectorR[_arrayIndx] > 0. ) {
 
      _vectorR[_nXtalsHit] = _vectorR[_arrayIndx];
      _xtalInd[_nXtalsHit] = _xtalInd[_arrayIndx];
      _xtalHits[_nXtalsHit] = _xtalHits[_arrayIndx];
      _xtalHitsMax[_nXtalsHit] = _xtalHitsMax[_arrayIndx];
      xRef[_arrayIndx] = _nXtalsHit;
 
      _nXtalsHit++;
    }
  }
 
  successful = _matrixM->reduce_Matrix( xRef );
 
  delete [] xRef;
 
  return successful;
}

Referenced by EmcBhaCalib::execute().

reset()

void EmcBhaCalibData::reset

(

)

Definition at line 249 of file EmcBhaCalibData.cxx.

{
  _nXtalsHit = 0;
  _matrixM->reset();
 
  for (int i=0;i<_nXtals;i++) {
    _vectorR[i]= 0.;
    _xtalInd[i] = i;
    _xtalHits[i] = 0;
    _xtalHitsMax[i] = 0;
  }
 
}

vectorR()

double & EmcBhaCalibData::vectorR ( int  ind )

inline

Definition at line 84 of file EmcBhaCalibData.h.

{ return (_vectorR[ind]); }

writeOut()

void EmcBhaCalibData::writeOut	(	ostream &	OutM,
		ostream &	outV
	)

Definition at line 130 of file EmcBhaCalibData.cxx.

                                                       {
 
  //write only non zero elements to file
  OutM.setf(ios::scientific,ios::floatfield);
  OutV.setf(ios::scientific,ios::floatfield);
  OutV.precision(12);
  OutM.precision(12);
 
  _matrixM->writeOut( OutM );
 
  _nXtalsHit = 0;
 
  double* vect_ele = _vectorR;
  double* vect_end = _vectorR + _nXtals;
  int* xtal_ind = _xtalInd;
  int* hits_max = _xtalHitsMax;
  int* hits = _xtalHits;
 
  while ( vect_ele < vect_end ) {
    if ( *vect_ele != 0.  ){
      _nXtalsHit++;
    }
    vect_ele++;
  }
 
  OutV << _nXtalsHit << " ";
 
  vect_ele = _vectorR;
 
  while ( vect_ele < vect_end ) {
    if ( *vect_ele != 0. ) {
      OutV << *xtal_ind << " " << *hits << " " << *hits_max << " " 
       << *vect_ele << " ";
    }
    xtal_ind++;
    vect_ele++;
    hits_max++;
    hits++;
  }
 
}

Referenced by EmcSelBhaEvent::OutputMV().

xtalHits()

int & EmcBhaCalibData::xtalHits ( int  ind )

inline

Definition at line 97 of file EmcBhaCalibData.h.

{ return (_xtalHits[ind]); }

xtalHitsDir()

int & EmcBhaCalibData::xtalHitsDir ( int  ind )

inline

Definition at line 99 of file EmcBhaCalibData.h.

{ return (_xtalHitsMax[ind]); }

xtalInd()

int EmcBhaCalibData::xtalInd ( int  ind )

inline

Definition at line 90 of file EmcBhaCalibData.h.

{ return (_xtalInd[ind]); }

Referenced by EmcBhaCalib::execute().

---

The documentation for this class was generated from the following files:

• EmcBhaCalibData.h
• EmcBhaCalibData.cxx