{
gROOT->Reset();
#include "Riostream.h"
#include "TFile.h"
#include "TH1.h"
#include "TH2.h"
#include "TNtuple.h"

Float_t CalibNTuple[10];
mFile          = new TFile("pedestal.root","RECREATE");
string varlist = "iwaf:wafer:ladder:rmsp:rmsn:strip:longstrip";
mNtuple         = new TNtuple("CalibrationNTuple","Calibration Ntuple",varlist.c_str());
for (int p=0;p<10;p++)
    {
      CalibNTuple[p]  = 0;
    }

gSystem->Load("St_base"); // needed for StDbModifier 
gSystem->Load("StUtilities");

// DB-specific libs
gSystem->Load("libStDb_Tables.so");
gSystem->Load("StDbLib");

TString path    = "/star/data06/SSD/pedestal_calibration/run7/";

TFile curr_root_file("/star/data06/SSD/pedestal_calibration/run7/ssdNoise.20070608.213326.root");
cout << " file found..."<<endl;
St_ssdNoise *aaa = curr_root_file.Get("ssdNoise"); //  TKey *eventKey = root_file.GetKey("ssdStripCalib");
ssdNoise_st *bbb = aaa->GetTable();

cout << " ReadStrip : printing few pedestal/noise values "<< endl;
Int_t  mNWaferPerLadder = 16;
Int_t  nMaxStrip        = 768;
Int_t iWaf              = 0;
Int_t iLad              = 0;

for (Int_t i = 0 ; i < aaa->GetNRows(); i++)
{
  iWaf = bbb[i].id-(bbb[i].id/mNWaferPerLadder)*mNWaferPerLadder;
  iLad = bbb[i].id/16;
  Int_t idWaf = bbb[i].id;
  CalibNTuple[0] = idWaf;
  CalibNTuple[1] = iWaf;
  CalibNTuple[2] = iLad;

  for(Int_t ii=0;ii<nMaxStrip;ii++)
    {
      CalibNTuple[3] = bbb[i].rmsp[ii];
      CalibNTuple[4] = bbb[i].rmsn[ii];
      CalibNTuple[5] = ii;
      CalibNTuple[6] = ii + (iWaf)*nMaxStrip;
      mNtuple->Fill(CalibNTuple);
    }
  if (idWaf==1)
    {
      for(Int_t nStrip=0;nStrip<10;nStrip++)
	{
	  cout<<"ReadStrip: iLad,idWaf,nStrip,rmsP,rmsN = "<<iLad+1
	      <<" "<<idWaf+1
	      <<" "<<nStrip
	      <<" "<<(int)(bbb[i].rmsp[nStrip])
	      <<" "<<(int)(bbb[i].rmsn[nStrip])<<endl;
	}
    }
}  
cout<<"ReadStrip: Number of rows in the table : "<< aaa->GetNRows()<<endl;
curr_root_file->Close();
delete aaa;


mFile->Write();
mFile->Close();

}