a4a39417c8e7901168fbda21fcd7facddca63fb0
[u/mrichter/AliRoot.git] / MUON / AliMUONFactoryV2.cxx
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 ////////////////////////////////////////////////////////////
17 //  Factory for muon chambers, segmentations and response //
18 ////////////////////////////////////////////////////////////
19
20 /* $Id$ */
21
22 #include "AliMUONFactoryV2.h"
23 #include "AliRun.h"
24 #include "AliSegmentation.h"
25 #include "AliLog.h"
26
27 #include "AliMUON.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 "AliMUONVGeometryDEIndexing.h"
36 #include "AliMUONSt12QuadrantSegmentation.h"
37 #include "AliMUONSt345SlatSegmentation.h"
38 #include "AliMUONTriggerSegmentation.h"
39 #include "AliMUONResponseTrigger.h"
40
41 ClassImp(AliMUONFactoryV2)
42
43 //__________________________________________________________________________
44   AliMUONFactoryV2::AliMUONFactoryV2(const char* name)
45     : AliMUONFactory(name),
46       fMUON(0),
47       fResponse0(0),
48       fDESegmentations(0)
49 {
50   // FactoryV2 inherite from Factory for switching in AliMUONv1::Init()
51   // to be changed when old segmentation will be removed.
52
53   fDESegmentations = new TObjArray();
54   fDESegmentations->SetOwner(kTRUE);
55 }
56
57 //__________________________________________________________________________
58   AliMUONFactoryV2::AliMUONFactoryV2()
59     : AliMUONFactory(),
60       fMUON(0),
61       fResponse0(0),
62       fDESegmentations(0)
63 {
64 // Default constructor
65 }
66
67 //__________________________________________________________________________
68 AliMUONFactoryV2::AliMUONFactoryV2(const AliMUONFactoryV2& rhs)
69   : AliMUONFactory(rhs)
70 {
71   // Protected copy constructor
72
73   AliFatal("Not implemented.");
74 }
75
76 //__________________________________________________________________________
77
78 AliMUONFactoryV2::~AliMUONFactoryV2()
79 {
80 // Destructor
81
82   delete fDESegmentations;
83 }
84
85 //__________________________________________________________________________
86 AliMUONFactoryV2&  AliMUONFactoryV2::operator=(const AliMUONFactoryV2& rhs)
87 {
88   // Protected assignement operator
89
90   if (this == &rhs) return *this;
91
92   AliFatal("Not implemented.");
93     
94   return *this;  
95 }    
96           
97 //__________________________________________________________________________
98 Bool_t AliMUONFactoryV2::IsGeometryDefined(Int_t ichamber)
99 {
100 // Return true, if det elements for the chamber with the given ichamber Id
101 // are defined in geometry (the geometry builder for this chamber was activated)
102
103   if ( ! fMUON ||
104        ! fMUON->Chamber(ichamber).GetGeometry() ||
105        ! fMUON->Chamber(ichamber).GetGeometry()->GetDEIndexing() ||
106        ! fMUON->Chamber(ichamber).GetGeometry()->GetDEIndexing()->GetNofDetElements() )
107        
108     return kFALSE;
109   
110   return kTRUE;
111 }  
112
113 //__________________________________________________________________________
114 void AliMUONFactoryV2::BuildCommon() 
115 {
116   //
117   // Construct the default response.
118   //
119
120   // Default response: 5 mm of gas
121   fResponse0 = new AliMUONResponseV0;
122   fResponse0->SetSqrtKx3AndDeriveKx2Kx4(0.7131); // sqrt(0.5085)
123   fResponse0->SetSqrtKy3AndDeriveKy2Ky4(0.7642); // sqrt(0.5840)
124   fResponse0->SetPitch(AliMUONConstants::Pitch()); // anode-cathode distance
125   fResponse0->SetSigmaIntegration(10.);
126   fResponse0->SetChargeSlope(10);
127   fResponse0->SetChargeSpread(0.18, 0.18);
128   fResponse0->SetMaxAdc(4096);
129   fResponse0->SetSaturation(3000);
130   fResponse0->SetZeroSuppression(6);
131 }       
132         
133 //__________________________________________________________________________
134 void AliMUONFactoryV2::BuildStation1() 
135 {
136   //--------------------------------------------------------
137   // Configuration for Chamber TC1/2  (Station 1) ----------           
138   //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
139
140
141   // Response for 4 mm of gas (station 1)
142   // automatic consistency with width of sensitive medium in CreateGeometry ????
143   AliMUONResponseV0* responseSt1 = new AliMUONResponseV0;
144   // Mathieson parameters from L.Kharmandarian's thesis, page 190
145   responseSt1->SetSqrtKx3AndDeriveKx2Kx4(0.7000); // sqrt(0.4900)
146   responseSt1->SetSqrtKy3AndDeriveKy2Ky4(0.7550); // sqrt(0.5700)
147   responseSt1->SetPitch(AliMUONConstants::PitchSt1()); // anode-cathode distance
148   responseSt1->SetSigmaIntegration(10.);
149   // ChargeSlope larger to compensate for the smaller anode-cathode distance
150   // and keep the same most probable ADC channel for mip's
151   responseSt1->SetChargeSlope(62.5); 
152   // assumed proportionality to anode-cathode distance for ChargeSpread
153   responseSt1->SetChargeSpread(0.144, 0.144);
154   responseSt1->SetMaxAdc(4096);
155   responseSt1->SetSaturation(3000);
156   responseSt1->SetZeroSuppression(6);
157
158   // Quadrant segmentations:
159   AliMUONSt12QuadrantSegmentation* bendSt1
160     = new AliMUONSt12QuadrantSegmentation(kStation1, kBendingPlane);
161   AliMUONSt12QuadrantSegmentation* nonbendSt1
162     = new AliMUONSt12QuadrantSegmentation(kStation1, kNonBendingPlane);
163   
164   // Add in the array (for safe deleting)  
165   fDESegmentations->Add(bendSt1);  
166   fDESegmentations->Add(nonbendSt1);  
167
168   AliMUONGeometrySegmentation* segmentation[2];
169
170   for (Int_t chamber = 0; chamber < 2; chamber++) {
171
172     segmentation[0] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
173     segmentation[1] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
174         
175     // id detection elt for chamber 1
176     Int_t id0 = (chamber+1)*100;
177
178     //--------------------------------------------------------
179     // Configuration for Chamber TC1/2  (Station 1) ----------           
180
181
182     fMUON->SetNsec(chamber,2);
183
184     // cathode 0
185     segmentation[0]->Add(id0,      bendSt1);
186     segmentation[0]->Add(id0 +  1, nonbendSt1);
187     segmentation[0]->Add(id0 + 50, bendSt1);
188     segmentation[0]->Add(id0 + 51, nonbendSt1); 
189     fMUON->SetSegmentationModel(chamber, 1, segmentation[0]);   
190
191     // cathode 1
192     segmentation[1]->Add(id0,      nonbendSt1);
193     segmentation[1]->Add(id0 +  1, bendSt1);
194     segmentation[1]->Add(id0 + 50, nonbendSt1);
195     segmentation[1]->Add(id0 + 51, bendSt1);
196     fMUON->SetSegmentationModel(chamber, 2, segmentation[1]);
197         
198     fMUON->SetResponseModel(chamber, responseSt1); // special response      
199     fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread
200         
201   }
202 }
203
204 //__________________________________________________________________________
205 void AliMUONFactoryV2::BuildStation2() 
206 {
207   //
208   //--------------------------------------------------------
209   // Configuration for Chamber TC3/4 (Station 2) -----------
210   ///^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
211
212
213   // Quadrant segmentations:
214   AliMUONSt12QuadrantSegmentation* bendSt2
215     = new AliMUONSt12QuadrantSegmentation(kStation2, kBendingPlane);
216   AliMUONSt12QuadrantSegmentation* nonbendSt2
217     = new AliMUONSt12QuadrantSegmentation(kStation2, kNonBendingPlane);
218
219   // Add in the array (for safe deleting)  
220   fDESegmentations->Add(bendSt2);  
221   fDESegmentations->Add(nonbendSt2);  
222
223   AliMUONGeometrySegmentation* segmentation[2];
224
225   for (Int_t chamber = 2; chamber < 4; chamber++) {
226
227     segmentation[0] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
228     segmentation[1] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
229         
230     // id detection elt for chamber 1
231     Int_t id0 = (chamber+1)*100;
232
233     //--------------------------------------------------------
234     // Configuration for Chamber TC3/4  (Station 2) ----------           
235
236
237     fMUON->SetNsec(chamber,2);
238
239     // cathode 0
240     segmentation[0]->Add(id0,      bendSt2);
241     segmentation[0]->Add(id0 +  1, nonbendSt2);
242     segmentation[0]->Add(id0 + 50, bendSt2);
243     segmentation[0]->Add(id0 + 51, nonbendSt2); 
244     fMUON->SetSegmentationModel(chamber, 1, segmentation[0]);   
245
246     // cathode 1
247     segmentation[1]->Add(id0,      nonbendSt2);
248     segmentation[1]->Add(id0 +  1, bendSt2);
249     segmentation[1]->Add(id0 + 50, nonbendSt2);
250     segmentation[1]->Add(id0 + 51, bendSt2);
251     fMUON->SetSegmentationModel(chamber, 2, segmentation[1]);
252         
253     fMUON->SetResponseModel(chamber, fResponse0); // normal response        
254     fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread
255         
256   }
257 }       
258         
259         
260 //__________________________________________________________________________
261 void AliMUONFactoryV2::BuildStation3() 
262 {
263   //--------------------------------------------------------
264   // Configuration for Chamber TC5/6  (Station 3) ----------          
265   //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
266
267   AliMUONGeometrySegmentation* segmentation[2];
268
269   //Slats Segmentations
270   AliMUONSt345SlatSegmentation *slatsegB[4]; // Types of segmentation for St3
271   AliMUONSt345SlatSegmentation *slatsegNB[4]; 
272   // Bending
273
274   Int_t ndiv[4] ={ 4, 4, 2, 1};  // densities zones 
275   for(Int_t i=0; i<4; i++) {
276     slatsegB[i] = new AliMUONSt345SlatSegmentation(1);
277     fDESegmentations->Add(slatsegB[i]);  
278     slatsegB[i]->SetPadSize(10.,0.5);
279     slatsegB[i]->SetPadDivision(ndiv);
280     slatsegB[i]->SetId(1); // Id elt ????
281     slatsegB[i]->SetDAnod(AliMUONConstants::Pitch());
282     slatsegNB[i] = new AliMUONSt345SlatSegmentation(0);
283     fDESegmentations->Add(slatsegNB[i]);  
284     slatsegNB[i]->SetPadSize(1./1.4,10.); // Nbending
285     slatsegNB[i]->SetPadDivision(ndiv);
286     slatsegNB[i]->SetId(1);
287     slatsegNB[i]->SetDAnod(AliMUONConstants::Pitch());
288   }
289
290   // Type 112200 for 500, 501, 508, 509, 510, 517 
291   // in Ch5 (similar for Ch6) for the futur official numbering
292   // Type 112200 for 503, 504, 505, 555, 554, 553 
293   // in Ch5 (similar for Ch6) actual numbering in the code to be changed in jan05
294   Int_t n0[4] = { 0, 2, 2, 0 };
295   slatsegB[0]->SetPcbBoards(n0);
296   slatsegB[0]->Init(0);
297   slatsegNB[0]->SetPcbBoards(n0);
298   slatsegNB[0]->Init(0);
299     
300   // Type 122200 for 502, 507, 511, 516 (similar in Ch6) 
301   // for future official numbering of ALICE
302   // Type 122200 for 502, 506, 556, 552 (similiarin Ch6) 
303   // for actual numbering in muon code to be changed in jan05
304   Int_t n1[4] = { 0, 1, 3, 0 }; 
305   slatsegB[1]->SetPcbBoards(n1);
306   slatsegB[1]->Init(0); 
307   slatsegNB[1]->SetPcbBoards(n1);
308   slatsegNB[1]->Init(0); 
309     
310   // Type 222000 for 503, 506, 512, 515 (similar in Ch6) 
311   // for future official numbering of ALICE
312   // Type 222000 for 501, 507, 557, 551 (similiarin Ch6) 
313   // for actual numbering in muon code to be changed in jan05
314   Int_t n2[4] = { 0, 0, 3, 0 };
315   slatsegB[2]->SetPcbBoards(n2);
316   slatsegB[2]->Init(0);
317   slatsegNB[2]->SetPcbBoards(n2);
318   slatsegNB[2]->Init(0);
319     
320   // Type 220000 for 504, 505, 513, 514 (similar in Ch6) 
321   // for future official numbering of ALICE
322   // Type 220000 for 500, 508, 558, 550 (similiarin Ch6) 
323   // for actual numbering in muon code to be changed in jan05
324   Int_t n3[4] = { 0, 0, 2, 0 };
325   slatsegB[3]->SetPcbBoards(n3);
326   slatsegB[3]->Init(0); 
327   slatsegNB[3]->SetPcbBoards(n3);
328   slatsegNB[3]->Init(0); 
329
330   for (Int_t chamber = 4; chamber < 6; chamber++) {
331
332     segmentation[0] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
333     segmentation[1] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
334         
335     // id detection elt for chamber 1
336     Int_t id0 = (chamber+1)*100;
337
338     //--------------------------------------------------------
339     // Configuration for Chamber TC3/4  (Station 2) ----------           
340
341
342     fMUON->SetNsec(chamber,2);
343
344     // cathode 0
345     // type 220000
346     segmentation[0]->Add(id0, slatsegB[3]);
347     segmentation[0]->Add(id0+ 8, slatsegB[3]);  
348     segmentation[0]->Add(id0+50, slatsegB[3]);  
349     segmentation[0]->Add(id0+58, slatsegB[3]);
350     // type 222000
351     segmentation[0]->Add(id0+ 1, slatsegB[2]);
352     segmentation[0]->Add(id0+ 7, slatsegB[2]);  
353     segmentation[0]->Add(id0+51, slatsegB[2]);  
354     segmentation[0]->Add(id0+57, slatsegB[2]);
355     // type 122200
356     segmentation[0]->Add(id0+ 2, slatsegB[1]);
357     segmentation[0]->Add(id0+ 6, slatsegB[1]);  
358     segmentation[0]->Add(id0+52, slatsegB[1]);  
359     segmentation[0]->Add(id0+56, slatsegB[1]);
360     // type 112200
361     segmentation[0]->Add(id0+ 3, slatsegB[0]);
362     segmentation[0]->Add(id0+ 4, slatsegB[0]);  
363     segmentation[0]->Add(id0+ 5, slatsegB[0]);  
364     segmentation[0]->Add(id0+53, slatsegB[0]);
365     segmentation[0]->Add(id0+54, slatsegB[0]);     
366     segmentation[0]->Add(id0+55, slatsegB[0]);
367     fMUON->SetSegmentationModel(chamber, 1, segmentation[0]);   
368
369     // cathode 1
370     // type 220000
371     segmentation[1]->Add(id0, slatsegNB[3]);
372     segmentation[1]->Add(id0+ 8, slatsegNB[3]);  
373     segmentation[1]->Add(id0+50, slatsegNB[3]);  
374     segmentation[1]->Add(id0+58, slatsegNB[3]);
375     // type 222000
376     segmentation[1]->Add(id0+ 1, slatsegNB[2]);
377     segmentation[1]->Add(id0+ 7, slatsegNB[2]);  
378     segmentation[1]->Add(id0+51, slatsegNB[2]);  
379     segmentation[1]->Add(id0+57, slatsegNB[2]);
380     // type 122200
381     segmentation[1]->Add(id0+ 2, slatsegNB[1]);
382     segmentation[1]->Add(id0+ 6, slatsegNB[1]);  
383     segmentation[1]->Add(id0+52, slatsegNB[1]);  
384     segmentation[1]->Add(id0+56, slatsegNB[1]);
385     // type 112200
386     segmentation[1]->Add(id0+ 3, slatsegNB[0]);
387     segmentation[1]->Add(id0+ 4, slatsegNB[0]);  
388     segmentation[1]->Add(id0+ 5, slatsegNB[0]);  
389     segmentation[1]->Add(id0+53, slatsegNB[0]);
390     segmentation[1]->Add(id0+54, slatsegNB[0]);     
391     segmentation[1]->Add(id0+55, slatsegNB[0]);
392     fMUON->SetSegmentationModel(chamber, 2, segmentation[1]);
393
394     fMUON->SetResponseModel(chamber, fResponse0);      
395     fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread
396   }
397 }
398
399         
400 //__________________________________________________________________________
401 void AliMUONFactoryV2::BuildStation4() 
402 {
403   //--------------------------------------------------------
404   // Configuration for Chamber TC7/8  (Station 4) ----------           
405   //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
406
407
408   AliMUONGeometrySegmentation* segmentation[2];
409
410   //Slats Segmentations
411   AliMUONSt345SlatSegmentation *slatsegB[7]; // Types of segmentation for St4
412   AliMUONSt345SlatSegmentation *slatsegNB[7]; 
413   // Bending
414
415   Int_t ndiv[4] ={ 4, 4, 2, 1};  // densities zones 
416   for(Int_t i = 0; i < 7; i++) {
417     slatsegB[i] = new AliMUONSt345SlatSegmentation(1);
418     fDESegmentations->Add(slatsegB[i]);  
419     slatsegB[i]->SetPadSize(10.,0.5);
420     slatsegB[i]->SetPadDivision(ndiv);
421     slatsegB[i]->SetId(1);
422     slatsegB[i]->SetDAnod(AliMUONConstants::Pitch());
423     slatsegNB[i] = new AliMUONSt345SlatSegmentation(0);
424     fDESegmentations->Add(slatsegNB[i]);  
425     slatsegNB[i]->SetPadSize(1./1.4,10.); 
426     slatsegNB[i]->SetPadDivision(ndiv);
427     slatsegNB[i]->SetId(1);
428     slatsegNB[i]->SetDAnod(AliMUONConstants::Pitch());
429   }
430
431   Int_t n4[4] = { 0, 1, 2, 2 };
432   slatsegB[0]->SetPcbBoards(n4);
433   slatsegB[0]->Init(0); // 0 detection element id
434   slatsegNB[0]->SetPcbBoards(n4);
435   slatsegNB[0]->Init(0); // 0 detection element id
436     
437   // Type 112233 for 701, 712, 714, 725 in Ch7 (similar for Ch8) 
438   // for the futur official numbering
439   // Type 112233 for 705, 707, 755, 757 in Ch7 (similar for Ch8) 
440   // actual numbering in the code to be changed in jan05
441   // Type 112233 for 901, 902, 911, 912, 914, 915, 924, 925 in Ch9 
442   // (similar for Ch10) for the futur official numbering
443   // Type 112233 for 904, 905, 907, 908, 954, 955, 957, 958 in Ch9 
444   // (similar for Ch10) actual numbering in the code to be changed in jan05
445   Int_t n5[4] = { 0, 2, 2, 2 };
446   slatsegB[1]->SetPcbBoards(n5);
447   slatsegB[1]->Init(0); // 0 detection element id
448   slatsegNB[1]->SetPcbBoards(n5);
449   slatsegNB[1]->Init(0); // 0 detection element id
450     
451   // Type 112230 for 702, 711, 715, 724 in Ch7 (similar for Ch8) 
452   // for the futur official numbering
453   // Type 112230 for 704, 708, 754, 758 in Ch7 (similar for Ch8) 
454   // actual numbering in the code to be changed in jan05
455   Int_t n6[4] = { 0, 2, 2, 1 };
456   slatsegB[2]->SetPcbBoards(n6);
457   slatsegB[2]->Init(0); // 0 detection element id
458   slatsegNB[2]->SetPcbBoards(n6);
459   slatsegNB[2]->Init(0); // 0 detection element id
460     
461   // Type 222330 for 703, 710, 716, 723 in Ch7 (similar for Ch8) 
462   // for the futur official numbering
463   // Type 222330 for 703, 709, 753, 759 in Ch7 (similar for Ch8) 
464   // actual numbering in the code to be changed in jan05
465   Int_t n7[4] = { 0, 0, 3, 2 };
466   slatsegB[3]->SetPcbBoards(n7);
467   slatsegB[3]->Init(0); // 0 detection element id
468   slatsegNB[3]->SetPcbBoards(n7);
469   slatsegNB[3]->Init(0); // 0 detection element id
470     
471   // Type 223300 for 704, 709, 717, 722 in Ch7 (similar for Ch8) 
472   // for the futur official numbering
473   // Type 223300 for 702, 710, 752, 760 in Ch7 (similar for Ch8) 
474   // actual numbering in the code to be changed in jan05
475   Int_t n8[4] = { 0, 0, 2, 2 };
476   slatsegB[4]->SetPcbBoards(n8);
477   slatsegB[4]->Init(0); // 0 detection element id
478   slatsegNB[4]->SetPcbBoards(n8);
479   slatsegNB[4]->Init(0); // 0 detection element id
480     
481   // Type 333000 for 705, 708, 718, 721 in Ch7 (similar for Ch8) 
482   // for the futur official numbering
483   // Type 333000 for 701, 711, 751, 761 in Ch7 (similar for Ch8) 
484   // actual numbering in the code to be changed in jan05
485   // Type 333000 for 906, 907, 919, 920 in Ch9 (similar for Ch10) 
486   // for the futur official numbering
487   // Type 333000 for 900, 912, 950, 962 in Ch9 (similar for Ch10) 
488   // actual numbering in the code to be changed in jan05
489   Int_t n9[4] = { 0, 0, 0, 3 };
490   slatsegB[5]->SetPcbBoards(n9);
491   slatsegB[5]->Init(0); // 0 detection element id
492   slatsegNB[5]->SetPcbBoards(n9);
493   slatsegNB[5]->Init(0); // 0 detection element id
494     
495   // Type 330000 for 706, 707, 719, 720 in Ch7 (similar for Ch8) 
496   // for the futur official numbering
497   // Type 330000 for 700, 712, 750, 762 in Ch7 (similar for Ch8) 
498   // actual numbering in the code to be changed in jan05
499   Int_t n10[4] = { 0, 0, 0, 2 };
500   slatsegB[6]->SetPcbBoards(n10);
501   slatsegB[6]->Init(0); // 0 detection element id
502   slatsegNB[6]->SetPcbBoards(n10);
503   slatsegNB[6]->Init(0); // 0 detection element id
504
505
506   for (Int_t chamber = 6; chamber < 8; chamber++) {
507
508     segmentation[0] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
509     segmentation[1] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
510         
511     // id detection elt for chamber 1
512     Int_t id0 = (chamber+1)*100;
513
514     //--------------------------------------------------------
515     // Configuration for Chamber TC6/7  (Station 4) ----------           
516
517
518     fMUON->SetNsec(chamber,2);
519
520     // cathode 0
521     // type 122330
522     segmentation[0]->Add(id0+ 6, slatsegB[0]);
523     segmentation[0]->Add(id0+56, slatsegB[0]);
524   
525     // type 112233
526     segmentation[0]->Add(id0+ 5, slatsegB[1]);
527     segmentation[0]->Add(id0+ 7, slatsegB[1]);  
528     segmentation[0]->Add(id0+55, slatsegB[1]);  
529     segmentation[0]->Add(id0+57, slatsegB[1]);
530    
531     // type 112230
532     segmentation[0]->Add(id0+ 4, slatsegB[2]);
533     segmentation[0]->Add(id0+ 8, slatsegB[2]);  
534     segmentation[0]->Add(id0+54, slatsegB[2]);  
535     segmentation[0]->Add(id0+58, slatsegB[2]);
536
537     // type 222330 
538     segmentation[0]->Add(id0+ 3, slatsegB[3]);
539     segmentation[0]->Add(id0+ 9, slatsegB[3]);  
540     segmentation[0]->Add(id0+53, slatsegB[3]);
541     segmentation[0]->Add(id0+59, slatsegB[3]);
542
543     // type 223300 
544     segmentation[0]->Add(id0+ 2, slatsegB[4]);
545     segmentation[0]->Add(id0+10, slatsegB[4]);  
546     segmentation[0]->Add(id0+52, slatsegB[4]);
547     segmentation[0]->Add(id0+60, slatsegB[4]);
548
549     // type 333000 
550     segmentation[0]->Add(id0+ 1, slatsegB[5]);
551     segmentation[0]->Add(id0+11, slatsegB[5]);  
552     segmentation[0]->Add(id0+51, slatsegB[5]);
553     segmentation[0]->Add(id0+61, slatsegB[5]);
554
555     // type 330000 
556     segmentation[0]->Add(id0   , slatsegB[6]);
557     segmentation[0]->Add(id0+12, slatsegB[6]);  
558     segmentation[0]->Add(id0+50, slatsegB[6]);
559     segmentation[0]->Add(id0+62, slatsegB[6]);
560     fMUON->SetSegmentationModel(chamber, 1, segmentation[0]);   
561
562     // cathode 1
563     // type 122330
564     segmentation[1]->Add(id0+ 6, slatsegNB[0]);
565     segmentation[1]->Add(id0+56, slatsegNB[0]);
566
567     // type 112233
568     segmentation[1]->Add(id0+ 5, slatsegNB[1]);
569     segmentation[1]->Add(id0+ 7, slatsegNB[1]);  
570     segmentation[1]->Add(id0+55, slatsegNB[1]);  
571     segmentation[1]->Add(id0+57, slatsegNB[1]);
572   
573     // type 112230
574     segmentation[1]->Add(id0+ 4, slatsegNB[2]);
575     segmentation[1]->Add(id0+ 8, slatsegNB[2]);  
576     segmentation[1]->Add(id0+54, slatsegNB[2]);  
577     segmentation[1]->Add(id0+58, slatsegNB[2]);
578
579     // type 222330 
580     segmentation[1]->Add(id0+ 3, slatsegNB[3]);
581     segmentation[1]->Add(id0+ 9, slatsegNB[3]);  
582     segmentation[1]->Add(id0+53, slatsegNB[3]);
583     segmentation[1]->Add(id0+59, slatsegNB[3]);
584
585     // type 223300 
586     segmentation[1]->Add(id0+ 2, slatsegNB[4]);
587     segmentation[1]->Add(id0+10, slatsegNB[4]);  
588     segmentation[1]->Add(id0+52, slatsegNB[4]);
589     segmentation[1]->Add(id0+60, slatsegNB[4]);
590
591     // type 333000 
592     segmentation[1]->Add(id0+ 1, slatsegNB[5]);
593     segmentation[1]->Add(id0+11, slatsegNB[5]);  
594     segmentation[1]->Add(id0+51, slatsegNB[5]);
595     segmentation[1]->Add(id0+61, slatsegNB[5]);
596
597     // type 330000 
598     segmentation[1]->Add(id0   , slatsegNB[6]);
599     segmentation[1]->Add(id0+12, slatsegNB[6]);  
600     segmentation[1]->Add(id0+50, slatsegNB[6]);
601     segmentation[1]->Add(id0+62, slatsegNB[6]);
602     fMUON->SetSegmentationModel(chamber, 2, segmentation[1]);
603
604     fMUON->SetResponseModel(chamber, fResponse0);      
605     fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread
606   }
607 }
608
609 //__________________________________________________________________________
610 void AliMUONFactoryV2::BuildStation5() 
611 {       
612   //--------------------------------------------------------
613   // Configuration for Chamber TC9/10  (Station 5) ---------           
614   //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
615
616   AliMUONGeometrySegmentation* segmentation[2];
617
618   //Slats Segmentations
619   AliMUONSt345SlatSegmentation *slatsegB[6]; // Types of segmentation for St5
620   AliMUONSt345SlatSegmentation *slatsegNB[6]; 
621   // Bending
622
623   Int_t ndiv[4] ={ 4, 4, 2, 1};  // densities zones 
624   for(Int_t i = 0; i < 6; i++) {
625     slatsegB[i] = new AliMUONSt345SlatSegmentation(1);
626     fDESegmentations->Add(slatsegB[i]);  
627     slatsegB[i]->SetPadSize(10.,0.5);
628     slatsegB[i]->SetPadDivision(ndiv);
629     slatsegB[i]->SetId(1);
630     slatsegB[i]->SetDAnod(AliMUONConstants::Pitch());
631     slatsegNB[i] = new AliMUONSt345SlatSegmentation(0);
632     fDESegmentations->Add(slatsegNB[i]);  
633     slatsegNB[i]->SetPadSize(1./1.4,10.); 
634     slatsegNB[i]->SetPadDivision(ndiv);
635     slatsegNB[i]->SetId(1);
636     slatsegNB[i]->SetDAnod(AliMUONConstants::Pitch());
637   }
638
639   // Type 122330 for 900, 913 in Ch9 (similar for Ch10) 
640   // for the futur official numbering
641   // Type 122330 for 906, 956 in Ch9 (similar for Ch10) 
642   // actual numbering in the code to be changed in jan05
643   Int_t n4[4] = { 0, 1, 2, 2 };
644   slatsegB[0]->SetPcbBoards(n4);
645   slatsegB[0]->Init(0); // 0 detection element id
646   slatsegNB[0]->SetPcbBoards(n4);
647   slatsegNB[0]->Init(0); // 0 detection element id
648     
649   // Type 112233 for 901, 902, 911, 912, 914, 915, 924, 925 in Ch9 
650   // (similar for Ch10) for the futur official numbering
651   // Type 112233 for 904, 905, 907, 908, 954, 955, 957, 958 in Ch9 
652   // (similar for Ch10) actual numbering in the code to be changed in jan05
653   Int_t n5[4] = { 0, 2, 2, 2 };
654   slatsegB[1]->SetPcbBoards(n5);
655   slatsegB[1]->Init(0); // 0 detection element id
656   slatsegNB[1]->SetPcbBoards(n5);
657   slatsegNB[1]->Init(0); // 0 detection element id
658
659   // Type 333000 for 906, 907, 919, 920 in Ch9 (similar for Ch10) 
660   // for the futur official numbering
661   // Type 333000 for 900, 912, 950, 962 in Ch9 (similar for Ch10) 
662   // actual numbering in the code to be changed in jan05
663   Int_t n9[4] = { 0, 0, 0, 3 };
664   slatsegB[2]->SetPcbBoards(n9);
665   slatsegB[2]->Init(0); // 0 detection element id
666   slatsegNB[2]->SetPcbBoards(n9);
667   slatsegNB[2]->Init(0); // 0 detection element id
668
669   // Type 222333 for 903, 910, 916, 923 in Ch9 (similar for Ch10) 
670   // for the futur official numbering
671   // Type 222333 for 903, 909, 953, 959 in Ch9 (similar for Ch10) 
672   // actual numbering in the code to be changed in jan05
673   Int_t n11[4] = { 0, 0, 3, 3 };
674   slatsegB[3]->SetPcbBoards(n11);
675   slatsegB[3]->Init(0); // 0 detection element id
676   slatsegNB[3]->SetPcbBoards(n11);
677   slatsegNB[3]->Init(0); // 0 detection element id
678   
679   // Type 223330 for 904, 909, 917, 922 in Ch9 (similar for Ch10) 
680   // for the futur official numbering
681   // Type 223330 for 902, 910, 952, 960 in Ch9 (similar for Ch10) 
682   // actual numbering in the code to be changed in jan05
683   Int_t n12[4] = { 0, 0, 2, 3 };
684   slatsegB[4]->SetPcbBoards(n12);
685   slatsegB[4]->Init(0); // 0 detection element id
686   slatsegNB[4]->SetPcbBoards(n12);
687   slatsegNB[4]->Init(0); // 0 detection element id
688     
689   // Type 333300 for 905, 908, 918, 921 in Ch9 (similar for Ch10) 
690   // for the futur official numbering
691   // Type 333300 for 901, 911, 951, 961 in Ch9 (similar for Ch10) 
692   // actual numbering in the code to be changed in jan05
693   Int_t n13[4] = { 0, 0, 0, 4 };
694   slatsegB[5]->SetPcbBoards(n13);
695   slatsegB[5]->Init(0); // 0 detection element id
696   slatsegNB[5]->SetPcbBoards(n13);
697   slatsegNB[5]->Init(0); // 0 detection element id
698
699   for (Int_t chamber = 8; chamber < 10; chamber++) {
700
701     segmentation[0] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
702     segmentation[1] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
703
704     // id detection elt for chamber 1
705     Int_t id0 = (chamber+1)*100;
706
707     //--------------------------------------------------------
708     // Configuration for Chamber TC8/9  (Station 5) ----------           
709
710
711     fMUON->SetNsec(chamber,2);
712
713     // cathode 0
714     // type 122330
715     segmentation[0]->Add(id0+ 6, slatsegB[0]);
716     segmentation[0]->Add(id0+56, slatsegB[0]);
717   
718     // type 112233
719     segmentation[0]->Add(id0+ 4, slatsegB[1]);
720     segmentation[0]->Add(id0+ 5, slatsegB[1]);
721     segmentation[0]->Add(id0+ 7, slatsegB[1]);  
722     segmentation[0]->Add(id0+ 8, slatsegB[1]);  
723     segmentation[0]->Add(id0+54, slatsegB[1]);  
724     segmentation[0]->Add(id0+55, slatsegB[1]);  
725     segmentation[0]->Add(id0+57, slatsegB[1]);
726     segmentation[0]->Add(id0+58, slatsegB[1]);
727
728     // type 333000 
729     segmentation[0]->Add(id0   , slatsegB[2]);
730     segmentation[0]->Add(id0+12, slatsegB[2]);  
731     segmentation[0]->Add(id0+50, slatsegB[2]);
732     segmentation[0]->Add(id0+62, slatsegB[2]);
733  
734     // type 222333 
735     segmentation[0]->Add(id0+ 3, slatsegB[3]);
736     segmentation[0]->Add(id0+ 9, slatsegB[3]);  
737     segmentation[0]->Add(id0+53, slatsegB[3]);
738     segmentation[0]->Add(id0+59, slatsegB[3]);
739  
740     // type 223330 
741     segmentation[0]->Add(id0+ 2, slatsegB[4]);
742     segmentation[0]->Add(id0+10, slatsegB[4]);  
743     segmentation[0]->Add(id0+52, slatsegB[4]);
744     segmentation[0]->Add(id0+60, slatsegB[4]);
745   
746     // type 333300 
747     segmentation[0]->Add(id0+ 1, slatsegB[5]);
748     segmentation[0]->Add(id0+11, slatsegB[5]);  
749     segmentation[0]->Add(id0+51, slatsegB[5]);
750     segmentation[0]->Add(id0+61, slatsegB[5]);
751     fMUON->SetSegmentationModel(chamber, 1, segmentation[0]);
752
753     // cathode 1
754     // type 122330
755     segmentation[1]->Add(id0+ 6, slatsegNB[0]);
756     segmentation[1]->Add(id0+56, slatsegNB[0]);
757   
758     // type 112233
759     segmentation[1]->Add(id0+ 4, slatsegNB[1]);
760     segmentation[1]->Add(id0+ 5, slatsegNB[1]);
761     segmentation[1]->Add(id0+ 7, slatsegNB[1]);  
762     segmentation[1]->Add(id0+ 8, slatsegNB[1]);  
763     segmentation[1]->Add(id0+54, slatsegNB[1]);  
764     segmentation[1]->Add(id0+55, slatsegNB[1]);  
765     segmentation[1]->Add(id0+57, slatsegNB[1]);
766     segmentation[1]->Add(id0+58, slatsegNB[1]);
767
768     // type 333000 
769     segmentation[1]->Add(id0   , slatsegNB[2]);
770     segmentation[1]->Add(id0+12, slatsegNB[2]);  
771     segmentation[1]->Add(id0+50, slatsegNB[2]);
772     segmentation[1]->Add(id0+62, slatsegNB[2]);
773  
774     // type 222333 
775     segmentation[1]->Add(id0+ 3, slatsegNB[3]);
776     segmentation[1]->Add(id0+ 9, slatsegNB[3]);  
777     segmentation[1]->Add(id0+53, slatsegNB[3]);
778     segmentation[1]->Add(id0+59, slatsegNB[3]);
779  
780     // type 223330 
781     segmentation[1]->Add(id0+ 2, slatsegNB[4]);
782     segmentation[1]->Add(id0+10, slatsegNB[4]);  
783     segmentation[1]->Add(id0+52, slatsegNB[4]);
784     segmentation[1]->Add(id0+60, slatsegNB[4]);
785   
786     // type 333300 
787     segmentation[1]->Add(id0+ 1, slatsegNB[5]);
788     segmentation[1]->Add(id0+11, slatsegNB[5]);  
789     segmentation[1]->Add(id0+51, slatsegNB[5]);
790     segmentation[1]->Add(id0+61, slatsegNB[5]);
791     fMUON->SetSegmentationModel(chamber, 2, segmentation[1]);
792         
793     fMUON->SetResponseModel(chamber, fResponse0);           
794     fMUON->Chamber(chamber).SetChargeCorrel(0.11); // 11% charge spread
795   }
796 }
797
798 //__________________________________________________________________________
799 void AliMUONFactoryV2::BuildStation6() 
800 {       
801  // Create Trigger geometry segmentation for given chamber and cathod
802
803  
804     AliMUONGeometrySegmentation *chamberSeg[2];
805 // Cluster-size off
806         AliMUONResponseTrigger* responseTrigger0 =  new AliMUONResponseTrigger;
807 // Cluster-size on  
808 //  AliMUONResponseTriggerV1* responseTrigger0 =  new AliMUONResponseTriggerV1;
809
810     for (Int_t chamber = 10; chamber < 14; chamber++) {
811
812       //Trigger Segmentation
813       AliMUONTriggerSegmentation *trigSegX[9]; 
814       AliMUONTriggerSegmentation *trigSegY[9]; 
815       for(Int_t i=0; i<9; i++) {
816         trigSegX[i] = new AliMUONTriggerSegmentation(1);
817         trigSegY[i] = new AliMUONTriggerSegmentation(0);
818         fDESegmentations->Add(trigSegX[i]);  
819         fDESegmentations->Add(trigSegY[i]);  
820         trigSegX[i]->SetLineNumber(9-i);    
821         trigSegY[i]->SetLineNumber(9-i);    
822       }
823
824       AliMUONChamber *iChamber, *iChamber1;
825       iChamber1 = &fMUON->Chamber(10);
826       iChamber  = &fMUON->Chamber(chamber);
827       Float_t zpos1= - iChamber1->Z();  
828       Float_t zpos = - iChamber->Z();        
829       Float_t zRatio = zpos / zpos1;
830
831       // init
832       Float_t stripWidth[3]={0.,0.,0.};     // 1.0625 2.125 4.25
833       Float_t stripLength[4]={0.,0.,0.,0.}; // 17. 34. 51. 68.
834       for (Int_t i=0; i<3; i++) 
835         stripWidth[i]=AliMUONTriggerConstants::StripWidth(i)*zRatio;
836       for (Int_t i=0; i<4; i++) 
837         stripLength[i]=AliMUONTriggerConstants::StripLength(i)*zRatio;
838       Int_t nStrip[7]={0,0,0,0,0,0,0};    
839       Float_t stripYsize[7]={0.,0.,0.,0.,0.,0.,0.};
840       Float_t stripXsize[7]={0.,0.,0.,0.,0.,0.,0.};
841
842       // chamber 8 0 cathode 0
843       for (Int_t i=0; i<7; i++) nStrip[i]=16;
844       for (Int_t i=0; i<7; i++) stripYsize[i]=stripWidth[2];
845       for (Int_t i=0; i<6; i++) stripXsize[i]=stripLength[1];
846       stripXsize[6]=stripLength[2];
847       trigSegX[8]->Init(0,nStrip,stripYsize,stripXsize,0.); 
848       trigSegX[0]->Init(0,nStrip,stripYsize,stripXsize,0.); 
849
850       // chamber 8 7 1 0 cathode 1
851       for (Int_t i=0; i<6; i++) nStrip[i]=8;
852       nStrip[6]=16;
853       for (Int_t i=0; i<7; i++) stripYsize[i]=stripLength[3];  
854       for (Int_t i=0; i<7; i++) stripXsize[i]=stripWidth[2];
855       trigSegY[8]->Init(0,nStrip,stripYsize,stripXsize,0.);  
856       trigSegY[7]->Init(0,nStrip,stripYsize,stripXsize,0.);
857       trigSegY[1]->Init(0,nStrip,stripYsize,stripXsize,0.);
858       trigSegY[0]->Init(0,nStrip,stripYsize,stripXsize,0.);
859  
860       // chamber 7 6 2 1 cathode 0
861       for (Int_t i=0; i<6; i++) nStrip[i]=32;
862       nStrip[6]=16;  
863       for (Int_t i=0; i<6; i++) stripYsize[i]=stripWidth[1];
864       stripYsize[6]=stripWidth[2];
865       for (Int_t i=0; i<6; i++) stripXsize[i]=stripLength[1];
866       stripXsize[6]=stripLength[2];
867       trigSegX[7]->Init(0,nStrip,stripYsize,stripXsize,0.);  
868       trigSegX[6]->Init(0,nStrip,stripYsize,stripXsize,0.);
869       trigSegX[2]->Init(0,nStrip,stripYsize,stripXsize,0.);  
870       trigSegX[1]->Init(0,nStrip,stripYsize,stripXsize,0.);
871
872       // chamber 6 2 cathode 1
873       for (Int_t i=0; i<5; i++) nStrip[i]=16;
874       for (Int_t i=5; i<6; i++) nStrip[i]=8;
875       nStrip[6]=16;
876       for (Int_t i=0; i<7; i++) stripYsize[i]=stripLength[3];
877       for (Int_t i=0; i<5; i++) stripXsize[i]=stripWidth[1];
878       for (Int_t i=5; i<7; i++) stripXsize[i]=stripWidth[2];
879       trigSegY[6]->Init(0,nStrip,stripYsize,stripXsize,0.);  
880       trigSegY[2]->Init(0,nStrip,stripYsize,stripXsize,0.);  
881
882       // chamber 5 3 cathode 0
883       nStrip[0]=48;
884       for (Int_t i=1; i<3; i++) nStrip[i]=64;
885       for (Int_t i=3; i<6; i++) nStrip[i]=32;
886       nStrip[6]=16;  
887       for (Int_t i=0; i<3; i++) stripYsize[i]=stripWidth[0];
888       for (Int_t i=3; i<6; i++) stripYsize[i]=stripWidth[1];
889       stripYsize[6]=stripWidth[2];
890       for (Int_t i=0; i<6; i++) stripXsize[i]=stripLength[1];
891       stripXsize[6]=stripLength[2];
892       trigSegX[5]->Init(0,nStrip,stripYsize,stripXsize,stripLength[0]);  
893       trigSegX[3]->Init(0,nStrip,stripYsize,stripXsize,0.);
894
895       // chamber 5 3 cathode 1
896       for (Int_t i=0; i<5; i++) nStrip[i]=16;
897       for (Int_t i=5; i<6; i++) nStrip[5]=8;  
898       nStrip[6]=16;  
899       stripYsize[0]=stripLength[2];
900       for (Int_t i=1; i<7; i++) stripYsize[i]=stripLength[3];
901       for (Int_t i=0; i<5; i++) stripXsize[i]=stripWidth[1];
902       for (Int_t i=5; i<7; i++) stripXsize[i]=stripWidth[2];
903       trigSegY[5]->Init(0,nStrip,stripYsize,stripXsize,stripLength[0]);  
904       trigSegY[3]->Init(0,nStrip,stripYsize,stripXsize,0.);
905
906       // chamber 4 cathode 0
907       nStrip[0]=0;
908       for (Int_t i=1; i<3; i++) nStrip[i]=64;  
909       for (Int_t i=3; i<6; i++) nStrip[i]=32;  
910       nStrip[6]=16;
911       stripYsize[0]=0.;
912       for (Int_t i=1; i<3; i++) stripYsize[i]=stripWidth[0];
913       for (Int_t i=3; i<6; i++) stripYsize[i]=stripWidth[1];
914       stripYsize[6]=stripWidth[2];
915       stripXsize[0]=0;  
916       stripXsize[1]=stripLength[0];  
917       for (Int_t i=2; i<6; i++) stripXsize[i]=stripLength[1];
918       stripXsize[6]=stripLength[2];
919       trigSegX[4]->Init(0,nStrip,stripYsize,stripXsize,0.);  
920
921       // chamber 4 cathode 1
922       nStrip[0]=0;  
923       nStrip[1]=8;  
924       for (Int_t i=2; i<5; i++) nStrip[i]=16;
925       for (Int_t i=5; i<6; i++) nStrip[i]=8;
926       nStrip[6]=16;
927       stripYsize[0]=0.;  
928       for (Int_t i=1; i<7; i++) stripYsize[i]=stripLength[3];
929       stripXsize[0]=0.;
930       for (Int_t i=1; i<5; i++) stripXsize[i]=stripWidth[1];
931       for (Int_t i=5; i<7; i++) stripXsize[i]=stripWidth[2];
932       trigSegY[4]->Init(0,nStrip,stripYsize,stripXsize,0.);
933
934       chamberSeg[0] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
935       chamberSeg[1] = new AliMUONGeometrySegmentation(fMUON->Chamber(chamber).GetGeometry());
936
937       fMUON->SetNsec(chamber,2);
938       Int_t icount=chamber-10;  // chamber counter (0 1 2 3)
939       Int_t id0=(10+icount+1)*100;
940
941
942       //  printf("in CreateTriggerSegmentation here 0 id0=%i \n",id0);  
943
944       for (Int_t i = 0; i < 9; i++) {       
945
946         // cathode 0
947         chamberSeg[0]->Add(id0+i,     trigSegX[i]);
948         chamberSeg[0]->Add(id0+50+i,  trigSegX[i]);
949         fMUON->SetSegmentationModel(chamber, 1, chamberSeg[0]);
950
951         // cathode 1
952         chamberSeg[1]->Add(id0+i,     trigSegY[i]);
953         chamberSeg[1]->Add(id0+50+i,  trigSegY[i]);
954         fMUON->SetSegmentationModel(chamber, 2, chamberSeg[1]);
955
956       }
957
958       fMUON->SetResponseModel(chamber, responseTrigger0);      
959       fMUON->Chamber(chamber).SetChargeCorrel(0); // same charge on cathodes
960
961   
962       //  printf("in CreateTriggerSegmentation here 1\n");  
963
964       if (!id0) {
965         AliWarning(Form("Segmentation for chamber %d is not yet defined",chamber));
966         return ;      
967       }
968     }
969 }       
970 //__________________________________________________________________________
971 void AliMUONFactoryV2::Build(AliMUON* where, const char* what) 
972 {
973   //
974   // Construct MUON from chambers, segmentation and responses
975   //
976
977   fMUON = where;
978   char tmp[20];
979   strcpy(tmp, what);
980
981   if (strcmp(tmp, "default")==0) {
982     // Set default parameters
983     fMUON->SetIshunt(0);
984     fMUON->SetMaxStepGas(0.1);
985     fMUON->SetMaxStepAlu(0.1);
986
987     // Build stations
988     BuildCommon();
989     if (IsGeometryDefined(0))  BuildStation1();
990     if (IsGeometryDefined(2))  BuildStation2();
991     if (IsGeometryDefined(4))  BuildStation3();
992     if (IsGeometryDefined(6))  BuildStation4();
993     if (IsGeometryDefined(8))  BuildStation5();
994     if (IsGeometryDefined(10)) BuildStation6();
995   } 
996   else
997     AliDebug(0,"Non default version of MUON selected. You have to construct yourself the MUON elements !!");
998 }
999
1000 //__________________________________________________________________________
1001 void AliMUONFactoryV2::BuildStation(AliMUON* where, Int_t stationNumber) 
1002 {
1003   //
1004   // Construct MUON from chambers, segmentation and responses
1005   //
1006   // Version 0
1007   //
1008   // First define the number of planes that are segmented (1 or 2) by a call
1009   // to SetNsec.
1010   // Then chose for each chamber (chamber plane) the segmentation
1011   // and response model.
1012   // They should be equal for the two chambers of each station. In a future
1013   // version this will be enforced.
1014   //
1015
1016   fMUON = where;
1017   if (!fResponse0) BuildCommon(); 
1018     
1019   switch (stationNumber) {    
1020   case 1:  BuildStation1(); break;
1021   case 2:  BuildStation2(); break;
1022   case 3:  BuildStation3(); break;
1023   case 4:  BuildStation4(); break;
1024   case 5:  BuildStation5(); break;
1025   case 6:  BuildStation6(); break;
1026     
1027   default: AliFatal("Wrong station number");
1028   }  
1029 }