#include <cstdio>
int FIntFixedPosAllPx(fits &datafile, vector<int16_t> &data, vector<int16_t> &data_offset, unsigned int &data_num, UInt_t data_px, vector<float> &drs_basemean, vector<float> &drs_gainmean, vector<float> &drs_triggeroffsetmean, UInt_t data_roi, UInt_t integration_size, UInt_t integration_delay, TH1F* spectrum[])
//Data, calibration, data_roi, data_num, data_px, threshold, two histograms
{
//	UInt_t integration_size = 5;
//	UInt_t integration_delay = 238;
	
//	float sample, last_sample, integral;
	float integral;
	UInt_t drs_calib_offset;

//**********************************************************************************	
	for (size_t i=0; (i<datafile.GetNumRows()) && (i<100000); i++)
	{
		datafile.GetRow(i);
		cout << "Event number: " << data_num << endl;
		
		//Iterate over the pixels
		for (int j=0; j<data_px; j++)
		{
//		size_t j=pixelnr; //Fix the Pixel to a SoftID
		
			//Integrate
			integral = 0;
			for(UInt_t l=integration_delay; l<integration_delay+integration_size; l++)
			{
				drs_calib_offset = (l+data_offset[j])%data_roi;
				integral+=(data[j*data_roi+l]*2000/4096.-drs_basemean[j*data_roi+drs_calib_offset]-drs_triggeroffsetmean[j*data_roi+l])/drs_gainmean[j*data_roi+drs_calib_offset]*1907.35;
			}
			spectrum[j]->Fill(integral);
//			for(Int_t l=-10; l<60; l++)
//			{
//				drs_calib_offset = (l+data_offset[j])%data_roi;
//				sample = (data[j*data_roi+l]*2000/4096.-drs_basemean[j*data_roi+drs_calib_offset]-drs_triggeroffsetmean[j*data_roi+l])/drs_gainmean[j*data_roi+drs_calib_offset]*1907.35;
//			}
//			k+=min_dist;
//			
//			drs_calib_offset = (start_sample-1+data_offset[j])%data_roi;
//			last_sample = (data[j*data_roi+start_sample-1]*2000/4096.-drs_basemean[j*data_roi+drs_calib_offset]-drs_triggeroffsetmean[j*data_roi+start_sample-1])/drs_gainmean[j*data_roi+drs_calib_offset]*1907.35;
//			for (UInt_t k=start_sample; k<end_sample; k++)
//			{
//				drs_calib_offset = (k+data_offset[j])%data_roi;
//				sample = (data[j*data_roi+k]*2000/4096.-drs_basemean[j*data_roi+drs_calib_offset]-drs_triggeroffsetmean[j*data_roi+k])/drs_gainmean[j*data_roi+drs_calib_offset]*1907.35;
//				if((last_sample<threshold)&&(sample>threshold))
//				{
//				}
//				//TBD: mistake that after the deadtime the last_sample must be set new...
//				//Solution: first process full pipeline
//				last_sample = sample;
//			}
		}
	}
	return 0;
}