From bc0c5e11c4d35c2fa99c2a6f655c0a03f7333b5b Mon Sep 17 00:00:00 2001
From: Vratislav Chudoba <chudoba@jinr.ru>
Date: Wed, 2 Aug 2017 17:00:17 +0300
Subject: [PATCH] Parameters for simulation in (75;90) modified.

---
 Be/BeAnalysis.cpp             |  4 +--
 macros/BeCorrPRC/ES1_7590.cxx | 61 ++++++++++++++++++++++-------------
 2 files changed, 41 insertions(+), 24 deletions(-)

diff --git a/Be/BeAnalysis.cpp b/Be/BeAnalysis.cpp
index 24e10be..68af630 100644
--- a/Be/BeAnalysis.cpp
+++ b/Be/BeAnalysis.cpp
@@ -198,7 +198,7 @@ void BeAnalysis::SetNoSimEvents(Long64_t *noSimEvents) {
 	if (!noSimEvents) {
 		for (Int_t i = 0; i < 6; i++) {
 			sEvents[i] = sMaxEvents;
-			cout << sEvents[i] << endl;
+//			cout << sEvents[i] << endl;
 		}
 		Warning("BeAnalysis::SetNoSimEvents", "Default numbers of simulated events were set.");
 		return;
@@ -206,7 +206,7 @@ void BeAnalysis::SetNoSimEvents(Long64_t *noSimEvents) {
 
 	for (Int_t i = 0; i < 6; i++) {
 		sEvents[i] = noSimEvents[i];
-		cout << sEvents[i] << endl;
+//		cout << sEvents[i] << endl;
 	}
 	return;
 
diff --git a/macros/BeCorrPRC/ES1_7590.cxx b/macros/BeCorrPRC/ES1_7590.cxx
index 4cb18b2..86abae7 100644
--- a/macros/BeCorrPRC/ES1_7590.cxx
+++ b/macros/BeCorrPRC/ES1_7590.cxx
@@ -1,6 +1,9 @@
 void ES1_7590()
 {
 
+	const Bool_t saveHists = 0;
+//	const Bool_t saveHists = 1;
+
 	TStopwatch sw;
 	sw.Start();
 
@@ -8,27 +11,30 @@ void ES1_7590()
 
 	ana.SetCMAngularRange(75, 90);
 
-	ana.SetFigures("figures/7590/", ".eps");
+//	ana.SetFigures("figures/7590/", ".eps");
 
 //	ana.SetChainsToDraw(0, 0, 0, 0, 1, 1);
 	ana.SetChainsToDraw();
-	Bool_t epsilonT[6] = {0,0,1,0,0};
-	Bool_t cosThetaT[6] = {0,0,1,0,0};
-	Bool_t epsilonY[6] = {0,0,1,0,0};
+	Bool_t epsilonT[6] = {0,1,0,0,0};
+	Bool_t cosThetaT[6] = {1,0,0,0,0};
+	Bool_t epsilonY[6] = {1,0,0,0,0};
 	Bool_t cosThetaY[6] = {0,0,1,0,0};
 	Bool_t thetaAT[6] = {0,0,1,0,0};
 
-	ana.SetSpectra();
-//	ana.SetEpsilonTintervals(epsilonT);
-	ana.SetEpsilonTintervals();
+	TString simEt[5] = {"", "f6BeIM>1.5 && f6BeIM<2.0", "", "", ""};
+	ana.SetSimCuts(simEt);
+
+//	ana.SetSpectra();
+	ana.SetEpsilonTintervals(epsilonT);
+//	ana.SetEpsilonTintervals();
 //	ana.SetCosThetaTkIntervals(cosThetaT);
-	ana.SetCosThetaTkIntervals(cosThetaT);
+//	ana.SetCosThetaTkIntervals();
 //	ana.SetEpsilonYintervals(epsilonY);
-	ana.SetEpsilonYintervals();
-	ana.SetCosThetaYkIntervals(cosThetaY);
-	ana.SetCosThetaYkIntervals();
-	ana.SetThetaATintervals(thetaAT);
-	ana.SetThetaATintervals();
+//	ana.SetEpsilonYintervals();
+//	ana.SetCosThetaYkIntervals(cosThetaY);
+//	ana.SetCosThetaYkIntervals();
+//	ana.SetThetaATintervals(thetaAT);
+//	ana.SetThetaATintervals();
 
 	ana.SetExpChain("../../../be/rootdata/correlations/v5_exp/Be.", 0, 40);
 
@@ -37,7 +43,8 @@ void ES1_7590()
 	ana.SetNoExpEvents(eEvents);
 
 	const Long64_t sMaxEvents = 3000000;
-	Long64_t sEvents[6] = {2050000, 2050000, 2050000, 2050000, 1900000, 1950000};
+//	Long64_t sEvents[6] = {2050000, 2050000, 2050000, 2050000, 1900000, 1950000};
+	Long64_t sEvents[6] = {2050000, 2050000, 2050000, 2050000, 2050000, 2050000};
 	ana.SetNoSimEvents(sEvents);
 
 	Int_t lowSimFile[6] = {0, 0, 0, 0, 0, 0};
@@ -63,13 +70,15 @@ void ES1_7590()
 	ana.OpenSimChains();
 
 
-	ana.SetRatiosGStoEX("sRatio>0.075 && sRatio<0.125",
-						"sRatio>0.076 && sRatio<0.126",
-						"sRatio>0.070 && sRatio<0.120",
-						"sRatio>0.070 && sRatio<0.120",
-						"sRatio>0.075 && sRatio<0.125",
-						"sRatio>0.076 && sRatio<0.126");
-	ana.SetCuts();
+	ana.SetRatiosGStoEX("sRatio>0.080 && sRatio<0.130",
+						"sRatio>0.080 && sRatio<0.130",
+//						"sRatio>0.070 && sRatio<0.120",
+//						"sRatio>0.070 && sRatio<0.120",
+						"sRatio>0.080 && sRatio<0.130",
+						"sRatio>0.080 && sRatio<0.130",
+						"sRatio>0.080 && sRatio<0.130",
+						"sRatio>0.080 && sRatio<0.130");
+//	ana.SetCuts();
 
 	Long64_t eEventsECuts[5][6] = {eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents,
 									eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents,
@@ -80,7 +89,8 @@ void ES1_7590()
 
 	const Long64_t sMaxEvents = 3000000;
 	Long64_t sEventsECuts[5][6] = {2020000, 2040000, 2130000, 2080000, 2060000, 2080000,
-									2000000, 2020000, 2000000, 2010000, 1990000, 1970000,
+//									2000000, 2020000, 2000000, 2010000, 1990000, 1970000,
+									2300000, 2300000, 2300000, 2300000, 2300000, 2300000,
 									2320000, 2300000, 2340000, 2310000, 2300000, 2290000,
 									2060000, 2070000, 2110000, 2110000, 2070000, 2060000,
 									2050000, 2050000, 2110000, 2110000, 2050000, 2050000};
@@ -114,6 +124,13 @@ void ES1_7590()
 	sw.Print();
 	sw.Continue();
 
+	if (saveHists) {
+		ana.SaveSpectrumHist(0);
+		ana.SaveEpsilonTHist(0, 5);
+		ana.SaveCosThetaTHist(0, 0);
+		ana.SaveEpsilonYHist(0, 0);
+	}
+
 	Info("ES1_7590.cxx", "Finished.");
 
 }
-- 
2.18.1