#include <cstdio>
int FPedestal(char *title_base, fits &datafile, vector<int16_t> &data, vector<int16_t> &data_offset, vector<float> &drs_basemean, vector<float> &drs_gainmean, vector<float> &drs_triggeroffsetmean, UInt_t data_roi, size_t pixelnr, float &pedestal_mean, float &pedestal_rms)
{
//	TCanvas *canv_base = new TCanvas( "canv_base", "Determine the baseline", 100, 520, 700, 500 );
	TH1F *baseline = new TH1F("baseline", title_base, 1000,-10,40);
//	canv_base->SetLogy();
	baseline->GetXaxis()->SetTitle("Sample value (mV)");
	baseline->GetYaxis()->SetTitle("Entries");
	
//-------------------------------------------
//Find the baseline
//-------------------------------------------
	for (size_t i=0; i<datafile.GetNumRows(); i++)
	{
		datafile.GetRow(i);
		
		//Iterate over the slices
		for (UInt_t k=0; k<data_roi; k++)
		{
			UInt_t drs_calib_offset = (k+data_offset[pixelnr])%data_roi;
			float sample = (data[pixelnr*data_roi+k]*2000/4096.-drs_basemean[pixelnr*data_roi+drs_calib_offset]-drs_triggeroffsetmean[pixelnr*data_roi+k])/drs_gainmean[pixelnr*data_roi+drs_calib_offset]*1907.35;
			baseline->Fill(sample);
		}
	}
	pedestal_mean = baseline->GetXaxis()->GetBinCenter(baseline->GetMaximumBin());
	pedestal_rms = baseline->GetRMS();
	std::cout << "Value of maximal probability: " << pedestal_mean << " +- " << pedestal_rms << std::endl;
	
//	canv_base->cd();
//	baseline->Draw();
//	canv_base->Modified();
//	canv_base->Update();
	
	return 0;
}