]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ESDCheck/AliPMDQATask.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Adding the new MUONcalib library (Ivana)
[u/mrichter/AliRoot.git] / ESDCheck / AliPMDQATask.cxx
CommitLineData
1dfe075f 1
2/**************************************************************************
3 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * *
5 * Author: The ALICE Off-line Project. *
6 * Contributors are mentioned in the code where appropriate. *
7 * *
8 * Permission to use, copy, modify and distribute this software and its *
9 * documentation strictly for non-commercial purposes is hereby granted *
10 * without fee, provided that the above copyright notice appears in all *
11 * copies and that both the copyright notice and this permission notice *
12 * appear in the supporting documentation. The authors make no claims *
13 * about the suitability of this software for any purpose. It is *
14 * provided "as is" without express or implied warranty. *
15 **************************************************************************/
0b28fd57 16
17/* $Id$ */
18
1dfe075f 19//_________________________________________________________________________
20// An analysis task to check the PMD data in simulated data
21//
22//*-- Yves Schutz
23//////////////////////////////////////////////////////////////////////////////
24
0b28fd57 25#include <TCanvas.h>
1dfe075f 26#include <TChain.h>
0b28fd57 27#include <TFile.h>
1dfe075f 28#include <TH1F.h>
29#include <TH2F.h>
1dfe075f 30#include <TLine.h>
0b28fd57 31#include <TROOT.h>
1dfe075f 32#include <TStyle.h>
1dfe075f 33
34#include "AliPMDQATask.h"
35#include "AliPMDUtility.h"
36#include "AliESD.h"
37#include "AliLog.h"
38
39//______________________________________________________________________________
40AliPMDQATask::AliPMDQATask(const char *name) :
41 AliAnalysisTask(name,""),
42 fChain(0),
43 fESD(0),
44 fhPMDP1(0),
45 fhPMDC2(0),
46 fhPMDP2(0),
47 fhPMDC3(0),
48 fhPMDP3(0),
49 fhPMDP4(0),
50 fhPMDC5(0),
51 fhPMDP5(0),
52 fhPMDCP0(0),
53 fhPMDCP1(0),
54 fhPMDCP2(0),
55 fhPMDCP3(0),
56 fhPMDCP4(0),
57 fhPMDSM1(0),
58 fhPMDSM2(0),
59 fhPMDSM3(0),
60 fhPMDSM4(0),
61 fhPMDSM5(0),
62 fhPMDSM6(0),
63 fhPMDSM7(0),
64 fhPMDSM8(0),
65 fhPMDSM9(0),
66 fhPMDSM10(0),
67 fhPMDSM11(0),
68 fhPMDSM12(0),
69 fhPMDSM13(0),
70 fhPMDSM14(0),
71 fhPMDSM15(0),
72 fhPMDSM16(0),
73 fhPMDSM17(0),
74 fhPMDSM18(0),
75 fhPMDSM19(0),
76 fhPMDSM20(0),
77 fhPMDSM21(0),
78 fhPMDSM22(0),
79 fhPMDSM23(0),
80 fhPMDSM24(0),
81 fhPMDSM (0)
82{
83 // Constructor.
84 // Input slot #0 works with an Ntuple
85 DefineInput(0, TChain::Class());
86 // Output slot #0 writes into a TH1 container
87 DefineOutput(0, TObjArray::Class()) ;
88}
89
90//______________________________________________________________________________
91AliPMDQATask::~AliPMDQATask()
92{
93 // dtor
94 fOutputContainer->Clear() ;
95 delete fOutputContainer ;
96
97 delete fhPMDP1 ;
98 delete fhPMDC2 ;
99 delete fhPMDP2 ;
100 delete fhPMDC3 ;
101 delete fhPMDP3 ;
102 delete fhPMDP4 ;
103 delete fhPMDC5 ;
104 delete fhPMDP5 ;
105 delete fhPMDCP0 ;
106 delete fhPMDCP1 ;
107 delete fhPMDCP2 ;
108 delete fhPMDCP3 ;
109 delete fhPMDCP4 ;
110 delete fhPMDSM1 ;
111 delete fhPMDSM2 ;
112 delete fhPMDSM3 ;
113 delete fhPMDSM4 ;
114 delete fhPMDSM5 ;
115 delete fhPMDSM6 ;
116 delete fhPMDSM7 ;
117 delete fhPMDSM8 ;
118 delete fhPMDSM9 ;
119 delete fhPMDSM10 ;
120 delete fhPMDSM11 ;
121 delete fhPMDSM12 ;
122 delete fhPMDSM13 ;
123 delete fhPMDSM14 ;
124 delete fhPMDSM15 ;
125 delete fhPMDSM16 ;
126 delete fhPMDSM17 ;
127 delete fhPMDSM18 ;
128 delete fhPMDSM19 ;
129 delete fhPMDSM20 ;
130 delete fhPMDSM21 ;
131 delete fhPMDSM22 ;
132 delete fhPMDSM23 ;
133 delete fhPMDSM24 ;
134 delete fhPMDSM ;
135
136}
137
138//______________________________________________________________________________
c52c2132 139void AliPMDQATask::ConnectInputData(const Option_t*)
1dfe075f 140{
141 // Initialisation of branch container and histograms
142
143 AliInfo(Form("*** Initialization of %s", GetName())) ;
144
145 // Get input data
146 fChain = dynamic_cast<TChain *>(GetInputData(0)) ;
147 if (!fChain) {
148 AliError(Form("Input 0 for %s not found\n", GetName()));
149 return ;
150 }
151
c52c2132 152 // One should first check if the branch address was taken by some other task
153 char ** address = (char **)GetBranchAddress(0, "ESD");
154 if (address) {
155 fESD = (AliESD*)(*address);
156 } else {
157 fESD = new AliESD();
158 SetBranchAddress(0, "ESD", &fESD);
1dfe075f 159 }
c52c2132 160}
161
162//________________________________________________________________________
163void AliPMDQATask::CreateOutputObjects()
164{
1dfe075f 165 // create histograms
166
167 fhPMDP1 = new TH2F("fhPMDP1","XY of Clusters",100,-100.,100.,100,-100.,100.);
168 fhPMDC2 = new TH1F("fhPMDC2","CPV PHI",200,-1,9);
169 fhPMDP2 = new TH1F("fhPMDP2","PRE PHI",200,-1,9);
170 fhPMDC3 = new TH1F("fhPMDC3","CPV Clus",30,0.,500.);
171 fhPMDP3 = new TH1F("fhPMDP3","PRE N-gammalike",20,0.,500.);
172 fhPMDP4 = new TH1F("fhPMDP4","PRE EDEP",30,0.,1000.);
173 fhPMDC5 = new TH1F("fhPMDC5","CPV n-cell",20,0.,100.);
174 fhPMDP5 = new TH1F("fhPMDP5","PMD n-cell",20,0.,100.);
175 fhPMDCP0 = new TH2F("fhPMDCP0","PRE CLUS Quad.1 vs 2",150,0.,300.,150,0.,300.);
176 fhPMDCP1 = new TH2F("fhPMDCP1","PRE CLUS Quad.3 vs 4",150,0.,300.,150,0.,300.);
177 fhPMDCP2 = new TH2F("fhPMDCP2","PRE EDEP Quad.3 vs 4",50,0.,300.,50,0.,300.);
178 fhPMDCP3 = new TH2F("fhPMDCP3","PRE EDEP vs Tot Clus ",10,0.,1000.,10,0.,300.);
179 fhPMDCP4 = new TH2F("fhPMDCP4","PRE Clus vs CPV Clus ",150,0.,200.,150,0.,200.);
180
181 fhPMDSM1 = new TH2F("fhPMDSM1","PRE Cluster XY",200,-100,100,200,-100,100);
182 fhPMDSM2 = new TH2F("fhPMDSM2","",999,-100.0,100.0,999,-100.0,100.0);
183 fhPMDSM3 = new TH2F("fhPMDSM3","",999,-100.0,100.0,999,-100.0,100.0);
184 fhPMDSM4 = new TH2F("fhPMDSM4","",999,-100.0,100.0,999,-100.0,100.0);
185 fhPMDSM5 = new TH2F("fhPMDSM5","",999,-100.0,100.0,999,-100.0,100.0);
186 fhPMDSM6 = new TH2F("fhPMDSM6","",999,-100.0,100.0,999,-100.0,100.0);
187 fhPMDSM7 = new TH2F("fhPMDSM7","",999,-100.0,100.0,999,-100.0,100.0);
188 fhPMDSM8 = new TH2F("fhPMDSM8","",999,-100.0,100.0,999,-100.0,100.0);
189 fhPMDSM9 = new TH2F("fhPMDSM9","",999,-100.0,100.0,999,-100.0,100.0);
190 fhPMDSM10 = new TH2F("fhPMDSM10","",999,-100.0,100.0,999,-100.0,100.0);
191 fhPMDSM11 = new TH2F("fhPMDSM11","",999,-100.0,100.0,999,-100.0,100.0);
192 fhPMDSM12 = new TH2F("fhPMDSM12","",999,-100.0,100.0,999,-100.0,100.0);
193 fhPMDSM13 = new TH2F("fhPMDSM13","",999,-100.0,100.0,999,-100.0,100.0);
194 fhPMDSM14 = new TH2F("fhPMDSM14","",999,-100.0,100.0,999,-100.0,100.0);
195 fhPMDSM15 = new TH2F("fhPMDSM15","",999,-100.0,100.0,999,-100.0,100.0);
196 fhPMDSM16 = new TH2F("fhPMDSM16","",999,-100.0,100.0,999,-100.0,100.0);
197 fhPMDSM17 = new TH2F("fhPMDSM17","",999,-100.0,100.0,999,-100.0,100.0);
198 fhPMDSM18 = new TH2F("fhPMDSM18","",999,-100.0,100.0,999,-100.0,100.0);
199 fhPMDSM19 = new TH2F("fhPMDSM19","",999,-100.0,100.0,999,-100.0,100.0);
200 fhPMDSM20 = new TH2F("fhPMDSM20","",999,-100.0,100.0,999,-100.0,100.0);
201 fhPMDSM21 = new TH2F("fhPMDSM21","",999,-100.0,100.0,999,-100.0,100.0);
202 fhPMDSM22 = new TH2F("fhPMDSM22","",999,-100.0,100.0,999,-100.0,100.0);
203 fhPMDSM23 = new TH2F("fhPMDSM23","",999,-100.0,100.0,999,-100.0,100.0);
204 fhPMDSM24 = new TH2F("fhPMDSM24","",999,-100.0,100.0,999,-100.0,100.0);
205 fhPMDSM = new TH1F("fhPMDSM","Plot of all 24 Super Modules",24,0,24);
206
207 // create output container
208
209 fOutputContainer = new TObjArray(38) ;
210 fOutputContainer->SetName("PMD") ;
211
212 fOutputContainer->AddAt(fhPMDP1, 0 );
213 fOutputContainer->AddAt(fhPMDC2, 1 );
214 fOutputContainer->AddAt(fhPMDP2, 2 );
215 fOutputContainer->AddAt(fhPMDC3, 3 );
216 fOutputContainer->AddAt(fhPMDP3, 4 );
217 fOutputContainer->AddAt(fhPMDP4, 5 );
218 fOutputContainer->AddAt(fhPMDC5, 6 );
219 fOutputContainer->AddAt(fhPMDP5, 7 );
220 fOutputContainer->AddAt(fhPMDCP0, 8 );
221 fOutputContainer->AddAt(fhPMDCP1, 9);
222 fOutputContainer->AddAt(fhPMDCP2, 10 );
223 fOutputContainer->AddAt(fhPMDCP3, 11 );
224 fOutputContainer->AddAt(fhPMDCP4, 12 );
225
226 fOutputContainer->AddAt(fhPMDSM1, 13 );
227 fOutputContainer->AddAt(fhPMDSM2, 14 );
228 fOutputContainer->AddAt(fhPMDSM3, 15 );
229 fOutputContainer->AddAt(fhPMDSM4, 16 );
230 fOutputContainer->AddAt(fhPMDSM5, 17 );
231 fOutputContainer->AddAt(fhPMDSM6, 18 );
232 fOutputContainer->AddAt(fhPMDSM7, 19 );
233 fOutputContainer->AddAt(fhPMDSM8, 20 );
234 fOutputContainer->AddAt(fhPMDSM9, 21 );
235 fOutputContainer->AddAt(fhPMDSM10, 22 );
236 fOutputContainer->AddAt(fhPMDSM11, 23 );
237 fOutputContainer->AddAt(fhPMDSM12, 24 );
238 fOutputContainer->AddAt(fhPMDSM13, 25 );
239 fOutputContainer->AddAt(fhPMDSM14, 26 );
240 fOutputContainer->AddAt(fhPMDSM15, 27 );
241 fOutputContainer->AddAt(fhPMDSM16, 28 );
242 fOutputContainer->AddAt(fhPMDSM17, 29 );
243 fOutputContainer->AddAt(fhPMDSM18, 30 );
244 fOutputContainer->AddAt(fhPMDSM19, 31 );
245 fOutputContainer->AddAt(fhPMDSM20, 32 );
246 fOutputContainer->AddAt(fhPMDSM21, 33 );
247 fOutputContainer->AddAt(fhPMDSM22, 34 );
248 fOutputContainer->AddAt(fhPMDSM23, 35 );
249 fOutputContainer->AddAt(fhPMDSM24, 36 );
250 fOutputContainer->AddAt(fhPMDSM, 37 );
251}
252
253//______________________________________________________________________________
254void AliPMDQATask::Exec(Option_t *)
255{
256 // Processing of one event
257
258 Long64_t entry = fChain->GetReadEntry() ;
259
260 if (!fESD) {
261 AliError("fESD is not connected to the input!") ;
262 return ;
263 }
264
265 if ( !((entry-1)%100) )
266 AliInfo(Form("%s ----> Processing event # %lld", (dynamic_cast<TChain *>(fChain))->GetFile()->GetName(), entry)) ;
267
268 // ************************ PMD
269
270 AliPMDUtility *cc = new AliPMDUtility();
271
272 Int_t smn;
273 Int_t n=0;
274 Float_t totCPVClus ;
275 Float_t totPREClus ;
276 Float_t totPREEdep ;
277 Float_t totCPVCell ;
278 Float_t totPRECell ;
279 Float_t preCluQUAD[4] ;
280 Float_t cpvCluQUAD[4] ;
281 Float_t preADCQUAD[4] ;
282 Float_t cpvADCQUAD[4] ;
283 Float_t preCelQUAD[4] ;
284 Float_t cpvCelQUAD[4] ;
285
286 Int_t npmdCl = fESD->GetNumberOfPmdTracks();
287
288 // ****** The loop over PMD clusters
289
290 for (Int_t kk = 0; kk < 4 ; kk++) {
291 cpvCluQUAD[kk] = 0.0 ;
292 preCluQUAD[kk] = 0.0 ;
293 cpvCelQUAD[kk] = 0.0 ;
294 preCelQUAD[kk] = 0.0 ;
295 preADCQUAD[kk] = 0.0 ;
296 }
297
298 while (npmdCl--) {
299
300 AliESDPmdTrack * pmdtr = fESD->GetPmdTrack(npmdCl);
301 Int_t det = pmdtr->GetDetector();
302 Float_t clsX = pmdtr->GetClusterX();
303 Float_t clsY = pmdtr->GetClusterY();
304 Float_t clsZ = pmdtr->GetClusterZ();
305 Float_t ncell = pmdtr->GetClusterCells();
306 Float_t adc = pmdtr->GetClusterADC();
307
308 cc->SetXYZ(clsX,clsY,clsZ);
309 cc->CalculateEta();
310 cc->CalculatePhi();
311 Float_t eta = cc->GetEta();
312 Float_t phi = cc->GetPhi();
313
314 // Calculating S.Module Number from Cluster .
315
316 CalculateSMN(clsX, clsY, smn);
317 if( det == 1)
318 {
319 if(smn >= 0 && smn <= 5) {
320 ++cpvCluQUAD[0] ;
321 cpvADCQUAD[0] =+ adc ;
322 cpvCelQUAD[0] =+ ncell ;
323 }
324 if(smn >= 6 && smn <=11) {
325 ++cpvCluQUAD[1] ;
326 cpvADCQUAD[1] =+ adc ;
327 cpvCelQUAD[1] =+ ncell ;
328 }
329 if(smn >=12 && smn <=17) {
330 ++cpvCluQUAD[2] ;
331 cpvADCQUAD[2] =+ adc ;
332 cpvCelQUAD[2] =+ ncell ;
333 }
334 if(smn >=18 && smn <=23) {
335 ++cpvCluQUAD[3] ;
336 cpvADCQUAD[3] =+ adc ;
337 cpvCelQUAD[3] =+ ncell ;
338 }
339
340 if(eta >= 2.3 && eta <= 3.5)
341 {
342 fhPMDC2->Fill(phi);
343 }
344 }
345 if( det == 0)
346 {
347 if(smn >= 0 && smn <= 5) {
348 ++preCluQUAD[0] ;
349 preADCQUAD[0] =+ adc ;
350 preCelQUAD[0] =+ ncell ;
351 }
352 if(smn >= 6 && smn <=11) {
353 ++preCluQUAD[1] ;
354 preADCQUAD[1] =+ adc ;
355 preCelQUAD[1] =+ ncell ;
356 }
357 if(smn >=12 && smn <=17) {
358 ++preCluQUAD[2] ;
359 preADCQUAD[2] =+ adc ;
360 preCelQUAD[2] =+ ncell ;
361 }
362 if(smn >=18 && smn <=23) {
363 ++preCluQUAD[3] ;
364 preADCQUAD[3] =+ adc ;
365 preCelQUAD[3] =+ ncell ;
366 }
367 if ( n <= 100 ) {
368 fhPMDSM->Fill(smn);
369 if(smn == 0) fhPMDSM1->Fill(-clsX,clsY);
370 if(smn == 0) fhPMDSM1->Fill(-clsX,clsY);
371 if(smn == 1) fhPMDSM2->Fill(-clsX,clsY);
372 if(smn == 2) fhPMDSM3->Fill(-clsX,clsY);
373 if(smn == 3) fhPMDSM4->Fill(-clsX,clsY);
374 if(smn == 4) fhPMDSM5->Fill(-clsX,clsY);
375 if(smn == 5) fhPMDSM6->Fill(-clsX,clsY);
376 if(smn == 6) fhPMDSM7->Fill(-clsX,clsY);
377 if(smn == 7) fhPMDSM8->Fill(-clsX,clsY);
378 if(smn == 8) fhPMDSM9->Fill(-clsX,clsY);
379 if(smn == 9) fhPMDSM10->Fill(-clsX,clsY);
380 if(smn ==10) fhPMDSM11->Fill(-clsX,clsY);
381 if(smn ==11) fhPMDSM12->Fill(-clsX,clsY);
382 if(smn ==12) fhPMDSM13->Fill(-clsX,clsY);
383 if(smn ==13) fhPMDSM14->Fill(-clsX,clsY);
384 if(smn ==14) fhPMDSM15->Fill(-clsX,clsY);
385 if(smn ==15) fhPMDSM16->Fill(-clsX,clsY);
386 if(smn ==16) fhPMDSM17->Fill(-clsX,clsY);
387 if(smn ==17) fhPMDSM18->Fill(-clsX,clsY);
388 if(smn ==18) fhPMDSM19->Fill(-clsX,clsY);
389 if(smn ==19) fhPMDSM20->Fill(-clsX,clsY);
390 if(smn ==20) fhPMDSM21->Fill(-clsX,clsY);
391 if(smn ==21) fhPMDSM22->Fill(-clsX,clsY);
392 if(smn ==22) fhPMDSM23->Fill(-clsX,clsY);
393 if(smn ==23) fhPMDSM24->Fill(-clsX,clsY);
394 }
395 if(eta >= 2.3 && eta <= 3.5)
396 {
397 fhPMDP2->Fill(phi);
398 }
399 fhPMDP1->Fill(clsX,clsY);
400 }
401 }
402 for (Int_t k = 0 ; k < 4 ; k++) {
403 totCPVClus =+ cpvCluQUAD [k] ;
404 totPREClus =+ preCluQUAD [k] ;
405 totCPVCell =+ cpvCelQUAD [k] ;
406 totPRECell =+ preCelQUAD [k] ;
407 totPREEdep =+ preADCQUAD [k] ;
408 }
409 Float_t totCPVpreClus = totPREClus + totCPVClus ;
410
411 // if(eta >= 2.3 && eta <= 3.5) {
412 fhPMDC3->Fill(totCPVClus);
413 fhPMDP3->Fill(totPREClus);
414 fhPMDP4->Fill(totPREEdep);
415 fhPMDP5->Fill(totPRECell);
416 fhPMDCP0->Fill(preCluQUAD[0],preCluQUAD[1]);
417 fhPMDCP1->Fill(preCluQUAD[2],preCluQUAD[3]);
418 fhPMDCP2->Fill(preADCQUAD[2],preADCQUAD[3]);
419 fhPMDCP3->Fill(totPREEdep,totCPVpreClus);
420 fhPMDCP4->Fill(totPREClus,totCPVClus);
421 // }
422 totCPVClus = 0.0;
423 totPREClus = 0.0;
424 totCPVCell = 0.0;
425 totPRECell = 0.0;
426 totPREEdep = 0.0;
427
428 PostData(0, fOutputContainer);
429}
430
431//______________________________________________________________________________
432void AliPMDQATask::Terminate(Option_t *)
433{
434 // Processing when the event loop is ended
c52c2132 435 fOutputContainer = (TObjArray*)GetOutputData(0);
1dfe075f 436
c52c2132 437 fhPMDP1 = (TH2F*)fOutputContainer->At(0);
438 fhPMDC2 = (TH1F*)fOutputContainer->At(1);
439 fhPMDP2 = (TH1F*)fOutputContainer->At(2);
440 fhPMDC3 = (TH1F*)fOutputContainer->At(3);
441 fhPMDP3 = (TH1F*)fOutputContainer->At(4);
442 fhPMDP4 = (TH1F*)fOutputContainer->At(5);
443 fhPMDC5 = (TH1F*)fOutputContainer->At(6);
444 fhPMDP5 = (TH1F*)fOutputContainer->At(7);
445 fhPMDCP0 = (TH2F*)fOutputContainer->At(8);
446 fhPMDCP1 = (TH2F*)fOutputContainer->At(9);
447 fhPMDCP2 = (TH2F*)fOutputContainer->At(10);
448 fhPMDCP3 = (TH2F*)fOutputContainer->At(11);
449 fhPMDCP4 = (TH2F*)fOutputContainer->At(12);
450
451 fhPMDSM1 = (TH2F*)fOutputContainer->At(13);
452 fhPMDSM2 = (TH2F*)fOutputContainer->At(14);
453 fhPMDSM3 = (TH2F*)fOutputContainer->At(15);
454 fhPMDSM4 = (TH2F*)fOutputContainer->At(16);
455 fhPMDSM5 = (TH2F*)fOutputContainer->At(17);
456 fhPMDSM6 = (TH2F*)fOutputContainer->At(18);
457 fhPMDSM7 = (TH2F*)fOutputContainer->At(19);
458 fhPMDSM8 = (TH2F*)fOutputContainer->At(20);
459 fhPMDSM9 = (TH2F*)fOutputContainer->At(21);
460 fhPMDSM10 = (TH2F*)fOutputContainer->At(22);
461 fhPMDSM11 = (TH2F*)fOutputContainer->At(23);
462 fhPMDSM12 = (TH2F*)fOutputContainer->At(24);
463 fhPMDSM13 = (TH2F*)fOutputContainer->At(25);
464 fhPMDSM14 = (TH2F*)fOutputContainer->At(26);
465 fhPMDSM15 = (TH2F*)fOutputContainer->At(27);
466 fhPMDSM16 = (TH2F*)fOutputContainer->At(28);
467 fhPMDSM17 = (TH2F*)fOutputContainer->At(29);
468 fhPMDSM18 = (TH2F*)fOutputContainer->At(30);
469 fhPMDSM19 = (TH2F*)fOutputContainer->At(31);
470 fhPMDSM20 = (TH2F*)fOutputContainer->At(32);
471 fhPMDSM21 = (TH2F*)fOutputContainer->At(33);
472 fhPMDSM22 = (TH2F*)fOutputContainer->At(34);
473 fhPMDSM23 = (TH2F*)fOutputContainer->At(35);
474 fhPMDSM24 = (TH2F*)fOutputContainer->At(36);
475 fhPMDSM = (TH1F*)fOutputContainer->At(37);
476
1dfe075f 477 gStyle->SetOptStat(110000);
478 gStyle->SetOptFit(1);
479
480 TCanvas *cPMD0 = new TCanvas("cPMD0","PMD ESD Test #1", 10,10, 600, 600);
481 cPMD0->Range(-100, -100,100 ,100 );
482 fhPMDSM1->SetMarkerColor(2);
483 fhPMDSM1->Draw();
484 fhPMDSM1->GetXaxis()->SetTitle("Cluster X");
485 fhPMDSM1->GetYaxis()->SetTitle("Cluster Y");
486 fhPMDSM2->SetMarkerColor(2);
487 fhPMDSM2->Draw("same");
488 fhPMDSM3->SetMarkerColor(2);
489 fhPMDSM3->Draw("same");
490 fhPMDSM4->SetMarkerColor(2);
491 fhPMDSM4->Draw("same");
492 fhPMDSM5->SetMarkerColor(2);
493 fhPMDSM5->Draw("same");
494 fhPMDSM6->SetMarkerColor(2);
495 fhPMDSM6->Draw("same");
496 fhPMDSM7->SetMarkerColor(4);
497 fhPMDSM7->Draw("same");
498 fhPMDSM8->SetMarkerColor(4);
499 fhPMDSM8->Draw("same");
500 fhPMDSM9->SetMarkerColor(4);
501 fhPMDSM9->Draw("same");
502 fhPMDSM10->SetMarkerColor(4);
503 fhPMDSM10->Draw("same");
504 fhPMDSM11->SetMarkerColor(4);
505 fhPMDSM11->Draw("same");
506 fhPMDSM12->SetMarkerColor(4);
507 fhPMDSM12->Draw("same");
508 fhPMDSM13->SetMarkerColor(6);
509 fhPMDSM13->Draw("same");
510 fhPMDSM14->SetMarkerColor(6);
511 fhPMDSM14->Draw("same");
512 fhPMDSM15->SetMarkerColor(6);
513 fhPMDSM15->Draw("same");
514 fhPMDSM16->SetMarkerColor(6);
515 fhPMDSM16->Draw("same");
516 fhPMDSM17->SetMarkerColor(6);
517 fhPMDSM17->Draw("same");
518 fhPMDSM18->SetMarkerColor(6);
519 fhPMDSM18->Draw("same");
520 fhPMDSM19->SetMarkerColor(8);
521 fhPMDSM19->Draw("same");
522 fhPMDSM20->SetMarkerColor(8);
523 fhPMDSM20->Draw("same");
524 fhPMDSM21->SetMarkerColor(8);
525 fhPMDSM21->Draw("same");
526 fhPMDSM22->SetMarkerColor(8);
527 fhPMDSM22->Draw("same");
528 fhPMDSM23->SetMarkerColor(8);
529 fhPMDSM23->Draw("same");
530 fhPMDSM24->SetMarkerColor(8);
531 fhPMDSM24->Draw("same");
532
533 DrawPMDBoundary();
534 DrawPMDBoundarySM1();
535 DrawPMDBoundarySM2();
536 DrawPMDBoundarySM3();
537 DrawPMDBoundarySM4();
538 cPMD0->Print("ClusterXY.eps");
539
540 TCanvas *cPMD1 = new TCanvas("cPMD1"," PMD ESD Test #2",10, 10, 600,600);
541 cPMD1->Divide(1,2);
542 cPMD1->cd(1);
543 cPMD1->SetFillColor(0);
544 fhPMDC2->SetLineColor(4);
545 fhPMDC2->Draw();
546 cPMD1->cd(2);
547 fhPMDP2->SetLineColor(2);
548 fhPMDP2->Draw();
549 cPMD1->Print("CPVPREphi.eps");
550
551 TCanvas *cPMD2 = new TCanvas("cPMD2","PMD ESD test #3",10, 10, 600, 600);
552 cPMD2->cd();
553 fhPMDSM->SetFillColor(2);
554 fhPMDSM->Draw();
555 cPMD2->Print("AllSMN.eps");
556
557 TCanvas *cPMD3 = new TCanvas("cPMD3", "PMD ESD test #4",10, 10, 600, 600);
558 cPMD3->Divide(2,2);
559 cPMD3->cd(1);
560 fhPMDCP0->SetMarkerColor(9);
561 fhPMDCP0->Draw();
562 cPMD3->cd(2);
563 fhPMDCP1->SetMarkerColor(6);
564 fhPMDCP1->Draw();
565 cPMD3->cd(3);
566 fhPMDP3->SetLineColor(2);
567 fhPMDP3->Draw();
568 cPMD3->cd(4);
569 fhPMDCP4->SetMarkerColor(3);
570 fhPMDCP4->Draw();
571 cPMD3->Print("CPVPREClus.eps");
572
573 TCanvas *cPMD4 = new TCanvas("cPMD4","PMD ESD test #5", 10, 10, 600, 600);
574 cPMD4->Divide(1,2);
575 cPMD4->cd(1);
576 fhPMDC3->SetLineColor(4);
577 fhPMDC3->Draw();
578 cPMD4->cd(2);
579 fhPMDP4->SetLineColor(2);
580 fhPMDP4->Draw();
581 cPMD4->Print("CPVPREAdc.eps");
582
583 char line[1024] ;
584 sprintf(line, ".!tar -zcvf %s.tar.gz *.eps", GetName()) ;
585 gROOT->ProcessLine(line);
586
587 AliInfo(Form("!!! All the eps files are in %s.tar.gz !!! \n", GetName())) ;
588
589}
590
591//______________________________________________________________________________
592void AliPMDQATask::CalculateSMN( Float_t clsX, Float_t clsY, Int_t & smn)
593{
594 Double_t xcon[96] = {75.133, 54.204, 53.254, 32.326, 31.376,10.447,
595 75.133, 54.204, 53.254, 32.326, 31.376,10.447,
596 75.133, 54.204, 53.254, 32.326, 31.376,10.447,
597 75.133, 54.204, 53.254, 32.326, 31.376,10.447,
598 -75.133, -54.204, -53.254, -32.326, -31.376,-10.447,
599 -75.133, -54.204, -53.254, -32.326, -31.376,-10.447,
600 -75.133, -54.204, -53.254, -32.326, -31.376,-10.447,
601 -75.133, -54.204, -53.254, -32.326, -31.376,-10.447,
602 9.167, -32.543, -33.493, -75.133,
603 9.167, -32.543, -33.493, -75.133,
604 9.167, -32.543, -33.493, -75.133,
605 9.167, -32.543, -33.493, -75.133,
606 9.167, -32.543, -33.493, -75.133,
607 9.167, -32.543, -33.493, -75.133,
608 -9.167, 32.543, 33.493, 75.133,
609 -9.167, 32.543, 33.493, 75.133,
610 -9.167, 32.543, 33.493, 75.133,
611 -9.167, 32.543, 33.493, 75.133,
612 -9.167, 32.543, 33.493, 75.133,
613 -9.167, 32.543, 33.493, 75.133};
614
615 Double_t ycon[96] = {86.475, 86.475, 86.475, 86.475, 86.475, 86.475,
616 38.225, 38.225, 38.225, 38.225, 38.225, 38.225,
617 37.325, 37.325, 37.325, 37.325, 37.325, 37.325,
618 -10.925, -10.925, -10.925, -10.925, -10.925, -10.925,
619 -86.475, -86.475, -86.475, -86.475, -86.475, -86.475,
620 -38.225, -38.225, -38.225, -38.225, -38.225, -38.225,
621 -37.325, -37.325, -37.325, -37.325, -37.325, -37.325,
622 10.925, 10.925, 10.925, 10.925, 10.925, 10.925,
623 86.475, 86.475, 86.475, 86.475,
624 62.225, 62.225, 62.225, 62.225,
625 61.325, 61.325, 61.325, 61.325,
626 37.075, 37.075, 37.075, 37.075,
627 36.175, 36.175, 36.175, 36.175,
628 11.925, 11.925, 11.925 , 11.925,
629 -86.475, -86.475, -86.475, -86.475,
630 -62.225, -62.225, -62.225, -62.225,
631 -61.325, -61.325, -61.325, -61.325,
632 -37.075, -37.075, -37.075, -37.075,
633 -36.175, -36.175, -36.175, -36.175,
634 -11.925, -11.925, -11.925 , -11.925 };
635
636 if((clsX <= xcon[0]) && (clsX >= xcon[1]) &&
637 (clsY <= ycon[0]) && (clsY >= ycon[6])) smn = 0 ;
638
639 else if((clsX <=xcon[2]) && (clsX >= xcon[3]) &&
640 (clsY <= ycon[1]) && (clsY >= ycon[7]))smn = 1 ;
641
642 else if((clsX <=xcon[4]) && (clsX >= xcon[5]) &&
643 (clsY <= ycon[3]) && (clsY >= ycon[8]))smn = 2 ;
644
645 else if((clsX <= xcon[0]) && (clsX >= xcon[1]) &&
646 (clsY <= ycon[12]) && (clsY >= ycon[18])) smn = 3 ;
647
648 else if((clsX <=xcon[2]) && (clsX >= xcon[3]) &&
649 (clsY <= ycon[12]) && (clsY >= ycon[18]))smn = 4 ;
650
651 else if((clsX <=xcon[4]) && (clsX >= xcon[5]) &&
652 (clsY <= ycon[12]) && (clsY >= ycon[18]))smn = 5 ;
653 //------------------------------------------------------------------
654 else if((clsX >= xcon[24]) && (clsX <= xcon[25]) &&
655 (clsY >= ycon[24]) && (clsY <= ycon[30])) smn = 6 ;
656
657 else if((clsX >=xcon[26]) && (clsX <= xcon[27]) &&
658 (clsY >= ycon[25]) && (clsY <= ycon[31]))smn = 7 ;
659
660 else if((clsX >=xcon[28]) && (clsX <= xcon[29]) &&
661 (clsY >= ycon[26]) && (clsY <= ycon[32]))smn = 8 ;
662
663 else if((clsX >= xcon[24]) && (clsX <= xcon[25]) &&
664 (clsY >= ycon[36]) && (clsY <= ycon[42])) smn = 9 ;
665
666 else if((clsX >=xcon[26]) && (clsX <= xcon[27]) &&
667 (clsY >= ycon[36]) && (clsY <= ycon[42]))smn = 10;
668
669 else if((clsX >=xcon[28]) && (clsX <= xcon[29]) &&
670 (clsY >= ycon[36]) && (clsY <= ycon[42]))smn = 11;
671 //------------------------------------------------------------------
672 else if((clsX <= xcon[48]) && (clsX >= xcon[49]) &&
673 (clsY <= ycon[48]) && (clsY >= ycon[52])) smn = 12 ;
674
675 else if((clsX <=xcon[50]) && (clsX >= xcon[51]) &&
676 (clsY <= ycon[48]) && (clsY >= ycon[52]))smn = 13 ;
677
678 else if((clsX <=xcon[48]) && (clsX >= xcon[49]) &&
679 (clsY <= ycon[56]) && (clsY >= ycon[60]))smn = 14 ;
680
681 else if((clsX <=xcon[50]) && (clsX >= xcon[51]) &&
682 (clsY <= ycon[56]) && (clsY >= ycon[60]))smn = 15 ;
683
684 else if((clsX <=xcon[48]) && (clsX >= xcon[49]) &&
685 (clsY <= ycon[64]) && (clsY >= ycon[68]))smn = 16 ;
686
687 else if((clsX <=xcon[50]) && (clsX >= xcon[51]) &&
688 (clsY <= ycon[64]) && (clsY >= ycon[68]))smn = 17 ;
689 //--------------------------------------------------------------
690 else if((clsX >= xcon[72]) && (clsX <= xcon[73]) &&
691 (clsY >= ycon[72]) && (clsY <= ycon[76])) smn = 18 ;
692
693 else if((clsX >=xcon[74]) && (clsX <= xcon[75]) &&
694 (clsY >= ycon[72]) && (clsY <= ycon[76]))smn = 19 ;
695
696 else if((clsX >=xcon[72]) && (clsX <= xcon[73]) &&
697 (clsY >= ycon[80]) && (clsY <= ycon[84]))smn = 20 ;
698
699 else if((clsX >=xcon[74]) && (clsX <= xcon[75]) &&
700 (clsY >= ycon[80]) && (clsY <= ycon[84]))smn = 21;
701
702 else if((clsX >= xcon[72]) && (clsX <= xcon[73]) &&
703 (clsY >= ycon[88]) && (clsY <= ycon[92])) smn = 22 ;
704
705 else if((clsX >=xcon[74]) && (clsX <= xcon[75]) &&
706 (clsY >= ycon[88]) && (clsY <= ycon[92]))smn = 23 ;
707 else smn = 111;
708 }
709
710//______________________________________________________________________________
711void AliPMDQATask::DrawPMDBoundary()
712{
713 // Draw PMD boundaries
714
715 gStyle->SetLineWidth(2);
716 gStyle->SetLineColor(2);
717 TLine * l;
718 l = new TLine(75.1333, 86.475, -75.1333, 86.475); l->Draw("same");
719 l = new TLine(-75.1333, 86.470,-75.1333, -86.475); l->Draw("same");
720 l = new TLine(-75.1333, -86.475,75.1333, -86.475); l->Draw("same");
721 l = new TLine(75.1333, -86.475,75.1333, 86.475); l->Draw("same");
722}
723
724//______________________________________________________________________________
725void AliPMDQATask::DrawPMDBoundarySM1()
726{
727 // Draw boundaries of Super Module 1
728
729 gStyle->SetLineWidth(1);
730 gStyle->SetLineColor(4);
731 TLine * l;
732 l = new TLine(-75.1333, 86.475, -10.447, 86.475); l->Draw("same");
733 l = new TLine(-10.447, 86.475, -10.446, -10.925); l->Draw("same");
734 l = new TLine(-10.446, -10.925, -75.1333,-10.925); l->Draw("same");
735 l = new TLine(-75.1333,-10.925, -75.1333, 86.475); l->Draw("same");
736}
737
738//______________________________________________________________________________
739void AliPMDQATask::DrawPMDBoundarySM2()
740{
741 // Draw boundaries of Super Module 2
742
743 gStyle->SetLineWidth(1);
744 gStyle->SetLineColor(4);
745 TLine * l;
746 l = new TLine(75.1333, -86.475, 10.446, -86.475); l->Draw("same");
747 l = new TLine(10.446, -86.475, 10.446, 10.925); l->Draw("same");
748 l = new TLine(10.446, 10.925, 75.1333, 10.925); l->Draw("same");
749 l = new TLine(75.1333, 10.925, 75.1333, -86.475); l->Draw("same");
750}
751
752
753//______________________________________________________________________________
754void AliPMDQATask::DrawPMDBoundarySM3()
755{
756 // Draw boundaries of Super Module 3
757
758 gStyle->SetLineWidth(1);
759 gStyle->SetLineColor(1);
760 TLine * l;
761 l = new TLine( -9.167, 86.475, 75.1333, 86.475); l->Draw("same");
762 l = new TLine(75.1333,86.475, 75.1333, 11.925); l->Draw("same");
763 l = new TLine(75.1333,11.925, -9.167, 11.925); l->Draw("same");
764 l = new TLine( -9.167, 11.925, -9.167, 86.475); l->Draw("same");
765}
766
767//______________________________________________________________________________
768void AliPMDQATask::DrawPMDBoundarySM4()
769{
770 // Draw boundaries of Super Module 4
771
772 gStyle->SetLineWidth(1);
773 gStyle->SetLineColor(1);
774 TLine * l;
775 l = new TLine(9.167, -86.475, -75.1333,-86.475); l->Draw("same");
776 l = new TLine(-75.1333,-86.475, -75.1333,-11.925); l->Draw("same");
777 l = new TLine(-75.1333,-11.925, 9.167, -11.925); l->Draw("same");
778 l = new TLine(9.167, -11.925, 9.167, -86.475); l->Draw("same");
779}