3 * @author Christian Holm Christensen <cholm@nbi.dk>
4 * @date Thu Nov 17 11:20:22 2011
6 * @brief Draw the cuts used in the analysis
8 * @ingroup pwg2_forward_scripts_qa
12 * Draw cuts used in analysis
14 * @param filename Input file name
16 * @ingroup pwg2_forward_scripts_qa
19 DrawCuts(const char* filename="forward.root")
21 gStyle->SetPalette(1);
23 gStyle->SetOptStat(0);
24 gStyle->SetTitleW(.4);
25 gStyle->SetTitleH(.1);
26 gStyle->SetTitleColor(0);
27 gStyle->SetTitleStyle(0);
28 gStyle->SetTitleBorderSize(0);
29 gStyle->SetTitleX(.6);
31 TFile* file = TFile::Open(filename, "READ");
33 Error("DrawCuts", "failed to open %s", filename);
37 TList* forward = static_cast<TList*>(file->Get("Forward"));
39 Error("DrawCuts", "List Forward not found in %s", filename);
43 TList* dc = static_cast<TList*>(forward->FindObject("fmdDensityCalculator"));
45 Error("DrawCuts", "List fmdDensityCalculator not found in Forward");
48 TList* sf = static_cast<TList*>(forward->FindObject("fmdSharingFilter"));
50 Error("DrawCuts", "List fmdSharingFilter not found in Forward");
53 TList* hc = static_cast<TList*>(forward->FindObject("fmdHistCollector"));
55 Error("DrawCuts", "List fmdHistCollector not found in Forward");
58 TH2* hC = static_cast<TH2*>(sf->FindObject("highCuts"));
60 Error("DrawCuts", "Histogram highCuts found in %s", sf->GetName());
63 TH2* lC = static_cast<TH2*>(dc->FindObject("lowCuts"));
65 Error("DrawCuts", "Histogram lowCuts found in %s", dc->GetName());
68 TH2* co = static_cast<TH2*>(hc->FindObject("coverage"));
70 Error("DrawCuts", "Histogram coverage found in %s", hc->GetName());
73 TCanvas* c = new TCanvas("cuts", "Cuts used in the analysis", 900, 700);
78 c->cd(1); hC->Draw("colz");
79 c->cd(2); lC->Draw("colz");
80 c->cd(3); co->Draw("colz");