28328eeb |
1 | /*************************************************************************** |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * |
3 | * * |
4 | * Author: The ALICE Off-line Project. * |
5 | * Contributors are mentioned in the code where appropriate. * |
6 | * * |
7 | * Permission to use, copy, modify and distribute this software and its * |
8 | * documentation strictly for non-commercial purposes is hereby granted * |
9 | * without fee, provided that the above copyright notice appears in all * |
10 | * copies and that both the copyright notice and this permission notice * |
11 | * appear in the supporting documentation. The authors make no claims * |
12 | * about the suitability of this software for any purpose. It is * |
13 | * provided "as is" without express or implied warranty. * |
14 | **************************************************************************/ |
15 | |
16 | #include <Riostream.h> |
17 | #include <TClonesArray.h> |
18 | #include <TTree.h> |
19 | #include <TBranch.h> |
20 | #include <TMath.h> |
95ec17e6 |
21 | #include <TString.h> |
22 | #include <TSystem.h> |
28328eeb |
23 | |
ecee2a1a |
24 | #include "AliLog.h" |
a4e4efaa |
25 | #include "AliRawDataHeader.h" |
cff75ba7 |
26 | #include "AliBitPacking.h" |
28328eeb |
27 | #include "AliPMDdigit.h" |
0a666212 |
28 | #include "AliPMDBlockHeader.h" |
29 | #include "AliPMDDspHeader.h" |
30 | #include "AliPMDPatchBusHeader.h" |
a4e4efaa |
31 | #include "AliPMDRawStream.h" |
28328eeb |
32 | #include "AliPMDDDLRawData.h" |
362c9d61 |
33 | #include "AliDAQ.h" |
28328eeb |
34 | |
35 | ClassImp(AliPMDDDLRawData) |
36 | |
37 | AliPMDDDLRawData::AliPMDDDLRawData(): |
38 | fDigits(new TClonesArray("AliPMDdigit", 1000)) |
39 | { |
40 | // Default Constructor |
41 | // |
a4e4efaa |
42 | |
28328eeb |
43 | } |
a48edddd |
44 | //____________________________________________________________________________ |
45 | AliPMDDDLRawData::AliPMDDDLRawData(const AliPMDDDLRawData& ddlraw): |
46 | TObject(ddlraw), |
47 | fDigits(ddlraw.fDigits) |
48 | { |
49 | //Copy Constructor |
50 | } |
51 | //____________________________________________________________________________ |
52 | AliPMDDDLRawData & AliPMDDDLRawData::operator=(const AliPMDDDLRawData& ddlraw) |
53 | { |
54 | //Assignment operator |
55 | if(this != &ddlraw) |
56 | { |
57 | fDigits = ddlraw.fDigits; |
58 | } |
59 | return *this; |
60 | } |
28328eeb |
61 | //____________________________________________________________________________ |
62 | |
63 | AliPMDDDLRawData::~AliPMDDDLRawData() |
64 | { |
65 | // Default Destructor |
66 | // |
67 | |
68 | } |
69 | |
70 | //____________________________________________________________________________ |
a4e4efaa |
71 | void AliPMDDDLRawData::WritePMDRawData(TTree *treeD) |
28328eeb |
72 | { |
73 | // write digits into raw data format |
74 | |
75 | ofstream outfile; |
76 | |
77 | TBranch *branch = treeD->GetBranch("PMDDigit"); |
ecee2a1a |
78 | if (!branch) |
79 | { |
80 | AliError("PMD Digit branch not found"); |
81 | return; |
82 | } |
28328eeb |
83 | branch->SetAddress(&fDigits); |
84 | |
ecee2a1a |
85 | Int_t nmodules = (Int_t) treeD->GetEntries(); |
86 | AliDebug(1,Form("Number of modules inside treeD = %d",nmodules)); |
28328eeb |
87 | |
362c9d61 |
88 | const Int_t kDDL = AliDAQ::NumberOfDdls("PMD"); |
28328eeb |
89 | Int_t modulePerDDL = 0; |
90 | |
91 | |
a4e4efaa |
92 | AliRawDataHeader header; |
28328eeb |
93 | UInt_t sizeRawData = 0; |
95ec17e6 |
94 | |
95 | const Int_t kSize = 1536; |
28328eeb |
96 | UInt_t buffer[kSize]; |
97 | |
95ec17e6 |
98 | UInt_t busPatch[50][1536]; |
99 | |
100 | Int_t contentsBus[50]; |
101 | |
28328eeb |
102 | Char_t filename[80]; |
103 | |
95ec17e6 |
104 | |
28328eeb |
105 | Int_t mmodule = 0; |
106 | for(Int_t iddl = 0; iddl < kDDL; iddl++) |
107 | { |
362c9d61 |
108 | strcpy(filename,AliDAQ::DdlFileName("PMD",iddl)); |
28328eeb |
109 | #ifndef __DECCXX |
110 | outfile.open(filename,ios::binary); |
111 | #else |
112 | outfile.open(filename); |
113 | #endif |
114 | |
115 | if (iddl < 4) |
116 | { |
117 | modulePerDDL = 6; |
118 | mmodule = 6*iddl; |
28328eeb |
119 | } |
120 | else if (iddl == 4) |
121 | { |
122 | modulePerDDL = 12; |
123 | mmodule = 24; |
28328eeb |
124 | } |
125 | else if (iddl == 5) |
126 | { |
127 | modulePerDDL = 12; |
a4e4efaa |
128 | mmodule = 30; |
28328eeb |
129 | } |
130 | |
28328eeb |
131 | |
95ec17e6 |
132 | |
a4e4efaa |
133 | // Write the Dummy Data Header into the file |
28328eeb |
134 | Int_t bHPosition = outfile.tellp(); |
a4e4efaa |
135 | outfile.write((char*)(&header),sizeof(header)); |
28328eeb |
136 | |
95ec17e6 |
137 | for (Int_t ibus = 0; ibus < 50; ibus++) |
138 | { |
139 | contentsBus[ibus] = 0; |
140 | for (Int_t ich = 0; ich < 1536; ich++) |
141 | { |
142 | busPatch[ibus][ich] = 0; |
143 | } |
144 | } |
28328eeb |
145 | |
146 | for(Int_t ium = 0; ium < modulePerDDL; ium++) |
147 | { |
95ec17e6 |
148 | if (iddl == 4 && ium == 6) mmodule = 42; |
a4e4efaa |
149 | |
28328eeb |
150 | // Extract energy deposition per cell and pack it |
151 | // in a 32-bit word and returns all the total words |
152 | // per one unit-module |
153 | |
95ec17e6 |
154 | GetUMDigitsData(treeD, mmodule, iddl, contentsBus, busPatch); |
28328eeb |
155 | mmodule++; |
95ec17e6 |
156 | } |
28328eeb |
157 | |
95ec17e6 |
158 | |
159 | |
160 | Int_t ij = 0; |
161 | Int_t dsp[10]; |
162 | Int_t dspBus[10]; |
163 | for (Int_t i = 0; i < 10; i++) |
164 | { |
165 | dsp[i] = 0; |
166 | dspBus[i] = 0; |
167 | for (Int_t ibus=0; ibus < 5; ibus++) |
168 | { |
169 | if (contentsBus[ij] > 0) |
170 | { |
171 | dsp[i] += contentsBus[ij]; |
172 | dspBus[i]++; |
173 | } |
174 | ij++; |
175 | } |
176 | // Add the patch Bus header to the DSP contents |
177 | dsp[i] += 4*dspBus[i]; |
178 | } |
179 | |
0a666212 |
180 | Int_t dspBlockARDL = 0; |
181 | Int_t dspBlockBRDL = 0; |
182 | |
95ec17e6 |
183 | for (Int_t i = 0; i < 5; i++) |
184 | { |
185 | Int_t ieven = 2*i; |
186 | Int_t iodd = 2*i + 1; |
187 | if (dsp[ieven] > 0) |
188 | { |
189 | dspBlockARDL += dsp[ieven]; |
0a666212 |
190 | Int_t remainder = dsp[ieven]%2; |
191 | if (remainder == 1) |
192 | { |
193 | dspBlockARDL++; |
194 | } |
95ec17e6 |
195 | } |
196 | if (dsp[iodd] > 0) |
197 | { |
198 | dspBlockBRDL += dsp[iodd]; |
0a666212 |
199 | Int_t remainder = dsp[ieven]%2; |
200 | if (remainder == 1) |
201 | { |
202 | dspBlockARDL++; |
203 | } |
95ec17e6 |
204 | } |
28328eeb |
205 | } |
95ec17e6 |
206 | |
207 | // Start writing the DDL file |
0a666212 |
208 | |
209 | AliPMDBlockHeader blHeader; |
210 | AliPMDDspHeader dspHeader; |
211 | AliPMDPatchBusHeader pbusHeader; |
212 | |
213 | const Int_t kblHLen = blHeader.GetHeaderLength(); |
214 | const Int_t kdspHLen = dspHeader.GetHeaderLength(); |
215 | const Int_t kpbusHLen = pbusHeader.GetHeaderLength(); |
216 | |
95ec17e6 |
217 | UInt_t dspRDL = 0; |
0a666212 |
218 | UInt_t dspBlockHeaderWord[8]; |
219 | UInt_t dspHeaderWord[10]; |
220 | UInt_t patchBusHeaderWord[4]; |
95ec17e6 |
221 | Int_t iskip[5]; |
222 | |
0a666212 |
223 | |
95ec17e6 |
224 | for (Int_t iblock = 0; iblock < 2; iblock++) |
225 | { |
226 | // DSP Block Header |
227 | |
0a666212 |
228 | for (Int_t i=0; i<kblHLen; i++) |
95ec17e6 |
229 | { |
0a666212 |
230 | dspBlockHeaderWord[i] = 0; |
231 | } |
232 | if (iblock == 0) |
233 | { |
234 | dspBlockHeaderWord[1] = (UInt_t) (dspBlockARDL + kblHLen); |
235 | dspBlockHeaderWord[2] = (UInt_t) dspBlockARDL; |
236 | } |
237 | else if (iblock == 1) |
238 | { |
239 | dspBlockHeaderWord[1] = (UInt_t) (dspBlockBRDL + kblHLen); |
240 | dspBlockHeaderWord[2] = (UInt_t) dspBlockBRDL; |
95ec17e6 |
241 | } |
242 | |
0a666212 |
243 | outfile.write((char*)dspBlockHeaderWord,kblHLen*sizeof(UInt_t)); |
95ec17e6 |
244 | |
245 | if (iblock == 0) |
246 | { |
247 | iskip[0] = 0; |
248 | iskip[1] = 10; |
249 | iskip[2] = 20; |
250 | iskip[3] = 30; |
251 | iskip[4] = 40; |
252 | } |
253 | else if (iblock == 1) |
254 | { |
255 | iskip[0] = 5; |
256 | iskip[1] = 15; |
257 | iskip[2] = 25; |
258 | iskip[3] = 35; |
259 | iskip[4] = 45; |
260 | } |
28328eeb |
261 | |
95ec17e6 |
262 | for (Int_t idsp = 0; idsp < 5; idsp++) |
263 | { |
264 | // DSP Header |
265 | Int_t dspno = 0; |
266 | if (iblock == 0) |
267 | { |
268 | dspno = 2*idsp; |
269 | dspRDL = (UInt_t) dsp[dspno]; |
270 | } |
271 | else if (iblock == 1) |
272 | { |
273 | dspno = 2*idsp + 1; |
274 | dspRDL = (UInt_t) dsp[dspno]; |
275 | } |
276 | |
0a666212 |
277 | for (Int_t i=0; i<kdspHLen; i++) |
95ec17e6 |
278 | { |
0a666212 |
279 | dspHeaderWord[i] = 0; |
95ec17e6 |
280 | } |
0a666212 |
281 | Int_t remainder = dspRDL%2; |
282 | if (remainder == 1) dspRDL++; |
283 | |
284 | dspHeaderWord[1] = dspRDL + kdspHLen; |
285 | dspHeaderWord[2] = dspRDL; |
286 | dspHeaderWord[3] = dspno; |
287 | if (remainder == 1) dspHeaderWord[8] = 1; // setting the padding word |
288 | outfile.write((char*)dspHeaderWord,kdspHLen*sizeof(UInt_t)); |
289 | |
95ec17e6 |
290 | for (Int_t ibus = 0; ibus < 5; ibus++) |
291 | { |
292 | // Patch Bus Header |
293 | Int_t busno = iskip[idsp] + ibus; |
294 | Int_t patchbusRDL = contentsBus[busno]; |
0a666212 |
295 | |
296 | if (patchbusRDL > 0) |
95ec17e6 |
297 | { |
0a666212 |
298 | patchBusHeaderWord[0] = 0; |
299 | patchBusHeaderWord[1] = (UInt_t) (patchbusRDL + kpbusHLen); |
300 | patchBusHeaderWord[2] = (UInt_t) patchbusRDL; |
301 | patchBusHeaderWord[3] = (UInt_t) busno; |
95ec17e6 |
302 | } |
0a666212 |
303 | else if (patchbusRDL == 0) |
304 | { |
305 | patchBusHeaderWord[0] = 0; |
306 | patchBusHeaderWord[1] = (UInt_t) kpbusHLen; |
307 | patchBusHeaderWord[2] = (UInt_t) 0; |
308 | patchBusHeaderWord[3] = (UInt_t) busno; |
309 | } |
310 | |
311 | |
312 | outfile.write((char*)patchBusHeaderWord,4*sizeof(UInt_t)); |
95ec17e6 |
313 | |
95ec17e6 |
314 | |
315 | for (Int_t iword = 0; iword < patchbusRDL; iword++) |
316 | { |
317 | buffer[iword] = busPatch[busno][iword]; |
318 | } |
319 | |
320 | outfile.write((char*)buffer,patchbusRDL*sizeof(UInt_t)); |
0a666212 |
321 | |
322 | } // End of patch bus loop |
323 | |
324 | |
325 | // Adding a padding word if the total words odd |
326 | if (remainder == 1) |
327 | { |
328 | UInt_t paddingWord = dspHeader.GetDefaultPaddingWord(); |
329 | outfile.write((char*)(&paddingWord),sizeof(UInt_t)); |
95ec17e6 |
330 | } |
331 | } |
332 | } |
0a666212 |
333 | |
a4e4efaa |
334 | // Write real data header |
335 | // take the pointer to the beginning of the data header |
28328eeb |
336 | // write the total number of words per ddl and bring the |
337 | // pointer to the current file position and close it |
338 | UInt_t cFPosition = outfile.tellp(); |
a4e4efaa |
339 | sizeRawData = cFPosition - bHPosition - sizeof(header); |
95ec17e6 |
340 | |
a4e4efaa |
341 | header.fSize = cFPosition - bHPosition; |
342 | header.SetAttribute(0); // valid data |
28328eeb |
343 | outfile.seekp(bHPosition); |
a4e4efaa |
344 | outfile.write((char*)(&header),sizeof(header)); |
28328eeb |
345 | outfile.seekp(cFPosition); |
346 | |
347 | outfile.close(); |
348 | } // DDL Loop over |
349 | |
350 | |
351 | } |
352 | //____________________________________________________________________________ |
95ec17e6 |
353 | void AliPMDDDLRawData::GetUMDigitsData(TTree *treeD, Int_t imodule, |
354 | Int_t ddlno, Int_t *contentsBus, |
355 | UInt_t busPatch[][1536]) |
28328eeb |
356 | { |
e54787da |
357 | // Retrives digits data UnitModule by UnitModule |
cff75ba7 |
358 | UInt_t baseword; |
28328eeb |
359 | UInt_t mcmno, chno; |
360 | UInt_t adc; |
28328eeb |
361 | Int_t det, smn, irow, icol; |
a98f9d26 |
362 | Int_t parity; |
28328eeb |
363 | |
95ec17e6 |
364 | const Int_t kMaxBus = 50; |
365 | Int_t totPatchBus, bPatchBus, ePatchBus; |
366 | Int_t ibus, totmcm, rows, cols, rowe, cole; |
367 | Int_t moduleno; |
368 | Int_t busno = 0; |
369 | Int_t patchBusNo[kMaxBus], mcmperBus[kMaxBus]; |
370 | Int_t startRowBus[kMaxBus], startColBus[kMaxBus]; |
371 | Int_t endRowBus[kMaxBus], endColBus[kMaxBus]; |
372 | |
373 | Int_t beginPatchBus = -1; |
374 | Int_t endPatchBus = -1; |
375 | for(Int_t i = 0; i < kMaxBus; i++) |
376 | { |
377 | patchBusNo[i] = -1; |
378 | mcmperBus[i] = -1; |
379 | startRowBus[i] = -1; |
380 | startColBus[i] = -1; |
381 | endRowBus[i] = -1; |
382 | endColBus[i] = -1; |
383 | } |
384 | Int_t modulePerDDL = 0; |
385 | if (ddlno < 4) |
386 | { |
387 | modulePerDDL = 6; |
388 | } |
389 | else if (ddlno == 4 || ddlno == 5) |
390 | { |
391 | modulePerDDL = 12; |
392 | } |
393 | |
394 | |
395 | |
396 | |
397 | |
398 | TString fileName(gSystem->Getenv("ALICE_ROOT")); |
399 | |
400 | if(ddlno == 0) |
401 | { |
402 | fileName += "/PMD/PMD_Mapping_ddl0.dat"; |
403 | } |
404 | else if(ddlno == 1) |
405 | { |
406 | fileName += "/PMD/PMD_Mapping_ddl1.dat"; |
407 | } |
408 | else if(ddlno == 2) |
409 | { |
410 | fileName += "/PMD/PMD_Mapping_ddl2.dat"; |
411 | } |
412 | else if(ddlno == 3) |
413 | { |
414 | fileName += "/PMD/PMD_Mapping_ddl3.dat"; |
415 | } |
416 | else if(ddlno == 4) |
417 | { |
418 | fileName += "/PMD/PMD_Mapping_ddl4.dat"; |
419 | } |
420 | else if(ddlno == 5) |
421 | { |
422 | fileName += "/PMD/PMD_Mapping_ddl5.dat"; |
423 | } |
424 | |
425 | ifstream infile; |
426 | infile.open(fileName.Data(), ios::in); // ascii file |
427 | if(!infile) |
428 | AliError(Form("Could not read the mapping file for DDL No = %d",ddlno)); |
429 | |
430 | for (Int_t im = 0; im < modulePerDDL; im++) |
431 | { |
432 | infile >> moduleno; |
433 | infile >> totPatchBus >> bPatchBus >> ePatchBus; |
434 | |
435 | if (moduleno == imodule) |
436 | { |
437 | beginPatchBus = bPatchBus; |
438 | endPatchBus = ePatchBus; |
439 | } |
440 | |
441 | for(Int_t i=0; i<totPatchBus; i++) |
442 | { |
443 | infile >> ibus >> totmcm >> rows >> rowe >> cols >> cole; |
444 | |
445 | if (moduleno == imodule) |
446 | { |
447 | patchBusNo[ibus] = ibus; |
448 | mcmperBus[ibus] = totmcm; |
449 | startRowBus[ibus] = rows; |
450 | startColBus[ibus] = cols; |
451 | endRowBus[ibus] = rowe; |
452 | endColBus[ibus] = cole; |
453 | } |
454 | } |
455 | |
456 | } |
457 | |
458 | infile.close(); |
459 | |
28328eeb |
460 | treeD->GetEntry(imodule); |
461 | Int_t nentries = fDigits->GetLast(); |
95ec17e6 |
462 | Int_t totword = nentries+1; |
28328eeb |
463 | |
0a666212 |
464 | AliPMDdigit *pmddigit; |
a48edddd |
465 | |
28328eeb |
466 | for (Int_t ient = 0; ient < totword; ient++) |
467 | { |
0a666212 |
468 | pmddigit = (AliPMDdigit*)fDigits->UncheckedAt(ient); |
28328eeb |
469 | |
0a666212 |
470 | det = pmddigit->GetDetector(); |
471 | smn = pmddigit->GetSMNumber(); |
472 | irow = pmddigit->GetRow(); |
473 | icol = pmddigit->GetColumn(); |
474 | adc = (UInt_t) pmddigit->GetADC(); |
28328eeb |
475 | |
e54787da |
476 | TransformS2H(smn,irow,icol); |
95ec17e6 |
477 | |
478 | GetMCMCh(ddlno, irow, icol, beginPatchBus, endPatchBus, |
479 | mcmperBus, startRowBus, startColBus, |
480 | endRowBus, endColBus, busno, mcmno, chno); |
28328eeb |
481 | |
28328eeb |
482 | baseword = 0; |
cff75ba7 |
483 | AliBitPacking::PackWord(adc,baseword,0,11); |
484 | AliBitPacking::PackWord(chno,baseword,12,17); |
485 | AliBitPacking::PackWord(mcmno,baseword,18,28); |
486 | AliBitPacking::PackWord(0,baseword,29,30); |
a98f9d26 |
487 | parity = ComputeParity(baseword); // generate the parity bit |
488 | AliBitPacking::PackWord(parity,baseword,31,31); |
95ec17e6 |
489 | |
490 | Int_t jj = contentsBus[busno]; |
491 | busPatch[busno][jj] = baseword; |
492 | |
493 | contentsBus[busno]++; |
28328eeb |
494 | } |
495 | |
e54787da |
496 | } |
95ec17e6 |
497 | |
e54787da |
498 | //____________________________________________________________________________ |
499 | void AliPMDDDLRawData::TransformS2H(Int_t smn, Int_t &irow, Int_t &icol) |
500 | { |
501 | // Does the Software to Hardware coordinate transformation |
502 | // |
95ec17e6 |
503 | |
e54787da |
504 | Int_t irownew = 0; |
505 | Int_t icolnew = 0; |
e54787da |
506 | |
507 | // First in digits we have all dimension 48x96 |
95ec17e6 |
508 | // Transform into the realistic one, i.e, For SM 0&1 96(row)x48(col) |
509 | // and for SM 2&3 48(row)x96(col) |
e54787da |
510 | // |
511 | if(smn < 12) |
512 | { |
95ec17e6 |
513 | irownew = icol; |
514 | icolnew = irow; |
e54787da |
515 | } |
516 | else if( smn >= 12 && smn < 24) |
517 | { |
95ec17e6 |
518 | irownew = irow; |
519 | icolnew = icol; |
e54787da |
520 | } |
521 | |
95ec17e6 |
522 | // In the new geometry always Geant (0,0) and Hardware (0,0) start |
523 | // from the top left corner |
524 | |
e54787da |
525 | irow = irownew; |
526 | icol = icolnew; |
28328eeb |
527 | |
528 | } |
e54787da |
529 | |
530 | |
28328eeb |
531 | //____________________________________________________________________________ |
e54787da |
532 | |
95ec17e6 |
533 | void AliPMDDDLRawData::GetMCMCh(Int_t ddlno, Int_t row, Int_t col, |
534 | Int_t beginPatchBus, Int_t endPatchBus, |
535 | Int_t *mcmperBus, |
536 | Int_t *startRowBus, Int_t *startColBus, |
537 | Int_t *endRowBus, Int_t *endColBus, |
538 | Int_t & busno, UInt_t &mcmno, UInt_t &chno) |
28328eeb |
539 | { |
9ac124a2 |
540 | // This converts row col to hardware channel number |
541 | // This is the final version of mapping supplied by Mriganka |
542 | |
543 | static const UInt_t kCh[16][4] = { {56, 58, 59, 57}, |
544 | {54, 62, 61, 53}, |
545 | {52, 60, 63, 51}, |
546 | {48, 50, 55, 49}, |
547 | {46, 40, 45, 47}, |
548 | {44, 32, 35, 43}, |
549 | {42, 34, 33, 41}, |
550 | {38, 36, 37, 39}, |
551 | {24, 26, 27, 25}, |
552 | {22, 30, 29, 21}, |
553 | {20, 28, 31, 19}, |
554 | {16, 18, 23, 17}, |
555 | {14, 8, 13, 15}, |
556 | {12, 0, 3, 11}, |
557 | {10, 2, 1, 9}, |
558 | {6, 4, 5, 7} }; |
28328eeb |
559 | |
95ec17e6 |
560 | Int_t irownew = row%16; |
561 | Int_t icolnew = col%4; |
562 | |
563 | chno = kCh[irownew][icolnew]; |
28328eeb |
564 | |
28328eeb |
565 | |
95ec17e6 |
566 | for (Int_t ibus = beginPatchBus; ibus <= endPatchBus; ibus++) |
567 | { |
568 | Int_t srow = startRowBus[ibus]; |
569 | Int_t erow = endRowBus[ibus]; |
570 | Int_t scol = startColBus[ibus]; |
571 | Int_t ecol = endColBus[ibus]; |
572 | Int_t tmcm = mcmperBus[ibus]; |
573 | if ((row >= srow && row <= erow) && (col >= scol && col <= ecol)) |
a4e4efaa |
574 | { |
95ec17e6 |
575 | busno = ibus; |
28328eeb |
576 | |
95ec17e6 |
577 | // Find out the MCM Number |
578 | // |
579 | |
580 | if (ddlno == 0 || ddlno == 1) |
581 | { |
582 | // PRE plane, SU Mod = 0, 1 |
583 | Int_t rowdiff = endRowBus[ibus] - startRowBus[ibus]; |
584 | if(rowdiff > 16) |
585 | { |
586 | Int_t midrow = srow + 16; |
587 | if(row >= srow && row < midrow) |
588 | { |
589 | mcmno = 12 + (col-scol)/4; |
590 | } |
591 | else if(row >= midrow && row < erow) |
592 | { |
593 | mcmno = (col-scol)/4; |
594 | } |
595 | } |
596 | else |
597 | { |
598 | mcmno = (col-scol)/4; |
599 | } |
600 | } // end of ddl 0 and 1 |
601 | else if (ddlno == 2) |
602 | { |
603 | // PRE plane, SU Mod = 2 |
604 | |
605 | Int_t icolnew = (col - scol)/4; |
606 | mcmno = tmcm - 1 - icolnew; |
607 | } |
608 | else if (ddlno == 3) |
609 | { |
610 | // PRE plane, SU Mod = 3 |
611 | |
612 | Int_t icolnew = (col - scol)/4; |
613 | mcmno = tmcm - 1 - icolnew; |
614 | } |
615 | else if (ddlno == 4) |
616 | { |
617 | // CPV plane, SU Mod = 0, 3 : ddl = 4 |
618 | |
619 | if(ibus <= 20) |
620 | { |
621 | Int_t rowdiff = endRowBus[ibus] - startRowBus[ibus]; |
622 | if(rowdiff > 16) |
623 | { |
624 | Int_t midrow = srow + 16; |
625 | if(row >= srow && row < midrow) |
626 | { |
627 | mcmno = 12 + (col-scol)/4; |
628 | } |
629 | else if(row >= midrow && row < erow) |
630 | { |
631 | mcmno = (col-scol)/4; |
632 | } |
633 | } |
634 | else |
635 | { |
636 | mcmno = (col-scol)/4; |
637 | } |
638 | } |
639 | else if (ibus > 20) |
640 | { |
641 | Int_t icolnew = (col - scol)/4; |
642 | mcmno = tmcm - 1 - icolnew; |
643 | } |
644 | } |
645 | else if (ddlno == 5) |
646 | { |
647 | // CPV plane, SU Mod = 2, 1 : ddl = 5 |
648 | |
649 | if(ibus <= 20) |
650 | { |
651 | Int_t rowdiff = endRowBus[ibus] - startRowBus[ibus]; |
652 | if(rowdiff > 16) |
653 | { |
654 | Int_t midrow = srow + 16; |
655 | if(row >= srow && row < midrow) |
656 | { |
657 | mcmno = 12 + (col-scol)/4; |
658 | } |
659 | else if(row >= midrow && row < erow) |
660 | { |
661 | mcmno = (col-scol)/4; |
662 | } |
663 | } |
664 | else |
665 | { |
666 | mcmno = (col-scol)/4; |
667 | } |
668 | } |
669 | else if (ibus > 20) |
670 | { |
671 | Int_t icolnew = (col - scol)/4; |
672 | mcmno = tmcm - 1 - icolnew; |
673 | } |
674 | } |
675 | } |
676 | } |
677 | |
28328eeb |
678 | } |
679 | //____________________________________________________________________________ |
680 | |
a98f9d26 |
681 | Int_t AliPMDDDLRawData::ComputeParity(UInt_t baseword) |
682 | { |
683 | // Generate the parity bit |
684 | |
685 | Int_t count = 0; |
686 | for(Int_t j=0; j<29; j++) |
687 | { |
688 | if (baseword & 0x01 ) count++; |
689 | baseword >>= 1; |
690 | } |
691 | Int_t parity = count%2; |
692 | return parity; |
693 | } |
694 | |
695 | //____________________________________________________________________________ |
696 | |