EVE/alice-macros/jetplane.C

   1 // $Id$
   2 // Main authors: Matevz Tadel & Alja Mrak-Tadel: 2006, 2007
   3
   4 /**************************************************************************
   5  * Copyright(c) 1998-2008, ALICE Experiment at CERN, all rights reserved. *
   6  * See http://aliceinfo.cern.ch/Offline/AliRoot/License.html for          *
   7  * full copyright notice.                                                 *
   8  **************************************************************************/
   9
  10 TEveViewer *gJPView   = 0;
  11 TEveScene  *gJPScene  = 0;
  12
  13 AliEveJetPlane* jetplane()
  14 {
  15   if (gJPView == 0)
  16   {
  17     TEveWindowSlot *slot    = 0;
  18     TEveBrowser    *browser = gEve->GetBrowser();
  19
  20     slot = TEveWindow::CreateWindowInTab(browser->GetTabRight());
  21     slot->MakeCurrent();
  22     gJPView  = gEve->SpawnNewViewer("JetPlane", "");
  23     gJPScene = gEve->SpawnNewScene("JetPlane", "Scene holding elements of the jet-plane view.");
  24     gJPView->AddScene(gJPScene);
  25
  26     gJPView->GetGLViewer()->SetCurrentCamera(TGLViewer::kCameraPerspXOY);
  27   }
  28
  29   AliAODEvent* aod = AliEveEventManager::AssertAOD();
  30
  31
  32   // We have event id everywhere now.
  33   Int_t iev = AliEveEventManager::GetMaster()->GetEventId();
  34
  35   AliEveJetPlane* jp = new AliEveJetPlane(iev);
  36   jp->SetPickable(kTRUE);
  37
  38   // Read Jets in current event
  39
  40   TClonesArray* jets = aod->GetJets();
  41   Int_t njets = jets->GetEntries();
  42   printf("Event: %5d Number of jets: %5d \n", iev, njets);
  43
  44   for (Int_t ij = 0; ij < njets; ij++)
  45   {
  46     AliAODJet jet = (AliAODJet) jets->At(ij);
  47     jp->AddJet(jet);
  48   }
  49
  50   // Read tracks in current event
  51
  52   TClonesArray* tracks = aod->GetTracks();
  53   Int_t ntracks = tracks->GetEntries();
  54   printf("Event: %5d Number of tracks: %5d \n", iev, ntracks);
  55
  56   for (Int_t ij = 0; ij < ntracks; ij++)
  57   {
  58     AliAODTrack track = (AliAODTrack) tracks->At(ij);
  59     jp->AddTrack(track);
  60   }
  61
  62   // Render Jet Plane
  63   gStyle->SetPalette(1, 0);
  64   // gEve->AddElement(jp);
  65   gJPScene->AddElement(jp);
  66
  67   gEve->Redraw3D();
  68
  69   return jp;
  70 }