Modifications needed by the HBT analysis (P.Skowronski)
[u/mrichter/AliRoot.git] / MUON / DrawTestAbso.C
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 }