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 **************************************************************************/
19 // Class AliMUONTriggerCircuit
21 // description of the global position of the trigger boards strip positions
22 // based on the "old" segmentation (to be replaced soon)
27 #include "AliMUONTriggerCircuit.h"
30 #include "AliMUONTriggerConstants.h"
31 #include "AliMUONGeometrySegmentation.h"
32 #include "AliMUONSegmentation.h"
33 #include "AliMUONConstants.h"
36 ClassImp(AliMUONTriggerCircuit)
38 //----------------------------------------------------------------------
39 AliMUONTriggerCircuit::AliMUONTriggerCircuit()
47 fOrMud[0]=fOrMud[1]=0;
50 for (Int_t j=0; j<32; j++) {
55 for (i=0; i<16; i++) { fXpos11[i]=0.; }
56 for (i=0; i<31; i++) { fYpos11[i]=0.; }
57 for (i=0; i<63; i++) { fYpos21[i]=0.; }
60 //----------------------------------------------------------------------
61 AliMUONTriggerCircuit::AliMUONTriggerCircuit(const AliMUONTriggerCircuit& theMUONTriggerCircuit)
62 : TObject(theMUONTriggerCircuit)
64 // Protected copy constructor
66 AliFatal("Not implemented.");
69 //----------------------------------------------------------------------
70 AliMUONTriggerCircuit &
71 AliMUONTriggerCircuit::operator=(const AliMUONTriggerCircuit& rhs)
73 // Protected assignement operator
75 if (this == &rhs) return *this;
77 AliFatal("Not implemented.");
82 //----------------------------------------------------------------------
83 void AliMUONTriggerCircuit::Init(Int_t iCircuit) {
84 // initialize circuit characteristics
85 fIdCircuit=AliMUONTriggerConstants::CircuitId(iCircuit);
96 //----------------------------------------------------------------------
97 Int_t AliMUONTriggerCircuit::CircuitNumber(Int_t idCircuit) const {
98 // returns circuit number iCircuit (0-234) corresponding to circuit idCircuit
100 for (Int_t i=0; i<234; i++) {
101 if (AliMUONTriggerConstants::CircuitId(i)==idCircuit) {
108 //----------------------------------------------------------------------
109 Int_t AliMUONTriggerCircuit::ModuleNumber(Int_t idModule) const {
110 // returns module number imod (from 0 to 63) corresponding to module idmodule
111 Int_t absidModule=TMath::Abs(idModule);
113 for (Int_t i=0; i<63; i++) {
114 if (AliMUONTriggerConstants::ModuleId(i)==absidModule) {
122 //----------------------------------------------------------------------
123 Int_t AliMUONTriggerCircuit::Module(Int_t idCircuit) const {
124 // returns ModuleId where Circuit idCircuit is sitting
125 return Int_t(idCircuit/10);
127 //----------------------------------------------------------------------
128 Int_t AliMUONTriggerCircuit::Position(Int_t idCircuit) const {
129 // returns position of idCircuit in correcponding Module
130 return TMath::Abs(idCircuit)-TMath::Abs(Module(idCircuit))*10;
133 //----------------------------------------------------------------------
134 void AliMUONTriggerCircuit::LoadX2() {
135 // initialize fX2m, fX2ud and fOrMud
137 Int_t idModule=Module(fIdCircuit); // corresponding module Id.
138 // and its number of X strips
139 Int_t nStrX=AliMUONTriggerConstants::NstripX(ModuleNumber(idModule));
140 // and its number of Y strips
141 Int_t nStrY=AliMUONTriggerConstants::NstripY(ModuleNumber(idModule));
142 Int_t iPosCircuit=Position(fIdCircuit); // position of circuit in module
144 // first step : look at lower part
145 if (iPosCircuit==1) { // need to scan lower module
146 if(idModule<91&&TMath::Abs(idModule)!=41&&idModule>-91) {
148 Int_t idModuleD=(TMath::Abs(idModule)+10)*(TMath::Abs(idModule)/idModule);
149 Int_t nStrD=AliMUONTriggerConstants::NstripY(ModuleNumber(idModuleD));
152 &&TMath::Abs(idModule)!=42&&TMath::Abs(idModule)!=52) {
153 if (nStrY==8) fX2m=1;
154 if (nStrD==8) fX2ud=1;
158 } else { // lower strips within same module
162 // second step : look at upper part
163 if ((iPosCircuit==1&&nStrX==16)||(iPosCircuit==2&&nStrX==32)||
164 (iPosCircuit==3&&nStrX==48)||(iPosCircuit==4&&nStrX==64)) {
165 if ((idModule>17||idModule<-17)&&TMath::Abs(idModule)!=61) {
167 Int_t idModuleU=(TMath::Abs(idModule)-10)*(TMath::Abs(idModule)/idModule);
168 Int_t nStrU=AliMUONTriggerConstants::NstripY(ModuleNumber(idModuleU));
171 &&TMath::Abs(idModule)!=62&&TMath::Abs(idModule)!=52) {
172 if (nStrY==8) fX2m=1;
173 if (nStrU==8) fX2ud=1;
177 } else { // upper strips within same module
182 //----------------------------------------------------------------------
183 void AliMUONTriggerCircuit::LoadXCode(){
184 // assign a Id. number to each X strip of current circuit
185 // Id.=(corresponding module Id.)*100+(Id. strip of module)
187 // first part : fill XMC11 XMC12 and strips 8 to 24 (middle) XMC21 XMC22
188 Int_t iStripCircMT1=0, iStripCircMT2=8;
189 Int_t idModule=Module(fIdCircuit); // corresponding module Id.
190 // and its number of strips
191 Int_t nStrX=AliMUONTriggerConstants::NstripX(ModuleNumber(idModule));
192 Int_t iPosCircuit=Position(fIdCircuit); // position of circuit in module
193 Int_t sign=TMath::Abs(idModule)/idModule; // left or right
196 for (istrip=(iPosCircuit-1)*16;
197 istrip<(iPosCircuit-1)*16+16; istrip++) {
199 fXcode[0][iStripCircMT1]=sign*(TMath::Abs(idModule)*100+istrip);
200 fXcode[1][iStripCircMT1]=sign*(TMath::Abs(idModule)*100+istrip);
201 fXcode[2][iStripCircMT2]=sign*(TMath::Abs(idModule)*100+istrip);
202 fXcode[3][iStripCircMT2]=sign*(TMath::Abs(idModule)*100+istrip);
208 // XMC21 XMC22 strips 0 to 7 and 24 to 31
209 Int_t idModuleD, idModuleU;
212 idModule=Module(fIdCircuit); // corresponding module Id.
213 // number of X strips
214 nStrX=AliMUONTriggerConstants::NstripX(ModuleNumber(idModule));
215 sign=TMath::Abs(idModule)/idModule;
217 // fill lower part (0 to 7)
218 if (iPosCircuit==1) { // need to scan lower module
219 if(idModule<91&&TMath::Abs(idModule)!=41&&idModule>-91) { // non-existing
220 idModuleD=sign*(TMath::Abs(idModule)+10); // lower module Id
221 // and its number of strips
222 nStrD=AliMUONTriggerConstants::NstripX(ModuleNumber(idModuleD));
225 for (istrip=nStrD-8; istrip<nStrD; istrip++) {
226 fXcode[2][iStripCircMT2]=sign*(TMath::Abs(idModuleD)*100+istrip);
227 fXcode[3][iStripCircMT2]=sign*(TMath::Abs(idModuleD)*100+istrip);
232 } else { // lower strips within same module
235 for (istrip=(iPosCircuit-1)*16-8;
236 istrip<(iPosCircuit-1)*16; istrip++) {
237 fXcode[2][iStripCircMT2]=sign*(TMath::Abs(idModule)*100+istrip);
238 fXcode[3][iStripCircMT2]=sign*(TMath::Abs(idModule)*100+istrip);
243 // fill upper part (24 to 31)
244 if ((iPosCircuit==1&&nStrX==16)||(iPosCircuit==2&&nStrX==32)||
245 (iPosCircuit==3&&nStrX==48)||(iPosCircuit==4&&nStrX==64)) {
246 if ((idModule>17||idModule<-17)&&TMath::Abs(idModule)!=61) {
247 idModuleU=sign*(TMath::Abs(idModule)-10); // upper module Id
248 // and its number of strips
249 nStrU=AliMUONTriggerConstants::NstripX(ModuleNumber(idModuleU));
252 for (istrip=0; istrip<8; istrip++) {
253 fXcode[2][iStripCircMT2]=sign*(TMath::Abs(idModuleU)*100+istrip);
254 fXcode[3][iStripCircMT2]=sign*(TMath::Abs(idModuleU)*100+istrip);
259 } else if ((iPosCircuit==1&&nStrX>16)||(iPosCircuit==2&&nStrX>32)||
260 (iPosCircuit==3&&nStrX>48)) { // upper strips within same mod.
263 for (istrip=(iPosCircuit-1)*16+16;
264 istrip<(iPosCircuit-1)*16+24; istrip++) {
265 fXcode[2][iStripCircMT2]=sign*(TMath::Abs(idModule)*100+istrip);
266 fXcode[3][iStripCircMT2]=sign*(TMath::Abs(idModule)*100+istrip);
272 //----------------------------------------------------------------------
273 void AliMUONTriggerCircuit::LoadYCode(){
274 // assign a Id. number to each Y strip of current circuit
275 // Id.=(corresponding module Id.)*100+(Id. strip of module)
276 // note : for Y plane fill only "central part" of circuit
277 // (upper and lower parts are filled in PreHandlingY of AliMUONTriggerDecision)
279 Int_t idModule=Module(fIdCircuit); // corresponding module Id.
280 // and its number of Y strips
281 Int_t nStrY=AliMUONTriggerConstants::NstripY(ModuleNumber(idModule));
282 Int_t sign=TMath::Abs(idModule)/idModule; // left or right
284 for (Int_t istrip=0; istrip<nStrY; istrip++) {
285 fYcode[0][istrip]=sign*(TMath::Abs(idModule)*100+istrip);
286 fYcode[1][istrip]=sign*(TMath::Abs(idModule)*100+istrip);
287 fYcode[2][istrip]=sign*(TMath::Abs(idModule)*100+istrip);
288 fYcode[3][istrip]=sign*(TMath::Abs(idModule)*100+istrip);
292 //----------------------------------------------------------------------
293 Float_t AliMUONTriggerCircuit::PtCal(Int_t istripX, Int_t idev, Int_t istripY){
294 // returns calculated pt for circuit/istripX/idev/istripY according
295 // to the formula of the TRD. Note : idev (input) is in [0+30]
297 // Int_t jdev = idev - 15; // jdev in [-15+15]
298 Int_t istripX2=istripX+idev+1; // find istripX2 using istripX and idev
300 Float_t yPosX1=fYpos11[istripX];
301 Float_t yPosX2=fYpos21[istripX2];
302 Float_t xPosY1=fXpos11[istripY];
305 Float_t z1=AliMUONConstants::DefaultChamberZ(10);
306 Float_t z2=AliMUONConstants::DefaultChamberZ(12);
307 Float_t thetaDev=(1./zf)*(yPosX1*z2-yPosX2*z1)/(z2-z1);
308 Float_t xf=xPosY1*zf/z1;
309 Float_t yf=yPosX2-((yPosX2-yPosX1)*(z2-zf))/(z2-z1);
310 return (3.*0.3/TMath::Abs(thetaDev)) * TMath::Sqrt(xf*xf+yf*yf)/zf;
312 //---------------------------------------------------------------------
313 //----------------------- New Segmentation ----------------------------
314 //---------------------------------------------------------------------
316 //---------------------------------------------------------------------
317 void AliMUONTriggerCircuit::LoadYPos2(){
318 // fill fYpos11 and fYpos21 -> y position of X declusterized strips
320 Int_t chamber, cathode;
321 Int_t code, idModule, idStrip, idSector;
322 Float_t x, y, z, width;
325 AliMUON *pMUON = (AliMUON*)gAlice->GetModule("MUON");
326 AliMUONGeometrySegmentation* segmentation;
332 = pMUON->GetSegmentation()->GetModuleSegmentation(chamber-1, cathode-1);
335 AliWarning("Segmentation not defined.");
339 for (istrip=0; istrip<16; istrip++) {
340 code=fXcode[0][istrip]; // decode current strip
341 idModule=Int_t(code/100); // corresponding module Id.
342 idDE = DetElemId(chamber, idModule);
343 idStrip=TMath::Abs(code-idModule*100); // corresp. strip number in module
344 idSector=segmentation->Sector(idDE, idModule, idStrip); // corresponding sector
345 width=segmentation->Dpy(idDE, idSector); // corresponding strip width
346 segmentation->GetPadC(idDE, idModule,idStrip,x,y,z); // get strip real position
349 if (istrip!=15) fYpos11[2*istrip+1]=y+width/2.;
356 = pMUON->GetSegmentation()->GetModuleSegmentation(chamber-1, cathode-1);
358 for (istrip=0; istrip<32; istrip++) {
359 code=fXcode[2][istrip]; // decode current strip
360 idModule=Int_t(code/100); // corresponding module Id.
361 idDE = DetElemId(chamber, idModule);
362 if (idModule == 0) continue;
363 idStrip=TMath::Abs(code-idModule*100); // corresp. strip number in module
364 idSector=segmentation->Sector(idDE, idModule, idStrip); // corresponding sector
365 width=segmentation->Dpy(idDE, idSector); // corresponding strip width
366 segmentation->GetPadC(idDE, idModule,idStrip,x,y,z); // get strip real position
368 // using idModule!=0 prevents to fill garbage in case of circuits
369 // in the first and last rows
372 if (istrip!=31) fYpos21[2*istrip+1]=y+width/2.;
377 //----------------------------------------------------------------------
378 void AliMUONTriggerCircuit::LoadXPos2(){
379 // fill fXpos11 -> x position of Y strips for the first plane only
380 // fXpos11 contains the x position of Y strip for the current circuit
381 // taking into account whether or nor not part(s) of the circuit
382 // (middle, up or down) has(have) 16 strips
389 AliMUON *pMUON = (AliMUON*)gAlice->GetModule("MUON");
390 AliMUONGeometrySegmentation* segmentation;
392 = pMUON->GetSegmentation()->GetModuleSegmentation(chamber-1, cathode-1);
395 AliWarning("Segmentation not defined.");
399 Int_t idModule=Module(fIdCircuit); // corresponding module Id.
400 // number of Y strips
401 idDE = DetElemId(chamber, idModule);
403 Int_t nStrY=AliMUONTriggerConstants::NstripY(ModuleNumber(idModule));
404 Int_t idSector=segmentation->Sector(idDE, idModule,0); // corresp. sector
405 Float_t width=segmentation->Dpx(idDE, idSector); // corresponding strip width
407 // first case : up middle and down parts have all 8 or 16 strip
408 if ((nStrY==16)||(nStrY==8&&fX2m==0&&fX2ud==0)) {
409 for (istrip=0; istrip<nStrY; istrip++) {
410 segmentation->GetPadC(idDE, idModule,istrip,x,y,z);
411 AliDebug(1,Form("idDE %d idModule %d istrip %d x,y,z=%e,%e,%e\n",
412 idDE,idModule,istrip,x,y,z));
415 // second case : mixing 8 and 16 strips within same circuit
417 for (istrip=0; istrip<nStrY; istrip++) {
418 if (nStrY!=8) { printf(" bug in LoadXpos \n");}
419 segmentation->GetPadC(idDE, idModule, istrip, x, y, z);
420 fXpos11[2*istrip]=x-width/4.;
421 fXpos11[2*istrip+1]=fXpos11[2*istrip]+width/2.;
426 //----------------------------------------------------------------------
427 //--- methods which return member data related info
428 //----------------------------------------------------------------------
429 Int_t AliMUONTriggerCircuit::GetIdCircuit() const {
430 // returns circuit Id
433 //----------------------------------------------------------------------
434 Int_t AliMUONTriggerCircuit::GetIdModule() const {
436 return Module(fIdCircuit);
438 //----------------------------------------------------------------------
439 Int_t AliMUONTriggerCircuit::GetNstripX() const {
440 // returns the number of X strips in the module where the circuit is sitting
441 return AliMUONTriggerConstants::NstripX(ModuleNumber(Module(fIdCircuit)));
443 //----------------------------------------------------------------------
444 Int_t AliMUONTriggerCircuit::GetNstripY() const {
445 // returns the number of Y strips in the module where the circuit is sitting
446 return AliMUONTriggerConstants::NstripY(ModuleNumber(Module(fIdCircuit)));
448 //----------------------------------------------------------------------
449 Int_t AliMUONTriggerCircuit::GetPosCircuit() const {
450 // returns the position of the circuit in its module
451 return Position(fIdCircuit);
453 //----------------------------------------------------------------------
454 Int_t AliMUONTriggerCircuit::GetIdCircuitD() const {
455 // returns the Id of the circuit down
456 Int_t idModule=Module(fIdCircuit);
457 Int_t idModuleD=(TMath::Abs(idModule)+10)*(TMath::Abs(idModule)/idModule);
458 return (TMath::Abs(idModuleD)*10+1)*(TMath::Abs(idModule)/idModule);
460 //----------------------------------------------------------------------
461 Int_t AliMUONTriggerCircuit::GetICircuitD() const {
462 // returns the number of the circuit down
463 Int_t idModule=Module(fIdCircuit);
464 Int_t idModuleD=(TMath::Abs(idModule)+10)*(TMath::Abs(idModule)/idModule);
466 (TMath::Abs(idModuleD)*10+1)*(TMath::Abs(idModule)/idModule);
467 return CircuitNumber(idCircuitD);
469 //----------------------------------------------------------------------
470 Int_t AliMUONTriggerCircuit::GetIdCircuitU() const {
471 // returns the Id of the circuit up
472 Int_t idModule=Module(fIdCircuit);
473 Int_t idModuleU=(TMath::Abs(idModule)-10)*(TMath::Abs(idModule)/idModule);
474 return (TMath::Abs(idModuleU)*10+1)*(TMath::Abs(idModule)/idModule);
476 //----------------------------------------------------------------------
477 Int_t AliMUONTriggerCircuit::GetICircuitU() const {
478 // returns the number of the circuit up
479 Int_t idModule=Module(fIdCircuit);
480 Int_t idModuleU=(TMath::Abs(idModule)-10)*(TMath::Abs(idModule)/idModule);
482 (TMath::Abs(idModuleU)*10+1)*(TMath::Abs(idModule)/idModule);
483 return CircuitNumber(idCircuitU);
485 //----------------------------------------------------------------------
486 Int_t AliMUONTriggerCircuit::GetX2m() const {
490 //----------------------------------------------------------------------
491 Int_t AliMUONTriggerCircuit::GetX2ud() const {
495 //----------------------------------------------------------------------
496 void AliMUONTriggerCircuit::GetOrMud(Int_t orMud[2]) const {
501 //----------------------------------------------------------------------
502 Int_t AliMUONTriggerCircuit::GetXcode(Int_t chamber, Int_t istrip) const {
503 // returns X code of circuit/chamber/istrip (warning : chamber in [0,3])
504 return fXcode[chamber][istrip];
506 //----------------------------------------------------------------------
507 Int_t AliMUONTriggerCircuit::GetYcode(Int_t chamber, Int_t istrip) const {
508 // returns Y code of circuit/chamber/istrip (warning : chamber in [0,3])
509 return fYcode[chamber][istrip];
511 //----------------------------------------------------------------------
512 Float_t AliMUONTriggerCircuit::GetY11Pos(Int_t istrip) const {
513 // returns Y position of X strip istrip in MC11
514 return fYpos11[istrip];
516 //----------------------------------------------------------------------
517 Float_t AliMUONTriggerCircuit::GetY21Pos(Int_t istrip) const {
518 // returns Y position of X strip istrip in MC21
519 return fYpos21[istrip];
521 //----------------------------------------------------------------------
522 Float_t AliMUONTriggerCircuit::GetX11Pos(Int_t istrip) const {
523 // returns X position of Y strip istrip in MC11
524 return fXpos11[istrip];
526 //----------------------------------------------------------------------
527 //--- end of methods which return member data related info
528 //----------------------------------------------------------------------
530 void dump(const char* what, const Float_t* array, Int_t size)
532 cout << what << " " << endl;
533 for ( Int_t i = 0; i < size; ++i )
535 cout << array[i] << " , ";
540 void dump(const char* what, const Int_t* array, Int_t size)
542 cout << what << " " << endl;
543 for ( Int_t i = 0; i < size; ++i )
545 cout << array[i] << " , ";
550 //_____________________________________________________________________________
552 AliMUONTriggerCircuit::Print(Option_t* ) const
554 // prints-out data member
555 cout << "IdCircuit " << fIdCircuit << " X2m,X2ud=" << fX2m << ","
557 for ( Int_t i = 0; i < 2; ++i )
559 cout << " OrMud[" << i << "]=" << fOrMud[i];
562 dump("Xpos11",fXpos11,16);
563 dump("Ypos11",fYpos11,31);
564 dump("Ypos21",fYpos21,63);
565 for ( Int_t i = 0; i < 4; ++i )
568 sprintf(s,"Xcode[%d]",i);
569 dump(s,fXcode[i],32);
570 sprintf(s,"Ycode[%d]",i);
571 dump(s,fYcode[i],32);
573 // Int_t fIdCircuit; // circuit Id number
574 // Int_t fX2m; // internal info needed by TriggerDecision
575 // Int_t fX2ud; // internal info needed by TriggerDecision
576 // Int_t fOrMud[2]; // internal info needed by TriggerDecision
577 // Int_t fXcode[4][32]; // code of X strips
578 // Int_t fYcode[4][32]; // code of Y strips
579 // Float_t fXpos11[16]; // X position of Y strips in MC11
580 // Float_t fYpos11[31]; // Y position of X strips in MC11
581 // Float_t fYpos21[63]; // Y position of X strips in MC21
585 Int_t AliMUONTriggerCircuit::DetElemId(Int_t ichamber, Int_t idModule)
587 // adpated to official numbering (09/20/05)
588 // returns the detection element Id for given chamber and module
589 // ichamber (from 11 to 14), idModule (from -97 to 97)
592 Int_t linenumber=Int_t(idModule/10);
593 switch (linenumber) // (from 1 to 9, from top to bottom)
651 return (ichamber*100)+itmp;