2 // Main authors: Matevz Tadel & Alja Mrak-Tadel: 2006, 2007
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 **************************************************************************/
11 #include <TEveManager.h>
12 #include <TEveQuadSet.h>
15 #include <TGeoManager.h>
17 #include <TEveTrans.h>
18 #include <TClonesArray.h>
21 #include <EveBase/AliEveEventManager.h>
23 #include <AliRunLoader.h>
24 #include <AliCluster.h>
25 #include <EMCAL/AliEMCALGeometry.h>
26 #include <EMCAL/AliEMCALDigit.h>
29 // #include <Riostream.h>
34 AliEveEventManager::AssertGeometry();
36 TGeoNode* node = gGeoManager->GetTopVolume()->FindNode("XEN1_1");
39 Int_t nModules = node->GetNdaughters();
41 TEveElementList* l = new TEveElementList("EMCAL");
42 l->SetTitle("Tooltip");
45 TGeoBBox* bbbox = (TGeoBBox*) node->GetDaughter(0) ->GetVolume()->GetShape();
46 TEveFrameBox* frame_big = new TEveFrameBox();
47 frame_big->SetFrameColorRGBA(200,200,0,50);
48 frame_big->SetAABoxCenterHalfSize(0, 0, 0, bbbox->GetDX(), bbbox->GetDY(), bbbox->GetDZ());
50 TEveFrameBox* frame_sml = 0x0;
53 TGeoBBox* sbbox = (TGeoBBox*) node->GetDaughter(10)->GetVolume()->GetShape();
54 frame_sml = new TEveFrameBox();
55 frame_sml->SetFrameColorRGBA(200,200,0,50);
56 frame_sml->SetAABoxCenterHalfSize(0, 0, 0, sbbox->GetDX(), sbbox->GetDY(), sbbox->GetDZ());
59 gStyle->SetPalette(1, 0);
60 TEveRGBAPalette* pal = new TEveRGBAPalette(0, 512);
61 pal->SetLimits(0, 1024);
63 TEveQuadSet* smodules[12];
64 memset(smodules,0,12*sizeof(TEveQuadSet*));
67 AliEMCALGeometry * geom = AliEMCALGeometry::GetInstance();
68 if (!geom) geom = AliEMCALGeometry::GetInstance("","");
70 for (Int_t sm=0; sm<nModules; ++sm)
72 TEveQuadSet* q = new TEveQuadSet(Form("SM %d", sm+1));
74 q->Reset(TEveQuadSet::kQT_RectangleYZFixedDimX, kFALSE, 32);
75 q->SetDefWidth (geom->GetPhiTileSize());
76 q->SetDefHeight(geom->GetEtaTileSize());
78 q->RefMainTrans().SetFrom(*node->GetDaughter(sm)->GetMatrix());
80 q->SetFrame(sm < 10 ? frame_big : frame_sml);
83 gEve->AddElement(q, l);
87 AliRunLoader* rl = AliEveEventManager::AssertRunLoader();
89 rl->LoadDigits("EMCAL");
90 TTree* dt = rl->GetTreeD("EMCAL", kFALSE);
93 TClonesArray *digits = 0;
94 dt->SetBranchAddress("EMCAL", &digits);
96 Int_t nEnt = digits->GetEntriesFast();
110 for (Int_t idig = 0; idig < nEnt; ++idig)
112 dig = static_cast<AliEMCALDigit *>(digits->At(idig));
115 id = dig->GetId() ; //cell (digit) label
116 amp = dig->GetAmp(); //amplitude in cell (digit)
117 time = dig->GetTime();//time of creation of digit after collision
119 AliDebugGeneral("emcal_digits", 5, Form("Cell ID %3d, Amplitude: %f", id, amp));
120 // cout<<"Cell ID "<<id<<" Amp "<<amp<<endl;//" time "<<time<<endl;
123 geom->GetCellIndex(id,iSupMod,iTower,iIphi,iIeta);
124 //Gives SuperModule and Tower numbers
125 geom->GetCellPhiEtaIndexInSModule(iSupMod,iTower,
126 iIphi, iIeta,iphi,ieta);
127 //Gives label of cell in eta-phi position per each supermodule
129 AliDebugGeneral("emcal_digits", 5, Form("SModule %3d; Tover %3d; Eta %3d; Phi %3d; Cell Eta %3d; Cell Phi %3d",
130 iSupMod, iTower, iIeta, iIphi, ieta, iphi));
131 // cout<< "SModule "<<iSupMod<<"; Tower "<<iTower
132 // <<"; Eta "<<iIeta<<"; Phi "<<iIphi
133 // <<"; Cell Eta "<<ieta<<"; Cell Phi "<<iphi<<endl;
135 geom->RelPosCellInSModule(id, x, y, z);
136 // cout << x <<" "<< y <<" "<< z <<endl;
137 AliDebugGeneral("emcal_digits", 5, Form("(x,y,z)=(%8.3f,%8.3f,%8.3f)", x, y, z));
139 TEveQuadSet* q = smodules[iSupMod];
142 q->QuadValue(TMath::Nint(amp));
143 q->QuadId(new AliEMCALDigit(*dig));
146 AliDebugGeneral("emcal_digits", 1, Form("Digit pointer 0x0"));
147 // cout<<"Digit pointer 0x0"<<endl;
151 rl->UnloadDigits("EMCAL");
154 for (Int_t sm = 0; sm < nModules; ++sm)
156 smodules[iSupMod]->RefitPlex();