]> git.uio.no Git - u/mrichter/AliRoot.git/blame - RICH/AliRICHDisplFast.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
digitization under AliSimulation fixed, all the information in chamber's local system...
[u/mrichter/AliRoot.git] / RICH / AliRICHDisplFast.cxx
CommitLineData
6b17b320 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
6b17b320 15#include "AliRICHDisplFast.h"
53fd478b 16#include "AliRICH.h"
17#include "AliRICHChamber.h"
6b17b320 18#include "AliRICHParam.h"
19#include <AliLoader.h>
6b17b320 20#include <TCanvas.h>
af3d25a6 21#include <TPolyLine.h>
6b17b320 22#include <TParticle.h>
23#include <TStyle.h>
6b17b320 24#include <TH2.h>
25#include <TMath.h>
6b17b320 26
aed240d4 27ClassImp(AliRICHDisplFast)
6b17b320 28
aed240d4 29//__________________________________________________________________________________________________
30void AliRICHDisplFast::Exec()
6b17b320 31{
af3d25a6 32 TH2F *pHitsH2 = new TH2F("pHitsH2" , "Event Display",165,-AliRICHParam::PcSizeX()/2,AliRICHParam::PcSizeX()/2,
33 144,-AliRICHParam::PcSizeY()/2,AliRICHParam::PcSizeY()/2);
aed240d4 34 pHitsH2->SetStats(kFALSE);
af3d25a6 35 TH2F *pDigitsH2 = new TH2F("pDigitsH2" ,"Event Display",165,-AliRICHParam::PcSizeX()/2,AliRICHParam::PcSizeX()/2,
36 144,-AliRICHParam::PcSizeY()/2,AliRICHParam::PcSizeY()/2);
37 TH2F *pClustersH2 = new TH2F("pClustersH2","Event Display",165,-AliRICHParam::PcSizeX()/2,AliRICHParam::PcSizeX()/2,
38 144,-AliRICHParam::PcSizeY()/2,AliRICHParam::PcSizeY()/2);
aed240d4 39
40 TCanvas *Display = new TCanvas("Display","RICH Display",0,0,600,600);
6b17b320 41
c94cc386 42 gStyle->SetPalette(1);
43
6b17b320 44 AliRICH *pRich = (AliRICH*)gAlice->GetDetector("RICH");
c94cc386 45 pRich->GetLoader()->LoadRecPoints();
6b17b320 46 pRich->GetLoader()->LoadDigits();
6b17b320 47 pRich->GetLoader()->LoadHits();
aed240d4 48 for(Int_t iEventN=0;iEventN<gAlice->GetEventsPerRun();iEventN++){//events Loop
49 pRich->GetLoader()->GetRunLoader()->GetEvent(iEventN);
50 pRich->GetLoader()->TreeD()->GetEntry(0);
51 pRich->GetLoader()->TreeR()->GetEntry(0);
6b17b320 52
c94cc386 53
aed240d4 54 Int_t nPrimaries = (Int_t)pRich->GetLoader()->TreeH()->GetEntries();
55 TObjArray * Hits = new TObjArray[nPrimaries];
6b17b320 56
aed240d4 57 for(Int_t i=0;i<nPrimaries;i++) {
58 pRich->GetLoader()->TreeH()->GetEntry(i);
59 Int_t nHits = pRich->Hits()->GetEntries();
60 for(Int_t k=0;k<nHits;k++) Hits[i].Add(pRich->Hits()->At(k));
61
6b17b320 62 }
6b17b320 63
aed240d4 64 for(Int_t iChamber=1;iChamber<=7;iChamber++){//modules loop
6b17b320 65
aed240d4 66 Int_t nDigits = pRich->Digits(iChamber)->GetEntries();
c94cc386 67 Int_t nClusters = pRich->Clusters(iChamber)->GetEntries();
6b17b320 68
aed240d4 69 pHitsH2->Reset(); pDigitsH2->Reset(); pClustersH2->Reset();
6b17b320 70
6b17b320 71
aed240d4 72 for(Int_t i=0;i<nPrimaries;i++){//prims loop
6b17b320 73 pRich->GetLoader()->TreeH()->GetEntry(i);
74 Int_t nHits = pRich->Hits()->GetEntries();
aed240d4 75 for(Int_t j=0;j<nHits;j++){//hits loop
6b17b320 76 AliRICHhit *pHit = (AliRICHhit*)Hits[i].At(j);
aed240d4 77 if(pHit->C()==iChamber){
6b17b320 78 TVector3 xyzhit(pHit->X(),pHit->Y(),pHit->Z());
3582c1f9 79 TVector2 hitlocal = pRich->C(iChamber)->Glob2Loc(xyzhit);
aed240d4 80 pHitsH2->Fill(hitlocal.X(),hitlocal.Y(),200);
81 }//if
82 }//hits loop
83 }//prims loop
c94cc386 84
aed240d4 85 for(Int_t j=0;j<nDigits;j++){//digits loop
86 AliRICHdigit *pDigit = (AliRICHdigit*)pRich->Digits(iChamber)->At(j);
3582c1f9 87 TVector2 x2=AliRICHParam::Pad2Loc(pDigit->X(),pDigit->Y());
88 pDigitsH2->Fill(x2.X(),x2.Y(),100);
aed240d4 89 }//digits loop
c94cc386 90
aed240d4 91 for(Int_t j=0;j<nClusters;j++){//clusters loop
92 AliRICHcluster *pCluster = (AliRICHcluster*)pRich->Clusters(iChamber)->At(j);
93 pClustersH2->Fill(pCluster->X(),pCluster->Y(),50);
94 }//clusters loop
95
96 pHitsH2->SetTitle(Form("event %i module %2i",iEventN,iChamber));
97 pHitsH2->SetMarkerColor(kRed); pHitsH2->SetMarkerStyle(29); pHitsH2->SetMarkerSize(0.4);
98 pHitsH2->Draw();
af3d25a6 99 DrawSectors();
6b17b320 100 Display->Update();
101 Display->Modified();
6b17b320 102 getchar();
aed240d4 103
104 pDigitsH2->SetMarkerColor(kGreen); pDigitsH2->SetMarkerStyle(29); pDigitsH2->SetMarkerSize(0.4);
105 pDigitsH2->Draw("same");
6b17b320 106 Display->Update();
aed240d4 107 Display->Modified();
6b17b320 108 getchar();
aed240d4 109
110 pClustersH2->SetMarkerColor(kBlue); pClustersH2->SetMarkerStyle(29); pClustersH2->SetMarkerSize(0.4);
111 pClustersH2->Draw("same");
c94cc386 112 Display->Update();
113 Display->Modified();
114 getchar();
aed240d4 115 }//modules loop
116 delete [] Hits;
117 }////events Loop
118 pRich->GetLoader()->UnloadRecPoints();
119 pRich->GetLoader()->UnloadDigits();
120 pRich->GetLoader()->UnloadHits();
121}//Exec()
122//__________________________________________________________________________________________________
af3d25a6 123void AliRICHDisplFast::DrawSectors()
124{
125 Double_t x1[5] = {-AliRICHParam::PcSizeX()/2,
126 -AliRICHParam::SectorSizeX()/2-AliRICHParam::DeadZone(),
127 -AliRICHParam::SectorSizeX()/2-AliRICHParam::DeadZone(),
128 -AliRICHParam::PcSizeX()/2,
129 -AliRICHParam::PcSizeX()/2};
130 Double_t y1[5] = {AliRICHParam::DeadZone()/2,
131 AliRICHParam::DeadZone()/2,
132 AliRICHParam::PcSizeY()/2,
133 AliRICHParam::PcSizeY()/2,
134 AliRICHParam::DeadZone()/2};
135 Double_t x2[5] = {-AliRICHParam::SectorSizeX()/2,
136 AliRICHParam::SectorSizeX()/2,
137 AliRICHParam::SectorSizeX()/2,
138 -AliRICHParam::SectorSizeX()/2,
139 -AliRICHParam::SectorSizeX()/2};
140 Double_t y2[5] = {AliRICHParam::DeadZone()/2,
141 AliRICHParam::DeadZone()/2,
142 AliRICHParam::PcSizeY()/2,
143 AliRICHParam::PcSizeY()/2,
144 AliRICHParam::DeadZone()/2};
145 Double_t x3[5] = { AliRICHParam::SectorSizeX()/2+AliRICHParam::DeadZone(),
146 AliRICHParam::PcSizeX()/2,
147 AliRICHParam::PcSizeX()/2,
148 AliRICHParam::SectorSizeX()/2+AliRICHParam::DeadZone(),
149 AliRICHParam::SectorSizeX()/2+AliRICHParam::DeadZone()};
150 Double_t y3[5] = {AliRICHParam::DeadZone()/2,
151 AliRICHParam::DeadZone()/2,
152 AliRICHParam::PcSizeY()/2,
153 AliRICHParam::PcSizeY()/2,
154 AliRICHParam::DeadZone()/2};
155 Double_t x4[5] = {-AliRICHParam::PcSizeX()/2,
156 -AliRICHParam::SectorSizeX()/2-AliRICHParam::DeadZone(),
157 -AliRICHParam::SectorSizeX()/2-AliRICHParam::DeadZone(),
158 -AliRICHParam::PcSizeX()/2,
159 -AliRICHParam::PcSizeX()/2};
160 Double_t y4[5] = {-AliRICHParam::PcSizeY()/2,
161 -AliRICHParam::PcSizeY()/2,
162 -AliRICHParam::DeadZone()/2,
163 -AliRICHParam::DeadZone()/2,
164 -AliRICHParam::PcSizeY()/2};
165 Double_t x5[5] = {-AliRICHParam::SectorSizeX()/2,
166 AliRICHParam::SectorSizeX()/2,
167 AliRICHParam::SectorSizeX()/2,
168 -AliRICHParam::SectorSizeX()/2,
169 -AliRICHParam::SectorSizeX()/2};
170 Double_t y5[5] = {-AliRICHParam::PcSizeY()/2,
171 -AliRICHParam::PcSizeY()/2,
172 -AliRICHParam::DeadZone()/2,
173 -AliRICHParam::DeadZone()/2,
174 -AliRICHParam::PcSizeY()/2};
175 Double_t x6[5] = { AliRICHParam::SectorSizeX()/2+AliRICHParam::DeadZone(),
176 AliRICHParam::PcSizeX()/2,
177 AliRICHParam::PcSizeX()/2,
178 AliRICHParam::SectorSizeX()/2+AliRICHParam::DeadZone(),
179 AliRICHParam::SectorSizeX()/2+AliRICHParam::DeadZone()};
180 Double_t y6[5] = {-AliRICHParam::PcSizeY()/2,
181 -AliRICHParam::PcSizeY()/2,
182 -AliRICHParam::DeadZone()/2,
183 -AliRICHParam::DeadZone()/2,
184 -AliRICHParam::PcSizeY()/2};
185 TPolyLine *sector1 = new TPolyLine(5,x1,y1);
186 TPolyLine *sector2 = new TPolyLine(5,x2,y2);
187 TPolyLine *sector3 = new TPolyLine(5,x3,y3);
188 TPolyLine *sector4 = new TPolyLine(5,x4,y4);
189 TPolyLine *sector5 = new TPolyLine(5,x5,y5);
190 TPolyLine *sector6 = new TPolyLine(5,x6,y6);
191 sector1->SetLineColor(21);
192 sector2->SetLineColor(21);
193 sector3->SetLineColor(21);
194 sector4->SetLineColor(21);
195 sector5->SetLineColor(21);
196 sector6->SetLineColor(21);
197 sector1->Draw();
198 sector2->Draw();
199 sector3->Draw();
200 sector4->Draw();
201 sector5->Draw();
202 sector6->Draw();
203}