tpc_clusters.C - added missing AliRunLoader::UnloadRecPoints(). New tpc_raw.C macro...
[u/mrichter/AliRoot.git] / 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 }