1 #if !defined(__CINT__) || defined(__MAKECINT__)
6 #include <TStopwatch.h>
10 #include <TLegendEntry.h>
13 #include "AliRawReaderDate.h"
14 #include "AliRawReaderRoot.h"
15 #include "AliITSRawStreamSDD.h"
16 #include "AliITSRawStreamSDDCompressed.h"
20 void CheckDataSizeSDD(TString datafil="12000188359004.16.root",
24 if(datafil.Contains("alien:")) TGrid::Connect("alien:");
25 printf("FILE: %s\n",datafil.Data());
27 gStyle->SetOptStat(0);
28 gStyle->SetTitleFont(42,"XY");
29 gStyle->SetLabelFont(42,"XYZ");
31 Double_t maxOcc=2500.;
32 TH2F* hCellsOnMod=new TH2F("hCellsOnMod","",260,239.5,499.5,500,0.,maxOcc);
33 TH2F* hCellsOnDDL=new TH2F("hCellsOnDDL","",24,-0.5,23.5,500,0.,maxOcc*12);
34 TH1F* hMaxOccMod=new TH1F("hMaxOccMod","",260,239.5,499.5);
35 TH1F* hMaxOccDDL=new TH1F("hMaxOccDDL","",24,-0.5,23.5);
37 AliITSDDLModuleMapSDD* dmap=new AliITSDDLModuleMapSDD();
43 if(datafil.Contains(".root")){
44 rd=new AliRawReaderRoot(datafil.Data(),iev);
46 rd=new AliRawReaderDate(datafil.Data(),iev);
49 Bool_t writtenoutput=kFALSE;
50 Int_t countMod[260],countDDL[24];
51 Int_t maxOccMod[260],maxOccDDL[24];
52 for(Int_t im=0; im<260; im++) maxOccMod[im]=0;
53 for(Int_t id=0; id<24; id++) maxOccDDL[id]=0;
57 printf("Event # %d\n",iev);
58 for(Int_t im=0; im<260; im++) countMod[im]=0;
59 for(Int_t id=0; id<24; id++) countDDL[id]=0;
61 UChar_t cdhAttr=AliITSRawStreamSDD::ReadBlockAttributes(rd);
62 UInt_t amSamplFreq=AliITSRawStreamSDD::ReadAMSamplFreqFromCDH(cdhAttr);
63 AliITSRawStream* s=AliITSRawStreamSDD::CreateRawStreamSDD(rd,cdhAttr);
65 printf("Use %s raw stream, sampling frequency %d MHz\n",s->ClassName(),amSamplFreq);
70 if(s->IsCompletedModule()==kFALSE && s->IsCompletedDDL()==kFALSE){
71 Int_t counts=s->GetSignal();
73 Int_t iDDL=rd->GetDDLID();
75 Int_t iMod=s->GetCarlosId();
76 Int_t modInd=dmap->GetModuleNumber(iDDL,iMod)-240;
77 if(modInd>=0 && modInd<260) ++countMod[modInd];
81 for(Int_t im=0; im<260; im++){
82 hCellsOnMod->Fill(im+240,countMod[im]);
83 if(countMod[im]>maxOccMod[im]) maxOccMod[im]=countMod[im];
85 for(Int_t id=0; id<24; id++){
86 hCellsOnDDL->Fill(id,countDDL[id]);
87 if(countDDL[id]>maxOccDDL[id]) maxOccDDL[id]=countDDL[id];
91 }while(rd->NextEvent()&&iev<=lastEv);
93 for(Int_t im=0; im<260; im++) hMaxOccMod->SetBinContent(im+1,maxOccMod[im]);
94 for(Int_t id=0; id<24; id++) hMaxOccDDL->SetBinContent(id+1,maxOccDDL[id]);
96 TProfile* hAveOccDDL=hCellsOnDDL->ProfileX();
97 TProfile* hAveOccMod=hCellsOnMod->ProfileX();
100 for(Int_t i=0; i<24;i++){
101 hOccDstDDL[i]=(TH1D*)hCellsOnDDL->ProjectionY(Form("hOccDstDDL%d",i),i+1,i+1);
102 printf("DDL %d Entries %.0f\n",i,hOccDstDDL[i]->Integral());
103 totEv=hOccDstDDL[i]->Integral();
105 Double_t xmax=hMaxOccDDL->GetMaximum()*1.1;
106 Int_t ddlForPlot[4]={9,15,17,21};
108 TCanvas* c1= new TCanvas("c1","Module occupancy",1100,750);
112 hCellsOnMod->GetXaxis()->SetTitle("Module Id");
113 hCellsOnMod->GetYaxis()->SetTitle("Number of cells on per event");
114 hCellsOnMod->GetYaxis()->SetTitleOffset(1.3);
115 hCellsOnMod->Draw("colz");
117 hAveOccMod->GetXaxis()->SetTitle("Module Id");
118 hAveOccMod->GetYaxis()->SetTitle("<Occupancy>");
119 hAveOccMod->GetYaxis()->SetTitleOffset(1.3);
122 hMaxOccMod->GetXaxis()->SetTitle("Module Id");
123 hMaxOccMod->GetYaxis()->SetTitle("Maximum Occupancy");
124 hMaxOccMod->GetYaxis()->SetTitleOffset(1.3);
127 TCanvas* c2= new TCanvas("c2","DDL occupancy",1100,750);
131 hCellsOnDDL->GetXaxis()->SetTitle("DDL Number");
132 hCellsOnDDL->GetYaxis()->SetTitle("Number of cells on per event");
133 hCellsOnDDL->GetYaxis()->SetTitleOffset(1.3);
134 hCellsOnDDL->Draw("colz");
135 TLatex* textev=new TLatex(0.15,0.8,Form("%d events",totEv));
140 TLegend* leg=new TLegend(0.6,0.6,0.89,0.89);
141 leg->SetFillColor(0);
142 for(Int_t i=0; i<4; i++){
143 Int_t iddl=ddlForPlot[i];
144 hOccDstDDL[iddl]->GetXaxis()->SetTitle("Occupancy");
145 hOccDstDDL[iddl]->GetXaxis()->SetRangeUser(0.,xmax);
146 hOccDstDDL[iddl]->SetLineColor(i+1);
147 if(i==0) hOccDstDDL[iddl]->Draw();
148 else hOccDstDDL[iddl]->Draw("same");
149 leg->AddEntry(hOccDstDDL[iddl],Form("DDL %d",iddl),"L")->SetTextColor(hOccDstDDL[iddl]->GetLineColor());
153 hAveOccDDL->GetXaxis()->SetTitle("DDL Number");
154 hAveOccDDL->GetYaxis()->SetTitle("<Occupancy>");
155 hAveOccDDL->GetYaxis()->SetTitleOffset(1.3);
158 hMaxOccDDL->GetXaxis()->SetTitle("DDL Number");
159 hMaxOccDDL->GetYaxis()->SetTitle("Maximum Occupancy");
160 hMaxOccDDL->GetYaxis()->SetTitleOffset(1.3);
163 TString outfilname=datafil.Data();
164 if(outfilname.Contains("alien")){
165 outfilname.Remove(0,outfilname.Length()-22);
167 outfilname.Prepend("DataSize_");
168 printf("%s\n",outfilname.Data());
169 TFile* outf=new TFile(outfilname.Data(),"recreate");
170 hCellsOnMod->Write();
172 hCellsOnDDL->Write();