Clean-up of include files
[u/mrichter/AliRoot.git] / MUON / AliMUONSt2GeometryBuilderV2.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 // $Id$
17
18 //-----------------------------------------------------------------------------
19 // Class AliMUONSt2GeometryBuilderV2
20 // ---------------------------------
21 // MUON Station2 coarse geometry construction class.
22 // Author: SANJOY PAL ,Prof. SUKALYAN CHATTOPADHAYAY  [SINP, KOLKATA]
23 //         &  Dr.SHAKEEL AHMAD (AMU), INDIA
24 //-----------------------------------------------------------------------------
25
26
27 #include <TVirtualMC.h>
28 #include <TGeoMatrix.h>
29 #include <Riostream.h>
30 #include <TArrayI.h>
31
32 #include "AliRun.h"
33 #include "AliLog.h"
34
35 #include "AliMUONSt2GeometryBuilderV2.h"
36 #include "AliMUON.h"
37 #include "AliMUONConstants.h"
38 #include "AliMUONGeometryModule.h"
39 #include "AliMUONGeometryEnvelopeStore.h"
40 #include "AliMUONConstants.h"
41
42 #define PI 3.14159
43
44 /// \cond CLASSIMP
45 ClassImp(AliMUONSt2GeometryBuilderV2)
46 /// \endcond
47
48 //______________________________________________________________________________
49 AliMUONSt2GeometryBuilderV2::AliMUONSt2GeometryBuilderV2(AliMUON* muon)
50  : AliMUONVGeometryBuilder(2, 2),
51    fMUON(muon)
52 {
53 /// Standard constructor
54
55 }
56
57 //______________________________________________________________________________
58 AliMUONSt2GeometryBuilderV2::AliMUONSt2GeometryBuilderV2()
59  : AliMUONVGeometryBuilder(),
60    fMUON(0)
61 {
62 /// Default constructor
63 }
64
65 //______________________________________________________________________________
66 AliMUONSt2GeometryBuilderV2::~AliMUONSt2GeometryBuilderV2() 
67 {
68 /// Destructor
69 }
70
71 //
72 // public methods
73 //
74
75 //______________________________________________________________________________
76 void AliMUONSt2GeometryBuilderV2::CreateGeometry()
77 {
78 /// Geometry construction
79
80 //
81 //********************************************************************
82 //                            Station 2                             **
83 //********************************************************************
84      // indices 1 and 2 for first and second chambers in the station
85      // iChamber (first chamber) kept for other quanties than Z,
86      // assumed to be the same in both chambers
87
88      // Get tracking medias Ids
89      Int_t *idtmed = fMUON->GetIdtmed()->GetArray()-1099;
90      Int_t idAir  = idtmed[1100]; // medium 1
91      Int_t idGas  = idtmed[1108]; // medium Ar-CO2 gas (80%+20%)
92      Int_t idPCB  = idtmed[1122]; // medium FR4
93      Int_t idCU   = idtmed[1110]; // medium copper
94      Int_t idRoha = idtmed[1113]; // medium roha cell
95      Int_t idPGF30= idtmed[1123]; // medium for Frame Eq.to Bakelite
96      //Int_t idScru = idtmed[1128]; // screw material - Stainless Steel(18%Cr,9%Ni,Fe)
97
98
99
100 /*########################################################################################
101     Create volume for one Quadrant
102 ##########################################################################################*/
103      Float_t tpar1[5];
104      tpar1[0] = 20.55;
105      tpar1[1] = 123.5;
106      tpar1[2] = 6.8/2;
107      tpar1[3] = -12.0;
108      tpar1[4] = 102.0;
109
110
111      gMC->Gsvolu("SQM3","TUBS", idAir, tpar1, 5);
112      gMC->Gsvolu("SQM4","TUBS", idAir, tpar1, 5);
113
114
115 //==================================================================================
116 //                                Plane      
117 //==================================================================================
118
119 //Thickness of variour parts
120   Float_t zCu   = 0.005;      // eff. cu in cathode pcb
121   Float_t zCbb  = 0.04;       //cathode pcb
122   Float_t zRoha = 2.5;        // Rhocell
123   Float_t zMeb = 0.04;      //Mech. exit board //0.08
124   Float_t zEeb = 0.04;     //Effective electronic readout board //0.02
125
126   //starting Z-positions of various parts--- in Plane-1
127   Float_t zposCu    = 0.25;   // 2.5 mm => gap between anode & chatode plane
128   Float_t zposCbb   = zposCu + zCu;
129   Float_t zposRoha  = zposCbb + zCbb;
130   Float_t zposMeb   = zposRoha + zRoha;
131   Float_t zposEeb   = zposMeb + zMeb;
132
133   Float_t zposCuBar   = zposCu +  zCu/2.; //for segment 0 & 6
134   Float_t zposCbbBar  = zposCbb + zCbb/2.;  
135   Float_t zposRohaBar = zposRoha + zRoha/2.;
136   Float_t zposMebBar  = zposMeb + zMeb/2.;
137   Float_t zposEebBar  = zposEeb + zEeb/2.;
138
139
140  //Cathode PCB + Copper sheet + Rohacell + mech exit board + eff. electronic exit board
141
142  //Segment-0 ~~~Horizantal box  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
143   Float_t bparH[3];
144   bparH[0] = 95.5/2.; // extension beyond 0 deg in x direction // 94.5
145   bparH[1] = 1.2/2.;  // extension beyond 0 deg in y direction (3.7[total extn] - 2.5[frame dim])
146   bparH[2] = zCu/2.;     //Thickness of Copper sheet
147   gMC->Gsvolu("SCU0L", "BOX", idCU, bparH, 3);
148   
149   bparH[2] = zCbb/2.; // thickness of cathode sheet in z direction
150   gMC->Gsvolu("SCB0L", "BOX", idPCB, bparH, 3);
151   
152   bparH[2] = zRoha/2.;     //Thickness of Roha cell
153   gMC->Gsvolu("SRH0L", "BOX", idRoha, bparH, 3);
154
155   bparH[2] = zMeb/2;       //Thickness of mechanical exit board
156   gMC->Gsvolu("SMB0L", "BOX", idPCB, bparH, 3);
157   
158   bparH[2] = zEeb/2;          //Thickness of effective electronic  exit board
159   gMC->Gsvolu("SEB0L", "BOX", idCU, bparH, 3);
160   
161   //Segment-1 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
162   Float_t pgpar[10]; // polygon
163   pgpar[0] = 0.;  // initial angle
164   pgpar[1] = 90.; // increment in angle starting from initial angle 
165   pgpar[2] = 5;   // number of side
166   pgpar[3] = 2.; // number of plane
167   pgpar[4] = -zCu/2.; // z-position of the first plane
168   pgpar[5] = 23.1; // innner radius first plane
169   pgpar[6] = 117.6;  // outer radious first plane
170   pgpar[7] = zCu/2.; // z-position of the second plane
171   pgpar[8] = pgpar[5];  // innner radius of second plane
172   pgpar[9] = pgpar[6];  // outer radious of second plane
173   gMC->Gsvolu("SCU1L", "PGON", idCU, pgpar, 10);
174
175   pgpar[4] = -zCbb/2.; // z-position of the first plane
176   pgpar[7] = zCbb/2.;  // Thickness of copper-sheet
177   gMC->Gsvolu("SCB1L", "PGON", idPCB, pgpar, 10);  
178   
179   pgpar[4] = -zRoha/2.;  // Thickness of Roha cell
180   pgpar[7] = zRoha/2.;  // Thickness of Roha cell
181   gMC->Gsvolu("SRH1L", "PGON", idRoha, pgpar, 10);
182        
183   pgpar[4] = -zMeb/2.;  // Thickness of mechanical exit board
184   pgpar[7] = zMeb/2.;  // Thickness of mechanical exit board
185   gMC->Gsvolu("SMB1L", "PGON", idPCB, pgpar, 10);
186   
187   pgpar[4] = -zEeb/2.;  // Thickness of electronic readout board
188   pgpar[7] = zEeb/2.;  // Thickness of electronic readout board
189   gMC->Gsvolu("SEB1L", "PGON", idCU, pgpar, 10);
190        
191
192 //Segment-2 - vertical box (simalar to horizontal bar as in Segment 0)~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
193        
194   Float_t bparV[3];
195   
196   bparV[0] = 1.0/2.;
197   bparV[1] = 95.5/2.; // 94.5
198   bparV[2] = zCu/2.;
199   gMC->Gsvolu("SCU2L", "BOX", idCU, bparV, 3);
200   
201   bparV[2] = zCbb/2.;
202   gMC->Gsvolu("SCB2L", "BOX", idPCB, bparV, 3);
203   
204   
205   bparV[2] = zRoha/2.;
206   gMC->Gsvolu("SRH2L", "BOX", idRoha, bparV, 3);
207   
208   bparV[2] = zMeb/2;
209   gMC->Gsvolu("SMB2L", "BOX", idPCB, bparV, 3);
210   
211   bparV[2] = zEeb/2;
212   gMC->Gsvolu("SEB2L", "BOX", idCU, bparV, 3);
213        
214   //....(Setting posion of Segment 0,1,2)..................................................................
215   
216   Float_t xposHorBox =  bparH[0] + 23.1; // 23.1 = 20.6(inner radius of qrd) + 2.5 (width of frame)
217   Float_t yposHorBox = -bparH[1];  
218   
219   Float_t xposVerBox = -bparV[0];
220   Float_t yposVerBox =  bparV[1] + 23.1;
221   
222   
223   //Positioning the PCB
224   
225   
226   // chamber 3
227   gMC->Gspos("SCB0L",1, "SQM3",xposHorBox,yposHorBox,zposCbbBar,0, "ONLY");
228   gMC->Gspos("SCB1L",1, "SQM3", 0.0,0.0,zposCbbBar,0, "ONLY");
229   gMC->Gspos("SCB2L",1, "SQM3",xposVerBox,yposVerBox,zposCbbBar,0, "ONLY");
230   
231   gMC->Gspos("SCB0L",2, "SQM3",xposHorBox,yposHorBox,-zposCbbBar,0, "ONLY");
232   gMC->Gspos("SCB1L",2, "SQM3", 0.0,0.0,-zposCbbBar,0, "ONLY");
233   gMC->Gspos("SCB2L",2, "SQM3",xposVerBox,yposVerBox,-zposCbbBar,0, "ONLY");
234
235   // chamber 4
236   gMC->Gspos("SCB0L",3, "SQM4",xposHorBox,yposHorBox,zposCbbBar,0, "ONLY");
237   gMC->Gspos("SCB1L",3, "SQM4", 0.0,0.0,zposCbbBar,0, "ONLY");
238   gMC->Gspos("SCB2L",3, "SQM4",xposVerBox,yposVerBox,zposCbbBar,0, "ONLY");
239   
240   gMC->Gspos("SCB0L",4, "SQM4",xposHorBox,yposHorBox,-zposCbbBar,0, "ONLY");
241   gMC->Gspos("SCB1L",4, "SQM4", 0.0,0.0,-zposCbbBar,0, "ONLY");
242   gMC->Gspos("SCB2L",4, "SQM4",xposVerBox,yposVerBox,-zposCbbBar,0, "ONLY");
243   
244
245   //Positioning Copper sheet
246   
247   // chamber 3
248
249   gMC->Gspos("SCU0L",1, "SQM3",xposHorBox,yposHorBox,zposCuBar,0, "ONLY");
250   gMC->Gspos("SCU1L",1, "SQM3", 0.0,0.0,zposCuBar,0, "ONLY");
251   gMC->Gspos("SCU2L",1, "SQM3",xposVerBox,yposVerBox,zposCuBar,0, "ONLY");
252   
253   gMC->Gspos("SCU0L",2, "SQM3",xposHorBox,yposHorBox,-zposCuBar,0, "ONLY");
254   gMC->Gspos("SCU1L",2, "SQM3", 0.0,0.0,-zposCuBar,0, "ONLY");
255   gMC->Gspos("SCU2L",2, "SQM3",xposVerBox,yposVerBox,-zposCuBar,0, "ONLY");
256
257   // chamber 4
258   gMC->Gspos("SCU0L",3, "SQM4",xposHorBox,yposHorBox,zposCuBar,0, "ONLY");
259   gMC->Gspos("SCU1L",3, "SQM4", 0.0,0.0,zposCuBar,0, "ONLY");
260   gMC->Gspos("SCU2L",3, "SQM4",xposVerBox,yposVerBox,zposCuBar,0, "ONLY");
261   
262   gMC->Gspos("SCU0L",4, "SQM4",xposHorBox,yposHorBox,-zposCuBar,0, "ONLY");
263   gMC->Gspos("SCU1L",4, "SQM4", 0.0,0.0,-zposCuBar,0, "ONLY");
264   gMC->Gspos("SCU2L",4, "SQM4",xposVerBox,yposVerBox,-zposCuBar,0, "ONLY");
265      
266   //Positioning Roha cell 
267   
268   // chamber 3
269   
270   gMC->Gspos("SRH0L",1, "SQM3",xposHorBox,yposHorBox,zposRohaBar,0, "ONLY");
271   gMC->Gspos("SRH1L",1, "SQM3", 0.0,0.0,zposRohaBar,0, "ONLY");
272   gMC->Gspos("SRH2L",1, "SQM3",xposVerBox,yposVerBox,zposRohaBar,0, "ONLY");
273   
274   gMC->Gspos("SRH0L",2, "SQM3",xposHorBox,yposHorBox,-zposRohaBar,0, "ONLY");
275   gMC->Gspos("SRH1L",2, "SQM3", 0.0,0.0,-zposRohaBar,0, "ONLY");
276   gMC->Gspos("SRH2L",2, "SQM3",xposVerBox,yposVerBox,-zposRohaBar,0, "ONLY");
277   
278   // chamber 4
279   
280   gMC->Gspos("SRH0L",3, "SQM4",xposHorBox,yposHorBox,zposRohaBar,0, "ONLY");
281   gMC->Gspos("SRH1L",3, "SQM4", 0.0,0.0,zposRohaBar,0, "ONLY");
282   gMC->Gspos("SRH2L",3, "SQM4",xposVerBox,yposVerBox,zposRohaBar,0, "ONLY");
283   
284   gMC->Gspos("SRH0L",4, "SQM4",xposHorBox,yposHorBox,-zposRohaBar,0, "ONLY");
285   gMC->Gspos("SRH1L",4, "SQM4", 0.0,0.0,-zposRohaBar,0, "ONLY");
286   gMC->Gspos("SRH2L",4, "SQM4",xposVerBox,yposVerBox,-zposRohaBar,0, "ONLY");
287
288   //Positioning of Mech. exit board
289   
290   // chamber 3
291   
292   gMC->Gspos("SMB0L",1, "SQM3",xposHorBox,yposHorBox,zposMebBar,0, "ONLY");
293   gMC->Gspos("SMB1L",1, "SQM3", 0.0,0.0,zposMebBar,0, "ONLY");
294   gMC->Gspos("SMB2L",1, "SQM3",xposVerBox,yposVerBox,zposMebBar,0, "ONLY");
295   
296   gMC->Gspos("SMB0L",2, "SQM3",xposHorBox,yposHorBox,-zposMebBar,0, "ONLY");
297   gMC->Gspos("SMB1L",2, "SQM3", 0.0,0.0,-zposMebBar,0, "ONLY");
298   gMC->Gspos("SMB2L",2, "SQM3",xposVerBox,yposVerBox,-zposMebBar,0, "ONLY");
299   // chamber 4
300   
301   gMC->Gspos("SMB0L",3, "SQM4",xposHorBox,yposHorBox,zposMebBar,0, "ONLY");
302   gMC->Gspos("SMB1L",3, "SQM4", 0.0,0.0,zposMebBar,0, "ONLY");
303   gMC->Gspos("SMB2L",3, "SQM4",xposVerBox,yposVerBox,zposMebBar,0, "ONLY");
304   
305   gMC->Gspos("SMB0L",4, "SQM4",xposHorBox,yposHorBox,-zposMebBar,0, "ONLY");
306   gMC->Gspos("SMB1L",4, "SQM4", 0.0,0.0,-zposMebBar,0, "ONLY");
307   gMC->Gspos("SMB2L",4, "SQM4",xposVerBox,yposVerBox,-zposMebBar,0, "ONLY");
308   
309   //Positioning of Electronic exit board
310   
311   
312   // chamber 3
313
314   gMC->Gspos("SEB0L",1, "SQM3",xposHorBox,yposHorBox,zposEebBar,0, "ONLY");
315   gMC->Gspos("SEB1L",1, "SQM3", 0.0,0.0,zposEebBar,0, "ONLY");
316   gMC->Gspos("SEB2L",1, "SQM3",xposVerBox,yposVerBox,zposEebBar,0, "ONLY");
317   
318   gMC->Gspos("SEB0L",2, "SQM3",xposHorBox,yposHorBox,-zposEebBar,0, "ONLY");
319   gMC->Gspos("SEB1L",2, "SQM3", 0.0,0.0,-zposEebBar,0, "ONLY");
320   gMC->Gspos("SEB2L",2, "SQM3",xposVerBox,yposVerBox,-zposEebBar,0, "ONLY");
321   // chamber 4
322
323   gMC->Gspos("SEB0L",3, "SQM4",xposHorBox,yposHorBox,zposEebBar,0, "ONLY");
324   gMC->Gspos("SEB1L",3, "SQM4", 0.0,0.0,zposEebBar,0, "ONLY");
325   gMC->Gspos("SEB2L",3, "SQM4",xposVerBox,yposVerBox,zposEebBar,0, "ONLY");
326   
327   gMC->Gspos("SEB0L",4, "SQM4",xposHorBox,yposHorBox,-zposEebBar,0, "ONLY");
328   gMC->Gspos("SEB1L",4, "SQM4", 0.0,0.0,-zposEebBar,0, "ONLY");
329   gMC->Gspos("SEB2L",4, "SQM4",xposVerBox,yposVerBox,-zposEebBar,0, "ONLY");
330   
331   
332   //----------------------------------------------------------------------
333   //                         Frames
334   //----------------------------------------------------------------------
335   //Frame-1
336
337   Float_t frame1[3] ;                
338   frame1[0] = 101.0/2.;             //100.6 = 94.5 + 2.5 + 3.6 
339   frame1[1] = 2.5/2.;               
340   frame1[2] = 5.0/2.;
341   
342   gMC->Gsvolu("SFRM1", "BOX", idPGF30, frame1, 3); //Frame - 1 // fill with pkk GF30
343   
344   Float_t arib1[3];
345   arib1[0] = frame1[0];
346   arib1[1] = 0.9/2.;
347   arib1[2] =(frame1[2]-0.95)/2.0;
348   
349   gMC->Gsvolu("SFRA1", "BOX", idAir, arib1, 3); // fill with air
350   
351   Float_t xposarib1 = 0;
352   Float_t yposarib1 = -frame1[1] + arib1[1];
353   Float_t zposarib1 = frame1[2] - arib1[2];
354   
355   gMC->Gspos("SFRA1",1, "SFRM1", xposarib1, yposarib1, zposarib1,0, "ONLY");  //replace pkk GF30 with air(b)
356   gMC->Gspos("SFRA1",2, "SFRM1", xposarib1, yposarib1, -zposarib1,0, "ONLY"); //replace pkk GF30 with air(nb)
357   
358   Float_t rrib1[3];
359   rrib1[0] = frame1[0];
360   rrib1[1] = 0.6/2.;
361   rrib1[2] =(frame1[2]-0.95)/2.0;
362   
363   gMC->Gsvolu("SFRR1", "BOX", idRoha, rrib1, 3); // fill with rohacell
364   
365   Float_t xposrrib1 = 0.0;
366   Float_t yposrrib1 = frame1[1] - rrib1[1];
367   Float_t zposrrib1 = frame1[2] - rrib1[2];
368   
369   gMC->Gspos("SFRR1",1, "SFRM1", xposrrib1, yposrrib1, zposrrib1,0, "ONLY");//replace pkk GF30 with rohacell
370   gMC->Gspos("SFRR1",2, "SFRM1", xposrrib1, yposrrib1, -zposrrib1,0, "ONLY");//replace pkk GF30 with rohacell
371   
372   
373   Float_t xposFr1 = frame1[0] + 20.6;
374   Float_t yposFr1 = -3.7 + frame1[1] ;
375   Float_t zposFr1 = 0.0;
376   
377   gMC->Gspos("SFRM1",1, "SQM3", xposFr1, yposFr1, zposFr1,0, "ONLY");// frame -1
378   gMC->Gspos("SFRM1",2, "SQM4", xposFr1, yposFr1, zposFr1,0, "ONLY");// frame -1
379   
380   
381   //......................................................................................
382   //Frame-2
383
384   Float_t frame2[3]; 
385   frame2[0] = 4.0/2.;
386   frame2[1] = 1.2/2.;
387   frame2[2] = 5.0/2;
388   
389   gMC->Gsvolu("SFRM2", "BOX", idPGF30, frame2, 3); //Frame - 2
390   
391   Float_t rrib2[3];
392   rrib2[0] = frame2[0]-1.0/2.0;
393   rrib2[1] = frame2[1];
394   rrib2[2] =(frame2[2]-0.95)/2.0;
395
396   gMC->Gsvolu("SFRR2", "BOX", idRoha, rrib2, 3);
397       
398   Float_t xposrrib2 = -1.0/2.0;
399   Float_t yposrrib2 = 0.0;
400   Float_t zposrrib2 = frame2[2] - rrib2[2];
401   
402   gMC->Gspos("SFRR2",1, "SFRM2", xposrrib2, yposrrib2, zposrrib2,0, "ONLY");//replace pkk GF30 with rohacell
403   gMC->Gspos("SFRR2",2, "SFRM2", xposrrib2, yposrrib2, -zposrrib2,0, "ONLY");//replace pkk GF30 with roha
404   
405
406   
407   Float_t xposFr2 = frame2[0] + 117.6;
408   Float_t yposFr2 = -frame2[1];
409   Float_t zposFr2 = 0.0;
410   
411   gMC->Gspos("SFRM2",1, "SQM3", xposFr2, yposFr2, zposFr2,0, "MANY");//global positing of frame in SQM3
412   gMC->Gspos("SFRM2",2, "SQM4", xposFr2, yposFr2, zposFr2,0, "MANY");//global positing of frame in SQM4
413   
414   
415   //......................................................................................
416   
417   //Frame-3
418   
419   Float_t pgparFr3[10];
420   pgparFr3[0] = 0.;
421   pgparFr3[1] = 90.;
422   pgparFr3[2] = 5;
423   pgparFr3[3] = 2.;
424   pgparFr3[4] = 0.;
425   pgparFr3[5] = 117.6;
426   pgparFr3[6] = 121.6;
427   pgparFr3[7] = pgparFr3[4] + 5.0;
428   pgparFr3[8] = pgparFr3[5];
429   pgparFr3[9] = pgparFr3[6];
430   
431   gMC->Gsvolu("SFRM3", "PGON", idPGF30, pgparFr3, 10);
432   
433   Float_t pgparRrib3[10];
434   pgparRrib3[0] = 0.;
435   pgparRrib3[1] = 90.;
436   pgparRrib3[2] = 5;
437   pgparRrib3[3] = 2.;
438   pgparRrib3[4] = 0.;
439   pgparRrib3[5] = 117.6;
440   pgparRrib3[6] = 120.6;
441   pgparRrib3[7] = pgparRrib3[4] +1.55 ;
442   pgparRrib3[8] = pgparRrib3[5];
443   pgparRrib3[9] = pgparRrib3[6];
444   
445   gMC->Gsvolu("SFRR3", "PGON", idRoha, pgparRrib3, 10);
446   
447   Float_t xposrrib3 = 0.0;
448   Float_t yposrrib3 = 0.0;
449   Float_t zposrrib3 = 0.0;
450
451   gMC->Gspos("SFRR3",1, "SFRM3", xposrrib3, yposrrib3, zposrrib3,0, "ONLY");
452   
453   zposrrib3 = 3.45;
454   
455   gMC->Gspos("SFRR3",2, "SFRM3", xposrrib3, yposrrib3, zposrrib3,0, "ONLY");
456   
457       
458   
459   Float_t xposFr3 = 0.0;
460   Float_t yposFr3 = 0.0;
461   Float_t zposFr3 = -frame1[2];
462       
463   gMC->Gspos("SFRM3",1, "SQM3", xposFr3, yposFr3, zposFr3,0, "ONLY");// frame -1
464   gMC->Gspos("SFRM3",2, "SQM4", xposFr3, yposFr3, zposFr3,0, "ONLY");// frame -1
465   
466   
467   //......................................................................................
468   //Frame-4
469   
470   Float_t frame4[3]; 
471   frame4[0] = 1.0/2.;
472   frame4[1] = 4.0/2.;
473   frame4[2] = frame1[2];
474   
475   gMC->Gsvolu("SFRM4", "BOX", idPGF30, frame4, 3); 
476   
477   Float_t rrib4[3];
478   rrib4[0] = frame4[0];
479   rrib4[1] = frame4[1]-1.0/2;
480   rrib4[2] =(frame4[2]-0.95)/2.0;
481
482   gMC->Gsvolu("SFRR4", "BOX", idRoha, rrib4, 3);
483   
484   Float_t xposrrib4 = 0.0;
485   Float_t yposrrib4 = -1.0/2;
486   Float_t zposrrib4 = frame4[2] - rrib4[2];
487   
488   gMC->Gspos("SFRR4",1, "SFRM4", xposrrib4, yposrrib4, zposrrib4,0, "ONLY");
489   gMC->Gspos("SFRR4",2, "SFRM4", xposrrib4, yposrrib4, -zposrrib4,0, "ONLY");
490   
491
492   
493   Float_t xposFr4 = -frame4[0];
494   Float_t yposFr4 = -frame4[1] + 117.6;
495   Float_t zposFr4 = 0.0;
496   
497   gMC->Gspos("SFRM4",1, "SQM3", xposFr4, yposFr4, zposFr4,0, "MANY");
498   gMC->Gspos("SFRM4",2, "SQM4", xposFr4, yposFr4, zposFr4,0, "MANY");
499   
500   
501   //......................................................................................
502   //Frame-5
503   Float_t frame5[3] ;               
504   frame5[0] = 2.7/2.;            
505   frame5[1] = 101.0/2.;              
506   frame5[2] = 5.0/2.;
507   
508   gMC->Gsvolu("SFRM5", "BOX", idPGF30, frame5, 3); //Frame - 1
509       
510   Float_t arib5[3];
511   arib5[0] = 0.9/2.0;
512   arib5[1] = frame5[1];
513   arib5[2] = (frame5[2]-0.95)/2.0;
514   
515   gMC->Gsvolu("SFRA5", "BOX", idAir, arib5, 3);
516   
517   Float_t xposarib5 = -frame5[0] + arib5[0];
518   Float_t yposarib5 = 0.0;
519   Float_t zposarib5 = frame5[2] - arib5[2];
520       
521   gMC->Gspos("SFRA5",1, "SFRM5", xposarib5, yposarib5, zposarib5,0, "ONLY");
522   gMC->Gspos("SFRA5",2, "SFRM5", xposarib5, yposarib5, -zposarib5,0, "ONLY");
523   
524   Float_t rrib5[3];
525   rrib5[0] = 0.8/2.0;
526   rrib5[1] = frame5[1];
527   rrib5[2] = (frame5[2]-0.95)/2.0;
528   
529   gMC->Gsvolu("SFRR5", "BOX", idRoha, rrib5, 3);
530   
531   Float_t xposrrib5 = frame5[0] - rrib5[0];
532   Float_t yposrrib5 = 0.0;
533   Float_t zposrrib5 = frame5[2] - rrib5[2];
534   
535   gMC->Gspos("SFRR5",1, "SFRM5", xposrrib5, yposrrib5, zposrrib5,0, "ONLY");
536   gMC->Gspos("SFRR5",2, "SFRM5", xposrrib5, yposrrib5, -zposrrib5,0, "ONLY");
537   
538   Float_t xposFr5 = -3.7 + frame5[0];
539   Float_t yposFr5 = frame5[1] + 20.6;
540   Float_t zposFr5 = 0.0;
541   
542   gMC->Gspos("SFRM5",1, "SQM3", xposFr5, yposFr5, zposFr5,0, "ONLY");
543   gMC->Gspos("SFRM5",2, "SQM4", xposFr5, yposFr5, zposFr5,0, "ONLY");
544   
545   //......................................................................................
546   //Frame -6 
547   
548   Float_t frame6[3]; 
549   frame6[0] = 1.0/2.;
550   frame6[1] = 2.5/2.;
551   frame6[2] = frame1[2];
552   
553   gMC->Gsvolu("SFRM6", "BOX", idPGF30, frame6, 3); 
554   
555   Float_t rrib6[3];
556   rrib6[0] = frame6[0];
557   rrib6[1] = 1.5/2.;
558   rrib6[2] =(frame2[2]-0.95)/2.0;
559
560   gMC->Gsvolu("SFRR6", "BOX", idRoha, rrib6, 3);
561   
562   Float_t xposrrib6 = 0.0;
563   Float_t yposrrib6 = 1.0/2.0;
564   Float_t zposrrib6 = frame6[2] - rrib6[2];
565
566   gMC->Gspos("SFRR6",1, "SFRM6", xposrrib6, yposrrib6, zposrrib6,0, "ONLY");
567   gMC->Gspos("SFRR6",2, "SFRM6", xposrrib6, yposrrib6, -zposrrib6,0, "ONLY");
568   
569   
570   
571   Float_t xposFr6 = -frame6[0];
572   Float_t yposFr6 = frame6[1] + 20.6;
573   Float_t zposFr6 = 0.0;
574
575   gMC->Gspos("SFRM6",1, "SQM3", xposFr6, yposFr6, zposFr6,0, "ONLY");
576   gMC->Gspos("SFRM6",2, "SQM4", xposFr6, yposFr6, zposFr6,0, "ONLY");
577   
578   
579 //......................................................................................
580 //Frame - 7 inner pgon
581   
582   Float_t pgparFr7[10];
583   pgparFr7[0] = 0.;
584   pgparFr7[1] = 90.;
585   pgparFr7[2] = 5;
586   pgparFr7[3] = 2.;
587   pgparFr7[4] = 0.;
588   pgparFr7[5] = 20.6;
589   pgparFr7[6] = 23.1;
590   pgparFr7[7] = pgparFr7[4] + 5.0;
591   pgparFr7[8] = pgparFr7[5];
592   pgparFr7[9] = pgparFr7[6];
593   
594   gMC->Gsvolu("SFRM7", "PGON", idPGF30, pgparFr7, 10);
595   
596   Float_t pgparRrib7[10];
597   pgparRrib7[0] = 0.;
598   pgparRrib7[1] = 90.;
599   pgparRrib7[2] = 5;
600   pgparRrib7[3] = 2.;
601   pgparRrib7[4] = 0.;
602   pgparRrib7[5] = 21.6;
603   pgparRrib7[6] = 23.1;
604   pgparRrib7[7] = pgparRrib7[4] + 1.55;
605   pgparRrib7[8] = pgparRrib7[5];
606   pgparRrib7[9] = pgparRrib7[6];
607   
608   gMC->Gsvolu("SFRR7", "PGON", idRoha, pgparRrib7, 10);
609   
610   Float_t xposrrib7 = 0.0;
611   Float_t yposrrib7 = 0.0;
612   Float_t zposrrib7 = 0.0;
613   
614   gMC->Gspos("SFRR7",1, "SFRM7", xposrrib7, yposrrib7, zposrrib7,0, "ONLY");
615   
616   zposrrib7 = 3.45;
617   
618   gMC->Gspos("SFRR7",2, "SFRM7", xposrrib7, yposrrib7, zposrrib7,0, "ONLY");
619   
620   Float_t xposFr7 = 0.0;
621   Float_t yposFr7 = 0.0;
622   Float_t zposFr7 = -frame1[2];
623   
624   gMC->Gspos("SFRM7",1, "SQM3", xposFr7, yposFr7, zposFr7,0, "ONLY");
625   gMC->Gspos("SFRM7",2, "SQM4", xposFr7, yposFr7, zposFr7,0, "ONLY");
626   
627   
628   //......................................................................................
629   //Frame - 8 
630
631   
632   Float_t frame8[3] ;
633   frame8[0] = 2.5/2.0;
634   frame8[1] = 1.2/2.0;
635   frame8[2] = frame1[2];
636   
637   gMC->Gsvolu("SFRM8", "BOX", idPGF30, frame8, 3); //Frame - 2
638       
639   Float_t rrib8[3];
640   rrib8[0] = frame8[0]-1.0/2;
641   rrib8[1] = frame8[1];
642   rrib8[2] =(frame8[2]-0.95)/2.0;
643   
644   gMC->Gsvolu("SFRR8", "BOX", idRoha, rrib8, 3);
645   
646   Float_t xposrrib8 = -1.0/2;
647   Float_t yposrrib8 = 0.0;
648   Float_t zposrrib8 = frame8[2] - rrib8[2];
649   
650   gMC->Gspos("SFRR8",1, "SFRM8", xposrrib8, yposrrib8, zposrrib8,0, "ONLY");
651   gMC->Gspos("SFRR8",2, "SFRM8", xposrrib8, yposrrib8, -zposrrib8,0, "ONLY");
652   
653   
654   
655   Float_t xposFr8 = frame8[0] + 20.6;
656   Float_t yposFr8 = -frame8[1];
657   Float_t zposFr8 = 0.0;
658   
659   gMC->Gspos("SFRM8",1, "SQM3", xposFr8, yposFr8, zposFr8,0, "ONLY");
660   gMC->Gspos("SFRM8",2, "SQM4", xposFr8, yposFr8, zposFr8,0, "ONLY");
661   
662   
663   
664
665   //^^^^^^^^^^^^^^^^^^^^^^^^^ Sensitive volumes ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
666   
667   Float_t zsenv = 0.5; // distance between two cathode plane
668   
669   //Segment-0 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
670   bparH[0] = 94.5/2.;
671   bparH[1] = 1.2/2.;
672   bparH[2] = zsenv/2.;
673   gMC->Gsvolu("SC3G0", "BOX", idGas, bparH, 3);
674   gMC->Gsvolu("SC4G0", "BOX", idGas, bparH, 3);
675   
676  //Segment-1 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
677   pgpar[0] = 0.;
678   pgpar[1] = 90.;
679   pgpar[2] = 5;
680   pgpar[3] = 2.;
681   pgpar[4] = -zsenv/2.;
682   pgpar[5] = 23.1;
683   pgpar[6] = 117.6;
684   pgpar[7] = zsenv/2.;
685   pgpar[8] = pgpar[5];
686   pgpar[9] = pgpar[6];
687   gMC->Gsvolu("SC3G1", "PGON", idGas, pgpar, 10);
688   gMC->Gsvolu("SC4G1", "PGON", idGas, pgpar, 10);
689   
690
691   //Segment-2 - vertical box ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
692   
693   bparV[0] = 1.0/2.;
694   bparV[1] = 95.5/2.;
695   bparV[2] = zsenv/2.;
696   gMC->Gsvolu("SC3G2", "BOX", idGas, bparV, 3);
697   gMC->Gsvolu("SC4G2", "BOX", idGas, bparV, 3);
698   
699   //...........................................................................................
700   
701   xposHorBox =  bparH[0] + 23.1;
702   yposHorBox = -bparH[1];
703   
704   xposVerBox = -bparV[0];
705   yposVerBox = bparV[1] + 23.1;
706   
707   gMC->Gspos("SC3G0",1, "SQM3", xposHorBox,yposHorBox,0.,0, "ONLY");
708   gMC->Gspos("SC3G1",1, "SQM3", 0.,0.,0.,0, "ONLY");
709   gMC->Gspos("SC3G2",1, "SQM3", xposVerBox,yposVerBox,0.,0, "ONLY");
710      
711   
712   gMC->Gspos("SC4G0",1, "SQM4", xposHorBox,yposHorBox,0.,0, "ONLY");
713   gMC->Gspos("SC4G1",1, "SQM4", 0.,0.,0.,0, "ONLY");
714   gMC->Gspos("SC4G2",1, "SQM4", xposVerBox,yposVerBox,0.,0, "ONLY");
715      
716   
717   //^^^^^^^^^^^^^^^^^^^^^^^^^ Sensitive volumes ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
718
719
720
721 //##################################################################################################
722 //   Positioning Quadrant  in chamber#3 and chamber#4
723 //##################################################################################################
724 /******Transformations for  Quadrant**********************************************
725       ||        I  => Quadrant I:   no rotation
726       ||
727   II. || I.    II  => Quadrant II:  Reflaction of Quadrant I in XZ plane
728       ||           => TGeoRotation("Qrot2",90.,180.,90.,90.,180.,0.);
729 =============
730       ||       III => Quadrant III: 180 degree rotation of Quadrant I in XY plane
731  III. || IV.       => TGeoRotation("Qrot3",90.,180.,90.,90.,180.,0.);
732       ||        IV => Quadrant IV:-180 degree rotation of Quadrant II in XY plane
733                    => TGeoRotation("Qrot4",90.,0.,90.,-90.,180.,0.);
734 **********************************************************************************************/
735
736  Int_t detElemId1 =  1;  // quadrant I
737  Int_t detElemId2 =  0;  // quadrant II
738  Int_t detElemId3 =  3;  // quadrant III
739  Int_t detElemId4 =  2;  // quadrant IV
740
741  //Float_t halfChamber = zCbb + zCu + zRoha + zMeb + zEeb + zsenv/2;
742 //   cout<<  "\n half_chamber \t" << half_chamber << endl;
743
744  Float_t halfChamber = tpar1[2];// tpar1[2] = 6.8/2;
745
746 // ------------------------------St2 Chamber3------------------------------------------------
747
748  //    GetEnvelopes(2)->AddEnvelope("S3M0", 300, true,TGeoTranslation(0.,0.,0.));
749     GetEnvelopes(2)->AddEnvelope("SQM3", 300+detElemId1, 1, TGeoTranslation( 0., 0., - halfChamber));
750     GetEnvelopes(2)->AddEnvelope("SQM3", 300+detElemId2, 2, TGeoTranslation( 0., 0., + halfChamber),
751                                  TGeoRotation("Qrot3",90.,180.,90.,90.,180.,0.));
752     GetEnvelopes(2)->AddEnvelope("SQM3", 300+detElemId3, 3, TGeoTranslation( 0., 0., - halfChamber),
753                                   TGeoRotation("Qrot3",90.,180.,90.,270.,0.,0.));
754     GetEnvelopes(2)->AddEnvelope("SQM3", 300+detElemId4, 4, TGeoTranslation( 0., 0., + halfChamber),
755                                   TGeoRotation("Qrot4",90.,0.,90.,-90.,180.,0.));
756
757 //--------------------------------St2 Chamber4-------------------------------------------------
758
759     GetEnvelopes(3)->AddEnvelope("SQM4", 400+detElemId1, 1, TGeoTranslation( 0., 0., - halfChamber));
760     GetEnvelopes(3)->AddEnvelope("SQM4", 400+detElemId2, 2, TGeoTranslation( 0., 0., + halfChamber),
761                                  TGeoRotation("Qrot2",90.,180.,90.,90.,180.,0.));
762     GetEnvelopes(3)->AddEnvelope("SQM4", 400+detElemId3, 3, TGeoTranslation( 0., 0., - halfChamber),
763                                   TGeoRotation("Qrot3",90.,180.,90.,270.,0.,0.));
764     GetEnvelopes(3)->AddEnvelope("SQM4", 400+detElemId4, 4, TGeoTranslation( 0., 0., + halfChamber),
765                                   TGeoRotation("Qrot4",90.,0.,90.,-90.,180.,0.));
766
767 //**********************************************************************************************
768
769 }
770 //______________________________________________________________________________
771 void AliMUONSt2GeometryBuilderV2::SetVolumes()
772 {
773 /// Defines the volumes for the station2 chambers.
774
775   if (gAlice->GetModule("SHIL")) {
776     SetMotherVolume(2, "YOUT1");
777     SetMotherVolume(3, "YOUT1");
778   }  
779
780   // Define chamber volumes as virtual
781   SetVolume(2, "SC03", true);
782   SetVolume(3, "SC04", true);
783 }
784
785 //______________________________________________________________________________
786 void AliMUONSt2GeometryBuilderV2::SetTransformations()
787 {
788 /// Defines the transformations for the station2 chambers.
789
790   Double_t zpos1 = - AliMUONConstants::DefaultChamberZ(2); 
791   SetTranslation(2, TGeoTranslation(0., 0., zpos1));
792
793   Double_t zpos2 = - AliMUONConstants::DefaultChamberZ(3); 
794   SetTranslation(3, TGeoTranslation(0., 0., zpos2));
795 }
796
797 //______________________________________________________________________________
798 void AliMUONSt2GeometryBuilderV2::SetSensitiveVolumes()
799 {
800 /// Defines the sensitive volumes for station2 chambers.
801
802   GetGeometry(2)->SetSensitiveVolume("SC3G0");
803   GetGeometry(2)->SetSensitiveVolume("SC3G1");
804   GetGeometry(2)->SetSensitiveVolume("SC3G2");
805
806   GetGeometry(3)->SetSensitiveVolume("SC4G0");
807   GetGeometry(3)->SetSensitiveVolume("SC4G1");
808   GetGeometry(3)->SetSensitiveVolume("SC4G2");
809 }
810