]> git.uio.no Git - u/mrichter/AliRoot.git/blame - MUON/AliMUONTriggerCircuit.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Removed unused include
[u/mrichter/AliRoot.git] / MUON / AliMUONTriggerCircuit.cxx
CommitLineData
8ff54bb6 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/* $Id$ */
17
18/// \class AliMUONTriggerCircuit
19/// Contains as data members the Y positions of the X declusturized strips and
20/// the X positions of the (doubled or not) Y strips.
21/// This is used to associate the global positions to the fired strips of the
22/// local trigger output (see AliMUONTrackReconstructor::MakeTriggerTrack)
23
24#include "AliMUONTriggerCircuit.h"
25#include "AliMUON.h"
26#include "AliMUONConstants.h"
27#include "AliMUONLocalTriggerBoard.h"
28#include "AliMUONTriggerCrateStore.h"
29#include "AliMUONTriggerCrate.h"
30#include "AliMUONGeometryTransformer.h"
31
32#include "AliMpTriggerSegmentation.h"
33#include "AliMpTrigger.h"
34#include "AliMpSlat.h"
35#include "AliMpPCB.h"
36#include "AliMpSegmentation.h"
37#include "AliMpVSegmentation.h"
38
39#include "AliRun.h"
40#include "AliLog.h"
41
42#include <TMath.h>
43
44/// \cond CLASSIMP
45ClassImp(AliMUONTriggerCircuit)
46/// \endcond
47
48//----------------------------------------------------------------------
49AliMUONTriggerCircuit::AliMUONTriggerCircuit()
50: TObject(),
51 fILocalBoard(0),
52 fTransformer(0x0)
53{
54/// Constructor
55
56 Int_t i;
57 for (i=0; i<16; i++) { fXpos11[i]=0.; }
58 for (i=0; i<31; i++) { fYpos11[i]=0.; }
59 for (i=0; i<63; i++) { fYpos21[i]=0.; }
60}
61
62//----------------------------------------------------------------------
63AliMUONTriggerCircuit::~AliMUONTriggerCircuit()
64{
65/// Destructor
66}
67
68//----------------------------------------------------------------------
69AliMUONTriggerCircuit::AliMUONTriggerCircuit(const AliMUONTriggerCircuit& circuit)
70 : TObject(circuit),
1cb43687 71 fILocalBoard(circuit.fILocalBoard),
72 fTransformer(0x0)
8ff54bb6 73{
74/// Copy constructor
75
76 for (Int_t i = 0; i < 16; ++i)
77 fXpos11[i] = circuit.fXpos11[i];
78
79 for (Int_t i = 0; i < 31; ++i)
80 fYpos11[i] = circuit.fYpos11[i];
81
82 for (Int_t i = 0; i < 63; ++i)
83 fYpos21[i] = circuit.fYpos21[i];
84
85}
86//----------------------------------------------------------------------
87AliMUONTriggerCircuit& AliMUONTriggerCircuit::operator=(const AliMUONTriggerCircuit& circuit)
88{
89/// Assignment operator
90
91 if (this == &circuit) return *this;
92
93 fILocalBoard = circuit.fILocalBoard;
94
95 for (Int_t i = 0; i < 16; ++i)
96 fXpos11[i] = circuit.fXpos11[i];
97
98 for (Int_t i = 0; i < 31; ++i)
99 fYpos11[i] = circuit.fYpos11[i];
100
101 for (Int_t i = 0; i < 63; ++i)
102 fYpos21[i] = circuit.fYpos21[i];
103
104 return *this;
105
106}
107//----------------------------------------------------------------------
108void AliMUONTriggerCircuit::Init(Int_t iCircuit, const AliMUONTriggerCrateStore& crates)
109{
110/// initialize circuit characteristics
111 fILocalBoard=iCircuit+1;//AliMUONTriggerConstants::CircuitId(iCircuit);
112
113 LoadXPos(crates);
114 LoadYPos(crates);
115
116}
117
118//---------------------------------------------------------------------
119void AliMUONTriggerCircuit::LoadYPos(const AliMUONTriggerCrateStore& crates)
120{
121/// fill fYpos11 and fYpos21 -> y position of X declusterized strips
122
123 const AliMUONLocalTriggerBoard* localBoard = crates.LocalBoard(fILocalBoard);
124
125 if (!localBoard)
126 {
127 AliError(Form("Did not get localboard %d",fILocalBoard));
128 return;
129 }
130 StdoutToAliDebug(1,localBoard->Print("CONF"););
131
132 Int_t ichamber = 0;
133 Int_t icathode = 0;
134
135 const AliMpVSegmentation* seg;
136
137 Int_t zeroDown = localBoard->GetSwitch(AliMUONLocalTriggerBoard::kZeroDown);
138 Int_t zeroUp = localBoard->GetSwitch(AliMUONLocalTriggerBoard::kZeroUp);
139
140 //--- first plane
141 ichamber = 10;
142
143 char side;
144 Int_t iline, icol;
145 DecodeBoardName(localBoard->GetName(),side,iline,icol);
146
147 Int_t detElemId = DetElemId(ichamber,side,iline);
148 seg = AliMpSegmentation::Instance()->GetMpSegmentation(detElemId, icathode);
149
150 Int_t iFirstStrip = FirstStrip(localBoard->GetName());
151 Int_t iLastStrip = iFirstStrip + 16;
152 Int_t iStripCircuit = 0;
153 FillXstrips(seg,detElemId,icol,
154 iFirstStrip,iLastStrip,iStripCircuit,fYpos11);
155
156 //--- second plane
157 ichamber = 12;
158
159 detElemId = DetElemId(ichamber,side,iline);
160 seg = AliMpSegmentation::Instance()->GetMpSegmentation(detElemId, icathode);
161
162 // second plane middle part
163 Int_t iFirstStripMiddle = FirstStrip(localBoard->GetName());
164 Int_t iLastStripMiddle = iFirstStrip + 16;
165 iStripCircuit=8;
166 FillXstrips(seg,detElemId,icol,
167 iFirstStripMiddle,iLastStripMiddle,iStripCircuit,fYpos21);
168
169 // second plane upper part
170 if (zeroUp == 0) { // something up
171 Int_t iFirstStripUp;
172 Int_t iLastStripUp;
173 Int_t icolUp=icol;
174 // check if we need to move to another detElemId
175 AliMpPad pad = seg->PadByIndices(AliMpIntPair(icol-1,iLastStripMiddle+1),kFALSE);
176 if (pad.IsValid()) { // upper strips within same detElemId
177 iFirstStripUp = iLastStripMiddle;
178 iLastStripUp = iFirstStripUp + 8;
179 // icolUp = icol;
180 } else { // upper strips in another detElemId
181 detElemId = DetElemId(ichamber,side,iline+1); // get detElemId
182 seg = AliMpSegmentation::Instance()->GetMpSegmentation(detElemId, icathode);
183
184 iFirstStripUp = 0;
185 iLastStripUp = iFirstStripUp + 8;
186 if (iline == 4) icolUp = icol - 1; // special case
187 if (iline == 5) icolUp = icol + 1; // special case
188 }
189
190 iStripCircuit=24;
191 FillXstrips(seg,detElemId,icolUp,
192 iFirstStripUp,iLastStripUp,iStripCircuit,fYpos21);
193
194 // fill strip between middle and upper part
195 fYpos21[47]=(fYpos21[46]+fYpos21[48])/2.;
196 } // end of something up
197
198 // restore current detElemId & segmentation
199 detElemId = DetElemId(ichamber,side,iline);
200 seg = AliMpSegmentation::Instance()->GetMpSegmentation(detElemId, icathode);
201
202 // second plane lower part
203 if (zeroDown == 0) { // something down
204 Int_t iFirstStripDo;
205 Int_t iLastStripDo;
206 Int_t icolDo=icol;
207
208 // check if we need to move to another detElemId
209 AliMpPad pad = seg->PadByIndices(AliMpIntPair(icol-1,iFirstStripMiddle-1),kFALSE);
210 if (pad.IsValid()) { // lower strips within same detElemId
211 iFirstStripDo = iFirstStripMiddle - 8;
212 iLastStripDo = iFirstStripDo + 8;
213 // icolDo = icol;
214 } else { // lower strips in another detElemId
215 detElemId = DetElemId(ichamber,side,iline-1); // get detElemId
216 seg = AliMpSegmentation::Instance()->GetMpSegmentation(detElemId, icathode);
217
218 // get iFirstStrip in this module
219 const AliMpTriggerSegmentation* trig = (AliMpTriggerSegmentation*)(seg);
220 const AliMpTrigger* t = trig->Slat();
221 const AliMpSlat* slat = t->GetLayer(0);
222 if (iline == 5) icolDo = icol + 1; // special case
223 if (iline == 6) icolDo = icol - 1; // special case
224 const AliMpPCB* pcb = slat->GetPCB(icolDo-1);
225 iFirstStripDo = (pcb->Iymax() + 1) - 8;
226 iLastStripDo = iFirstStripDo + 8;
227 }
228
229 iStripCircuit=0;
230 FillXstrips(seg,detElemId,icolDo,
231 iFirstStripDo,iLastStripDo,iStripCircuit,fYpos21);
232
233 // fill strip between middle and upper part
234 fYpos21[15]=(fYpos21[14]+fYpos21[16])/2.;
235 } // end of something down
236
237}
238
239//----------------------------------------------------------------------
240void AliMUONTriggerCircuit::LoadXPos(const AliMUONTriggerCrateStore& crates)
241{
242/// fill fXpos11 -> x position of Y strips for the first plane only
243/// fXpos11 contains the x position of Y strip for the current circuit
244/// taking into account whether or nor not part(s) of the circuit
245/// (middle, up or down) has(have) 16 strips (handdled by means of switchs)
246
247 const AliMUONLocalTriggerBoard* localBoard = crates.LocalBoard(fILocalBoard);
248
249 if (!localBoard)
250 {
251 AliError(Form("Did not get localboard %d",fILocalBoard));
252 return;
253 }
254 StdoutToAliDebug(1,localBoard->Print("CONF"););
255
256 Int_t ichamber = 10;
257 Int_t icathode = 1;
258
259 Int_t x2u = localBoard->GetSwitch(AliMUONLocalTriggerBoard::kX2u);
260 Int_t x2m = localBoard->GetSwitch(AliMUONLocalTriggerBoard::kX2m);
261 Int_t x2d = localBoard->GetSwitch(AliMUONLocalTriggerBoard::kX2d);
262 Int_t zeroAllYLSB = localBoard->GetSwitch(AliMUONLocalTriggerBoard::kZeroAllYLSB);
263 Int_t iStripCircuit = 0;
264 Int_t iFirstStrip = 0;
265 Int_t iLastStrip = 0;
266 Bool_t doubling = kFALSE;
267
268 const AliMpVSegmentation* seg;
269
270 char side;
271 Int_t iline, icol;
272
273 DecodeBoardName(localBoard->GetName(),side,iline,icol);
274
275 Int_t detElemId=DetElemId(ichamber,side,iline); // get detElemId
276 seg = AliMpSegmentation::Instance()->GetMpSegmentation(detElemId, icathode);
277
278 // check if one needs a strip doubling or not
279 if ( (x2u == 1 || x2m == 1 || x2d == 1) && x2m == 1) doubling = kTRUE;
280
281 // check if one starts at strip = 0 or 8 (boards 26-29 and 143-146)
282 if (zeroAllYLSB == 1) iStripCircuit = 8;
283
284 // get iFirstStrip in this module
285 const AliMpTriggerSegmentation* trig = (AliMpTriggerSegmentation*)(seg);
286 const AliMpTrigger* t = trig->Slat();
287 const AliMpSlat* slat = t->GetLayer(0);
288
289 const AliMpPCB* pcb = slat->GetPCB(icol-1);
290 iFirstStrip = pcb->Ixmin();
291
292 if (doubling) iLastStrip = iFirstStrip + 8;
293 else iLastStrip = iFirstStrip + 16;
294
295 FillYstrips(seg,detElemId,
296 iFirstStrip,iLastStrip,iStripCircuit,doubling);
297}
298
299//----------------------------------------------------------------------
300void AliMUONTriggerCircuit::FillYstrips(
301 const AliMpVSegmentation* seg,
302 const Int_t detElemId,
303 const Int_t iFirstStrip, const Int_t iLastStrip, Int_t liStripCircuit,
304 const Bool_t doubling)
305{
306/// fill
307 Double_t xyGlobal[4]={0.,0.,0.,0.};
308 for (Int_t istrip=iFirstStrip; istrip<iLastStrip; istrip++) {
309 AliMpPad pad = seg->PadByIndices(AliMpIntPair(istrip,0),kTRUE);
310 if ( !pad.IsValid() )
311 {
312 StdoutToAliError(cout << "Pad not found in seg " << endl;
313 seg->Print();
314 cout << " ix,iy=" << istrip << "," << 0 << endl;
315 );
316 }
317 XYGlobal(detElemId,pad,xyGlobal);
318
319 if (!doubling) {
320 fXpos11[liStripCircuit]=xyGlobal[0];
321 } else if (doubling) {
322 fXpos11[2*liStripCircuit]=TMath::Sign(1.,xyGlobal[0]) *
323 (TMath::Abs(xyGlobal[0]) - xyGlobal[2]/2.);
324 fXpos11[2*liStripCircuit+1]=TMath::Sign(1.,xyGlobal[0]) *
325 (TMath::Abs(fXpos11[2*liStripCircuit]) + xyGlobal[2]);
326 }
327 liStripCircuit++;
328 }
329}
330
331//----------------------------------------------------------------------
332void AliMUONTriggerCircuit::FillXstrips(
333 const AliMpVSegmentation* seg,
334 const Int_t detElemId, const Int_t icol,
335 const Int_t iFirstStrip, const Int_t iLastStrip,
336 Int_t liStripCircuit, Float_t *tab)
337{
338/// fill
339 Double_t xyGlobal[4]={0.,0.,0.,0.};
340 for (Int_t istrip=iFirstStrip; istrip<iLastStrip; istrip++) {
341 AliMpPad pad = seg->PadByIndices(AliMpIntPair(icol-1,istrip),kTRUE);
342 if ( !pad.IsValid() )
343 {
344 StdoutToAliError(cout << "Pad not found in seg " << endl;
345 seg->Print();
346 cout << " ix,iy=" << icol-1 << "," << istrip << endl;
347 );
348 }
349
350 XYGlobal(detElemId,pad,xyGlobal);
351
352 tab[2*liStripCircuit]=xyGlobal[1];
353 if (istrip!=(iLastStrip-1)) tab[2*liStripCircuit+1]=xyGlobal[1]+xyGlobal[3];
354 liStripCircuit++;
355 }
356}
357
358//----------------------------------------------------------------------
359//--- methods which return member data related info
360//----------------------------------------------------------------------
361Float_t AliMUONTriggerCircuit::GetY11Pos(Int_t istrip) const {
362/// returns Y position of X strip istrip in MC11
363 return fYpos11[istrip];
364}
365//----------------------------------------------------------------------
366Float_t AliMUONTriggerCircuit::GetY21Pos(Int_t istrip) const {
367/// returns Y position of X strip istrip in MC21
368 return fYpos21[istrip];
369}
370//----------------------------------------------------------------------
371Float_t AliMUONTriggerCircuit::GetX11Pos(Int_t istrip) const {
372/// returns X position of Y strip istrip in MC11
373 return fXpos11[istrip];
374}
375//----------------------------------------------------------------------
376//--- end of methods which return member data related info
377//----------------------------------------------------------------------
378/* removed tmp
379void AliMUONTriggerCircuit::dump(const char* what, const Float_t* array, Int_t size)
380{
381 cout << what << " " << endl;
382 for ( Int_t i = 0; i < size; ++i )
383 {
384 cout << array[i] << " , ";
385 }
386 cout << endl;
387}
388
389void AliMUONTriggerCircuit::dump(const char* what, const Int_t* array, Int_t size)
390{
391 cout << what << " " << endl;
392 for ( Int_t i = 0; i < size; ++i )
393 {
394 cout << array[i] << " , ";
395 }
396 cout << endl;
397}
398
399//_____________________________________________________________________________
400void AliMUONTriggerCircuit::Print(Option_t* ) const
401{
402 cout << "IdCircuit " << fILocalBoard << " X2m,X2ud=" << fX2m << ","
403 << fX2ud;
404 for ( Int_t i = 0; i < 2; ++i )
405 {
406 cout << " OrMud[" << i << "]=" << fOrMud[i];
407 }
408 cout << endl;
409 dump("Xpos11",fXpos11,16);
410 dump("Ypos11",fYpos11,31);
411 dump("Ypos21",fYpos21,63);
412 for ( Int_t i = 0; i < 4; ++i )
413 {
414 char s[80];
415 sprintf(s,"Xcode[%d]",i);
416 dump(s,fXcode[i],32);
417 sprintf(s,"Ycode[%d]",i);
418 dump(s,fYcode[i],32);
419 }
420 // Int_t fILocalBoard; // circuit Id number
421 // Int_t fX2m; // internal info needed by TriggerDecision
422 // Int_t fX2ud; // internal info needed by TriggerDecision
423 // Int_t fOrMud[2]; // internal info needed by TriggerDecision
424 // Int_t fXcode[4][32]; // code of X strips
425 // Int_t fYcode[4][32]; // code of Y strips
426 // Float_t fXpos11[16]; // X position of Y strips in MC11
427 // Float_t fYpos11[31]; // Y position of X strips in MC11
428 // Float_t fYpos21[63]; // Y position of X strips in MC21
429
430}
431removed tmp*/
432
433//----------------------------------------------------------------------
434Int_t AliMUONTriggerCircuit::DetElemId(Int_t ichamber, char side, Int_t iline)
435{
436/// returns detection element Id for chamber iChamber, side side and line iline
437 Int_t itmp=0;
438 if ( side == 'R' ) { // right side
439 switch (iline) // (from 1 to 9, from bottom to top)
440 {
441 case 1:
442 itmp = 14;
443 break;
444 case 2:
445 itmp = 15;
446 break;
447 case 3:
448 itmp = 16;
449 break;
450 case 4:
451 itmp = 17;
452 break;
453 case 5:
454 itmp = 0;
455 break;
456 case 6:
457 itmp = 1;
458 break;
459 case 7:
460 itmp = 2;
461 break;
462 case 8:
463 itmp = 3;
464 break;
465 case 9:
466 itmp = 4;
467 break;
468 }
469 } else if ( side == 'L' ) { // left side
470 switch (iline) // (from 1 to 9, from bottom to top)
471 {
472 case 1:
473 itmp = 13;
474 break;
475 case 2:
476 itmp = 12;
477 break;
478 case 3:
479 itmp = 11;
480 break;
481 case 4:
482 itmp = 10;
483 break;
484 case 5:
485 itmp = 9;
486 break;
487 case 6:
488 itmp = 8;
489 break;
490 case 7:
491 itmp = 7;
492 break;
493 case 8:
494 itmp = 6;
495 break;
496 case 9:
497 itmp = 5;
498 break;
499 }
500 }
501 return ((ichamber+1)*100)+itmp;
502}
503
504//----------------------------------------------------------------------
505Int_t
506AliMUONTriggerCircuit::DetElemId(Int_t iChamber, const char* boardName)
507{
508/// returns detection element Id for chamber iChamber and board boardName
509 char side;
510 Int_t iline;
511 Int_t icol;
512
513 DecodeBoardName(boardName, side, iline, icol);
514
515 return DetElemId(iChamber,side,iline);
516}
517
518//----------------------------------------------------------------------
519void
520AliMUONTriggerCircuit::DecodeBoardName(const char* boardName,
521 char& side,
522 Int_t& iLine,
523 Int_t& iCol)
524{
525/// get side, line and col from board boardName
526/// note: icol = icol -1 for iline = 5 w.r.t other ilines
527 side = boardName[0];
528 iLine = boardName[4] - '0';
529 iCol = boardName[2] - '0';
530 if ( iLine == 5 ) --iCol;
531}
532
533//----------------------------------------------------------------------
534Int_t
535AliMUONTriggerCircuit::FirstStrip(const char* boardName)
536{
537/// returns the first strip from mapping for board boardName
538/// take care of special case for boards RC1L6B12 & LC1L6B12
539 Int_t iFirstStrip = -1;
540 Int_t boardNumber = atoi(boardName+6);
541 char side;
542 Int_t iline, icol;
543 DecodeBoardName(boardName,side,iline,icol);
544 switch (boardNumber)
545 {
546 case 12:
547 iFirstStrip = 0;
548 break;
549 case 34:
550 iFirstStrip = 16;
551 break;
552 case 56:
553 iFirstStrip = 32;
554 break;
555 case 78:
556 iFirstStrip = 48;
557 break;
558 }
559 if (icol == 1 && iline == 6) iFirstStrip = iFirstStrip + 16; // special case
560 return iFirstStrip;
561}
562
563//----------------------------------------------------------------------
564void AliMUONTriggerCircuit::XYGlobal(
565 Int_t detElemId, const AliMpPad& pad,
566 Double_t xyGlobal[4])
567{
568/// returns pad x & y positions and x & y pad dimensions in global coordinates
569/// note: no need for transformation for pad dimensions
570
571 // get the pad position and dimensions
572 Double_t xl1 = pad.Position().X();
573 Double_t yl1 = pad.Position().Y();
574 xyGlobal[2] = pad.Dimensions().X(); // half size!
575 xyGlobal[3] = pad.Dimensions().Y(); // half size!
576 Double_t zg1=0;
577
578 // positions from local to global
579 fTransformer->Local2Global(detElemId, xl1, yl1, 0,
580 xyGlobal[0], xyGlobal[1], zg1);
581}