EVE/alice-macros/acorde_raw.C

   1 // $Id: geom_acorde.C 23412 2008-01-18 21:04:54Z mtadel $
   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 TString acorde_module_path(Int_t module);
  11
  12
  13 void acorde_raw()
  14 {
  15   AliEveEventManager::AssertGeometry();
  16
  17   AliRawReader       * reader = AliEveEventManager::AssertRawReader();
  18   AliACORDERawStream * stream = new AliACORDERawStream(reader);
  19
  20   stream->Reset();
  21   stream->Next();
  22
  23   UInt_t dy[4];
  24   dy[0] = stream->GetWord(0);
  25   dy[1] = stream->GetWord(1);
  26   dy[2] = stream->GetWord(2);
  27   dy[3] = stream->GetWord(3);
  28
  29   printf ("ACORDE event 0x%08x 0x%08x 0x%08x 0x%08x\n", dy[0], dy[1], dy[2], dy[3]);
  30
  31   TEveElementList* acorde = new TEveElementList("ACORDE Raw");
  32
  33   gEve->AddElement(acorde);
  34
  35   Int_t shape_offset = TEveGeoShape::Class()->GetDataMemberOffset("fShape");
  36
  37   for (Int_t module=0; module < 60; ++module)
  38   {
  39     TString path = acorde_module_path(module);
  40     // printf("%2d - %s\n", i, path.Data());
  41
  42     if ( ! gGeoManager->cd(path))
  43     {
  44       Warning("acorde_raw", "Module id=%d, path='%s' not found.\n", module, path.Data());
  45       continue;
  46     }
  47
  48     TEveGeoShape* eg_shape = new TEveGeoShape(Form("Module %d", module));
  49
  50     eg_shape->RefMainTrans().SetFrom(* gGeoManager->GetCurrentMatrix());
  51
  52     // @@NEWROOT@@ Temporary hack.
  53     // Hack to set shape pointer, no interface available in TEveGeoShape.
  54     * (TGeoShape**) (((char*)eg_shape) + shape_offset) = gGeoManager->GetCurrentVolume()->GetShape();
  55
  56     // From Matevz:
  57     // Here check ctate and assign color, I do it partially for now.
  58     Int_t  word_idx = module / 30;
  59     Int_t  bit_idx  = module % 30;
  60     Bool_t val      = (dy[word_idx] & (1 << bit_idx)) != 0;
  61     //printf("Module %2d: word_idx = %d, bit_idx = %2d => val = %d\n",
  62     //       module, word_idx, bit_idx, val);
  63     if (val)
  64       eg_shape->SetMainColor((Color_t) 2);
  65     else
  66       eg_shape->SetMainColor((Color_t) 4);
  67     eg_shape->StampColorSelection();
  68
  69     acorde->AddElement(eg_shape);
  70   }
  71
  72   gEve->Redraw3D();
  73 }
  74
  75 //==============================================================================
  76 //==============================================================================
  77
  78 TString acorde_module_path(Int_t module)
  79 {
  80   if (module < 0 || module > 59)
  81   {
  82     Error("acorde_module_path", "module %d out of range.", module);
  83     return "";
  84   }
  85
  86   TGeoPNEntry* pne = gGeoManager->GetAlignableEntry(Form("ACORDE/Array%d", module + 1));
  87   if(!pne) return "missing_pne";
  88
  89   return Form("%s/ACORDE2_5", pne->GetTitle());
  90 }