1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
18 //-----------------------------------------------------------------------------
19 /// \class AliMUONTriggerCircuit
20 /// Contains as data members the Y positions of the X declusturized strips and
21 /// the X positions of the (doubled or not) Y strips.
22 /// This is used to associate the global positions to the fired strips of the
23 /// local trigger output (see AliMUONTrackReconstructor::MakeTriggerTrack)
25 /// \author Philippe Crochet (LPCCFd)
26 //-----------------------------------------------------------------------------
28 #include "AliMUONTriggerCircuit.h"
29 #include "AliMUONConstants.h"
30 #include "AliMUONGeometryTransformer.h"
32 #include "AliMpTrigger.h"
33 #include "AliMpSlat.h"
35 #include "AliMpSegmentation.h"
36 #include "AliMpVSegmentation.h"
37 #include "AliMpCathodType.h"
38 #include "AliMpDDLStore.h"
39 #include "AliMpLocalBoard.h"
40 #include "AliMpConstants.h"
47 #include <Riostream.h>
50 ClassImp(AliMUONTriggerCircuit)
53 //----------------------------------------------------------------------
54 AliMUONTriggerCircuit::AliMUONTriggerCircuit(const AliMUONGeometryTransformer* transformer)
56 fkTransformer(transformer),
59 fCurrentLocalBoard(0x0)
63 for (Int_t i = 1; i < AliMpConstants::NofLocalBoards()+1; ++i) {
69 for (Int_t i = 1; i < AliMpConstants::NofLocalBoards()+1; ++i) { // board begins at 1
71 AliMpLocalBoard* localBoard = AliMpDDLStore::Instance()->GetLocalBoard(i);
75 AliError(Form("Did not get localboard %d",i));
86 //----------------------------------------------------------------------
87 AliMUONTriggerCircuit::~AliMUONTriggerCircuit()
90 for (Int_t i = 1; i < AliMpConstants::NofLocalBoards()+1; ++i) {
98 //----------------------------------------------------------------------
99 AliMUONTriggerCircuit::AliMUONTriggerCircuit(const AliMUONTriggerCircuit& circuit)
101 fkTransformer(circuit.fkTransformer), // do not copy, just pointed to
102 fkCurrentSeg(circuit.fkCurrentSeg),
103 fCurrentDetElem(circuit.fCurrentDetElem),
104 fCurrentLocalBoard(circuit.fCurrentLocalBoard)
108 for (Int_t i = 1; i < AliMpConstants::NofLocalBoards()+1; ++i) {
109 fXpos11[i] = circuit.fXpos11[i];
110 fYpos11[i] = circuit.fYpos11[i];
111 fYpos21[i] = circuit.fYpos21[i];
115 //----------------------------------------------------------------------
116 AliMUONTriggerCircuit& AliMUONTriggerCircuit::operator=(const AliMUONTriggerCircuit& circuit)
118 /// Assignment operator
120 if (this == &circuit) return *this;
122 fkTransformer = circuit.fkTransformer;
123 fkCurrentSeg = circuit.fkCurrentSeg;
124 fCurrentDetElem = circuit.fCurrentDetElem;
125 fCurrentLocalBoard = circuit.fCurrentLocalBoard;
127 for (Int_t i = 1; i < AliMpConstants::NofLocalBoards()+1; ++i) {
128 fXpos11[i] = circuit.fXpos11[i];
129 fYpos11[i] = circuit.fYpos11[i];
130 fYpos21[i] = circuit.fYpos21[i];
137 //---------------------------------------------------------------------
138 void AliMUONTriggerCircuit::LoadYPos(AliMpLocalBoard* const localBoard)
140 /// fill fYpos11 and fYpos21 -> y position of X declusterized strips
142 fCurrentLocalBoard = localBoard->GetId();
146 Int_t zeroDown = localBoard->GetSwitch(AliMpLocalBoard::kZeroDown);
147 Int_t zeroUp = localBoard->GetSwitch(AliMpLocalBoard::kZeroUp);
149 Int_t iline = localBoard->GetPosition().GetFirst();
150 Int_t icol = localBoard->GetPosition().GetSecond();
151 if ( iline == 5 ) --icol;
155 fCurrentDetElem = AliMpDDLStore::Instance()->GetDEfromLocalBoard(fCurrentLocalBoard, ichamber);
157 fkCurrentSeg = AliMpSegmentation::Instance()
158 ->GetMpSegmentation(fCurrentDetElem, AliMp::GetCathodType(icathode));
160 Int_t iFirstStrip = FirstStrip(localBoard);
161 Int_t iLastStrip = iFirstStrip + 16;
162 Int_t iStripCircuit = 0;
164 FillXstrips(icol, iFirstStrip, iLastStrip,
165 iStripCircuit, fYpos11[fCurrentLocalBoard]);
169 fCurrentDetElem = AliMpDDLStore::Instance()->GetDEfromLocalBoard(fCurrentLocalBoard, ichamber);
171 fkCurrentSeg = AliMpSegmentation::Instance()
172 ->GetMpSegmentation(fCurrentDetElem, AliMp::GetCathodType(icathode));
174 // second plane middle part
175 Int_t iFirstStripMiddle = FirstStrip(localBoard);
176 Int_t iLastStripMiddle = iFirstStrip + 16;
179 FillXstrips(icol, iFirstStripMiddle, iLastStripMiddle,
180 iStripCircuit, fYpos21[fCurrentLocalBoard]);
182 // second plane upper part
183 if (zeroUp == 0) { // something up
188 // check if we need to move to another detElemId
189 AliMpPad pad = fkCurrentSeg->PadByIndices(AliMpIntPair(icol-1,iLastStripMiddle+1),kFALSE);
191 if (pad.IsValid()) { // upper strips within same detElemId
192 iFirstStripUp = iLastStripMiddle;
193 iLastStripUp = iFirstStripUp + 8;
195 } else { // upper strips in another detElemId
196 fCurrentDetElem = AliMpDDLStore::Instance()->
197 GetNextDEfromLocalBoard(fCurrentLocalBoard, ichamber);
199 fkCurrentSeg = AliMpSegmentation::Instance()
200 ->GetMpSegmentation(fCurrentDetElem, AliMp::GetCathodType(icathode));
203 iLastStripUp = iFirstStripUp + 8;
204 if (iline == 4) icolUp = icol - 1; // special case
205 if (iline == 5) icolUp = icol + 1; // special case
209 FillXstrips(icolUp, iFirstStripUp, iLastStripUp,
210 iStripCircuit, fYpos21[fCurrentLocalBoard]);
212 // fill strip between middle and upper part
213 fYpos21[fCurrentLocalBoard][47] = (fYpos21[fCurrentLocalBoard][46] +
214 fYpos21[fCurrentLocalBoard][48])/2.;
215 } // end of something up
217 // restore current detElemId & segmentation
218 fCurrentDetElem = AliMpDDLStore::Instance()->GetDEfromLocalBoard(fCurrentLocalBoard, ichamber);
219 fkCurrentSeg = AliMpSegmentation::Instance()
220 ->GetMpSegmentation(fCurrentDetElem, AliMp::GetCathodType(icathode));
222 // second plane lower part
223 if (zeroDown == 0) { // something down
228 // check if we need to move to another detElemId
229 AliMpPad pad = fkCurrentSeg->PadByIndices(AliMpIntPair(icol-1,iFirstStripMiddle-1),kFALSE);
230 if (pad.IsValid()) { // lower strips within same detElemId
231 iFirstStripDo = iFirstStripMiddle - 8;
232 iLastStripDo = iFirstStripDo + 8;
234 } else { // lower strips in another detElemId
235 fCurrentDetElem = AliMpDDLStore::Instance()
236 ->GetPreviousDEfromLocalBoard(fCurrentLocalBoard, ichamber);
238 fkCurrentSeg = AliMpSegmentation::Instance()
239 ->GetMpSegmentation(fCurrentDetElem, AliMp::GetCathodType(icathode));
241 // get iFirstStrip in this module
242 const AliMpTrigger* t = AliMpSegmentation::Instance()->GetTrigger(fkCurrentSeg);
243 const AliMpSlat* slat = t->GetLayer(0);
245 if (iline == 5) icolDo = icol + 1; // special case
246 if (iline == 6) icolDo = icol - 1; // special case
248 const AliMpPCB* pcb = slat->GetPCB(icolDo-1);
249 iFirstStripDo = (pcb->Iymax() + 1) - 8;
250 iLastStripDo = iFirstStripDo + 8;
254 FillXstrips(icolDo, iFirstStripDo, iLastStripDo,
255 iStripCircuit, fYpos21[fCurrentLocalBoard]);
257 // fill strip between middle and upper part
258 fYpos21[fCurrentLocalBoard][15] = (fYpos21[fCurrentLocalBoard][14] +
259 fYpos21[fCurrentLocalBoard][16])/2.;
260 } // end of something down
264 //----------------------------------------------------------------------
265 void AliMUONTriggerCircuit::FillXstrips(const Int_t icol,
266 const Int_t iFirstStrip, const Int_t iLastStrip,
267 Int_t liStripCircuit, TArrayF& ypos)
270 Double_t xyGlobal[2] = {0.};
271 for (Int_t istrip = iFirstStrip; istrip < iLastStrip; ++istrip) {
273 AliMpPad pad = fkCurrentSeg->PadByIndices(AliMpIntPair(icol-1,istrip),kTRUE);
274 if ( !pad.IsValid() ) {
275 StdoutToAliError(cout << "Pad not found in seg " << endl;
276 fkCurrentSeg->Print();
277 cout << " ix,iy=" << icol-1 << "," << istrip << endl;
280 Float_t yDim = pad.Dimensions().Y(); // half size!
282 XYGlobal(pad,xyGlobal);
284 ypos[2*liStripCircuit] = xyGlobal[1];
285 if (istrip != (iLastStrip - 1)) ypos[2*liStripCircuit+1] = xyGlobal[1] + yDim;
291 //----------------------------------------------------------------------
292 void AliMUONTriggerCircuit::LoadXPos(AliMpLocalBoard* const localBoard)
294 /// fill fXpos11 -> x position of Y strips for the first plane only
295 /// fXpos11 contains the x position of Y strip for the current circuit
296 /// taking into account whether or nor not part(s) of the circuit
297 /// (middle, up or down) has(have) 16 strips (handdled by means of switchs)
299 fCurrentLocalBoard = localBoard->GetId();
304 Int_t x2u = localBoard->GetSwitch(AliMpLocalBoard::kX2u);
305 Int_t x2m = localBoard->GetSwitch(AliMpLocalBoard::kX2m);
306 Int_t x2d = localBoard->GetSwitch(AliMpLocalBoard::kX2d);
307 Int_t zeroAllYLSB = localBoard->GetSwitch(AliMpLocalBoard::kZeroAllYLSB);
309 Int_t iStripCircuit = 0;
310 Int_t iFirstStrip = 0;
311 Int_t iLastStrip = 0;
312 Bool_t doubling = kFALSE;
314 Int_t iline = localBoard->GetPosition().GetFirst();
315 Int_t icol = localBoard->GetPosition().GetSecond();
316 if ( iline == 5 ) --icol;
318 fCurrentDetElem = AliMpDDLStore::Instance()->GetDEfromLocalBoard(fCurrentLocalBoard, ichamber);
320 fkCurrentSeg = AliMpSegmentation::Instance()
321 ->GetMpSegmentation(fCurrentDetElem, AliMp::GetCathodType(icathode));
323 // check if one needs a strip doubling or not
324 if ( (x2u || x2m || x2d ) && x2m ) doubling = kTRUE;
327 // check if one starts at strip = 0 or 8 (boards 26-29 and 143-146)
328 if (zeroAllYLSB) iStripCircuit = 8;
330 // get iFirstStrip in this module
331 const AliMpTrigger* t = AliMpSegmentation::Instance()->GetTrigger(fkCurrentSeg);
332 const AliMpSlat* slat = t->GetLayer(0);
334 const AliMpPCB* pcb = slat->GetPCB(icol-1);
335 iFirstStrip = pcb->Ixmin();
338 if (doubling || zeroAllYLSB == 1) iLastStrip = iFirstStrip + 8;
339 else iLastStrip = iFirstStrip + 16;
341 FillYstrips(iFirstStrip, iLastStrip, iStripCircuit, doubling);
344 //----------------------------------------------------------------------
345 void AliMUONTriggerCircuit::FillYstrips(const Int_t iFirstStrip, const Int_t iLastStrip,
346 Int_t liStripCircuit,
347 const Bool_t doubling)
350 Double_t xyGlobal[2] = {0.};
352 for (Int_t istrip = iFirstStrip; istrip < iLastStrip; ++istrip) {
354 AliMpPad pad = fkCurrentSeg->PadByIndices(AliMpIntPair(istrip,0),kTRUE);
356 if ( !pad.IsValid() )
358 StdoutToAliError(cout << "Pad not found in seg " << endl;
359 fkCurrentSeg->Print();
360 cout << " ix,iy=" << istrip << "," << 0 << endl;
363 Float_t xDim = pad.Dimensions().X(); // half size!
365 XYGlobal(pad,xyGlobal);
368 fXpos11[fCurrentLocalBoard].AddAt(xyGlobal[0], liStripCircuit);
369 } else if (doubling) {
371 fXpos11[fCurrentLocalBoard].AddAt(TMath::Sign(1.,xyGlobal[0]) *
372 (TMath::Abs(xyGlobal[0]) - xDim/2.), 2*liStripCircuit);
374 fXpos11[fCurrentLocalBoard].AddAt(TMath::Sign(1.,xyGlobal[0]) *
375 (TMath::Abs(fXpos11[fCurrentLocalBoard][2*liStripCircuit]) + xDim),
376 2*liStripCircuit + 1);
383 //----------------------------------------------------------------------
384 void AliMUONTriggerCircuit::XYGlobal(const AliMpPad& pad,
387 /// returns pad x & y positions and x & y pad dimensions in global coordinates
388 /// note: no need for transformation for pad dimensions
390 // get the pad position and dimensions
391 Double_t xl1 = pad.Position().X();
392 Double_t yl1 = pad.Position().Y();
395 // positions from local to global
396 fkTransformer->Local2Global(fCurrentDetElem, xl1, yl1, 0,
397 xyGlobal[0], xyGlobal[1], zg1);
401 //----------------------------------------------------------------------
402 //--- methods which return member data related info
403 //----------------------------------------------------------------------
404 Float_t AliMUONTriggerCircuit::GetY11Pos(Int_t localBoardId, Int_t istrip) const
406 /// returns Y position of X strip istrip in MC11
407 return fYpos11[localBoardId][istrip];
409 //----------------------------------------------------------------------
410 Float_t AliMUONTriggerCircuit::GetY21Pos(Int_t localBoardId, Int_t istrip) const
412 /// returns Y position of X strip istrip in MC21
413 return fYpos21[localBoardId][istrip];
415 //----------------------------------------------------------------------
416 Float_t AliMUONTriggerCircuit::GetX11Pos(Int_t localBoardId, Int_t istrip) const
418 /// returns X position of Y strip istrip in MC11
419 return fXpos11[localBoardId][istrip];
422 //----------------------------------------------------------------------
423 Int_t AliMUONTriggerCircuit::FirstStrip(AliMpLocalBoard* localBoard)
425 /// returns the first strip from mapping for board boardName
426 /// take care of special case for boards RC1L6B12 & LC1L6B12
427 Int_t iFirstStrip = -1;
428 Int_t boardNumber = atoi(localBoard->GetName()+6);
430 Int_t iline = localBoard->GetPosition().GetFirst();
431 Int_t icol = localBoard->GetPosition().GetSecond();
432 if ( iline == 5 ) --icol;
449 if (icol == 1 && iline == 6) iFirstStrip += 16; // special case
453 //----------------------------------------------------------------------
454 Float_t AliMUONTriggerCircuit::PtCal(Int_t localBoardId, Int_t istripX, Int_t idev, Int_t istripY) const{
455 /// returns calculated pt for circuit/istripX/idev/istripY according
456 /// to the formula of the TRD. Note : idev (input) is in [0+30]
458 // Int_t jdev = idev - 15; // jdev in [-15+15]
459 Int_t istripX2=istripX+idev+1; // find istripX2 using istripX and idev
461 Float_t yPosX1=fYpos11[localBoardId][istripX];
462 Float_t yPosX2=fYpos21[localBoardId][istripX2];
463 Float_t xPosY1=fXpos11[localBoardId][istripY];
465 // Z distance between IP and center of dipole
466 Float_t zf= TMath::Abs(0.5 *(AliMUONConstants::CoilZ() + AliMUONConstants::YokeZ()));
467 Float_t z1=AliMUONConstants::DefaultChamberZ(10);
468 Float_t z2=AliMUONConstants::DefaultChamberZ(12);
469 Float_t thetaDev=(1./zf)*(yPosX1*z2-yPosX2*z1)/(z2-z1);
470 Float_t xf=xPosY1*zf/z1;
471 Float_t yf=yPosX2-((yPosX2-yPosX1)*(z2-zf))/(z2-z1);
472 return (3.*0.3/TMath::Abs(thetaDev)) * TMath::Sqrt(xf*xf+yf*yf)/zf;