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

---
 macros/BeCorrPRC/ES1_4560.cxx | 59 ++++++++++++++++++++++++-----------
 1 file changed, 41 insertions(+), 18 deletions(-)

diff --git a/macros/BeCorrPRC/ES1_4560.cxx b/macros/BeCorrPRC/ES1_4560.cxx
index fd4848d..26557f2 100644
--- a/macros/BeCorrPRC/ES1_4560.cxx
+++ b/macros/BeCorrPRC/ES1_4560.cxx
@@ -1,6 +1,14 @@
 void ES1_4560()
 {
 
+	/* accessible sRatio in simulated beams - (0.05,0.15)
+	 *
+	 *
+	 */
+
+	const Bool_t saveHists = 0;
+//	const Bool_t saveHists = 1;
+
 	TStopwatch sw;
 	sw.Start();
 
@@ -8,35 +16,41 @@ void ES1_4560()
 
 	ana.SetCMAngularRange(45, 60);
 
-	ana.SetFigures("figures/4560/", ".eps");
+//	ana.SetFigures("figures/4560/", ".eps");
 
-	//	ana.SetChainsToDraw(0, 0, 0, 0, 1, 1);
+//	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 epsilonT[6] = {0,1,0,0,0};
+	Bool_t cosThetaT[6] = {1,0,0,0,0};
 	Bool_t epsilonY[6] = {0,0,1,0,0};
 	Bool_t cosThetaY[6] = {0,0,1,0,0};
-	Bool_t thetaAT[6] = {0,0,1,0,0};
+	Bool_t thetaAT[6] = {0,1,0,0,0};
+
+	TString simEt[5] = {"", "f6BeIM>1.6 && f6BeIM<2.1", "", "", ""};
+//	TString simEt[5];
+	ana.SetSimCuts(simEt);
 
-	ana.SetSpectra();
+//	ana.SetNoBinsSpectra(200);
+
+//	ana.SetSpectra();
 	ana.SetEpsilonTintervals(epsilonT);
-	ana.SetEpsilonTintervals();
-	ana.SetCosThetaTkIntervals(cosThetaT);
-	ana.SetCosThetaTkIntervals();
-	ana.SetEpsilonYintervals(epsilonY);
-	ana.SetEpsilonYintervals();
-	ana.SetCosThetaYkIntervals(cosThetaY);
-	ana.SetCosThetaYkIntervals();
-	ana.SetThetaATintervals(thetaAT);
-	ana.SetThetaATintervals();
+//	ana.SetEpsilonTintervals();
+//	ana.SetCosThetaTkIntervals(cosThetaT);
+//	ana.SetCosThetaTkIntervals();
+//	ana.SetEpsilonYintervals(epsilonY);
+//	ana.SetEpsilonYintervals();
+//	ana.SetCosThetaYkIntervals(cosThetaY);
+//	ana.SetCosThetaYkIntervals();
+//	ana.SetThetaATintervals(thetaAT);
+//	ana.SetThetaATintervals();
 
 
 
 	ana.SetExpChain("../../../be/rootdata/correlations/v5_exp/Be.", 0, 40);
 
 	const Long64_t eMaxEvents = 5000000;
-	Long64_t eEvents[6] = {3150000, 3150000, 3150000, 3150000, 3150000, 3150000};
+	Long64_t eEvents[6] = {3150000, 3150000, 3150000, 3150000, 3150000, /*1115000*/ 3150000};
 	ana.SetNoExpEvents(eEvents);
 
 	const Long64_t sMaxEvents = 3000000;
@@ -72,6 +86,7 @@ void ES1_4560()
 			"sRatio>0.05 && sRatio<0.095",
 			"sRatio>0.05 && sRatio<0.095",
 			"sRatio>0.05 && sRatio<0.095");
+//			"sRatio>0.1 && sRatio<0.115");
 
 //	ana.SetRatiosGStoEX("sRatio>0.051 && sRatio<0.101",
 //			"sRatio>0.050 && sRatio<0.100",
@@ -79,10 +94,11 @@ void ES1_4560()
 //			"sRatio>0.048 && sRatio<0.098",
 //			"sRatio>0.051 && sRatio<0.101",
 //			"sRatio>0.051 && sRatio<0.101");
-	ana.SetCuts();
+//	ana.SetCuts();
 
 	Long64_t eEventsECuts[5][6] = {2900000, 2900000, 2900000, 2900000, 2900000, 2900000,
-										3400000, 3400000, 3400000, 3500000, 3400000, 3400000,
+//										3400000, 3400000, 3400000, 3500000, 3400000, 3400000,
+										2100000, 2100000, 2100000, 2100000, 2100000, 2100000,
 										2400000, 2400000, 2400000, 2400000, 2400000, 2400000,
 										2900000, 3000000, 2900000, 3000000, 2900000, 3000000,
 										3200000, 3300000, 3300000, 3400000, 3300000, 3400000};
@@ -126,6 +142,13 @@ void ES1_4560()
 	sw.Print();
 	sw.Continue();
 
+	if (saveHists) {
+//		ana.SaveSpectrumHist(0);
+		ana.SaveEpsilonTHist(0, 0);
+		ana.SaveCosThetaTHist(0, 3);
+	}
+
+
 	Info("ES1_4560.cxx", "Finished.");
 
 	return;
-- 
2.18.1