#include <iostream>
#include <fstream>
#include <cmath>
#include <TApplication.h>
#include <TCanvas.h>
#include <TH1F.h>

#define NOFCHANNELS 5	

using namespace std; 


int spark() {

TCanvas * c1 = new TCanvas("c", "c", 1200, 1200);

ifstream ifs("trial_3.txt");  				// construct the interface to the file
string line;

int npoints=0;
	
	while(!ifs.eof()){					//count number of lines
		getline(ifs,line);
		npoints++;
	}

ifs.close();

TH1F* histo[NOFCHANNELS];	
char name[50];

	for(int i=0;i<NOFCHANNELS;i++){
		
		sprintf(name,"channel%d",i+1);					//fill the char with injj00_i	
		histo[i] = new TH1F(name,name,npoints,0,npoints);
	}

ifstream ifs("trial_3.txt");  

	while(!ifs.eof()) {		
       		
		int bin = 1;
		float entries[NOFCHANNELS]=0;
	
		line.clear();		
		
		getline(ifs, line);				//get and store the entire line of the file

        	const char* cline = line.c_str();		// convert string to char
		
		sscanf(cline, "%f %f %f %f %f", &entries[0], &entries[1] ,&entries[2],&entries[3],&entries[4]);	//read and store data from string
		
		histo[0]->SetBinContent(bin,entries[0]);		
		/*for(int i=0;i<NOFCHANNELS;i++){
		
			histo[i]->SetBinContent(bin,entries[i]);
		}*/	
		
		bin++;
	}

/*TSpectrum* peaks = new TSpectrum();

double sigma = 2.;
double treshold = 20;
int position[1];

peaks -> Search(histo, sigma, "", treshold);
peaks -> Print >> position;
histo-> GetBinValue
*/
histo[0]->Draw("C");
/*histo[1]->Draw("Csame");
histo[2]->Draw("Csame");
histo[3]->Draw("Csame");
histo[4]->Draw("Csame");
*/
c1->Draw();


}