From 7a0fcc974c352f0cbbe36688dc8c2fbe0d004875 Mon Sep 17 00:00:00 2001
From: Vratislav Chudoba <chudoba@jinr.ru>
Date: Wed, 2 Aug 2017 16:03:52 +0300
Subject: [PATCH] Figures of GS correlations implemented in figCorrGS.cxx

---
 macros/BeCorrPRC/figCorrGS.cxx | 125 +++++++++++++++++++++++++++++++++
 1 file changed, 125 insertions(+)
 create mode 100755 macros/BeCorrPRC/figCorrGS.cxx

diff --git a/macros/BeCorrPRC/figCorrGS.cxx b/macros/BeCorrPRC/figCorrGS.cxx
new file mode 100755
index 0000000..33546c3
--- /dev/null
+++ b/macros/BeCorrPRC/figCorrGS.cxx
@@ -0,0 +1,125 @@
+//#include "TFile.h"
+//#include "TCanvas.h"
+//#include "TH1F.h"
+
+void figCorrGS()
+{
+
+	TFile *fr1 = new TFile("figures/PRCepsilonT.root", "READ");
+	TFile *fr2 = new TFile("figures/PRCepsilonT90120.root", "READ");
+
+	TH1F *hExp[4];
+	TH1F *hSim[4];
+	TH1F *hInput[4];
+
+	TH1F *hCurr = 0;
+
+	hCurr = (TH1F*)fr1->Get("heeT0_0_4560");
+	hExp[0] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("heeT0_5_6075");
+	hExp[1] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("heeT0_5_7590");
+	hExp[2] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr2->Get("heeT0_0_90120");
+	hExp[3] = new TH1F(*hCurr);
+
+	hCurr = (TH1F*)fr1->Get("hseT0_0_4560");
+	hSim[0] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("hseT0_5_6075");
+	hSim[1] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("hseT0_5_7590");
+	hSim[2] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr2->Get("hseT0_0_90120");
+	hSim[3] = new TH1F(*hCurr);
+
+	hCurr = (TH1F*)fr1->Get("hieT0_0_4560");
+	hInput[0] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("hieT0_5_6075");
+	hInput[1] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("hieT0_5_7590");
+	hInput[2] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr2->Get("hieT0_0_90120");
+	hInput[3] = new TH1F(*hCurr);
+
+
+
+	TCanvas *c1 = new TCanvas();
+	c1->ToggleEditor();
+//	c1->SetWindowSize(800, 1200);
+	c1->SetWindowSize(800, 800);
+	c1->Divide(2, 2, .001, .001);
+
+
+	for (Int_t i = 0; i < 4; i++) {
+		c1->cd(i+1);
+
+//		cout << i << endl;
+
+		hSim[i]->SetTitle("");
+		hSim[i]->SetLineColor(kGray+1);
+		hSim[i]->SetFillColor(kGray+1);
+		hSim[i]->GetXaxis()->SetNdivisions(405);
+		hSim[i]->GetYaxis()->SetTitle("");
+		if (i == 0 || i == 1) {
+			hSim[i]->GetYaxis()->SetNdivisions(1005);
+		}
+		hInput[i]->SetLineWidth(2);
+
+		hSim[i]->Draw("");
+		hExp[i]->Draw("E same");
+		hInput[i]->Draw("same");
+	}
+
+//	c1->SaveAs("figures/GSepsilonT.eps");
+
+//	return;
+
+	TCanvas *c2 = new TCanvas();
+	c2->ToggleEditor();
+	c2->SetWindowSize(800, 427);
+	c2->Divide(2, 1, .001, .001);
+
+	TFile *fr3 = new TFile("figures/PRCepsilonY.root");
+	TFile *fr4 = new TFile("figures/PRCcosThetaT.root");
+
+	TH1F *hCorrExp[2];
+	TH1F *hCorrSim[2];
+	TH1F *hCorrInput[2];
+
+	hCurr = (TH1F*)fr3->Get("heeY0_0");
+	hCorrExp[0] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr4->Get("hecoskT0_0_7590");
+	hCorrExp[1] = new TH1F(*hCurr);
+
+	hCurr = (TH1F*)fr3->Get("hseY0_0");
+	hCorrSim[0] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr4->Get("hscoskT0_0_7590");
+	hCorrSim[1] = new TH1F(*hCurr);
+
+	hCurr = (TH1F*)fr3->Get("hieY0_0");
+	hCorrInput[0] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr4->Get("hicoskT0_0_7590");
+	hCorrInput[1] = new TH1F(*hCurr);
+
+	for (Int_t i = 0; i < 2; i++) {
+		c2->cd(i+1);
+
+		hCorrSim[i]->SetTitle("");
+		hCorrSim[i]->SetLineColor(kGray+1);
+		hCorrSim[i]->SetFillColor(kGray+1);
+		hCorrSim[i]->GetXaxis()->SetNdivisions(505);
+		hCorrSim[i]->GetYaxis()->SetTitle("");
+		if (i == 1) {
+			hCorrSim[i]->GetXaxis()->SetTitleSize(0.08);
+			hCorrSim[i]->GetXaxis()->SetTitleOffset(0.73);
+		}
+		hCorrInput[i]->SetLineWidth(2);
+
+		hCorrSim[i]->Draw("");
+		hCorrExp[i]->Draw("E same");
+		hCorrInput[i]->Draw("same");
+	}
+
+	c2->SaveAs("figures/GSotherCorrs.eps");
+
+}
-- 
2.18.1