diff --git a/data/rawDataDSR4/LOGBOOK.md b/data/rawDataDSR4/LOGBOOK.md
index 0bfe215678b77d2967a68411c68f4177806f71db..4e38a9ff8cd3753a064e006512ba5b0d98899f99 100644
--- a/data/rawDataDSR4/LOGBOOK.md
+++ b/data/rawDataDSR4/LOGBOOK.md
@@ -2,20 +2,22 @@ Raw data. Both in binary and root format.
 
 Logbook
 
+PMT amplification: HV = -900 V, threshold ~ -10 mV
+
 File name		| Num of events | Comment
-------------------------|---------------|---------------------------------------------------------------------|
-NeuRad_test_071216_1.dat| 500           | collimated gamma source; ch0 and ch1 have coincidence; trigger ch0  |
-NeuRad_test_071216_2.dat| 5000          | collimated gamma source; ch0 and ch1 have coincidence; trigger ch0  |
-NeuRad_test_071216_3.dat| 1039		| neutron source Cf						      |
-NeuRad_test_071216_4.dat| 2000		| source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back	      |
-NeuRad_test_071216_5.dat| 2000		| source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back	      |
-NeuRad_test_071216_6.dat| 2000		| source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back	      |
-NeuRad_test_071216_7.dat| 2000		| source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back	      |
-------------------------|---------------|---------------------------------------------------------------------|
+------------------------|---------------|--------------------------------------------------------------------------------------------------------------------|
+NeuRad_test_071216_1.dat| 500           | collimated gamma source; ch0 and ch1 have coincidence; trigger ch0; 2 PMT H9500; periphery PMT pixels              |
+NeuRad_test_071216_2.dat| 5000          | collimated gamma source; ch0 and ch1 have coincidence; trigger ch0; 2 PMT H9500; periphery PMT pixels              |
+NeuRad_test_071216_3.dat| 1039		| non-collimated neutron source Cf;  2 PMT H9500; periphery PMT pixels				                     |
+NeuRad_test_071216_4.dat| 2000		| source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back; internal PMT pixels    			             |
+NeuRad_test_071216_5.dat| 2000		| source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back; internal PMT pixels   			             |
+NeuRad_test_071216_6.dat| 2000		| source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back; internal PMT pixels			                     |
+NeuRad_test_071216_7.dat| 2000		| source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back; internal PMT pixels			                     |
+------------------------|---------------|--------------------------------------------------------------------------------------------------------------------|
 
-Change of PMT amplification: increasing voltage up to 1100 mV, threshold 25 mV
+Change of PMT amplification: increasing HV up to -1100 V, threshold -25 mV
 
-------------------------|---------------|---------------------------------------------------------------------|
-NeuRad_test_071216_8.dat| 10000		| source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back 	      |
-NeuRad_test_081216_1.dat| 2204		| trigger ch2	      				?????		      |
-NeuRad_test_081216_2.dat| 10000		| collimated gamma source; two channels only: ch0 and ch1	      |
+------------------------|---------------|--------------------------------------------------------------------------------------------------------------------|
+NeuRad_test_071216_8.dat| 10000		| source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back; internal PMT pixels	 !!!!                                |
+NeuRad_test_081216_1.dat| 2204		| trigger ch2; 	     internal PMT pixels 				?????		                             |
+NeuRad_test_081216_2.dat| 10000		| collimated gamma source; two channels only: ch0 and ch1; changed PMT H9500 -ch0, H10966 - ch1 ; internal PMT pixels|
diff --git a/macros/histCorrNew.C b/macros/histCorrNew.C
index a1ead1c14ec49de4d1c7a09fc2a47196dc573a8e..eea620929bdd129d704ea9e9e7ead6dfd81cd7b4 100644
--- a/macros/histCorrNew.C
+++ b/macros/histCorrNew.C
@@ -4,7 +4,7 @@
 
 void histCorrNew (
 
-		const char *filename = "analysis_08_2.root",
+		const char *filename = "analysis_07_8.root",
 //		const char *foldername = "report",
 		const char *foldername = "8_2",
 //		const char *ext = ".eps",
@@ -64,7 +64,7 @@ void histCorrNew (
 
 		c1->Update();
 
-		c1->Print(Form("../macros/picsDRS4/file%s/time_front_corr1%s", foldername, ext));
+	//	c1->Print(Form("../macros/picsDRS4/file%s/time_front_corr1%s", foldername, ext));
 	}
 
 	if(can2) {
@@ -103,7 +103,7 @@ void histCorrNew (
 		h6->GetYaxis()->SetTitle("Max_ampl_{0} [V]");
 		h6->GetYaxis()->CenterTitle();
 
-		c2->Print(Form("../macros/picsDRS4/file%s/time_front_corr2%s", foldername, ext));
+	//	c2->Print(Form("../macros/picsDRS4/file%s/time_front_corr2%s", foldername, ext));
 	}
 
 	if(can3) {
@@ -111,22 +111,66 @@ void histCorrNew (
 		c3->Divide(2,2);
 
 		c3->cd(1);
-		TH1F *h7 = new TH1F("h7", "tau TF", 100, -5 , 25);
-		t->Draw("Ach0.fTimeFront - Ach1.fTimeFront >> h7","Ach0.fTimeFront>100. && Ach0.fTimeFront < 150 && Ach1.fTimeFront>100 && Ach1.fTimeFront < 150 && Ach0.fEdgeSlope > 0.01 && Ach1.fEdgeSlope > 0.01","");
-		h7->Fit("gaus", "","",7,14);
+		TH1F *h10 = new TH1F("h10", "delta t (front edge fit)", 100, -5 , 25);
+		t->Draw("Ach0.fTimeFront - Ach1.fTimeFront >> h10","Ach0.fTimeFront>100. && Ach0.fTimeFront < 150 && Ach1.fTimeFront>100 && Ach1.fTimeFront < 150","");
+		h10->Fit("gaus", "","",-5,5);
 		gStyle->SetOptFit();
+		h10->GetXaxis()->SetTitle("delta t [ns]");
+		h10->GetXaxis()->CenterTitle();
 
 		c3->cd(2);
-		TH1F *h8 = new TH1F("h8", "tau TF", 100, -5 , 25);
-		t->Draw("Ach0.fTimeFront - Ach1.fTimeFront >> h8","Ach0.fTimeFront>100. && Ach0.fTimeFront < 150 && Ach1.fTimeFront>100 && Ach1.fTimeFront < 150 && Ach0.fChargeTF > 0.05e-9 && Ach1.fChargeTF > 0.025e-9","");
-		h8->Fit("gaus","","",7,14);
+		TH1F *h7 = new TH1F("h7", "delta t (front edge fit) with edge slope correction", 100, -5 , 25);
+	//	t->Draw("Ach0.fTimeFront - Ach1.fTimeFront >> h7",
+	//	"Ach0.fTimeFront>100. && Ach0.fTimeFront < 150 && Ach1.fTimeFront>100 && Ach1.fTimeFront < 150 && Ach0.fEdgeSlope > 0.01 && Ach1.fEdgeSlope > 0.01","");
+		t->Draw("Ach0.fTimeFront - Ach1.fTimeFront >> h7","Ach0.fTimeFront>100. && Ach0.fTimeFront < 150 && Ach1.fTimeFront>100 && Ach1.fTimeFront < 150 && Ach0.fEdgeSlope > 0.005 && Ach1.fEdgeSlope > 0.001","");
+		h7->Fit("gaus", "","",-5,5);
 		gStyle->SetOptFit();
+		h7->GetXaxis()->SetTitle("delta t [ns]");
+		h7->GetXaxis()->CenterTitle();
+
+		c3->cd(3);
+		TH1F *h8 = new TH1F("h8", "delta t (front edge fit) with charge correction", 100, -5 , 25);
+		t->Draw("Ach0.fTimeFront - Ach1.fTimeFront >> h8",
+		"Ach0.fTimeFront>100. && Ach0.fTimeFront < 150 && Ach1.fTimeFront>100 && Ach1.fTimeFront < 150 && Ach0.fChargeTF > 0.05e-9 && Ach1.fChargeTF > 0.025e-9","");
+		h8->Fit("gaus","","",-5,5);
+		gStyle->SetOptFit();
+		h8->GetXaxis()->SetTitle("delta t [ns]");
+		h8->GetXaxis()->CenterTitle();
+
+		c3->cd(4);
+		TH1F *h9 = new TH1F("h9", "delta t (front edge fit) with rise time correction", 100, -5 , 25);
+		t->Draw("Ach0.fTimeFront - Ach1.fTimeFront >> h9",
+		"Ach0.fTimeFront> 100 && Ach0.fTimeFront<150 && Ach1.fTimeFront> 100 && Ach1.fTimeFront<150 && (Ach0.fTime90 - Ach0.fTime10) > 1.5 && (Ach0.fTime90 - Ach0.fTime10)<3 && (Ach1.fTime90 - Ach1.fTime10) > 1.5 && (Ach1.fTime90 - Ach1.fTime10)<3","");
+		h9->Fit("gaus","","",-5,5);
+		gStyle->SetOptFit();
+		h9->GetXaxis()->SetTitle("delta t [ns]");
+		h9->GetXaxis()->CenterTitle();
 
-		c3->Print(Form("../macros/picsDRS4/file%s/time_front_corr3%s", foldername, ext));
-
-
+	//	c3->Print(Form("../macros/picsDRS4/file%s/time_front_corr3%s", foldername, ext));
 	}
 
+	if(can4) {
+		TCanvas *c4 = new TCanvas ("c4", "Correlation pictures for collimated data: time difference");
+		c4->Divide(2,2);
+
+		c4->cd(1);
+		TH1F *h11 = new TH1F("h11", "rise time", 100, 0 , 5);
+		t->Draw("Ach0.fTime90 - Ach0.fTime10 >> h11","","");
+		//gStyle->SetLineColor(kRed);
+
+		c4->cd(2);
+		TH1F *h12 = new TH1F("h12", "rise time", 100, 0 , 5);
+		t->Draw("Ach1.fTime90 - Ach1.fTime10 >> h12","","");
+	//	h10->Fit("gaus", "","",7,14);
+	//	gStyle->SetOptFit();
+	//	h10->GetXaxis()->SetTitle("delta t [ns]");
+	//	h10->GetXaxis()->CenterTitle();
+
+		c4->cd(3);
+		TH1F *h13 = new TH1F("h13", "rise time", 100, -10, 10);
+		t->Draw("( Ach0.fTime10 + (Ach0.fTime90 - Ach0.fTime10)*0.5 ) - ( Ach1.fTime10 + (Ach1.fTime90 - Ach1.fTime10)*0.5 )>> h13","Ach0.fEdgeSlope > 0.03 && Ach1.fEdgeSlope > 0.006","");
+		h13->Fit("gaus","","",-5,5);
+		gStyle->SetOptFit();
 
-
+	}
 }