1 ///////////////////////////////////////////////////////////////////////////////
3 // Forward Multiplicity Detector version 0 //
7 <img src="picts/AliFMDv0Class.gif">
10 <font size=+2 color=red>
11 <p>The responsible person for this module is
12 <a href="mailto:Valeri.Kondratiev@cern.ch">Valeri Kondratiev</a>.
18 ///////////////////////////////////////////////////////////////////////////////
26 //_____________________________________________________________________________
27 AliFMDv0::AliFMDv0() : AliFMD()
30 // Default constructor for FMD version 0
34 //_____________________________________________________________________________
35 AliFMDv0::AliFMDv0(const char *name, const char *title)
39 // Standard constructor for FMD version 0
43 //___________________________________________
44 void AliFMDv0::CreateGeometry()
47 // Creation of the geometry of the FMD version 0
51 <img src="picts/AliFMDv0Tree.gif">
56 <img src="picts/AliFMDv0.gif">
64 Int_t *idtmed = fIdtmed->GetArray()-899;
66 // ********************************************************
67 // DEFINE RIGHT DISK#3 OF FMD
68 // ********************************************************
79 gMC->Gsvolu("R1R3", "TUBE", idtmed[900], par, 3);
80 gMC->Gspos("R1R3", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
85 gMC->Gsvolu("R2R3", "TUBE", idtmed[901], par, 3);
86 gMC->Gspos("R2R3", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
91 gMC->Gsvolu("R3R3", "TUBE", idtmed[902], par, 3);
92 gMC->Gspos("R3R3", 1, "ALIC", 0., 0., z + 1.525, 0, "ONLY");
93 // ***********************************************************
94 // DEFINE LEFT DISK#3 OF FMD
95 // ***********************************************************
106 gMC->Gsvolu("R1L3", "TUBE", idtmed[900], par, 3);
107 gMC->Gspos("R1L3", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
112 gMC->Gsvolu("R2L3", "TUBE", idtmed[901], par, 3);
113 gMC->Gspos("R2L3", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
116 par[1] = rout + 5.65;
118 gMC->Gsvolu("R3L3", "TUBE", idtmed[902], par, 3);
119 gMC->Gspos("R3L3", 1, "ALIC", 0., 0., z - 1.525, 0, "ONLY");
120 // ********************************************************
121 // DEFINE RIGHT DISK#2 OF FMD
122 // ********************************************************
133 gMC->Gsvolu("R1R2", "TUBE", idtmed[900], par, 3);
134 gMC->Gspos("R1R2", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
139 gMC->Gsvolu("R2R2", "TUBE", idtmed[901], par, 3);
140 gMC->Gspos("R2R2", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
143 par[1] = rout + 5.65;
145 gMC->Gsvolu("R3R2", "TUBE", idtmed[902], par, 3);
146 gMC->Gspos("R3R2", 1, "ALIC", 0., 0., z + 1.525, 0, "ONLY");
147 // ***********************************************************
148 // DEFINE LEFT DISK#2 OF FMD
149 // ***********************************************************
160 gMC->Gsvolu("R1L2", "TUBE", idtmed[900], par, 3);
161 gMC->Gspos("R1L2", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
166 gMC->Gsvolu("R2L2", "TUBE", idtmed[901], par, 3);
167 gMC->Gspos("R2L2", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
170 par[1] = rout + 5.65;
172 gMC->Gsvolu("R3L2", "TUBE", idtmed[902], par, 3);
173 gMC->Gspos("R3L2", 1, "ALIC", 0., 0., z - 1.525, 0, "ONLY");
174 // ********************************************************
175 // DEFINE RIGHT DISK#1 OF FMD
176 // ********************************************************
187 gMC->Gsvolu("R1R1", "TUBE", idtmed[900], par, 3);
188 gMC->Gspos("R1R1", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
193 gMC->Gsvolu("R2R1", "TUBE", idtmed[901], par, 3);
194 gMC->Gspos("R2R1", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
197 par[1] = rout + 5.65;
199 gMC->Gsvolu("R3R1", "TUBE", idtmed[902], par, 3);
200 gMC->Gspos("R3R1", 1, "ALIC", 0., 0., z + 1.525, 0, "ONLY");
201 // ***********************************************************
202 // DEFINE LEFT DISK#1 OF FMD
203 // ***********************************************************
214 gMC->Gsvolu("R1L1", "TUBE", idtmed[900], par, 3);
215 gMC->Gspos("R1L1", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
220 gMC->Gsvolu("R2L1", "TUBE", idtmed[901], par, 3);
221 gMC->Gspos("R2L1", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
224 par[1] = rout + 5.65;
226 gMC->Gsvolu("R3L1", "TUBE", idtmed[902], par, 3);
227 gMC->Gspos("R3L1", 1, "ALIC", 0., 0., z - 1.525, 0, "ONLY");
228 // ***********************************************************
229 // DEFINE LEFT DISK#4 OF FMD
230 // ***********************************************************
241 gMC->Gsvolu("R1L4", "TUBE", idtmed[900], par, 3);
242 gMC->Gspos("R1L4", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
247 gMC->Gsvolu("R2L4", "TUBE", idtmed[901], par, 3);
248 gMC->Gspos("R2L4", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
251 par[1] = rout + 5.65;
253 gMC->Gsvolu("R3L4", "TUBE", idtmed[902], par, 3);
254 gMC->Gspos("R3L4", 1, "ALIC", 0., 0., z - 1.525, 0, "ONLY");
257 //_____________________________________________________________________________
258 void AliFMDv0::DrawModule()
261 // Draw a shaded view of the FMD version 0
264 // Set everything unseen
265 gMC->Gsatt("*", "seen", -1);
267 // Set ALIC mother transparent
268 gMC->Gsatt("ALIC","SEEN",0);
270 // Set the volumes visible
271 gMC->Gsatt("R1R3","SEEN",1);
272 gMC->Gsatt("R2R3","SEEN",1);
273 gMC->Gsatt("R3R3","SEEN",1);
274 gMC->Gsatt("R1L3","SEEN",1);
275 gMC->Gsatt("R2L3","SEEN",1);
276 gMC->Gsatt("R3L3","SEEN",1);
277 gMC->Gsatt("R1R2","SEEN",1);
278 gMC->Gsatt("R2R2","SEEN",1);
279 gMC->Gsatt("R3R2","SEEN",1);
280 gMC->Gsatt("R1L2","SEEN",1);
281 gMC->Gsatt("R2L2","SEEN",1);
282 gMC->Gsatt("R3L2","SEEN",1);
283 gMC->Gsatt("R1R1","SEEN",1);
284 gMC->Gsatt("R2R1","SEEN",1);
285 gMC->Gsatt("R3R1","SEEN",1);
286 gMC->Gsatt("R1L1","SEEN",1);
287 gMC->Gsatt("R2L1","SEEN",1);
288 gMC->Gsatt("R3L1","SEEN",1);
289 gMC->Gsatt("R1L4","SEEN",1);
290 gMC->Gsatt("R2L4","SEEN",1);
291 gMC->Gsatt("R3L4","SEEN",1);
293 gMC->Gdopt("hide", "on");
294 gMC->Gdopt("shad", "on");
295 gMC->Gsatt("*", "fill", 7);
296 gMC->SetClipBox(".");
297 gMC->SetClipBox("*", 0, 1000, -1000, 1000, -1000, 1000);
299 gMC->Gdraw("alic", 40, 30, 0, 6, 9, .08, .08);
300 gMC->Gdhead(1111, "Forward Multiplicity Detector version 0");
301 gMC->Gdman(13, 9, "MAN");
304 //_____________________________________________________________________________
305 void AliFMDv0::CreateMaterials()
308 // Create Materials for version 0 of FMD
311 // Material for ring #1
312 Float_t ar1[8] = { 55.8,58.7,52.,47.9,16.,28.,207.2,27. };
313 Float_t zr1[8] = { 26.,28.,24.,22.,8.,14.,82.,13. };
314 Float_t wr1[8] = { .27,.081,.054,.045,.18,.25,.06,.06 };
315 // Material for ring #2
316 Float_t ar2[3] = { 55.8,27.,16. };
317 Float_t zr2[3] = { 26.,13.,8. };
318 Float_t wr2[3] = { .35,.34,.31 };
319 // Material for ring #3
320 Float_t ar3[3] = { 28.,27.,16. };
321 Float_t zr3[3] = { 14.,13.,8. };
322 Float_t wr3[3] = { .37,.33,.3 };
324 Float_t epsil, stmin, deemax, tmaxfd, stemax;
326 Int_t ISXFLD = gAlice->Field()->Integ();
327 Float_t SXMGMX = gAlice->Field()->Max();
331 AliMixture(1, "FMD_R1$", ar1, zr1, 2.69, 8, wr1);
335 AliMixture(2, "FMD_R2$", ar2, zr2, 2.63, 3, wr2);
339 AliMixture(3, "FMD_R3$", ar3, zr3, 3.15, 3, wr3);
340 // *******************************************************
341 // Defines tracking media parameters.
342 // *******************************************************
343 epsil = .001; // Tracking precision, DLS
344 stemax = -1.; // Maximum displacement for multiple scattering
345 tmaxfd = -20.; // Maximum angle due to field deflection
346 deemax = -.3; // Maximum fractional energy loss, DLS
348 // ********************************************************
349 AliMedium(1, "FMD_R1_L3 ", 1, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
350 AliMedium(2, "FMD_R2_L3 ", 2, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
351 AliMedium(3, "FMD_R3_L3 ", 3, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);