diff --git a/macros/BeCorrPRC/ES1_4560.cxx b/macros/BeCorrPRC/ES1_4560.cxx index fd4848d78d8af39cdfafa7103c86e7aa409e7b9d..26557f253aec88cace446717365ae687422f3300 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;