/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Emc/EmcCalib/EmcBhaCalib/EmcBhaCalib-00-00-34/EmcBhaCalib/EmcBhaCalibData.h

Go to the documentation of this file.

//--------------------------------------------------------------------------
// Environment:
//      This software was developed for the BESIII collaboration.  If you
//      use all or part of it, please give an appropriate acknowledgement.
//
// Copyright Information:
//      Copyright (C) 2005           IHEP
//
//------------------------------------------------------------------------

#ifndef EMCBHACALIBDATA_H
#define EMCBHACALIBDATA_H

//-------------
// C Headers --
//-------------
extern "C" {
}

//---------------
// C++ Headers --
//---------------
#include <iostream>
#include <fstream>

//----------------------
// Base Class Headers --
//----------------------

//-------------------------------
// Collaborating Class Headers --
//-------------------------------
#include "EmcBhaCalib/EmcLSSMatrix.h"

using namespace std;

//------------------------------------
// Collaborating Class Declarations --
//------------------------------------
//              ---------------------
//              -- Class Interface --
//              ---------------------

class EmcBhaCalibData { 

public:

  // Constructor
  EmcBhaCalibData(int numxtals);

  EmcBhaCalibData(int numxtals, int MsgFlag);

  // Destructor
  ~EmcBhaCalibData( );

  // access
  // calibration matrix element
  double& matrixMEle(int row, int column) { return ((*_matrixM)(row,column)); }

  // access
  // pointer to calibration matrix
  EmcLSSMatrix* getMatrixM() { return _matrixM; }

  // calibration vector 
  double& vectorR(int ind) { return (_vectorR[ind]); }

  // pointer to calibration vector 
  double* getVectorR() { return _vectorR; }

  // array of xtal indices 
  int xtalInd(int ind) { return (_xtalInd[ind]); }

  // number of xtals 
  int& nXtals() { return _nXtals; }
  // number of xtals hit 
  int& nXtalsHit() { return _nXtalsHit; }
  // array with number of hits per xtal 
  int& xtalHits(int ind) { return (_xtalHits[ind]); }
  // array with number of direct hits per xtal 
  int& xtalHitsDir(int ind) { return (_xtalHitsMax[ind]); }

  //reduce matrix and vector to an array containing only non zero elements,
  //   needed for matrix inversion package slap 
  bool reduce();

  //reset all data 
  void reset();

  // write matrix and vector to files 
  void writeOut(ostream& OutM, ostream& outV);

  // read in matrix and vector from files 
  void readIn(istream& InM, istream& InV);

  // output vector element 
  void printVec(int number);

 
private:

  // Data members
  // calibration matrix 
  EmcLSSMatrix* _matrixM;
  // calibration vector 
  double* _vectorR;
  // array of xtal indices 
  int* _xtalInd;

  // number of xtals 
  int _nXtals;
  // number of xtals hit 
  int _nXtalsHit;
  // array with number of hits per xtal 
  int* _xtalHits;
  // array with number of direct hits per xtal 
  int* _xtalHitsMax;
  // message output flag, 
  int _MsgFlag;
};

#endif // EMCBHACALIBDATA_H

---