From 260c021d1f8c8f21665c53306a7742665bd97d24 Mon Sep 17 00:00:00 2001 From: Vratislav Chudoba Date: Wed, 2 Aug 2017 17:50:06 +0300 Subject: [PATCH] Parameters for simulation in (90;120) modified. --- macros/BeCorrPRC/ES1_90120.cxx | 56 ++++++++++++++++++++-------------- 1 file changed, 33 insertions(+), 23 deletions(-) diff --git a/macros/BeCorrPRC/ES1_90120.cxx b/macros/BeCorrPRC/ES1_90120.cxx index edf065a..90e2bec 100644 --- a/macros/BeCorrPRC/ES1_90120.cxx +++ b/macros/BeCorrPRC/ES1_90120.cxx @@ -1,8 +1,8 @@ void ES1_90120() { - //open chain with Be events only and draw E_Be - //order: al0, noal0, al180, noal180 +// const Bool_t saveHists = 0; + const Bool_t saveHists = 0; TStopwatch sw; sw.Start(); @@ -11,27 +11,30 @@ void ES1_90120() ana.SetCMAngularRange(90, 120); -// ana.SetFigures("figures/", ".eps"); +// ana.SetFigures("figures/90120/", ".eps"); // ana.SetChainsToDraw(0, 0, 0, 0, 1, 1); ana.SetChainsToDraw(); - Bool_t epsilonT[6] = {0,0,1,0,0}; + Bool_t epsilonT[6] = {0,1,0,0,0}; Bool_t cosThetaT[6] = {0,0,1,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}; - ana.SetSpectra(); + TString simEt[5] = {"", "f6BeIM>1.40 && f6BeIM<1.95", "", "", ""}; + ana.SetSimCuts(simEt); + +// 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); @@ -41,7 +44,8 @@ void ES1_90120() ana.SetNoExpEvents(eEvents); const Long64_t sMaxEvents = 3000000; - Long64_t sEvents[6] = {2280000, 2290000, 2100000, 2300000, 2330000, 2300000}; +// Long64_t sEvents[6] = {2280000, 2290000, 2100000, 2300000, 2330000, 2300000}; + Long64_t sEvents[6] = {2200000, 2200000, 2200000, 2200000, 2200000, 2200000}; ana.SetNoSimEvents(sEvents); Int_t lowSimFile[6] = {0, 0, 0, 0, 0, 0}; @@ -66,13 +70,13 @@ void ES1_90120() ana.OpenSimChains(); - ana.SetRatiosGStoEX("sRatio>0.108 && sRatio<0.158", - "sRatio>0.104 && sRatio<0.154", - "sRatio>0.0109 && sRatio<0.159", - "sRatio>0.103 && sRatio<0.153", - "sRatio>0.101 && sRatio<0.151", - "sRatio>0.104 && sRatio<0.154"); - ana.SetCuts(); + ana.SetRatiosGStoEX("sRatio>0.115 && sRatio<0.165", + "sRatio>0.115 && sRatio<0.165", + "sRatio>0.115 && sRatio<0.165", + "sRatio>0.115 && sRatio<0.165", + "sRatio>0.115 && sRatio<0.165", + "sRatio>0.115 && sRatio<0.165"); +// ana.SetCuts(); Long64_t eEventsECuts[5][6] = {eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, @@ -82,7 +86,8 @@ void ES1_90120() ana.ExpEventsECuts(eEventsECuts); Long64_t sEventsECuts[5][6] = {2300000, 2310000, 2130000, 2330000, 2320000, 2310000, - 2260000, 2270000, 2010000, 2280000, 2280000, 2280000, +// 2260000, 2270000, 2010000, 2280000, 2280000, 2280000, + 2280000, 2280000, 2280000, 2280000, 2280000, 2280000, sMaxEvents, sMaxEvents, 2290000, sMaxEvents, sMaxEvents, sMaxEvents, 3270000, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, 2320000, 2330000, 2100000, 2340000, 2370000, 2340000}; @@ -116,6 +121,11 @@ void ES1_90120() sw.Print(); sw.Continue(); + if (saveHists) { + ana.SaveSpectrumHist(0); + ana.SaveEpsilonTHist(0, 0); + } + Info("ES1_90120.cxx", "Finished."); return; -- 2.18.1