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">
60 AliMC* pMC = AliMC::GetMC();
65 Int_t *idtmed = fIdtmed->GetArray()-899;
67 // ********************************************************
68 // DEFINE RIGHT DISK#3 OF FMD
69 // ********************************************************
80 pMC->Gsvolu("R1R3", "TUBE", idtmed[900], par, 3);
81 pMC->Gspos("R1R3", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
86 pMC->Gsvolu("R2R3", "TUBE", idtmed[901], par, 3);
87 pMC->Gspos("R2R3", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
92 pMC->Gsvolu("R3R3", "TUBE", idtmed[902], par, 3);
93 pMC->Gspos("R3R3", 1, "ALIC", 0., 0., z + 1.525, 0, "ONLY");
94 // ***********************************************************
95 // DEFINE LEFT DISK#3 OF FMD
96 // ***********************************************************
107 pMC->Gsvolu("R1L3", "TUBE", idtmed[900], par, 3);
108 pMC->Gspos("R1L3", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
113 pMC->Gsvolu("R2L3", "TUBE", idtmed[901], par, 3);
114 pMC->Gspos("R2L3", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
117 par[1] = rout + 5.65;
119 pMC->Gsvolu("R3L3", "TUBE", idtmed[902], par, 3);
120 pMC->Gspos("R3L3", 1, "ALIC", 0., 0., z - 1.525, 0, "ONLY");
121 // ********************************************************
122 // DEFINE RIGHT DISK#2 OF FMD
123 // ********************************************************
134 pMC->Gsvolu("R1R2", "TUBE", idtmed[900], par, 3);
135 pMC->Gspos("R1R2", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
140 pMC->Gsvolu("R2R2", "TUBE", idtmed[901], par, 3);
141 pMC->Gspos("R2R2", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
144 par[1] = rout + 5.65;
146 pMC->Gsvolu("R3R2", "TUBE", idtmed[902], par, 3);
147 pMC->Gspos("R3R2", 1, "ALIC", 0., 0., z + 1.525, 0, "ONLY");
148 // ***********************************************************
149 // DEFINE LEFT DISK#2 OF FMD
150 // ***********************************************************
161 pMC->Gsvolu("R1L2", "TUBE", idtmed[900], par, 3);
162 pMC->Gspos("R1L2", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
167 pMC->Gsvolu("R2L2", "TUBE", idtmed[901], par, 3);
168 pMC->Gspos("R2L2", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
171 par[1] = rout + 5.65;
173 pMC->Gsvolu("R3L2", "TUBE", idtmed[902], par, 3);
174 pMC->Gspos("R3L2", 1, "ALIC", 0., 0., z - 1.525, 0, "ONLY");
175 // ********************************************************
176 // DEFINE RIGHT DISK#1 OF FMD
177 // ********************************************************
188 pMC->Gsvolu("R1R1", "TUBE", idtmed[900], par, 3);
189 pMC->Gspos("R1R1", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
194 pMC->Gsvolu("R2R1", "TUBE", idtmed[901], par, 3);
195 pMC->Gspos("R2R1", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
198 par[1] = rout + 5.65;
200 pMC->Gsvolu("R3R1", "TUBE", idtmed[902], par, 3);
201 pMC->Gspos("R3R1", 1, "ALIC", 0., 0., z + 1.525, 0, "ONLY");
202 // ***********************************************************
203 // DEFINE LEFT DISK#1 OF FMD
204 // ***********************************************************
215 pMC->Gsvolu("R1L1", "TUBE", idtmed[900], par, 3);
216 pMC->Gspos("R1L1", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
221 pMC->Gsvolu("R2L1", "TUBE", idtmed[901], par, 3);
222 pMC->Gspos("R2L1", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
225 par[1] = rout + 5.65;
227 pMC->Gsvolu("R3L1", "TUBE", idtmed[902], par, 3);
228 pMC->Gspos("R3L1", 1, "ALIC", 0., 0., z - 1.525, 0, "ONLY");
229 // ***********************************************************
230 // DEFINE LEFT DISK#4 OF FMD
231 // ***********************************************************
242 pMC->Gsvolu("R1L4", "TUBE", idtmed[900], par, 3);
243 pMC->Gspos("R1L4", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
248 pMC->Gsvolu("R2L4", "TUBE", idtmed[901], par, 3);
249 pMC->Gspos("R2L4", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
252 par[1] = rout + 5.65;
254 pMC->Gsvolu("R3L4", "TUBE", idtmed[902], par, 3);
255 pMC->Gspos("R3L4", 1, "ALIC", 0., 0., z - 1.525, 0, "ONLY");
258 //_____________________________________________________________________________
259 void AliFMDv0::DrawModule()
262 // Draw a shaded view of the FMD version 0
265 AliMC* pMC = AliMC::GetMC();
267 // Set everything unseen
268 pMC->Gsatt("*", "seen", -1);
270 // Set ALIC mother transparent
271 pMC->Gsatt("ALIC","SEEN",0);
273 // Set the volumes visible
274 pMC->Gsatt("R1R3","SEEN",1);
275 pMC->Gsatt("R2R3","SEEN",1);
276 pMC->Gsatt("R3R3","SEEN",1);
277 pMC->Gsatt("R1L3","SEEN",1);
278 pMC->Gsatt("R2L3","SEEN",1);
279 pMC->Gsatt("R3L3","SEEN",1);
280 pMC->Gsatt("R1R2","SEEN",1);
281 pMC->Gsatt("R2R2","SEEN",1);
282 pMC->Gsatt("R3R2","SEEN",1);
283 pMC->Gsatt("R1L2","SEEN",1);
284 pMC->Gsatt("R2L2","SEEN",1);
285 pMC->Gsatt("R3L2","SEEN",1);
286 pMC->Gsatt("R1R1","SEEN",1);
287 pMC->Gsatt("R2R1","SEEN",1);
288 pMC->Gsatt("R3R1","SEEN",1);
289 pMC->Gsatt("R1L1","SEEN",1);
290 pMC->Gsatt("R2L1","SEEN",1);
291 pMC->Gsatt("R3L1","SEEN",1);
292 pMC->Gsatt("R1L4","SEEN",1);
293 pMC->Gsatt("R2L4","SEEN",1);
294 pMC->Gsatt("R3L4","SEEN",1);
296 pMC->Gdopt("hide", "on");
297 pMC->Gdopt("shad", "on");
298 pMC->Gsatt("*", "fill", 7);
299 pMC->SetClipBox(".");
300 pMC->SetClipBox("*", 0, 1000, -1000, 1000, -1000, 1000);
302 pMC->Gdraw("alic", 40, 30, 0, 6, 9, .08, .08);
303 pMC->Gdhead(1111, "Forward Multiplicity Detector version 0");
304 pMC->Gdman(13, 9, "MAN");
307 //_____________________________________________________________________________
308 void AliFMDv0::CreateMaterials()
311 // Create Materials for version 0 of FMD
314 // Material for ring #1
315 Float_t ar1[8] = { 55.8,58.7,52.,47.9,16.,28.,207.2,27. };
316 Float_t zr1[8] = { 26.,28.,24.,22.,8.,14.,82.,13. };
317 Float_t wr1[8] = { .27,.081,.054,.045,.18,.25,.06,.06 };
318 // Material for ring #2
319 Float_t ar2[3] = { 55.8,27.,16. };
320 Float_t zr2[3] = { 26.,13.,8. };
321 Float_t wr2[3] = { .35,.34,.31 };
322 // Material for ring #3
323 Float_t ar3[3] = { 28.,27.,16. };
324 Float_t zr3[3] = { 14.,13.,8. };
325 Float_t wr3[3] = { .37,.33,.3 };
327 Float_t epsil, stmin, deemax, tmaxfd, stemax;
329 Int_t ISXFLD = gAlice->Field()->Integ();
330 Float_t SXMGMX = gAlice->Field()->Max();
334 AliMixture(1, "FMD_R1$", ar1, zr1, 2.69, 8, wr1);
338 AliMixture(2, "FMD_R2$", ar2, zr2, 2.63, 3, wr2);
342 AliMixture(3, "FMD_R3$", ar3, zr3, 3.15, 3, wr3);
343 // *******************************************************
344 // Defines tracking media parameters.
345 // *******************************************************
346 epsil = .001; // Tracking precision, DLS
347 stemax = -1.; // Maximum displacement for multiple scattering
348 tmaxfd = -20.; // Maximum angle due to field deflection
349 deemax = -.3; // Maximum fractional energy loss, DLS
351 // ********************************************************
352 AliMedium(1, "FMD_R1_L3 ", 1, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
353 AliMedium(2, "FMD_R2_L3 ", 2, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
354 AliMedium(3, "FMD_R3_L3 ", 3, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);