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 **************************************************************************/
16 ////////////////////////////////////////////////////////////
17 // Factory for muon chambers, segmentations and response //
18 ////////////////////////////////////////////////////////////
22 #include "AliMUONFactoryV2.h"
24 #include "AliSegmentation.h"
28 #include "AliMUONConstants.h"
29 #include "AliMUONTriggerConstants.h"
30 #include "AliMUONChamber.h"
31 #include "AliMUONResponseV0.h"
32 #include "AliMUONGeometryModule.h"
33 #include "AliMUONGeometryStore.h"
34 #include "AliMUONGeometrySegmentation.h"
35 #include "AliMUONSt12QuadrantSegmentation.h"
36 #include "AliMUONSt345SlatSegmentation.h"
37 #include "AliMUONTriggerSegmentation.h"
38 #include "AliMUONResponseTrigger.h"
40 ClassImp(AliMUONFactoryV2)
42 //__________________________________________________________________________
43 AliMUONFactoryV2::AliMUONFactoryV2(const char* name)
44 : AliMUONFactory(name),
49 // FactoryV2 inherite from Factory for switching in AliMUONv1::Init()
50 // to be changed when old segmentation will be removed.
52 fDESegmentations = new TObjArray();
53 fDESegmentations->SetOwner(kTRUE);
56 //__________________________________________________________________________
57 AliMUONFactoryV2::AliMUONFactoryV2()
63 // Default constructor
66 //__________________________________________________________________________
67 AliMUONFactoryV2::AliMUONFactoryV2(const AliMUONFactoryV2& rhs)
70 // Protected copy constructor
72 AliFatal("Not implemented.");
75 //__________________________________________________________________________
77 AliMUONFactoryV2::~AliMUONFactoryV2()
81 delete fDESegmentations;
84 //__________________________________________________________________________
85 AliMUONFactoryV2& AliMUONFactoryV2::operator=(const AliMUONFactoryV2& rhs)
87 // Protected assignement operator
89 if (this == &rhs) return *this;
91 AliFatal("Not implemented.");
96 //__________________________________________________________________________
97 void AliMUONFactoryV2::BuildCommon()
100 // Construct the default response.
103 // Default response: 5 mm of gas
104 fResponse0 = new AliMUONResponseV0;
105 fResponse0->SetSqrtKx3AndDeriveKx2Kx4(0.7131); // sqrt(0.5085)
106 fResponse0->SetSqrtKy3AndDeriveKy2Ky4(0.7642); // sqrt(0.5840)
107 fResponse0->SetPitch(AliMUONConstants::Pitch()); // anode-cathode distance
108 fResponse0->SetSigmaIntegration(10.);
109 fResponse0->SetChargeSlope(10);
110 fResponse0->SetChargeSpread(0.18, 0.18);
111 fResponse0->SetMaxAdc(4096);
112 fResponse0->SetSaturation(3000);
113 fResponse0->SetZeroSuppression(6);
116 //__________________________________________________________________________
117 void AliMUONFactoryV2::BuildStation1()
119 //--------------------------------------------------------
120 // Configuration for Chamber TC1/2 (Station 1) ----------
121 //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
124 // Response for 4 mm of gas (station 1)
125 // automatic consistency with width of sensitive medium in CreateGeometry ????
126 AliMUONResponseV0* responseSt1 = new AliMUONResponseV0;
127 // Mathieson parameters from L.Kharmandarian's thesis, page 190
128 responseSt1->SetSqrtKx3AndDeriveKx2Kx4(0.7000); // sqrt(0.4900)
129 responseSt1->SetSqrtKy3AndDeriveKy2Ky4(0.7550); // sqrt(0.5700)
130 responseSt1->SetPitch(AliMUONConstants::PitchSt1()); // anode-cathode distance
131 responseSt1->SetSigmaIntegration(10.);
132 // ChargeSlope larger to compensate for the smaller anode-cathode distance
133 // and keep the same most probable ADC channel for mip's
134 responseSt1->SetChargeSlope(62.5);
135 // assumed proportionality to anode-cathode distance for ChargeSpread
136 responseSt1->SetChargeSpread(0.144, 0.144);
137 responseSt1->SetMaxAdc(4096);
138 responseSt1->SetSaturation(3000);
139 responseSt1->SetZeroSuppression(6);
141 // Quadrant segmentations:
142 AliMUONSt12QuadrantSegmentation* bendSt1
143 = new AliMUONSt12QuadrantSegmentation(kStation1, kBendingPlane);
144 AliMUONSt12QuadrantSegmentation* nonbendSt1
145 = new AliMUONSt12QuadrantSegmentation(kStation1, kNonBendingPlane);
147 // Add in the array (for safe deleting)
148 fDESegmentations->Add(bendSt1);
149 fDESegmentations->Add(nonbendSt1);
151 AliMUONGeometrySegmentation* segmentation[2];
153 for (Int_t chamber = 0; chamber < 2; chamber++) {
155 segmentation[0] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
156 segmentation[1] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
158 // id detection elt for chamber 1
159 Int_t id0 = (chamber+1)*100;
161 //--------------------------------------------------------
162 // Configuration for Chamber TC1/2 (Station 1) ----------
165 fMUON->SetNsec(chamber,2);
168 segmentation[0]->Add(id0, bendSt1);
169 segmentation[0]->Add(id0 + 1, nonbendSt1);
170 segmentation[0]->Add(id0 + 50, bendSt1);
171 segmentation[0]->Add(id0 + 51, nonbendSt1);
172 fMUON->SetSegmentationModel(chamber, 1, segmentation[0]);
175 segmentation[1]->Add(id0, nonbendSt1);
176 segmentation[1]->Add(id0 + 1, bendSt1);
177 segmentation[1]->Add(id0 + 50, nonbendSt1);
178 segmentation[1]->Add(id0 + 51, bendSt1);
179 fMUON->SetSegmentationModel(chamber, 2, segmentation[1]);
181 fMUON->SetResponseModel(chamber, responseSt1); // special response
182 fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread
187 //__________________________________________________________________________
188 void AliMUONFactoryV2::BuildStation2()
191 //--------------------------------------------------------
192 // Configuration for Chamber TC3/4 (Station 2) -----------
193 ///^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
196 // Quadrant segmentations:
197 AliMUONSt12QuadrantSegmentation* bendSt2
198 = new AliMUONSt12QuadrantSegmentation(kStation2, kBendingPlane);
199 AliMUONSt12QuadrantSegmentation* nonbendSt2
200 = new AliMUONSt12QuadrantSegmentation(kStation2, kNonBendingPlane);
202 // Add in the array (for safe deleting)
203 fDESegmentations->Add(bendSt2);
204 fDESegmentations->Add(nonbendSt2);
206 AliMUONGeometrySegmentation* segmentation[2];
208 for (Int_t chamber = 2; chamber < 4; chamber++) {
210 segmentation[0] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
211 segmentation[1] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
213 // id detection elt for chamber 1
214 Int_t id0 = (chamber+1)*100;
216 //--------------------------------------------------------
217 // Configuration for Chamber TC3/4 (Station 2) ----------
220 fMUON->SetNsec(chamber,2);
223 segmentation[0]->Add(id0, bendSt2);
224 segmentation[0]->Add(id0 + 1, nonbendSt2);
225 segmentation[0]->Add(id0 + 50, bendSt2);
226 segmentation[0]->Add(id0 + 51, nonbendSt2);
227 fMUON->SetSegmentationModel(chamber, 1, segmentation[0]);
230 segmentation[1]->Add(id0, nonbendSt2);
231 segmentation[1]->Add(id0 + 1, bendSt2);
232 segmentation[1]->Add(id0 + 50, nonbendSt2);
233 segmentation[1]->Add(id0 + 51, bendSt2);
234 fMUON->SetSegmentationModel(chamber, 2, segmentation[1]);
236 fMUON->SetResponseModel(chamber, fResponse0); // normal response
237 fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread
243 //__________________________________________________________________________
244 void AliMUONFactoryV2::BuildStation3()
246 //--------------------------------------------------------
247 // Configuration for Chamber TC5/6 (Station 3) ----------
248 //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
250 AliMUONGeometrySegmentation* segmentation[2];
252 //Slats Segmentations
253 AliMUONSt345SlatSegmentation *slatsegB[4]; // Types of segmentation for St3
254 AliMUONSt345SlatSegmentation *slatsegNB[4];
257 Int_t ndiv[4] ={ 4, 4, 2, 1}; // densities zones
258 for(Int_t i=0; i<4; i++) {
259 slatsegB[i] = new AliMUONSt345SlatSegmentation(1);
260 fDESegmentations->Add(slatsegB[i]);
261 slatsegB[i]->SetPadSize(10.,0.5);
262 slatsegB[i]->SetPadDivision(ndiv);
263 slatsegB[i]->SetId(1); // Id elt ????
264 slatsegB[i]->SetDAnod(AliMUONConstants::Pitch());
265 slatsegNB[i] = new AliMUONSt345SlatSegmentation(0);
266 fDESegmentations->Add(slatsegNB[i]);
267 slatsegNB[i]->SetPadSize(0.713,10.); // Nbending
268 slatsegNB[i]->SetPadDivision(ndiv);
269 slatsegNB[i]->SetId(1);
270 slatsegNB[i]->SetDAnod(AliMUONConstants::Pitch());
273 // Type 112200 for 500, 501, 508, 509, 510, 517
274 // in Ch5 (similar for Ch6) for the futur official numbering
275 // Type 112200 for 503, 504, 505, 555, 554, 553
276 // in Ch5 (similar for Ch6) actual numbering in the code to be changed in jan05
277 Int_t n0[4] = { 0, 2, 2, 0 };
278 slatsegB[0]->SetPcbBoards(n0);
279 slatsegB[0]->Init(0);
280 slatsegNB[0]->SetPcbBoards(n0);
281 slatsegNB[0]->Init(0);
283 // Type 122200 for 502, 507, 511, 516 (similar in Ch6)
284 // for future official numbering of ALICE
285 // Type 122200 for 502, 506, 556, 552 (similiarin Ch6)
286 // for actual numbering in muon code to be changed in jan05
287 Int_t n1[4] = { 0, 1, 3, 0 };
288 slatsegB[1]->SetPcbBoards(n1);
289 slatsegB[1]->Init(0);
290 slatsegNB[1]->SetPcbBoards(n1);
291 slatsegNB[1]->Init(0);
293 // Type 222000 for 503, 506, 512, 515 (similar in Ch6)
294 // for future official numbering of ALICE
295 // Type 222000 for 501, 507, 557, 551 (similiarin Ch6)
296 // for actual numbering in muon code to be changed in jan05
297 Int_t n2[4] = { 0, 0, 3, 0 };
298 slatsegB[2]->SetPcbBoards(n2);
299 slatsegB[2]->Init(0);
300 slatsegNB[2]->SetPcbBoards(n2);
301 slatsegNB[2]->Init(0);
303 // Type 220000 for 504, 505, 513, 514 (similar in Ch6)
304 // for future official numbering of ALICE
305 // Type 220000 for 500, 508, 558, 550 (similiarin Ch6)
306 // for actual numbering in muon code to be changed in jan05
307 Int_t n3[4] = { 0, 0, 2, 0 };
308 slatsegB[3]->SetPcbBoards(n3);
309 slatsegB[3]->Init(0);
310 slatsegNB[3]->SetPcbBoards(n3);
311 slatsegNB[3]->Init(0);
313 for (Int_t chamber = 4; chamber < 6; chamber++) {
315 segmentation[0] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
316 segmentation[1] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
318 // id detection elt for chamber 1
319 Int_t id0 = (chamber+1)*100;
321 //--------------------------------------------------------
322 // Configuration for Chamber TC3/4 (Station 2) ----------
325 fMUON->SetNsec(chamber,2);
329 segmentation[0]->Add(id0, slatsegB[3]);
330 segmentation[0]->Add(id0+ 8, slatsegB[3]);
331 segmentation[0]->Add(id0+50, slatsegB[3]);
332 segmentation[0]->Add(id0+58, slatsegB[3]);
334 segmentation[0]->Add(id0+ 1, slatsegB[2]);
335 segmentation[0]->Add(id0+ 7, slatsegB[2]);
336 segmentation[0]->Add(id0+51, slatsegB[2]);
337 segmentation[0]->Add(id0+57, slatsegB[2]);
339 segmentation[0]->Add(id0+ 2, slatsegB[1]);
340 segmentation[0]->Add(id0+ 6, slatsegB[1]);
341 segmentation[0]->Add(id0+52, slatsegB[1]);
342 segmentation[0]->Add(id0+56, slatsegB[1]);
344 segmentation[0]->Add(id0+ 3, slatsegB[0]);
345 segmentation[0]->Add(id0+ 4, slatsegB[0]);
346 segmentation[0]->Add(id0+ 5, slatsegB[0]);
347 segmentation[0]->Add(id0+53, slatsegB[0]);
348 segmentation[0]->Add(id0+54, slatsegB[0]);
349 segmentation[0]->Add(id0+55, slatsegB[0]);
350 fMUON->SetSegmentationModel(chamber, 1, segmentation[0]);
354 segmentation[1]->Add(id0, slatsegNB[3]);
355 segmentation[1]->Add(id0+ 8, slatsegNB[3]);
356 segmentation[1]->Add(id0+50, slatsegNB[3]);
357 segmentation[1]->Add(id0+58, slatsegNB[3]);
359 segmentation[1]->Add(id0+ 1, slatsegNB[2]);
360 segmentation[1]->Add(id0+ 7, slatsegNB[2]);
361 segmentation[1]->Add(id0+51, slatsegNB[2]);
362 segmentation[1]->Add(id0+57, slatsegNB[2]);
364 segmentation[1]->Add(id0+ 2, slatsegNB[1]);
365 segmentation[1]->Add(id0+ 6, slatsegNB[1]);
366 segmentation[1]->Add(id0+52, slatsegNB[1]);
367 segmentation[1]->Add(id0+56, slatsegNB[1]);
369 segmentation[1]->Add(id0+ 3, slatsegNB[0]);
370 segmentation[1]->Add(id0+ 4, slatsegNB[0]);
371 segmentation[1]->Add(id0+ 5, slatsegNB[0]);
372 segmentation[1]->Add(id0+53, slatsegNB[0]);
373 segmentation[1]->Add(id0+54, slatsegNB[0]);
374 segmentation[1]->Add(id0+55, slatsegNB[0]);
375 fMUON->SetSegmentationModel(chamber, 2, segmentation[1]);
377 fMUON->SetResponseModel(chamber, fResponse0);
378 fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread
383 //__________________________________________________________________________
384 void AliMUONFactoryV2::BuildStation4()
386 //--------------------------------------------------------
387 // Configuration for Chamber TC7/8 (Station 4) ----------
388 //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
391 AliMUONGeometrySegmentation* segmentation[2];
393 //Slats Segmentations
394 AliMUONSt345SlatSegmentation *slatsegB[7]; // Types of segmentation for St4
395 AliMUONSt345SlatSegmentation *slatsegNB[7];
398 Int_t ndiv[4] ={ 4, 4, 2, 1}; // densities zones
399 for(Int_t i = 0; i < 7; i++) {
400 slatsegB[i] = new AliMUONSt345SlatSegmentation(1);
401 fDESegmentations->Add(slatsegB[i]);
402 slatsegB[i]->SetPadSize(10.,0.5);
403 slatsegB[i]->SetPadDivision(ndiv);
404 slatsegB[i]->SetId(1);
405 slatsegB[i]->SetDAnod(AliMUONConstants::Pitch());
406 slatsegNB[i] = new AliMUONSt345SlatSegmentation(0);
407 fDESegmentations->Add(slatsegNB[i]);
408 slatsegNB[i]->SetPadSize(0.713,10.);
409 slatsegNB[i]->SetPadDivision(ndiv);
410 slatsegNB[i]->SetId(1);
411 slatsegNB[i]->SetDAnod(AliMUONConstants::Pitch());
414 Int_t n4[4] = { 0, 1, 2, 2 };
415 slatsegB[0]->SetPcbBoards(n4);
416 slatsegB[0]->Init(0); // 0 detection element id
417 slatsegNB[0]->SetPcbBoards(n4);
418 slatsegNB[0]->Init(0); // 0 detection element id
420 // Type 112233 for 701, 712, 714, 725 in Ch7 (similar for Ch8)
421 // for the futur official numbering
422 // Type 112233 for 705, 707, 755, 757 in Ch7 (similar for Ch8)
423 // actual numbering in the code to be changed in jan05
424 // Type 112233 for 901, 902, 911, 912, 914, 915, 924, 925 in Ch9
425 // (similar for Ch10) for the futur official numbering
426 // Type 112233 for 904, 905, 907, 908, 954, 955, 957, 958 in Ch9
427 // (similar for Ch10) actual numbering in the code to be changed in jan05
428 Int_t n5[4] = { 0, 2, 2, 2 };
429 slatsegB[1]->SetPcbBoards(n5);
430 slatsegB[1]->Init(0); // 0 detection element id
431 slatsegNB[1]->SetPcbBoards(n5);
432 slatsegNB[1]->Init(0); // 0 detection element id
434 // Type 112230 for 702, 711, 715, 724 in Ch7 (similar for Ch8)
435 // for the futur official numbering
436 // Type 112230 for 704, 708, 754, 758 in Ch7 (similar for Ch8)
437 // actual numbering in the code to be changed in jan05
438 Int_t n6[4] = { 0, 2, 2, 1 };
439 slatsegB[2]->SetPcbBoards(n6);
440 slatsegB[2]->Init(0); // 0 detection element id
441 slatsegNB[2]->SetPcbBoards(n6);
442 slatsegNB[2]->Init(0); // 0 detection element id
444 // Type 222330 for 703, 710, 716, 723 in Ch7 (similar for Ch8)
445 // for the futur official numbering
446 // Type 222330 for 703, 709, 753, 759 in Ch7 (similar for Ch8)
447 // actual numbering in the code to be changed in jan05
448 Int_t n7[4] = { 0, 0, 3, 2 };
449 slatsegB[3]->SetPcbBoards(n7);
450 slatsegB[3]->Init(0); // 0 detection element id
451 slatsegNB[3]->SetPcbBoards(n7);
452 slatsegNB[4]->Init(0); // 0 detection element id
454 // Type 223300 for 704, 709, 717, 722 in Ch7 (similar for Ch8)
455 // for the futur official numbering
456 // Type 223300 for 702, 710, 752, 760 in Ch7 (similar for Ch8)
457 // actual numbering in the code to be changed in jan05
458 Int_t n8[4] = { 0, 0, 2, 2 };
459 slatsegB[4]->SetPcbBoards(n8);
460 slatsegB[4]->Init(0); // 0 detection element id
461 slatsegNB[4]->SetPcbBoards(n8);
462 slatsegNB[4]->Init(0); // 0 detection element id
464 // Type 333000 for 705, 708, 718, 721 in Ch7 (similar for Ch8)
465 // for the futur official numbering
466 // Type 333000 for 701, 711, 751, 761 in Ch7 (similar for Ch8)
467 // actual numbering in the code to be changed in jan05
468 // Type 333000 for 906, 907, 919, 920 in Ch9 (similar for Ch10)
469 // for the futur official numbering
470 // Type 333000 for 900, 912, 950, 962 in Ch9 (similar for Ch10)
471 // actual numbering in the code to be changed in jan05
472 Int_t n9[4] = { 0, 0, 0, 3 };
473 slatsegB[5]->SetPcbBoards(n9);
474 slatsegB[5]->Init(0); // 0 detection element id
475 slatsegNB[5]->SetPcbBoards(n9);
476 slatsegNB[5]->Init(0); // 0 detection element id
478 // Type 330000 for 706, 707, 719, 720 in Ch7 (similar for Ch8)
479 // for the futur official numbering
480 // Type 330000 for 700, 712, 750, 762 in Ch7 (similar for Ch8)
481 // actual numbering in the code to be changed in jan05
482 Int_t n10[4] = { 0, 0, 0, 2 };
483 slatsegB[6]->SetPcbBoards(n10);
484 slatsegB[6]->Init(0); // 0 detection element id
485 slatsegNB[6]->SetPcbBoards(n10);
486 slatsegNB[6]->Init(0); // 0 detection element id
489 for (Int_t chamber = 6; chamber < 8; chamber++) {
491 segmentation[0] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
492 segmentation[1] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
494 // id detection elt for chamber 1
495 Int_t id0 = (chamber+1)*100;
497 //--------------------------------------------------------
498 // Configuration for Chamber TC6/7 (Station 4) ----------
501 fMUON->SetNsec(chamber,2);
505 segmentation[0]->Add(id0+ 6, slatsegB[0]);
506 segmentation[0]->Add(id0+56, slatsegB[0]);
509 segmentation[0]->Add(id0+ 5, slatsegB[1]);
510 segmentation[0]->Add(id0+ 7, slatsegB[1]);
511 segmentation[0]->Add(id0+55, slatsegB[1]);
512 segmentation[0]->Add(id0+57, slatsegB[1]);
515 segmentation[0]->Add(id0+ 4, slatsegB[2]);
516 segmentation[0]->Add(id0+ 8, slatsegB[2]);
517 segmentation[0]->Add(id0+54, slatsegB[2]);
518 segmentation[0]->Add(id0+58, slatsegB[2]);
521 segmentation[0]->Add(id0+ 3, slatsegB[3]);
522 segmentation[0]->Add(id0+ 9, slatsegB[3]);
523 segmentation[0]->Add(id0+53, slatsegB[3]);
524 segmentation[0]->Add(id0+59, slatsegB[3]);
527 segmentation[0]->Add(id0+ 2, slatsegB[4]);
528 segmentation[0]->Add(id0+10, slatsegB[4]);
529 segmentation[0]->Add(id0+52, slatsegB[4]);
530 segmentation[0]->Add(id0+60, slatsegB[4]);
533 segmentation[0]->Add(id0+ 1, slatsegB[5]);
534 segmentation[0]->Add(id0+11, slatsegB[5]);
535 segmentation[0]->Add(id0+51, slatsegB[5]);
536 segmentation[0]->Add(id0+61, slatsegB[5]);
539 segmentation[0]->Add(id0 , slatsegB[6]);
540 segmentation[0]->Add(id0+12, slatsegB[6]);
541 segmentation[0]->Add(id0+50, slatsegB[6]);
542 segmentation[0]->Add(id0+62, slatsegB[6]);
543 fMUON->SetSegmentationModel(chamber, 1, segmentation[0]);
547 segmentation[1]->Add(id0+ 6, slatsegNB[0]);
548 segmentation[1]->Add(id0+56, slatsegNB[0]);
551 segmentation[1]->Add(id0+ 5, slatsegNB[1]);
552 segmentation[1]->Add(id0+ 7, slatsegNB[1]);
553 segmentation[1]->Add(id0+55, slatsegNB[1]);
554 segmentation[1]->Add(id0+57, slatsegNB[1]);
557 segmentation[1]->Add(id0+ 4, slatsegNB[2]);
558 segmentation[1]->Add(id0+ 8, slatsegNB[2]);
559 segmentation[1]->Add(id0+54, slatsegNB[2]);
560 segmentation[1]->Add(id0+58, slatsegNB[2]);
563 segmentation[1]->Add(id0+ 3, slatsegNB[3]);
564 segmentation[1]->Add(id0+ 9, slatsegNB[3]);
565 segmentation[1]->Add(id0+53, slatsegNB[3]);
566 segmentation[1]->Add(id0+59, slatsegNB[3]);
569 segmentation[1]->Add(id0+ 2, slatsegNB[4]);
570 segmentation[1]->Add(id0+10, slatsegNB[4]);
571 segmentation[1]->Add(id0+52, slatsegNB[4]);
572 segmentation[1]->Add(id0+60, slatsegNB[4]);
575 segmentation[1]->Add(id0+ 1, slatsegNB[5]);
576 segmentation[1]->Add(id0+11, slatsegNB[5]);
577 segmentation[1]->Add(id0+51, slatsegNB[5]);
578 segmentation[1]->Add(id0+61, slatsegNB[5]);
581 segmentation[1]->Add(id0 , slatsegNB[6]);
582 segmentation[1]->Add(id0+12, slatsegNB[6]);
583 segmentation[1]->Add(id0+50, slatsegNB[6]);
584 segmentation[1]->Add(id0+62, slatsegNB[6]);
585 fMUON->SetSegmentationModel(chamber, 2, segmentation[1]);
587 fMUON->SetResponseModel(chamber, fResponse0);
588 fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread
592 //__________________________________________________________________________
593 void AliMUONFactoryV2::BuildStation5()
595 //--------------------------------------------------------
596 // Configuration for Chamber TC9/10 (Station 5) ---------
597 //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
599 AliMUONGeometrySegmentation* segmentation[2];
601 //Slats Segmentations
602 AliMUONSt345SlatSegmentation *slatsegB[6]; // Types of segmentation for St5
603 AliMUONSt345SlatSegmentation *slatsegNB[6];
606 Int_t ndiv[4] ={ 4, 4, 2, 1}; // densities zones
607 for(Int_t i = 0; i < 6; i++) {
608 slatsegB[i] = new AliMUONSt345SlatSegmentation(1);
609 fDESegmentations->Add(slatsegB[i]);
610 slatsegB[i]->SetPadSize(10.,0.5);
611 slatsegB[i]->SetPadDivision(ndiv);
612 slatsegB[i]->SetId(1);
613 slatsegB[i]->SetDAnod(AliMUONConstants::Pitch());
614 slatsegNB[i] = new AliMUONSt345SlatSegmentation(0);
615 fDESegmentations->Add(slatsegNB[i]);
616 slatsegNB[i]->SetPadSize(0.713,10.);
617 slatsegNB[i]->SetPadDivision(ndiv);
618 slatsegNB[i]->SetId(1);
619 slatsegNB[i]->SetDAnod(AliMUONConstants::Pitch());
622 // Type 122330 for 900, 913 in Ch9 (similar for Ch10)
623 // for the futur official numbering
624 // Type 122330 for 906, 956 in Ch9 (similar for Ch10)
625 // actual numbering in the code to be changed in jan05
626 Int_t n4[4] = { 0, 1, 2, 2 };
627 slatsegB[0]->SetPcbBoards(n4);
628 slatsegB[0]->Init(0); // 0 detection element id
629 slatsegNB[0]->SetPcbBoards(n4);
630 slatsegNB[0]->Init(0); // 0 detection element id
632 // Type 112233 for 901, 902, 911, 912, 914, 915, 924, 925 in Ch9
633 // (similar for Ch10) for the futur official numbering
634 // Type 112233 for 904, 905, 907, 908, 954, 955, 957, 958 in Ch9
635 // (similar for Ch10) actual numbering in the code to be changed in jan05
636 Int_t n5[4] = { 0, 2, 2, 2 };
637 slatsegB[1]->SetPcbBoards(n5);
638 slatsegB[1]->Init(0); // 0 detection element id
639 slatsegNB[1]->SetPcbBoards(n5);
640 slatsegNB[1]->Init(0); // 0 detection element id
642 // Type 333000 for 906, 907, 919, 920 in Ch9 (similar for Ch10)
643 // for the futur official numbering
644 // Type 333000 for 900, 912, 950, 962 in Ch9 (similar for Ch10)
645 // actual numbering in the code to be changed in jan05
646 Int_t n9[4] = { 0, 0, 0, 3 };
647 slatsegB[2]->SetPcbBoards(n9);
648 slatsegB[2]->Init(0); // 0 detection element id
649 slatsegNB[2]->SetPcbBoards(n9);
650 slatsegNB[2]->Init(0); // 0 detection element id
652 // Type 222333 for 903, 910, 916, 923 in Ch9 (similar for Ch10)
653 // for the futur official numbering
654 // Type 222333 for 903, 909, 953, 959 in Ch9 (similar for Ch10)
655 // actual numbering in the code to be changed in jan05
656 Int_t n11[4] = { 0, 0, 3, 3 };
657 slatsegB[3]->SetPcbBoards(n11);
658 slatsegB[3]->Init(0); // 0 detection element id
659 slatsegNB[3]->SetPcbBoards(n11);
660 slatsegNB[3]->Init(0); // 0 detection element id
662 // Type 223330 for 904, 909, 917, 922 in Ch9 (similar for Ch10)
663 // for the futur official numbering
664 // Type 223330 for 902, 910, 952, 960 in Ch9 (similar for Ch10)
665 // actual numbering in the code to be changed in jan05
666 Int_t n12[4] = { 0, 0, 2, 3 };
667 slatsegB[4]->SetPcbBoards(n12);
668 slatsegB[4]->Init(0); // 0 detection element id
669 slatsegNB[4]->SetPcbBoards(n12);
670 slatsegNB[4]->Init(0); // 0 detection element id
672 // Type 333300 for 905, 908, 918, 921 in Ch9 (similar for Ch10)
673 // for the futur official numbering
674 // Type 333300 for 901, 911, 951, 961 in Ch9 (similar for Ch10)
675 // actual numbering in the code to be changed in jan05
676 Int_t n13[4] = { 0, 0, 0, 4 };
677 slatsegB[5]->SetPcbBoards(n13);
678 slatsegB[5]->Init(0); // 0 detection element id
679 slatsegNB[5]->SetPcbBoards(n13);
680 slatsegNB[5]->Init(0); // 0 detection element id
682 for (Int_t chamber = 8; chamber < 10; chamber++) {
684 segmentation[0] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
685 segmentation[1] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
687 // id detection elt for chamber 1
688 Int_t id0 = (chamber+1)*100;
690 //--------------------------------------------------------
691 // Configuration for Chamber TC8/9 (Station 5) ----------
694 fMUON->SetNsec(chamber,2);
698 segmentation[0]->Add(id0+ 6, slatsegB[0]);
699 segmentation[0]->Add(id0+56, slatsegB[0]);
702 segmentation[0]->Add(id0+ 4, slatsegB[1]);
703 segmentation[0]->Add(id0+ 5, slatsegB[1]);
704 segmentation[0]->Add(id0+ 7, slatsegB[1]);
705 segmentation[0]->Add(id0+ 8, slatsegB[1]);
706 segmentation[0]->Add(id0+54, slatsegB[1]);
707 segmentation[0]->Add(id0+55, slatsegB[1]);
708 segmentation[0]->Add(id0+57, slatsegB[1]);
709 segmentation[0]->Add(id0+58, slatsegB[1]);
712 segmentation[0]->Add(id0 , slatsegB[2]);
713 segmentation[0]->Add(id0+12, slatsegB[2]);
714 segmentation[0]->Add(id0+50, slatsegB[2]);
715 segmentation[0]->Add(id0+62, slatsegB[2]);
718 segmentation[0]->Add(id0+ 3, slatsegB[3]);
719 segmentation[0]->Add(id0+ 9, slatsegB[3]);
720 segmentation[0]->Add(id0+53, slatsegB[3]);
721 segmentation[0]->Add(id0+59, slatsegB[3]);
724 segmentation[0]->Add(id0+ 2, slatsegB[4]);
725 segmentation[0]->Add(id0+10, slatsegB[4]);
726 segmentation[0]->Add(id0+52, slatsegB[4]);
727 segmentation[0]->Add(id0+60, slatsegB[4]);
730 segmentation[0]->Add(id0+ 1, slatsegB[5]);
731 segmentation[0]->Add(id0+11, slatsegB[5]);
732 segmentation[0]->Add(id0+51, slatsegB[5]);
733 segmentation[0]->Add(id0+61, slatsegB[5]);
734 fMUON->SetSegmentationModel(chamber, 1, segmentation[0]);
738 segmentation[1]->Add(id0+ 6, slatsegNB[0]);
739 segmentation[1]->Add(id0+56, slatsegNB[0]);
742 segmentation[1]->Add(id0+ 4, slatsegNB[1]);
743 segmentation[1]->Add(id0+ 5, slatsegNB[1]);
744 segmentation[1]->Add(id0+ 7, slatsegNB[1]);
745 segmentation[1]->Add(id0+ 8, slatsegNB[1]);
746 segmentation[1]->Add(id0+54, slatsegNB[1]);
747 segmentation[1]->Add(id0+55, slatsegNB[1]);
748 segmentation[1]->Add(id0+57, slatsegNB[1]);
749 segmentation[1]->Add(id0+58, slatsegNB[1]);
752 segmentation[1]->Add(id0 , slatsegNB[2]);
753 segmentation[1]->Add(id0+12, slatsegNB[2]);
754 segmentation[1]->Add(id0+50, slatsegNB[2]);
755 segmentation[1]->Add(id0+62, slatsegNB[2]);
758 segmentation[1]->Add(id0+ 3, slatsegNB[3]);
759 segmentation[1]->Add(id0+ 9, slatsegNB[3]);
760 segmentation[1]->Add(id0+53, slatsegNB[3]);
761 segmentation[1]->Add(id0+59, slatsegNB[3]);
764 segmentation[1]->Add(id0+ 2, slatsegNB[4]);
765 segmentation[1]->Add(id0+10, slatsegNB[4]);
766 segmentation[1]->Add(id0+52, slatsegNB[4]);
767 segmentation[1]->Add(id0+60, slatsegNB[4]);
770 segmentation[1]->Add(id0+ 1, slatsegNB[5]);
771 segmentation[1]->Add(id0+11, slatsegNB[5]);
772 segmentation[1]->Add(id0+51, slatsegNB[5]);
773 segmentation[1]->Add(id0+61, slatsegNB[5]);
774 fMUON->SetSegmentationModel(chamber, 2, segmentation[1]);
776 fMUON->SetResponseModel(chamber, fResponse0);
777 fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread
781 //__________________________________________________________________________
782 void AliMUONFactoryV2::BuildStation6()
784 // Create Trigger geometry segmentation for given chamber and cathod
787 AliMUONGeometrySegmentation *chamberSeg[2];
789 AliMUONResponseTrigger* responseTrigger0 = new AliMUONResponseTrigger;
791 // AliMUONResponseTriggerV1* responseTrigger0 = new AliMUONResponseTriggerV1;
793 for (Int_t chamber = 10; chamber < 14; chamber++) {
795 //Trigger Segmentation
796 AliMUONTriggerSegmentation *trigSegX[9];
797 AliMUONTriggerSegmentation *trigSegY[9];
798 for(Int_t i=0; i<9; i++) {
799 trigSegX[i] = new AliMUONTriggerSegmentation(1);
800 trigSegY[i] = new AliMUONTriggerSegmentation(0);
801 fDESegmentations->Add(trigSegX[i]);
802 fDESegmentations->Add(trigSegY[i]);
803 trigSegX[i]->SetLineNumber(9-i);
804 trigSegY[i]->SetLineNumber(9-i);
807 AliMUONChamber *iChamber, *iChamber1;
808 iChamber1 = &fMUON->Chamber(10);
809 iChamber = &fMUON->Chamber(chamber);
810 Float_t zpos1= - iChamber1->Z();
811 Float_t zpos = - iChamber->Z();
812 Float_t zRatio = zpos / zpos1;
815 Float_t stripWidth[3]={0.,0.,0.}; // 1.0625 2.125 4.25
816 Float_t stripLength[4]={0.,0.,0.,0.}; // 17. 34. 51. 68.
817 for (Int_t i=0; i<3; i++)
818 stripWidth[i]=AliMUONTriggerConstants::StripWidth(i)*zRatio;
819 for (Int_t i=0; i<4; i++)
820 stripLength[i]=AliMUONTriggerConstants::StripLength(i)*zRatio;
821 Int_t nStrip[7]={0,0,0,0,0,0,0};
822 Float_t stripYsize[7]={0.,0.,0.,0.,0.,0.,0.};
823 Float_t stripXsize[7]={0.,0.,0.,0.,0.,0.,0.};
825 // chamber 8 0 cathode 0
826 for (Int_t i=0; i<7; i++) nStrip[i]=16;
827 for (Int_t i=0; i<7; i++) stripYsize[i]=stripWidth[2];
828 for (Int_t i=0; i<6; i++) stripXsize[i]=stripLength[1];
829 stripXsize[6]=stripLength[2];
830 trigSegX[8]->Init(0,nStrip,stripYsize,stripXsize,0.);
831 trigSegX[0]->Init(0,nStrip,stripYsize,stripXsize,0.);
833 // chamber 8 7 1 0 cathode 1
834 for (Int_t i=0; i<6; i++) nStrip[i]=8;
836 for (Int_t i=0; i<7; i++) stripYsize[i]=stripLength[3];
837 for (Int_t i=0; i<7; i++) stripXsize[i]=stripWidth[2];
838 trigSegY[8]->Init(0,nStrip,stripYsize,stripXsize,0.);
839 trigSegY[7]->Init(0,nStrip,stripYsize,stripXsize,0.);
840 trigSegY[1]->Init(0,nStrip,stripYsize,stripXsize,0.);
841 trigSegY[0]->Init(0,nStrip,stripYsize,stripXsize,0.);
843 // chamber 7 6 2 1 cathode 0
844 for (Int_t i=0; i<6; i++) nStrip[i]=32;
846 for (Int_t i=0; i<6; i++) stripYsize[i]=stripWidth[1];
847 stripYsize[6]=stripWidth[2];
848 for (Int_t i=0; i<6; i++) stripXsize[i]=stripLength[1];
849 stripXsize[6]=stripLength[2];
850 trigSegX[7]->Init(0,nStrip,stripYsize,stripXsize,0.);
851 trigSegX[6]->Init(0,nStrip,stripYsize,stripXsize,0.);
852 trigSegX[2]->Init(0,nStrip,stripYsize,stripXsize,0.);
853 trigSegX[1]->Init(0,nStrip,stripYsize,stripXsize,0.);
855 // chamber 6 2 cathode 1
856 for (Int_t i=0; i<5; i++) nStrip[i]=16;
857 for (Int_t i=5; i<6; i++) nStrip[i]=8;
859 for (Int_t i=0; i<7; i++) stripYsize[i]=stripLength[3];
860 for (Int_t i=0; i<5; i++) stripXsize[i]=stripWidth[1];
861 for (Int_t i=5; i<7; i++) stripXsize[i]=stripWidth[2];
862 trigSegY[6]->Init(0,nStrip,stripYsize,stripXsize,0.);
863 trigSegY[2]->Init(0,nStrip,stripYsize,stripXsize,0.);
865 // chamber 5 3 cathode 0
867 for (Int_t i=1; i<3; i++) nStrip[i]=64;
868 for (Int_t i=3; i<6; i++) nStrip[i]=32;
870 for (Int_t i=0; i<3; i++) stripYsize[i]=stripWidth[0];
871 for (Int_t i=3; i<6; i++) stripYsize[i]=stripWidth[1];
872 stripYsize[6]=stripWidth[2];
873 for (Int_t i=0; i<6; i++) stripXsize[i]=stripLength[1];
874 stripXsize[6]=stripLength[2];
875 trigSegX[5]->Init(0,nStrip,stripYsize,stripXsize,stripLength[0]);
876 trigSegX[3]->Init(0,nStrip,stripYsize,stripXsize,0.);
878 // chamber 5 3 cathode 1
879 for (Int_t i=0; i<5; i++) nStrip[i]=16;
880 for (Int_t i=5; i<6; i++) nStrip[5]=8;
882 stripYsize[0]=stripLength[2];
883 for (Int_t i=1; i<8; i++) stripYsize[i]=stripLength[3];
884 for (Int_t i=0; i<5; i++) stripXsize[i]=stripWidth[1];
885 for (Int_t i=5; i<7; i++) stripXsize[i]=stripWidth[2];
886 trigSegY[5]->Init(0,nStrip,stripYsize,stripXsize,stripLength[0]);
887 trigSegY[3]->Init(0,nStrip,stripYsize,stripXsize,0.);
889 // chamber 4 cathode 0
891 for (Int_t i=1; i<3; i++) nStrip[i]=64;
892 for (Int_t i=3; i<6; i++) nStrip[i]=32;
895 for (Int_t i=1; i<3; i++) stripYsize[i]=stripWidth[0];
896 for (Int_t i=3; i<6; i++) stripYsize[i]=stripWidth[1];
897 stripYsize[6]=stripWidth[2];
899 stripXsize[1]=stripLength[0];
900 for (Int_t i=2; i<6; i++) stripXsize[i]=stripLength[1];
901 stripXsize[6]=stripLength[2];
902 trigSegX[4]->Init(0,nStrip,stripYsize,stripXsize,0.);
904 // chamber 4 cathode 1
907 for (Int_t i=2; i<5; i++) nStrip[i]=16;
908 for (Int_t i=5; i<6; i++) nStrip[i]=8;
911 for (Int_t i=1; i<7; i++) stripYsize[i]=stripLength[3];
913 for (Int_t i=1; i<5; i++) stripXsize[i]=stripWidth[1];
914 for (Int_t i=5; i<7; i++) stripXsize[i]=stripWidth[2];
915 trigSegY[4]->Init(0,nStrip,stripYsize,stripXsize,0.);
917 chamberSeg[0] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
918 chamberSeg[1] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
920 fMUON->SetNsec(chamber,2);
921 Int_t icount=chamber-10; // chamber counter (0 1 2 3)
922 Int_t id0=(10+icount+1)*100;
925 // printf("in CreateTriggerSegmentation here 0 id0=%i \n",id0);
927 for (Int_t i = 0; i < 9; i++) {
930 chamberSeg[0]->Add(id0+i, trigSegX[i]);
931 chamberSeg[0]->Add(id0+50+i, trigSegX[i]);
932 fMUON->SetSegmentationModel(chamber, 1, chamberSeg[0]);
935 chamberSeg[1]->Add(id0+i, trigSegY[i]);
936 chamberSeg[1]->Add(id0+50+i, trigSegY[i]);
937 fMUON->SetSegmentationModel(chamber, 2, chamberSeg[1]);
941 fMUON->SetResponseModel(chamber, responseTrigger0);
942 fMUON->Chamber(chamber).SetChargeCorrel(0); // same charge on cathodes
945 // printf("in CreateTriggerSegmentation here 1\n");
948 AliWarning(Form("Segmentation for chamber %d is not yet defined",chamber));
953 //__________________________________________________________________________
954 void AliMUONFactoryV2::Build(AliMUON* where, const char* what)
957 // Construct MUON from chambers, segmentation and responses
964 if (strcmp(tmp, "default")==0) {
965 // Set default parameters
967 fMUON->SetMaxStepGas(0.1);
968 fMUON->SetMaxStepAlu(0.1);
970 // Build all stations
980 AliDebug(0,"Non default version of MUON selected. You have to construct yourself the MUON elements !!");
983 //__________________________________________________________________________
984 void AliMUONFactoryV2::BuildStation(AliMUON* where, Int_t stationNumber)
987 // Construct MUON from chambers, segmentation and responses
991 // First define the number of planes that are segmented (1 or 2) by a call
993 // Then chose for each chamber (chamber plane) the segmentation
994 // and response model.
995 // They should be equal for the two chambers of each station. In a future
996 // version this will be enforced.
1000 if (!fResponse0) BuildCommon();
1002 switch (stationNumber) {
1003 case 1: BuildStation1(); break;
1004 case 2: BuildStation2(); break;
1005 case 3: BuildStation3(); break;
1006 case 4: BuildStation4(); break;
1007 case 5: BuildStation5(); break;
1008 case 6: BuildStation6(); break;
1010 default: AliFatal("Wrong station number");