#define nsensor 12
#include "Riostream.h"
#include "TMatrixD.h"
#include "TVectorD.h"
#include "TDecompLU.h"
#include "TDecompSVD.h"

class ManyMatrix{
public:
  ManyMatrix();
 
private:
  TMatrixD* RcurPtr[nsensor];
  TMatrixD* m_tmpPtr;
  TVectorD* r0curPtr[nsensor];
  TVectorD* v_tmpPtr;

};

// Constructor
ManyMatrix::ManyMatrix()
{

  for (Int_t i=0; i<nsensor; i++){
    RcurPtr[i] = new TMatrixD(3,3);
    m_tmpPtr=RcurPtr[i];
    m_tmpPtr->UnitMatrix();
    TMatrixD &m_tmp = *m_tmpPtr;
    m_tmp(0,0)=1; m_tmp(1,1)=1; m_tmp(2,2)=1;
  }

  for (Int_t i=0; i<nsensor; i++){
    r0curPtr[i] = new TVectorD(3); 
    v_tmpPtr=r0curPtr[i]; 
    TVectorD &v_tmp = *v_tmpPtr;
    v_tmp(2)=100.;
  }

  m_tmpPtr->Print();
  v_tmpPtr->Print();

}