Transition PWG0 -> PWGUD

[u/mrichter/AliRoot.git] / PWGUD / selectors / trigger / show.C

diff --git a/PWGUD/selectors/trigger/show.C b/PWGUD/selectors/trigger/show.C
new file mode 100644 (file)
index 0000000..4a8a1d3
--- /dev/null
+++ b/PWGUD/selectors/trigger/show.C
@@ -0,0 +1,60 @@
+void showAll(const char* fileName = "trigger.root")
+{
+  TH1* fStats[100];
+
+  TFile::Open(fileName);
+  
+  Int_t count = 0;
+  while (1)
+  {
+    hist = (TH1*) gFile->Get(Form("fStats_%d", count));
+    if (!hist)
+      break;
+    
+    c = new TCanvas;
+    hist->Rebin(10);
+    hist->Draw();
+    c->SaveAs(Form("trigger_%d.png", count));
+    Printf("%s: %d", hist->GetTitle(), (Int_t) hist->GetEntries());
+    
+    fStats[count] = hist;
+    count++;
+  }
+
+  Int_t base = 1;
+  if (fStats[0])
+    base = (Int_t) fStats[0]->Integral(0, fStats[0]->GetNbinsX()+1);
+
+  for (Int_t i=0; i<count; i++)
+    if (fStats[i])
+    {
+      c->cd(i+1);
+      fStats[i]->Draw();
+      Printf("%s: %d triggers | %f %% of all triggered", fStats[i]->GetTitle(), (UInt_t) fStats[i]->Integral(0, fStats[i]->GetNbinsX()+1), 100.0 * fStats[i]->Integral(0, fStats[i]->GetNbinsX()+1) / base);
+    }
+
+}
+
+void show(const char* fileName = "trigger.root")
+{
+  TFile::Open(fileName);
+  
+  c = new TCanvas("c", "c", 1500, 1000);
+  c->Divide(4, 3);
+  
+  for (Int_t i=0; i<4; i++)
+  {
+    for (Int_t j=6; j<=8; j++)
+    {
+      hist = (TH1*) gFile->Get(Form("fStats_%d_%d", i, j));
+      c->cd((j-6)*4+i+1);
+      hist->Rebin(2);
+      hist->Rebin(2);
+      hist->Rebin(2);
+      hist->Draw();
+    }
+  }  
+      
+  c->SaveAs("trigger_rates.png");
+}
+