tpc_clusters.C - added missing AliRunLoader::UnloadRecPoints(). New tpc_raw.C macro...
[u/mrichter/AliRoot.git] / EVE / alice-macros / acorde_raw.C
CommitLineData
c2c4b7a2 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
10TString acorde_module_path(Int_t module);
11
12
13void acorde_raw()
14{
6226bf2b 15 AliEveEventManager::AssertGeometry();
647814a2 16
c2c4b7a2 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
c2c4b7a2 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
647814a2 52 // @@NEWROOT@@ Temporary hack.
c2c4b7a2 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
78TString acorde_module_path(Int_t module)
79{
c2c4b7a2 80 if (module < 0 || module > 59)
81 {
82 Error("acorde_module_path", "module %d out of range.", module);
83 return "";
84 }
85
d3c0a260 86 TGeoPNEntry* pne = gGeoManager->GetAlignableEntry(Form("ACORDE/Array%d", module + 1));
87 if(!pne) return "missing_pne";
c2c4b7a2 88
d3c0a260 89 return Form("%s/ACORDE2_5", pne->GetTitle());
c2c4b7a2 90}