{
	gSystem->Load("CsICalib_C.so");

//	CsICalib cal;
	CsICalib cal("parameters/SQ23Alpha.par");


//	cal.PrintCuts();

//return;

	TCanvas *c1 = new TCanvas("c1", "Plain");	
	TCanvas *c2 = new TCanvas("c2", "Plain");
	TCanvas *c3 = new TCanvas("c3", "Plain");
	TCanvas *c4 = new TCanvas("c4", "Plain");

	cal.DrawBeam(c4, 4, "SQ23");

return;

//	cal.ReadParFile("SQ23Alpha.par");
//	cal.OpenTrees();

//	cal.ReadClbParameters("calSQ23Alpha.clb");
//	cal.DrawEnergyDeposite("SQ23", c4, 3, "same");

//return;

//	cal.DrawVariable("SQ23", 1, c1);
//	cal.DrawVariable("TDC", 2, c2, 0, 700);
//	cal.DrawVariable("QDC", 3, c3);
//return;

//	cal.DrawVariableCut("SQ23", 0, c1, "cutAlpha16", "cutAlpha16Amp");
//	cal.DrawVariableCut("SQ23", 0, c2, "cutAlpha20", "cutAlpha20Amp");
//	cal.DrawVariableCut("SQ23", 0, c3, "cutAlpha25", "cutAlpha25Amp");
//	cal.DrawVariableCut("SQ23", 0, c4, "cutAlpha30", "cutAlpha30Amp");

//	cal.GetPeakMean("SQ23", 0, c1, "cutAlpha16Amp", minAlpha16, maxAlpha16);
//	cal.GetPeakMean("SQ23", 1, c2, "cutAlpha20Amp", minAlpha20, maxAlpha20);
//	cal.GetPeakMean("SQ23", 2, c3, "cutAlpha25Amp", minAlpha25, maxAlpha25);
//	cal.GetPeakMean("SQ23", 3, c4, "cutAlpha30Amp", minAlpha30, maxAlpha30);

	cal.GetPeakMean("SQ23", 0, 0, c1, "cutAlpha16Amp");
	cal.GetPeakMean("SQ23", 1, 1, c2, "cutAlpha20Amp");
	cal.GetPeakMean("SQ23", 2, 2, c3, "cutAlpha25Amp");
	cal.GetPeakMean("SQ23", 3, 3, c4, "cutAlpha30Amp");

	TCanvas *cCal = new TCanvas("cCal", "calibration Alpha");
	cal.Calibrate(cCal);
//	cal.SaveClbGraphs("gSQ23Alpha.root", "RECREATE");
//	cal.WriteClbParameters("calSQ23Alpha.clb");

return;
}