]> git.uio.no Git - u/mrichter/AliRoot.git/blame - EVE/alice-macros/vzero_digits.C
Adding more information to the debug output
[u/mrichter/AliRoot.git] / EVE / alice-macros / vzero_digits.C
CommitLineData
d810d0de 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 *
51346b82 7 * full copyright notice. *
d810d0de 8 **************************************************************************/
63192539 9
16718cdc 10// fV0CHeight1 = 2.5; // height of cell 1
11// fV0CHeight2 = 4.4; // height of cell 2
12// fV0CHeight3 = 7.4; // height of cell 3
13// fV0CHeight4 = 12.5; // height of cell 4
14// fV0CRMin = 4.6; // inner radius of box
15//
16// fV0AR0 = 4.2; // Radius of hole
17// fV0AR1 = 7.6; // Maximun radius of 1st cell
18// fV0AR2 = 13.8; // Maximun radius of 2nd cell
19// fV0AR3 = 22.7; // Maximun radius of 3rd cell
20// fV0AR4 = 41.3; // Maximun radius of 4th cell
63192539 21
22void vzero_digits()
23{
24 static const Float_t RadC[] = { 4.6, 7.1, 11.5, 18.9, 31.4 };
25 static const Float_t RadA[] = { 4.2, 7.6, 13.8, 22.7, 41.4 };
26 static const Float_t RadEps = 0.4;
27 static const Float_t PhiEps = 0.025;
28 static const Float_t PhiStp = TMath::TwoPi()/8.0;
29
30 gStyle->SetPalette(1, 0);
31
d810d0de 32 AliRunLoader* rl = AliEveEventManager::AssertRunLoader();
63192539 33 rl->LoadDigits("VZERO");
34
35 TTree* dt = rl->GetTreeD("VZERO", false);
36 TClonesArray* dca = 0;
37 dt->SetBranchAddress("VZERODigit", &dca);
38 dt->GetEntry(0);
39
40 Float_t v[12];
41
84aff7a4 42 TEveQuadSet* modR = new TEveQuadSet("V0R");
43 modR->Reset(TEveQuadSet::kQT_FreeQuad, kFALSE, 32);
63192539 44
84aff7a4 45 TEveQuadSet* modL = new TEveQuadSet("V0L");
46 modL->Reset(TEveQuadSet::kQT_FreeQuad, kFALSE, 32);
51346b82 47
63192539 48 Int_t numEntr = dca->GetEntriesFast();
49 for (Int_t entr=0; entr<numEntr; ++entr)
50 {
51 AliVZEROdigit* d = (AliVZEROdigit*) dca->UncheckedAt(entr);
52 Int_t i = d->PMNumber();
53
d810d0de 54 if (i < 32) // AliEveV0 Right
63192539 55 {
84aff7a4 56 TEveQuadSet* module = modR;
63192539 57 Int_t ri = i / 8;
58 Int_t pi = i % 8;
59 Float_t minR = RadC[ri] + RadEps, maxR = RadC[ri+1] - RadEps;
60 Float_t minP = pi*PhiStp + PhiEps, maxP = (pi+1)*PhiStp - PhiEps;
61
62 v[ 0] = minR*TMath::Cos(minP); v[ 1] = minR*TMath::Sin(minP); v[ 2] = 0;
63 v[ 3] = maxR*TMath::Cos(minP); v[ 4] = maxR*TMath::Sin(minP); v[ 5] = 0;
64 v[ 6] = maxR*TMath::Cos(maxP); v[ 7] = maxR*TMath::Sin(maxP); v[ 8] = 0;
65 v[ 9] = minR*TMath::Cos(maxP); v[10] = minR*TMath::Sin(maxP); v[11] = 0;
66
67 module->AddQuad(v);
68 module->QuadValue(d->ADC());
69 module->QuadId(d);
70 }
d810d0de 71 else // AliEveV0 Left
63192539 72 {
84aff7a4 73 TEveQuadSet* module = modL;
63192539 74 Int_t ri = (i-32) / 8;
75 Int_t pi = i % 8;
76 Float_t minR = RadA[ri] + RadEps, maxR = RadA[ri+1] - RadEps;
77 Float_t minP = pi*PhiStp + PhiEps, maxP = (pi+1)*PhiStp - PhiEps;
78
79 v[ 0] = minR*TMath::Cos(minP); v[ 1] = minR*TMath::Sin(minP); v[ 2] = 0;
80 v[ 3] = maxR*TMath::Cos(minP); v[ 4] = maxR*TMath::Sin(minP); v[ 5] = 0;
81 v[ 6] = maxR*TMath::Cos(maxP); v[ 7] = maxR*TMath::Sin(maxP); v[ 8] = 0;
82 v[ 9] = minR*TMath::Cos(maxP); v[10] = minR*TMath::Sin(maxP); v[11] = 0;
83
84 module->AddQuad(v);
85 module->QuadValue(d->ADC());
86 module->QuadId(d);
87 }
88 }
89
a15e6d7d 90 modL->RefMainTrans().SetPos(0, 0, 324);
91 modR->RefMainTrans().SetPos(0, 0, -84);
63192539 92
84aff7a4 93 gEve->AddElement(modL);
94 gEve->AddElement(modR);
63192539 95
84aff7a4 96 gEve->Redraw3D();
63192539 97}