EmcBhaCalib Class Reference

#include <EmcBhaCalib.h>

List of all members.

Public Member Functions
	EmcBhaCalib (const std::string &name, ISvcLocator *pSvcLocator)
	EmcBhaCalib (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	execute ()
StatusCode	execute ()
StatusCode	finalize ()
StatusCode	finalize ()
virtual void	help ()
virtual void	help ()
StatusCode	initialize ()
StatusCode	initialize ()
virtual	~EmcBhaCalib ()
virtual	~EmcBhaCalib ()
Private Member Functions
void	initCalibConst ()
void	initCalibConst ()
void	ntupleOut ()
void	ntupleOut ()
bool	prepareConstants ()
bool	prepareConstants ()
void	printInfo (std::string xtalHitsDirFile)
void	printInfo (std::string xtalHitsDirFile)
bool	readInFromDB ()
bool	readInFromDB ()
bool	readInFromFile ()
bool	readInFromFile ()
bool	solveEqua ()
bool	solveEqua ()
void	writeOut ()
void	writeOut ()
void	writeOutConst ()
void	writeOutConst ()
Private Attributes
NTuple::Item< float >	calibConst
NTuple::Item< float >	calibConst
double	DigiCalibConst
NTuple::Item< float >	err
NTuple::Item< float >	err
NTuple::Item< float >	gi0
NTuple::Item< float >	gi0
NTuple::Item< long >	ixtal
NTuple::Item< long >	ixtal
double *	m_absoluteConstants
double *	m_absoluteConstants
bool	m_askForMatrixInversion
double *	m_calibConst
double *	m_calibConst
double *	m_calibConstUnred
double *	m_calibConstUnred
bool	m_calibrationSuccessful
double	m_convCrit
int	m_dirHitsCut
IEmcCalibConstSvc *	m_emcCalibConstSvc
IEmcCalibConstSvc *	m_emcCalibConstSvc
std::string	m_equationSolver
std::string	m_fileDir
std::string	m_fileExt
bool	m_fitPolynom
int	m_MsgFlag
EmcBhaCalibData *	m_myCalibData
EmcBhaCalibData *	m_myCalibData
long int	m_nrNonZeros
int	m_nrXtalsEnoughHits
double *	m_oldConstants
double *	m_oldConstants
bool	m_readDataFromDB
std::string	m_runNumberFile
NTuple::Tuple *	m_tuple1
NTuple::Tuple *	m_tuple1
bool	m_writeToFile
NTuple::Item< long >	nhitxtal
NTuple::Item< long >	nhitxtal

---

Detailed Description

class EmcBhaCalib Algorithm - performs calibration of EMC Xtals with Bhabha events and Chi^2 fit, the resulting system of linear equations of the fit is solved with the SLAP Algebra package

The matrix and the vector can either be fetched from file(s). If they are read from files the file names are supposed to have the name CalibVector.dat and calibMatrix.dat with an extension at the beginning. Only these extensions (usually the runnumbers) need to be listed in the read-in-file that is specified with runNumberFile.

Author:

Chunxiu Liu (originator/contributor etc.);

---

Constructor & Destructor Documentation

EmcBhaCalib::EmcBhaCalib (  const std::string &  name, ISvcLocator *  pSvcLocator )

:Algorithm(name, pSvcLocator), m_dirHitsCut(200), m_convCrit(0.000001), m_askForMatrixInversion(true), m_fitPolynom(true), //no used now m_writeToFile(true), m_readDataFromDB(false), m_equationSolver("GMR"), m_fileExt(""), m_fileDir("/ihepbatch/besdata/public/liucx/Calib/"), m_nrNonZeros(0), m_nrXtalsEnoughHits(0), m_runNumberFile("runnumbers.dat"), m_MsgFlag(0) { // Declare the properties declareProperty("equationSolver", m_equationSolver); declareProperty("dirHitsCut", m_dirHitsCut); declareProperty("writeToFile", m_writeToFile); declareProperty("fileExt", m_fileExt); declareProperty("fileDir", m_fileDir); declareProperty("readDataFromDB", m_readDataFromDB); declareProperty("askForMatrixInversion", m_askForMatrixInversion); declareProperty("fitPolynom", m_fitPolynom); declareProperty("convCrit", m_convCrit); declareProperty("runNumberFile", m_runNumberFile); declareProperty("MsgFlag", m_MsgFlag); m_calibConst = new double[6240]; m_calibConstUnred = new double[6240]; m_absoluteConstants = new double[6240]; m_oldConstants = new double[6240]; //ntuple m_tuple1=0; }

EmcBhaCalib::~EmcBhaCalib (   )  [virtual]

{ if ( 0 != m_calibConst) { delete [] m_calibConst; m_calibConst = 0; } if ( 0 != m_calibConstUnred) { delete [] m_calibConstUnred; m_calibConstUnred = 0; } if ( 0 != m_absoluteConstants) { delete [] m_absoluteConstants; m_absoluteConstants = 0; } if ( 0 != m_oldConstants) { delete [] m_oldConstants; m_oldConstants = 0; } if ( 0 != m_myCalibData) { delete m_myCalibData; m_myCalibData = 0; } }

EmcBhaCalib::EmcBhaCalib (  const std::string &  name, ISvcLocator *  pSvcLocator )

virtual EmcBhaCalib::~EmcBhaCalib (   )  [virtual]

---

Member Function Documentation

StatusCode EmcBhaCalib::execute (   )

StatusCode EmcBhaCalib::execute (   )

{ MsgStream log(msgSvc(), name()); log << MSG::DEBUG << "in execute()" << endreq; //read in accumulated matrix/matrices if ( m_readDataFromDB ) { m_calibrationSuccessful = readInFromDB(); log << MSG::INFO << "read in accumulated matrix from DB"<<endreq; } else { m_calibrationSuccessful = readInFromFile(); log << MSG::INFO <<"read in accumulated matrix from file"<<endreq; } if ( m_calibrationSuccessful == true ) { //ask first if to do a matrix inversion if (m_askForMatrixInversion==true) { m_calibrationSuccessful = false; log << MSG::INFO << " Well, we have " << m_nrXtalsEnoughHits << " crystals whith more " << "than " << m_dirHitsCut << " direct hits. Do you want do " << "a matrix inversion ? [N]" << endreq; m_calibrationSuccessful = true; } if ( m_calibrationSuccessful == true ) { //write added matrix and vector to files if ( m_writeToFile == true) { writeOut(); } cout <<"writeout"<<endl; m_myCalibData->printVec(2); //reduce to continious array of only non zeros elements for SLAP if ( m_calibrationSuccessful = m_myCalibData->reduce() ) { cout<<"m_calibrationSuccessful="<<m_calibrationSuccessful<<endl; if ( m_myCalibData->nXtalsHit() != m_myCalibData->nXtals() ){ log << MSG::INFO << " Reduced number of Xtals for calibration: " << m_myCalibData->nXtalsHit() << endreq; } cout<<"m_myCalibData->nXtalsHit()="<<m_myCalibData->nXtalsHit() <<"m_myCalibData->nXtals()="<<m_myCalibData->nXtals()<<endl; m_myCalibData->printVec(2); //for ( int i = 0; i< m_myCalibData->nXtals(); i++ ) { // cout<<m_calibConst[i]<<endl; //} //solve matrix equation if ( m_calibrationSuccessful = solveEqua() ) { // for ( int i = 0; i< m_myCalibData->nXtals(); i++ ) { //cout<<m_calibConst[i]<<endl; // } for (int i=0; i<m_myCalibData->nXtalsHit(); i++){ //fill an array of constants for storage in database m_calibConstUnred[m_myCalibData->xtalInd(i)] = m_calibConst[i]; } //do validation, fit polynom if necessary, fill CalList prepareConstants(); //if wanted write constants to plain ASCII file if ( m_writeToFile==true){ writeOutConst(); } ntupleOut(); } } else { log << MSG::WARNING << " Reduction of the Emc Bhabha calibration matrix FAILED !" << endl << " Will NOT perform Emc Bhabha calibration !" << endreq; return( StatusCode::FAILURE); } } } else { log << MSG::WARNING << " ERROR in reading in Emc Bhabha calibration data !" << endl << " Will NOT perform Emc Bhabha calibration !" << endreq; return( StatusCode::FAILURE); } return StatusCode::SUCCESS; }

StatusCode EmcBhaCalib::finalize (   )

StatusCode EmcBhaCalib::finalize (   )

{ MsgStream log(msgSvc(), name()); log << MSG::INFO << "in endRun()" << endreq; return StatusCode::SUCCESS; }

virtual void EmcBhaCalib::help (   )  [virtual]

void EmcBhaCalib::help (   )  [virtual]

{ MsgStream log(msgSvc(), name()); log << MSG::INFO<< "Performs the Chi square fit of the accumulated " << endreq; }

void EmcBhaCalib::initCalibConst (   )  [private]

void EmcBhaCalib::initCalibConst (   )  [private]

{ MsgStream log(msgSvc(), name()); for ( int i = 0; i< m_myCalibData->nXtals(); i++ ) { m_calibConst[i] = 20.; } ifstream inConst("EmcCalib.Const"); if ( !inConst ) { log << MSG::VERBOSE << " No starting values for calibration constants found ! " << "(EmcCalib.Const)" << endl << "Set them to 1.0 ! " << endreq; } else { log << MSG::VERBOSE << " Read in starting values of calibration constants !" << endreq; int nConst; inConst >> nConst; log << MSG::VERBOSE << " Have starting values for " << nConst << " Xtals !" << endl << "Set the others to 1.0 ! " << endmsg; int numberXtal; for ( int i = 0; i< nConst; i++ ) { inConst >> numberXtal >> m_calibConst[numberXtal]; } } int nConstEmc; nConstEmc= m_emcCalibConstSvc -> getDigiCalibConstNo() ; log << MSG::VERBOSE << " Have starting values for " << nConstEmc << " Xtals !" << endl << "Set the others to 1.0 ! " << endmsg; int numberXtal; //cout<<"ind"<<" "<<"getThetaIndex(ind)"<<" " // <<"getPhiIndex(ind)"<<"getEmcID"<<endl; for ( int i = 0; i< nConstEmc; i++ ) { numberXtal=i; m_calibConst[numberXtal]= m_emcCalibConstSvc -> getDigiCalibConst(i); //cout<<"m_calibConst[numberXtal]"<<numberXtal //<<" "<<m_calibConst[numberXtal]<<endl; //cout<<numberXtal<<" " //<<m_emcCalibConstSvc->getThetaIndex(numberXtal)<<" " //<<m_emcCalibConstSvc->getPhiIndex(numberXtal) //<<" "<<m_emcCalibConstSvc->getEmcID(numberXtal)<<endl; } }

StatusCode EmcBhaCalib::initialize (   )

StatusCode EmcBhaCalib::initialize (   )

{ MsgStream log(msgSvc(), name()); log << MSG::INFO << "in initialize()" << endreq; m_myCalibData = new EmcBhaCalibData(6240, m_MsgFlag); m_calibrationSuccessful = false; StatusCode status1; NTuplePtr nt1(ntupleSvc(),"FILE302/n1"); if ( nt1 ) m_tuple1 = nt1; else { m_tuple1=ntupleSvc()->book("FILE302/n1",CLID_ColumnWiseTuple,"Calib"); if( m_tuple1 ) { status1 = m_tuple1->addItem ("ixtal",ixtal); status1 = m_tuple1->addItem ("gi0",gi0); status1 = m_tuple1->addItem ("calibConst",calibConst); status1 = m_tuple1->addItem ("err",err); status1 = m_tuple1->addItem ("nhitxtal",nhitxtal); } else { // did not manage to book the N tuple.... log << MSG::ERROR <<"Cannot book N-tuple:" << long(m_tuple1) << endmsg; return StatusCode::FAILURE; } } // use EmcCalibConstSvc StatusCode scCalib; scCalib = Gaudi::svcLocator() -> service("EmcCalibConstSvc", m_emcCalibConstSvc); if( scCalib != StatusCode::SUCCESS){ log << MSG::ERROR << "can not use EmcCalibConstSvc" << endreq; } else { std::cout << "Test EmcCalibConstSvc DigiCalibConst(0)= " << m_emcCalibConstSvc -> getDigiCalibConst(0) << std::endl; } //init starting values for calibration constants from file or set to 1 initCalibConst(); return StatusCode::SUCCESS; }

void EmcBhaCalib::ntupleOut (   )  [private]

void EmcBhaCalib::ntupleOut (   )  [private]

{ //readout inputConst from file, only for MC data std::string InconstFile = m_fileDir; InconstFile += m_fileExt; InconstFile += "CalibInput.dat"; double inputConst[6240]; ifstream inputConstfile(InconstFile.c_str()); int ind=0; while( !inputConstfile.eof() ) { inputConstfile>>inputConst[ind]; ind++; } inputConstfile.close(); for (int i=0; i < m_myCalibData->nXtalsHit(); i++) { int xtalIndex=m_myCalibData->xtalInd(i); gi0 = 1.0/inputConst[xtalIndex]; err = (m_calibConstUnred[xtalIndex]-gi0)/gi0; calibConst=m_calibConstUnred[xtalIndex]; ixtal=xtalIndex; nhitxtal=m_myCalibData->xtalHitsDir(xtalIndex); m_tuple1->write(); } }

bool EmcBhaCalib::prepareConstants (   )  [private]

bool EmcBhaCalib::prepareConstants (   )  [private]

{ bool successfull=false; //the old constant //float theCalConst = 1.0; int chanIndex; double digiCalibConst[6240]; for ( int i = 0; i< m_myCalibData->nXtals(); i++ ) { chanIndex=m_myCalibData->xtalInd(i); //xtalInd(i) array of xtal indices //the old constant m_oldConstants[chanIndex]=m_calibConst[chanIndex]; //---- get the new absolute constant ---- //set it to 0 if we have not enough hits -> we'll keep the old constant if ( m_myCalibData->xtalHitsDir(chanIndex) < m_dirHitsCut ) m_absoluteConstants[chanIndex] = m_oldConstants[chanIndex]; else m_absoluteConstants[chanIndex] = m_oldConstants[chanIndex] * m_calibConstUnred[chanIndex]; } double DigiConst[6240]; for(int ind=0; ind<m_myCalibData->nXtals(); ind++ ) { digiCalibConst[ind] =m_calibConstUnred[ind] ; DigiCalibConst=digiCalibConst[ind]; DigiConst[ind]=m_absoluteConstants[ind]; } //define tree fill the absolute digicalibconsts into the root-file TTree* constr= new TTree("DigiCalibConst","DigiCalibConst"); constr->Branch("DigiCalibConst",DigiConst,"DigiConst[6240]/D"); constr->Fill(); double aa[6240]; constr->SetBranchAddress("DigiCalibConst", aa); constr->GetEntry(0); TFile fconst("EmcCalibConst.root", "recreate"); constr->Write(); fconst.Close(); delete constr; // end tree successfull=true; return successfull; }

void EmcBhaCalib::printInfo (  std::string  xtalHitsDirFile  )  [private]

void EmcBhaCalib::printInfo (  std::string  xtalHitsDirFile  )  [private]

{ ofstream xtalHitsDirOut; std::string xtalHitsDirFile = m_fileDir; xtalHitsDirFile += m_fileExt; xtalHitsDirFile += fileName; xtalHitsDirOut.open(xtalHitsDirFile.c_str()); //nXtalsHit() is number of xtals hit for (int i=0; i < m_myCalibData->nXtalsHit(); i++) { if ( m_myCalibData->xtalHitsDir(i) > m_dirHitsCut ) { int chanindex=m_myCalibData->xtalInd(i); xtalHitsDirOut<<chanindex<<endl; } } xtalHitsDirOut.close(); }

bool EmcBhaCalib::readInFromDB (   )  [private]

bool EmcBhaCalib::readInFromDB (   )  [private]

{ bool success = true; return success; }

bool EmcBhaCalib::readInFromFile (   )  [private]

bool EmcBhaCalib::readInFromFile (   )  [private]

{ MsgStream log(msgSvc(), name()); log << MSG::INFO << " Read Bhabha calibration data from files !" << endreq; std::string runNumberFile = m_fileDir; runNumberFile += m_fileExt; runNumberFile += m_runNumberFile; bool success = false; log << MSG::INFO << " Get file names from input file : " << runNumberFile << endreq; std::string vectorFile; std::string matrixFile; log << MSG::INFO << "Runnumber non-zeros xtals" << endreq; ifstream datafile(runNumberFile.c_str()); //cout<<"datafile="<<runNumberFile.c_str()<<""<<datafile.good()<<endl; if (datafile.good() > 0 ) { while( !datafile.eof() ) { ifstream vectorIn; ifstream matrixIn; std::string num; datafile >> num; if ( num.length() > 0 ) { vectorFile = m_fileDir; vectorFile += m_fileExt; vectorFile += num; vectorFile += "CalibVector.dat"; matrixFile = m_fileDir; matrixFile += m_fileExt; matrixFile += num; matrixFile += "CalibMatrix.dat"; vectorIn.open(vectorFile.c_str()); matrixIn.open(matrixFile.c_str()); if ( vectorIn.good() > 0 && matrixIn.good() > 0 ) { m_myCalibData->readIn(matrixIn,vectorIn); log << MSG::INFO << num << " " << m_myCalibData->getMatrixM()->num_nonZeros() << " " << m_myCalibData->nXtalsHit() << endreq; success = true; } else { if ( !vectorIn ) log << MSG::ERROR << " ERROR: Vector input file " << vectorFile << " doesn't exist !" << endreq; if ( !matrixIn ) log << MSG::ERROR << " ERROR: Matrix input file " << matrixFile << " doesn't exist !" << endreq; } vectorIn.close(); matrixIn.close(); } } } if ( success == true) { for (int i=0; i < m_myCalibData->nXtals(); i++) { if ( m_myCalibData->xtalHitsDir(i) > m_dirHitsCut ) { m_nrXtalsEnoughHits++; } } m_nrNonZeros = m_myCalibData->getMatrixM()->num_nonZeros(); log << MSG::INFO<< " Final matrix : " << "Number of non zero elements: " << m_nrNonZeros << " " << m_myCalibData->nXtalsHit() << endreq; } return success; }

bool EmcBhaCalib::solveEqua (   )  [private]

bool EmcBhaCalib::solveEqua (   )  [private]

{ MsgStream log(msgSvc(), name()); //----------------------------------------------------- // solve system of equations with SLAP package //----------------------------------------------------- log << MSG::VERBOSE << " Solve Matrix equation with method " << m_equationSolver << endl << "Xtals hit: " << m_myCalibData->nXtalsHit() << endl << "Nr of non Zeros: " << m_nrNonZeros << endreq; //input parameters for SLAP long int isym = 1; //matrix M is symmetric long int nsave = 20; long int itol = 0; //type of convergence criterion double convCriterion = m_convCrit; //convergence crtiterion long int maxIter = 1000; //maximum number of iterations long int errUnit = 6; //unit on which to write error estimation // at each iteration //working arrays needed by SLAP long int lengthDoubleWork; double* doubleWork; long int lengthIntWork; long int* intWork; //output parameters of slap long int iters=0; //number of iterations required double errorv=1000; //error estimate long int errorFlag; //sparse Ax=b solver. //uses the generalized minimum residual method //The preconditioner is diagonal scaling. if ( m_equationSolver == "GMR" ) { cout<<m_equationSolver<<endl; //workspaces lengthDoubleWork = (1 + m_myCalibData->nXtals()*(nsave+7) + nsave*(nsave+3)) + 50000; doubleWork = new double[lengthDoubleWork]; lengthIntWork = 50; intWork = new long int [lengthIntWork]; dsdgmr_ ( &(m_myCalibData->nXtalsHit()), m_myCalibData->getVectorR(), m_calibConst, &m_nrNonZeros, m_myCalibData->getMatrixM()->row(), m_myCalibData->getMatrixM()->column(), m_myCalibData->getMatrixM()->matrix(), &isym, &nsave, &itol, &convCriterion, &maxIter, &iters, &errorv, &errorFlag, &errUnit, doubleWork, &lengthDoubleWork, intWork, &lengthIntWork ); log << MSG::VERBOSE << " Error flag: " << errorFlag << endreq; if ( errorFlag < 0 ) errorFlag = labs(errorFlag) + 2; switch ( errorFlag ) { case 0: log << MSG::VERBOSE << " all went well" << endreq; break; case 1: log << MSG::ERROR << " insufficient storage allocated for _doubleWork or _intWork" << endreq; break; case 2: log << MSG::ERROR << " method failed to reduce norm of current residual" << endreq; break; case 3: log << MSG::ERROR << " insufficient length for _doubleWork" << endreq; break; case 4: log << MSG::ERROR << " inconsistent _itol and values" << endreq; break; default: log << MSG::ERROR << " unknown flag" << endreq; } log << MSG::VERBOSE << " Integer workspace used = " << intWork[8] << endl << " Double workspace used = " << intWork[9] << endreq; } //Routine to solve a symmetric positive definite linear //system Ax = b using the Preconditioned Conjugate //Gradient method. The preconditioner is diagonal scaling. if ( m_equationSolver == "PCG" ) { cout<<m_equationSolver<<endl; itol = 1; //workspaces lengthDoubleWork = 5 *m_myCalibData->nXtals() + 50000; doubleWork = new double[lengthDoubleWork]; lengthIntWork = 50; intWork = new long int [lengthIntWork]; dsdcg_( &(m_myCalibData->nXtalsHit()), m_myCalibData->getVectorR(), m_calibConst, &m_nrNonZeros, m_myCalibData->getMatrixM()->row(), m_myCalibData->getMatrixM()->column(), m_myCalibData->getMatrixM()->matrix(), &isym, &itol, &convCriterion, &maxIter, &iters, &errorv, &errorFlag, &errUnit, doubleWork, &lengthDoubleWork, intWork, &lengthIntWork ); switch ( errorFlag ) { case 0: log << MSG::VERBOSE << "all went well" << endreq; break; case 1: log << MSG::ERROR << " insufficient storage allocated for WORK or IWORK" << endreq; break; case 2: log << MSG::ERROR << " method failed to to converge in maxIter steps." << endreq; break; case 3:log << MSG::ERROR << " Error in user input. Check input value of _nXtal,_itol." << endreq; break; case 4:log << MSG::ERROR << " User error tolerance set too tight. " << "Reset to 500.0*D1MACH(3). Iteration proceeded." << endreq; break; case 5:log << MSG::ERROR << " Preconditioning matrix, M, is not Positive Definite. " << endreq; break; case 6: log << MSG::ERROR << " Matrix A is not Positive Definite." << endreq; break; default: log << MSG::ERROR << " unknown flag" << endreq; } log << MSG::VERBOSE << " Integer workspace used = " << intWork[9] << endl << "Double workspace used = " << intWork[10] << endreq; } log << MSG::VERBOSE << "------ Calibration fit statistics ------- " << endl << "maximum number of iterations =" << maxIter << endl << " number of iterations =" << iters << endl << "error estimate of error in final solution =" << errorv << endreq; if ( 0 != doubleWork) delete [] doubleWork; if ( 0 != intWork) delete [] intWork; if ( errorFlag != 0 ) return false; else return true; return true; }

void EmcBhaCalib::writeOut (   )  [private]

void EmcBhaCalib::writeOut (   )  [private]

{ ofstream vectorOut; std::string vectorFile = m_fileDir; vectorFile += m_fileExt; vectorFile += "allCalibVector.dat"; vectorOut.open(vectorFile.c_str()); ofstream matrixOut; std::string matrixFile = m_fileDir; matrixFile += m_fileExt; matrixFile += "allCalibMatrix.dat"; matrixOut.open(matrixFile.c_str()); MsgStream log(msgSvc(), name()); log << MSG::VERBOSE << " Write out files " << vectorFile << " " << matrixFile << endreq; m_myCalibData->writeOut(matrixOut,vectorOut); vectorOut.close(); matrixOut.close(); }

void EmcBhaCalib::writeOutConst (   )  [private]

void EmcBhaCalib::writeOutConst (   )  [private]

{ ofstream constOut; std::string constFile = m_fileDir; constFile += m_fileExt; constFile += "CalibConst.dat"; constOut.open(constFile.c_str()); constOut << "#crystalIndex relative-constant absolute-constant" << endl; //output constants to file for ( int i=0; i < m_myCalibData->nXtals(); i++) { long Xtal_Index = i; if(m_calibConstUnred[Xtal_Index]>0){ constOut << Xtal_Index << " " << m_calibConstUnred[Xtal_Index] << " " << m_absoluteConstants[Xtal_Index] << endl; } } constOut.close(); }

---

Member Data Documentation

NTuple::Item<float> EmcBhaCalib::calibConst [private]

NTuple::Item<float> EmcBhaCalib::calibConst [private]

double EmcBhaCalib::DigiCalibConst [private]

NTuple::Item<float> EmcBhaCalib::err [private]

NTuple::Item<float> EmcBhaCalib::err [private]

NTuple::Item<float> EmcBhaCalib::gi0 [private]

NTuple::Item<float> EmcBhaCalib::gi0 [private]

NTuple::Item<long> EmcBhaCalib::ixtal [private]

NTuple::Item<long> EmcBhaCalib::ixtal [private]

double* EmcBhaCalib::m_absoluteConstants [private]

double* EmcBhaCalib::m_absoluteConstants [private]

bool EmcBhaCalib::m_askForMatrixInversion [private]

double* EmcBhaCalib::m_calibConst [private]

double* EmcBhaCalib::m_calibConst [private]

double* EmcBhaCalib::m_calibConstUnred [private]

double* EmcBhaCalib::m_calibConstUnred [private]

bool EmcBhaCalib::m_calibrationSuccessful [private]

double EmcBhaCalib::m_convCrit [private]

int EmcBhaCalib::m_dirHitsCut [private]

IEmcCalibConstSvc* EmcBhaCalib::m_emcCalibConstSvc [private]

IEmcCalibConstSvc* EmcBhaCalib::m_emcCalibConstSvc [private]

std::string EmcBhaCalib::m_equationSolver [private]

std::string EmcBhaCalib::m_fileDir [private]

std::string EmcBhaCalib::m_fileExt [private]

bool EmcBhaCalib::m_fitPolynom [private]

int EmcBhaCalib::m_MsgFlag [private]

EmcBhaCalibData* EmcBhaCalib::m_myCalibData [private]

EmcBhaCalibData* EmcBhaCalib::m_myCalibData [private]

long int EmcBhaCalib::m_nrNonZeros [private]

int EmcBhaCalib::m_nrXtalsEnoughHits [private]

double* EmcBhaCalib::m_oldConstants [private]

double* EmcBhaCalib::m_oldConstants [private]

bool EmcBhaCalib::m_readDataFromDB [private]

std::string EmcBhaCalib::m_runNumberFile [private]

NTuple::Tuple* EmcBhaCalib::m_tuple1 [private]

NTuple::Tuple* EmcBhaCalib::m_tuple1 [private]

bool EmcBhaCalib::m_writeToFile [private]

NTuple::Item<long> EmcBhaCalib::nhitxtal [private]

NTuple::Item<long> EmcBhaCalib::nhitxtal [private]

---

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

• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/Emc/EmcCalib/EmcBhaCalib/EmcBhaCalib-00-00-18/EmcBhaCalib/EmcBhaCalib.h
• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/InstallArea/include/EmcBhaCalib/EmcBhaCalib/EmcBhaCalib.h
• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/Emc/EmcCalib/EmcBhaCalib/EmcBhaCalib-00-00-18/src/EmcBhaCalib.cxx

---