diff --git a/macros/histCorr.C b/macros/histCorr.C
index 8c2ccde970c8488cea8ac134fe83248e4a755f7c..8f7e008c8c79c0a1139311a2cba276e980302ecd 100644
--- a/macros/histCorr.C
+++ b/macros/histCorr.C
@@ -1,22 +1,28 @@
 #include "TH2F.h"
 #include "TH3F.h"
 #include "TStyle.h"
+#include "TMath.h"
+#include "TF1.h"
 
 void histCorr(
 
 
-//		const char *filename = "Neurad_7_08_NEW.root",
-		const char *filename1 = "analysis_08_2.root",
-		const char *filename = "analysis_07_8.root",
-		const char *foldername = "report",
-//		const char *foldername = "8_2",
-//		const char *ext = ".eps",
-		const char *ext = ".gif",
+//	const char *filename = "Neurad_7_08_NEW.root",
+	const char *filename1 = "analysis_08_2.root",
+	const char *filename = "analysis_07_8.root",
+	const char *foldername = "report",
+//	const char *foldername = "8_2",
+//	const char *ext = ".eps",
+	const char *ext = ".gif",
+
+	const Bool_t can1 = 1, const Bool_t can2 = 0, // delta tLed and delta tcfd // dvumerki
+	const Bool_t can3 = 0, const Bool_t can4 = 0,
+	const Bool_t can5 = 0, const Bool_t can6 = 0) // tLEd and CFD for colollinated and NONcollimated data
 
-		const Bool_t can1 = 1, const Bool_t can2 = 1,
-		const Bool_t can3 = 1, const Bool_t can4 = 1,
-		const Bool_t can5 = 1, const Bool_t can6 = 1)
 {
+	Double_t rMin[4];
+	Double_t rMax[4];
+
 	using std::cout;
 	using std::endl;
 	gSystem->Load("../libData.so");
@@ -37,44 +43,66 @@ void histCorr(
 		c1->Divide(2,2);
 		c1->cd(1);
 
+		TF1 *g1 = new TF1("g1","gaus",0,20);
+		g1->SetParLimits(0,1e+2,1e+4);
+		g1->SetParLimits(1,-100,100);
+		g1->SetParLimits(2,0.1,10);
+
+
+		rMin[0]= -4.8; rMax[0] = 0.8;
+		rMin[1]= -7.6; rMax[1] = -0.8;
+		rMin[2]= 9.2; rMax[2] = 12.8;
+		rMin[3]= 6.8; rMax[3] = 12.4;
+
 		TH1F *h1 = new TH1F("h1", "delta tau_CFD for NONcollimated data",1000	, -200., 200);
 
-		//t->Draw("Ach0.fTimeCFD - Ach3.fTimeCFD >> h1","(Ach0.fTimeCFD - Ach3.fTimeCFD)>0 && (Ach0.fTimeCFD - Ach3.fTimeCFD)<8","");
-		t->Draw("Ach0.fTimeCFD - Ach3.fTimeCFD >> h1","","");
-		h1->GetXaxis()->SetRangeUser(-65, 8);
+		t->Draw("Ach0.fTimeCFD - Ach3.fTimeCFD >> h1","(Ach0.fTimeCFD - Ach3.fTimeCFD)>-20 && (Ach0.fTimeCFD - Ach3.fTimeCFD)<10","");
+		h1->GetXaxis()->SetRangeUser(-20, 10);
 		h1->GetXaxis()->SetTitle("delta tauCFD [ns]");
 		h1->GetXaxis()->CenterTitle();
+		g1->SetRange(rMin[0],rMax[0]);
+		h1->Fit("g1","R");
+		cout<<"sigma= "<<g1->GetParameter(2)<<endl;
 
 		c1->cd(2);
 		TH1F *h2 = new TH1F("h2", "delta tau_LED for NONcollimated data",1000, -200., 200);
 
-		//t->Draw("Ach0.fTimeLED - Ach3.fTimeLED >> h2","(Ach0.fTimeLED - Ach3.fTimeLED)>0 && (Ach0.fTimeLED - Ach3.fTimeLED)<8","");
-		t->Draw("Ach0.fTimeLED - Ach3.fTimeLED >> h2","","");
-		h2->GetXaxis()->SetRangeUser(-65, 8);
+		t->Draw("Ach0.fTimeLED - Ach3.fTimeLED >> h2","(Ach0.fTimeLED - Ach3.fTimeLED)>-20 && (Ach0.fTimeLED - Ach3.fTimeLED)<10","");
+		//t->Draw("Ach0.fTimeLED - Ach3.fTimeLED >> h2","","");
+		h2->GetXaxis()->SetRangeUser(-20, 10);
 		h2->GetXaxis()->SetTitle("delta tauLED [ns]");
 		h2->GetXaxis()->CenterTitle();
+		g1->SetRange(rMin[1],rMax[1]);
+		h2->Fit("g1","R");
+		cout<<"sigma= "<<g1->GetParameter(2)<<endl;
 
 		c1->cd(3);
 
-		TH1F *h3 = new TH1F("h3", "delta tau_CFD for collimated data",1000	, 0, 200);
+		TH1F *h3 = new TH1F("h3", "delta tau_CFD for collimated data",1000	, -200, 200);
 
-		t1->Draw("Ach0.fTimeCFD - Ach1.fTimeCFD >> h3","(Ach0.fTimeCFD - Ach1.fTimeCFD)>0 && (Ach0.fTimeCFD - Ach1.fTimeCFD)<150","");
+		t1->Draw("Ach0.fTimeCFD - Ach1.fTimeCFD >> h3","(Ach0.fTimeCFD - Ach1.fTimeCFD)>0 && (Ach0.fTimeCFD - Ach1.fTimeCFD)<20","");
 		//t1->Draw("Ach0.fTimeCFD - Ach1.fTimeCFD >> h3","","");
-		h3->GetXaxis()->SetRangeUser(0, 150);
+		h3->GetXaxis()->SetRangeUser(0, 20);
 		h3->GetXaxis()->SetTitle("delta tauCFD [ns]");
 		h3->GetXaxis()->CenterTitle();
+		g1->SetRange(rMin[2],rMax[2]);
+		h3->Fit("g1","R");
+		cout<<"sigma= "<<g1->GetParameter(2)<<endl;
 
 		c1->cd(4);
 		TH1F *h4 = new TH1F("h4", "delta tau_LED for collimated data",1000	, -200., 200);
 
-		t1->Draw("Ach0.fTimeLED - Ach1.fTimeLED >> h4","(Ach0.fTimeLED - Ach1.fTimeLED)>-2 && (Ach0.fTimeLED - Ach1.fTimeLED)<15","");
+		t1->Draw("Ach0.fTimeLED - Ach1.fTimeLED >> h4","(Ach0.fTimeLED - Ach1.fTimeLED)>0 && (Ach0.fTimeLED - Ach1.fTimeLED)<20","");
 		//t1->Draw("Ach0.fTimeLED - Ach1.fTimeLED >> h4","","");
-		h4->GetXaxis()->SetRangeUser(-2,15);
+		h4->GetXaxis()->SetRangeUser(0,20);
 		h4->GetXaxis()->SetTitle("delta tauLED [ns]");
 		h4->GetXaxis()->CenterTitle();
+		g1->SetRange(rMin[3],rMax[3]);
+		h4->Fit("g1","R");
+		cout<<"sigma= "<<g1->GetParameter(2)<<endl;
 
 		c1->Update();
-		//c1->Print(Form("../macros/picsDRS4/%s/deltaT%s", foldername, ext));
+		c1->Print(Form("../macros/picsDRS4/%s/deltaT%s", foldername, ext));
 	}
 
 	if (can2) {
@@ -208,45 +236,78 @@ void histCorr(
 		c5->Divide(2,2);
 		c5->cd(1);
 
+		rMin[0]= 138.4; rMax[0] = 140.4;
+		rMin[1]= 137; rMax[1] = 138.4;
+		rMin[2]= 136.4; rMax[2] = 142.4;
+		rMin[3]= 135.2; rMax[3] = 141.4;
+
+
+//		TF1 *gaus5 = new TF1("gaus5","TMath:gaus(0)");
+		TF1 *g5 = new TF1("g5","gaus",135,140);
+		g5->SetParLimits(0,1e+2,1e+4);
+		g5->SetParLimits(1,10,1000);
+		g5->SetParLimits(2,0.1,10);
+
 		TH1F *h12 = new TH1F("h12", "tau_CFD triggered channel for NONcollimated data",1000, 0., 200);
 
-		t->Draw("Ach0.fTimeCFD >> h12","Ach0.fTimeCFD>132 && Ach0.fTimeCFD<142","");
-		h12->GetXaxis()->SetRangeUser(120, 150);
+		t->Draw("Ach0.fTimeCFD >> h12","Ach0.fTimeCFD>132 && Ach0.fTimeCFD<150","");
+		h12->GetXaxis()->SetRangeUser(132, 150);
 		h12->GetXaxis()->SetTitle("tauCFD_{0} [ns]");
 		h12->GetXaxis()->CenterTitle();
+		g5->SetRange(rMin[0],rMax[0]);
+		h12->Fit("g5","R");
+		cout<<"sigma= "<<g5->GetParameter(2)<<endl;
 
 		c5->cd(2);
 		TH1F *h13 = new TH1F("h13", "tau LED triggered channel for NONcollimated data",1000, 0., 200);
 
-		t->Draw("Ach0.fTimeLED >> h13","Ach0.fTimeLED>132 && Ach0.fTimeLED<142","");
-		h13->GetXaxis()->SetRangeUser(120, 150);
+		t->Draw("Ach0.fTimeLED >> h13","Ach0.fTimeLED>132 && Ach0.fTimeLED<150","");
+		h13->GetXaxis()->SetRangeUser(132, 150);
 		h13->GetXaxis()->SetTitle("tauLED_{0} [ns]");
 		h13->GetXaxis()->CenterTitle();
+		g5->SetRange(rMin[1],rMax[1]);
+		h13->Fit("g5","R");
+		cout<<"sigma= "<<g5->GetParameter(2)<<endl;
 
 		c5->cd(3);
-
 		TH1F *h14 = new TH1F("h14", "tau_CFD nontriggered channel for NONcollimated data",1000	, 0., 200);
 
-		t->Draw("Ach1.fTimeCFD >> h14","Ach1.fTimeCFD > 120 && Ach1.fTimeCFD<180","");
-		h14->GetXaxis()->SetRangeUser(120, 150);
+		t->Draw("Ach1.fTimeCFD >> h14","Ach1.fTimeCFD > 132 && Ach1.fTimeCFD<150","");
+		h14->GetXaxis()->SetRangeUser(132, 150);
 		h14->GetXaxis()->SetTitle("tauCFD_{0} [ns]");
 		h14->GetXaxis()->CenterTitle();
+		g5->SetRange(rMin[2],rMax[2]);
+		h14->Fit("g5","R");
+		cout<<"sigma= "<<g5->GetParameter(2)<<endl;
 
 		c5->cd(4);
 		TH1F *h15 = new TH1F("h15", "tau LED nontriggered channel for NONcollimated data",1000	, 0., 200);
 
-		t->Draw("Ach1.fTimeLED >> h15","Ach1.fTimeLED>120 && Ach1.fTimeLED<180","");
-		h15->GetXaxis()->SetRangeUser(120, 	150);
+		t->Draw("Ach1.fTimeLED >> h15","Ach1.fTimeLED>132 && Ach1.fTimeLED<150","");
+		h15->GetXaxis()->SetRangeUser(132, 	150);
 		h15->GetXaxis()->SetTitle("tauLED_{0} [ns]");
 		h15->GetXaxis()->CenterTitle();
-
+		g5->SetRange(rMin[3],rMax[3]);
+		h15->Fit("g5","R");
+		cout<<"sigma= "<<g5->GetParameter(2)<<endl;
 		c5->Update();
+
 		c5->Print(Form("../macros/picsDRS4/%s/NONcollimatedT%s", foldername, ext));
 	}
 
 	if (can6) {
 		TCanvas *c6 = new TCanvas("c6", "Changing CFD parameters");
 		c6->Divide(2,2);
+		TF1 *g6 = new TF1("g6","gaus",135,140);
+		//g6->SetParLimits(0,1e+2,1e+4);
+		//g6->SetParLimits(1,10,1000);
+		//g6->SetParLimits(2,0.1,10);
+
+		rMin[0]= 138.8; rMax[0] = 140.2;
+		rMin[1]= 137.4; rMax[1] = 138.4;
+		rMin[2]= 125.6; rMax[2] = 131.6;
+		rMin[3]= 125.8; rMax[3] = 131.4;
+
 		c6->cd(1);
 
 		TH1F *h16 = new TH1F("h16", "tau_CFD triggered channel for collimated data",1000, 0., 200);
@@ -255,6 +316,10 @@ void histCorr(
 		h16->GetXaxis()->SetRangeUser(120, 150);
 		h16->GetXaxis()->SetTitle("tauCFD_{0} [ns]");
 		h16->GetXaxis()->CenterTitle();
+		g6->SetRange(rMin[0],rMax[0]);
+		h16->Fit("g6","R");
+		cout<<"sigma= "<<g6->GetParameter(2)<<endl;
+
 
 		c6->cd(2);
 		TH1F *h17 = new TH1F("h17", "tau LED triggered channel for collimated data",1000, 0., 200);
@@ -263,6 +328,9 @@ void histCorr(
 		h17->GetXaxis()->SetRangeUser(120, 150);
 		h17->GetXaxis()->SetTitle("tauLED_{0} [ns]");
 		h17->GetXaxis()->CenterTitle();
+		g6->SetRange(rMin[1],rMax[1]);
+		h17->Fit("g6","R");
+		cout<<"sigma= "<<g6->GetParameter(2)<<endl;
 
 		c6->cd(3);
 
@@ -272,6 +340,9 @@ void histCorr(
 		h18->GetXaxis()->SetRangeUser(120, 150);
 		h18->GetXaxis()->SetTitle("tauCFD_{0} [ns]");
 		h18->GetXaxis()->CenterTitle();
+		g6->SetRange(rMin[2],rMax[2]);
+		h18->Fit("g6","R");
+		cout<<"sigma= "<<g6->GetParameter(2)<<endl;
 
 		c6->cd(4);
 		TH1F *h19 = new TH1F("h19", "tau LED nontriggered channel for collimated data",1000	, 0., 200);
@@ -280,6 +351,9 @@ void histCorr(
 		h19->GetXaxis()->SetRangeUser(120, 	150);
 		h19->GetXaxis()->SetTitle("tauLED_{0} [ns]");
 		h19->GetXaxis()->CenterTitle();
+		g6->SetRange(rMin[3],rMax[3]);
+		h19->Fit("g6","R");
+		cout<<"sigma= "<<g6->GetParameter(2)<<endl;
 
 		c6->Update();
 		c6->Print(Form("../macros/picsDRS4/%s/collimatedT%s", foldername, ext));