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

---
 macros/BeCorrPRC/ES1_6075.cxx | 62 +++++++++++++++++++++--------------
 1 file changed, 38 insertions(+), 24 deletions(-)

diff --git a/macros/BeCorrPRC/ES1_6075.cxx b/macros/BeCorrPRC/ES1_6075.cxx
index 8f0fb8a..76d86a9 100644
--- a/macros/BeCorrPRC/ES1_6075.cxx
+++ b/macros/BeCorrPRC/ES1_6075.cxx
@@ -1,6 +1,9 @@
 void ES1_6075()
 {
 
+	const Bool_t saveHists = 0;
+//	const Bool_t saveHists = 1;
+
 	TStopwatch sw;
 	sw.Start();
 
@@ -8,31 +11,34 @@ void ES1_6075()
 
 	ana.SetCMAngularRange(60, 75);
 
-	ana.SetFigures("figures/6075/", ".eps");
+//	ana.SetFigures("figures/6075/", ".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();
+//	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.SetNoBinsSpectra(200);
 
 	ana.SetExpChain("../../../be/rootdata/correlations/v5_exp/Be.", 0, 40);
 
@@ -41,7 +47,8 @@ void ES1_6075()
 	ana.SetNoExpEvents(eEvents);
 
 	const Long64_t sMaxEvents = 3000000;
-	Long64_t sEvents[6] = {sMaxEvents, 2700000, sMaxEvents, 2700000, 2690000, 2710000};
+//	Long64_t sEvents[6] = {sMaxEvents, 2700000, sMaxEvents, 2700000, 2690000, 2710000};
+	Long64_t sEvents[6] = {sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents};
 	ana.SetNoSimEvents(sEvents);
 
 	Int_t lowSimFile[6] = {0, 0, 0, 0, 0, 0};
@@ -65,23 +72,25 @@ void ES1_6075()
 	ana.OpenExpChain();
 	ana.OpenSimChains();
 
-	ana.SetRatiosGStoEX("sRatio>0.051 && sRatio<0.101",
-			"sRatio>0.050 && sRatio<0.100",
-			"sRatio>0.047 && sRatio<0.097",
-			"sRatio>0.048 && sRatio<0.098",
-			"sRatio>0.051 && sRatio<0.101",
-			"sRatio>0.051 && sRatio<0.101");
-	ana.SetCuts();
+	ana.SetRatiosGStoEX("sRatio>0.060 && sRatio<0.110",
+			"sRatio>0.060 && sRatio<0.110",
+			"sRatio>0.060 && sRatio<0.110",
+			"sRatio>0.060 && sRatio<0.110",
+			"sRatio>0.060 && sRatio<0.110",
+			"sRatio>0.060 && sRatio<0.110");
+//	ana.SetCuts();
 
 	Long64_t eEventsECuts[5][6] = {3950000, 3950000, 3800000, 3850000, 3950000, 3900000,
-									eMaxEvents, eMaxEvents, 4100000, eMaxEvents, eMaxEvents, eMaxEvents,
+//									eMaxEvents, eMaxEvents, 4100000, eMaxEvents, eMaxEvents, eMaxEvents,
+									3500000, 3500000, 3500000, 3500000, 3500000, 3500000,
 									3050000, 3150000, 3200000, 3200000, 3250000, 3350000,
 									3850000, 3950000, 3900000, 3900000, 4000000, 4000000,
 									4100000, 4100000, 4000000, 4100000, eMaxEvents, eMaxEvents};
 	ana.ExpEventsECuts(eEventsECuts);
 
 	Long64_t sEventsECuts[5][6] = {sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents,
-										sMaxEvents, 2630000, sMaxEvents, sMaxEvents, 2660000, 2660000,
+//										sMaxEvents, 2630000, sMaxEvents, sMaxEvents, 2660000, 2660000,
+										sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, 2660000, 2660000,
 										sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents,
 										sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents,
 										sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents};
@@ -118,6 +127,11 @@ void ES1_6075()
 	sw.Print();
 	sw.Continue();
 
+	if (saveHists) {
+		ana.SaveSpectrumHist(0);
+		ana.SaveEpsilonTHist(0, 5);
+	}
+
 	Info("ES1_6075.cxx", "Finished.");
 
 	return;
-- 
2.18.1