Commit bc0c5e11 authored by Vratislav Chudoba's avatar Vratislav Chudoba

Parameters for simulation in (75;90) modified.

parent 7a0fcc97
...@@ -198,7 +198,7 @@ void BeAnalysis::SetNoSimEvents(Long64_t *noSimEvents) { ...@@ -198,7 +198,7 @@ void BeAnalysis::SetNoSimEvents(Long64_t *noSimEvents) {
if (!noSimEvents) { if (!noSimEvents) {
for (Int_t i = 0; i < 6; i++) { for (Int_t i = 0; i < 6; i++) {
sEvents[i] = sMaxEvents; sEvents[i] = sMaxEvents;
cout << sEvents[i] << endl; // cout << sEvents[i] << endl;
} }
Warning("BeAnalysis::SetNoSimEvents", "Default numbers of simulated events were set."); Warning("BeAnalysis::SetNoSimEvents", "Default numbers of simulated events were set.");
return; return;
...@@ -206,7 +206,7 @@ void BeAnalysis::SetNoSimEvents(Long64_t *noSimEvents) { ...@@ -206,7 +206,7 @@ void BeAnalysis::SetNoSimEvents(Long64_t *noSimEvents) {
for (Int_t i = 0; i < 6; i++) { for (Int_t i = 0; i < 6; i++) {
sEvents[i] = noSimEvents[i]; sEvents[i] = noSimEvents[i];
cout << sEvents[i] << endl; // cout << sEvents[i] << endl;
} }
return; return;
......
void ES1_7590() void ES1_7590()
{ {
const Bool_t saveHists = 0;
// const Bool_t saveHists = 1;
TStopwatch sw; TStopwatch sw;
sw.Start(); sw.Start();
...@@ -8,27 +11,30 @@ void ES1_7590() ...@@ -8,27 +11,30 @@ void ES1_7590()
ana.SetCMAngularRange(75, 90); ana.SetCMAngularRange(75, 90);
ana.SetFigures("figures/7590/", ".eps"); // ana.SetFigures("figures/7590/", ".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] = {1,0,0,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,0,1,0,0};
ana.SetSpectra(); TString simEt[5] = {"", "f6BeIM>1.5 && f6BeIM<2.0", "", "", ""};
// ana.SetEpsilonTintervals(epsilonT); ana.SetSimCuts(simEt);
ana.SetEpsilonTintervals();
// ana.SetSpectra();
ana.SetEpsilonTintervals(epsilonT);
// ana.SetEpsilonTintervals();
// ana.SetCosThetaTkIntervals(cosThetaT); // ana.SetCosThetaTkIntervals(cosThetaT);
ana.SetCosThetaTkIntervals(cosThetaT); // 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);
...@@ -37,7 +43,8 @@ void ES1_7590() ...@@ -37,7 +43,8 @@ void ES1_7590()
ana.SetNoExpEvents(eEvents); ana.SetNoExpEvents(eEvents);
const Long64_t sMaxEvents = 3000000; const Long64_t sMaxEvents = 3000000;
Long64_t sEvents[6] = {2050000, 2050000, 2050000, 2050000, 1900000, 1950000}; // Long64_t sEvents[6] = {2050000, 2050000, 2050000, 2050000, 1900000, 1950000};
Long64_t sEvents[6] = {2050000, 2050000, 2050000, 2050000, 2050000, 2050000};
ana.SetNoSimEvents(sEvents); ana.SetNoSimEvents(sEvents);
Int_t lowSimFile[6] = {0, 0, 0, 0, 0, 0}; Int_t lowSimFile[6] = {0, 0, 0, 0, 0, 0};
...@@ -63,13 +70,15 @@ void ES1_7590() ...@@ -63,13 +70,15 @@ void ES1_7590()
ana.OpenSimChains(); ana.OpenSimChains();
ana.SetRatiosGStoEX("sRatio>0.075 && sRatio<0.125", ana.SetRatiosGStoEX("sRatio>0.080 && sRatio<0.130",
"sRatio>0.076 && sRatio<0.126", "sRatio>0.080 && sRatio<0.130",
"sRatio>0.070 && sRatio<0.120", // "sRatio>0.070 && sRatio<0.120",
"sRatio>0.070 && sRatio<0.120", // "sRatio>0.070 && sRatio<0.120",
"sRatio>0.075 && sRatio<0.125", "sRatio>0.080 && sRatio<0.130",
"sRatio>0.076 && sRatio<0.126"); "sRatio>0.080 && sRatio<0.130",
ana.SetCuts(); "sRatio>0.080 && sRatio<0.130",
"sRatio>0.080 && sRatio<0.130");
// ana.SetCuts();
Long64_t eEventsECuts[5][6] = {eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, Long64_t eEventsECuts[5][6] = {eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents,
eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents, eMaxEvents,
...@@ -80,7 +89,8 @@ void ES1_7590() ...@@ -80,7 +89,8 @@ void ES1_7590()
const Long64_t sMaxEvents = 3000000; const Long64_t sMaxEvents = 3000000;
Long64_t sEventsECuts[5][6] = {2020000, 2040000, 2130000, 2080000, 2060000, 2080000, Long64_t sEventsECuts[5][6] = {2020000, 2040000, 2130000, 2080000, 2060000, 2080000,
2000000, 2020000, 2000000, 2010000, 1990000, 1970000, // 2000000, 2020000, 2000000, 2010000, 1990000, 1970000,
2300000, 2300000, 2300000, 2300000, 2300000, 2300000,
2320000, 2300000, 2340000, 2310000, 2300000, 2290000, 2320000, 2300000, 2340000, 2310000, 2300000, 2290000,
2060000, 2070000, 2110000, 2110000, 2070000, 2060000, 2060000, 2070000, 2110000, 2110000, 2070000, 2060000,
2050000, 2050000, 2110000, 2110000, 2050000, 2050000}; 2050000, 2050000, 2110000, 2110000, 2050000, 2050000};
...@@ -114,6 +124,13 @@ void ES1_7590() ...@@ -114,6 +124,13 @@ void ES1_7590()
sw.Print(); sw.Print();
sw.Continue(); sw.Continue();
if (saveHists) {
ana.SaveSpectrumHist(0);
ana.SaveEpsilonTHist(0, 5);
ana.SaveCosThetaTHist(0, 0);
ana.SaveEpsilonYHist(0, 0);
}
Info("ES1_7590.cxx", "Finished."); Info("ES1_7590.cxx", "Finished.");
} }
......
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