]> git.uio.no Git - u/mrichter/AliRoot.git/blame_incremental - TRD/AliTRDfeeParam.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame_incremental
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (non-incremental) | history | HEAD
New macro to keep track of timing performances of the segmentation methods (Laurent)
[u/mrichter/AliRoot.git] / TRD / AliTRDfeeParam.cxx
... / ...
CommitLineData
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 **************************************************************************/
16
17/* $Id$ */
18
19////////////////////////////////////////////////////////////////////////////
20// //
21// TRD front end electronics parameters class //
22// Contains all FEE (MCM, TRAP, PASA) related //
23// parameters, constants, and mapping. //
24// //
25// New release on 2007/08/17: //
26// The default raw data version (now fRAWversion ) is set to 3 //
27// in the constructor because version 3 raw data read and write //
28// are fully debugged. //
29// //
30// Author: //
31// Ken Oyama (oyama@physi.uni-heidelberg.de) //
32// //
33////////////////////////////////////////////////////////////////////////////
34
35#include <TMath.h>
36
37#include "AliLog.h"
38
39#include "AliTRDfeeParam.h"
40#include "AliTRDgeometry.h"
41#include "AliTRDCommonParam.h"
42
43ClassImp(AliTRDfeeParam)
44
45AliTRDfeeParam *AliTRDfeeParam::fgInstance = 0;
46Bool_t AliTRDfeeParam::fgTerminated = kFALSE;
47
48//_____________________________________________________________________________
49AliTRDfeeParam* AliTRDfeeParam::Instance()
50{
51 //
52 // Instance constructor
53 //
54
55 if (fgTerminated != kFALSE) {
56 return 0;
57 }
58
59 if (fgInstance == 0) {
60 fgInstance = new AliTRDfeeParam();
61 }
62
63 return fgInstance;
64
65}
66
67//_____________________________________________________________________________
68void AliTRDfeeParam::Terminate()
69{
70 //
71 // Terminate the class and release memory
72 //
73
74 fgTerminated = kTRUE;
75
76 if (fgInstance != 0) {
77 delete fgInstance;
78 fgInstance = 0;
79 }
80
81}
82
83//_____________________________________________________________________________
84AliTRDfeeParam::AliTRDfeeParam()
85 :TObject()
86 ,fCP(0)
87 ,fTFnExp(1)
88 ,fTFr1(0)
89 ,fTFr2(0)
90 ,fTFc1(0)
91 ,fTFc2(0)
92 ,fEBsglIndThr(5)
93 ,fEBsumIndThr(5)
94 ,fEBindLUT(0xF0)
95 ,fEBignoreNeighbour(0)
96 ,fRAWversion(3)
97 ,fRAWstoreRaw(kTRUE)
98{
99 //
100 // Default constructor
101 //
102
103 // PASA V.4
104 if (fTFnExp == 1) {
105 fTFr1 = 1.1563;
106 fTFr2 = 0.1299;
107 fTFc1 = 0.0657;
108 fTFc2 = 0.0000;
109 }
110 else if (fTFnExp == 2) {
111 fTFr1 = 1.1563;
112 fTFr2 = 0.1299;
113 fTFc1 = 0.1141;
114 fTFc2 = 0.6241;
115 }
116
117 fCP = AliTRDCommonParam::Instance();
118
119}
120
121//_____________________________________________________________________________
122AliTRDfeeParam::AliTRDfeeParam(TRootIoCtor *)
123 :TObject()
124 ,fCP(0)
125 ,fTFnExp(1)
126 ,fTFr1(0)
127 ,fTFr2(0)
128 ,fTFc1(0)
129 ,fTFc2(0)
130 ,fEBsglIndThr(0)
131 ,fEBsumIndThr(0)
132 ,fEBindLUT(0)
133 ,fEBignoreNeighbour(0)
134 ,fRAWversion(0)
135 ,fRAWstoreRaw(0)
136{
137 //
138 // IO constructor
139 //
140
141}
142
143//_____________________________________________________________________________
144AliTRDfeeParam::AliTRDfeeParam(const AliTRDfeeParam &p)
145 :TObject(p)
146 ,fCP(p.fCP)
147 ,fTFnExp(p.fTFnExp)
148 ,fTFr1(p.fTFr1)
149 ,fTFr2(p.fTFr2)
150 ,fTFc1(p.fTFc1)
151 ,fTFc2(p.fTFc2)
152 ,fEBsglIndThr(p.fEBsglIndThr)
153 ,fEBsumIndThr(p.fEBsumIndThr)
154 ,fEBindLUT(p.fEBindLUT)
155 ,fEBignoreNeighbour (p.fEBignoreNeighbour)
156 ,fRAWversion(p.fRAWversion)
157 ,fRAWstoreRaw(p.fRAWstoreRaw)
158{
159 //
160 // AliTRDfeeParam copy constructor
161 //
162
163}
164
165//_____________________________________________________________________________
166AliTRDfeeParam::~AliTRDfeeParam()
167{
168 //
169 // AliTRDfeeParam destructor
170 //
171
172}
173
174//_____________________________________________________________________________
175AliTRDfeeParam &AliTRDfeeParam::operator=(const AliTRDfeeParam &p)
176{
177 //
178 // Assignment operator
179 //
180
181 if (this != &p) {
182 ((AliTRDfeeParam &) p).Copy(*this);
183 }
184
185 return *this;
186
187}
188
189//_____________________________________________________________________________
190void AliTRDfeeParam::Copy(TObject &p) const
191{
192 //
193 // Copy function
194 //
195
196 ((AliTRDfeeParam &) p).fCP = fCP;
197 ((AliTRDfeeParam &) p).fTFnExp = fTFnExp;
198 ((AliTRDfeeParam &) p).fTFr1 = fTFr1;
199 ((AliTRDfeeParam &) p).fTFr2 = fTFr2;
200 ((AliTRDfeeParam &) p).fTFc1 = fTFc1;
201 ((AliTRDfeeParam &) p).fTFc2 = fTFc2;
202 ((AliTRDfeeParam &) p).fEBsglIndThr = fEBsglIndThr;
203 ((AliTRDfeeParam &) p).fEBsumIndThr = fEBsumIndThr;
204 ((AliTRDfeeParam &) p).fEBindLUT = fEBindLUT;
205 ((AliTRDfeeParam &) p).fEBignoreNeighbour = fEBignoreNeighbour;
206 ((AliTRDfeeParam &) p).fRAWversion = fRAWversion;
207 ((AliTRDfeeParam &) p).fRAWstoreRaw = fRAWstoreRaw;
208
209 TObject::Copy(p);
210
211}
212
213//_____________________________________________________________________________
214Int_t AliTRDfeeParam::GetPadRowFromMCM(Int_t irob, Int_t imcm) const
215{
216 //
217 // Return on which pad row this mcm sits
218 //
219
220 return fgkNmcmRobInRow*(irob/2) + imcm/fgkNmcmRobInCol;
221
222}
223
224//_____________________________________________________________________________
225Int_t AliTRDfeeParam::GetPadColFromADC(Int_t irob, Int_t imcm, Int_t iadc) const
226{
227 //
228 // Return which pad is connected to this adc channel.
229 //
230 // Return virtual pad number even if ADC is outside chamber
231 // to keep compatibility of data processing at the edge MCM.
232 // User has to check that this is in the chamber if it is essential.
233 // Return -100 if iadc is invalid.
234 //
235 // Caution: ADC ordering in the online data is opposite to the pad column ordering.
236 // And it is not one-by-one correspondence. Precise drawing can be found in:
237 // http://wiki.kip.uni-heidelberg.de/ti/TRD/index.php/Image:ROB_MCM_numbering.pdf
238 //
239
240 if (iadc < 0 || iadc > 19 ) return -100;
241 Int_t mcmcol = imcm%fgkNmcmRobInCol + GetRobSide(irob)*fgkNmcmRobInCol; // MCM column number on ROC [0..7]
242 Int_t padcol = mcmcol*fgkNcolMcm + fgkNcolMcm + 1 - iadc;
243 // if( padcol < 0 || padcol >= fgkNcol ) return -1; // thisi s commented because of reson above KO
244
245 return padcol;
246
247}
248
249//_____________________________________________________________________________
250Int_t AliTRDfeeParam::GetMCMfromPad(Int_t irow, Int_t icol) const
251{
252 //
253 // Return on which MCM this pad is directry connected.
254 // Return -1 for error.
255 //
256
257 if ( irow < 0 || icol < 0 || irow > fgkNrowC1 || icol > fgkNcol ) return -1;
258
259 return (icol%(fgkNcol/2))/fgkNcolMcm + fgkNmcmRobInCol*(irow%fgkNmcmRobInRow);
260
261}
262
263//_____________________________________________________________________________
264Int_t AliTRDfeeParam::GetROBfromPad(Int_t irow, Int_t icol) const
265{
266 //
267 // Return on which rob this pad is
268 //
269
270 return (irow/fgkNmcmRobInRow)*2 + GetColSide(icol);
271
272}
273
274//_____________________________________________________________________________
275Int_t AliTRDfeeParam::GetRobSide(Int_t irob) const
276{
277 //
278 // Return on which side this rob sits (A side = 0, B side = 1)
279 //
280
281 if ( irob < 0 || irob >= fgkNrobC1 ) return -1;
282
283 return irob%2;
284
285}
286
287//_____________________________________________________________________________
288Int_t AliTRDfeeParam::GetColSide(Int_t icol) const
289{
290 //
291 // Return on which side this column sits (A side = 0, B side = 1)
292 //
293
294 if ( icol < 0 || icol >= fgkNcol ) return -1;
295
296 return icol/(fgkNcol/2);
297
298}
299
300//_____________________________________________________________________________
301//void AliTRDfeeParam::GetFilterParam( Float_t &r1, Float_t &r2, Float_t &c1
302// , Float_t &c2, Float_t &ped ) const
303//{
304 //
305 // Return current filter parameter
306 //
307
308 // r1 = fR1;
309 //r2 = fR2;
310 //c1 = fC1;
311 //c2 = fC2;
312 //ped = fPedestal;
313//};
314
315//_____________________________________________________________________________
316void AliTRDfeeParam::SetEBsglIndThr(Int_t val)
317{
318 //
319 // Set Event Buffer Sngle Indicator Threshold (EBIS in TRAP conf).
320 // Timebin is indicated if ADC value >= val.
321 //
322
323 if( val >= 0 && val <= 1023 ) {
324 fEBsglIndThr = val;
325 } else {
326 AliError(Form("EBsglIndThr value %d is out of range, keep previously set value (%d).",
327 val, fEBsglIndThr));
328 }
329
330}
331
332//_____________________________________________________________________________
333void AliTRDfeeParam::SetEBsumIndThr(Int_t val)
334{
335 //
336 // Set Event Buffer Sum Indicator Threshold (EBIT in TRAP conf).
337 // Timebin is indicated if ADC sum value >= val.
338 //
339
340 if( val >= 0 && val <= 4095 ) {
341 fEBsumIndThr = val;
342 }
343 else {
344 AliError(Form("EBsumIndThr value %d is out of range, keep previously set value (%d).",
345 val, fEBsumIndThr));
346 }
347
348}
349
350//_____________________________________________________________________________
351void AliTRDfeeParam::SetEBindLUT(Int_t val)
352{
353 //
354 // Set Event Buffer Indicator Look-Up Table (EBIL in TRAP conf).
355 // 8 bits value forms lookup table for combination of three criterions.
356 //
357
358 if( val >= 0 && val <= 255 ) {
359 fEBindLUT = val;
360 }
361 else {
362 AliError(Form("EBindLUT value %d is out of range, keep previously set value (%d).",
363 val, fEBindLUT));
364 }
365
366}
367
368//_____________________________________________________________________________
369void AliTRDfeeParam::SetEBignoreNeighbour(Int_t val)
370{
371 //
372 // Set Event Buffer Indicator Neighbor Sensitivity. (EBIN in TRAP conf).
373 // If 0, take account of neigbor's values.
374 //
375
376 if( val >= 0 && val <= 1 ) {
377 fEBignoreNeighbour = val;
378 }
379 else {
380 AliError(Form("EBignoreNeighbour value %d is out of range, keep previously set value (%d).",
381 val, fEBignoreNeighbour));
382 }
383}
384
385//_____________________________________________________________________________
386void AliTRDfeeParam::SetRAWversion( Int_t rawver )
387{
388 //
389 // Set raw data version (major number only)
390 // Maximum available number is preset in fgkMaxRAWversion
391 //
392
393 if( rawver >= 0 && rawver <= fgkMaxRAWversion ) {
394 fRAWversion = rawver;
395 }
396 else {
397 AliError(Form("Raw version is out of range: %d",rawver));
398 }
399
400}
401
402//_____________________________________________________________________________
403void AliTRDfeeParam::SetXenon()
404{
405 //
406 // Sets the filter parameters for the Xenon gas mixture
407 //
408
409 fTFnExp = 1;
410 fTFr1 = 1.1563;
411 fTFr2 = 0.1299;
412 fTFc1 = 0.0657;
413 fTFc2 = 0.0000;
414
415}
416
417//_____________________________________________________________________________
418void AliTRDfeeParam::SetArgon()
419{
420 //
421 // Sets the filter parameters for the Argon gas mixture
422 //
423
424 fTFnExp = 2;
425 fTFr1 = 6.0;
426 fTFr2 = 0.62;
427 fTFc1 = 0.0087;
428 fTFc2 = 0.07;
429
430}