Removing AliMCProcess and AliMC
[u/mrichter/AliRoot.git] / STRUCT / AliMAG.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 $Log$
18 Revision 1.13  2002/11/06 15:20:32  hristov
19 Corrected logics (I.Hrivnacova)
20
21 Revision 1.12  2002/11/01 10:59:49  hristov
22 Correcting geometry (corresponds to version 1.9)
23
24 Revision 1.11  2002/10/17 16:33:07  hristov
25 Corrected check of the MC type
26
27 Revision 1.10  2002/10/17 16:14:05  hristov
28 MAG geometry with resolved MANY (I.Hrivnacova)
29
30 Revision 1.9  2001/05/16 14:57:22  alibrary
31 New files for folders and Stack
32
33 Revision 1.8  2000/10/02 21:28:15  fca
34 Removal of useless dependecies via forward declarations
35
36 Revision 1.7  2000/06/11 12:35:56  morsch
37 Coding rule violations corrected
38
39 Revision 1.6  1999/09/29 09:24:30  fca
40 Introduction of the Copyright and cvs Log
41
42 */
43
44 ///////////////////////////////////////////////////////////////////////////////
45 //                                                                           //
46 //  L3 Magnet                                                                //
47 //                                                                           //
48 //Begin_Html
49 /*
50 <img src="picts/AliMAGClass.gif">
51 </pre>
52 <br clear=left>
53 <font size=+2 color=red>
54 <p>The responsible person for this module is
55 <a href="mailto:andreas.morsch@cern.ch">Andreas Morsch</a>.
56 </font>
57 <pre>
58 */
59 //End_Html
60 //                                                                           //
61 ///////////////////////////////////////////////////////////////////////////////
62  
63 #include "AliMAG.h"
64 #include "AliRun.h"
65 #include "AliMagF.h"
66  
67 ClassImp(AliMAG)
68  
69 //_____________________________________________________________________________
70 AliMAG::AliMAG()
71 {
72   //
73   // Default constructor for L3 magnet
74   //
75 }
76  
77 //_____________________________________________________________________________
78 AliMAG::AliMAG(const char *name, const char *title)
79   : AliModule(name,title)
80 {
81   //
82   // Standard constructor for L3 magnet
83   //
84   //Begin_Html
85   /*
86     <img src="picts/aliMAG.gif">
87   */
88   //End_Html
89   
90   SetMarkerColor(7);
91   SetMarkerStyle(2);
92   SetMarkerSize(0.4);
93 }
94
95 //_____________________________________________________________________________
96 void AliMAG::CreateGeometry()
97 {
98   //
99   // Create geometry for L3 magnet
100   //
101   //Begin_Html
102   /*
103     <img src="picts/mag.gif">
104   */
105   //End_Html
106   //Begin_Html
107   /*
108     <img src="picts/tree_mag.gif">
109   */
110   //End_Html
111   
112   Int_t *idtmed = fIdtmed->GetArray()-299;
113   
114   Float_t dpar[13];
115   Int_t idrotm[399];
116   Float_t par[10];
117   
118   // ANGLE POLAIRE MAXIMUM 
119   
120   //  Define Mother 
121   
122   par[0] = 22.5;
123   par[1] = 360.;
124   par[2] = 8.;
125   par[3] = 2.;
126   par[4] = -600.;
127   par[5] = 580.;
128   par[6] = 790.;
129   par[7] = 600.;
130   par[8] = 580.;
131   par[9] = 790.;
132   gMC->Gsvolu("L3MO", "PGON", idtmed[334], par, 10);
133   gMC->Gspos("L3MO", 1, "ALIC", 0., -30., 0., 0, "ONLY");
134   
135   // Define coils 
136   
137   par[5] = 585.;
138   par[6] = 690.;
139   par[8] = 585.;
140   par[9] = 690.;
141   gMC->Gsvolu("L3CO", "PGON", idtmed[328], par, 10);
142   gMC->Gspos("L3CO", 1, "L3MO", 0., 0., 0., 0, "ONLY");
143   
144   par[5] = 580.;
145   par[6] = 585.;
146   par[8] = 580.;
147   par[9] = 585.;
148   gMC->Gsvolu("L3C1", "PGON", idtmed[308], par, 10);
149   gMC->Gspos("L3C1", 1, "L3MO", 0., 0., 0., 0, "ONLY");
150   
151   // Define yoke 
152   
153   par[5] = 690.;
154   par[6] = 790.;
155   par[8] = 690.;
156   par[9] = 790.;
157   gMC->Gsvolu("L3YO", "PGON", idtmed[329], par, 10);
158   gMC->Gspos("L3YO", 1, "L3MO", 0., 0., 0., 0, "ONLY");
159   
160   // Define the return yoke of L3 (DOOR) 
161   
162   par[4] = 600.;
163   par[5] = 0.;
164   par[6] = 790.;
165   par[7] = 700.;
166   par[8] = par[5];
167   par[9] = par[6];
168   gMC->Gsvolu("L3DO", "PGON", idtmed[334], par, 10);
169   gMC->Gsvolu("L3DX", "PGON", idtmed[334], par, 10);
170   
171   par[4] = 610.;
172   par[5] = 0.;
173   par[6] = 790.;
174   par[7] = 700.;
175   par[8] = par[5];
176   par[9] = par[6];
177   gMC->Gsvolu("L3FR", "PGON", idtmed[329], par, 10);
178   gMC->Gsvolu("L3FX", "PGON", idtmed[329], par, 10);
179   
180   // INNER LAYER 
181   
182   par[4] = 600.;
183   par[7] = 610.;
184   gMC->Gsvolu("L3IR", "PGON", idtmed[309], par, 10);
185   gMC->Gsvolu("L3IX", "PGON", idtmed[309], par, 10);
186   
187   //     DOOR OPENING 
188   
189   dpar[0] = 22.5;
190   dpar[1] = 360;
191   dpar[2] = 8.;
192   dpar[3] = 3.;
193   dpar[4] = 610.;
194   dpar[5] = 0.;
195   dpar[6] = 163.5;
196   dpar[7] = 670.;
197   dpar[8] = dpar[5];
198   dpar[9] = dpar[6];
199   dpar[10] = 700.;
200   dpar[11] = dpar[5];
201   dpar[12] = dpar[6] + 50.;
202   gMC->Gsvolu("L3O1", "PGON", idtmed[314], dpar, 13);
203   gMC->Gsvolu("L3O3", "PGON", idtmed[314], dpar, 13);
204   par[4] = 600.;
205   par[5] = 0.;
206   par[6] = 163.5;
207   par[7] = 610.;
208   par[8] = 0.;
209   par[9] = 163.5;
210   gMC->Gsvolu("L3O2", "PGON", idtmed[314], par, 10);
211   gMC->Gsvolu("L3O4", "PGON", idtmed[314], par, 10);
212   
213   //     THE DOOR OPENING HAS TO BE PLACED WITH 'MANY' SINCE THE REGION 
214   //     WILL CONTAIN A MUON CHAMBER, BEAM PIPE AND BEAM SHIELD 
215   //     PLACED WITH 'ONLY'. 
216   
217   AliMatrix(idrotm[300], 90., 0., 90., 90., 180., 0.);
218
219   gMC->Gspos("L3DO", 1, "ALIC", 0., -30., 0., 0, "MANY");
220   gMC->Gspos("L3FR", 1, "L3DO", 0., 0., 0., 0, "MANY");
221   gMC->Gspos("L3IR", 1, "L3DO", 0., 0., 0., 0, "MANY");
222   
223   gMC->Gspos("L3DX", 1, "ALIC", 0., -30., 0., idrotm[300], "MANY");
224   gMC->Gspos("L3FX", 1, "L3DX", 0., 0., 0., 0, "MANY");
225   gMC->Gspos("L3IX", 1, "L3DX", 0., 0., 0., 0, "MANY");
226
227   if(!strcmp(gMC->GetName(),"TGeant3")) {
228
229     gMC->Gspos("L3O1", 1, "L3FR", 0., 30., 0., 0, "MANY");
230     gMC->Gspos("L3O2", 1, "L3IR", 0., 30., 0., 0, "MANY");
231     gMC->Gspos("L3O3", 1, "L3FX", 0., 30., 0., 0, "MANY");
232     gMC->Gspos("L3O4", 1, "L3IX", 0., 30., 0., 0, "MANY");
233
234   }
235   else {
236     gMC->Gspos("L3O1", 1, "ALIC", 0., 0., 0., 0, "MANY");
237     gMC->Gspos("L3O2", 1, "ALIC", 0., 0., 0., 0, "MANY");
238     gMC->Gsbool("L3O1", "L3DO");
239     gMC->Gsbool("L3O2", "L3DO");
240   
241     gMC->Gspos("L3O3", 1, "ALIC", 0., 0., 0., idrotm[300], "MANY");
242     gMC->Gspos("L3O4", 1, "ALIC", 0., 0., 0., idrotm[300], "MANY");
243     gMC->Gsbool("L3O3", "L3DX");
244     gMC->Gsbool("L3O4", "L3DX");
245   }      
246 }
247
248 //_____________________________________________________________________________
249 void AliMAG::CreateMaterials()
250 {
251   //
252   // Create materials for L3 magnet
253   //
254   
255   Int_t   isxfld = gAlice->Field()->Integ();
256   Float_t sxmgmx = gAlice->Field()->Max();
257   
258   Float_t epsil, stmin, deemax, tmaxfd, stemax;
259   // --- Define the various materials for GEANT --- 
260   
261   //     Aluminum 
262   AliMaterial(9, "Al$", 26.98, 13., 2.7, 8.9, 37.2);
263   AliMaterial(29, "Al$", 26.98, 13., 2.7, 8.9, 37.2);
264   
265   //     Iron 
266   AliMaterial(10, "Fe$", 55.85, 26., 7.87, 1.76, 17.1);
267   AliMaterial(30, "Fe$", 55.85, 26., 7.87, 1.76, 17.1);
268   
269   //     Air 
270   AliMaterial(15, "Air$", 14.61, 7.3, .001205, 30420., 67500);
271   AliMaterial(35, "Air$", 14.61, 7.3, .001205, 30420., 67500);
272   
273   // **************** 
274   //     Defines tracking media parameters. 
275   //     Les valeurs sont commentees pour laisser le defaut 
276   //     a GEANT (version 3-21, page CONS200), f.m. 
277   epsil  = .001;  // Tracking precision, 
278   stemax = -1.;   // Maximum displacement for multiple scat 
279   tmaxfd = -20.;  // Maximum angle due to field deflection 
280   deemax = -.3;   // Maximum fractional energy loss, DLS 
281   stmin  = -.8;
282   // *************** 
283   
284   //    IRON 
285   
286   AliMedium(10, "FE_C0             ", 10, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
287   AliMedium(30, "FE_C1             ", 30, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
288   
289   //     ALUMINUM 
290
291   AliMedium(9, "ALU_C0            ",  9, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
292   AliMedium(29, "ALU_C1            ", 29, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
293   
294   //     AIR 
295   
296   AliMedium(15, "AIR_C0            ", 15, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
297   AliMedium(35, "AIR_C1            ", 35, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
298 }
299
300 //_____________________________________________________________________________
301 void AliMAG::DrawModule()
302 {
303   //
304   // Draw a shaded view of the L3 magnet
305   //
306 }
307
308 //_____________________________________________________________________________
309 void AliMAG::Init()
310 {
311   //
312   // Initialise L3 magnet after it has been built
313   Int_t i;
314   //
315   if(fDebug) {
316     printf("\n%s: ",ClassName());
317     for(i=0;i<35;i++) printf("*");
318     printf(" MAG_INIT ");
319     for(i=0;i<35;i++) printf("*");
320     printf("\n%s: ",ClassName());
321     //
322     // Here the MAG initialisation code (if any!)
323     for(i=0;i<80;i++) printf("*");
324     printf("\n");
325   }
326 }
327