b9f1715d71b7958070eef8cb3ffa80ca44f2a53b
[u/mrichter/AliRoot.git] / STRUCT / AliBODY.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 //                                                                           //
20 //  Alice external volume                                                    //
21 //  This class contains the description of the Alice external volume         //
22 //                                                                           //
23 //Begin_Html
24 /*
25 <img src="picts/AliBODYClass.gif">
26 </pre>
27 <br clear=left>
28 <font size=+2 color=red>
29 <p>The responsible person for this module is
30 <a href="mailto:andreas.morsch@cern.ch">Andreas Morsch</a>.
31 </font>
32 <pre>
33 */
34 //End_Html
35 //                                                                           //
36 //                                                                           //
37 ///////////////////////////////////////////////////////////////////////////////
38
39 #include <TGeoGlobalMagField.h>
40 #include <TVirtualMC.h>
41
42 #include "AliBODY.h"
43 #include "AliMagF.h"
44 #include "AliRun.h"
45
46 ClassImp(AliBODY)
47  
48 //_____________________________________________________________________________
49 AliBODY::AliBODY()
50 {
51   //
52   // Default constructor
53   //
54 }
55  
56 //_____________________________________________________________________________
57 AliBODY::AliBODY(const char *name, const char *title)
58        : AliModule(name,title)
59 {
60   //
61   // Standard constructor of the Alice external volume
62   //
63   //PH  SetMarkerColor(7);
64   //PH  SetMarkerStyle(2);
65   //PH  SetMarkerSize(0.4);
66 }
67  
68 //_____________________________________________________________________________
69 void AliBODY::CreateGeometry()
70 {
71   //
72   // Create the geometry of the Alice external body
73   //
74   //Begin_Html
75   /*
76     <img src="picts/AliBODYTree.gif">
77   */
78   //End_Html
79   //
80   // If the ZDC is present we have an asymmetric box
81   // made by a four sides polygone
82   //  
83   //Begin_Html
84   /*
85     <img src="picts/AliBODYLarge.gif">
86   */
87   //End_Html
88   //
89   // If the ZDC is not present make just a BOX
90   //
91   //Begin_Html
92   /*
93     <img src="picts/AliBODYSmall.gif">
94   */
95   //End_Html
96
97   Float_t dALIC[10];
98   Int_t *idtmed = fIdtmed->GetArray()+1;
99   //
100   if(gAlice->GetModule("ZDC")) {
101     //
102     // If the ZDC is present we have an asymmetric box
103     // made by a four sides polygone
104     //
105     dALIC[0]=45;
106     dALIC[1]=360;
107     dALIC[2]=4;
108     dALIC[3]=2;
109
110     dALIC[4]=-15000;
111     dALIC[5]=0;
112     dALIC[6]=2000;
113
114     dALIC[7]=15000;
115     dALIC[8]=0;
116     dALIC[9]=2000;
117     gMC->Gsvolu("ALIC","PGON",idtmed[1],dALIC,10);
118   } else if ( gAlice->GetModule("ACORDE")) {
119     //
120     // If the Cosmic Ray Trigger  is present we need a large box
121     // 
122     //
123     dALIC[0]=13000.;
124     dALIC[1]=5000.;
125     dALIC[2]=13000.;
126     gMC->Gsvolu("ALIC","BOX ",idtmed[1],dALIC,3);
127       
128   } else {
129     //
130     // If the ZDC and ACORDE are not present make just a BOX
131     //
132     dALIC[0]=2000;
133     dALIC[1]=2000;
134     dALIC[2]=3000;
135     gMC->Gsvolu("ALIC","BOX ",idtmed[1],dALIC,3);
136   }
137 }
138  
139 //_____________________________________________________________________________
140 void AliBODY::CreateMaterials()
141 {
142 // Create materials and media
143   Int_t isxfld = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Integ();
144   Float_t sxmgmx = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Max();
145   
146   // AIR
147
148   Float_t aAir[4]={12.0107,14.0067,15.9994,39.948};
149   Float_t zAir[4]={6.,7.,8.,18.};
150   Float_t wAir[4]={0.000124,0.755267,0.231781,0.012827};
151   Float_t dAir = 1.20479E-3;
152   Float_t dAir1 = 1.20479E-10;
153   //
154   AliMixture(1,"Vacuum  $",aAir,zAir,dAir1,4,wAir);
155   AliMixture(2,"Air     $",aAir,zAir,dAir,4,wAir);
156   AliMaterial(3,"Be      $", 9.01,4 ,1.848   ,35.30,36.70);
157   //
158   AliMedium(1,"Vacuum  $",1,0,isxfld,sxmgmx,10,1,0.1,0.1,10);
159   AliMedium(2,"Air     $",2,0,isxfld,sxmgmx,10,-1,-0.1,0.1 ,-10);
160   AliMedium(3,"Be pipe $",3,0,isxfld,sxmgmx,10,0.1,0.1,0.01,0.01);
161 }
162  
163 //_____________________________________________________________________________
164 void AliBODY::DrawModule() const
165 {
166   //
167   // Draw a view of the Alice outside box
168   //
169   // Set everything unseen
170   gMC->Gsatt("*", "seen", -1);
171   // 
172   // Set ALIC mother visible
173   gMC->Gsatt("ALIC","SEEN",1);
174   //
175   // Set the volumes visible
176   //
177   gMC->Gdopt("hide","off");
178   if(gAlice->GetModule("ZDC")) {
179     //
180     // ZDC is present
181     //
182     gMC->DefaultRange();
183     gMC->Gdraw("alic", 40, 30, 0, 15, 10, .0014, .0014);
184     gMC->Gdhead(1111, "Aice Main body with Zero Degree Calorimeter");
185   } else {
186     //
187     // ZDC is not present
188     //
189     gMC->Gdraw("alic", 40, 30, 0, 10, 9, .0027, .0027);
190     gMC->Gdhead(1111, "Aice Main body");
191   }
192   gMC->Gdman(18, 4, "MAN");
193 }
194  
195