calibration1.cxx 1.06 KB
Newer Older
Kostyleva D.A's avatar
Kostyleva D.A committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
void calibration1()
{
	gSystem->Load("/home/dariak/AculUtils/libAculData.so");
//	gSystem->Load("/home/vratik/workspace/AculUtils/libAculData.so");

	AculCalibration cal;
	cal.SetInputParameters("/home/dariak/AculUtils/parforcal.par"); //from .par
	cal.PrintInputParameters();
	cal.CalculateCalibParameters("/home/dariak/AculUtils/exp1016/clb_raw_files/sq22/clb07_0001.root", "SQ22", "AnalysisxTree", 100, 1500);

//	    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);
	}*/
}