]>
Commit | Line | Data |
---|---|---|
91282711 | 1 | #if !defined(__CINT__) || defined(__MAKECINT__) |
2 | #include <TH2F.h> | |
3 | #include <TCanvas.h> | |
4 | #include <TGraphErrors.h> | |
5 | #include <TStopwatch.h> | |
6 | #include <TStyle.h> | |
7 | #include <TLatex.h> | |
8 | #include <TFile.h> | |
9 | #include <TGrid.h> | |
10 | #include <TF1.h> | |
11 | #include <TLine.h> | |
12 | #include "AliRawReader.h" | |
13 | #include "AliRawReaderDate.h" | |
14 | #include "AliRawReaderRoot.h" | |
15 | #include "AliITSOnlineSDDInjectors.h" | |
16 | #include "AliITSRawStreamSDD.h" | |
17 | #include "AliITSDDLModuleMapSDD.h" | |
18 | #endif | |
19 | ||
20 | // Macro for the analysis of PULSER runs (equivalent to ITSSDDINJda.cxx) | |
21 | // Two functions named AnalyzeSDDInjectorsAllModules: | |
22 | // The first is for analyzing a local raw data file and takes as agrument the file name. | |
23 | // The second is for running on ALIEN | |
24 | // All DDLs are analyzed, the argument nDDL selects the DDL to be plotted | |
25 | // Origin: F. Prino (prino@to.infn.it) | |
26 | ||
27 | ||
28 | void AnalyzeSDDInjectorsAllMod(Char_t *datafil, Int_t nDDL, Int_t firstEv=10, Int_t lastEv=15){ | |
29 | Int_t eqOffset = 256; | |
30 | Int_t DDLid_range = 24; | |
31 | // Int_t eqOffset = 100; | |
32 | // Int_t DDLid_range = 1; | |
33 | const Int_t kTotDDL=24; | |
34 | const Int_t kModPerDDL=12; | |
35 | const Int_t kSides=2; | |
36 | ||
37 | AliITSDDLModuleMapSDD* dmap=new AliITSDDLModuleMapSDD(); | |
38 | dmap->SetJun08Map(); | |
39 | ||
40 | TH2F** histo = new TH2F*[kTotDDL*kModPerDDL*kSides]; | |
41 | Int_t nWrittenEv[kTotDDL*kModPerDDL*kSides]; | |
42 | TGraphErrors** gvel = new TGraphErrors*[kTotDDL*kModPerDDL*kSides]; | |
43 | AliITSOnlineSDDInjectors **anal=new AliITSOnlineSDDInjectors*[kTotDDL*kModPerDDL*kSides]; | |
44 | ||
45 | Char_t hisnam[20]; | |
46 | for(Int_t iddl=0; iddl<kTotDDL;iddl++){ | |
47 | for(Int_t imod=0; imod<kModPerDDL;imod++){ | |
48 | for(Int_t isid=0;isid<kSides;isid++){ | |
49 | Int_t index=kSides*(kModPerDDL*iddl+imod)+isid; | |
50 | sprintf(hisnam,"h%02dc%02ds%d",iddl,imod,isid); | |
51 | histo[index]=new TH2F(hisnam,"",256,-0.5,255.5,256,-0.5,255.5); | |
52 | anal[index]=new AliITSOnlineSDDInjectors(iddl,imod,isid); | |
53 | nWrittenEv[index]=0; | |
54 | } | |
55 | } | |
56 | } | |
57 | TGraph *gvvsmod0=new TGraph(0); | |
58 | TGraph *gvvsmod1=new TGraph(0); | |
59 | ||
60 | TCanvas* c0 = new TCanvas("c0","",900,900); | |
61 | gStyle->SetPalette(1); | |
62 | TCanvas* c1 = new TCanvas("c1","",900,900); | |
63 | Char_t text[50]; | |
64 | ||
65 | Int_t iev=firstEv; | |
66 | AliRawReader *rd; | |
67 | if(strstr(datafil,".root")!=0){ | |
68 | rd=new AliRawReaderRoot(datafil,iev); | |
69 | }else{ | |
70 | rd=new AliRawReaderDate(datafil,iev); | |
71 | } | |
72 | ||
73 | Char_t gname[15]; | |
74 | TF1 *funz=new TF1("funz","[0]+[1]*x+[2]*x*x+[3]*x*x*x",0.,255.); | |
75 | TLatex *t0=new TLatex(); | |
76 | t0->SetNDC(); | |
77 | t0->SetTextSize(0.06); | |
78 | t0->SetTextColor(4); | |
79 | ||
80 | do{ | |
81 | c0->Clear(); | |
82 | c0->Divide(4,6,0.001,0.001); | |
83 | c1->Clear(); | |
84 | c1->Divide(4,6,0.001,0.001); | |
85 | printf("Event # %d\n",iev); | |
86 | UInt_t timeSt=rd->GetTimestamp(); | |
87 | //rd->SelectEvents(7); | |
88 | rd->SelectEquipment(17,eqOffset,eqOffset+DDLid_range); | |
89 | rd->Reset(); | |
90 | for(Int_t iddl=0; iddl<kTotDDL;iddl++){ | |
91 | for(Int_t imod=0; imod<kModPerDDL;imod++){ | |
92 | for(Int_t isid=0;isid<kSides;isid++){ | |
93 | Int_t index=kSides*(kModPerDDL*iddl+imod)+isid; | |
94 | histo[index]->Reset(); | |
95 | } | |
96 | } | |
97 | } | |
98 | ||
99 | AliITSRawStreamSDD s(rd); | |
100 | rd->SelectEquipment(17,eqOffset,eqOffset+DDLid_range); | |
101 | while(s.Next()){ | |
102 | Int_t iDDL=rd->GetDDLID(); | |
103 | Int_t iCarlos=s.GetCarlosId(); | |
104 | if(iDDL>=0 && iDDL<kTotDDL && s.IsCompletedModule()==kFALSE){ | |
105 | Int_t index=kSides*(kModPerDDL*iDDL+iCarlos)+s.GetChannel(); | |
106 | histo[index]->Fill(s.GetCoord2(),s.GetCoord1(),s.GetSignal()); | |
107 | } | |
108 | } | |
109 | ||
110 | for(Int_t iddl=0; iddl<kTotDDL;iddl++){ | |
111 | for(Int_t imod=0; imod<kModPerDDL;imod++){ | |
112 | for(Int_t isid=0;isid<kSides;isid++){ | |
113 | Int_t index=kSides*(kModPerDDL*iddl+imod)+isid; | |
114 | anal[index]->AnalyzeEvent(histo[index]); | |
115 | anal[index]->WriteToASCII(iev,timeSt,nWrittenEv[index]); | |
116 | nWrittenEv[index]++; | |
117 | if(iev==firstEv && anal[index]->GetInjPadStatus(16)>=6){ | |
118 | Float_t vel=anal[index]->GetDriftSpeed(16); | |
119 | Int_t iMod=dmap->GetModuleNumber(iddl,imod); | |
120 | if(isid==0) gvvsmod0->SetPoint(gvvsmod0->GetN(),(Float_t)iMod,vel); | |
121 | if(isid==1) gvvsmod1->SetPoint(gvvsmod1->GetN(),(Float_t)iMod,vel); | |
122 | } | |
123 | if(iddl==nDDL){ | |
124 | Int_t index2=kSides*imod+isid; | |
125 | c0->cd(index2+1); | |
126 | histo[index]->SetMaximum(100.); | |
127 | histo[index]->DrawCopy("colz"); | |
128 | sprintf(text,"DDL %d channel %d Side %d",nDDL,imod,isid); | |
129 | t0->DrawLatex(0.15,0.92,text); | |
130 | c0->Update(); | |
131 | c1->cd(index2+1); | |
132 | gvel[index]=anal[index]->GetDriftSpeedGraph(); | |
133 | gvel[index]->SetMarkerStyle(20); | |
134 | gvel[index]->SetTitle(""); | |
135 | sprintf(gname,"gvel%dev%d",index,iev); | |
136 | gvel[index]->SetName(gname); | |
137 | // gvel[index]->SetMinimum(0); | |
138 | //gvel[index]->SetMaximum(200); | |
139 | ||
140 | gvel[index]->GetXaxis()->SetLimits(0,256); | |
141 | gvel[index]->GetXaxis()->SetTitle("Anode"); | |
142 | gvel[index]->GetYaxis()->SetTitle("Drift vel."); | |
143 | gvel[index]->GetXaxis()->SetTitleSize(0.07); | |
144 | gvel[index]->GetYaxis()->SetTitleSize(0.07); | |
145 | gvel[index]->GetXaxis()->SetTitleOffset(0.6); | |
146 | gvel[index]->GetYaxis()->SetTitleOffset(0.6); | |
147 | if(gvel[index]->GetN()>0) gvel[index]->Draw("AP"); | |
148 | Float_t *param=anal[index]->GetDriftSpeedFitParam(); | |
149 | funz->SetParameters(param[0],param[1],param[2],param[3]); | |
150 | funz->SetLineColor(2); | |
151 | funz->DrawCopy("LSAME"); | |
152 | t0->DrawLatex(0.15,0.92,text); | |
153 | c1->Update(); | |
154 | } | |
155 | } | |
156 | } | |
157 | } | |
158 | iev++; | |
159 | printf(" --- OK\n"); | |
160 | }while(rd->NextEvent()&&iev<=lastEv); | |
161 | ||
162 | TCanvas* c8=new TCanvas("c8"); | |
163 | gvvsmod0->SetTitle(""); | |
164 | gvvsmod1->SetTitle(""); | |
165 | ||
166 | gvvsmod0->SetMarkerStyle(20); | |
167 | gvvsmod1->SetMarkerStyle(21); | |
168 | gvvsmod1->SetMarkerColor(2); | |
169 | gvvsmod0->Draw("AP"); | |
170 | gvvsmod0->SetMinimum(6.2); | |
171 | gvvsmod0->SetMaximum(7.2); | |
172 | gvvsmod0->GetXaxis()->SetTitle("Module Number"); | |
173 | gvvsmod0->GetYaxis()->SetTitle("Vdrift at injector pad 16"); | |
174 | gvvsmod1->Draw("PSAME"); | |
175 | TLatex* tleft=new TLatex(0.7,0.82,"Side 0"); | |
176 | tleft->SetNDC(); | |
177 | tleft->SetTextColor(1); | |
178 | tleft->Draw(); | |
179 | TLatex* tright=new TLatex(0.7,0.75,"Side 1"); | |
180 | tright->SetNDC(); | |
181 | tright->SetTextColor(2); | |
182 | tright->Draw(); | |
183 | ||
184 | TLine *lin=new TLine(323,6.2,323,7.2); | |
185 | lin->SetLineColor(4); | |
186 | lin->Draw(); | |
187 | c8->Update(); | |
188 | } | |
189 | ||
190 | void AnalyzeSDDInjectorsAllMod(Int_t nrun, Int_t n2, Int_t nDDL=0, Int_t firstEv=10, Int_t lastEv=15){ | |
191 | TGrid::Connect("alien:",0,0,"t"); | |
192 | Char_t filnam[200]; | |
193 | sprintf(filnam,"alien:///alice/data/2008/LHC08c_SDD/%09d/raw/08%09d%03d.10.root",nrun,nrun,n2); | |
194 | printf("Open file %s\n",filnam); | |
195 | AnalyzeSDDInjectorsAllMod(filnam,nDDL,firstEv,lastEv); | |
196 | } | |
197 | ||
198 | ||
199 |