macros/lego.C

   1 #if !defined(__CINT__) || defined(__MAKECINT__)
   2 #include <TROOT.h>
   3 #include <TFile.h>
   4 #include <TCanvas.h>
   5 #include <TH2F.h>
   6 #endif
   7
   8 void lego() {
   9 // macro to visualize the lego plots generated by gAlive->RunLego
  10
  11    gROOT->Reset();
  12    TFile *file = new TFile("galice.root");
  13
  14    Float_t theta = 10;
  15    Float_t phi   = 170;
  16    Int_t ncont   = 50;
  17
  18    TCanvas *cgcm2 = new TCanvas("cgcm2","gcm2",200,100,600,400);
  19    cgcm2->SetTheta(theta);
  20    cgcm2->SetPhi(phi);
  21    TH2F *hgcm2 = (TH2F*)file->Get("hgcm2");
  22    hgcm2->SetFillColor(2);
  23    hgcm2->SetMaximum(1);
  24    hgcm2->SetContour(ncont);
  25    hgcm2->SetMaximum(50);
  26    hgcm2->Draw("lego2sphe");
  27
  28    TCanvas *cabso = new TCanvas("cabso","abso",100,50,600,400);
  29    cabso->SetTheta(theta);
  30    cabso->SetPhi(phi);
  31    TH2F *habso = (TH2F*)file->Get("habso");
  32    habso->SetFillColor(2);
  33    habso->SetMaximum(1);
  34    habso->SetContour(ncont);
  35    habso->SetMaximum(1);
  36    habso->Draw("lego2sphe");
  37
  38    TCanvas *cradl = new TCanvas("cradl","radl",10,10,600,400);
  39    cradl->SetTheta(theta);
  40    cradl->SetPhi(phi);
  41    TH2F *hradl = (TH2F*)file->Get("hradl");
  42    hradl->SetFillColor(2);
  43    hradl->SetMaximum(1);
  44    hradl->SetContour(ncont);
  45    hradl->SetMaximum(5);
  46    hradl->Draw("lego2sphe");
  47 }