diff --git a/data/rawDataDSR4/LOGBOOK.md b/data/rawDataDSR4/LOGBOOK.md
new file mode 100644
index 0000000000000000000000000000000000000000..d4ac709a96e1e45d025c5c47149333c657e7d51e
--- /dev/null
+++ b/data/rawDataDSR4/LOGBOOK.md
@@ -0,0 +1,21 @@
+Raw data. Both in binary and root format. 
+
+Logbook
+
+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|		| 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	      |
+------------------------|---------------|---------------------------------------------------------------------|
+
+Change of PMT amplification: increasing voltage up to 1100 mV, threshold 25 mV
+
+------------------------|---------------|---------------------------------------------------------------------|
+NeuRad_test_071216_8.dat| 10000		| source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back 	      |
+NeuRad_test_081216_1.dat|		| trigger ch2	      				?????		      |
+NeuRad_test_081216_2.dat| 10000		| collimated gamma source; two channels only: ch0 and ch1	      |
diff --git a/dataClasses/AEvent.cpp b/dataClasses/AEvent.cpp
index e9d13e504b9ea570b0ea5c5766cd396739c37f72..7d23ddd74bc3e018ec7dabd672d419d83a0fa5c3 100644
--- a/dataClasses/AEvent.cpp
+++ b/dataClasses/AEvent.cpp
@@ -60,18 +60,7 @@ void AEvent::ProcessEvent() {
 		fAmpPos[j] = fAmpPos[j] - fZeroLevel;
 	}
 
-	Double_t maxAmp = 0.;
-	Double_t maxAmpT = 0.;
-
-	maxAmp = fAmpPos[0];
-	for(Int_t j=0; j < fNPoints; j++) {
-		if(fAmpPos[j] > maxAmp) {
-			maxAmp = fAmpPos[j];
-			maxAmpT = fTime[j];
-		}
-	}
-	fAmpMax = maxAmp;
-	fTimeAmpMax = maxAmpT;
+	SetMaxAmplitudes();
 
 	SetGraphs();
 	FindFrontProperties();
@@ -267,4 +256,19 @@ Double_t AEvent::GetT_90() {
 
 Double_t AEvent::GetEdgeSlope() {
 		return EdgeSlope;
+
+void AEvent::SetMaxAmplitudes() {
+	Double_t maxAmp = 0.;
+	Double_t maxAmpT = 0.;
+
+	maxAmp = fAmpPos[0];
+	for(Int_t j=0; j < fNPoints; j++) {
+		if(fAmpPos[j] > maxAmp) {
+			maxAmp = fAmpPos[j];
+			maxAmpT = fTime[j];
+		}
+	}
+	fAmpMax = maxAmp;
+	fTimeAmpMax = maxAmpT;
+
 }
diff --git a/dataClasses/AEvent.h b/dataClasses/AEvent.h
index 172c31105b7cb7f66c92250e683f60ef204997d6..968e9e2455d577ff6fb41231d09990902b6ae7c9 100644
--- a/dataClasses/AEvent.h
+++ b/dataClasses/AEvent.h
@@ -111,6 +111,7 @@ public:
 
 private:
 	void Init();
+	void SetMaxAmplitudes();
 	void SetGraphs();
 	void SetCFD();
 };
diff --git a/macros/histCorr.C b/macros/histCorr.C
new file mode 100644
index 0000000000000000000000000000000000000000..ca55de8eefe5144524ae5c8f47f4019989dee804
--- /dev/null
+++ b/macros/histCorr.C
@@ -0,0 +1,56 @@
+//#include "TCanvas.h"
+
+void histCorr()
+{
+	gSystem->Load("../libData.so");
+	TFile *f = new TFile("../data/dataDSR4/analysis_07_8.root");
+	TTree *t = (TTree*)f->Get("atree");
+	TCanvas *c1 = new TCanvas("canvas", "Correlation pictures");
+	c1->Divide(2,2);
+
+//	TH2F *h1 = new TH2F("h1", "Time difference between two channels vs an integral of the first channel",500,-150.,150.,500,0.,5e-9);
+	TH2F *h1 = new TH2F("h1", "Time difference between two channels vs an integral of the first channel",125,-50.,50.,125,0.,5e-9);
+	TH2F *h2 = new TH2F("h2", "Time difference between two channels vs an integral of the second channel",125,-50.,50.,125,0e-9,5e-9);
+//	TH2F *h3 = new TH2F("h3", "Integral of the first channel vs integral of the second channel",500,0.,5e-9,500,0.,5e-9);
+	TH2F *h3 = new TH2F("h3", "Integral of the first channel vs integral of the second channel",500,0.,2e-9,500,0.,2e-9);
+	TH2F *h4 = new TH2F("h4", "Time difference between two channels vs amplitude of the first channel",100,10,1e-10,100,-300,300);
+
+	c1->cd(1);
+	t->Draw("Ach0.fChargeCFD : ( Ach0.fTimeCFD - Ach2.fTimeCFD ) >> h1","","col");
+	//x
+	h1->GetXaxis()->SetTitle("\\Delta \\tau_{0-2}  [ns]");
+	h1->GetXaxis()->CenterTitle();
+	//y
+	h1->GetYaxis()->SetTitle("Q_{0} [C]");
+	h1->GetYaxis()->CenterTitle();
+
+	c1->cd(2);
+	t->Draw("Ach2.fChargeCFD : ( Ach0.fTimeCFD - Ach2.fTimeCFD ) >> h2","","col");
+	h2->GetXaxis()->SetTitle("\\Delta \\tau_{0-2}  [ns]");
+	h2->GetXaxis()->CenterTitle();
+	//y
+	h2->GetYaxis()->SetTitle("Q_{2} [C]");
+	h2->GetYaxis()->CenterTitle();
+
+	c1->cd(3);
+	t->Draw("Ach0.fChargeCFD:Ach2.fChargeCFD >> h3","","col");
+	//x
+	h3->GetXaxis()->SetTitle("Q_{2} [C]");
+	h3->GetXaxis()->CenterTitle();
+	//y
+	h3->GetYaxis()->SetTitle("Q_{0} [C]");
+	h3->GetYaxis()->CenterTitle();
+
+	c1->cd(4);
+
+	//TCanvas *c2 = new TCanvas("Correlations 2");
+	//c2->Divide(1,2);
+
+}
+/*Time difference between two channels vs an integral of the first channel,
+2. Time difference between two channels vs an integral of the second channel,
+3. Integral of the first channel vs integral of the second channel,
+4. Time difference between two channels vs amplitude of the first channel,
+5. Time difference between two channels vs amplitude of the second channel,
+6. Amplitude of the first channel vs amplitude of the second channel
+*/