]>
Commit | Line | Data |
---|---|---|
91b112bd | 1 | DrawPi0Flow(const TString filename = "Pi0Flow_000167920.root") |
8e6cc59c | 2 | { |
b8bea077 | 3 | /* $Id$ */ |
4 | ||
8e6cc59c | 5 | gStyle->SetOptStat(0); |
91b112bd | 6 | TFile * f = new TFile(filename) ; |
8b4ce769 | 7 | //TList *histoList = (TList*)f->Get("PHOSPi0Flow"); |
8 | TList *histoList = (TList*)f->Get("PHOSPi0Flow/PHOSPi0FlowCoutput1"); // lego train | |
8e6cc59c | 9 | |
91b112bd | 10 | TH1F *hev = (TH1F*)histoList->FindObject("hTotSelEvents") ; |
11 | TH2F* hZvertex = (TH2F*)histoList->FindObject("hZvertex"); | |
12 | TH2F* hCenTrack = (TH2F*)histoList->FindObject("hCenTrack"); | |
13 | TH2F *hCenPHOS = (TH2F*)histoList->FindObject("hCenPHOS") ; | |
14 | TH2F* phiRPV0A = (TH2F*)histoList->FindObject("phiRPV0A"); | |
15 | TH2F* phiRPV0C = (TH2F*)histoList->FindObject("phiRPV0C"); | |
16 | TH2F* hCellNXZM1 = (TH2F*)histoList->FindObject("hCellNXZM1"); | |
17 | TH2F* hCellNXZM2 = (TH2F*)histoList->FindObject("hCellNXZM2"); | |
18 | TH2F* hCellNXZM3 = (TH2F*)histoList->FindObject("hCellNXZM3"); | |
19 | TH2F *hPi0All_cen0 = (TH2F*)histoList->FindObject("hPi0All_cen0"); | |
8e6cc59c | 20 | TH2F *hPi0Allcore_cen0 = (TH2F*)histoList->FindObject("hPi0Allcore_cen0"); |
21 | ||
91b112bd | 22 | //----------------------------------------------------------------------------- |
8e6cc59c | 23 | TCanvas *c1 = new TCanvas("c1","Event selection"); |
8b4ce769 | 24 | hev->GetXaxis()->SetRangeUser(0,5); |
91b112bd | 25 | hev->GetXaxis()->SetBinLabel(1,"Total"); |
8b4ce769 | 26 | hev->GetXaxis()->SetBinLabel(2,"Has Vertex"); |
27 | hev->GetXaxis()->SetBinLabel(3,"abs(z_vertex) < 10."); | |
28 | hev->GetXaxis()->SetBinLabel(4,"Has Centrality"); | |
29 | hev->GetXaxis()->SetBinLabel(5,"C. upper edge"); | |
30 | hev->GetXaxis()->SetBinLabel(6,"C. lower edge"); | |
8e6cc59c | 31 | hev->SetYTitle("N_{events}"); |
8b4ce769 | 32 | hev->GetYaxis()->SetTitleOffset(1.2); |
8e6cc59c | 33 | hev->Draw(); |
34 | c1->Print("PHOS_EvSel.eps"); | |
35 | ||
91b112bd | 36 | //----------------------------------------------------------------------------- |
37 | TCanvas *c2 = new TCanvas("c2","Track multiplicity vs centrality"); | |
38 | hCenTrack->SetXTitle("centrality (%)"); | |
39 | hCenTrack->SetYTitle("Number of tracks"); | |
40 | hCenTrack->Draw("colz"); | |
41 | c2->Print("TrackMultCentrality.eps"); | |
42 | ||
43 | //----------------------------------------------------------------------------- | |
44 | TCanvas *c3 = new TCanvas("c3","PHOS multiplicity vs centrality"); | |
8e6cc59c | 45 | hCenPHOS->SetXTitle("centrality (%)"); |
46 | hCenPHOS->SetYTitle("Number of PHOS clusters"); | |
47 | hCenPHOS->Draw("colz"); | |
91b112bd | 48 | c3->Print("PHOSMultCentrality.eps"); |
49 | ||
50 | //----------------------------------------------------------------------------- | |
51 | TCanvas *c4 = new TCanvas("c4","RP phi"); | |
52 | phiRPV0A1 = phiRPV0A->ProjectionX(); | |
53 | phiRPV0C1 = phiRPV0C->ProjectionX(); | |
54 | phiRPV0A1->SetLineColor(kRed); | |
55 | phiRPV0C1->SetLineColor(kBlue); | |
56 | phiRPV0A1->SetXTitle("RP #phi"); | |
57 | phiRPV0C1->SetXTitle("RP #phi"); | |
58 | phiRPV0A1->Draw(); | |
59 | phiRPV0C1->Draw("same"); | |
60 | leg = new TLegend(0.7,0.8,0.89,0.89); | |
61 | leg->SetFillColor(kWhite); | |
62 | leg->SetBorderSize(1); | |
63 | leg->AddEntry(phiRPV0A1,"V0A","l"); | |
64 | leg->AddEntry(phiRPV0C1,"V0C","l"); | |
65 | leg->Draw(); | |
66 | c4->Print("V0RPphi.eps"); | |
67 | ||
68 | //----------------------------------------------------------------------------- | |
69 | TCanvas *c5 = new TCanvas("c5","XZ in M1"); | |
70 | c5->SetLogz(); | |
71 | hCellNXZM1->SetTitle("Cell occupancy in module 4"); | |
72 | hCellNXZM1->SetXTitle("X (cells)"); | |
73 | hCellNXZM1->SetYTitle("X (cells)"); | |
74 | hCellNXZM1->Draw("colz"); | |
75 | c5->Print("PHOS_XYM4.eps"); | |
8e6cc59c | 76 | |
91b112bd | 77 | //----------------------------------------------------------------------------- |
78 | TCanvas *c6 = new TCanvas("c6","XZ in M2"); | |
79 | c6->SetLogz(); | |
80 | hCellNXZM2->SetTitle("Cell occupancy in module 3"); | |
81 | hCellNXZM2->SetXTitle("X (cells)"); | |
82 | hCellNXZM2->SetYTitle("X (cells)"); | |
83 | hCellNXZM2->Draw("colz"); | |
84 | c6->Print("PHOS_XYM3.eps"); | |
85 | ||
86 | //----------------------------------------------------------------------------- | |
87 | TCanvas *c7 = new TCanvas("c7","XZ in M3"); | |
88 | c7->SetLogz(); | |
89 | hCellNXZM3->SetTitle("Cell occupancy in module 2"); | |
90 | hCellNXZM3->SetXTitle("X (cells)"); | |
91 | hCellNXZM3->SetYTitle("X (cells)"); | |
92 | hCellNXZM3->Draw("colz"); | |
93 | c7->Print("PHOS_XYM2.eps"); | |
94 | ||
95 | //----------------------------------------------------------------------------- | |
96 | TCanvas *c8 = new TCanvas("c8","gg mass vs pt, PID=All"); | |
97 | c8->SetLogz(); | |
8e6cc59c | 98 | hPi0All_cen0->SetTitle("PID: All"); |
99 | hPi0All_cen0->SetXTitle("M_{#gamma#gamma} (GeV/c^{2})"); | |
100 | hPi0All_cen0->SetYTitle("p_{T} (GeV/c)"); | |
101 | hPi0All_cen0->Draw("colz"); | |
102 | c3->Print("PHOS_MggAll.eps"); | |
103 | ||
91b112bd | 104 | //----------------------------------------------------------------------------- |
105 | TCanvas *c9 = new TCanvas("c9","gg mass vs pt, PID=Allcore"); | |
106 | c9->SetLogz(); | |
8e6cc59c | 107 | hPi0Allcore_cen0->SetTitle("PID: Allcore"); |
108 | hPi0Allcore_cen0->SetXTitle("M_{#gamma#gamma} (GeV/c^{2})"); | |
109 | hPi0Allcore_cen0->SetYTitle("p_{T} (GeV/c)"); | |
110 | hPi0Allcore_cen0->Draw("colz"); | |
91b112bd | 111 | c9->Print("PHOS_MggAllcore.eps"); |
8e6cc59c | 112 | |
113 | } |