macros/lego.C

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