]>
Commit | Line | Data |
---|---|---|
b6340d07 | 1 | void DrawTestAbso(Text_t *FileName = "MUONTestAbso") |
2 | { | |
3 | // FileName : name of the input file (.root) | |
4 | // Draw histograms to test the absorber correction (mean energy loss + Branson correction) | |
5 | // MUONTestAbso.root can be produced by the mean of the macro MUONTestAbso.C | |
6 | ||
7 | ||
8 | strcat(FileName,".root"); | |
9 | TFile *hfile = (TFile*)gROOT->FindObject(FileName); | |
10 | if (hfile) hfile->Close(); | |
11 | TFile *fill = new TFile(FileName); | |
12 | ||
13 | gROOT->Reset(); | |
14 | gStyle->SetOptStat(11111111); | |
15 | gStyle->SetOptFit(1); | |
16 | ||
17 | ||
18 | TCanvas *c1 = new TCanvas("c1",FileName, 250, 20, 800, 820); | |
19 | c1->GetFrame()->SetBorderSize(6); | |
20 | c1->GetFrame()->SetBorderMode(-1); | |
21 | ||
22 | c1->Divide(1,2); | |
23 | ||
24 | TPad *pad1 = new TPad("pad1", "The pad with the Histo", 0.01, 0.01, 0.99, 0.99, 21); | |
25 | c1->cd(1); | |
26 | pad1->Draw(); | |
27 | pad1->cd(); | |
28 | pad1->SetGridx(); | |
29 | pad1->SetGridy(); | |
30 | pad1->SetFillColor(22); | |
31 | pad1->GetFrame()->SetBorderMode(-1); | |
32 | hDeltaP1->Draw(); | |
33 | ||
34 | TPad *pad2 = new TPad("pad2", "The pad with the Histo", 0.01, 0.01, 0.99, 0.99, 21); | |
35 | c1->cd(2); | |
36 | pad2->Draw(); | |
37 | pad2->cd(); | |
38 | pad2->SetGridx(); | |
39 | pad2->SetGridy(); | |
40 | pad2->SetFillColor(22); | |
41 | pad2->GetFrame()->SetBorderMode(-1); | |
42 | hDeltaAngle1->Draw(); | |
43 | ||
44 | ||
45 | TCanvas *c2 = new TCanvas("c2",FileName, 200, 10, 800, 820); | |
46 | c2->GetFrame()->SetBorderSize(6); | |
47 | c2->GetFrame()->SetBorderMode(-1); | |
48 | ||
49 | c2->Divide(1,2); | |
50 | ||
51 | TPad *pad3 = new TPad("pad3", "The pad with the Histo", 0.01, 0.01, 0.99, 0.99, 21); | |
52 | c2->cd(1); | |
53 | pad3->Draw(); | |
54 | pad3->cd(); | |
55 | pad3->SetGridx(); | |
56 | pad3->SetGridy(); | |
57 | pad3->SetFillColor(22); | |
58 | pad3->GetFrame()->SetBorderMode(-1); | |
59 | hDeltaP2->Draw(); | |
60 | ||
61 | TPad *pad4 = new TPad("pad4", "The pad with the Histo", 0.01, 0.01, 0.99, 0.99, 21); | |
62 | c2->cd(2); | |
63 | pad4->Draw(); | |
64 | pad4->cd(); | |
65 | pad4->SetGridx(); | |
66 | pad4->SetGridy(); | |
67 | pad4->SetFillColor(22); | |
68 | pad4->GetFrame()->SetBorderMode(-1); | |
69 | hDeltaAngle2->Draw(); | |
70 | ||
71 | TCanvas *c3 = new TCanvas("c3",FileName, 150, 5, 800, 700); | |
72 | ||
73 | c3->GetFrame()->SetBorderSize(6); | |
74 | c3->GetFrame()->SetBorderMode(-1); | |
75 | ||
76 | TPad *pad5 = new TPad("pad5", "The pad with the Histo", 0.01, 0.01, 0.99, 0.99, 21); | |
77 | pad5->Draw(); | |
78 | pad5->cd(); | |
79 | pad5->SetGridx(); | |
80 | pad5->SetGridy(); | |
81 | pad5->SetFillColor(22); | |
82 | pad5->GetFrame()->SetBorderMode(-1); | |
83 | ||
84 | g1= new TF1("g1","gaus",9.3,9.8) ; // 9.25 | |
85 | hInvMassRes->Fit("g1","RQ"); | |
86 | hInvMassRes->GetXaxis()->SetTitleFont(20); | |
87 | hInvMassRes->SetXTitle("Mass (GeV/c^2!)"); | |
88 | hInvMassRes->SetFillColor(63); | |
89 | hInvMassRes->Draw(); | |
90 | ||
91 | } |