Commit d178d5c3 authored by Vratislav Chudoba's avatar Vratislav Chudoba

Parameters for simulation in (45;60) modified.

parent 260c021d
void ES1_4560() void ES1_4560()
{ {
/* accessible sRatio in simulated beams - (0.05,0.15)
*
*
*/
const Bool_t saveHists = 0;
// const Bool_t saveHists = 1;
TStopwatch sw; TStopwatch sw;
sw.Start(); sw.Start();
...@@ -8,35 +16,41 @@ void ES1_4560() ...@@ -8,35 +16,41 @@ void ES1_4560()
ana.SetCMAngularRange(45, 60); 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(); 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 cosThetaT[6] = {1,0,0,0,0};
Bool_t epsilonY[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 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(epsilonT);
ana.SetEpsilonTintervals(); // ana.SetEpsilonTintervals();
ana.SetCosThetaTkIntervals(cosThetaT); // ana.SetCosThetaTkIntervals(cosThetaT);
ana.SetCosThetaTkIntervals(); // ana.SetCosThetaTkIntervals();
ana.SetEpsilonYintervals(epsilonY); // ana.SetEpsilonYintervals(epsilonY);
ana.SetEpsilonYintervals(); // ana.SetEpsilonYintervals();
ana.SetCosThetaYkIntervals(cosThetaY); // ana.SetCosThetaYkIntervals(cosThetaY);
ana.SetCosThetaYkIntervals(); // ana.SetCosThetaYkIntervals();
ana.SetThetaATintervals(thetaAT); // ana.SetThetaATintervals(thetaAT);
ana.SetThetaATintervals(); // ana.SetThetaATintervals();
ana.SetExpChain("../../../be/rootdata/correlations/v5_exp/Be.", 0, 40); ana.SetExpChain("../../../be/rootdata/correlations/v5_exp/Be.", 0, 40);
const Long64_t eMaxEvents = 5000000; 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); ana.SetNoExpEvents(eEvents);
const Long64_t sMaxEvents = 3000000; const Long64_t sMaxEvents = 3000000;
...@@ -72,6 +86,7 @@ void ES1_4560() ...@@ -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.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", // ana.SetRatiosGStoEX("sRatio>0.051 && sRatio<0.101",
// "sRatio>0.050 && sRatio<0.100", // "sRatio>0.050 && sRatio<0.100",
...@@ -79,10 +94,11 @@ void ES1_4560() ...@@ -79,10 +94,11 @@ void ES1_4560()
// "sRatio>0.048 && sRatio<0.098", // "sRatio>0.048 && sRatio<0.098",
// "sRatio>0.051 && sRatio<0.101", // "sRatio>0.051 && sRatio<0.101",
// "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, 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, 2400000, 2400000, 2400000, 2400000, 2400000, 2400000,
2900000, 3000000, 2900000, 3000000, 2900000, 3000000, 2900000, 3000000, 2900000, 3000000, 2900000, 3000000,
3200000, 3300000, 3300000, 3400000, 3300000, 3400000}; 3200000, 3300000, 3300000, 3400000, 3300000, 3400000};
...@@ -126,6 +142,13 @@ void ES1_4560() ...@@ -126,6 +142,13 @@ void ES1_4560()
sw.Print(); sw.Print();
sw.Continue(); sw.Continue();
if (saveHists) {
// ana.SaveSpectrumHist(0);
ana.SaveEpsilonTHist(0, 0);
ana.SaveCosThetaTHist(0, 3);
}
Info("ES1_4560.cxx", "Finished."); Info("ES1_4560.cxx", "Finished.");
return; return;
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment