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(TString filename, Int_t firstEv=0, Int_t lastEv=5){
19 Bool_t writtenoutput=kFALSE;
20 AliITSDDLModuleMapSDD* ddlmap=new AliITSDDLModuleMapSDD();
21 ddlmap->SetJun09Map();
23 TH2F* hzphi3=new TH2F("hzphi3","Layer 3",1536,-0.5,1535.5,3584,-0.5,3584.5);
24 TH2F* hzphi4=new TH2F("hzphi4","Layer 4",2048,-0.5,2047.5,5632,-0.5,5631.5);
26 TLine** lA3=new TLine*[5];
27 for(Int_t ilin=0;ilin<5;ilin++){
28 lA3[ilin]=new TLine((ilin+1)*256,0,(ilin+1)*256,3584.5);
29 lA3[ilin]->SetLineColor(kGray);
30 lA3[ilin]->SetLineStyle(2);
32 TLine** lT3=new TLine*[13];
33 for(Int_t ilin=0;ilin<13;ilin++){
34 lT3[ilin]=new TLine(0,(ilin+1)*256,1535.5,(ilin+1)*256);
35 lT3[ilin]->SetLineColor(kGray);
36 lT3[ilin]->SetLineStyle(2);
39 TLine** lA4=new TLine*[7];
40 for(Int_t ilin=0;ilin<7;ilin++){
41 lA4[ilin]=new TLine((ilin+1)*256,0,(ilin+1)*256,5631.5);
42 lA4[ilin]->SetLineColor(kGray);
43 lA4[ilin]->SetLineStyle(2);
45 TLine** lT4=new TLine*[21];
46 for(Int_t ilin=0;ilin<21;ilin++){
47 lT4[ilin]=new TLine(0,(ilin+1)*256,2047.5,(ilin+1)*256);
48 lT4[ilin]->SetLineColor(kGray);
49 lT4[ilin]->SetLineStyle(2);
57 if(filename.Contains(".root")){
58 rd=new AliRawReaderRoot(filename.Data(),iev);
60 rd=new AliRawReaderDate(filename.Data(),iev);
63 TStopwatch *evtime=new TStopwatch();
64 TCanvas* c0 = new TCanvas("cd0","c0",800,800);
65 gStyle->SetPalette(1);
68 c0->Divide(1,2,0.001,0.001);
71 printf("Event # %d\n",iev);
76 UChar_t cdhAttr=AliITSRawStreamSDD::ReadBlockAttributes(rd);
77 UInt_t amSamplFreq=AliITSRawStreamSDD::ReadAMSamplFreqFromCDH(cdhAttr);
78 AliITSRawStream* s=AliITSRawStreamSDD::CreateRawStreamSDD(rd,cdhAttr);
80 printf("Use %s raw stream, sampling frequency %d MHz\n",s->ClassName(),amSamplFreq);
87 if(s->IsCompletedModule()==kFALSE && s->IsCompletedDDL()==kFALSE){
89 Int_t modID=ddlmap->GetModuleNumber(rd->GetDDLID(),s->GetCarlosId());
90 AliITSgeomTGeo::GetModuleId(modID,lay,lad,det);
91 Int_t iz=s->GetCoord1()+256*(det-1);
92 Int_t iphi=s->GetCoord2()+256*(lad-1)+128*s->GetChannel();
94 hzphi3->SetBinContent(iz+1,iphi+1,s->GetSignal());
96 hzphi4->SetBinContent(iz+1,iphi+1,s->GetSignal());
101 printf("**** Event=%d \n",iev);
107 hzphi3->Draw("colz");
108 for(Int_t ilin=0;ilin<5;ilin++) lA3[ilin]->Draw("same");
109 for(Int_t ilin=0;ilin<13;ilin++) lT3[ilin]->Draw("same");
110 hzphi3->GetXaxis()->SetTitle("Z (anode)");
111 hzphi3->GetYaxis()->SetTitle("PHI (time bin)");
114 hzphi4->Draw("colz");
115 for(Int_t ilin=0;ilin<7;ilin++) lA4[ilin]->Draw("same");
116 for(Int_t ilin=0;ilin<21;ilin++) lT4[ilin]->Draw("same");
117 hzphi4->GetXaxis()->SetTitle("Z (anode)");
118 hzphi4->GetYaxis()->SetTitle("PHI (time bin)");
120 }while(rd->NextEvent()&&iev<=lastEv);
124 void DisplaySDDRawData(Int_t nrun, Int_t n2, Int_t chunk=10, Int_t year=2009, Char_t* dir="LHC09b",
126 Int_t lastEv=21*3+1){
127 TGrid::Connect("alien:",0,0,"t");
128 TString filnam(Form("alien:///alice/data/%d/%s/%09d/raw/%02d%09d%03d.%02d.root",year,dir,nrun,year-2000,nrun,n2,chunk));
129 printf("Open file %s\n",filnam.Data());
130 DisplaySDDRawData(filnam,firstEv,lastEv);