FairRoot/PandaRoot
overall_perf_full.C
Go to the documentation of this file.
1 
2 // ############### TOY #################
3 
4 // *** MASS cut only
5 // 2.40 3.77 4.50 5.50
6 // ---------------------------------------------
7 // e^+e^- | 79.41 80.28 80.43 80.88
8 // \phi | 88.32 91.16 92.08 91.81
9 // \eta_c | - 87.41 89.24 89.14
10 // J/\psi(1) | - 79.88 80.55 80.84
11 // J/\psi(2) | - 81.49 82.01 82.78
12 // D^0 | - 85.16 87.36 87.04
13 // D^\pm | - 81.87 85.55 85.54
14 // D_s | - - 89.13 90.38
15 // \Lambda | 91.12 91.12 90.87 91.05
16 // \Lambda_c | - - - 89.06
17 // DPM | 3.46 11.70 15.27 21.93
18 //
19 
20 // *** intermediate solution
21 // 2.40 3.77 4.50 5.50
22 // ---------------------------------------------
23 // e^+e^- | 79.41 80.27 78.92 79.24
24 // \phi | 71.31 79.88 73.34 84.64
25 // \eta_c | - 43.07 41.60 52.76
26 // J/\psi(1) | - 79.22 78.20 77.71
27 // J/\psi(2) | - 80.94 79.92 79.91
28 // D^0 | - 50.52 38.94 55.53
29 // D^\pm | - 29.65 28.90 46.91
30 // D_s | - - 41.90 54.91
31 // \Lambda | 45.98 67.65 76.23 72.77
32 // \Lambda_c | - - - 39.52
33 // DPM | 0.10 0.16 0.18 0.12
34 
35 // *** FACTOR 1000
36 // 2.40 3.77 4.50 5.50
37 // ---------------------------------------------
38 // e^+e^- | 79.41 80.27 78.92 79.24
39 // \phi | 71.31 79.81 73.32 84.64
40 // \eta_c | - 29.33 39.63 52.70
41 // J/\psi(ee) | - 79.22 78.15 77.71
42 // J/\psi(\mu\mu) | - 80.94 79.86 79.91
43 // D^0 | - 45.08 21.55 55.50
44 // D^\pm | - 15.59 17.58 46.86
45 // D_s | - - 24.30 54.83
46 // \Lambda | 45.98 67.65 76.23 72.77
47 // \Lambda_c | - - - 36.13
48 // DPM | 0.10 0.09 0.09 0.09
49 
50 
51 // *** High EFFICIENCY >90% after combinatorics
52 // 2.40 3.77 4.50 5.50
53 // ---------------------------------------------
54 // e^+e^- | 79.41 80.27 78.92 79.24
55 // \phi | 80.84 86.55 83.84 85.83
56 // \eta_c | - 69.03 77.04 71.28
57 // J/\psi(ee) | - 79.17 78.28 77.87
58 // J/\psi(\mu\mu) | - 80.90 79.97 80.02
59 // D^0 | - 76.64 79.19 78.52
60 // D^\pm | - 71.74 76.10 73.67
61 // D_s | - - 74.75 74.42
62 // \Lambda | 84.72 83.03 83.77 81.75
63 // \Lambda_c | - - - 72.39
64 // DPM | 0.08 0.80 1.49 1.11
65 //
66 
67 // ############### FULL #################
68 
69 int cols[10] = {1,2,4,6,kOrange+5,9,kGreen+2, kBlue+3,kCyan-2,kMagenta+2};
70 
71 // *** setup style
72 void setStyle()
73 {
74  gStyle->SetPadTopMargin(0.10);
75  gStyle->SetPadBottomMargin(0.16);
76  gStyle->SetPadLeftMargin(0.12);
77  gStyle->SetPadRightMargin(0.05);
78 
79  gStyle->SetLabelSize(0.05,"X");
80  gStyle->SetLabelSize(0.05,"Y");
81 
82  gStyle->SetTitleH(0.07);
83  gStyle->SetTitleX(0.08);
84  gStyle->SetTitleOffset(1.3,"x");
85  gStyle->SetTitleXSize(0.055);
86 }
87 void renorm(double* a, int n=4)
88 {
89  for (int i=0;i<n;++i) a[i]/=100.;
90 }
91 
92 void config_graph(TGraph *g, int m=20, TString name="")
93 {
94  g->SetMarkerStyle(m);
95  g->SetMarkerColor(1);
96  g->SetMarkerSize(1.2);
97  g->SetLineWidth(2);
98  g->SetLineColor(cols[m-20]);
99  g->SetMarkerColor(cols[m-20]);
100  g->SetName(name);
101 }
102 
103 int overall_perf_full()
104 {
105  setStyle();
106 
107  TCanvas *c1=new TCanvas("c1","c1",10,10,700,500);
108  TCanvas *c2=new TCanvas("c2","c2",900,10,700,500);
109  TCanvas *c3=new TCanvas("c3","c3",900,500,700,500);
110 
111  TH1F *h=new TH1F("h","Full MC - Background fraction",100,2.0, 6.0);
112  h->SetMaximum(1);
113  h->SetMinimum(5*1e-4);
114  h->SetStats(0);
115  h->SetXTitle("#sqrt{s} [GeV]");
116 
117  TH1F *h2=new TH1F("h","Full MC - Efficiency - mass cut only",100,2.0, 6.0);
118  h2->SetMaximum(0.9);
119  h2->SetMinimum(0.0);
120  h2->SetStats(0);
121  h2->SetXTitle("#sqrt{s} [GeV]");
122 
123  TH1F *h3=new TH1F("h","Full MC - Efficiency - all cuts (high suppression)",100,2.0, 6.0);
124  h3->SetMaximum(0.9);
125  h3->SetMinimum(0.0);
126  h3->SetStats(0);
127  h3->SetXTitle("#sqrt{s} [GeV]");
128 
129  double E[4]={2.4, 3.77, 4.5, 5.5};
130  double E2[3]={ 3.77, 4.5, 5.5};
131  double E3[2]={4.5, 5.5};
132  double E4[1]={ 5.5};
133 
134  // ****
135  // *** Combinatorics and mass cut only
136  // ****
137 // & e^+e^- & \phi & \eta_c & J/\psi(e)&J/\psi(mu)& D^0 & D^\pm & D_s & \Lambda & \Lambda_c& DPM \\\hline
138 // 2.4 & 50.81 & 35.57 & - & - & - & - & - & - & 19.69 & - & 3.16\\
139 // 3.77 & 43.30 & 43.81 & 42.95 & 43.35 & 56.35 & 49.56 & 40.88 & - & 21.06 & - & 16.06\\
140 // 4.5 & 45.48 & 43.43 & 52.69 & 42.90 & 55.17 & 54.86 & 50.39 & 53.07 & 21.73 & - & 26.48\\
141 // 5.5 & 36.53 & 41.36 & 58.02 & 47.62 & 58.14 & 57.50 & 55.08 & 58.10 & 23.05 & 60.63 & 40.11\\
142 
143  double ee1[4] = { 50.81 , 43.30 , 45.48 , 36.53 };
144  double phi1[4] = { 35.57 , 43.81 , 43.43 , 41.36 };
145  double etac1[3] = { 42.95 , 52.69 , 58.02 };
146  double j11[3] = { 43.35 , 42.90 , 47.62 };
147  double j21[3] = { 56.35 , 55.17 , 58.14 };
148  double d01[3] = { 49.56 , 54.86 , 57.50 };
149  double d1[3] = { 40.88 , 50.39 , 55.08 };
150  double ds1[2] = { 53.07 , 58.10 };
151  double lam1[4] = { 19.69 , 21.06 , 21.73 , 23.05 };
152  double lamc1[1] = { 60.63 };
153  double dpm1[4] = { 3.16 , 16.06 , 26.48 , 40.11 };
154 
155  // ****
156  // *** values for high efficiency (90%)
157  // ****
158 /* & e^+e^- & \phi & \eta_c & J/\psi(e)&J/\psi(mu)& D^0 & D^\pm & D_s & \Lambda & \Lambda_c& DPM \\\hline
159  2.4 & 50.55 & 31.03 & - & - & - & - & - & - & 17.07 & - & 0.24\\
160  3.77 & 40.86 & 36.92 & 21.51 & 38.90 & 53.15 & 40.75 & 27.80 & - & 17.81 & - & 0.80\\
161  4.5 & 44.05 & 31.70 & 35.23 & 36.23 & 50.22 & 45.91 & 37.08 & 39.15 & 17.93 & - & 1.97\\
162  5.5 & 33.17 & 32.93 & 39.68 & 33.57 & 48.87 & 48.32 & 40.79 & 45.34 & 17.94 & 43.40 & 14.08\\*/
163 // double ee2[4] = { 50.55 , 40.86 , 44.05 , 33.17 };
164 // double phi2[4] = { 31.03 , 36.92 , 31.70 , 32.93 };
165 // double etac2[3] = { 21.51 , 35.23 , 39.68 };
166 // double j12[3] = { 38.90 , 36.23 , 33.57 };
167 // double j22[3] = { 53.15 , 50.22 , 48.87 };
168 // double d02[3] = { 40.75 , 45.91 , 48.32 };
169 // double d2[3] = { 27.80 , 37.08 , 40.79 };
170 // double ds2[2] = { 39.15 , 45.34 };
171 // double lam2[4] = { 17.07 , 17.81 , 17.93 , 17.94 };
172 // double lamc2[1] = { 43.40 };
173 // double dpm2[4] = { 0.24 , 0.80 , 1.97 , 14.08 };
174 
175 
176  // ****
177  // *** values for forced factor 1000 suppression
178  // ****
179 /* e^+e^- \phi \eta_cJ/\psi(ee)J/\psi(\mu\mu) D^0 D^\pm D_s \Lambda \Lambda_c DPM
180 2.4 | 50.59 30.43 - - - - - - 13.93 - 0.09
181 3.77 | 40.04 32.49 7.82 38.67 53.09 34.47 20.49 - 10.95 - 0.08
182 4.5 | 44.05 26.09 11.23 36.17 50.20 21.44 13.70 17.78 9.21 - 0.09
183 5.5 | 33.37 25.39 6.72 30.70 48.56 22.73 14.81 12.74 10.57 6.99 0.10 */
184  double ee2[4] = { 50.59 , 40.04 , 44.05 , 33.37 };
185  double phi2[4] = { 30.43 , 32.49 , 26.09 , 25.39 };
186  double etac2[3] = { 7.82 , 11.23 , 6.72 };
187  double j12[3] = { 38.67 , 36.17 , 30.70 };
188  double j22[3] = { 53.09 , 50.20 , 48.56 };
189  double d02[3] = { 34.47 , 21.44 , 22.73 };
190  double d2[3] = { 20.49 , 13.70 , 14.81 };
191  double ds2[2] = { 17.78 , 12.74 };
192  double lam2[4] = { 13.93 , 10.95 , 9.21 , 10.57 };
193  double lamc2[1] = { 6.99 };
194  double dpm2[4] = { 0.09 , 0.08 , 0.09 , 0.10 };
195 
196 
197  renorm(dpm1);renorm(dpm2);renorm(ee1);renorm(ee2);
198  renorm(phi1);renorm(phi2);renorm(etac1,3);renorm(etac2,3);
199  renorm(j11,3);renorm(j12,3);renorm(j21,3);renorm(j22,3);
200  renorm(d01,3);renorm(d02,3);renorm(d1,3);renorm(d2,3);
201  renorm(ds1,2);renorm(ds2,2);renorm(lam1);renorm(lam2);
202  renorm(lamc1,1);renorm(lamc2,1);
203 
204  TGraph *gdpm1=new TGraph(4,E,dpm1);
205  TGraph *gdpm2=new TGraph(4,E,dpm2);
206 
207  TGraph *gee1=new TGraph(4,E,ee1);
208  TGraph *gee2=new TGraph(4,E,ee2);
209 
210  TGraph *gphi1=new TGraph(4,E,phi1);
211  TGraph *gphi2=new TGraph(4,E,phi2);
212 
213  TGraph *get1=new TGraph(3,E2,etac1);
214  TGraph *get2=new TGraph(3,E2,etac2);
215 
216  TGraph *gj11=new TGraph(3,E2,j11);
217  TGraph *gj12=new TGraph(3,E2,j12);
218 
219  TGraph *gj21=new TGraph(3,E2,j21);
220  TGraph *gj22=new TGraph(3,E2,j22);
221 
222  TGraph *gd01=new TGraph(3,E2,d01);
223  TGraph *gd02=new TGraph(3,E2,d02);
224 
225  TGraph *gd1=new TGraph(3,E2,d1);
226  TGraph *gd2=new TGraph(3,E2,d2);
227 
228  TGraph *gds1=new TGraph(2,E3,ds1);
229  TGraph *gds2=new TGraph(2,E3,ds2);
230 
231  TGraph *glam1=new TGraph(4,E,lam1);
232  TGraph *glam2=new TGraph(4,E,lam2);
233 
234  TGraph *glc1=new TGraph(1,E4,lamc1);
235  TGraph *glc2=new TGraph(1,E4,lamc2);
236 
237 
238  config_graph(gdpm1,20, "gdpm1");config_graph(gdpm2,21, "gdpm2");
239  config_graph(gee1,20, "gee1");config_graph(gee2,20, "gee2");
240  config_graph(gphi1,21, "gphi1");config_graph(gphi2,21, "gphi2");
241  config_graph(gj11,22 , "gj11");config_graph(gj12,22, "gj12");
242  config_graph(gj21,23, "gj21");config_graph(gj22,23, "gj22");
243  config_graph(get1,24, "get1");config_graph(get2,24, "get2");
244  config_graph(gd01,25, "gd01");config_graph(gd02,25, "gd02");
245  config_graph(gd1,26, "gd1");config_graph(gd2,26, "gd2");
246  config_graph(gds1,27, "gds1");config_graph(gds2,27, "gds2");
247  config_graph(glam1,28, "glam1");config_graph(glam2,28, "glam2");
248  config_graph(glc1,29, "glc1");config_graph(glc2,29, "glc2");
249 
250  glc1->SetMarkerSize(2);glc2->SetMarkerSize(2);
251 
252  (c1->cd())->SetLogy(1);
253  gPad->SetGridx();
254  gPad->SetGridy();
255 
256  h->Draw();
257  gdpm1->Draw("LP same");
258  gdpm2->Draw("LP same");
259  TLegend *leg1 = new TLegend(0.12,0.75,0.55,0.9);
260  //TLegend *leg1 = new TLegend(0.12,0.83,0.58,0.9);
261  leg1->AddEntry("gdpm1","mass cut only","lp");
262  leg1->AddEntry("gdpm2","all cuts (high efficiency)","lp");
263  leg1->Draw();
264 
265  c1->SaveAs("figfull/Background_sum_supr1000_full.gif");
266  c1->SaveAs("figfull/Background_sum_supr1000_full.pdf");
267 
268 
269  c2->cd();
270  gPad->SetGridx();
271  gPad->SetGridy();
272  h2->Draw();
273  gee1->Draw("LP same");
274  gphi1->Draw("LP same");
275  get1->Draw("LP same");
276  gj11->Draw("LP same");
277  gj21->Draw("LP same");
278  gd01->Draw("LP same");
279  gd1->Draw("LP same");
280  gds1->Draw("LP same");
281  glam1->Draw("LP same");
282  glc1->Draw("LP same");
283  TLegend *leg2 = new TLegend(0.12,0.68,0.5,0.9);
284 // TLegend *leg2 = new TLegend(0.12,0.16,0.5,0.40);
285  leg2->SetNColumns(3);
286  leg2->AddEntry("glam1","#Lambda","lp");
287  leg2->AddEntry("glc1","#Lambda_{c}","lp");
288  leg2->AddEntry("gphi1","#phi","lp");
289  leg2->AddEntry("gj11","J/#psi(2e)","lp");
290  leg2->AddEntry("gj21","J/#psi(2#mu)","lp");
291  leg2->AddEntry("get1","#eta_{c}","lp");
292  leg2->AddEntry("gd01","D^{0}","lp");
293  leg2->AddEntry("gd1","D^{#pm}","lp");
294  leg2->AddEntry("gds1","D_{s}","lp");
295  leg2->AddEntry("gee1","e+e-","lp");
296  leg2->Draw();
297 /* c2->SaveAs("figfull/Efficiency_sum_mass_only_full.gif");
298  c2->SaveAs("figfull/Efficiency_sum_mass_only_full.pdf");*/
299 
300  c3->cd();
301  gPad->SetGridx();
302  gPad->SetGridy();
303  h3->Draw();
304  gee2->Draw("LP same");
305  gphi2->Draw("LP same");
306  get2->Draw("LP same");
307  gj12->Draw("LP same");
308  gj22->Draw("LP same");
309  gd02->Draw("LP same");
310  gd2->Draw("LP same");
311  gds2->Draw("LP same");
312  glam2->Draw("LP same");
313  glc2->Draw("LP same");
314 
315  TLegend *leg3 = new TLegend(0.12,0.68,0.5,0.9);
316 // TLegend *leg3 = new TLegend(0.12,0.16,0.45,0.45);
317  leg3->SetNColumns(3);
318 /* TLegend *leg3 = new TLegend(0.12,0.66,0.5,0.9);
319  leg3->SetNColumns(3);*/
320  leg3->AddEntry("glam2","#Lambda","lp");
321  leg3->AddEntry("glc2","#Lambda_{c}","lp");
322  leg3->AddEntry("gphi2","#phi","lp");
323  leg3->AddEntry("gj12","J/#psi(2e)","lp");
324  leg3->AddEntry("gj22","J/#psi(2#mu)","lp");
325  leg3->AddEntry("get2","#eta_{c}","lp");
326  leg3->AddEntry("gd02","D^{0}","lp");
327  leg3->AddEntry("gd2","D^{#pm}","lp");
328  leg3->AddEntry("gds2","D_{s}","lp");
329  leg3->AddEntry("gee2","e+e-","lp");
330  leg3->Draw();
331  c3->SaveAs("figfull/Efficiency_sum_opt_supr1000_full.gif");
332  c3->SaveAs("figfull/Efficiency_sum_opt_supr1000_full.pdf");
333  return 0;
334 }
335 
336 
337 
338 
339 
340 
341 
config_graph
void config_graph(TGraph *g, int m=0, TString name="")
Definition: evaltrig.C:73
setStyle
void setStyle()
Definition: crosstag.C:39
SetLogy
gPad SetLogy()
i
Int_t i
Definition: run_full.C:25
m
__m128 m
Definition: P4_F32vec4.h:28
n
int n
Definition: HypStatDecay.cxx:21
g
TFile * g
Definition: GammaSpectraAnalysis.C:22
c2
c2
Definition: plot_dirc.C:39
a
Int_t a
Definition: anaLmdDigi.C:126
cols
int cols[10]
Definition: evaltrig.C:69
c1
c1
Definition: plot_dirc.C:35
h3
TH1F * h3
Definition: anal_point_fwendcap.C:83
renorm
void renorm(double *a, int n=4)
Definition: overall_perf_full.C:87
name
TString name
Definition: createRootGeoFile.C:82
c3
c3
Definition: plot_dirc.C:50
h2
TH1F * h2
Definition: anal_hit_digi_cluster_fwendcap.C:122
overall_perf_full
int overall_perf_full()
Definition: overall_perf_full.C:103
Generated on Wed Apr 3 2019 10:02:48 for FairRoot/PandaRoot by   doxygen 1.8.5