FairRoot/PandaRoot
draw_bands.C
Go to the documentation of this file.
1 {
2 
3  TString fundedx = "[0] * (1./x)**2 + [1] * TMath::Log(x) + [2]";
4 
5  // TString funsig = "pol4(3)";
6  TString funsig = "([3] * (1./x)**2 + [4] * TMath::Log(x) + [5])";
7 
8  // ====================================
9 
10  TCanvas *c = new TCanvas("c", "dedx_p");
11  c->Divide(3, 2);
12 
13  TFile infile("dedx_out.root", "READ");
14 
15  TH2F *h1 = (TH2F*) infile.Get("hdedx_p_e");
16  TH2F *h2 = (TH2F*) infile.Get("hdedx_p_mu");
17  TH2F *h3 = (TH2F*) infile.Get("hdedx_p_pi");
18  TH2F *h4 = (TH2F*) infile.Get("hdedx_p_k");
19  TH2F *h5 = (TH2F*) infile.Get("hdedx_p_p");
20 
21  gStyle->SetPalette(1, 0);
22  c->cd(1); h1->Draw("colz");
23  c->cd(2); h2->Draw("colz");
24  c->cd(3); h3->Draw("colz");
25  c->cd(4); h4->Draw("colz");
26  c->cd(5); h5->Draw("colz");
27  c->cd(6);
28  h1->Draw();
29  h2->Draw("SAME");
30  h3->Draw("SAME");
31  h4->Draw("SAME");
32  h5->Draw("SAME");
33 
34  TString funsum = fundedx + " + " + funsig;
35  TString fundiff = fundedx + " - " + funsig;
36 
37  TF1 *fdedx1 = new TF1("fdedx1", fundedx, 0., 1.5);
38  fdedx1->SetParameters(-0.0051344, 0.163231, 7.30319);
39  TF1 *fsum1 = new TF1("fsum1", funsum, 0., 1.5);
40  fsum1->SetParameters(-0.0051344, 0.163231, 7.30319, 0.00428134, -0.00359152, 0.655438);
41  TF1 *fdiff1 = new TF1("fdiff1", fundiff, 0., 1.5);
42  fdiff1->SetParameters(-0.0051344, 0.163231, 7.30319, 0.00428134, -0.00359152, 0.655438);
43  c->cd(1);
44  fdedx1->SetLineColor(kOrange);
45  fdedx1->Draw("SAME");
46  fsum1->Draw("SAME");
47  fdiff1->Draw("SAME");
48 
49  TF1 *fdedx2 = new TF1("fdedx2", fundedx, 0., 1.5);
50  fdedx2->SetParameters(0.0737857, 0.643465, 5.42333);
51  TF1 *fsum2 = new TF1("fsum2", funsum, 0., 1.5);
52  fsum2->SetParameters(0.0737857, 0.643465, 5.42333, 0.0193587, -0.0174329, 0.539616);
53  TF1 *fdiff2 = new TF1("fdiff2", fundiff, 0., 1.5);
54  fdiff2->SetParameters(0.0737857, 0.643465, 5.42333, 0.0193587, -0.0174329, 0.539616);
55  c->cd(2);
56  fdedx2->SetLineColor(kCyan);
57  fdedx2->Draw("SAME");
58  fsum2->Draw("SAME");
59  fdiff2->Draw("SAME");
60 
61  TF1 *fdedx3 = new TF1("fdedx3", fundedx, 0., 1.5);
62  fdedx3->SetParameters(0.124558, 0.585125, 5.20178);
63  TF1 *fsum3 = new TF1("fsum3", funsum, 0., 1.5);
64  fsum3->SetParameters(0.124558, 0.585125, 5.20178, 0.0395417, -0.0443629, 0.502671);
65  TF1 *fdiff3 = new TF1("fdiff3", fundiff, 0., 1.5);
66  fdiff3->SetParameters(0.124558, 0.585125, 5.20178, 0.0395417, -0.0443629, 0.502671);
67  c->cd(3);
68  fdedx3->SetLineColor(kGreen + 3);
69  fdedx3->Draw("SAME");
70  fsum3->Draw("SAME");
71  fdiff3->Draw("SAME");
72 
73  TF1 *fdedx4 = new TF1("fdedx4", fundedx, 0., 1.5);
74  fdedx4->SetParameters(1.71069, 2.21786, 3.82597);
75  TF1 *fsum4 = new TF1("fsum4", funsum, 0., 1.5);
76  fsum4->SetParameters(1.71069, 2.21786, 3.82597, 0.260654, 0.385838, 0.307763);
77  TF1 *fdiff4 = new TF1("fdiff4", fundiff, 0., 1.5);
78  fdiff4->SetParameters(1.71069, 2.21786, 3.82597, 0.260654, 0.385838, 0.307763);
79  c->cd(4);
80  fdedx4->SetLineColor(kViolet + 2);
81  fdedx4->Draw("SAME");
82  fsum4->Draw("SAME");
83  fdiff4->Draw("SAME");
84 
85  TF1 *fdedx5 = new TF1("fdedx5", fundedx, 0., 1.5);
86  fdedx5->SetParameters(4.26561, -0.152697, 4.03156);
87  TF1 *fsum5 = new TF1("fsum5", funsum, 0., 1.5);
88  fsum5->SetParameters(4.26561, -0.152697, 4.03156, 0.338277, -0.664876, 0.656844);
89  TF1 *fdiff5 = new TF1("fdiff5", fundiff, 0., 1.5);
90  fdiff5->SetParameters(4.26561, -0.152697, 4.03156, 0.338277, -0.664876, 0.656844);
91 
92  c->cd(5);
93  fdedx5->SetLineColor(kGray);
94  fdedx5->Draw("SAME");
95  fsum5->SetLineColor(kYellow);
96  fsum5->Draw("SAME");
97  fdiff5->SetLineColor(kYellow);
98  fdiff5->Draw("SAME");
99 
100  c->cd(6);
101  fdedx1->Draw("SAME");
102  fdedx2->Draw("SAME");
103  fdedx3->Draw("SAME");
104  fdedx4->Draw("SAME");
105  fdedx5->Draw("SAME");
106 
107 }
TF1 * fdedx3
Definition: draw_bands.C:61
TF1 * fdiff4
Definition: draw_bands.C:77
TString funsum
Definition: dedx_bands.C:148
TH1F * h4
TF1 * fdiff1
Definition: draw_bands.C:41
TF1 * fsum2
Definition: draw_bands.C:51
TF1 * fdedx1
Definition: draw_bands.C:37
TH2F * h5
Definition: draw_bands.C:19
TString fundiff
Definition: dedx_bands.C:151
TF1 * fsum4
Definition: draw_bands.C:75
TString fundedx
Definition: dedx_bands.C:26
TF1 * fdiff5
Definition: draw_bands.C:89
TF1 * fdiff2
Definition: draw_bands.C:53
TF1 * fsum1
Definition: draw_bands.C:39
TF1 * fsum3
Definition: draw_bands.C:63
TF1 * fdedx2
Definition: draw_bands.C:49
TF1 * fdiff3
Definition: draw_bands.C:65
TH1F * h3
TF1 * fdedx5
Definition: draw_bands.C:85
TF1 * fsum5
Definition: draw_bands.C:87
TString funsig
Definition: dedx_bands.C:36
TF1 * fdedx4
Definition: draw_bands.C:73
TFile infile("dedx_out.root","READ")