1 #if !defined(__CINT__) || defined(__MAKECINT__)
4 #include <TStopwatch.h>
8 #include "AliRawReaderDate.h"
9 #include "AliRawReaderRoot.h"
10 #include "AliITSRawStreamSDD.h"
13 // Macro for a z-phi event display of the SDD Raw Data
14 // Origin: F. Prino, prino@to.infn.it
16 void DisplaySDDRawData(Char_t *datafil, Int_t firstEv=0, Int_t lastEv=5){
19 AliITSDDLModuleMapSDD* ddlmap=new AliITSDDLModuleMapSDD();
20 ddlmap->SetJun08Map();
22 TH2F* hzphi3=new TH2F("hzphi3","Layer 3",1536,-0.5,1535.5,3584,-0.5,3584.5);
23 TH2F* hzphi4=new TH2F("hzphi4","Layer 4",2048,-0.5,2047.5,5632,-0.5,5631.5);
25 TLine** lA3=new TLine*[5];
26 for(Int_t ilin=0;ilin<5;ilin++){
27 lA3[ilin]=new TLine((ilin+1)*256,0,(ilin+1)*256,3584.5);
28 lA3[ilin]->SetLineColor(kGray);
29 lA3[ilin]->SetLineStyle(2);
31 TLine** lT3=new TLine*[13];
32 for(Int_t ilin=0;ilin<13;ilin++){
33 lT3[ilin]=new TLine(0,(ilin+1)*256,1535.5,(ilin+1)*256);
34 lT3[ilin]->SetLineColor(kGray);
35 lT3[ilin]->SetLineStyle(2);
38 TLine** lA4=new TLine*[7];
39 for(Int_t ilin=0;ilin<7;ilin++){
40 lA4[ilin]=new TLine((ilin+1)*256,0,(ilin+1)*256,5631.5);
41 lA4[ilin]->SetLineColor(kGray);
42 lA4[ilin]->SetLineStyle(2);
44 TLine** lT4=new TLine*[21];
45 for(Int_t ilin=0;ilin<21;ilin++){
46 lT4[ilin]=new TLine(0,(ilin+1)*256,2047.5,(ilin+1)*256);
47 lT4[ilin]->SetLineColor(kGray);
48 lT4[ilin]->SetLineStyle(2);
56 if(strstr(datafil,".root")!=0){
57 rd=new AliRawReaderRoot(datafil,iev);
59 rd=new AliRawReaderDate(datafil,iev);
61 TStopwatch *evtime=new TStopwatch();
62 TCanvas* c0 = new TCanvas("cd0","c0",800,800);
63 gStyle->SetPalette(1);
67 c0->Divide(1,2,0.001,0.001);
70 printf("Event # %d\n",iev);
74 AliITSRawStreamSDD s(rd);
79 if(s.IsCompletedModule()==kFALSE && s.IsCompletedDDL()==kFALSE){
81 Int_t modID=ddlmap->GetModuleNumber(rd->GetDDLID(),s.GetCarlosId());
82 AliITSgeomTGeo::GetModuleId(modID,lay,lad,det);
83 Int_t iz=s.GetCoord1()+256*(det-1);
84 Int_t iphi=s.GetCoord2()+256*(lad-1)+128*s.GetChannel();
86 hzphi3->SetBinContent(iz+1,iphi+1,s.GetSignal());
88 hzphi4->SetBinContent(iz+1,iphi+1,s.GetSignal());
94 printf("**** Event=%d ID=%d\n",iev,idev);
100 hzphi3->Draw("colz");
101 for(Int_t ilin=0;ilin<5;ilin++) lA3[ilin]->Draw("same");
102 for(Int_t ilin=0;ilin<13;ilin++) lT3[ilin]->Draw("same");
103 hzphi3->GetXaxis()->SetTitle("Z (anode)");
104 hzphi3->GetYaxis()->SetTitle("PHI (time bin)");
107 hzphi4->Draw("colz");
108 for(Int_t ilin=0;ilin<7;ilin++) lA4[ilin]->Draw("same");
109 for(Int_t ilin=0;ilin<21;ilin++) lT4[ilin]->Draw("same");
110 hzphi4->GetXaxis()->SetTitle("Z (anode)");
111 hzphi4->GetYaxis()->SetTitle("PHI (time bin)");
113 }while(rd->NextEvent()&&iev<=lastEv);
117 void DisplaySDDRawData(Int_t nrun, Int_t n2, Int_t nchunk=10, Int_t firstEv=1, Int_t lastEv=20){
118 TGrid::Connect("alien:",0,0,"t");
120 sprintf(filnam,"alien:///alice/data/2008/LHC08c/%09d/raw/08%09d%03d.%02d.root",nrun,nrun,n2,nchunk);
121 printf("Open file %s\n",filnam);
122 DisplaySDDRawData(filnam,firstEv,lastEv);