plot_outliers.C
Go to the documentation of this file.
1 void Cut(double x)
2 {
3  TGraph* g = new TGraph;
4  g->SetPoint(0, x, 0);
5  g->SetPoint(1, x, 1e6);
6  g->SetLineColor(kRed);
7  g->SetLineStyle(7);
8  g->SetLineWidth(2);
9  g->Draw("l same");
10 }
11 
12 void Grid()
13 {
14  for(int b = 1; b < 22; ++b){
15  TLine* l = new TLine(32*b, 0, 32*b, 32*12);
16  l->SetLineColor(kGray);
17  l->Draw();
18  }
19 
20  for(int m = 1; m < 12; ++m){
21  TLine* l = new TLine(0, 32*m, 64*11, 32*m);
22  l->SetLineColor(kGray);
23  l->Draw();
24  }
25 }
26 
28 {
29  gROOT->ForceStyle();
30 
31  TFile* f = new TFile("outliers_hist.root");
32 
33  TH1* hq = f->Get("find/qual");
34  hq->SetTitle("Far detector data");
35  hq->Draw();
36 
37  Cut(0.2);
38 
39  gPad->Print("plots/quality.eps");
40  gPad->SetLogy();
41  gPad->Print("plots/quality_log.eps");
42 
43 
44  new TCanvas;
45  TH1* hp = f->Get("find/pe");
46  hp->SetTitle("Far detector data");
47  hp->Draw();
48 
49  Cut(4.75);
50  Cut(11);
51 
52  gPad->Print("plots/pe.eps");
53  gPad->SetLogy();
54  gPad->Print("plots/pe_log.eps");
55 
56 
57  new TCanvas;
58  TH2* hex = f->Get("find/emptyX");
59  hex->SetTitle("Empty cells - X view");
60  hex->Draw("col");
61  Grid();
62  gPad->Print("plots/map_empty_x.eps");
63 
64  new TCanvas;
65  TH2* hey = f->Get("find/emptyY");
66  hey->SetTitle("Empty cells - Y view");
67  hey->Draw("col");
68  Grid();
69  gPad->Print("plots/map_empty_y.eps");
70 
71  new TCanvas;
72  TH2* hbx = f->Get("find/badfitX");
73  hbx->SetTitle("Bad fits - X view");
74  hbx->Draw("col");
75  Grid();
76  gPad->Print("plots/map_bad_fit_x.eps");
77 
78  new TCanvas;
79  TH2* hby = f->Get("find/badfitY");
80  hby->SetTitle("Bad fits - Y view");
81  hby->Draw("col");
82  Grid();
83  gPad->Print("plots/map_bad_fit_y.eps");
84 
85  new TCanvas;
86  TH2* hlx = f->Get("find/lowpeX");
87  hlx->SetTitle("Low PE cells - X view");
88  hlx->Draw("col");
89  Grid();
90  gPad->Print("plots/map_low_pe_x.eps");
91 
92  new TCanvas;
93  TH2* hly = f->Get("find/lowpeY");
94  hly->SetTitle("Low PE cells - Y view");
95  hly->Draw("col");
96  Grid();
97  gPad->Print("plots/map_low_pe_y.eps");
98 
99  new TCanvas;
100  TH2* hhx = f->Get("find/highpeX");
101  hhx->SetTitle("High PE cells - X view");
102  hhx->Draw("col");
103  Grid();
104  gPad->Print("plots/map_high_pe_x.eps");
105 
106  new TCanvas;
107  TH2* hhy = f->Get("find/highpeY");
108  hhy->SetTitle("High PE cells - Y view");
109  hhy->Draw("col");
110  Grid();
111  gPad->Print("plots/map_high_pe_y.eps");
112 }
enum BeamMode kRed
void Cut(double x)
Definition: plot_outliers.C:1
void plot_outliers()
Definition: plot_outliers.C:27
const hit & b
Definition: hits.cxx:21
void Grid()
Definition: plot_outliers.C:12