#include <TSystem.h>
#include <iostream>

using namespace std;

void calibration1()
{
//	gSystem->Load("/home/dariak/AculUtils/libAculData.so");
	gSystem->Load("/home/vratik/workspace/AculUtils/libTELoss.so");
	gSystem->Load("/home/vratik/workspace/AculUtils/libAculData.so");

	AculCalibration cal;
//	cal.SetParFileName("/home/dariak/AculUtils/parforcal.par");	//from .par
	cal.SetParFileName("/home/vratik/workspace/AculUtils/parforcal.par");	//from .par
	cal.Init();
//return;
	cal.PrintInputParameters();

	TCanvas *c1 = new TCanvas();
	cal.ShowEnergySpectra("SQ22[]E.root", c1, 32);

//	return;
//	cal.CalculateCalibParameters("/home/dariak/AculUtils/exp1016/clb_raw_files/sq22/clb07_0001.root", "SQ22", "AnalysisxTree", 100, 1500);
//	cal.CalculateCalibParameters("/home/vratik/workspace/AculUtils/exp1016/clb_raw_files/clb07_0001.root", "SQ22", "AnalysisxTree", 100, 1500);
return;
//	    CalculateCalibParameters(const char* inputfile, const char* block, const char* treename, Int_t lowerchannel = 0, Int_t upperchannel = 4095, Int_t nEBins = 1000, Int_t lowersubaddress = 0, Int_t uppersubaddress = ADDRESSNUMBER-1);

	cal.PrintCalibrationParameters();
	TFile fr("/home/dariak/AculUtils/exp1016/SQ22[]E.root");
	cout << fr.GetListOfKeys()->GetEntries() << " histograms" << endl;
	TCanvas *c2 = new TCanvas();
	cal.FindEnergyPeaks(c2,"SQ22[]E.root","outputenergrun07.txt");

/*	TList *histList = fr.GetListOfKeys();
	TH1 *hWork = 0;
	for (Int_t i = 0; i < 32; i++) 
	{
		fr.GetObject(histList->At(i)->GetName(), hWork);		
		cal->PeaksFitting(hWork);
	}*/
}