Various updates for the production scripts

[u/mrichter/AliRoot.git] / PWGLF / FORWARD / analysis2 / sim / PlotSysInfo.C

diff --git a/PWGLF/FORWARD/analysis2/sim/PlotSysInfo.C b/PWGLF/FORWARD/analysis2/sim/PlotSysInfo.C
new file mode 100644 (file)
index 0000000..99ccd98
--- /dev/null
+++ b/PWGLF/FORWARD/analysis2/sim/PlotSysInfo.C
@@ -0,0 +1,182 @@
+TVirtualPad* MakeCanvas(const char* title="")
+{
+  static Int_t cId = 0;
+  TCanvas* c = new TCanvas(Form("c%02d", ++cId), title);
+  c->SetTopMargin(0.02);
+  c->SetRightMargin(0.02);
+  c->cd();
+  return c;
+}
+
+Double_t SumUsage(TTree*      tree, 
+                 const char* exp, 
+                 const char* cut, 
+                 bool        draw=false)
+{
+  //
+  // return sum of usage
+  //
+  if (draw) MakeCanvas(Form("%s [%s]", exp, cut));
+  
+  Int_t  entries = tree->Draw(exp, cut, (draw ? "" : "goff"));
+  if (entries==0) return 0;
+
+  Double_t mean = TMath::Mean(entries, tree->GetV1());
+  return entries * mean;
+}
+Double_t TopUsage(TTree*      tree, 
+                 const char* exp, 
+                 const char* cut, 
+                 Int_t       order)
+{
+  Int_t entries = tree->Draw(exp, cut, "goff");
+  if (entries <= 1 || !tree->GetV1()) return -10000;
+  
+  TArrayI index(entries);
+  TMath::Sort(entries, tree->GetV1(), index.fArray);
+
+  Int_t    oindex = TMath::Min(order, entries);
+  Double_t value  = tree->GetV1()[index[oindex-1]];
+  
+  return value;
+}
+  
+TH1* ExtractHist(TTree*      tree,
+                const char* exp, 
+                const char* cut,
+                const char* name,
+                const char* xtitle="", 
+                const char* ytitle="",
+                Bool_t      draw=false)
+{
+  tree->Draw(Form("%s>>tmpa", exp), cut, "GOFF");
+  if (!tree->GetHistogram()) return 0;
+
+  TH1* ret = static_cast<TH1*>(tree->GetHistogram()->Clone(name));
+  delete tree->GetHistogram();
+  ret->SetXTitle(xtitle);
+  ret->SetYTitle(ytitle);
+  ret->SetMarkerStyle(22);
+  ret->SetMarkerSize(1);
+  if (draw) {
+    if (draw) MakeCanvas(name);
+    ret->Draw();
+  }
+  return ret;
+}
+                
+                
+void Print(std::ostream& o,
+          const char*   name, 
+          Double_t      dT, 
+          Double_t      dVM, 
+          Double_t      alldT, 
+          Double_t      alldVM)
+{
+  o << name << "\t" 
+    << dT   << "\t" 
+    << dVM  << "\t" 
+    << 100*(alldT  > 0 ? dT  / alldT  : 0) << "\t" 
+    << 100*(alldVM > 0 ? dVM / alldVM : 0) << std::endl;
+}
+
+const char* dets[] = {"ITS", 
+                     "TPC", 
+                     "TRD", 
+                     "TOF", 
+                     "PHOS", 
+                     "HMPID", 
+                     "EMCAL", 
+                     "MUON", 
+                     "FMD", 
+                     "ZDC", 
+                     "PMD", 
+                     "T0", 
+                     "VZERO", 
+                     "ACORDE", 
+                     "HLT",
+                     0 };
+
+
+void
+Plot1SysInfo(const char* file, UShort_t draw=0x1)
+{
+  // gROOT->LoadMacro("$ALICE_ROOT/../master-src/macros/PlotSys.C+");
+
+  // --- Create output file and tree ---------------------------------
+  TString rootOut(file); 
+  rootOut.ReplaceAll(".log", ".root");
+  rootOut.ReplaceAll(".foo", ".root");
+  Info("", "Writing to ROOT file %s", rootOut.Data());
+  TFile* out  = TFile::Open(rootOut, "RECREATE");
+  TTree* tree = AliSysInfo::MakeTree(file);
+
+
+  
+  // --- Create ASCII output ------------------------------------------
+  TString sumOut(rootOut); 
+  sumOut.ReplaceAll(".root", ".sum");
+  Info("", "Writing to ASCII file %s", sumOut.Data());
+  std::ofstream ascii(sumOut.Data());
+  ascii << "Det/C:sumDt/F:sumDvm/F:fracDt/F:fracDvm/F" << std::endl;
+
+  // --- Get global stuff ---------------------------------------------
+  const char* all     = "id0>=0&&id2>=0";
+  Double_t sumdTAll   = SumUsage(tree, "deltaT",  all, draw & 0x1);
+  Double_t sumdVMAll  = SumUsage(tree, "deltaVM", all, draw & 0x1);
+  Double_t topdT      = TopUsage(tree, "deltaT",  "id2<3", 20);
+  Double_t topdVM     = TopUsage(tree, "deltaVM", "", 20);
+  TCut     cutT("cutDT", Form("deltaT > %f", topdT));
+  TCut     cutVM("cutVM", Form("deltaVM > %f", topdVM));
+  
+  ExtractHist(tree, "deltaVM:sname", "1"+cutVM,
+             "DVMvsName","","#DeltaVM [MB]", draw&0x4);
+  ExtractHist(tree, "VM:sname", "id2<3"+cutVM, 
+             "VMvsName", "", "VM [MB]", draw&0x4);
+  ExtractHist(tree, "VM:T", "deltaVM>1", 
+             "VMvsTime", "Time [sec]", "VM [MB]", draw&0x4);
+  ExtractHist(tree, "deltaT:sname","id2<3"+cutT,
+             "CPUvsName","","#DeltaT [sec]", draw&0x4);
+  
+
+
+  Print(ascii, "all", sumdTAll, sumdVMAll, sumdTAll, sumdVMAll);
+
+  
+
+  // --- Loop over detetors ------------------------------------------
+  const char** pdet = dets;
+  Int_t        idet = 0;
+  while (*pdet) { 
+    TString  cut    = Form("id0==%d && id2 >= 0", idet);
+    Double_t sumdT  = SumUsage(tree, "deltaT",  cut, draw & 0x2);
+    Double_t sumdVM = SumUsage(tree, "deltaVM", cut, draw & 0x2);
+    Print(ascii, *pdet, sumdT, sumdVM, sumdTAll, sumdVMAll);
+    
+#if 0
+    TString cut2    = Form("id0==%d",idet);
+    ExtractHist(tree, "deltaVM:sname", cut2.Data()+cutVM,
+               Form("DVMvsName_%02d", idet), "", "#DeltaVM [MB]", draw&0x8);
+    ExtractHist(tree, "VM:sname",      cut2.Data()+cutVM,
+               Form("VMvsName_%02d", idet), "", "VM [MB]", draw&0x8);
+    ExtractHist(tree, "deltaT:sname",  cut2.Data()+cutT, 
+               Form("CPUvsName_%02d", idet),"", "#DeltaT [sec]", draw&0x8);
+#endif
+
+    pdet++;
+    idet++;
+  }
+  ascii.close();
+
+  new TBrowser;
+}
+
+  
+
+PlotSysInfo(ULong_t pid=431808952)
+{
+  Plot1SysInfo(Form("%d_simwatch.log", pid));
+  Plot1SysInfo(Form("%d_recowatch.log", pid));
+}
+
+