Change of argument list of methods GetPadCxy, GetPadIxy, SetHit and FirstPad
[u/mrichter/AliRoot.git] / MUON / AliMUONTriggerCircuit.cxx
CommitLineData
a9e2aefa 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/*
aac5be80 16$Log$
a9e2aefa 17*/
aac5be80 18
a9e2aefa 19#include "AliRun.h"
20#include "AliMUON.h"
21#include "AliMUONPoints.h"
22#include "AliMUONTriggerCircuit.h"
aac5be80 23#include "AliMUONTriggerConstants.h"
a9e2aefa 24#include "AliMUONSegmentation.h"
25#include "AliMUONResponse.h"
26#include "AliMUONChamber.h"
27#include "TMath.h"
28#include "iostream.h"
29
30ClassImp(AliMUONTriggerCircuit)
31
a9e2aefa 32//----------------------------------------------------------------------
33AliMUONTriggerCircuit::AliMUONTriggerCircuit()
34{
35// Constructor
36 fidCircuit=0;
37 fx2m=0;
38 fx2ud=0;
39 fOrMud[0]=fOrMud[1]=0;
aac5be80 40 Int_t i;
41 for (i=0; i<4; i++) {
a9e2aefa 42 for (Int_t j=0; j<32; j++) {
43 fXcode[i][j]=0;
44 fYcode[i][j]=0;
45 }
46 }
aac5be80 47 for (i=0; i<16; i++) { fXpos11[i]=0.; }
48 for (i=0; i<31; i++) { fYpos11[i]=0.; }
49 for (i=0; i<63; i++) { fYpos21[i]=0.; }
a9e2aefa 50}
51
52//----------------------------------------------------------------------
53AliMUONTriggerCircuit::AliMUONTriggerCircuit(const AliMUONTriggerCircuit& MUONTriggerCircuit)
54{
55// Dummy copy constructor
56}
57
58//----------------------------------------------------------------------
59AliMUONTriggerCircuit & AliMUONTriggerCircuit::operator=(const AliMUONTriggerCircuit& MUONTriggerCircuit)
60{
61// Dummy assignment operator
62 return *this;
63}
64
65//----------------------------------------------------------------------
66void AliMUONTriggerCircuit::Init(Int_t iCircuit) {
67// initialize circuit characteristics
aac5be80 68 fidCircuit=AliMUONTriggerConstants::CircuitId(iCircuit);
a9e2aefa 69 LoadX2();
70 LoadXCode();
71 LoadYCode();
72 LoadXPos();
73 LoadYPos();
74}
75
76//----------------------------------------------------------------------
77Int_t AliMUONTriggerCircuit::CircuitNumber(Int_t idCircuit){
78// returns circuit number iCircuit (0-234) corresponding to circuit idCircuit
79 Int_t iCircuit=0;
80 for (Int_t i=0; i<234; i++) {
aac5be80 81 if (AliMUONTriggerConstants::CircuitId(i)==idCircuit) {
a9e2aefa 82 iCircuit=i;
83 break;
84 }
85 }
86 return iCircuit;
87}
88//----------------------------------------------------------------------
89Int_t AliMUONTriggerCircuit::ModuleNumber(Int_t idModule){
90// returns module number imod (from 0 to 63) corresponding to module idmodule
91 Int_t absidModule=TMath::Abs(idModule);
92 Int_t iModule=0;
93 for (Int_t i=0; i<63; i++) {
aac5be80 94 if (AliMUONTriggerConstants::ModuleId(i)==absidModule) {
a9e2aefa 95 iModule=i;
96 break;
97 }
98 }
99 return iModule;
100}
101
102//----------------------------------------------------------------------
103Int_t AliMUONTriggerCircuit::Module(Int_t idCircuit) {
104// returns ModuleId where Circuit idCircuit is sitting
105 return Int_t(idCircuit/10);
106}
107//----------------------------------------------------------------------
108Int_t AliMUONTriggerCircuit::Position(Int_t idCircuit) {
109// returns position of idCircuit in correcponding Module
110 return TMath::Abs(idCircuit)-TMath::Abs(Module(idCircuit))*10;
111}
112
113//----------------------------------------------------------------------
114void AliMUONTriggerCircuit::LoadX2() {
115// initialize fx2m, fx2ud and fOrMud
116
117 Int_t idModule=Module(fidCircuit); // corresponding module Id.
aac5be80 118// and its number of X strips
119 Int_t nStrX=AliMUONTriggerConstants::NstripX(ModuleNumber(idModule));
120// and its number of Y strips
121 Int_t nStrY=AliMUONTriggerConstants::NstripY(ModuleNumber(idModule));
a9e2aefa 122 Int_t iPosCircuit=Position(fidCircuit); // position of circuit in module
123
124// first step : look at lower part
125 if (iPosCircuit==1) { // need to scan lower module
126 if(idModule<91&&TMath::Abs(idModule)!=41&&idModule>-91) {
127 fOrMud[0]=1;
128 Int_t idModuleD=(TMath::Abs(idModule)+10)*(TMath::Abs(idModule)/idModule);
aac5be80 129 Int_t nStrD=AliMUONTriggerConstants::NstripY(ModuleNumber(idModuleD));
a9e2aefa 130
131 if (nStrY!=nStrD
132 &&TMath::Abs(idModule)!=42&&TMath::Abs(idModule)!=52) {
133 if (nStrY==8) fx2m=1;
134 if (nStrD==8) fx2ud=1;
135 }
136 }
137
138 } else { // lower strips within same module
139 fOrMud[0]=0;
140 }
141
142// second step : look at upper part
143 if ((iPosCircuit==1&&nStrX==16)||(iPosCircuit==2&&nStrX==32)||
144 (iPosCircuit==3&&nStrX==48)||(iPosCircuit==4&&nStrX==64)) {
145 if ((idModule>17||idModule<-17)&&TMath::Abs(idModule)!=61) {
146 fOrMud[1]=1;
147 Int_t idModuleU=(TMath::Abs(idModule)-10)*(TMath::Abs(idModule)/idModule);
aac5be80 148 Int_t nStrU=AliMUONTriggerConstants::NstripY(ModuleNumber(idModuleU));
a9e2aefa 149
150 if (nStrY!=nStrU
151 &&TMath::Abs(idModule)!=62&&TMath::Abs(idModule)!=52) {
152 if (nStrY==8) fx2m=1;
153 if (nStrU==8) fx2ud=1;
154 }
155 }
156
157 } else { // upper strips within same module
158 fOrMud[1]=0;
159 }
160}
161
162//----------------------------------------------------------------------
163void AliMUONTriggerCircuit::LoadXCode(){
164// assign a Id. number to each X strip of current circuit
165// Id.=(corresponding module Id.)*100+(Id. strip of module)
166
167// first part : fill XMC11 XMC12 and strips 8 to 24 (middle) XMC21 XMC22
168 Int_t iStripCircMT1=0, iStripCircMT2=8;
169 Int_t idModule=Module(fidCircuit); // corresponding module Id.
aac5be80 170// and its number of strips
171 Int_t nStrX=AliMUONTriggerConstants::NstripX(ModuleNumber(idModule));
a9e2aefa 172 Int_t iPosCircuit=Position(fidCircuit); // position of circuit in module
173 Int_t sign=TMath::Abs(idModule)/idModule; // left or right
aac5be80 174 Int_t istrip;
a9e2aefa 175
aac5be80 176 for (istrip=(iPosCircuit-1)*16;
a9e2aefa 177 istrip<(iPosCircuit-1)*16+16; istrip++) {
178
179 fXcode[0][iStripCircMT1]=sign*(TMath::Abs(idModule)*100+istrip);
180 fXcode[1][iStripCircMT1]=sign*(TMath::Abs(idModule)*100+istrip);
181 fXcode[2][iStripCircMT2]=sign*(TMath::Abs(idModule)*100+istrip);
182 fXcode[3][iStripCircMT2]=sign*(TMath::Abs(idModule)*100+istrip);
183 iStripCircMT1++;
184 iStripCircMT2++;
185 }
186
187// second part
188// XMC21 XMC22 strips 0 to 7 and 24 to 31
189 Int_t idModuleD, idModuleU;
190 Int_t nStrD, nStrU;
191
192 idModule=Module(fidCircuit); // corresponding module Id.
aac5be80 193// number of X strips
194 nStrX=AliMUONTriggerConstants::NstripX(ModuleNumber(idModule));
a9e2aefa 195 sign=TMath::Abs(idModule)/idModule;
196
197// fill lower part (0 to 7)
198 if (iPosCircuit==1) { // need to scan lower module
199 if(idModule<91&&TMath::Abs(idModule)!=41&&idModule>-91) { // non-existing
200 idModuleD=sign*(TMath::Abs(idModule)+10); // lower module Id
aac5be80 201// and its number of strips
202 nStrD=AliMUONTriggerConstants::NstripX(ModuleNumber(idModuleD));
a9e2aefa 203
204 iStripCircMT2=0;
aac5be80 205 for (istrip=nStrD-8; istrip<nStrD; istrip++) {
a9e2aefa 206 fXcode[2][iStripCircMT2]=sign*(TMath::Abs(idModuleD)*100+istrip);
207 fXcode[3][iStripCircMT2]=sign*(TMath::Abs(idModuleD)*100+istrip);
208 iStripCircMT2++;
209 }
210 }
211
212 } else { // lower strips within same module
213
214 iStripCircMT2=0;
aac5be80 215 for (istrip=(iPosCircuit-1)*16-8;
a9e2aefa 216 istrip<(iPosCircuit-1)*16; istrip++) {
217 fXcode[2][iStripCircMT2]=sign*(TMath::Abs(idModule)*100+istrip);
218 fXcode[3][iStripCircMT2]=sign*(TMath::Abs(idModule)*100+istrip);
219 iStripCircMT2++;
220 }
221 }
222
223// fill upper part (24 to 31)
224 if ((iPosCircuit==1&&nStrX==16)||(iPosCircuit==2&&nStrX==32)||
225 (iPosCircuit==3&&nStrX==48)||(iPosCircuit==4&&nStrX==64)) {
226 if ((idModule>17||idModule<-17)&&TMath::Abs(idModule)!=61) {
227 idModuleU=sign*(TMath::Abs(idModule)-10); // upper module Id
aac5be80 228// and its number of strips
229 nStrU=AliMUONTriggerConstants::NstripX(ModuleNumber(idModuleU));
a9e2aefa 230
231 iStripCircMT2=24;
aac5be80 232 for (istrip=0; istrip<8; istrip++) {
a9e2aefa 233 fXcode[2][iStripCircMT2]=sign*(TMath::Abs(idModuleU)*100+istrip);
234 fXcode[3][iStripCircMT2]=sign*(TMath::Abs(idModuleU)*100+istrip);
235 iStripCircMT2++;
236 }
237 }
238
239 } else if ((iPosCircuit==1&&nStrX>16)||(iPosCircuit==2&&nStrX>32)||
240 (iPosCircuit==3&&nStrX>48)) { // upper strips within same mod.
241
242 iStripCircMT2=24;
aac5be80 243 for (istrip=(iPosCircuit-1)*16+16;
a9e2aefa 244 istrip<(iPosCircuit-1)*16+24; istrip++) {
245 fXcode[2][iStripCircMT2]=sign*(TMath::Abs(idModule)*100+istrip);
246 fXcode[3][iStripCircMT2]=sign*(TMath::Abs(idModule)*100+istrip);
247 iStripCircMT2++;
248 }
249 }
250}
251
252//----------------------------------------------------------------------
253void AliMUONTriggerCircuit::LoadYCode(){
254// assign a Id. number to each Y strip of current circuit
255// Id.=(corresponding module Id.)*100+(Id. strip of module)
256// note : for Y plane fill only "central part" of circuit
257// (upper and lower parts are filled in PreHandlingY of AliMUONTriggerDecision)
258
259 Int_t idModule=Module(fidCircuit); // corresponding module Id.
aac5be80 260// and its number of Y strips
261 Int_t nStrY=AliMUONTriggerConstants::NstripY(ModuleNumber(idModule));
a9e2aefa 262 Int_t sign=TMath::Abs(idModule)/idModule; // left or right
263
264 for (Int_t istrip=0; istrip<nStrY; istrip++) {
265 fYcode[0][istrip]=sign*(TMath::Abs(idModule)*100+istrip);
266 fYcode[1][istrip]=sign*(TMath::Abs(idModule)*100+istrip);
267 fYcode[2][istrip]=sign*(TMath::Abs(idModule)*100+istrip);
268 fYcode[3][istrip]=sign*(TMath::Abs(idModule)*100+istrip);
269 }
270}
271
272//----------------------------------------------------------------------
273void AliMUONTriggerCircuit::LoadYPos(){
274// fill fYpos11 and fYpos21 -> y position of X declusterized strips
275
276 Int_t chamber, cathode;
277 Int_t code, idModule, idStrip, idSector;
278 Float_t x, y, width;
aac5be80 279 Int_t istrip;
a9e2aefa 280
281 AliMUON *pMUON = (AliMUON*)gAlice->GetModule("MUON");
282 AliMUONChamber* iChamber;
283 AliMUONSegmentation* segmentation;
284
285// first plane (11)
286 chamber=11;
287 cathode=1;
288 iChamber = &(pMUON->Chamber(chamber-1));
289 segmentation=iChamber->SegmentationModel(cathode);
290
aac5be80 291 for (istrip=0; istrip<16; istrip++) {
a9e2aefa 292 code=fXcode[0][istrip]; // decode current strip
293 idModule=Int_t(code/100); // corresponding module Id.
294 idStrip=TMath::Abs(code-idModule*100); // corresp. strip number in module
295 idSector=segmentation->Sector(idModule,idStrip); // corresponding sector
296 width=segmentation->Dpy(idSector); // corresponding strip width
297 segmentation->GetPadCxy(idModule,idStrip,x,y); // get strip real position
298
299 fYpos11[2*istrip]=y;
300 if (istrip!=15) fYpos11[2*istrip+1]=y+width/2.;
301 }
302
303// second plane (21)
304 chamber=13;
305 cathode=1;
306 iChamber = &(pMUON->Chamber(chamber-1));
307 segmentation=iChamber->SegmentationModel(cathode);
308
aac5be80 309 for (istrip=0; istrip<32; istrip++) {
a9e2aefa 310 code=fXcode[2][istrip]; // decode current strip
311 idModule=Int_t(code/100); // corresponding module Id.
312 idStrip=TMath::Abs(code-idModule*100); // corresp. strip number in module
313 idSector=segmentation->Sector(idModule,idStrip); // corresponding sector
314 width=segmentation->Dpy(idSector); // corresponding strip width
315 segmentation->GetPadCxy(idModule,idStrip,x,y); // get strip real position
316
317// using idModule!=0 prevents to fill garbage in case of circuits
318// in the first and last rows
319 if (idModule!=0) {
320 fYpos21[2*istrip]=y;
321 if (istrip!=31) fYpos21[2*istrip+1]=y+width/2.;
322 }
323 }
324}
325
326//----------------------------------------------------------------------
327void AliMUONTriggerCircuit::LoadXPos(){
328// fill fXpos11 -> x position of Y strips for the first plane only
329// fXpos11 contains the x position of Y strip for the current circuit
330// taking into account whether or nor not part(s) of the circuit
331// (middle, up or down) has(have) 16 strips
332
333 Float_t x, y;
aac5be80 334 Int_t istrip;
a9e2aefa 335
336 Int_t chamber=11;
337 Int_t cathode=2;
338 AliMUON *pMUON = (AliMUON*)gAlice->GetModule("MUON");
339 AliMUONChamber* iChamber;
340 AliMUONSegmentation* segmentation;
341 iChamber = &(pMUON->Chamber(chamber-1));
342 segmentation=iChamber->SegmentationModel(cathode);
343
344 Int_t idModule=Module(fidCircuit); // corresponding module Id.
aac5be80 345// number of Y strips
346 Int_t nStrY=AliMUONTriggerConstants::NstripY(ModuleNumber(idModule));
a9e2aefa 347 Int_t idSector=segmentation->Sector(idModule,0); // corresp. sector
348 Float_t width=segmentation->Dpx(idSector); // corresponding strip width
349
350// first case : up middle and down parts have all 8 or 16 strip
351 if ((nStrY==16)||(nStrY==8&&fx2m==0&&fx2ud==0)) {
aac5be80 352 for (istrip=0; istrip<nStrY; istrip++) {
a9e2aefa 353 segmentation->GetPadCxy(idModule,istrip,x,y);
354 fXpos11[istrip]=x;
355 }
356// second case : mixing 8 and 16 strips within same circuit
357 } else {
aac5be80 358 for (istrip=0; istrip<nStrY; istrip++) {
a9e2aefa 359 if (nStrY!=8) { cout << " bug in LoadXpos " << "\n";}
360 segmentation->GetPadCxy(idModule,istrip,x,y);
361 fXpos11[2*istrip]=x-width/4.;
362 fXpos11[2*istrip+1]=fXpos11[2*istrip]+width/2.;
363 }
364 }
365}
366
367//----------------------------------------------------------------------
368Float_t AliMUONTriggerCircuit::PtCal(Int_t istripX, Int_t idev, Int_t istripY){
369// returns calculated pt for circuit/istripX/idev/istripY according
370// to the formula of the TRD. Note : idev (input) is in [0+30]
371
372 // Int_t jdev = idev - 15; // jdev in [-15+15]
373 Int_t istripX2=istripX+idev+1; // find istripX2 using istripX and idev
374
375 Float_t yPosX1=fYpos11[istripX];
376 Float_t yPosX2=fYpos21[istripX2];
377 Float_t xPosY1=fXpos11[istripY];
378
379 Float_t zf=975., z1=1603.5, z2=1703.5;
380 Float_t thetaDev=(1./zf)*(yPosX1*z2-yPosX2*z1)/(z2-z1);
381 Float_t xf=xPosY1*zf/z1;
382 Float_t yf=yPosX2-((yPosX2-yPosX1)*(z2-zf))/(z2-z1);
383 return (3.*0.3/TMath::Abs(thetaDev)) * TMath::Sqrt(xf*xf+yf*yf)/zf;
384}
385
386//----------------------------------------------------------------------
387//--- methods which return member data related info
388//----------------------------------------------------------------------
389Int_t AliMUONTriggerCircuit::GetIdCircuit(){
390// returns circuit Id
391 return fidCircuit;
392}
393//----------------------------------------------------------------------
394Int_t AliMUONTriggerCircuit::GetIdModule(){
395// returns module Id
396 return Module(fidCircuit);
397}
398//----------------------------------------------------------------------
399Int_t AliMUONTriggerCircuit::GetNstripX() {
400// returns the number of X strips in the module where the circuit is sitting
aac5be80 401 return AliMUONTriggerConstants::NstripX(ModuleNumber(Module(fidCircuit)));
a9e2aefa 402}
403//----------------------------------------------------------------------
404Int_t AliMUONTriggerCircuit::GetNstripY() {
405// returns the number of Y strips in the module where the circuit is sitting
aac5be80 406 return AliMUONTriggerConstants::NstripY(ModuleNumber(Module(fidCircuit)));
a9e2aefa 407}
408//----------------------------------------------------------------------
409Int_t AliMUONTriggerCircuit::GetPosCircuit() {
410// returns the position of the circuit in its module
411 return Position(fidCircuit);
412}
413//----------------------------------------------------------------------
414Int_t AliMUONTriggerCircuit::GetIdCircuitD(){
415// returns the Id of the circuit down
416 Int_t idModule=Module(fidCircuit);
417 Int_t idModuleD=(TMath::Abs(idModule)+10)*(TMath::Abs(idModule)/idModule);
418 return (TMath::Abs(idModuleD)*10+1)*(TMath::Abs(idModule)/idModule);
419}
420//----------------------------------------------------------------------
421Int_t AliMUONTriggerCircuit::GetICircuitD(){
422// returns the number of the circuit down
423 Int_t idModule=Module(fidCircuit);
424 Int_t idModuleD=(TMath::Abs(idModule)+10)*(TMath::Abs(idModule)/idModule);
425 Int_t idCircuitD=
426 (TMath::Abs(idModuleD)*10+1)*(TMath::Abs(idModule)/idModule);
427 return CircuitNumber(idCircuitD);
428}
429//----------------------------------------------------------------------
430Int_t AliMUONTriggerCircuit::GetIdCircuitU(){
431// returns the Id of the circuit up
432 Int_t idModule=Module(fidCircuit);
433 Int_t idModuleU=(TMath::Abs(idModule)-10)*(TMath::Abs(idModule)/idModule);
434 return (TMath::Abs(idModuleU)*10+1)*(TMath::Abs(idModule)/idModule);
435}
436//----------------------------------------------------------------------
437Int_t AliMUONTriggerCircuit::GetICircuitU(){
438// returns the number of the circuit up
439 Int_t idModule=Module(fidCircuit);
440 Int_t idModuleU=(TMath::Abs(idModule)-10)*(TMath::Abs(idModule)/idModule);
441 Int_t idCircuitU=
442 (TMath::Abs(idModuleU)*10+1)*(TMath::Abs(idModule)/idModule);
443 return CircuitNumber(idCircuitU);
444}
445//----------------------------------------------------------------------
446Int_t AliMUONTriggerCircuit::GetX2m(){
447// returns fx2m
448 return fx2m;
449}
450//----------------------------------------------------------------------
451Int_t AliMUONTriggerCircuit::GetX2ud(){
452// returns fx2ud
453 return fx2ud;
454}
455//----------------------------------------------------------------------
456void AliMUONTriggerCircuit::GetOrMud(Int_t orMud[2]){
457// returns fOrMud
458 orMud[0]=fOrMud[0];
459 orMud[1]=fOrMud[1];
460}
461//----------------------------------------------------------------------
462Int_t AliMUONTriggerCircuit::GetXcode(Int_t chamber, Int_t istrip){
463// returns X code of circuit/chamber/istrip (warning : chamber in [0,3])
464 return fXcode[chamber][istrip];
465}
466//----------------------------------------------------------------------
467Int_t AliMUONTriggerCircuit::GetYcode(Int_t chamber, Int_t istrip){
468// returns Y code of circuit/chamber/istrip (warning : chamber in [0,3])
469 return fYcode[chamber][istrip];
470}
471//----------------------------------------------------------------------
472Float_t AliMUONTriggerCircuit::GetY11Pos(Int_t istrip){
473// returns Y position of X strip istrip in MC11
474 return fYpos11[istrip];
475}
476//----------------------------------------------------------------------
477Float_t AliMUONTriggerCircuit::GetY21Pos(Int_t istrip){
478// returns Y position of X strip istrip in MC21
479 return fYpos21[istrip];
480}
481//----------------------------------------------------------------------
482Float_t AliMUONTriggerCircuit::GetX11Pos(Int_t istrip){
483// returns X position of Y strip istrip in MC11
484 return fXpos11[istrip];
485}
486//----------------------------------------------------------------------
487//--- end of methods which return member data related info
488//----------------------------------------------------------------------
489
490
491
492
493
494
495