EVE/alice-macros/emcal_hits.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 TEvePointSet*
  11 emcal_hits(const char *varexp    = "fX:fY:fZ",
  12            const char *selection = "",
  13            TEveElement* cont = 0)
  14 {
  15   AliRunLoader* rl =  AliEveEventManager::AssertRunLoader();
  16   rl->LoadHits("EMCAL");
  17
  18   TTree* ht = rl->GetTreeH("EMCAL", false);
  19
  20   TEvePointSet* points = new TEvePointSet(Form("EMCAL Hits '%s'", selection));
  21
  22   TEvePointSelector ps(ht, points, varexp, selection);
  23   ps.Select();
  24
  25   if (points->Size() == 0 && gEve->GetKeepEmptyCont() == kFALSE) {
  26     Warning("emcal_hits", Form("No hits match '%s'", selection));
  27     delete points;
  28     return 0;
  29   }
  30
  31   points->SetName(Form("EMCAL Hits"));
  32   const TString viz_tag("SIM Hits EMCAL");
  33   points->ApplyVizTag(viz_tag, "Hits");
  34
  35   points->SetTitle(Form("N=%d", points->Size()));
  36   points->SetMarkerSize(.5);
  37   points->SetMarkerColor(2);
  38
  39   gEve->AddElement(points, cont);
  40   gEve->Redraw3D();
  41
  42   return points;
  43 }