]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ITS/AliITSv4.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Medium array is now private for each module
[u/mrichter/AliRoot.git] / ITS / AliITSv4.cxx
CommitLineData
fe4da5cc 1///////////////////////////////////////////////////////////////////////////////
2// //
3// Inner Traking System version 4 //
4// //
5//Begin_Html
6/*
1439f98e 7<img src="picts/AliITSv4Class.gif">
fe4da5cc 8</pre>
9<br clear=left>
10<font size=+2 color=red>
11<p>The responsible person for this module is
12<a href="mailto:roberto.barbera@ct.infn.it">Roberto Barbera</a>.
13</font>
14<pre>
15*/
16//End_Html
17// //
18///////////////////////////////////////////////////////////////////////////////
19
20#include <TMath.h>
21#include <TRandom.h>
22#include <TVector.h>
23#include "AliITSv4.h"
24#include "AliRun.h"
25#include "AliCallf77.h"
26
27#include "TGeant3.h"
28
29# define its_geo4 its_geo4_
30extern "C" void type_of_call its_geo4(Int_t&, Int_t*);
31
32ClassImp(AliITSv4)
33
34//_____________________________________________________________________________
35AliITSv4::AliITSv4() : AliITS()
36{
37 //
38 // Default constructor for the ITS
39 //
40 fMinorVersion=1;
41}
42
43//_____________________________________________________________________________
44AliITSv4::AliITSv4(const char *name, const char *title)
45 : AliITS(name, title)
46{
47 //
48 // Standard constructor
49 //
50 fMinorVersion=1;
51}
52
53//_____________________________________________________________________________
54void AliITSv4::CreateGeometry()
55{
56 //
57 // Create geometry for its version 4
58 //
59 AliMC* pMC = AliMC::GetMC();
60#ifndef NEW
61 Int_t ivers=4*10+fMinorVersion;
ad51aeb0 62 its_geo4(ivers, fIdtmed->GetArray()-200;
fe4da5cc 63#else
64// *** DEFINITION OF THE GEOMETRY OF THE ITS ***
65// *** RB APR-1990 CATANIA ***
66// ORIGIN : ROBERTO BARBERA
67
68 Float_t xx[13] = { 0.,0.,-4.824,-4.833,-22.167,-22.585,-28.07,-27.626,
69 -38.139,-19.749,-13.449,-14.726,0. };
70 Float_t yy[13] = { 0.,34.028,34.064,
71 32.594,29.984,30.914,29.18,27.777,22.522,-7.918,-4.769,-2.216,0. };
72 Float_t xbeg[12] = { 0.,-.497,-4.827,-5.425,-22.395,-23.009,-27.919,
73 -28.161,-37.859,-19.337,-13.673,-13.773 };
74 Float_t ybeg[12] = { .43,34.032,33.564,
75 32.505,30.493,30.781,28.702,27.509,22.058,-7.711,-4.322,-2.073 };
76 Float_t xend[12] = { 0.,-4.323,-4.83,-21.643,-22.395,-27.593,-27.807,
77 -37.655,-19.988,-13.897,-14.294,-.425 };
78 Float_t yend[12] = { 33.531,34.06,33.192,
79 30.064,30.493,29.33,28.347,22.764,-7.523,-4.992,-3.077,-.064 };
80 Float_t xarc[12] = { -19.56,-14.12,-13.669,-.5,-.5,-4.327,-5.529,-21.746,
81 -22.858,-27.442,-28.474,-37.431 };
82 Float_t yarc[12] = { -7.264,-4.545,-2.765,
83 .431,33.531,33.56,33.197,30.756,30.304,28.854,28.136,22.316 };
84 Float_t rarc[12] = { .5,.5,.7,.5,.5,.5,.7,.7,.5,.5,.7,.5 };
85 Float_t rr = 4.08332;
86 Float_t tteta = 63.;
87 Float_t pphi = -27.;
88 Float_t gteta = 89.4;
89
90 static Double_t biga, bigb;
91 Float_t dcei[3], dela[3], dchi[3], dpcb[3], darc[5], dfra[10], dcer[3], dkap[3], dpla[3], xccc, yccc, aphi, dcop[3], dtra[3], dsil[3],
92 atheta1213, dbus[3], dtub[3], dwat[3], depx[3], dits[3], atheta1112, atheta1011, dsup[3], xtra[8], ytra[8], ztra[8], dsrv[3];
93 static Double_t biga1, bigb1;
94 Float_t runo, xpos, ypos, zpos, rtwo, aphi1, aphi2, dtra1[3], dtra2[3], dtra3[3], dtra4[3], dbox1[3], dbox2[3];
95 Int_t jbox1, jbox2;
96 Float_t xtra1[6], ytra1[6], ztra1[6];
97 Int_t i;
98 Float_t xpos1, ypos1;
99 Int_t j;
100 Float_t angle, dcone[5], dtube[3], dpgon[10];
101 Float_t rzero, xzero, yzero;
102 static Double_t coeffa, coeffb, coeffc;
103 Int_t idrotm[5250];
104 Float_t atheta, offset;
105 Float_t offset1, offset2, dgh[15];
106 extern Double_t atg_(Float_t *, real *);
107 Float_t xcc, ycc, sep, atheta12, atheta23, atheta34, atheta45, atheta56, atheta67, atheta78, atheta89, yos, xxm, dal1[3], dal2[3];
108
109 static Double_t xcc1, ycc1, xcc2, ycc2;
110 Float_t atheta910;
111
ad51aeb0 112 Int_t *idtmed = fIdtmed->GetArray()-199;
fe4da5cc 113
114#define natra_ref(a_0,a_1) &natra[(a_1)*4 + a_0 - 4]
115#define natra1_ref(a_0,a_1) &natra1[(a_1)*4 + a_0 - 4]
116#define natra2_ref(a_0,a_1) &natra2[(a_1)*4 + a_0 - 4]
117#define natra3_ref(a_0,a_1) &natra3[(a_1)*4 + a_0 - 4]
118#define natra4_ref(a_0,a_1) &natra4[(a_1)*4 + a_0 - 4]
119
120
121// --- Define ghost volume containing the whole ITS and fill it with air
122// or vacuum
123
124// 89.4 = 90.-0.6
125 dgh[0] = 0.;
126 dgh[1] = 360.;
127 dgh[2] = 4.;
128 dgh[3] = -70.;
129 dgh[4] = 49.999;
130 dgh[5] = 49.999;
131 dgh[6] = -25.;
132 dgh[7] = 3.;
133 dgh[8] = 49.999;
134 dgh[9] = 25.;
135 dgh[10] = 3.;
136 dgh[11] = 49.999;
137 dgh[12] = 70.;
138 dgh[13] = 49.999;
139 dgh[14] = 49.999;
140 geant->Gsvolu("ITSV", "PCON", idtmed[275], dgh, 15);
141
142// --- Place the ghost volume in its mother volume (ALIC) and make it
143// invisible
144
145 geant->Gspos("ITSV", 1, "ALIC", 0., 0., 0., 0, "ONLY");
146 geant->Gsatt("ITSV", "SEEN", 0);
147
148//************************************************************************
149//* *
150//* P I X E L S *
151//* =========== *
152//* *
153//************************************************************************
154
155// GOTO 2345 ! skip ITS layer no. 1 and 2
156
157//--- Define ghost volume containing the Pixel Detectors and fill it with air
158// or vacuum
159
160 xxm = 1.0444222222222224;
161 dgh[0] = 0.;
162 dgh[1] = 360.;
163 dgh[2] = 4.;
164 dgh[3] = -25. - 5.9900000000000002 / xxm;
165 dgh[4] = 9.;
166 dgh[5] = 9.;
167 dgh[6] = -25.;
168 dgh[7] = 3.01;
169 dgh[8] = 9.;
170 dgh[9] = 25.;
171 dgh[10] = 3.01;
172 dgh[11] = 9.;
173 dgh[12] = 5.9900000000000002 / xxm + 25.;
174 dgh[13] = 9.;
175 dgh[14] = 9.;
176 geant->Gsvolu("IT12", "PCON", idtmed[275], dgh, 15);
177
178// --- Place the ghost volume in its mother volume (ITSV) and make it
179// invisible
180
181 geant->Gspos("IT12", 1, "ITSV", 0., 0., 0., 0, "ONLY");
182 geant->Gsatt("IT12", "SEEN", 0);
183
184// --- Define a ghost volume containing a single element of layer #1
185// and fill it with air or vacuum
186
187 dbox1[0] = .022499999999999999;
188 dbox1[1] = .79;
189 dbox1[2] = 12.67;
190 geant->Gsvolu("IPV1", "BOX ", idtmed[203], dbox1, 3);
191 ;
192
193//--- Divide each element of layer #1 in three ladders along the beam direction
194
195 geant->Gsdvn("IPB1", "IPV1", 3, 3);
196
197// --- Make the ghost volumes invisible
198
199 geant->Gsatt("IPV1", "SEEN", 0);
200 geant->Gsatt("IPB1", "SEEN", 0);
201
202// --- Define a volume containing the chip of pixels (silicon, layer #1)
203
204 dchi[0] = .005;
205 dchi[1] = .79;
206 dchi[2] = dbox1[2] / 3.;
207 geant->Gsvolu("ICH1", "BOX ", idtmed[200], dchi, 3);
208
209// --- Define a volume containing the bus of pixels (silicon, layer #1)
210
211 dbus[0] = .01;
212 dbus[1] = .64;
213 dbus[2] = 4.19;
214 geant->Gsvolu("IBU1", "BOX ", idtmed[201], dbus, 3);
215
216// --- Define a volume containing the sensitive part of pixels
217// (silicon, layer #1)
218
219 dits[0] = .0075;
220 dits[1] = .64;
221 dits[2] = 4.19;
222 geant->Gsvolu("ITS1", "BOX ", idtmed[199], dits, 3);
223
224// --- Place the chip into its mother (IPB1)
225
226 xpos = dbox1[0] - dchi[0];
227 ypos = 0.;
228 zpos = 0.;
229 geant->Gspos("ICH1", 1, "IPB1", xpos, ypos, zpos, 0, "ONLY");
230
231// --- Place the sensitive volume into its mother (IPB1)
232
233 xpos = dbox1[0] - dchi[0] * 2. - dits[0];
234 ypos = dchi[1] - dits[1];
235 zpos = -(dchi[2] - dits[2]);
236 geant->Gspos("ITS1", 1, "IPB1", xpos, ypos, zpos, 0, "ONLY");
237
238// --- Place the bus into its mother (IPB1)
239
240 xpos = dbox1[0] - dchi[0] * 2. - dits[0] * 2. - dbus[0];
241 ypos = dchi[1] - dbus[1];
242 zpos = -(dchi[2] - dbus[2]);
243 geant->Gspos("IBU1", 1, "IPB1", xpos, ypos, zpos, 0, "ONLY");
244
245// --- Define a ghost volume containing a single element of layer #2
246// and fill it with air or vacuum
247
248 dbox2[0] = .022499999999999999;
249 dbox2[1] = .79;
250 dbox2[2] = 16.91;
251 geant->Gsvolu("IPV2", "BOX ", idtmed[203], dbox2, 3);
252 ;
253
254//--- Divide each element of layer #2 in four ladders along the beam direction
255
256 geant->Gsdvn("IPB2", "IPV2", 4, 3);
257
258// --- Make the ghost volumes invisible
259
260 geant->Gsatt("IPV2", "SEEN", 0);
261 geant->Gsatt("IPB2", "SEEN", 0);
262
263// --- Define a volume containing the chip of pixels (silicon, layer #2)
264
265 dchi[0] = .005;
266 dchi[1] = .79;
267 dchi[2] = dbox2[2] / 4.;
268 geant->Gsvolu("ICH2", "BOX ", idtmed[200], dchi, 3);
269
270// --- Define a volume containing the bus of pixels (silicon, layer #2)
271
272 dbus[0] = .01;
273 dbus[1] = .64;
274 dbus[2] = 4.19;
275 geant->Gsvolu("IBU2", "BOX ", idtmed[201], dbus, 3);
276
277// --- Define a volume containing the sensitive part of pixels
278// (silicon, layer #2)
279
280 dits[0] = .0075;
281 dits[1] = .64;
282 dits[2] = 4.19;
283 geant->Gsvolu("ITS2", "BOX ", idtmed[199], dits, 3);
284
285// --- Place the chip into its mother (IPB2)
286
287 xpos = dbox1[0] - dbus[0] * 2. - dits[0] * 2. - dchi[0];
288 ypos = 0.;
289 zpos = 0.;
290 geant->Gspos("ICH2", 1, "IPB2", xpos, ypos, zpos, 0, "ONLY");
291
292// --- Place the sensitive volume into its mother (IPB2)
293
294 xpos = dbox1[0] - dbus[0] * 2. - dits[0];
295 ypos = -(dchi[1] - dits[1]);
296 zpos = -(dchi[2] - dits[2]);
297 geant->Gspos("ITS2", 1, "IPB2", xpos, ypos, zpos, 0, "ONLY");
298
299// --- Place the bus into its mother (IPB2)
300
301 xpos = dbox1[0] - dbus[0];
302 ypos = -(dchi[1] - dbus[1]);
303 zpos = -(dchi[2] - dbus[2]);
304 geant->Gspos("IBU2", 1, "IPB2", xpos, ypos, zpos, 0, "ONLY");
305
306// --- Define a generic segment of an element of the mechanical support
307
308 dsup[0] = 0.;
309 dsup[1] = 0.;
310 dsup[2] = 0.;
311 geant->Gsvolu("SPIX", "BOX ", idtmed[202], dsup, 0);
312
313// --- Define a generic arc of an element of the mechanical support
314
315 darc[0] = 0.;
316 darc[1] = 0.;
317 darc[2] = 0.;
318 geant->Gsvolu("SARC", "TUBS", idtmed[202], darc, 0);
319
320// --- Define the mechanical supports of layers #1 and #2 and place the
321// elements of the layers in it
322
323 jbox1 = 0;
324// counter over the number of elements of layer #1 (
325 jbox2 = 0;
326
327// counter over the number of elements of layer #2 (
328 for (i = 1; i <= 10; ++i) {
329
330// --- Place part # 1-2 (see sketch)
331
332// number of carbon fiber supports (see sketch)
333 offset1 = -27.;
334 dsup[0] = .01;
335 dsup[1] = TMath::Sqrt((xend[0] - xbeg[0]) * (xend[0] - xbeg[0]) + (yend[0] - ybeg[0]) * (yend[0] - ybeg[0])) / 20.;
336 dsup[2] = 25.;
337 xcc = (xx[0] + xx[1]) / 20.;
338 ycc = (yy[0] + yy[1]) / 20.;
339 xccc = (xbeg[0] + xend[0]) / 20.;
340 yccc = (ybeg[0] + yend[0]) / 20.;
341 if (xx[0] == xx[1]) {
342 offset2 = 0.;
343 } else {
344 r2 = xx[1] - xx[0];
345 offset2 = atg_(&yy[1] - yy[0], &r2) * raddeg - 90.;
346 }
347 aphi = (pphi + (i-1) * 36.) * degrad;
348 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
349 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
350 xpos1 = xccc * TMath::Cos(aphi) - yccc * TMath::Sin(aphi) + xzero;
351 ypos1 = xccc * TMath::Sin(aphi) + yccc * TMath::Cos(aphi) + yzero;
352 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
353 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
354 zpos = 0.;
355 atheta12 = (i-1) * 36. + offset1 + offset2 - gteta;
356 AliMatrix(idrotm[(i-1) * 13 + 1100], 90., atheta12, 90., atheta12 + 90., 0., 0.);
357 geant->Gsposp("SPIX", (i-1) * 13 + 1, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1100], "ONLY", dsup, 3);
358
359// --- Place part # 2-3 (see sketch)
360
361 offset1 = -27.;
362 dsup[0] = .01;
363 dsup[1] = TMath::Sqrt((xend[1] - xbeg[1]) * (xend[1] - xbeg[1]) + (yend[1] - ybeg[1]) * (yend[1] - ybeg[1])) / 20.;
364 dsup[2] = 25.;
365 xcc = (xx[1] + xx[2]) / 20.;
366 ycc = (yy[1] + yy[2]) / 20.;
367 xccc = (xbeg[1] + xend[1]) / 20.;
368 yccc = (ybeg[1] + yend[1]) / 20.;
369 if (xx[1] == xx[2]) {
370 offset2 = 0.;
371 } else {
372 r2 = xx[2] - xx[1];
373 offset2 = atg_(&yy[2] - yy[1], &r2) * raddeg - 90.;
374 }
375 aphi = (pphi + (i-1) * 36.) * degrad;
376 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
377 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
378 xpos1 = xccc * TMath::Cos(aphi) - yccc * TMath::Sin(aphi) + xzero;
379 ypos1 = xccc * TMath::Sin(aphi) + yccc * TMath::Cos(aphi) + yzero;
380 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
381 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
382 zpos = 0.;
383 atheta23 = (i-1) * 36. + offset1 + offset2 - gteta;
384 AliMatrix(idrotm[(i-1) * 13 + 1101], 90., atheta23, 90., atheta23 + 90., 0., 0.);
385 geant->Gsposp("SPIX", (i-1) * 13 + 2, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1101], "ONLY", dsup, 3);
386
387// --- Place an element of layer #2
388
389 biga = (yy[2] - yy[1]) / (xx[2] - xx[1]);
390 bigb = (xx[2] * yy[1] - xx[1] * yy[2]) / (xx[2] - xx[1]) / 10.;
391 coeffa = biga * biga + 1.;
392 coeffb = biga * bigb - biga * ycc - xcc;
393 coeffc = xcc * xcc + ycc * ycc - ycc * 2. * bigb + bigb * bigb - .0022212369;
394 xcc1 = (-coeffb + TMath::Sqrt(coeffb * coeffb - coeffa * coeffc)) / coeffa;
395 ycc1 = biga * xcc1 + bigb;
396 biga1 = -1. / biga;
397 bigb1 = xcc1 / biga + ycc1;
398 coeffa = biga1 * biga1 + 1.;
399 coeffb = biga1 * bigb1 - biga1 * ycc1 - xcc1;
400 coeffc = xcc1 * xcc1 + ycc1 * ycc1 - ycc1 * 2. * bigb1 + bigb1 * bigb1 - (dsup[0] + dbox2[0]) * (dsup[0] + dbox2[0]);
401 xcc2 = (-coeffb + TMath::Sqrt(coeffb * coeffb - coeffa * coeffc)) / coeffa;
402 ycc2 = biga1 * xcc2 + bigb1;
403 xpos1 = xcc2 * TMath::Cos(aphi) - ycc2 * TMath::Sin(aphi) + xzero;
404 ypos1 = xcc2 * TMath::Sin(aphi) + ycc2 * TMath::Cos(aphi) + yzero;
405 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
406 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
407 zpos = 0.;
408 ++jbox2;
409 geant->Gspos("IPV2", jbox2, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1101], "ONLY");
410
411// --- Place part # 3-4 (see sketch)
412
413 offset1 = -27.;
414 dsup[0] = .01;
415 dsup[1] = TMath::Sqrt((xend[2] - xbeg[2]) * (xend[2] - xbeg[2]) + (yend[2] - ybeg[2]) * (yend[2] - ybeg[2])) / 20.;
416 dsup[2] = 25.;
417 xcc = (xx[1] + xx[2]) / 20.;
418 ycc = (yy[1] + yy[2]) / 20.;
419 xccc = (xbeg[2] + xend[2]) / 20.;
420 yccc = (ybeg[2] + yend[2]) / 20.;
421 if (xx[2] == xx[3]) {
422 offset2 = 0.;
423 } else {
424 r2 = xx[3] - xx[2];
425 offset2 = atg_(&yy[3] - yy[2], &r2) * raddeg - 90.;
426 }
427 aphi = (pphi + (i-1) * 36.) * degrad;
428 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
429 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
430 xpos1 = xccc * TMath::Cos(aphi) - yccc * TMath::Sin(aphi) + xzero;
431 ypos1 = xccc * TMath::Sin(aphi) + yccc * TMath::Cos(aphi) + yzero;
432 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
433 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
434 zpos = 0.;
435 atheta34 = (i-1) * 36. + offset1 + offset2 - gteta;
436 AliMatrix(idrotm[(i-1) * 13 + 1102], 90., atheta34, 90., atheta34 + 90., 0., 0.);
437 geant->Gsposp("SPIX", (i-1) * 13 + 3, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1102], "ONLY", dsup, 3);
438
439// --- Place part # 4-5 (see sketch)
440
441 offset1 = -27.;
442 dsup[0] = .01;
443 dsup[1] = TMath::Sqrt((xend[3] - xbeg[3]) * (xend[3] - xbeg[3]) + (yend[3] - ybeg[3]) * (yend[3] - ybeg[3])) / 20.;
444 dsup[2] = 25.;
445 xcc = (xx[3] + xx[4]) / 20.;
446 ycc = (yy[3] + yy[4]) / 20.;
447 xccc = (xbeg[3] + xend[3]) / 20.;
448 yccc = (ybeg[3] + yend[3]) / 20.;
449 if (xx[3] == xx[4]) {
450 offset2 = 0.;
451 } else {
452 r2 = xx[4] - xx[3];
453 offset2 = atg_(&yy[4] - yy[3], &r2) * raddeg - 90.;
454 }
455 aphi = (pphi + (i-1) * 36.) * degrad;
456 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
457 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
458 xpos1 = xccc * TMath::Cos(aphi) - yccc * TMath::Sin(aphi) + xzero;
459 ypos1 = xccc * TMath::Sin(aphi) + yccc * TMath::Cos(aphi) + yzero;
460 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
461 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
462 zpos = 0.;
463 atheta45 = (i-1) * 36. + offset1 + offset2 - gteta;
464 AliMatrix(idrotm[(i-1) * 13 + 1103], 90., atheta45, 90., atheta45 + 90., 0., 0.);
465 geant->Gsposp("SPIX", (i-1) * 13 + 4, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1103], "ONLY", dsup, 3);
466
467// --- Place an element of layer #2
468
469 biga = (yy[4] - yy[3]) / (xx[4] - xx[3]);
470 bigb = (xx[4] * yy[3] - xx[3] * yy[4]) / (xx[4] - xx[3]) / 10.;
471 coeffa = biga * biga + 1.;
472 coeffb = biga * bigb - biga * ycc - xcc;
473 coeffc = xcc * xcc + ycc * ycc - ycc * 2. * bigb + bigb * bigb + 0.;
474 xcc1 = xccc;
475 ycc1 = yccc;
476 biga1 = -1. / biga;
477 bigb1 = xcc1 / biga + ycc1;
478 coeffa = biga1 * biga1 + 1.;
479 coeffb = biga1 * bigb1 - biga1 * ycc1 - xcc1;
480 coeffc = xcc1 * xcc1 + ycc1 * ycc1 - ycc1 * 2. * bigb1 + bigb1 * bigb1 - (dsup[0] + dbox2[0]) * (dsup[0] + dbox2[0]);
481 xcc2 = (-coeffb - TMath::Sqrt(coeffb * coeffb - coeffa * coeffc)) / coeffa;
482 ycc2 = biga1 * xcc2 + bigb1;
483 xpos1 = xcc2 * TMath::Cos(aphi) - ycc2 * TMath::Sin(aphi) + xzero;
484 ypos1 = xcc2 * TMath::Sin(aphi) + ycc2 * TMath::Cos(aphi) + yzero;
485 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
486 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
487 zpos = 0.;
488 ++jbox2;
489 geant->Gspos("IPV2", jbox2, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1103], "ONLY");
490
491// --- Place part # 5-6 (see sketch)
492
493 offset1 = -27.;
494 dsup[0] = .01;
495 dsup[1] = TMath::Sqrt((xend[4] - xbeg[4]) * (xend[4] - xbeg[4]) + (yend[4] - ybeg[4]) * (yend[4] - ybeg[4])) / 20.;
496 dsup[2] = 25.;
497 xcc = (xx[4] + xx[5]) / 20.;
498 ycc = (yy[4] + yy[5]) / 20.;
499 xccc = (xbeg[4] + xend[4]) / 20.;
500 yccc = (ybeg[4] + yend[4]) / 20.;
501 if (xx[4] == xx[5]) {
502 offset2 = 0.;
503 } else {
504 r2 = xx[5] - xx[4];
505 offset2 = atg_(&yy[5] - yy[4], &r2) * raddeg - 90.;
506 }
507 aphi = (pphi + (i-1) * 36.) * degrad;
508 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
509 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
510 xpos1 = xccc * TMath::Cos(aphi) - yccc * TMath::Sin(aphi) + xzero;
511 ypos1 = xccc * TMath::Sin(aphi) + yccc * TMath::Cos(aphi) + yzero;
512 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
513 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
514 zpos = 0.;
515 atheta56 = (i-1) * 36. + offset1 + offset2 - gteta;
516 AliMatrix(idrotm[(i-1) * 13 + 1104], 90., atheta56, 90., atheta56 + 90., 0., 0.);
517 geant->Gsposp("SPIX", (i-1) * 13 + 5, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1104], "ONLY", dsup, 3);
518
519// --- Place part # 6-7 (see sketch)
520
521 offset1 = -27.;
522 dsup[0] = .01;
523 dsup[1] = TMath::Sqrt((xend[5] - xbeg[5]) * (xend[5] - xbeg[5]) + (yend[5] - ybeg[5]) * (yend[5] - ybeg[5])) / 20.;
524 dsup[2] = 25.;
525 xcc = (xx[5] + xx[6]) / 20.;
526 ycc = (yy[5] + yy[6]) / 20.;
527 xccc = (xbeg[5] + xend[5]) / 20.;
528 yccc = (ybeg[5] + yend[5]) / 20.;
529 if (xx[5] == xx[6]) {
530 offset2 = 0.;
531 } else {
532 r2 = xx[6] - xx[5];
533 offset2 = atg_(&yy[6] - yy[5], &r2) * raddeg - 90.;
534 }
535 aphi = (pphi + (i-1) * 36.) * degrad;
536 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
537 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
538 xpos1 = xccc * TMath::Cos(aphi) - yccc * TMath::Sin(aphi) + xzero;
539 ypos1 = xccc * TMath::Sin(aphi) + yccc * TMath::Cos(aphi) + yzero;
540 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
541 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
542 zpos = 0.;
543 atheta67 = (i-1) * 36. + offset1 + offset2 - gteta;
544 AliMatrix(idrotm[(i-1) * 13 + 1105], 90., atheta67, 90., atheta67 + 90., 0., 0.);
545 geant->Gsposp("SPIX", (i-1) * 13 + 6, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1105], "ONLY", dsup, 3);
546
547// --- Place an element of layer #2
548
549 biga = (yy[6] - yy[5]) / (xx[6] - xx[5]);
550 bigb = (xx[6] * yy[5] - xx[5] * yy[6]) / (xx[6] - xx[5]) / 10.;
551 coeffa = biga * biga + 1.;
552 coeffb = biga * bigb - biga * ycc - xcc;
553 coeffc = xcc * xcc + ycc * ycc - ycc * 2. * bigb + bigb * bigb + 0.;
554 xcc1 = xccc;
555 ycc1 = yccc;
556 biga1 = -1. / biga;
557 bigb1 = xcc1 / biga + ycc1;
558 coeffa = biga1 * biga1 + 1.;
559 coeffb = biga1 * bigb1 - biga1 * ycc1 - xcc1;
560 coeffc = xcc1 * xcc1 + ycc1 * ycc1 - ycc1 * 2. * bigb1 + bigb1 * bigb1 - (dsup[0] + dbox2[0]) * (dsup[0] + dbox2[0]);
561 xcc2 = (-coeffb - TMath::Sqrt(coeffb * coeffb - coeffa * coeffc)) / coeffa;
562 ycc2 = biga1 * xcc2 + bigb1;
563 xpos1 = xcc2 * TMath::Cos(aphi) - ycc2 * TMath::Sin(aphi) + xzero;
564 ypos1 = xcc2 * TMath::Sin(aphi) + ycc2 * TMath::Cos(aphi) + yzero;
565 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
566 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
567 zpos = 0.;
568 ++jbox2;
569 geant->Gspos("IPV2", jbox2, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1105], "ONLY");
570
571// --- Place part # 7-8 (see sketch)
572
573 offset1 = -27.;
574 dsup[0] = .01;
575 dsup[1] = TMath::Sqrt((xend[6] - xbeg[6]) * (xend[6] - xbeg[6]) + (yend[6] - ybeg[6]) * (yend[6] - ybeg[6])) / 20.;
576 dsup[2] = 25.;
577 xcc = (xx[6] + xx[7]) / 20.;
578 ycc = (yy[6] + yy[7]) / 20.;
579 xccc = (xbeg[6] + xend[6]) / 20.;
580 yccc = (ybeg[6] + yend[6]) / 20.;
581 if (xx[6] == xx[7]) {
582 offset2 = 0.;
583 } else {
584 r2 = xx[7] - xx[6];
585 offset2 = atg_(&yy[7] - yy[6], &r2) * raddeg - 90.;
586 }
587 aphi = (pphi + (i-1) * 36.) * degrad;
588 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
589 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
590 xpos1 = xccc * TMath::Cos(aphi) - yccc * TMath::Sin(aphi) + xzero;
591 ypos1 = xccc * TMath::Sin(aphi) + yccc * TMath::Cos(aphi) + yzero;
592 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
593 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
594 zpos = 0.;
595 atheta78 = (i-1) * 36. + offset1 + offset2 - gteta;
596 AliMatrix(idrotm[(i-1) * 13 + 1106], 90., atheta78, 90., atheta78 + 90., 0., 0.);
597 geant->Gsposp("SPIX", (i-1) * 13 + 7, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1106], "ONLY", dsup, 3);
598
599// --- Place part # 8-9 (see sketch)
600
601 offset1 = -27.;
602 dsup[0] = .01;
603 dsup[1] = TMath::Sqrt((xend[7] - xbeg[7]) * (xend[7] - xbeg[7]) + (yend[7] - ybeg[7]) * (yend[7] - ybeg[7])) / 20.;
604 dsup[2] = 25.;
605 xcc = (xx[7] + xx[8]) / 20.;
606 ycc = (yy[7] + yy[8]) / 20.;
607 xccc = (xbeg[7] + xend[7]) / 20.;
608 yccc = (ybeg[7] + yend[7]) / 20.;
609 if (xx[1] == xx[2]) {
610 offset2 = 0.;
611 } else {
612 r2 = xx[8] - xx[7];
613 offset2 = atg_(&yy[8] - yy[7], &r2) * raddeg - 90.;
614 }
615 aphi = (pphi + (i-1) * 36.) * degrad;
616 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
617 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
618 xpos1 = xccc * TMath::Cos(aphi) - yccc * TMath::Sin(aphi) + xzero;
619 ypos1 = xccc * TMath::Sin(aphi) + yccc * TMath::Cos(aphi) + yzero;
620 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
621 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
622 zpos = 0.;
623 atheta89 = (i-1) * 36. + offset1 + offset2 - gteta;
624 AliMatrix(idrotm[(i-1) * 13 + 1107], 90., atheta89, 90., atheta89 + 90., 0., 0.);
625 geant->Gsposp("SPIX", (i-1) * 13 + 8, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1107], "ONLY", dsup, 3);
626
627// --- Place an element of layer #2
628
629 biga = (yy[8] - yy[7]) / (xx[8] - xx[7]);
630 bigb = (xx[8] * yy[7] - xx[7] * yy[8]) / (xx[8] - xx[7]) / 10.;
631 coeffa = biga * biga + 1.;
632 coeffb = biga * bigb - biga * ycc - xcc;
633 coeffc = xcc * xcc + ycc * ycc - ycc * 2. * bigb + bigb * bigb - .055526209599999994;
634 xcc1 = (-coeffb - TMath::Sqrt(coeffb * coeffb - coeffa * coeffc)) / coeffa;
635 ycc1 = biga * xcc1 + bigb;
636 biga1 = -1. / biga;
637 bigb1 = xcc1 / biga + ycc1;
638 coeffa = biga1 * biga1 + 1.;
639 coeffb = biga1 * bigb1 - biga1 * ycc1 - xcc1;
640 coeffc = xcc1 * xcc1 + ycc1 * ycc1 - ycc1 * 2. * bigb1 + bigb1 * bigb1 - (dsup[0] + dbox2[0]) * (dsup[0] + dbox2[0]);
641 xcc2 = (-coeffb - TMath::Sqrt(coeffb * coeffb - coeffa * coeffc)) / coeffa;
642 ycc2 = biga1 * xcc2 + bigb1;
643 xpos1 = xcc2 * TMath::Cos(aphi) - ycc2 * TMath::Sin(aphi) + xzero;
644 ypos1 = xcc2 * TMath::Sin(aphi) + ycc2 * TMath::Cos(aphi) + yzero;
645 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
646 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
647 zpos = 0.;
648 ++jbox2;
649 geant->Gspos("IPV2", jbox2, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1107], "ONLY");
650
651// --- Place part # 9-10 (see sketch)
652
653 offset1 = -27.;
654 dsup[0] = .01;
655 dsup[1] = TMath::Sqrt((xend[8] - xbeg[8]) * (xend[8] - xbeg[8]) + (yend[8] - ybeg[8]) * (yend[8] - ybeg[8])) / 20.;
656 dsup[2] = 25.;
657 xcc = (xx[8] + xx[9]) / 20.;
658 ycc = (yy[8] + yy[9]) / 20.;
659 xccc = (xbeg[8] + xend[8]) / 20.;
660 yccc = (ybeg[8] + yend[8]) / 20.;
661 if (xx[8] == xx[9]) {
662 offset2 = 0.;
663 } else {
664 r2 = xx[9] - xx[8];
665 offset2 = atg_(&yy[9] - yy[8], &r2) * raddeg - 90.;
666 }
667 aphi = (pphi + (i-1) * 36.) * degrad;
668 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
669 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
670 xpos1 = xccc * TMath::Cos(aphi) - yccc * TMath::Sin(aphi) + xzero;
671 ypos1 = xccc * TMath::Sin(aphi) + yccc * TMath::Cos(aphi) + yzero;
672 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
673 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
674 zpos = 0.;
675 atheta910 = (i-1) * 36. + offset1 + offset2 - gteta;
676 AliMatrix(idrotm[(i-1) * 13 + 1108], 90., atheta910, 90., atheta910 + 90., 0., 0.);
677 geant->Gsposp("SPIX", (i-1) * 13 + 9, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1108], "ONLY", dsup, 3);
678
679// --- Place part # 12-13 (see sketch)
680
681 offset1 = -27.;
682 dsup[0] = .01;
683 dsup[1] = TMath::Sqrt((xend[11] - xbeg[11]) * (xend[11] - xbeg[11]) + (yend[11] - ybeg[11]) * (yend[11] - ybeg[11])) / 20.;
684 dsup[2] = 25.;
685 xcc = (xx[11] + xx[12]) / 20.;
686 ycc = (yy[11] + yy[12]) / 20.;
687 xccc = (xbeg[11] + xend[11]) / 20.;
688 yccc = (ybeg[11] + yend[11]) / 20.;
689 if (xx[11] == xx[12]) {
690 offset2 = 0.;
691 } else {
692 r2 = xx[12] - xx[11];
693 offset2 = atg_(&yy[12] - yy[11], &r2) * raddeg - 90.;
694 }
695 aphi = (pphi + (i-1) * 36.) * degrad;
696 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
697 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
698 xpos1 = xccc * TMath::Cos(aphi) - yccc * TMath::Sin(aphi) + xzero;
699 ypos1 = xccc * TMath::Sin(aphi) + yccc * TMath::Cos(aphi) + yzero;
700 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
701 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
702 zpos = 0.;
703 atheta1213 = (i-1) * 36. + offset1 + offset2 - gteta;
704 AliMatrix(idrotm[(i-1) * 13 + 1111], 270., atheta1213, 90., atheta1213 + 270., 0., 0.);
705 geant->Gsposp("SPIX", (i-1) * 13 + 12, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1111], "ONLY", dsup, 3);
706
707// --- Place an element of layer #1
708
709 biga = (yy[12] - yy[11]) / (xx[12] - xx[11]);
710 bigb = (xx[12] * yy[11] - xx[11] * yy[12]) / (xx[12] - xx[11]) / 10.;
711 coeffa = biga * biga + 1.;
712 coeffb = biga * bigb - biga * ycc - xcc;
713 coeffc = xcc * xcc + ycc * ycc - ycc * 2. * bigb + bigb * bigb - .015792948900000003;
714 xcc1 = (-coeffb + TMath::Sqrt(coeffb * coeffb - coeffa * coeffc)) / coeffa;
715 ycc1 = biga * xcc1 + bigb;
716 biga1 = -1. / biga;
717 bigb1 = xcc1 / biga + ycc1;
718 coeffa = biga1 * biga1 + 1.;
719 coeffb = biga1 * bigb1 - biga1 * ycc1 - xcc1;
720 coeffc = xcc1 * xcc1 + ycc1 * ycc1 - ycc1 * 2. * bigb1 + bigb1 * bigb1 - (dsup[0] + dbox1[0]) * (dsup[0] + dbox1[0]);
721 xcc2 = (-coeffb + TMath::Sqrt(coeffb * coeffb - coeffa * coeffc)) / coeffa;
722 ycc2 = biga1 * xcc2 + bigb1;
723 xpos1 = xcc2 * TMath::Cos(aphi) - ycc2 * TMath::Sin(aphi) + xzero;
724 ypos1 = xcc2 * TMath::Sin(aphi) + ycc2 * TMath::Cos(aphi) + yzero;
725 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
726 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
727 zpos = 0.;
728 ++jbox1;
729 geant->Gspos("IPV1", jbox1, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1111], "ONLY");
730
731// --- Place part # 11-12 (see sketch)
732
733 offset1 = -27.;
734 dsup[0] = .01;
735 dsup[1] = TMath::Sqrt((xend[10] - xbeg[10]) * (xend[10] - xbeg[10]) + (yend[10] - ybeg[10]) * (yend[10] - ybeg[10])) / 20.;
736 dsup[2] = 25.;
737 xcc = (xx[10] + xx[11]) / 20.;
738 ycc = (yy[10] + yy[11]) / 20.;
739 xccc = (xbeg[10] + xend[10]) / 20.;
740 yccc = (ybeg[10] + yend[10]) / 20.;
741 if (xx[10] == xx[11]) {
742 offset2 = 0.;
743 } else {
744 r2 = xx[11] - xx[10];
745 offset2 = atg_(&yy[11] - yy[10], &r2) * raddeg - 90.;
746 }
747 aphi = (pphi + (i-1) * 36.) * degrad;
748 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
749 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
750 xpos1 = xccc * TMath::Cos(aphi) - yccc * TMath::Sin(aphi) + xzero;
751 ypos1 = xccc * TMath::Sin(aphi) + yccc * TMath::Cos(aphi) + yzero;
752 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
753 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
754 zpos = 0.;
755 atheta1112 = (i-1) * 36. + offset1 + offset2 - gteta;
756 AliMatrix(idrotm[(i-1) * 13 + 1110], 90., atheta1112, 90.,atheta1112 + 90., 0., 0.);
757 geant->Gsposp("SPIX", (i-1) * 13 + 11, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1110], "ONLY", dsup, 3);
758
759// --- Place part # 10-11 (see sketch)
760
761 offset1 = -27.;
762 dsup[0] = .01;
763 dsup[1] = TMath::Sqrt((xend[9] - xbeg[9]) * (xend[9] - xbeg[9]) + (yend[9] - ybeg[9]) * (yend[9] - ybeg[9])) / 20.;
764 dsup[2] = 25.;
765 xcc = (xx[9] + xx[10]) / 20.;
766 ycc = (yy[9] + yy[10]) / 20.;
767 xccc = (xbeg[9] + xend[9]) / 20.;
768 yccc = (ybeg[9] + yend[9]) / 20.;
769 if (xx[9] == xx[10]) {
770 offset2 = 0.;
771 } else {
772 r2 = xx[10] - xx[9];
773 offset2 = atg_(&yy[10] - yy[9], &r2) * raddeg - 90.;
774 }
775 aphi = (pphi + (i-1) * 36.) * degrad;
776 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
777 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
778 xpos1 = xccc * TMath::Cos(aphi) - yccc * TMath::Sin(aphi) + xzero;
779 ypos1 = xccc * TMath::Sin(aphi) + yccc * TMath::Cos(aphi) + yzero;
780 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
781 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
782 zpos = 0.;
783 atheta1011 = (i-1) * 36. + offset1 + offset2 - gteta;
784 AliMatrix(idrotm[(i-1) * 13 + 1109], 270., atheta1011, 90., atheta1011 + 270., 0., 0.);
785 geant->Gsposp("SPIX", (i-1) * 13 + 10, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1109], "ONLY", dsup, 3);
786
787// --- Place an element of layer #1
788
789 biga = (yy[10] - yy[9]) / (xx[10] - xx[9]);
790 bigb = (xx[10] * yy[9] - xx[9] * yy[10]) / (xx[10] - xx[9]) / 10.;
791 coeffa = biga * biga + 1.;
792 coeffb = biga * bigb - biga * ycc - xcc;
793 coeffc = xcc * xcc + ycc * ycc - ycc * 2. * bigb + bigb * bigb - 9.9999999999999995e-7;
794 xcc1 = (-coeffb - TMath::Sqrt(coeffb * coeffb - coeffa * coeffc)) / coeffa;
795 ycc1 = biga * xcc1 + bigb;
796 biga1 = -1. / biga;
797 bigb1 = xcc1 / biga + ycc1;
798 coeffa = biga1 * biga1 + 1.;
799 coeffb = biga1 * bigb1 - biga1 * ycc1 - xcc1;
800 coeffc = xcc1 * xcc1 + ycc1 * ycc1 - ycc1 * 2. * bigb1 + bigb1 * bigb1 - (dsup[0] + dbox1[0]) * (dsup[0] + dbox1[0]);
801 xcc2 = (-coeffb + TMath::Sqrt(coeffb * coeffb - coeffa * coeffc)) / coeffa;
802 ycc2 = biga1 * xcc2 + bigb1;
803 xpos1 = xcc2 * TMath::Cos(aphi) - ycc2 * TMath::Sin(aphi) + xzero;
804 ypos1 = xcc2 * TMath::Sin(aphi) + ycc2 * TMath::Cos(aphi) + yzero;
805 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
806 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
807 zpos = 0.;
808 ++jbox1;
809 geant->Gspos("IPV1", jbox1, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1109], "ONLY");
810
811// --- Place arc # 5 (between part 1-2 and part 2-3) (see sketch)
812
813 darc[0] = rarc[4] / 10. - .02;
814 darc[1] = rarc[4] / 10.;
815 darc[2] = 25.;
816 darc[3] = atheta12 + 90. - (i-1) * 36.;
817 darc[4] = atheta23 + 90. - (i-1) * 36.;
818 xcc = xarc[4] / 10.;
819 ycc = yarc[4] / 10.;
820 aphi = (pphi + (i-1) * 36.) * degrad;
821 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
822 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
823 xpos1 = xcc * TMath::Cos(aphi) - ycc * TMath::Sin(aphi) + xzero;
824 ypos1 = xcc * TMath::Sin(aphi) + ycc * TMath::Cos(aphi) + yzero;
825 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
826 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
827 zpos = 0.;
828 geant->Gsposp("SARC", (i-1) * 13 + 5, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1104], "ONLY", darc, 5);
829
830// --- Place arc # 6 (between part 2-3 and part 3-4) (see sketch)
831
832 darc[0] = rarc[5] / 10. - .02;
833 darc[1] = rarc[5] / 10.;
834 darc[2] = 25.;
835 darc[3] = atheta23 - (i-1) * 36.;
836 darc[4] = atheta34 - (i-1) * 36.;
837 xcc = xarc[5] / 10.;
838 ycc = yarc[5] / 10.;
839 aphi = (pphi + (i-1) * 36.) * degrad;
840 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
841 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
842 xpos1 = xcc * TMath::Cos(aphi) - ycc * TMath::Sin(aphi) + xzero;
843 ypos1 = xcc * TMath::Sin(aphi) + ycc * TMath::Cos(aphi) + yzero;
844 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
845 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
846 zpos = 0.;
847 geant->Gsposp("SARC", (i-1) * 13 + 6, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1105], "ONLY", darc, 5);
848
849// --- Place arc # 7 (between part 3-4 and part 4-5) (see sketch)
850
851 darc[0] = rarc[6] / 10. - .02;
852 darc[1] = rarc[6] / 10.;
853 darc[2] = 25.;
854 darc[3] = atheta45 + 90. - (i-1) * 36.;
855 darc[4] = atheta34 + 90. - (i-1) * 36.;
856 xcc = xarc[6] / 10.;
857 ycc = yarc[6] / 10.;
858 aphi = (pphi + (i-1) * 36.) * degrad;
859 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
860 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
861 xpos1 = xcc * TMath::Cos(aphi) - ycc * TMath::Sin(aphi) + xzero;
862 ypos1 = xcc * TMath::Sin(aphi) + ycc * TMath::Cos(aphi) + yzero;
863 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
864 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
865 zpos = 0.;
866 geant->Gsposp("SARC", (i-1) * 13 + 7, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1106], "ONLY", darc, 5);
867
868// --- Place arc # 8 (between part 4-5 and part 5-6) (see sketch)
869
870 darc[0] = rarc[7] / 10. - .02;
871 darc[1] = rarc[7] / 10.;
872 darc[2] = 25.;
873 darc[3] = atheta56 + 180. - (i-1) * 36.;
874 darc[4] = atheta45 + 180. - (i-1) * 36.;
875 xcc = xarc[7] / 10.;
876 ycc = yarc[7] / 10.;
877 aphi = (pphi + (i-1) * 36.) * degrad;
878 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
879 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
880 xpos1 = xcc * TMath::Cos(aphi) - ycc * TMath::Sin(aphi) + xzero;
881 ypos1 = xcc * TMath::Sin(aphi) + ycc * TMath::Cos(aphi) + yzero;
882 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
883 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
884 zpos = 0.;
885 geant->Gsposp("SARC", (i-1) * 13 + 8, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1107], "ONLY", darc, 5);
886
887// --- Place arc # 9 (between part 5-6 and part 6-7) (see sketch)
888
889 darc[0] = rarc[8] / 10. - .02;
890 darc[1] = rarc[8] / 10.;
891 darc[2] = 25.;
892 darc[3] = atheta56 - 90. - (i-1) * 36.;
893 darc[4] = atheta67 - 90. - (i-1) * 36.;
894 xcc = xarc[8] / 10.;
895 ycc = yarc[8] / 10.;
896 aphi = (pphi + (i-1) * 36.) * degrad;
897 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
898 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
899 xpos1 = xcc * TMath::Cos(aphi) - ycc * TMath::Sin(aphi) + xzero;
900 ypos1 = xcc * TMath::Sin(aphi) + ycc * TMath::Cos(aphi) + yzero;
901 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
902 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
903 zpos = 0.;
904 geant->Gsposp("SARC", (i-1) * 13 + 9, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1108], "ONLY", darc, 5);
905
906// --- Place arc # 10 (between part 6-7 and part 7-8) (see sketch)
907
908 darc[0] = rarc[9] / 10. - .02;
909 darc[1] = rarc[9] / 10.;
910 darc[2] = 25.;
911 darc[3] = atheta67 - (i-1) * 36.;
912 darc[4] = atheta78 - (i-1) * 36.;
913 xcc = xarc[9] / 10.;
914 ycc = yarc[9] / 10.;
915 aphi = (pphi + (i-1) * 36.) * degrad;
916 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
917 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
918 xpos1 = xcc * TMath::Cos(aphi) - ycc * TMath::Sin(aphi) + xzero;
919 ypos1 = xcc * TMath::Sin(aphi) + ycc * TMath::Cos(aphi) + yzero;
920 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
921 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
922 zpos = 0.;
923 geant->Gsposp("SARC", (i-1) * 13 + 10, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1109], "ONLY", darc, 5);
924
925// --- Place arc # 11 (between part 7-8 and part 8-9) (see sketch)
926
927 darc[0] = rarc[10] / 10. - .02;
928 darc[1] = rarc[10] / 10.;
929 darc[2] = 25.;
930 darc[3] = atheta89 - 90. - (i-1) * 36.;
931 darc[4] = atheta78 - 90. - (i-1) * 36.;
932 xcc = xarc[10] / 10.;
933 ycc = yarc[10] / 10.;
934 aphi = (pphi + (i-1) * 36.) * degrad;
935 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
936 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
937 xpos1 = xcc * TMath::Cos(aphi) - ycc * TMath::Sin(aphi) + xzero;
938 ypos1 = xcc * TMath::Sin(aphi) + ycc * TMath::Cos(aphi) + yzero;
939 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
940 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
941 zpos = 0.;
942 geant->Gsposp("SARC", (i-1) * 13 + 11, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1110], "ONLY", darc, 5);
943
944// --- Place arc # 12 (between part 8-9 and part 9-10) (see sketch)
945
946 darc[0] = rarc[11] / 10. - .02;
947 darc[1] = rarc[11] / 10.;
948 darc[2] = 25.;
949 darc[3] = atheta89 - (i-1) * 36.;
950 darc[4] = atheta910 - (i-1) * 36.;
951 xcc = xarc[11] / 10.;
952 ycc = yarc[11] / 10.;
953 aphi = (pphi + (i-1) * 36.) * degrad;
954 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
955 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
956 xpos1 = xcc * TMath::Cos(aphi) - ycc * TMath::Sin(aphi) + xzero;
957 ypos1 = xcc * TMath::Sin(aphi) + ycc * TMath::Cos(aphi) + yzero;
958 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
959 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
960 zpos = 0.;
961 geant->Gsposp("SARC", (i-1) * 13 + 12, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1111], "ONLY", darc, 5);
962
963// --- Place arc # 1 (between part 9-10 and part 10-11) (see sketch)
964
965 darc[0] = rarc[0] / 10. - .02;
966 darc[1] = rarc[0] / 10.;
967 darc[2] = 25.;
968 darc[3] = atheta1011 + 10. - (i-1) * 36.;
969 darc[4] = atheta910 - 120. - (i-1) * 36.;
970 xcc = xarc[0] / 10.;
971 ycc = yarc[0] / 10.;
972 aphi = (pphi + (i-1) * 36.) * degrad;
973 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
974 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
975 xpos1 = xcc * TMath::Cos(aphi) - ycc * TMath::Sin(aphi) + xzero;
976 ypos1 = xcc * TMath::Sin(aphi) + ycc * TMath::Cos(aphi) + yzero;
977 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
978 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
979 zpos = 0.;
980 geant->Gsposp("SARC", (i-1) * 13 + 1, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1100], "ONLY", darc, 5);
981
982// --- Place arc # 2 (between part 10-11 and part 11-12) (see sketch)
983
984 darc[0] = rarc[1] / 10. - .02;
985 darc[1] = rarc[1] / 10.;
986 darc[2] = 25.;
987 darc[3] = atheta1011 - (i-1) * 36.;
988 darc[4] = atheta1112 + 45. - (i-1) * 36.;
989 xcc = xarc[1] / 10.;
990 ycc = yarc[1] / 10.;
991 aphi = (pphi + (i-1) * 36.) * degrad;
992 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
993 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
994 xpos1 = xcc * TMath::Cos(aphi) - ycc * TMath::Sin(aphi) + xzero;
995 ypos1 = xcc * TMath::Sin(aphi) + ycc * TMath::Cos(aphi) + yzero;
996 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
997 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
998 zpos = 0.;
999 geant->Gsposp("SARC", (i-1) * 13 + 2, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1101], "ONLY", darc, 5);
1000
1001// --- Place arc # 3 (between part 11-12 and part 12-13) (see sketch)
1002
1003 darc[0] = rarc[2] / 10. - .02;
1004 darc[1] = rarc[2] / 10.;
1005 darc[2] = 25.;
1006 darc[3] = atheta1213 + 115. - (i-1) * 36.;
1007 darc[4] = atheta1112 + 115. - (i-1) * 36.;
1008 xcc = xarc[2] / 10.;
1009 ycc = yarc[2] / 10.;
1010 aphi = (pphi + (i-1) * 36.) * degrad;
1011 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
1012 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
1013 xpos1 = xcc * TMath::Cos(aphi) - ycc * TMath::Sin(aphi) + xzero;
1014 ypos1 = xcc * TMath::Sin(aphi) + ycc * TMath::Cos(aphi) + yzero;
1015 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
1016 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
1017 zpos = 0.;
1018 geant->Gsposp("SARC", (i-1) * 13 + 3, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1102], "ONLY", darc, 5);
1019
1020// --- Place arc # 4 (between part 12-13 and part 1-2) (see sketch)
1021
1022 darc[0] = rarc[3] / 10. - .02;
1023 darc[1] = rarc[3] / 10.;
1024 darc[2] = 25.;
1025 darc[3] = atheta1213 - (i-1) * 36.;
1026 darc[4] = atheta12 - (i-1) * 36.;
1027 xcc = xarc[3] / 10.;
1028 ycc = yarc[3] / 10.;
1029 aphi = (pphi + (i-1) * 36.) * degrad;
1030 xzero = rr * TMath::Cos((tteta + (i-1) * 36.) * degrad);
1031 yzero = rr * TMath::Sin((tteta + (i-1) * 36.) * degrad);
1032 xpos1 = xcc * TMath::Cos(aphi) - ycc * TMath::Sin(aphi) + xzero;
1033 ypos1 = xcc * TMath::Sin(aphi) + ycc * TMath::Cos(aphi) + yzero;
1034 xpos = xpos1 * TMath::Cos(gteta * degrad) + ypos1 * TMath::Sin(gteta *degrad);
1035 ypos = -xpos1 * TMath::Sin(gteta * degrad) + ypos1 * TMath::Cos(gteta * degrad);
1036 zpos = 0.;
1037 geant->Gsposp("SARC", (i-1) * 13 + 4, "IT12", xpos, ypos, zpos, idrotm[(i-1) * 13 + 1103], "ONLY", darc, 5);
1038
1039 }
1040
1041//************************************************************************
1042//* *
1043//* D R I F T S *
1044//* =========== *
1045//* *
1046//************************************************************************
1047
1048// --- Define a ghost volume containing the Silicon Drift Detectors
1049// (layer #3 and #4) and fill it with air or vacuum
1050
1051 xxm = 1.0444222222222224;
1052 dgh[0] = 0.;
1053 dgh[1] = 360.;
1054 dgh[2] = 4.;
1055 dgh[3] = -25. - 5.9900000000000002 / xxm -
1056 .0099999999999997868 / xxm - 17.990000000000002 / xxm;
1057 dgh[4] = 27.;
1058 dgh[5] = 27.;
1059 dgh[6] = -25. - 5.9900000000000002 / xxm - .0099999999999997868 / xxm;
1060 dgh[7] = 9.01;
1061 dgh[8] = 27.;
1062 dgh[9] = 5.9900000000000002 / xxm + 25. + .0099999999999997868 / xxm;
1063 dgh[10] = 9.01;
1064 dgh[11] = 27.;
1065 dgh[12] = 5.9900000000000002 / xxm + 25. +
1066 .0099999999999997868 / xxm + 17.990000000000002 / xxm;
1067 dgh[13] = 27.;
1068 dgh[14] = 27.;
1069 geant->Gsvolu("IT34", "PCON", idtmed[275], dgh, 15);
1070
1071// --- Place the ghost volume in its mother volume (ITSV) and make it
1072// invisible
1073
1074 geant->Gspos("IT34", 1, "ITSV", 0., 0., 0., 0, "ONLY");
1075 geant->Gsatt("IT34", "SEEN", 0);
1076
1077// --- Layer #3
1078
1079// GOTO 3456 ! skip ITS layer no. 3
1080
1081//--- Define a ghost volume containing a single ladder of layer #3 (with the
1082// smaller lenght of ribs) and fill it with air or vacuum
1083
1084 dbox1[0] = .61539999999999995;
1085 dbox1[1] = 3.85;
1086// the widest element is the sensitive element
1087 dbox1[2] = 35.850000000000001;
1088// 7.5 cm is the lenght
1089 geant->Gsvolu("IDV1", "BOX ", idtmed[228], dbox1, 3);
1090
1091// --- Make the ghost volume invisible
1092
1093 geant->Gsatt("IDV1", "SEEN", 0);
1094
1095// --- Define a volume containing the sensitive part of drifts
1096// (silicon, layer #3)
1097
1098 dits[0] = .0172;
1099// see material budget report by G. Feofilov
1100 dits[1] = 3.85;
1101 dits[2] = 4.35;
1102 geant->Gsvolu("ITS3", "BOX ", idtmed[224], dits, 3);
1103
1104//--- Define the part of the (smaller) rib between two sensitive parts made of
1105// carbon (layer #3)
1106
1107 dsup[0] = .5 - dits[0];
1108 dsup[1] = .01;
1109 dsup[2] = 35.850000000000001;
1110// 7.5 cm is the lenght
1111 geant->Gsvolu("IR11", "BOX ", idtmed[227], dsup, 3);
1112
1113//--- Define the first part of the (smaller) rib between two sensitive parts
1114// made of aluminum (layer #3)
1115
1116 dal1[0] = .5 - dits[0];
1117 dal1[1] = 4.8000000000000001e-4;
1118 dal1[2] = 35.850000000000001;
1119// 7.5 cm is the lenght
1120 geant->Gsvolu("IR12", "BOX ", idtmed[230], dal1, 3);
1121
1122//--- Define the part of the (smaller) rib between two sensitive parts made of
1123// kapton (layer #3)
1124
1125 dkap[0] = .5 - dits[0];
1126 dkap[1] = .01585;
1127 dkap[2] = 35.850000000000001;
1128// 7.5 cm is the lenght
1129 geant->Gsvolu("IR13", "BOX ", idtmed[236], dkap, 3);
1130
1131//--- Define the second part of the (smaller) rib between two sensitive parts
1132// made of aluminum (layer #3)
1133
1134 dal2[0] = .5 - dits[0];
1135 dal2[1] = .0013500000000000001;
1136 dal2[2] = 35.850000000000001;
1137// 7.5 cm is the lenght
1138 geant->Gsvolu("IR14", "BOX ", idtmed[230], dal2, 3);
1139
1140// --- Define the part of the (smaller) rib between two sensitive parts
1141// made of silicon (the electronics) (layer #3)
1142
1143 dchi[0] = .5 - dits[0];
1144 dchi[1] = .0035500000000000002;
1145 dchi[2] = 35.850000000000001;
1146// 7.5 cm is the lenght
1147 geant->Gsvolu("IR15", "BOX ", idtmed[225], dal2, 3);
1148
1149// --- Define the part of the (smaller) rib between two sensitive parts
1150// made of water (the cooler) (layer #3)
1151
1152 dwat[0] = .5 - dits[0];
1153 dwat[1] = .0046499999999999996;
1154 dwat[2] = 35.850000000000001;
1155// 7.5 cm is the lenght
1156 geant->Gsvolu("IR16", "BOX ", idtmed[231], dwat, 3);
1157
1158//--- Define the third part of the (smaller) rib between two sensitive parts
1159// made of aluminum (the cooling tubes) (layer #3)
1160
1161 dtub[0] = .5 - dits[0];
1162 dtub[1] = 6.7000000000000002e-4;
1163 dtub[2] = 35.850000000000001;
1164// 7.5 cm is the lenght
1165 geant->Gsvolu("IR17", "BOX ", idtmed[230], dtub, 3);
1166
1167// --- Define the part of the end-ladder stuff made of PCB (layer #3)
1168
1169 dpcb[0] = .03;
1170// twice the foreseen thickness
1171 dpcb[1] = 3.5;
1172 dpcb[2] = 7.5;
1173 geant->Gsvolu("IEL1", "BOX ", idtmed[233], dpcb, 3);
1174
1175// --- Define the part of the end-ladder stuff made of copper (layer #3)
1176
1177 dcop[0] = .0252;
1178// twice the foreseen thickness
1179 dcop[1] = 3.5;
1180 dcop[2] = 7.5;
1181 geant->Gsvolu("IEL2", "BOX ", idtmed[234], dcop, 3);
1182
1183// --- Define the part of the end-ladder stuff made of ceramics (layer #3)
1184 */
1185
1186 dcer[0] = .04;
1187// twice the foreseen thickness
1188 dcer[1] = 3.5;
1189 dcer[2] = 7.5;
1190 geant->Gsvolu("IEL3", "BOX ", idtmed[235], dcer, 3);
1191
1192// --- Define the part of the end-ladder stuff made of silicon (layer #3)
1193
1194 dsil[0] = .003;
1195// twice the foreseen thickness
1196 dsil[1] = 3.5;
1197 dsil[2] = 7.5;
1198 geant->Gsvolu("IEL4", "BOX ", idtmed[226], dsil, 3);
1199
1200//--- Place the sensitive part of the drifts (smaller ribs) into its mother
1201// (IDV1)
1202
1203 ypos = 0.;
1204 for (j = 1; j <= 5; ++j) {
1205// odd elements are up and even elements are down
1206 if (j == 1) {
1207 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0] * 2. - dits[0];
1208 zpos = 0. - dits[2] + 1. - dits[2] * 2. - .1 - dits[2];
1209 } else if (j == 2) {
1210 xpos = -dbox1[0] + dits[0];
1211 zpos = 0. - dits[2] + 1. - dits[2];
1212 } else if (j == 3) {
1213 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0] * 2. - dits[0];
1214 zpos = 0.;
1215 } else if (j == 4) {
1216 xpos = -dbox1[0] + dits[0];
1217 zpos = dits[2] + 0. - 1. + dits[2];
1218 } else if (j == 5) {
1219 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0] * 2. - dits[0];
1220 zpos = dits[2] + 0. - 1. + dits[2] * 2. + .1 + dits[2];
1221 }
1222 geant->Gspos("ITS3", j, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1223 }
1224
1225// --- Place the smaller ribs into their mother (IDV1)
1226
1227// --- Right ribs (just a matter of convention)
1228
1229 xpos = .5 - dbox1[0] + dits[0];
1230 zpos = 0.;
1231
1232// --- Carbon
1233
1234 ypos = 2.81;
1235 geant->Gspos("IR11", 1, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1236
1237// --- Aluminum #1
1238
1239 ypos = dsup[1] + 2.81 + dal1[1];
1240 geant->Gspos("IR12", 1, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1241
1242// --- Kapton
1243
1244 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1];
1245 geant->Gspos("IR13", 1, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1246
1247// --- Aluminum #2
1248
1249 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1];
1250 geant->Gspos("IR14", 1, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1251
1252// --- Silicon (chip)
1253
1254 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1] * 2. + dchi[1];
1255 geant->Gspos("IR15", 1, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1256
1257// --- Water
1258
1259 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1] * 2. + dchi[1] * 2. + dwat[1];
1260 geant->Gspos("IR16", 1, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1261
1262// --- Aluminum #3
1263
1264 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1] * 2. + dchi[1] * 2. + dwat[1] * 2.
1265 + dtub[1];
1266 geant->Gspos("IR17", 1, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1267
1268// --- Right ribs (just a matter of convention)
1269
1270// --- Carbon
1271
1272 ypos = -2.81;
1273 geant->Gspos("IR11", 2, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1274
1275// --- Aluminum #1
1276
1277 ypos = -(dsup[1] + 2.81 + dal1[1]);
1278 geant->Gspos("IR12", 2, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1279
1280// --- Kapton
1281
1282 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1]);
1283 geant->Gspos("IR13", 2, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1284
1285// --- Aluminum #2
1286
1287 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] *
1288 2. + dal2[1]);
1289 geant->Gspos("IR14", 2, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1290
1291// --- Silicon (chip)
1292
1293 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] *
1294 2. + dal2[1] * 2. + dchi[1]);
1295 geant->Gspos("IR15", 2, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1296
1297// --- Water
1298
1299 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] *
1300 2. + dal2[1] * 2. + dchi[1] * 2. + dwat[1]);
1301 geant->Gspos("IR16", 2, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1302
1303// --- Aluminum #3
1304
1305 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] *
1306 2. + dal2[1] * 2. + dchi[1] * 2. + dwat[1] * 2. + dtub[1]);
1307 geant->Gspos("IR17", 2, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1308
1309// --- Place the end-ladder stuff into its mother (IDV1)
1310
1311
1312// --- Negative-Z end-ladder
1313
1314 ypos = 0.;
1315 zpos = -28.350000000000001;
1316
1317// --- PCB
1318
1319 xpos = dbox1[0] - dpcb[0];
1320 geant->Gspos("IEL1", 1, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1321
1322// --- Copper
1323
1324 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0];
1325 geant->Gspos("IEL2", 1, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1326
1327// --- Ceramics
1328
1329 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0];
1330 geant->Gspos("IEL3", 1, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1331
1332// --- Silicon (bus)
1333
1334 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0];
1335 geant->Gspos("IEL4", 1, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1336
1337// --- Positive-Z end-ladder
1338
1339 ypos = 0.;
1340 zpos = 28.350000000000001;
1341
1342// --- PCB
1343
1344 xpos = dbox1[0] - dpcb[0];
1345 geant->Gspos("IEL1", 2, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1346
1347// --- Copper
1348
1349 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0];
1350 geant->Gspos("IEL2", 2, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1351
1352// --- Ceramics
1353
1354 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0];
1355 geant->Gspos("IEL3", 2, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1356
1357// --- Silicon (bus)
1358
1359 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0];
1360 geant->Gspos("IEL4", 2, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1361
1362//--- Define a ghost volume containing a single ladder of layer #3 (with the
1363// larger lenght of ribs) and fill it with air or vacuum
1364
1365 dbox2[0] = .76540000000000008;
1366 dbox2[1] = 3.85;
1367// the widest element is the sensitive element
1368 dbox2[2] = 35.850000000000001;
1369// 7.5 cm is the lenght
1370 geant->Gsvolu("IDV2", "BOX ", idtmed[228], dbox2, 3);
1371
1372// --- Make the ghost volume invisible
1373
1374 geant->Gsatt("IDV2", "SEEN", 0);
1375
1376//--- Define the part of the (larger) rib between two sensitive parts madeof
1377// carbon (layer #3)
1378
1379 dsup[0] = .65 - dits[0];
1380 dsup[1] = .01;
1381 dsup[2] = 35.850000000000001;
1382// 7.5 cm is the lenght
1383 geant->Gsvolu("IR21", "BOX ", idtmed[227], dsup, 3);
1384
1385//--- Define the first part of the (larger) rib between two sensitive parts
1386// made of aluminum (layer #3)
1387
1388 dal1[0] = .65 - dits[0];
1389 dal1[1] = 4.8000000000000001e-4;
1390 dal1[2] = 35.850000000000001;
1391// 7.5 cm is the lenght
1392 geant->Gsvolu("IR22", "BOX ", idtmed[230], dal1, 3);
1393
1394//--- Define the part of the (larger) rib between two sensitive parts madeof
1395// kapton (layer #3)
1396
1397 dkap[0] = .65 - dits[0];
1398 dkap[1] = .01585;
1399 dkap[2] = 35.850000000000001;
1400// 7.5 cm is the lenght
1401 geant->Gsvolu("IR23", "BOX ", idtmed[236], dkap, 3);
1402
1403//--- Define the second part of the (larger) rib between two sensitive parts
1404// made of aluminum (layer #3)
1405
1406 dal2[0] = .65 - dits[0];
1407 dal2[1] = .0013500000000000001;
1408 dal2[2] = 35.850000000000001;
1409// 7.5 cm is the lenght
1410 geant->Gsvolu("IR24", "BOX ", idtmed[230], dal2, 3);
1411
1412// --- Define the part of the (larger) rib between two sensitive parts
1413// made of silicon (the electronics) (layer #3)
1414
1415 dchi[0] = .65 - dits[0];
1416 dchi[1] = .0035500000000000002;
1417 dchi[2] = 35.850000000000001;
1418// 7.5 cm is the lenght
1419 geant->Gsvolu("IR25", "BOX ", idtmed[225], dal2, 3);
1420
1421// --- Define the part of the (larger) rib between two sensitive parts
1422// made of water (the cooler) (layer #3)
1423
1424 dwat[0] = .65 - dits[0];
1425 dwat[1] = .0046499999999999996;
1426 dwat[2] = 35.850000000000001;
1427// 7.5 cm is the lenght
1428 geant->Gsvolu("IR26", "BOX ", idtmed[231], dwat, 3);
1429
1430//--- Define the third part of the (larger) rib between two sensitive parts
1431// made of aluminum (the cooling tubes) (layer #3)
1432
1433 dtub[0] = .65 - dits[0];
1434 dtub[1] = 6.7000000000000002e-4;
1435 dtub[2] = 35.850000000000001;
1436// 7.5 cm is the lenght
1437 geant->Gsvolu("IR27", "BOX ", idtmed[230], dtub, 3);
1438
1439//--- Place the sensitive part of the drifts (smaller ribs) into its mother
1440// (IDV2)
1441
1442 ypos = 0.;
1443 for (j = 1; j <= 5; ++j) {
1444// odd element are up and even elements are down
1445 if (j == 1) {
1446 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0] * 2. - dits[0];
1447 zpos = 0. - dits[2] + 1. - dits[2] * 2. - .1 - dits[2];
1448 } else if (j == 2) {
1449 xpos = -dbox2[0] + dits[0];
1450 zpos = 0. - dits[2] + 1. - dits[2];
1451 } else if (j == 3) {
1452 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0] * 2. - dits[0];
1453 zpos = 0.;
1454 } else if (j == 4) {
1455 xpos = -dbox2[0] + dits[0];
1456 zpos = dits[2] + 0. - 1. + dits[2];
1457 } else if (j == 5) {
1458 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0] * 2. - dits[0];
1459 zpos = dits[2] + 0. - 1. + dits[2] * 2. + .1 + dits[2];
1460 }
1461 geant->Gspos("ITS3", j, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1462 }
1463
1464// --- Place the larger ribs into their mother (IDV2)
1465
1466
1467// --- Right ribs (just a matter of convention)
1468
1469 xpos = .65 - dbox2[0] + dits[0];
1470 zpos = 0.;
1471
1472// --- Carbon
1473
1474 ypos = 2.81;
1475 geant->Gspos("IR21", 1, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1476
1477// --- Aluminum #1
1478
1479 ypos = dsup[1] + 2.81 + dal1[1];
1480 geant->Gspos("IR22", 1, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1481
1482// --- Kapton
1483
1484 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1];
1485 geant->Gspos("IR23", 1, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1486
1487// --- Aluminum #2
1488
1489 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1];
1490 geant->Gspos("IR24", 1, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1491
1492// --- Silicon (chip)
1493
1494 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1] * 2. + dchi[1];
1495 geant->Gspos("IR25", 1, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1496
1497// --- Water
1498
1499 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1] * 2. + dchi[1] * 2. + dwat[1];
1500 geant->Gspos("IR26", 1, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1501
1502// --- Aluminum #3
1503
1504 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1] * 2. + dchi[1] * 2. + dwat[1] * 2.
1505 + dtub[1];
1506 geant->Gspos("IR27", 1, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1507
1508// --- Right ribs (just a matter of convention)
1509
1510// --- Carbon
1511
1512 ypos = -2.81;
1513 geant->Gspos("IR21", 2, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1514
1515// --- Aluminum #1
1516
1517 ypos = -(dsup[1] + 2.81 + dal1[1]);
1518 geant->Gspos("IR22", 2, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1519
1520// --- Kapton
1521
1522 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1]);
1523 geant->Gspos("IR23", 2, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1524
1525// --- Aluminum #2
1526
1527 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] *
1528 2. + dal2[1]);
1529 geant->Gspos("IR24", 2, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1530
1531// --- Silicon (chip)
1532
1533 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] *
1534 2. + dal2[1] * 2. + dchi[1]);
1535 geant->Gspos("IR25", 2, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1536
1537// --- Water
1538
1539 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] *
1540 2. + dal2[1] * 2. + dchi[1] * 2. + dwat[1]);
1541 geant->Gspos("IR26", 2, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1542
1543// --- Aluminum #3
1544
1545 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] *
1546 2. + dal2[1] * 2. + dchi[1] * 2. + dwat[1] * 2. + dtub[1]);
1547 geant->Gspos("IR27", 2, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1548
1549// --- Place the end-ladder stuff into its mother (IDV1)
1550
1551
1552// --- Negative-Z end-ladder
1553
1554 ypos = 0.;
1555 zpos = -28.350000000000001;
1556
1557// --- PCB
1558
1559 xpos = dbox2[0] - dpcb[0];
1560 geant->Gspos("IEL1", 3, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1561
1562// --- Copper
1563
1564 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0];
1565 geant->Gspos("IEL2", 3, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1566
1567// --- Ceramics
1568
1569 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0];
1570 geant->Gspos("IEL3", 3, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1571
1572// --- Silicon (bus)
1573
1574 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0];
1575 geant->Gspos("IEL4", 3, "IDV1", xpos, ypos, zpos, 0, "ONLY");
1576
1577// --- Positive-Z end-ladder
1578
1579 yos = 0.;
1580 zpos = 28.350000000000001;
1581
1582// --- PCB
1583
1584 xpos = dbox2[0] - dpcb[0];
1585 geant->Gspos("IEL1", 4, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1586
1587// --- Copper
1588
1589 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0];
1590 geant->Gspos("IEL2", 4, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1591
1592// --- Ceramics
1593
1594 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0];
1595 geant->Gspos("IEL3", 4, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1596
1597// --- Silicon (bus)
1598
1599 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0];
1600 geant->Gspos("IEL4", 4, "IDV2", xpos, ypos, zpos, 0, "ONLY");
1601
1602//--- Place the ghost volumes containing the drift ladders of layer #3 in their
1603// mother volume (IT34)
1604// Odd elements have large ribs and even elements have small ribs
1605
1606 for (i = 1; i <= 12; ++i) {
1607 atheta = (i-1) * 30.;
1608 AliMatrix(idrotm[i+1299], 90., atheta, 90., atheta + 90., 0.,0.);
1609 if (i % 2 == 0) {
1610 rzero = 14.;
1611 xpos = rzero * TMath::Cos((i-1) * twopi / 12.);
1612 ypos = rzero * TMath::Sin((i-1) * twopi / 12.);
1613 zpos = 0.;
1614 geant->Gspos("IDV1", i, "IT34", xpos, ypos, zpos, idrotm[i+1299], "ONLY");
1615 } else {
1616 rzero = 13.85;
1617 xpos = rzero * TMath::Cos((i-1) * twopi / 12.);
1618 ypos = rzero * TMath::Sin((i-1) * twopi / 12.);
1619 zpos = 0.;
1620 geant->Gspos("IDV2", i, "IT34", xpos, ypos, zpos, idrotm[i+1299], "ONLY");
1621 }
1622 }
1623
1624
1625// --- Layer #4
1626
1627// GOTO 4567 ! skip ITS layer no. 4
1628
1629//--- Define a ghost volume containing a single ladder of layer #4 (with the
1630// smaller lenght of ribs) and fill it with air or vacuum
1631
1632 dbox1[0] = .61539999999999995;
1633 dbox1[1] = 3.5;
1634// the widest element is the end-ladder stuff
1635 dbox1[2] = 43.449999999999996;
1636// 7.5 cm is the lenght
1637 geant->Gsvolu("IDV3", "BOX ", idtmed[228], dbox1, 3);
1638
1639// --- Make the ghost volume invisible
1640
1641 geant->Gsatt("IDV3", "SEEN", 0);
1642
1643// --- Define a volume containing the sensitive part of drifts
1644// (silicon, layer #4)
1645
1646 dits[0] = .0172;
1647// see material budget report by G. Feofilov
1648 dits[1] = 3.125;
1649 dits[2] = 4.35;
1650 geant->Gsvolu("ITS4", "BOX ", idtmed[224], dits, 3);
1651
1652//--- Define the part of the (smaller) rib between two sensitive parts made of
1653// carbon (layer #4)
1654
1655 dsup[0] = .5 - dits[0];
1656 dsup[1] = .01;
1657 dsup[2] = 43.449999999999996;
1658// 7.5 cm is the lengh
1659 geant->Gsvolu("IR31", "BOX ", idtmed[227], dsup, 3);
1660
1661//--- Define the first part of the (smaller) rib between two sensitive parts
1662// made of aluminum (layer #4)
1663
1664 dal1[0] = .5 - dits[0];
1665 dal1[1] = 4.8000000000000001e-4;
1666 dal1[2] = 43.449999999999996;
1667// 7.5 cm is the lengh
1668 geant->Gsvolu("IR32", "BOX ", idtmed[230], dal1, 3);
1669
1670//--- Define the part of the (smaller) rib between two sensitive parts made of
1671// kapton (layer #4)
1672
1673 dkap[0] = .5 - dits[0];
1674 dkap[1] = .01585;
1675 dkap[2] = 43.449999999999996;
1676// 7.5 cm is the lengh
1677 geant->Gsvolu("IR33", "BOX ", idtmed[236], dkap, 3);
1678
1679//--- Define the second part of the (smaller) rib between two sensitive parts
1680// made of aluminum (layer #4)
1681
1682 dal2[0] = .5 - dits[0];
1683 dal2[1] = .0013500000000000001;
1684 dal2[2] = 43.449999999999996;
1685// 7.5 cm is the lengh
1686 geant->Gsvolu("IR34", "BOX ", idtmed[230], dal2, 3);
1687
1688// --- Define the part of the (smaller) rib between two sensitive parts
1689// made of silicon (the electronics) (layer #4)
1690
1691 dchi[0] = .5 - dits[0];
1692 dchi[1] = .0035500000000000002;
1693 dchi[2] = 43.449999999999996;
1694// 7.5 cm is the lengh
1695 geant->Gsvolu("IR35", "BOX ", idtmed[225], dal2, 3);
1696
1697// --- Define the part of the (smaller) rib between two sensitive parts
1698// made of water (the cooler) (layer #4)
1699
1700 dwat[0] = .5 - dits[0];
1701 dwat[1] = .0046499999999999996;
1702 dwat[2] = 43.449999999999996;
1703// 7.5 cm is the lenght
1704 geant->Gsvolu("IR36", "BOX ", idtmed[231], dwat, 3);
1705
1706//--- Define the third part of the (smaller) rib between two sensitive parts
1707// made of aluminum (the cooling tubes) (layer #4)
1708
1709 dtub[0] = .5 - dits[0];
1710 dtub[1] = 6.7000000000000002e-4;
1711 dtub[2] = 43.449999999999996;
1712// 7.5 cm is the lengh
1713 geant->Gsvolu("IR37", "BOX ", idtmed[230], dtub, 3);
1714
1715// --- Define the part of the end-ladder stuff made of PCB (layer #4)
1716
1717 dpcb[0] = .03;
1718// twice the foreseen thickness
1719 dpcb[1] = 3.5;
1720 dpcb[2] = 7.5;
1721 geant->Gsvolu("IEL5", "BOX ", idtmed[233], dpcb, 3);
1722
1723// --- Define the part of the end-ladder stuff made of copper (layer #4)
1724
1725 dcop[0] = .0252;
1726// twice the foreseen thickness
1727 dcop[1] = 3.5;
1728 dcop[2] = 7.5;
1729 geant->Gsvolu("IEL6", "BOX ", idtmed[234], dcop, 3);
1730
1731// --- Define the part of the end-ladder stuff made of ceramics (layer #4)
1732 */
1733
1734 dcer[0] = .04;
1735// twice the foreseen thickness
1736 dcer[1] = 3.5;
1737 dcer[2] = 7.5;
1738 geant->Gsvolu("IEL7", "BOX ", idtmed[235], dcer, 3);
1739
1740// --- Define the part of the end-ladder stuff made of silicon (layer #4)
1741
1742 dsil[0] = .003;
1743// twice the foreseen thickness
1744 dsil[1] = 3.5;
1745 dsil[2] = 7.5;
1746 geant->Gsvolu("IEL8", "BOX ", idtmed[226], dsil, 3);
1747
1748//--- Place the sensitive part of the drifts (smaller ribs) into its mother
1749// (IDV3)
1750
1751 ypos = 0.;
1752 for (j = 1; j <= 7; ++j) {
1753// odd elements are down and even elements are up
1754 if (j == 1) {
1755 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0] * 2. - dits[0];
1756 zpos = 0. - dits[2] + .7 - dits[2] * 2. + 0. - dits[2] * 2. + 1.3 - dits[2];
1757 } else if (j == 2) {
1758 xpos = -dbox1[0] + dits[0];
1759 zpos = 0. - dits[2] + .7 - dits[2] * 2. + 0. - dits[2];
1760 } else if (j == 3) {
1761 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0] * 2. - dits[0];
1762 zpos = 0. - dits[2] + .7 - dits[2];
1763 } else if (j == 4) {
1764 xpos = -dbox1[0] + dits[0];
1765 zpos = 0.;
1766 } else if (j == 5) {
1767 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0] * 2. - dits[0];
1768 zpos = dits[2] + 0. - .7 + dits[2];
1769 } else if (j == 6) {
1770 xpos = -dbox1[0] + dits[0];
1771 zpos = dits[2] + 0. - .7 + dits[2] * 2. + 0. + dits[2];
1772 } else if (j == 7) {
1773 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0] * 2. - dits[0];
1774 zpos = dits[2] + 0. - .7 + dits[2] * 2. + 0. + dits[2] * 2. - 1.3 + dits[2];
1775 }
1776 geant->Gspos("ITS4", j, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1777 }
1778
1779// --- Place the smaller ribs into their mother (IDV3)
1780
1781// --- Right ribs (just a matter of convention)
1782
1783 xpos = .5 - dbox1[0] + dits[0];
1784 zpos = 0.;
1785
1786// --- Carbon
1787
1788 ypos = 2.81;
1789 geant->Gspos("IR31", 1, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1790
1791// --- Aluminum #1
1792
1793 ypos = dsup[1] + 2.81 + dal1[1];
1794 geant->Gspos("IR32", 1, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1795
1796// --- Kapton
1797
1798 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1];
1799 geant->Gspos("IR33", 1, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1800
1801// --- Aluminum #2
1802
1803 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1];
1804 geant->Gspos("IR34", 1, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1805
1806// --- Silicon (chip)
1807
1808 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1] * 2. + dchi[1];
1809 geant->Gspos("IR35", 1, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1810
1811// --- Water
1812
1813 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1] * 2. + dchi[1] * 2. + dwat[1];
1814 geant->Gspos("IR36", 1, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1815
1816// --- Aluminum #3
1817
1818 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1] * 2. + dchi[1] * 2. + dwat[1] * 2.
1819 + dtub[1];
1820 geant->Gspos("IR37", 1, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1821
1822// --- Right ribs (just a matter of convention)
1823
1824// --- Carbon
1825
1826 ypos = -2.81;
1827 geant->Gspos("IR31", 2, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1828
1829// --- Aluminum #1
1830
1831 ypos = -(dsup[1] + 2.81 + dal1[1]);
1832 geant->Gspos("IR32", 2, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1833
1834// --- Kapton
1835
1836 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1]);
1837 geant->Gspos("IR33", 2, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1838
1839// --- Aluminum #2
1840
1841 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1]);
1842 geant->Gspos("IR34", 2, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1843
1844// --- Silicon (chip)
1845
1846 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1] * 2. + dchi[1]);
1847 geant->Gspos("IR35", 2, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1848
1849// --- Water
1850
1851 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1] * 2. + dchi[1] * 2. + dwat[1]);
1852 geant->Gspos("IR36", 2, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1853
1854// --- Aluminum #3
1855
1856 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] *
1857 2. + dal2[1] * 2. + dchi[1] * 2. + dwat[1] * 2. + dtub[1]);
1858 geant->Gspos("IR37", 2, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1859
1860// --- Place the end-ladder stuff into its mother (IDV1)
1861
1862
1863// --- Negative-Z end-ladder
1864
1865 ypos = 0.;
1866 zpos = -35.949999999999996;
1867
1868// --- PCB
1869
1870 xpos = dbox1[0] - dpcb[0];
1871 geant->Gspos("IEL5", 1, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1872
1873// --- Copper
1874
1875 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0];
1876 geant->Gspos("IEL6", 1, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1877
1878// --- Ceramics
1879
1880 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0];
1881 geant->Gspos("IEL7", 1, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1882
1883// --- Silicon (bus)
1884
1885 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0];
1886 geant->Gspos("IEL8", 1, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1887
1888// --- Positive-Z end-ladder
1889
1890 ypos = 0.;
1891 zpos = 35.949999999999996;
1892
1893// --- PCB
1894
1895 xpos = dbox1[0] - dpcb[0];
1896 geant->Gspos("IEL5", 2, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1897
1898// --- Copper
1899
1900 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0];
1901 geant->Gspos("IEL6", 2, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1902
1903// --- Ceramics
1904
1905 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0];
1906 geant->Gspos("IEL7", 2, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1907
1908// --- Silicon (bus)
1909
1910 xpos = dbox1[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0];
1911 geant->Gspos("IEL8", 2, "IDV3", xpos, ypos, zpos, 0, "ONLY");
1912
1913//--- Define a ghost volume containing a single ladder of layer #4 (with the
1914// larger lenght of ribs) and fill it with air or vacuum
1915
1916 dbox2[0] = .76540000000000008;
1917 dbox2[1] = 3.5;
1918// the widest element is the end-ladder stuff
1919 dbox2[2] = 43.449999999999996;
1920// 7.5 cm is the lenght
1921 geant->Gsvolu("IDV4", "BOX ", idtmed[228], dbox2, 3);
1922
1923// --- Make the ghost volume invisible
1924
1925 geant->Gsatt("IDV4", "SEEN", 0);
1926
1927//--- Define the part of the (larger) rib between two sensitive parts madeof
1928// carbon (layer #4)
1929
1930 dsup[0] = .65 - dits[0];
1931 dsup[1] = .01;
1932 dsup[2] = 43.449999999999996;
1933// 7.5 cm is the lengh
1934 geant->Gsvolu("IR41", "BOX ", idtmed[227], dsup, 3);
1935
1936//--- Define the first part of the (larger) rib between two sensitive parts
1937// made of aluminum (layer #4)
1938
1939 dal1[0] = .65 - dits[0];
1940 dal1[1] = 4.8000000000000001e-4;
1941 dal1[2] = 43.449999999999996;
1942// 7.5 cm is the lengh
1943 geant->Gsvolu("IR42", "BOX ", idtmed[230], dal1, 3);
1944
1945//--- Define the part of the (larger) rib between two sensitive parts madeof
1946// kapton (layer #4)
1947
1948 dkap[0] = .65 - dits[0];
1949 dkap[1] = .01585;
1950 dkap[2] = 43.449999999999996;
1951// 7.5 cm is the lengh
1952 geant->Gsvolu("IR43", "BOX ", idtmed[236], dkap, 3);
1953
1954//--- Define the second part of the (larger) rib between two sensitive parts
1955// made of aluminum (layer #4)
1956
1957 dal2[0] = .65 - dits[0];
1958 dal2[1] = .0013500000000000001;
1959 dal2[2] = 43.449999999999996;
1960// 7.5 cm is the lengh
1961 geant->Gsvolu("IR44", "BOX ", idtmed[230], dal2, 3);
1962
1963// --- Define the part of the (larger) rib between two sensitive parts
1964// made of silicon (the electronics) (layer #4)
1965
1966 dchi[0] = .65 - dits[0];
1967 dchi[1] = .0035500000000000002;
1968 dchi[2] = 43.449999999999996;
1969// 7.5 cm is the lengh
1970 geant->Gsvolu("IR45", "BOX ", idtmed[225], dal2, 3);
1971
1972// --- Define the part of the (larger) rib between two sensitive parts
1973// made of water (the cooler) (layer #4)
1974
1975 dwat[0] = .65 - dits[0];
1976 dwat[1] = .0046499999999999996;
1977 dwat[2] = 43.449999999999996;
1978// 7.5 cm is the lengh
1979 geant->Gsvolu("IR46", "BOX ", idtmed[231], dwat, 3);
1980
1981//--- Define the third part of the (larger) rib between two sensitive parts
1982// made of aluminum (the cooling tubes) (layer #4)
1983
1984 dtub[0] = .65 - dits[0];
1985 dtub[1] = 6.7000000000000002e-4;
1986 dtub[2] = 43.449999999999996;
1987// 7.5 cm is the lengh
1988 geant->Gsvolu("IR47", "BOX ", idtmed[230], dtub, 3);
1989
1990//--- Place the sensitive part of the drifts (smaller ribs) into its mother
1991// (IDV4)
1992
1993 ypos = 0.;
1994 for (j = 1; j <= 7; ++j) {
1995// odd elements are down and even elements are up
1996 if (j == 1) {
1997 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0] * 2. - dits[0];
1998 zpos = 0. - dits[2] + .7 - dits[2] * 2. + 0. - dits[2] * 2. + 1.3 - dits[2];
1999 } else if (j == 2) {
2000 xpos = -dbox2[0] + dits[0];
2001 zpos = 0. - dits[2] + .7 - dits[2] * 2. + 0. - dits[2];
2002 } else if (j == 3) {
2003 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0] * 2. - dits[0];
2004 zpos = 0. - dits[2] + .7 - dits[2];
2005 } else if (j == 4) {
2006 xpos = -dbox2[0] + dits[0];
2007 zpos = 0.;
2008 } else if (j == 5) {
2009 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0] * 2. - dits[0];
2010 zpos = dits[2] + 0. - .7 + dits[2];
2011 } else if (j == 6) {
2012 xpos = -dbox2[0] + dits[0];
2013 zpos = dits[2] + 0. - .7 + dits[2] * 2. + 0. + dits[2];
2014 } else if (j == 7) {
2015 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0] * 2. - dits[0];
2016 zpos = dits[2] + 0. - .7 + dits[2] * 2. + 0. + dits[2] * 2. - 1.3 + dits[2];
2017 }
2018 geant->Gspos("ITS4", j, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2019 }
2020
2021// --- Place the larger ribs into their mother (IDV4)
2022
2023
2024// --- Right ribs (just a matter of convention)
2025
2026 xpos = .65 - dbox2[0] + dits[0];
2027 zpos = 0.;
2028
2029// --- Carbon
2030
2031 ypos = 2.81;
2032 geant->Gspos("IR41", 1, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2033
2034// --- Aluminum #1
2035
2036 ypos = dsup[1] + 2.81 + dal1[1];
2037 geant->Gspos("IR42", 1, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2038
2039// --- Kapton
2040
2041 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1];
2042 geant->Gspos("IR43", 1, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2043
2044// --- Aluminum #2
2045
2046 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1];
2047 geant->Gspos("IR44", 1, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2048
2049// --- Silicon (chip)
2050
2051 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1] * 2. + dchi[1];
2052 geant->Gspos("IR45", 1, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2053
2054// --- Water
2055
2056 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1] * 2. + dchi[1] * 2. + dwat[1];
2057 geant->Gspos("IR46", 1, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2058
2059// --- Aluminum #3
2060
2061 ypos = dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] * 2. + dal2[1] * 2. + dchi[1] * 2. + dwat[1] * 2.
2062 + dtub[1];
2063 geant->Gspos("IR47", 1, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2064
2065// --- Right ribs (just a matter of convention)
2066
2067// --- Carbon
2068
2069 ypos = -2.81;
2070 geant->Gspos("IR41", 2, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2071
2072// --- Aluminum #1
2073
2074 ypos = -(dsup[1] + 2.81 + dal1[1]);
2075 geant->Gspos("IR42", 2, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2076
2077// --- Kapton
2078
2079 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1]);
2080 geant->Gspos("IR43", 2, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2081
2082// --- Aluminum #2
2083
2084 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] *
2085 2. + dal2[1]);
2086 geant->Gspos("IR44", 2, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2087
2088// --- Silicon (chip)
2089
2090 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] *
2091 2. + dal2[1] * 2. + dchi[1]);
2092 geant->Gspos("IR45", 2, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2093
2094// --- Water
2095
2096 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] *
2097 2. + dal2[1] * 2. + dchi[1] * 2. + dwat[1]);
2098 geant->Gspos("IR46", 2, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2099
2100// --- Aluminum #3
2101
2102 ypos = -(dsup[1] + 2.81 + dal1[1] * 2. + dkap[1] *
2103 2. + dal2[1] * 2. + dchi[1] * 2. + dwat[1] * 2. + dtub[1]);
2104 geant->Gspos("IR47", 2, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2105
2106// --- Place the end-ladder stuff into its mother (IDV1)
2107
2108
2109// --- Negative-Z end-ladder
2110
2111 ypos = 0.;
2112 zpos = -35.949999999999996;
2113
2114// --- PCB
2115
2116 xpos = dbox2[0] - dpcb[0];
2117 geant->Gspos("IEL5", 3, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2118
2119// --- Copper
2120
2121 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0];
2122 geant->Gspos("IEL6", 3, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2123
2124// --- Ceramics
2125
2126 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0];
2127 geant->Gspos("IEL7", 3, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2128
2129// --- Silicon (bus)
2130
2131 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0];
2132 geant->Gspos("IEL8", 3, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2133
2134// --- Positive-Z end-ladder
2135
2136 yos = 0.;
2137 zpos = 35.949999999999996;
2138
2139// --- PCB
2140
2141 xpos = dbox2[0] - dpcb[0];
2142 geant->Gspos("IEL5", 4, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2143
2144// --- Copper
2145
2146 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0];
2147 geant->Gspos("IEL6", 4, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2148
2149// --- Ceramics
2150
2151 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0];
2152 geant->Gspos("IEL7", 4, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2153
2154// --- Silicon (bus)
2155
2156 xpos = dbox2[0] - dpcb[0] * 2. - dcop[0] * 2. - dcer[0] * 2. - dsil[0];
2157 geant->Gspos("IEL8", 4, "IDV4", xpos, ypos, zpos, 0, "ONLY");
2158
2159//--- Place the ghost volumes containing the drift ladders of layer #4 in their
2160// mother volume (IT34)
2161// Odd elements have large ribs and even elements have small ribs
2162
2163 for (i = 1; i <= 24; ++i) {
2164 atheta = (i-1) * 15.;
2165 AliMatrix(idrotm[i+1399], 90., atheta, 90., atheta + 90., 0.,0.);
2166 if (i % 2 == 0) {
2167 rzero = 23.5;
2168 xpos = rzero * TMath::Cos((i-1) * twopi / 24.);
2169 ypos = rzero * TMath::Sin((i-1) * twopi / 24.);
2170 zpos = 0.;
2171 geant->Gspos("IDV3", i, "IT34", xpos, ypos, zpos, idrotm[i+1399], "ONLY");
2172 } else {
2173 rzero = 23.399999999999999;
2174 xpos = rzero * TMath::Cos((i-1) * twopi / 24.);
2175 ypos = rzero * TMath::Sin((i-1) * twopi / 24.);
2176 zpos = 0.;
2177 geant->Gspos("IDV4", i, "IT34", xpos, ypos, zpos, idrotm[i+1399], "ONLY");
2178 }
2179 }
2180
2181
2182//************************************************************************
2183//* *
2184//* S T R I P S *
2185//* =========== *
2186//* *
2187//************************************************************************
2188
2189// --- Define SSD with the 35+39 lay-out
2190
2191 if (fMinorVersion == 0 || fMinorVersion == 1 || fMinorVersion ==2) {
2192
2193//--- Define ghost volume containing the Strip Detectors and fill it with air
2194// or vacuum
2195
2196 xxm = 1.0444222222222224;
2197 dgh[0] = 0.;
2198 dgh[1] = 360.;
2199 dgh[2] = 4.;
2200 dgh[3] = -25. - 5.9900000000000002 / xxm -
2201 .0099999999999997868 / xxm - 17.990000000000002 / xxm - 10. / xxm - 12.997999999999998 / xxm;
2202 dgh[4] = 49.998;
2203 dgh[5] = 49.998;
2204 dgh[6] = -25. - 5.9900000000000002 / xxm -
2205 .0099999999999997868 / xxm - 17.990000000000002 / xxm - 10. / xxm;
2206 dgh[7] = 37.;
2207 dgh[8] = 49.998;
2208 dgh[9] = 5.9900000000000002 / xxm + 25. +
2209 .0099999999999997868 / xxm + 17.990000000000002 / xxm + 10. / xxm;
2210 dgh[10] = 37.;
2211 dgh[11] = 49.998;
2212 dgh[12] = 5.9900000000000002 / xxm + 25. +
2213 .0099999999999997868 / xxm + 17.990000000000002 / xxm + 10. / xxm + 12.997999999999998 / xxm;
2214 dgh[13] = 49.998;
2215 dgh[14] = 49.998;
2216 geant->Gsvolu("IT56", "PCON", idtmed[275], dgh, 15);
2217 geant->Gspos("IT56", 1, "ITSV", 0., 0., 0., 0, "ONLY");
2218 geant->Gsatt("IT56", "SEEN", 0);
2219
2220// --- Layer #5
2221
2222// GOTO 5678 ! skip ITS layer no. 5
2223
2224//--- Define a ghost volume containing a single ladder of layer #5 andfill
2225// it with air or vacuum
2226
2227 dbox1[0] = .044999999999999998;
2228 dbox1[1] = 3.75;
2229 dbox1[2] = 45.109999999999999;
2230 geant->Gsvolu("ISV1", "BOX ", idtmed[253], dbox1, 3);
2231
2232// --- Make the ghost volume invisible
2233
2234 geant->Gsatt("ISV1", "SEEN", 0);
2235
2236// --- Define a ghost volume containing the electronics and cooling of
2237// a single ladder of layer #5 and fill it with air or vacuum
2238
2239 dsrv[0] = (TMath::Sqrt(3.) / 2. * 4.2 + .47 + .05) / 2.;
2240 dsrv[1] = 3.75;
2241 dsrv[2] = 45.109999999999999;
2242 geant->Gsvolu("SSV1", "BOX ", idtmed[253], dsrv, 3);
2243
2244// --- Make the ghost volume invisible
2245
2246 geant->Gsatt("SSV1", "SEEN", 0);
2247
2248// --- Define a ghost volume containing the end-ladder stuff of
2249// a single ladder of layer #5 and fill it with air or vacuum
2250
2251 dela[0] = 2.;
2252 dela[1] = 3.5;
2253 dela[2] = 4.;
2254 geant->Gsvolu("ELL5", "BOX ", idtmed[253], dela, 3);
2255
2256// --- Make the ghost volume invisible
2257
2258 geant->Gsatt("ELL5", "SEEN", 0);
2259
2260// --- Define a volume containing the sensitive part of the strips
2261// (silicon, layer #5)
2262
2263 dits[0] = .015;
2264 dits[1] = 3.75;
2265 dits[2] = 2.1;
2266 geant->Gsvolu("ITS5", "BOX ", idtmed[249], dits, 3);
2267
2268// --- Define a volume containing the electronics of the strips
2269// (silicon, layer #5)
2270
2271 dchi[0] = .02;
2272 dchi[1] = 3.4;
2273 dchi[2] = .525;
2274 geant->Gsvolu("SCH5", "BOX ", idtmed[250], dchi, 3);
2275
2276// --- Define the cooling tubes (aluminum, layer #5)
2277
2278 dtub[0] = .09;
2279 dtub[1] = dtub[0] + .01;
2280 dtub[2] = 45.109999999999999;
2281 geant->Gsvolu("STB5", "TUBE", idtmed[255], dtub, 3);
2282
2283// --- Define the cooling fluid (water or freon, layer #5)
2284
2285 dwat[0] = 0.;
2286 dwat[1] = .09;
2287 dwat[2] = 45.109999999999999;
2288 geant->Gsvolu("SWT5", "TUBE", idtmed[256], dwat, 3);
2289// CALL GSVOLU('SWT5','TUBE',IDTMED(258),DWAT,3,IOUT) ! freon
2290
2291//--- Define the (triangular) element of the heat bridge (carbon, layer #5)
2292
2293// water
2294 dfra[0] = 120.;
2295 dfra[1] = 360.;
2296 dfra[2] = 3.;
2297 dfra[3] = 2.;
2298 dfra[4] = -.015;
2299 dfra[5] = TMath::Sqrt(3.) * 4.2 / 6.;
2300 dfra[6] = dfra[5] + .03;
2301 dfra[7] = .015;
2302 dfra[8] = dfra[5];
2303 dfra[9] = dfra[6];
2304 geant->Gsvolu("SFR5", "PGON", idtmed[252], dfra, 10);
2305
2306// --- Define the element connecting the triangles of the heat bridge
2307// (carbon, layer #5)
2308
2309 dcei[0] = 0.;
2310 dcei[1] = .03;
2311 dcei[2] = 45.109999999999999;
2312 geant->Gsvolu("SCE5", "TUBE", idtmed[252], dcei, 3);
2313
2314// --- Define the part of the end-ladder stuff made of plastic (G10FR4)
2315// (layer #5)
2316
2317 dpla[0] = .089285714285714288;
2318 dpla[1] = 3.5;
2319 dpla[2] = 4.;
2320 geant->Gsvolu("EPL5", "BOX ", idtmed[262], dpla, 3);
2321
2322// --- Define the part of the end-ladder stuff made of copper (layer #5)
2323
2324 dcop[0] = .017857142857142856;
2325 dcop[1] = 3.5;
2326 dcop[2] = 4.;
2327 geant->Gsvolu("ECU5", "BOX ", idtmed[259], dcop, 3);
2328
2329// --- Define the part of the end-ladder stuff made of epoxy (layer #5)
2330
2331 depx[0] = .26785714285714285;
2332 depx[1] = 3.5;
2333 depx[2] = 4.;
2334 geant->Gsvolu("EPX5", "BOX ", idtmed[262], depx, 3);
2335
2336// --- Define the part of the end-ladder stuff made of silicon (bus)
2337// (layer #5)
2338
2339 dsil[0] = .17857142857142858;
2340 dsil[1] = 3.5;
2341 dsil[2] = 4.;
2342 geant->Gsvolu("ESI5", "BOX ", idtmed[251], dsil, 3);
2343
2344// --- Place the end-ladder stuff into its mother (ELL5)
2345
2346 sep = (4. - (dpla[0] + dcop[0] + depx[0] + dsil[0]) * 2.) / 3.;
2347 ypos = 0.;
2348 zpos = 0.;
2349
2350// --- Plastic
2351
2352 xpos = -dela[0] + dpla[0];
2353 geant->Gspos("EPL5", 1, "ELL5", xpos, ypos, zpos, 0, "ONLY");
2354
2355// --- Copper
2356
2357 xpos = -dela[0] + dpla[0] * 2. + sep + dcop[0];
2358 geant->Gspos("ECU5", 1, "ELL5", xpos, ypos, zpos, 0, "ONLY");
2359
2360// --- Epoxy
2361
2362 xpos = -dela[0] + dpla[0] * 2. + sep + dcop[0] * 2. + sep + depx[0];
2363 geant->Gspos("EPX5", 1, "ELL5", xpos, ypos, zpos, 0, "ONLY");
2364
2365// --- Silicon (bus)
2366
2367 xpos = -dela[0] + dpla[0] * 2. + sep + dcop[0] * 2. + sep + depx[0] * 2. + sep + dsil[0];
2368 geant->Gspos("ESI5", 1, "ELL5", xpos, ypos, zpos, 0, "ONLY");
2369
2370// --- Place the sensitive part of the strips into its mother (ISV1)
2371
2372 ypos = 0.;
2373 for (j = 1; j <= 23; ++j) {
2374 if (j % 2 == 0) xpos = dbox1[0] - dits[0];
2375 else xpos = -dbox1[0] + dits[0];
2376 zpos = ((j - 1) - 11.) * 3.91;
2377 geant->Gspos("ITS5", j, "ISV1", xpos, ypos, zpos, 0, "ONLY");
2378 }
2379
2380// --- Place the electronics of the strips into its mother (SSV1)
2381
2382 ypos = 0.;
2383 for (j = 1; j <= 23; ++j) {
2384 if (j % 2 == 0) xpos = -dsrv[0] + .28;
2385 else xpos = -dsrv[0] + .28 - dits[0] * 2. - .03;
2386 zpos = ((j - 1) - 11.) * 3.91 + .85;
2387 geant->Gspos("SCH5", j, "SSV1", xpos, ypos, zpos, 0, "ONLY");
2388 }
2389
2390//--- Place the cooling tubes and the cooling fluid into their mother (SSV1)
2391
2392 xpos = -dsrv[0] + .41;
2393 zpos = 0.;
2394
2395// --- Left tube (just a matter of convention)
2396
2397 ypos = -2.3500000000000001;
2398 geant->Gspos("STB5", 1, "SSV1", xpos, ypos, zpos, 0, "ONLY");
2399 geant->Gspos("SWT5", 1, "SSV1", xpos, ypos, zpos, 0, "ONLY");
2400
2401// --- Right tube (just a matter of convention)
2402
2403 ypos = 2.3500000000000001;
2404 geant->Gspos("STB5", 2, "SSV1", xpos, ypos, zpos, 0, "ONLY");
2405 geant->Gspos("SWT5", 2, "SSV1", xpos, ypos, zpos, 0, "ONLY");
2406
2407// --- Place the heat bridge elements into their mother (SSV1)
2408
2409 xpos = -dsrv[0] + .47 + TMath::Sqrt(3.) / 6. * 4.2;
2410 ypos = 0.;
2411 for (j = 1; j <= 24; ++j) {
2412 zpos = ((j - 1) - 11.) * 3.91 - 2.1000000000000001;
2413 geant->Gspos("SFR5", j, "SSV1", xpos, ypos, zpos, 0, "ONLY");
2414 }
2415
2416// --- Place the elements connecting the triangles of the heat bridge
2417// into their mother (SSV1)
2418
2419 zpos = 0.;
2420
2421// --- Left element (just a matter of convention)
2422
2423 xpos = -dsrv[0] + .47;
2424 ypos = -2.1150000000000002;
2425 geant->Gspos("SCE5", 1, "SSV1", xpos, ypos, zpos, 0, "ONLY");
2426
2427// --- Right element
2428
2429 xpos = -dsrv[0] + .47;
2430 ypos = 2.1150000000000002;
2431 geant->Gspos("SCE5", 2, "SSV1", xpos, ypos, zpos, 0, "ONLY");
2432
2433// --- Top element
2434
2435 xpos = -dsrv[0] + .47 + TMath::Sqrt(3.) / 2. * 4.2 + .015;
2436 ypos = 0.;
2437 geant->Gspos("SCE5", 3, "SSV1", xpos, ypos, zpos, 0, "ONLY");
2438
2439// --- Place the ghost volumes containing the strip ladders (ISV1),
2440// electronics/cooling (SSV1) and end-ladder stuff (ELL5) of layer #5 in
2441// their mother volume (IT56)
2442
2443 offset1 = atg_(.9, 40.);
2444 offset2 = 5.2;
2445 rzero = dbox1[0] + 40.;
2446 runo = dbox1[0] * 2. + 40. + dsrv[0];
2447 rtwo = dbox1[0] * 2. + 40. + dela[0];
2448 for (i = 1; i <= 35; ++i) {
2449 atheta = (i-1) * twopi * raddeg / 35. + offset2;
2450 AliMatrix(idrotm[i+1499], 90., atheta, 90., atheta + 90., 0., 0.);
2451
2452// --- Strip ladders
2453
2454 xpos = rzero * TMath::Cos((i-1) * twopi / 35. + offset1);
2455 ypos = rzero * TMath::Sin((i-1) * twopi / 35. + offset1);
2456 zpos = 0.;
2457 geant->Gspos("ISV1", i, "IT56", xpos, ypos, zpos, idrotm[i+1499], "ONLY");
2458
2459// --- Electronics/cooling
2460
2461 xpos = runo * TMath::Cos((i-1) * twopi / 35. + offset1);
2462 ypos = runo * TMath::Sin((i-1) * twopi / 35. + offset1);
2463 zpos = 0.;
2464 geant->Gspos("SSV1", i, "IT56", xpos, ypos, zpos, idrotm[i+1499], "ONLY");
2465
2466// --- End-ladders (nagative-Z and positive-Z)
2467
2468 xpos = rtwo * TMath::Cos((i-1) * twopi / 35. + offset1);
2469 ypos = rtwo * TMath::Sin((i-1) * twopi / 35. + offset1);
2470 zpos = -(dbox1[2] + dela[2] + 6.);
2471 geant->Gspos("ELL5", i, "IT56", xpos, ypos, zpos, idrotm[i+1499], "ONLY");
2472 zpos = dbox1[2] + dela[2] + 6.;
2473 geant->Gspos("ELL5", i + 35, "IT56", xpos, ypos, zpos, idrotm[i+1499], "ONLY");
2474 }
2475
2476
2477// --- Layer #6
2478
2479// GOTO 5778 ! skip ITS layer no. 6
2480
2481//--- Define a ghost volume containing a single ladder of layer #6 andfill
2482// it with air or vacuum
2483
2484 dbox2[0] = .044999999999999998;
2485 dbox2[1] = 3.75;
2486 dbox2[2] = 50.975000000000001;
2487 geant->Gsvolu("ISV2", "BOX ", idtmed[253], dbox2, 3);
2488
2489// --- Make the ghost volume invisible
2490
2491 geant->Gsatt("ISV2", "SEEN", 0);
2492
2493// --- Define a ghost volume containing the electronics and cooling of
2494// a single ladder of layer #6 and fill it with air or vacuum
2495
2496 dsrv[0] = (TMath::Sqrt(3.) / 2. * 4.2 + .47 + .05) / 2.;
2497 dsrv[1] = 3.75;
2498 dsrv[2] = 50.975000000000001;
2499 geant->Gsvolu("SSV2", "BOX ", idtmed[253], dsrv, 3);
2500
2501// --- Make the ghost volume invisible
2502
2503 geant->Gsatt("SSV2", "SEEN", 0);
2504
2505// --- Define a ghost volume containing the end-ladder stuff of
2506// a single ladder of layer #6 and fill it with air or vacuum
2507
2508 dela[0] = 2.;
2509 dela[1] = 3.5;
2510 dela[2] = 4.;
2511 geant->Gsvolu("ELL6", "BOX ", idtmed[253], dela, 3);
2512
2513// --- Make the ghost volume invisible
2514
2515 geant->Gsatt("ELL6", "SEEN", 0);
2516
2517// --- Define a volume containing the sensitive part of the strips
2518// (silicon, layer #6)
2519
2520 dits[0] = .015;
2521 dits[1] = 3.75;
2522 dits[2] = 2.1;
2523 geant->Gsvolu("ITS6", "BOX ", idtmed[249], dits, 3);
2524
2525// --- Define a volume containing the electronics of the strips
2526// (silicon, layer #6)
2527
2528 dchi[0] = .02;
2529 dchi[1] = 3.4;
2530 dchi[2] = .525;
2531 geant->Gsvolu("SCH6", "BOX ", idtmed[250], dchi, 3);
2532
2533// --- Define the cooling tubes (aluminum, layer #6)
2534
2535 dtub[0] = .09;
2536 dtub[1] = dtub[0] + .01;
2537 dtub[2] = 50.975000000000001;
2538 geant->Gsvolu("STB6", "TUBE", idtmed[255], dtub, 3);
2539
2540// --- Define the cooling fluid (water or freon, layer #6)
2541
2542 dwat[0] = 0.;
2543 dwat[1] = .09;
2544 dwat[2] = 50.975000000000001;
2545 geant->Gsvolu("SWT6", "TUBE", idtmed[256], dwat, 3);
2546// CALL GSVOLU('SWT6','TUBE',IDTMED(258),DWAT,3,IOUT) ! freon
2547
2548//--- Define the (triangular) element of the heat bridge (carbon, layer #6)
2549
2550// water
2551 dfra[0] = 120.;
2552 dfra[1] = 360.;
2553 dfra[2] = 3.;
2554 dfra[3] = 2.;
2555 dfra[4] = -.015;
2556 dfra[5] = TMath::Sqrt(3.) * 4.2 / 6.;
2557 dfra[6] = dfra[5] + .03;
2558 dfra[7] = .015;
2559 dfra[8] = dfra[5];
2560 dfra[9] = dfra[6];
2561 geant->Gsvolu("SFR6", "PGON", idtmed[252], dfra, 10);
2562
2563// --- Define the element connecting the triangles of the heat bridge
2564// (carbon, layer #6)
2565
2566 dcei[0] = 0.;
2567 dcei[1] = .03;
2568 dcei[2] = 50.975000000000001;
2569 geant->Gsvolu("SCE6", "TUBE", idtmed[252], dcei, 3);
2570
2571// --- Define the part of the end-ladder stuff made of plastic (G10FR4)
2572// (layer #6)
2573
2574 dpla[0] = .089285714285714288;
2575 dpla[1] = 3.5;
2576 dpla[2] = 4.;
2577 geant->Gsvolu("EPL6", "BOX ", idtmed[262], dpla, 3);
2578
2579// --- Define the part of the end-ladder stuff made of copper (layer #6)
2580
2581 dcop[0] = .017857142857142856;
2582 dcop[1] = 3.5;
2583 dcop[2] = 4.;
2584 geant->Gsvolu("ECU6", "BOX ", idtmed[259], dcop, 3);
2585
2586// --- Define the part of the end-ladder stuff made of epoxy (layer #6)
2587
2588 depx[0] = .26785714285714285;
2589 depx[1] = 3.5;
2590 depx[2] = 4.;
2591 geant->Gsvolu("EPX6", "BOX ", idtmed[262], depx, 3);
2592
2593// --- Define the part of the end-ladder stuff made of silicon (bus)
2594// (layer #6)
2595
2596 dsil[0] = .17857142857142858;
2597 dsil[1] = 3.5;
2598 dsil[2] = 4.;
2599 geant->Gsvolu("ESI6", "BOX ", idtmed[251], dsil, 3);
2600
2601// --- Place the end-ladder stuff into its mother (ELL5)
2602
2603 sep = (4. - (dpla[0] + dcop[0] + depx[0] + dsil[0]) * 2.) / 3.;
2604 ypos = 0.;
2605 zpos = 0.;
2606
2607// --- Plastic
2608
2609 xpos = -dela[0] + dpla[0];
2610 geant->Gspos("EPL6", 1, "ELL6", xpos, ypos, zpos, 0, "ONLY");
2611
2612// --- Copper
2613
2614 xpos = -dela[0] + dpla[0] * 2. + sep + dcop[0];
2615 geant->Gspos("ECU6", 1, "ELL6", xpos, ypos, zpos, 0, "ONLY");
2616
2617// --- Epoxy
2618
2619 xpos = -dela[0] + dpla[0] * 2. + sep + dcop[0] * 2. + sep + depx[0];
2620 geant->Gspos("EPX6", 1, "ELL6", xpos, ypos, zpos, 0, "ONLY");
2621
2622// --- Silicon (bus)
2623
2624 xpos = -dela[0] + dpla[0] * 2. + sep + dcop[0] * 2. + sep + depx[0] * 2. + sep + dsil[0];
2625 geant->Gspos("ESI6", 1, "ELL6", xpos, ypos, zpos, 0, "ONLY");
2626
2627// --- Place the sensitive part of the strips into its mother (ISV2)
2628
2629 ypos = 0.;
2630 for (j = 1; j <= 26; ++j) {
2631 if (j % 2 == 0) xpos = dbox2[0] - dits[0];
2632 else xpos = -dbox2[0] + dits[0];
2633 zpos = ((j - 1) - 12.) * 3.91 - 1.96;
2634 geant->Gspos("ITS6", j, "ISV2", xpos, ypos, zpos, 0, "ONLY");
2635 }
2636
2637// --- Place the electronics of the strips into its mother (SSV2)
2638
2639 ypos = 0.;
2640 for (j = 1; j <= 26; ++j) {
2641 if (j % 2 == 0) xpos = -dsrv[0] + .28;
2642 else xpos = -dsrv[0] + .28 - dits[0] * 2. - .03;
2643 zpos = ((j - 1) - 12.) * 3.91 - 1.96 + .85;
2644 geant->Gspos("SCH5", j, "SSV2", xpos, ypos, zpos, 0, "ONLY");
2645 }
2646
2647//--- Place the cooling tubes and the cooling fluid into their mother (SSV1)
2648
2649 xpos = -dsrv[0] + .41;
2650 zpos = 0.;
2651
2652// --- Left tube (just a matter of convention)
2653
2654 ypos = -2.3500000000000001;
2655 geant->Gspos("STB6", 1, "SSV2", xpos, ypos, zpos, 0, "ONLY");
2656 geant->Gspos("SWT6", 1, "SSV2", xpos, ypos, zpos, 0, "ONLY");
2657
2658// --- Right tube (just a matter of convention)
2659
2660 ypos = 2.3500000000000001;
2661 geant->Gspos("STB6", 2, "SSV2", xpos, ypos, zpos, 0, "ONLY");
2662 geant->Gspos("SWT6", 2, "SSV2", xpos, ypos, zpos, 0, "ONLY");
2663
2664// --- Place the heat bridge elements into their mother (SSV2)
2665
2666 xpos = -dsrv[0] + .47 + TMath::Sqrt(3.) / 6. * 4.2;
2667 ypos = 0.;
2668 for (j = 1; j <= 27; ++j) {
2669 zpos = ((j - 1) - 12.) * 3.91 - 1.96 - 2.1000000000000001;
2670 geant->Gspos("SFR6", j, "SSV2", xpos, ypos, zpos, 0, "ONLY");
2671 }
2672
2673// --- Place the elements connecting the triangles of the heat bridge
2674// into their mother (SSV2)
2675
2676 zpos = 0.;
2677
2678// --- Left element (just a matter of convention)
2679
2680 xpos = -dsrv[0] + .47;
2681 ypos = -2.1150000000000002;
2682 geant->Gspos("SCE6", 1, "SSV2", xpos, ypos, zpos, 0, "ONLY");
2683
2684// --- Right element
2685
2686 xpos = -dsrv[0] + .47;
2687 ypos = 2.1150000000000002;
2688 geant->Gspos("SCE6", 2, "SSV2", xpos, ypos, zpos, 0, "ONLY");
2689
2690// --- Top element
2691
2692 xpos = -dsrv[0] + .47 + TMath::Sqrt(3.) / 2. * 4.2 + .015;
2693 ypos = 0.;
2694 geant->Gspos("SCE6", 3, "SSV2", xpos, ypos, zpos, 0, "ONLY");
2695
2696// --- Place the ghost volumes containing the strip ladders (ISV2),
2697// electronics/cooling (SSV2) and end-ladder stuff (ELL6) of layer #6 in
2698// their mother volume (IT56)
2699
2700 offset1 = atg_(1., 45.);
2701 offset2 = 5.2;
2702 rzero = dbox2[0] + 45.;
2703 runo = dbox2[0] * 2. + 45. + dsrv[0];
2704 rtwo = dbox2[0] * 2. + 45. + dela[0];
2705 for (i = 1; i <= 39; ++i) {
2706 atheta = (i-1) * twopi * raddeg / 39. + offset2;
2707 AliMatrix(idrotm[i+1599], 90., atheta, 90., atheta + 90., 0., 0.);
2708
2709// --- Strip ladders
2710
2711 xpos = rzero * TMath::Cos((i-1) * twopi / 39. + offset1);
2712 ypos = rzero * TMath::Sin((i-1) * twopi / 39. + offset1);
2713 zpos = 0.;
2714 geant->Gspos("ISV2", i, "IT56", xpos, ypos, zpos, idrotm[i+1599], "ONLY");
2715
2716// --- Electronics/cooling
2717
2718 xpos = runo * TMath::Cos((i-1) * twopi / 39. + offset1);
2719 ypos = runo * TMath::Sin((i-1) * twopi / 39. + offset1);
2720 zpos = 0.;
2721 geant->Gspos("SSV2", i, "IT56", xpos, ypos, zpos, idrotm[i+1599], "ONLY");
2722
2723// --- End-ladders (nagative-Z and positive-Z)
2724
2725 xpos = rtwo * TMath::Cos((i-1) * twopi / 39. + offset1);
2726 ypos = rtwo * TMath::Sin((i-1) * twopi / 39. + offset1);
2727 zpos = -(dbox2[2] + dela[2] + 6.);
2728 geant->Gspos("ELL6", i, "IT56", xpos, ypos, zpos, idrotm[i+1599], "ONLY");
2729 zpos = dbox2[2] + dela[2] + 6.;
2730 geant->Gspos("ELL6", i + 39, "IT56", xpos, ypos, zpos, idrotm[i+1599], "ONLY");
2731 }
2732
2733
2734 }
2735
2736// --- Define SSD with the 32+36 lay-out
2737
2738 if (fMinorVersion == 3 || fMinorVersion == 4 || fMinorVersion == 5) {
2739
2740//--- Define ghost volume containing the Strip Detectors and fill it with air
2741// or vacuum
2742
2743 xxm = 1.0444222222222224;
2744 dgh[0] = 0.;
2745 dgh[1] = 360.;
2746 dgh[2] = 4.;
2747 dgh[3] = -25. - 5.9900000000000002 / xxm -
2748 .0099999999999997868 / xxm - 17.990000000000002 / xxm - 9. / xxm - 13.997999999999998 / xxm;
2749 dgh[4] = 49.998;
2750 dgh[5] = 49.998;
2751 dgh[6] = -25. - 5.9900000000000002 / xxm -
2752 .0099999999999997868 / xxm - 17.990000000000002 / xxm - 9. / xxm;
2753 dgh[7] = 36.;
2754 dgh[8] = 49.998;
2755 dgh[9] = 5.9900000000000002 / xxm + 25. +
2756 .0099999999999997868 / xxm + 17.990000000000002 / xxm + 9. / xxm;
2757 dgh[10] = 36.;
2758 dgh[11] = 49.998;
2759 dgh[12] = 5.9900000000000002 / xxm + 25. +
2760 .0099999999999997868 / xxm + 17.990000000000002 / xxm + 9. / xxm + 13.997999999999998 / xxm;
2761 dgh[13] = 49.998;
2762 dgh[14] = 49.998;
2763 geant->Gsvolu("IT56", "PCON", idtmed[275], dgh, 15);
2764 geant->Gspos("IT56", 1, "ITSV", 0., 0., 0., 0, "ONLY");
2765 geant->Gsatt("IT56", "SEEN", 0);
2766
2767// --- Layer #5
2768
2769// GOTO 6678 ! skip ITS layer no. 5
2770
2771//--- Define a ghost volume containing a single ladder of layer #5 andfill
2772// it with air or vacuum
2773
2774 dbox1[0] = .044999999999999998;
2775 dbox1[1] = 3.75;
2776 dbox1[2] = 43.155000000000001;
2777 geant->Gsvolu("ISV1", "BOX ", idtmed[253], dbox1, 3);
2778
2779// --- Make the ghost volume invisible
2780
2781 geant->Gsatt("ISV1", "SEEN", 0);
2782
2783// --- Define a ghost volume containing the electronics and cooling of
2784// a single ladder of layer #5 and fill it with air or vacuum
2785
2786 dsrv[0] = (TMath::Sqrt(3.) / 2. * 4.2 + .47 + .05) / 2.;
2787 dsrv[1] = 3.75;
2788 dsrv[2] = 43.155000000000001;
2789 geant->Gsvolu("SSV1", "BOX ", idtmed[253], dsrv, 3);
2790
2791// --- Make the ghost volume invisible
2792
2793 geant->Gsatt("SSV1", "SEEN", 0);
2794
2795// --- Define a ghost volume containing the end-ladder stuff of
2796// a single ladder of layer #5 and fill it with air or vacuum
2797
2798 dela[0] = 2.;
2799 dela[1] = 3.5;
2800 dela[2] = 4.;
2801 geant->Gsvolu("ELL5", "BOX ", idtmed[253], dela, 3);
2802
2803// --- Make the ghost volume invisible
2804
2805 geant->Gsatt("ELL5", "SEEN", 0);
2806
2807// --- Define a volume containing the sensitive part of the strips
2808// (silicon, layer #5)
2809
2810 dits[0] = .015;
2811 dits[1] = 3.75;
2812 dits[2] = 2.1;
2813 geant->Gsvolu("ITS5", "BOX ", idtmed[249], dits, 3);
2814
2815// --- Define a volume containing the electronics of the strips
2816// (silicon, layer #5)
2817
2818 dchi[0] = .02;
2819 dchi[1] = 3.4;
2820 dchi[2] = .525;
2821 geant->Gsvolu("SCH5", "BOX ", idtmed[250], dchi, 3);
2822
2823// --- Define the cooling tubes (aluminum, layer #5)
2824
2825 dtub[0] = .09;
2826 dtub[1] = dtub[0] + .01;
2827 dtub[2] = 43.155000000000001;
2828 geant->Gsvolu("STB5", "TUBE", idtmed[255], dtub, 3);
2829
2830// --- Define the cooling fluid (water or freon, layer #5)
2831
2832 dwat[0] = 0.;
2833 dwat[1] = .09;
2834 dwat[2] = 43.155000000000001;
2835 geant->Gsvolu("SWT5", "TUBE", idtmed[256], dwat, 3);
2836// CALL GSVOLU('SWT5','TUBE',IDTMED(258),DWAT,3,IOUT) ! freon
2837
2838//--- Define the (triangular) element of the heat bridge (carbon, layer #5)
2839
2840// water
2841 dfra[0] = 120.;
2842 dfra[1] = 360.;
2843 dfra[2] = 3.;
2844 dfra[3] = 2.;
2845 dfra[4] = -.015;
2846 dfra[5] = TMath::Sqrt(3.) * 4.2 / 6.;
2847 dfra[6] = dfra[5] + .03;
2848 dfra[7] = .015;
2849 dfra[8] = dfra[5];
2850 dfra[9] = dfra[6];
2851 geant->Gsvolu("SFR5", "PGON", idtmed[252], dfra, 10);
2852
2853// --- Define the element connecting the triangles of the heat bridge
2854// (carbon, layer #5)
2855
2856 dcei[0] = 0.;
2857 dcei[1] = .03;
2858 dcei[2] = 43.155000000000001;
2859 geant->Gsvolu("SCE5", "TUBE", idtmed[252], dcei, 3);
2860
2861// --- Define the part of the end-ladder stuff made of plastic (G10FR4)
2862// (layer #5)
2863
2864 dpla[0] = .089285714285714288;
2865 dpla[1] = 3.5;
2866 dpla[2] = 4.;
2867 geant->Gsvolu("EPL5", "BOX ", idtmed[262], dpla, 3);
2868
2869// --- Define the part of the end-ladder stuff made of copper (layer #5)
2870
2871 dcop[0] = .017857142857142856;
2872 dcop[1] = 3.5;
2873 dcop[2] = 4.;
2874 geant->Gsvolu("ECU5", "BOX ", idtmed[259], dcop, 3);
2875
2876// --- Define the part of the end-ladder stuff made of epoxy (layer #5)
2877
2878 depx[0] = .26785714285714285;
2879 depx[1] = 3.5;
2880 depx[2] = 4.;
2881 geant->Gsvolu("EPX5", "BOX ", idtmed[262], depx, 3);
2882
2883// --- Define the part of the end-ladder stuff made of silicon (bus)
2884// (layer #5)
2885
2886 dsil[0] = .17857142857142858;
2887 dsil[1] = 3.5;
2888 dsil[2] = 4.;
2889 geant->Gsvolu("ESI5", "BOX ", idtmed[251], dsil, 3);
2890
2891// --- Place the end-ladder stuff into its mother (ELL5)
2892
2893 sep = (4. - (dpla[0] + dcop[0] + depx[0] + dsil[0]) * 2.) / 3.;
2894 ypos = 0.;
2895 zpos = 0.;
2896
2897// --- Plastic
2898
2899 xpos = -dela[0] + dpla[0];
2900 geant->Gspos("EPL5", 1, "ELL5", xpos, ypos, zpos, 0, "ONLY");
2901
2902// --- Copper
2903
2904 xpos = -dela[0] + dpla[0] * 2. + sep + dcop[0];
2905 geant->Gspos("ECU5", 1, "ELL5", xpos, ypos, zpos, 0, "ONLY");
2906
2907// --- Epoxy
2908
2909 xpos = -dela[0] + dpla[0] * 2. + sep + dcop[0] * 2. + sep + depx[0];
2910 geant->Gspos("EPX5", 1, "ELL5", xpos, ypos, zpos, 0, "ONLY");
2911
2912// --- Silicon (bus)
2913
2914 xpos = -dela[0] + dpla[0] * 2. + sep + dcop[0] * 2. + sep + depx[0] * 2. + sep + dsil[0];
2915 geant->Gspos("ESI5", 1, "ELL5", xpos, ypos, zpos, 0, "ONLY");
2916
2917// --- Place the sensitive part of the strips into its mother (ISV1)
2918
2919 ypos = 0.;
2920 for (j = 1; j <= 22; ++j) {
2921 if (j % 2 == 0) xpos = dbox1[0] - dits[0];
2922 else xpos = -dbox1[0] + dits[0];
2923 zpos = ((j - 1) - 10.) * 3.91 - 1.96;
2924 geant->Gspos("ITS5", j, "ISV1", xpos, ypos, zpos, 0, "ONLY");
2925 }
2926
2927// --- Place the electronics of the strips into its mother (SSV1)
2928
2929 ypos = 0.;
2930 for (j = 1; j <= 22; ++j) {
2931 if (j % 2 == 0) xpos = -dsrv[0] + .28;
2932 else xpos = -dsrv[0] + .28 - dits[0] * 2. - .03;
2933 zpos = ((j - 1) - 10.) * 3.91 - 1.96 + .85;
2934 geant->Gspos("SCH5", j, "SSV1", xpos, ypos, zpos, 0, "ONLY");
2935 }
2936
2937//--- Place the cooling tubes and the cooling fluid into their mother (SSV1)
2938
2939 xpos = -dsrv[0] + .41;
2940 zpos = 0.;
2941
2942// --- Left tube (just a matter of convention)
2943
2944 ypos = -2.3500000000000001;
2945 geant->Gspos("STB5", 1, "SSV1", xpos, ypos, zpos, 0, "ONLY");
2946 geant->Gspos("SWT5", 1, "SSV1", xpos, ypos, zpos, 0, "ONLY");
2947
2948// --- Right tube (just a matter of convention)
2949
2950 ypos = 2.3500000000000001;
2951 geant->Gspos("STB5", 2, "SSV1", xpos, ypos, zpos, 0, "ONLY");
2952 geant->Gspos("SWT5", 2, "SSV1", xpos, ypos, zpos, 0, "ONLY");
2953
2954// --- Place the heat bridge elements into their mother (SSV1)
2955
2956 xpos = -dsrv[0] + .47 + TMath::Sqrt(3.) / 6. * 4.2;
2957 ypos = 0.;
2958 for (j = 1; j <= 23; ++j) {
2959 zpos = ((j - 1) - 10.) * 3.91 - 1.96 - 2.1000000000000001;
2960 geant->Gspos("SFR5", j, "SSV1", xpos, ypos, zpos, 0, "ONLY");
2961 }
2962
2963// --- Place the elements connecting the triangles of the heat bridge
2964// into their mother (SSV1)
2965
2966 zpos = 0.;
2967
2968// --- Left element (just a matter of convention)
2969
2970 xpos = -dsrv[0] + .47;
2971 ypos = -2.1150000000000002;
2972 geant->Gspos("SCE5", 1, "SSV1", xpos, ypos, zpos, 0, "ONLY");
2973
2974// --- Right element
2975
2976 xpos = -dsrv[0] + .47;
2977 ypos = 2.1150000000000002;
2978 geant->Gspos("SCE5", 2, "SSV1", xpos, ypos, zpos, 0, "ONLY");
2979
2980// --- Top element
2981
2982 xpos = -dsrv[0] + .47 + TMath::Sqrt(3.) / 2. * 4.2 + .015;
2983 ypos = 0.;
2984 geant->Gspos("SCE5", 3, "SSV1", xpos, ypos, zpos, 0, "ONLY");
2985
2986// --- Place the ghost volumes containing the strip ladders (ISV1),
2987// electronics/cooling (SSV1) and end-ladder stuff (ELL5) of layer #5 in
2988// their mother volume (IT56)
2989
2990 offset1 = atg_(.8, 36.6);
2991 offset2 = 5.2;
2992 rzero = dbox1[0] + 36.6;
2993 runo = dbox1[0] * 2. + 36.6 + dsrv[0];
2994 rtwo = dbox1[0] * 2. + 36.6 + dela[0];
2995 for (i = 1; i <= 32; ++i) {
2996 atheta = (i-1) * twopi * raddeg / 32. + offset2;
2997 AliMatrix(idrotm[i+1499], 90., atheta, 90., atheta + 90., 0., 0.);
2998
2999// --- Strip ladders
3000
3001 xpos = rzero * TMath::Cos((i-1) * twopi / 32. + offset1);
3002 ypos = rzero * TMath::Sin((i-1) * twopi / 32. + offset1);
3003 zpos = 0.;
3004 geant->Gspos("ISV1", i, "IT56", xpos, ypos, zpos, idrotm[i+1499], "ONLY");
3005
3006// --- Electronics/cooling
3007
3008 xpos = runo * TMath::Cos((i-1) * twopi / 32. + offset1);
3009 ypos = runo * TMath::Sin((i-1) * twopi / 32. + offset1);
3010 zpos = 0.;
3011 geant->Gspos("SSV1", i, "IT56", xpos, ypos, zpos, idrotm[i+1499], "ONLY");
3012
3013// --- End-ladders (nagative-Z and positive-Z)
3014
3015 xpos = rtwo * TMath::Cos((i-1) * twopi / 32. + offset1);
3016 ypos = rtwo * TMath::Sin((i-1) * twopi / 32. + offset1);
3017 zpos = -(dbox1[2] + dela[2] + 6.);
3018 geant->Gspos("ELL5", i, "IT56", xpos, ypos, zpos, idrotm[i+1499], "ONLY");
3019 zpos = dbox1[2] + dela[2] + 6.;
3020 geant->Gspos("ELL5", i + 35, "IT56", xpos, ypos, zpos, idrotm[i+1499], "ONLY");
3021 }
3022
3023
3024// --- Layer #6
3025
3026// GOTO 6778 ! skip ITS layer no. 6
3027
3028//--- Define a ghost volume containing a single ladder of layer #6 andfill
3029// it with air or vacuum
3030
3031 dbox2[0] = .044999999999999998;
3032 dbox2[1] = 3.75;
3033 dbox2[2] = 47.064999999999998;
3034 geant->Gsvolu("ISV2", "BOX ", idtmed[253], dbox2, 3);
3035
3036// --- Make the ghost volume invisible
3037
3038 geant->Gsatt("ISV2", "SEEN", 0);
3039
3040// --- Define a ghost volume containing the electronics and cooling of
3041// a single ladder of layer #6 and fill it with air or vacuum
3042
3043 dsrv[0] = (TMath::Sqrt(3.) / 2. * 4.2 + .47 + .05) / 2.;
3044 dsrv[1] = 3.75;
3045 dsrv[2] = 47.064999999999998;
3046 geant->Gsvolu("SSV2", "BOX ", idtmed[253], dsrv, 3);
3047
3048// --- Make the ghost volume invisible
3049
3050 geant->Gsatt("SSV2", "SEEN", 0);
3051
3052// --- Define a ghost volume containing the end-ladder stuff of
3053// a single ladder of layer #6 and fill it with air or vacuum
3054
3055 dela[0] = 2.;
3056 dela[1] = 3.5;
3057 dela[2] = 4.;
3058 geant->Gsvolu("ELL6", "BOX ", idtmed[253], dela, 3);
3059
3060// --- Make the ghost volume invisible
3061
3062 geant->Gsatt("ELL6", "SEEN", 0);
3063
3064// --- Define a volume containing the sensitive part of the strips
3065// (silicon, layer #6)
3066
3067 dits[0] = .015;
3068 dits[1] = 3.75;
3069 dits[2] = 2.1;
3070 geant->Gsvolu("ITS6", "BOX ", idtmed[249], dits, 3);
3071
3072// --- Define a volume containing the electronics of the strips
3073// (silicon, layer #6)
3074
3075 dchi[0] = .02;
3076 dchi[1] = 3.4;
3077 dchi[2] = .525;
3078 geant->Gsvolu("SCH6", "BOX ", idtmed[250], dchi, 3);
3079
3080// --- Define the cooling tubes (aluminum, layer #6)
3081
3082 dtub[0] = .09;
3083 dtub[1] = dtub[0] + .01;
3084 dtub[2] = 47.064999999999998;
3085 geant->Gsvolu("STB6", "TUBE", idtmed[255], dtub, 3);
3086
3087// --- Define the cooling fluid (water or freon, layer #6)
3088
3089 dwat[0] = 0.;
3090 dwat[1] = .09;
3091 dwat[2] = 47.064999999999998;
3092 geant->Gsvolu("SWT6", "TUBE", idtmed[256], dwat, 3);
3093// CALL GSVOLU('SWT6','TUBE',IDTMED(258),DWAT,3,IOUT) ! freon
3094
3095//--- Define the (triangular) element of the heat bridge (carbon, layer #6)
3096
3097// water
3098 dfra[0] = 120.;
3099 dfra[1] = 360.;
3100 dfra[2] = 3.;
3101 dfra[3] = 2.;
3102 dfra[4] = -.015;
3103 dfra[5] = TMath::Sqrt(3.) * 4.2 / 6.;
3104 dfra[6] = dfra[5] + .03;
3105 dfra[7] = .015;
3106 dfra[8] = dfra[5];
3107 dfra[9] = dfra[6];
3108 geant->Gsvolu("SFR6", "PGON", idtmed[252], dfra, 10);
3109
3110// --- Define the element connecting the triangles of the heat bridge
3111// (carbon, layer #6)
3112
3113 dcei[0] = 0.;
3114 dcei[1] = .03;
3115 dcei[2] = 47.064999999999998;
3116 geant->Gsvolu("SCE6", "TUBE", idtmed[252], dcei, 3);
3117
3118// --- Define the part of the end-ladder stuff made of plastic (G10FR4)
3119// (layer #6)
3120
3121 dpla[0] = .089285714285714288;
3122 dpla[1] = 3.5;
3123 dpla[2] = 4.;
3124 geant->Gsvolu("EPL6", "BOX ", idtmed[262], dpla, 3);
3125
3126// --- Define the part of the end-ladder stuff made of copper (layer #6)
3127
3128 dcop[0] = .017857142857142856;
3129 dcop[1] = 3.5;
3130 dcop[2] = 4.;
3131 geant->Gsvolu("ECU6", "BOX ", idtmed[259], dcop, 3);
3132
3133// --- Define the part of the end-ladder stuff made of epoxy (layer #6)
3134
3135 depx[0] = .26785714285714285;
3136 depx[1] = 3.5;
3137 depx[2] = 4.;
3138 geant->Gsvolu("EPX6", "BOX ", idtmed[262], depx, 3);
3139
3140// --- Define the part of the end-ladder stuff made of silicon (bus)
3141// (layer #6)
3142
3143 dsil[0] = .17857142857142858;
3144 dsil[1] = 3.5;
3145 dsil[2] = 4.;
3146 geant->Gsvolu("ESI6", "BOX ", idtmed[251], dsil, 3);
3147
3148// --- Place the end-ladder stuff into its mother (ELL5)
3149
3150 sep= (4. - (dpla[0] + dcop[0] + depx[0] + dsil[0]) * 2.) / 3.;
3151 ypos = 0.;
3152 zpos = 0.;
3153
3154// --- Plastic
3155
3156 xpos = -dela[0] + dpla[0];
3157 geant->Gspos("EPL6", 1, "ELL6", xpos, ypos, zpos, 0, "ONLY");
3158
3159// --- Copper
3160
3161 xpos = -dela[0] + dpla[0] * 2. + sep + dcop[0];
3162 geant->Gspos("ECU6", 1, "ELL6", xpos, ypos, zpos, 0, "ONLY");
3163
3164// --- Epoxy
3165
3166 xpos = -dela[0] + dpla[0] * 2. + sep + dcop[0] * 2. + sep + depx[0];
3167 geant->Gspos("EPX6", 1, "ELL6", xpos, ypos, zpos, 0, "ONLY");
3168
3169// --- Silicon (bus)
3170
3171 xpos = -dela[0] + dpla[0] * 2. + sep + dcop[0] * 2. + sep + depx[0] * 2. + sep + dsil[0];
3172 geant->Gspos("ESI6", 1, "ELL6", xpos, ypos, zpos, 0, "ONLY");
3173
3174// --- Place the sensitive part of the strips into its mother (ISV2)
3175
3176 ypos = 0.;
3177 for (j = 1; j <= 24; ++j) {
3178 if (j % 2 == 0) xpos = -dbox2[0] + dits[0];
3179 else xpos = dbox2[0] - dits[0];
3180 zpos = ((j - 1) - 11.) * 3.91 - 1.96;
3181 geant->Gspos("ITS6", j, "ISV2", xpos, ypos, zpos, 0, "ONLY");
3182 }
3183
3184// --- Place the electronics of the strips into its mother (SSV2)
3185
3186 ypos = 0.;
3187 for (j = 1; j <= 24; ++j) {
3188 if (j % 2 == 0) xpos = -dsrv[0] + .28 - dits[0] * 2. - .03;
3189 else xpos = -dsrv[0] + .28;
3190 zpos = ((j - 1) - 11.) * 3.91 - 1.96 + .85;
3191 geant->Gspos("SCH5", j, "SSV1", xpos, ypos, zpos, 0, "ONLY");
3192 }
3193
3194//--- Place the cooling tubes and the cooling fluid into their mother (SSV2)
3195
3196 xpos = -dsrv[0] + .41;
3197 zpos = 0.;
3198
3199// --- Left tube (just a matter of convention)
3200
3201 ypos = -2.3500000000000001;
3202 geant->Gspos("STB6", 1, "SSV2", xpos, ypos, zpos, 0, "ONLY");
3203 geant->Gspos("SWT6", 1, "SSV2", xpos, ypos, zpos, 0, "ONLY");
3204
3205// --- Right tube (just a matter of convention)
3206
3207 ypos = 2.3500000000000001;
3208 geant->Gspos("STB6", 2, "SSV2", xpos, ypos, zpos, 0, "ONLY");
3209 geant->Gspos("SWT6", 2, "SSV2", xpos, ypos, zpos, 0, "ONLY");
3210
3211// --- Place the heat bridge elements into their mother (SSV2)
3212
3213 xpos = -dsrv[0] + .47 + TMath::Sqrt(3.) / 6. * 4.2;
3214 ypos = 0.;
3215 for (j = 1; j <= 25; ++j) {
3216 zpos = ((j - 1) - 11.) * 3.91 - 1.96 - 2.1000000000000001;
3217 geant->Gspos("SFR6", j, "SSV2", xpos, ypos, zpos, 0, "ONLY");
3218 }
3219
3220// --- Place the elements connecting the triangles of the heat bridge
3221// into their mother (SSV2)
3222
3223 zpos = 0.;
3224
3225// --- Left element (just a matter of convention)
3226
3227 xpos = -dsrv[0] + .47;
3228 ypos = -2.1150000000000002;
3229 geant->Gspos("SCE6", 1, "SSV2", xpos, ypos, zpos, 0, "ONLY");
3230
3231// --- Right element
3232
3233 xpos = -dsrv[0] + .47;
3234 ypos = 2.1150000000000002;
3235 geant->Gspos("SCE6", 2, "SSV2", xpos, ypos, zpos, 0, "ONLY");
3236
3237// --- Top element
3238
3239 xpos = -dsrv[0] + .47 + TMath::Sqrt(3.) / 2. * 4.2 + .015;
3240 ypos = 0.;
3241 geant->Gspos("SCE6", 3, "SSV2", xpos, ypos, zpos, 0, "ONLY");
3242
3243// --- Place the ghost volumes containing the strip ladders (ISV2),
3244// electronics/cooling (SSV2) and end-ladder stuff (ELL6) of layer #6 in
3245// their mother volume (IT56)
3246
3247 offset1 = atg_(.9, 41.2);
3248 offset2 = 5.2;
3249 rzero = dbox2[0] + 41.2;
3250 runo = dbox2[0] * 2. + 41.2 + dsrv[0];
3251 rtwo = dbox2[0] * 2. + 41.2 + dela[0];
3252 for (i = 1; i <= 36; ++i) {
3253 atheta = (i-1) * twopi * raddeg / 36. + offset2;
3254 AliMatrix(idrotm[i+1599], 90., atheta, 90., atheta + 90., 0., 0.);
3255
3256// --- Strip ladders
3257
3258 xpos = rzero * TMath::Cos((i-1) * twopi / 36. + offset1);
3259 ypos = rzero * TMath::Sin((i-1) * twopi / 36. + offset1);
3260 zpos = 0.;
3261 geant->Gspos("ISV2", i, "IT56", xpos, ypos, zpos, idrotm[i+1599], "ONLY");
3262
3263// --- Electronics/cooling
3264
3265 xpos = runo * TMath::Cos((i-1) * twopi / 36. + offset1);
3266 ypos = runo * TMath::Sin((i-1) * twopi / 36. + offset1);
3267 zpos = 0.;
3268 geant->Gspos("SSV2", i, "IT56", xpos, ypos, zpos, idrotm[i+1599], "ONLY");
3269
3270// --- End-ladders (nagative-Z and positive-Z)
3271
3272 xpos = rtwo * TMath::Cos((i-1) * twopi / 36. + offset1);
3273 ypos = rtwo * TMath::Sin((i-1) * twopi / 36. + offset1);
3274 zpos = -(dbox2[2] + dela[2] + 6.);
3275 geant->Gspos("ELL6", i, "IT56", xpos, ypos, zpos, idrotm[i+1599], "ONLY");
3276 zpos = dbox2[2] + dela[2] + 6.;
3277 geant->Gspos("ELL6", i + 39, "IT56", xpos, ypos, zpos, idrotm[i+1599], "ONLY");
3278 }
3279
3280
3281 }
3282
3283//************************************************************************
3284//* *
3285//* E N D - C A P S A N D F R A M E S *
3286//* ========================================= *
3287//* *
3288//************************************************************************
3289
3290// --- Define a dummy cylinder for multiple scattering tests
3291
3292// GOTO 7890 ! skip dummy cylinder for multiple scatteringtests
3293
3294// DITS(1)=49.
3295// DITS(2)=DITS(1)+0.1
3296// DITS(3)=60.3
3297// CALL GSVOLU('ITST','TUBE',IDTMED(255),DITS,3,IOUT)
3298// CALL GSPOS('ITST',1,'ITSV',0.,0.,0.,0,'ONLY')
3299// 7890 CONTINUE
3300
3301// --- The 0.74% X0 outer wall (C) of the gas vessel at r=50cm ---
3302
3303 if (fMinorVersion == 0 || fMinorVersion == 3) {
3304
3305// GOTO 8901 ! skip outer wall
3306
3307 dits[0] = 49.9;
3308 dits[1] = dits[0] + .06926;
3309 dits[2] = dpcb[2] * 2. + 62.7 - 10.5;
3310// old value 60.3
3311 geant->Gsvolu("ITSG", "TUBE", idtmed[274], dits, 3);
3312 geant->Gspos("ITSG", 1, "ITSV", 0., 0., 0., 0, "ONLY");
3313 } else {
3314 goto L8901;
3315 }
3316L8901:
3317// --- The frame between the end-caps (octagonal lay-out) ---
3318
3319// GOTO 9012 ! skip octagonal frame
3320
3321 if (fMinorVersion == 1) {
3322
3323 rzero = 34.;
3324 dtra[0] = .92;
3325 dtra[1] = 1.;
3326 dtra[2] = dpcb[2] * 2. + 50.5 - 10.5;
3327 dtra1[0] = .92;
3328 dtra1[1] = 1.;
3329 dtra1[2] = TMath::Sqrt(dtra[2] * dtra[2] + 518.90999999999985) / 2.;
3330 angle = 45.;
3331 offset = angle / 2.;
3332 for (i = 1; i <= 8; ++i) {
3333 xtra[i - 1] = rzero * TMath::Cos((i-1) * angle * degrad);
3334 ytra[i - 1] = rzero * TMath::Sin((i-1) * angle * degrad);
3335 ztra[i - 1] = 0.;
3336 geant->Gsvolu(natra_ref(0, i), "TUBE", idtmed[274], dtra, 3);
3337 geant->Gspos(natra_ref(0, i), 1, "ITSV", xtra[i - 1], ytra[i - 1], ztra[i - 1], 0, "ONLY");
3338 }
3339
3340 atheta = 22.5;
3341 aphi1 = TMath::ACos(dtra[2] / TMath::Sqrt(dtra[2] * dtra[2] + (50.5 / cos(.48869219055841229) * (50.5 / cos(.48869219055841229))
3342 - 2550.25))) * raddeg;
3343 aphi2 = 180. - aphi1;
3344 xpos = (xtra[0] + xtra[1]) / 2.;
3345 ypos = (ytra[0] + ytra[1]) / 2.;
3346 zpos = dtra[2] / 2.;
3347 geant->Gsvolu(natra1_ref(0, 1), "TUBE", idtmed[274], dtra1, 3);
3348 r2 = atheta + 90.;
3349 r3 = atheta + 90.;
3350 AliMatrix(idrotm[5100], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3351 geant->Gspos(natra1_ref(0, 1), 1, "ITSV", xpos, ypos, zpos, idrotm[5100], "ONLY");
3352 zpos = -dtra[2] / 2.;
3353 geant->Gsvolu(natra1_ref(0, 2), "TUBE", idtmed[274], dtra1, 3);
3354 r2 = atheta + 90.;
3355 r3 = atheta + 90.;
3356 AliMatrix(idrotm[5101], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3357 geant->Gspos(natra1_ref(0, 2), 1, "ITSV", xpos, ypos, zpos, idrotm[5101], "ONLY");
3358
3359 atheta = 67.5;
3360 aphi2 = TMath::ACos(dtra[2] / TMath::Sqrt(dtra[2] * dtra[2] + (50.5 / cos(.48869219055841229) * (50.5 / cos(.48869219055841229))
3361 - 2550.25))) * raddeg;
3362 aphi1 = 180. - aphi2;
3363 xpos = (xtra[1] + xtra[2]) / 2.;
3364 ypos = (ytra[1] + ytra[2]) / 2.;
3365 zpos = dtra[2] / 2.;
3366 geant->Gsvolu(natra1_ref(0, 3), "TUBE", idtmed[274], dtra1, 3);
3367 r2 = atheta + 90.;
3368 r3 = atheta + 90.;
3369 AliMatrix(idrotm[5102], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3370 geant->Gspos(natra1_ref(0, 3), 1, "ITSV", xpos, ypos, zpos, idrotm[5102], "ONLY");
3371 zpos = -dtra[2] / 2.;
3372 geant->Gsvolu(natra1_ref(0, 4), "TUBE", idtmed[274], dtra1, 3);
3373 r2 = atheta + 90.;
3374 r3 = atheta + 90.;
3375 AliMatrix(idrotm[5103], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3376 geant->Gspos(natra1_ref(0, 4), 1, "ITSV", xpos, ypos, zpos, idrotm[5103], "ONLY");
3377
3378 atheta = 112.5;
3379 aphi1 = TMath::ACos(dtra[2] / TMath::Sqrt(dtra[2] * dtra[2] + (50.5 / cos(.48869219055841229) * (50.5 / cos(.48869219055841229))
3380 - 2550.25))) * raddeg;
3381 aphi2 = 180. - aphi1;
3382 xpos = (xtra[2] + xtra[3]) / 2.;
3383 ypos = (ytra[2] + ytra[3]) / 2.;
3384 zpos = dtra[2] / 2.;
3385 geant->Gsvolu(natra1_ref(0, 5), "TUBE", idtmed[274], dtra1, 3);
3386 r2 = atheta + 90.;
3387 r3 = atheta + 90.;
3388 AliMatrix(idrotm[5104], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3389 geant->Gspos(natra1_ref(0, 5), 1, "ITSV", xpos, ypos, zpos, idrotm[5104], "ONLY");
3390 zpos = -dtra[2] / 2.;
3391 geant->Gsvolu(natra1_ref(0, 6), "TUBE", idtmed[274], dtra1, 3);
3392 r2 = atheta + 90.;
3393 r3 = atheta + 90.;
3394 AliMatrix(idrotm[5105], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3395 geant->Gspos(natra1_ref(0, 6), 1, "ITSV", xpos, ypos, zpos, idrotm[5105], "ONLY");
3396
3397 atheta = 157.5;
3398 aphi2 = TMath::ACos(dtra[2] / TMath::Sqrt(dtra[2] * dtra[2] + (50.5 / cos(.48869219055841229) * (50.5 / cos(.48869219055841229))
3399 - 2550.25))) * raddeg;
3400 aphi1 = 180. - aphi2;
3401 xpos = (xtra[3] + xtra[4]) / 2.;
3402 ypos = (ytra[3] + ytra[4]) / 2.;
3403 zpos = dtra[2] / 2.;
3404 geant->Gsvolu(natra1_ref(0, 7), "TUBE", idtmed[274], dtra1, 3);
3405 r2 = atheta + 90.;
3406 r3 = atheta + 90.;
3407 AliMatrix(idrotm[5106], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3408 geant->Gspos(natra1_ref(0, 7), 1, "ITSV", xpos, ypos, zpos, idrotm[5106], "ONLY");
3409 zpos = -dtra[2] / 2.;
3410 geant->Gsvolu(natra1_ref(0, 8), "TUBE", idtmed[274], dtra1, 3);
3411 r2 = atheta + 90.;
3412 r3 = atheta + 90.;
3413 AliMatrix(idrotm[5107], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3414 geant->Gspos(natra1_ref(0, 8), 1, "ITSV", xpos, ypos, zpos, idrotm[5107], "ONLY");
3415
3416 atheta = 22.5;
3417 aphi2 = TMath::ACos(dtra[2] / TMath::Sqrt(dtra[2] * dtra[2] + (50.5 / cos(.48869219055841229) * (50.5 / cos(.48869219055841229))
3418 - 2550.25))) * raddeg;
3419 aphi1 = 180. - aphi2;
3420 xpos = (xtra[4] + xtra[5]) / 2.;
3421 ypos = (ytra[4] + ytra[5]) / 2.;
3422 zpos = dtra[2] / 2.;
3423 geant->Gsvolu(natra1_ref(0, 9), "TUBE", idtmed[274], dtra1, 3);
3424 r2 = atheta + 90.;
3425 r3 = atheta + 90.;
3426 AliMatrix(idrotm[5108], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3427 geant->Gspos(natra1_ref(0, 9), 1, "ITSV", xpos, ypos, zpos, idrotm[5108], "ONLY");
3428 zpos = -dtra[2] / 2.;
3429 geant->Gsvolu(natra1_ref(0, 10), "TUBE", idtmed[274], dtra1, 3);
3430 r2 = atheta + 90.;
3431 r3 = atheta + 90.;
3432 AliMatrix(idrotm[5109], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3433 geant->Gspos(natra1_ref(0, 10), 1, "ITSV", xpos, ypos, zpos, idrotm[5109], "ONLY");
3434
3435 atheta = 67.5;
3436 aphi1 = TMath::ACos(dtra[2] / TMath::Sqrt(dtra[2] * dtra[2] + (50.5 / cos(.48869219055841229) * (50.5 / cos(.48869219055841229))
3437 - 2550.25))) * raddeg;
3438 aphi2 = 180. - aphi1;
3439 xpos = (xtra[5] + xtra[6]) / 2.;
3440 ypos = (ytra[5] + ytra[6]) / 2.;
3441 zpos = dtra[2] / 2.;
3442 geant->Gsvolu(natra1_ref(0, 11), "TUBE", idtmed[274], dtra1, 3);
3443 r2 = atheta + 90.;
3444 r3 = atheta + 90.;
3445 AliMatrix(idrotm[5110], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3446 geant->Gspos(natra1_ref(0, 11), 1, "ITSV", xpos, ypos, zpos, idrotm[5110], "ONLY");
3447 zpos = -dtra[2] / 2.;
3448 geant->Gsvolu(natra1_ref(0, 12), "TUBE", idtmed[274], dtra1, 3);
3449 r2 = atheta + 90.;
3450 r3 = atheta + 90.;
3451 AliMatrix(idrotm[5111], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3452 geant->Gspos(natra1_ref(0, 12), 1, "ITSV", xpos, ypos, zpos, idrotm[5111], "ONLY");
3453
3454 atheta = 112.5;
3455 aphi2 = TMath::ACos(dtra[2] / TMath::Sqrt(dtra[2] * dtra[2] + (50.5 / cos(.48869219055841229) * (50.5 / cos(.48869219055841229))
3456 - 2550.25))) * raddeg;
3457 aphi1 = 180. - aphi2;
3458 xpos = (xtra[6] + xtra[7]) / 2.;
3459 ypos = (ytra[6] + ytra[7]) / 2.;
3460 zpos = dtra[2] / 2.;
3461 geant->Gsvolu(natra1_ref(0, 13), "TUBE", idtmed[274], dtra1, 3);
3462 r2 = atheta + 90.;
3463 r3 = atheta + 90.;
3464 AliMatrix(idrotm[5112], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3465 geant->Gspos(natra1_ref(0, 13), 1, "ITSV", xpos, ypos, zpos, idrotm[5112], "ONLY");
3466 zpos = -dtra[2] / 2.;
3467 geant->Gsvolu(natra1_ref(0, 14), "TUBE", idtmed[274], dtra1, 3);
3468 r2 = atheta + 90.;
3469 r3 = atheta + 90.;
3470 AliMatrix(idrotm[5113], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3471 geant->Gspos(natra1_ref(0, 14), 1, "ITSV", xpos, ypos, zpos, idrotm[5113], "ONLY");
3472
3473 atheta = 157.5;
3474 aphi1 = TMath::ACos(dtra[2] / TMath::Sqrt(dtra[2] * dtra[2] + (50.5 / cos(.48869219055841229) * (50.5 / cos(.48869219055841229))
3475 - 2550.25))) * raddeg;
3476 aphi2 = 180. - aphi1;
3477 xpos = (xtra[7] + xtra[0]) / 2.;
3478 ypos = (ytra[7] + ytra[0]) / 2.;
3479 zpos = dtra[2] / 2.;
3480 geant->Gsvolu(natra1_ref(0, 15), "TUBE", idtmed[274], dtra1, 3);
3481 r2 = atheta + 90.;
3482 r3 = atheta + 90.;
3483 AliMatrix(idrotm[5114], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3484 geant->Gspos(natra1_ref(0, 15), 1, "ITSV", xpos, ypos, zpos, idrotm[5114], "ONLY");
3485 zpos = -dtra[2] / 2.;
3486 geant->Gsvolu(natra1_ref(0, 16), "TUBE", idtmed[274], dtra1, 3);
3487 r2 = atheta + 90.;
3488 r3 = atheta + 90.;
3489 AliMatrix(idrotm[5115], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3490 geant->Gspos(natra1_ref(0, 16), 1, "ITSV", xpos, ypos, zpos, idrotm[5115], "ONLY");
3491
3492
3493 } else if (fMinorVersion == 4) {
3494
3495
3496 rzero = 34.;
3497 dtra[0] = .92;
3498 dtra[1] = 1.;
3499 dtra[2] = dpcb[2] * 2. + 50.5 - 10.5;
3500 dtra1[0] = .92;
3501 dtra1[1] = 1.;
3502 dtra1[2] = TMath::Sqrt(dtra[2] * dtra[2] + 518.91) / 2.;
3503 angle = 45.;
3504 offset = angle / 2.;
3505 for (i = 1; i <= 8; ++i) {
3506 xtra[i - 1] = rzero * TMath::Cos((i-1) * angle * degrad);
3507 ytra[i - 1] = rzero * TMath::Sin((i-1) * angle * degrad);
3508 ztra[i - 1] = 0.;
3509 geant->Gsvolu(natra_ref(0, i), "TUBE", idtmed[274], dtra, 3);
3510 geant->Gspos(natra_ref(0, i), 1, "ITSV", xtra[i - 1], ytra[i - 1], ztra[i - 1], 0, "ONLY");
3511 }
3512
3513 atheta = 22.5;
3514 aphi1 = TMath::ACos(dtra[2] / TMath::Sqrt(dtra[2] * dtra[2] + (50.5 / cos(.48869219055841229) * (50.5 / cos(.48869219055841229))
3515 - 2550.25))) * raddeg;
3516 aphi2 = 180. - aphi1;
3517 xpos = (xtra[0] + xtra[1]) / 2.;
3518 ypos = (ytra[0] + ytra[1]) / 2.;
3519 zpos = dtra[2] / 2.;
3520 geant->Gsvolu(natra1_ref(0, 1), "TUBE", idtmed[274], dtra1, 3);
3521 r2 = atheta + 90.;
3522 r3 = atheta + 90.;
3523 AliMatrix(idrotm[5100], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3524 geant->Gspos(natra1_ref(0, 1), 1, "ITSV", xpos, ypos, zpos, idrotm[5100], "ONLY");
3525 zpos = -dtra[2] / 2.;
3526 geant->Gsvolu(natra1_ref(0, 2), "TUBE", idtmed[274], dtra1, 3);
3527 r2 = atheta + 90.;
3528 r3 = atheta + 90.;
3529 AliMatrix(idrotm[5101], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3530 geant->Gspos(natra1_ref(0, 2), 1, "ITSV", xpos, ypos, zpos, idrotm[5101], "ONLY");
3531
3532 atheta = 67.5;
3533 aphi2 = TMath::ACos(dtra[2] / TMath::Sqrt(dtra[2] * dtra[2] + (50.5 / cos(.48869219055841229) * (50.5 / cos(.48869219055841229))
3534 - 2550.25))) * raddeg;
3535 aphi1 = 180. - aphi2;
3536 xpos = (xtra[1] + xtra[2]) / 2.;
3537 ypos = (ytra[1] + ytra[2]) / 2.;
3538 zpos = dtra[2] / 2.;
3539 geant->Gsvolu(natra1_ref(0, 3), "TUBE", idtmed[274], dtra1, 3);
3540 r2 = atheta + 90.;
3541 r3 = atheta + 90.;
3542 AliMatrix(idrotm[5102], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3543 geant->Gspos(natra1_ref(0, 3), 1, "ITSV", xpos, ypos, zpos, idrotm[5102], "ONLY");
3544 zpos = -dtra[2] / 2.;
3545 geant->Gsvolu(natra1_ref(0, 4), "TUBE", idtmed[274], dtra1, 3);
3546 r2 = atheta + 90.;
3547 r3 = atheta + 90.;
3548 AliMatrix(idrotm[5103], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3549 geant->Gspos(natra1_ref(0, 4), 1, "ITSV", xpos, ypos, zpos, idrotm[5103], "ONLY");
3550
3551 atheta = 112.5;
3552 aphi1 = TMath::ACos(dtra[2] / TMath::Sqrt(dtra[2] * dtra[2] + (50.5 / cos(.48869219055841229) * (50.5 / cos(.48869219055841229))
3553 - 2550.25))) * raddeg;
3554 aphi2 = 180. - aphi1;
3555 xpos = (xtra[2] + xtra[3]) / 2.;
3556 ypos = (ytra[2] + ytra[3]) / 2.;
3557 zpos = dtra[2] / 2.;
3558 geant->Gsvolu(natra1_ref(0, 5), "TUBE", idtmed[274], dtra1, 3);
3559 r2 = atheta + 90.;
3560 r3 = atheta + 90.;
3561 AliMatrix(idrotm[5104], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3562 geant->Gspos(natra1_ref(0, 5), 1, "ITSV", xpos, ypos, zpos, idrotm[5104], "ONLY");
3563 zpos = -dtra[2] / 2.;
3564 geant->Gsvolu(natra1_ref(0, 6), "TUBE", idtmed[274], dtra1, 3);
3565 r2 = atheta + 90.;
3566 r3 = atheta + 90.;
3567 AliMatrix(idrotm[5105], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3568 geant->Gspos(natra1_ref(0, 6), 1, "ITSV", xpos, ypos, zpos, idrotm[5105], "ONLY");
3569
3570 atheta = 157.5;
3571 aphi2 = TMath::ACos(dtra[2] / TMath::Sqrt(dtra[2] * dtra[2] + (50.5 / cos(.48869219055841229) * (50.5 / cos(.48869219055841229))
3572 - 2550.25))) * raddeg;
3573 aphi1 = 180. - aphi2;
3574 xpos = (xtra[3] + xtra[4]) / 2.;
3575 ypos = (ytra[3] + ytra[4]) / 2.;
3576 zpos = dtra[2] / 2.;
3577 geant->Gsvolu(natra1_ref(0, 7), "TUBE", idtmed[274], dtra1, 3);
3578 r2 = atheta + 90.;
3579 r3 = atheta + 90.;
3580 AliMatrix(idrotm[5106], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3581 geant->Gspos(natra1_ref(0, 7), 1, "ITSV", xpos, ypos, zpos, idrotm[5106], "ONLY");
3582 zpos = -dtra[2] / 2.;
3583 geant->Gsvolu(natra1_ref(0, 8), "TUBE", idtmed[274], dtra1, 3);
3584 r2 = atheta + 90.;
3585 r3 = atheta + 90.;
3586 AliMatrix(idrotm[5107], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3587 geant->Gspos(natra1_ref(0, 8), 1, "ITSV", xpos, ypos, zpos, idrotm[5107], "ONLY");
3588
3589 atheta = 22.5;
3590 aphi2 = TMath::ACos(dtra[2] / TMath::Sqrt(dtra[2] * dtra[2] + (50.5 / cos(.48869219055841229) * (50.5 / cos(.48869219055841229))
3591 - 2550.25))) * raddeg;
3592 aphi1 = 180. - aphi2;
3593 xpos = (xtra[4] + xtra[5]) / 2.;
3594 ypos = (ytra[4] + ytra[5]) / 2.;
3595 zpos = dtra[2] / 2.;
3596 geant->Gsvolu(natra1_ref(0, 9), "TUBE", idtmed[274], dtra1, 3);
3597 r2 = atheta + 90.;
3598 r3 = atheta + 90.;
3599 AliMatrix(idrotm[5108], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3600 geant->Gspos(natra1_ref(0, 9), 1, "ITSV", xpos, ypos, zpos, idrotm[5108], "ONLY");
3601 zpos = -dtra[2] / 2.;
3602 geant->Gsvolu(natra1_ref(0, 10), "TUBE", idtmed[274], dtra1, 3);
3603 r2 = atheta + 90.;
3604 r3 = atheta + 90.;
3605 AliMatrix(idrotm[5109], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3606 geant->Gspos(natra1_ref(0, 10), 1, "ITSV", xpos, ypos, zpos, idrotm[5109], "ONLY");
3607
3608 atheta = 67.5;
3609 aphi1 = TMath::ACos(dtra[2] / TMath::Sqrt(dtra[2] * dtra[2] + (50.5 / cos(.48869219055841229) * (50.5 / cos(.48869219055841229))
3610 - 2550.25))) * raddeg;
3611 aphi2 = 180. - aphi1;
3612 xpos = (xtra[5] + xtra[6]) / 2.;
3613 ypos = (ytra[5] + ytra[6]) / 2.;
3614 zpos = dtra[2] / 2.;
3615 geant->Gsvolu(natra1_ref(0, 11), "TUBE", idtmed[274], dtra1, 3);
3616 r2 = atheta + 90.;
3617 r3 = atheta + 90.;
3618 AliMatrix(idrotm[5110], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3619 geant->Gspos(natra1_ref(0, 11), 1, "ITSV", xpos, ypos, zpos, idrotm[5110], "ONLY");
3620 zpos = -dtra[2] / 2.;
3621 geant->Gsvolu(natra1_ref(0, 12), "TUBE", idtmed[274], dtra1, 3);
3622 r2 = atheta + 90.;
3623 r3 = atheta + 90.;
3624 AliMatrix(idrotm[5111], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3625 geant->Gspos(natra1_ref(0, 12), 1, "ITSV", xpos, ypos, zpos, idrotm[5111], "ONLY");
3626
3627 atheta = 112.5;
3628 aphi2 = TMath::ACos(dtra[2] / TMath::Sqrt(dtra[2] * dtra[2] + (50.5 / cos(.48869219055841229) * (50.5 / cos(.48869219055841229))
3629 - 2550.25))) * raddeg;
3630 aphi1 = 180. - aphi2;
3631 xpos = (xtra[6] + xtra[7]) / 2.;
3632 ypos = (ytra[6] + ytra[7]) / 2.;
3633 zpos = dtra[2] / 2.;
3634 geant->Gsvolu(natra1_ref(0, 13), "TUBE", idtmed[274], dtra1, 3);
3635 r2 = atheta + 90.;
3636 r3 = atheta + 90.;
3637 AliMatrix(idrotm[5112], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3638 geant->Gspos(natra1_ref(0, 13), 1, "ITSV", xpos, ypos, zpos, idrotm[5112], "ONLY");
3639 zpos = -dtra[2] / 2.;
3640 geant->Gsvolu(natra1_ref(0, 14), "TUBE", idtmed[274], dtra1, 3);
3641 r2 = atheta + 90.;
3642 r3 = atheta + 90.;
3643 AliMatrix(idrotm[5113], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3644 geant->Gspos(natra1_ref(0, 14), 1, "ITSV", xpos, ypos, zpos, idrotm[5113], "ONLY");
3645
3646 atheta = 157.5;
3647 aphi1 = TMath::ACos(dtra[2] / TMath::Sqrt(dtra[2] * dtra[2] + (50.5 / cos(.48869219055841229) * (50.5 / cos(.48869219055841229))
3648 - 2550.25))) * raddeg;
3649 aphi2 = 180. - aphi1;
3650 xpos = (xtra[7] + xtra[0]) / 2.;
3651 ypos = (ytra[7] + ytra[0]) / 2.;
3652 zpos = dtra[2] / 2.;
3653 geant->Gsvolu(natra1_ref(0, 15), "TUBE", idtmed[274], dtra1, 3);
3654 r2 = atheta + 90.;
3655 r3 = atheta + 90.;
3656 AliMatrix(idrotm[5114], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3657 geant->Gspos(natra1_ref(0, 15), 1, "ITSV", xpos, ypos, zpos, idrotm[5114], "ONLY");
3658 zpos = -dtra[2] / 2.;
3659 geant->Gsvolu(natra1_ref(0, 16), "TUBE", idtmed[274], dtra1, 3);
3660 r2 = atheta + 90.;
3661 r3 = atheta + 90.;
3662 AliMatrix(idrotm[5115], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3663 geant->Gspos(natra1_ref(0, 16), 1, "ITSV", xpos, ypos, zpos, idrotm[5115], "ONLY");
3664 } else {
3665 goto L9012;
3666 }
3667
3668L9012:
3669
3670// --- The frame between the end-caps (hexagonal lay-out) ---
3671
3672// GOTO 9123 ! skip hexagonal frame
3673
3674 if (fMinorVersion == 2) {
3675
3676 rzero = 33.5;
3677 dtra2[0] = .92;
3678 dtra2[1] = 1.;
3679 dtra2[2] = dpcb[2] * 2. + 50. - 10.5;
3680 dtra3[0] = .92;
3681 dtra3[1] = 1.;
3682 dtra3[2] = 16.75;
3683 dtra4[0] = .92;
3684 dtra4[1] = 1.;
3685 dtra4[2] = TMath::Sqrt(dtra2[2] * dtra2[2] + 1088.0099999999998) / 2.;
3686 angle = 60.;
3687 offset = angle / 2.;
3688 for (i = 1; i <= 6; ++i) {
3689 xtra1[i - 1] = rzero * TMath::Cos(((i-1) * angle + offset) *degrad);
3690 ytra1[i - 1] = rzero * TMath::Sin(((i-1) * angle + offset) *degrad);
3691 ztra1[i - 1] = 0.;
3692 geant->Gsvolu(natra2_ref(0, i), "TUBE", idtmed[274], dtra2, 3);
3693 geant->Gspos(natra2_ref(0, i), 1, "ITSV", xtra1[i - 1], ytra1[i - 1], ztra1[i - 1], 0, "ONLY");
3694 }
3695
3696 atheta = 60.;
3697 aphi = 90.;
3698 xpos = (xtra1[0] + xtra1[1]) / 2.;
3699 ypos = (ytra1[0] + ytra1[1]) / 2.;
3700 zpos = 0.;
3701 geant->Gsvolu(natra3_ref(0, 1), "TUBE", idtmed[274], dtra3, 3);
3702 r2 = atheta + 90.;
3703 r3 = atheta + 90.;
3704 AliMatrix(idrotm[5200], 90., atheta, aphi + 90., r2, aphi, r3);
3705 geant->Gspos(natra3_ref(0, 1), 1, "ITSV", xpos, ypos, zpos, idrotm[5200], "ONLY");
3706
3707 atheta = 120.;
3708 aphi = 90.;
3709 xpos = (xtra1[1] + xtra1[2]) / 2.;
3710 ypos = (ytra1[1] + ytra1[2]) / 2.;
3711 zpos = 0.;
3712 geant->Gsvolu(natra3_ref(0, 2), "TUBE", idtmed[274], dtra3, 3);
3713 r2 = atheta + 90.;
3714 r3 = atheta + 90.;
3715 AliMatrix(idrotm[5201], 90., atheta, aphi + 90., r2, aphi, r3);
3716 geant->Gspos(natra3_ref(0, 2), 1, "ITSV", xpos, ypos, zpos, idrotm[5201], "ONLY");
3717
3718 atheta = 180.;
3719 aphi = 90.;
3720 xpos = (xtra1[2] + xtra1[3]) / 2.;
3721 ypos = (ytra1[2] + ytra1[3]) / 2.;
3722 zpos = 0.;
3723 geant->Gsvolu(natra3_ref(0, 3), "TUBE", idtmed[274], dtra3, 3);
3724 r2 = atheta + 90.;
3725 r3 = atheta + 90.;
3726 AliMatrix(idrotm[5202], 90., atheta, aphi + 90., r2, aphi, r3);
3727 geant->Gspos(natra3_ref(0, 3), 1, "ITSV", xpos, ypos, zpos, idrotm[5202], "ONLY");
3728
3729 atheta = 60.;
3730 aphi = 90.;
3731 xpos = (xtra1[3] + xtra1[4]) / 2.;
3732 ypos = (ytra1[3] + ytra1[4]) / 2.;
3733 zpos = 0.;
3734 geant->Gsvolu(natra3_ref(0, 4), "TUBE", idtmed[274], dtra3, 3);
3735 r2 = atheta + 90.;
3736 r3 = atheta + 90.;
3737 AliMatrix(idrotm[5203], 90., atheta, aphi + 90., r2, aphi, r3);
3738 geant->Gspos(natra3_ref(0, 4), 1, "ITSV", xpos, ypos, zpos, idrotm[5203], "ONLY");
3739
3740 atheta = 120.;
3741 aphi = 90.;
3742 xpos = (xtra1[4] + xtra1[5]) / 2.;
3743 ypos = (ytra1[4] + ytra1[5]) / 2.;
3744 zpos = 0.;
3745 geant->Gsvolu(natra3_ref(0, 5), "TUBE", idtmed[274], dtra3, 3);
3746 r2 = atheta + 90.;
3747 r3 = atheta + 90.;
3748 AliMatrix(idrotm[5204], 90., atheta, aphi + 90., r2, aphi, r3);
3749 geant->Gspos(natra3_ref(0, 5), 1, "ITSV", xpos, ypos, zpos, idrotm[5204], "ONLY");
3750
3751 atheta = 180.;
3752 aphi = 90.;
3753 xpos = (xtra1[5] + xtra1[0]) / 2.;
3754 ypos = (ytra1[5] + ytra1[0]) / 2.;
3755 zpos = 0.;
3756 geant->Gsvolu(natra3_ref(0, 6), "TUBE", idtmed[274], dtra3, 3);
3757 r2 = atheta + 90.;
3758 r3 = atheta + 90.;
3759 AliMatrix(idrotm[5205], 90., atheta, aphi + 90., r2, aphi, r3);
3760 geant->Gspos(natra3_ref(0, 6), 1, "ITSV", xpos, ypos, zpos, idrotm[5205], "ONLY");
3761
3762 atheta = 60.;
3763 aphi2 = TMath::ACos(dtra2[2] / TMath::Sqrt(dtra2[2] * dtra2[2] + (50. / cos(.59341194567807209) * (50. / cos(.59341194567807209))
3764 - 2500.))) * raddeg;
3765 aphi1 = 180. - aphi2;
3766 xpos = (xtra1[0] + xtra1[1]) / 2.;
3767 ypos = (ytra1[0] + ytra1[1]) / 2.;
3768 zpos = dtra2[2] / 2.;
3769 geant->Gsvolu(natra4_ref(0, 1), "TUBE", idtmed[274], dtra4, 3);
3770 r2 = atheta + 90.;
3771 r3 = atheta + 90.;
3772 AliMatrix(idrotm[5210], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3773 geant->Gspos(natra4_ref(0, 1), 1, "ITSV", xpos, ypos, zpos, idrotm[5210], "ONLY");
3774 zpos = -dtra2[2] / 2.;
3775 geant->Gsvolu(natra4_ref(0, 2), "TUBE", idtmed[274], dtra4, 3);
3776 r2 = atheta + 90.;
3777 r3 = atheta + 90.;
3778 AliMatrix(idrotm[5211], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3779 geant->Gspos(natra4_ref(0, 2), 1, "ITSV", xpos, ypos, zpos, idrotm[5211], "ONLY");
3780
3781 atheta = 120.;
3782 aphi1 = TMath::ACos(dtra2[2] / TMath::Sqrt(dtra2[2] * dtra2[2] + (50. / cos(.59341194567807209) * (50. / cos(.59341194567807209))
3783 - 2500.))) * raddeg;
3784 aphi2 = 180. - aphi1;
3785 xpos = (xtra1[1] + xtra1[2]) / 2.;
3786 ypos = (ytra1[1] + ytra1[2]) / 2.;
3787 zpos = dtra2[2] / 2.;
3788 geant->Gsvolu(natra4_ref(0, 3), "TUBE", idtmed[274], dtra4, 3);
3789 r2 = atheta + 90.;
3790 r3 = atheta + 90.;
3791 AliMatrix(idrotm[5212], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3792 geant->Gspos(natra4_ref(0, 3), 1, "ITSV", xpos, ypos, zpos, idrotm[5212], "ONLY");
3793 zpos = -dtra2[2] / 2.;
3794 geant->Gsvolu(natra4_ref(0, 4), "TUBE", idtmed[274], dtra4, 3);
3795 r2 = atheta + 90.;
3796 r3 = atheta + 90.;
3797 AliMatrix(idrotm[5213], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3798 geant->Gspos(natra4_ref(0, 4), 1, "ITSV", xpos, ypos, zpos, idrotm[5213], "ONLY");
3799
3800 atheta = 180.;
3801 aphi2 = TMath::ACos(dtra2[2] / TMath::Sqrt(dtra2[2] * dtra2[2] + (50. / cos(.59341194567807209) * (50. / cos(.59341194567807209))
3802 - 2500.))) * raddeg;
3803 aphi1 = 180. - aphi2;
3804 xpos = (xtra1[2] + xtra1[3]) / 2.;
3805 ypos = (ytra1[2] + ytra1[3]) / 2.;
3806 zpos = dtra2[2] / 2.;
3807 geant->Gsvolu(natra4_ref(0, 5), "TUBE", idtmed[274], dtra4, 3);
3808 r2 = atheta + 90.;
3809 r3 = atheta + 90.;
3810 AliMatrix(idrotm[5214], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3811 geant->Gspos(natra4_ref(0, 5), 1, "ITSV", xpos, ypos, zpos, idrotm[5214], "ONLY");
3812 zpos = -dtra2[2] / 2.;
3813 geant->Gsvolu(natra4_ref(0, 6), "TUBE", idtmed[274], dtra4, 3);
3814 r2 = atheta + 90.;
3815 r3 = atheta + 90.;
3816 AliMatrix(idrotm[5215], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3817 geant->Gspos(natra4_ref(0, 6), 1, "ITSV", xpos, ypos, zpos, idrotm[5215], "ONLY");
3818 atheta = 180.;
3819 aphi1 = TMath::ACos(dtra2[2] / TMath::Sqrt(dtra2[2] * dtra2[2] + (50. / cos(.59341194567807209) * (50. / cos(.59341194567807209))
3820 - 2500.))) * raddeg;
3821 aphi2 = 180. - aphi1;
3822 xpos = (xtra1[2] + xtra1[3]) / 2.;
3823 ypos = (ytra1[2] + ytra1[3]) / 2.;
3824 zpos = dtra2[2] / 2.;
3825 geant->Gsvolu(natra4_ref(0, 7), "TUBE", idtmed[274], dtra4, 3);
3826 r2 = atheta + 90.;
3827 r3 = atheta + 90.;
3828 AliMatrix(idrotm[5216], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3829 geant->Gspos(natra4_ref(0, 7), 1, "ITSV", xpos, ypos, zpos, idrotm[5216], "ONLY");
3830 zpos = -dtra2[2] / 2.;
3831 geant->Gsvolu(natra4_ref(0, 8), "TUBE", idtmed[274], dtra4, 3);
3832 r2 = atheta + 90.;
3833 r3 = atheta + 90.;
3834 AliMatrix(idrotm[5217], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3835 geant->Gspos(natra4_ref(0, 8), 1, "ITSV", xpos, ypos, zpos, idrotm[5217], "ONLY");
3836
3837 atheta = 60.;
3838 aphi2 = TMath::ACos(dtra2[2] / TMath::Sqrt(dtra2[2] * dtra2[2] + (50. / cos(.59341194567807209) * (50. / cos(.59341194567807209))
3839 - 2500.))) * raddeg;
3840 aphi1 = 180. - aphi2;
3841 xpos = (xtra1[3] + xtra1[4]) / 2.;
3842 ypos = (ytra1[3] + ytra1[4]) / 2.;
3843 zpos = dtra2[2] / 2.;
3844 geant->Gsvolu(natra4_ref(0, 9), "TUBE", idtmed[274], dtra4, 3);
3845 r2 = atheta + 90.;
3846 r3 = atheta + 90.;
3847 AliMatrix(idrotm[5218], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3848 geant->Gspos(natra4_ref(0, 9), 1, "ITSV", xpos, ypos, zpos, idrotm[5218], "ONLY");
3849 zpos = -dtra2[2] / 2.;
3850 geant->Gsvolu(natra4_ref(0, 10), "TUBE", idtmed[274], dtra4, 3);
3851 r2 = atheta + 90.;
3852 r3 = atheta + 90.;
3853 AliMatrix(idrotm[5219], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3854 geant->Gspos(natra4_ref(0, 10), 1, "ITSV", xpos, ypos, zpos, idrotm[5219], "ONLY");
3855
3856 atheta = 120.;
3857 aphi1 = TMath::ACos(dtra2[2] / TMath::Sqrt(dtra2[2] * dtra2[2] + (50. / cos(.59341194567807209) * (50. / cos(.59341194567807209))
3858 - 2500.))) * raddeg;
3859 aphi2 = 180. - aphi1;
3860 xpos = (xtra1[4] + xtra1[5]) / 2.;
3861 ypos = (ytra1[4] + ytra1[5]) / 2.;
3862 zpos = dtra2[2] / 2.;
3863 geant->Gsvolu(natra4_ref(0, 11), "TUBE", idtmed[274], dtra4, 3);
3864 r2 = atheta + 90.;
3865 r3 = atheta + 90.;
3866 AliMatrix(idrotm[5220], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3867 geant->Gspos(natra4_ref(0, 11), 1, "ITSV", xpos, ypos, zpos, idrotm[5220], "ONLY");
3868 zpos = -dtra2[2] / 2.;
3869 geant->Gsvolu(natra4_ref(0, 12), "TUBE", idtmed[274], dtra4, 3);
3870 r2 = atheta + 90.;
3871 r3 = atheta + 90.;
3872 AliMatrix(idrotm[5221], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3873 geant->Gspos(natra4_ref(0, 12), 1, "ITSV", xpos, ypos, zpos, idrotm[5221], "ONLY");
3874
3875 atheta = 180.;
3876 aphi2 = TMath::ACos(dtra2[2] / TMath::Sqrt(dtra2[2] * dtra2[2] + (50. / cos(.59341194567807209) * (50. / cos(.59341194567807209))
3877 - 2500.))) * raddeg;
3878 aphi1 = 180. - aphi2;
3879 xpos = (xtra1[5] + xtra1[0]) / 2.;
3880 ypos = (ytra1[5] + ytra1[0]) / 2.;
3881 zpos = dtra2[2] / 2.;
3882 geant->Gsvolu(natra4_ref(0, 13), "TUBE", idtmed[274], dtra4, 3);
3883 r2 = atheta + 90.;
3884 r3 = atheta + 90.;
3885 AliMatrix(idrotm[5222], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3886 geant->Gspos(natra4_ref(0, 13), 1, "ITSV", xpos, ypos, zpos, idrotm[5222], "ONLY");
3887 zpos = -dtra2[2] / 2.;
3888 geant->Gsvolu(natra4_ref(0, 14), "TUBE", idtmed[274], dtra4, 3);
3889 r2 = atheta + 90.;
3890 r3 = atheta + 90.;
3891 AliMatrix(idrotm[5223], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3892 geant->Gspos(natra4_ref(0, 14), 1, "ITSV", xpos, ypos, zpos, idrotm[5223], "ONLY");
3893 atheta = 180.;
3894 aphi1 = TMath::ACos(dtra2[2] / TMath::Sqrt(dtra2[2] * dtra2[2] + (50. / cos(.59341194567807209) * (50. / cos(.59341194567807209))
3895 - 2500.))) * raddeg;
3896 aphi2 = 180. - aphi1;
3897 xpos = (xtra1[5] + xtra1[0]) / 2.;
3898 ypos = (ytra1[5] + ytra1[0]) / 2.;
3899 zpos = dtra2[2] / 2.;
3900 geant->Gsvolu(natra4_ref(0, 15), "TUBE", idtmed[274], dtra4, 3);
3901 r2 = atheta + 90.;
3902 r3 = atheta + 90.;
3903 AliMatrix(idrotm[5224], 90., atheta, aphi1 + 90., r2, aphi1, r3);
3904 geant->Gspos(natra4_ref(0, 15), 1, "ITSV", xpos, ypos, zpos, idrotm[5224], "ONLY");
3905 zpos = -dtra2[2] / 2.;
3906 geant->Gsvolu(natra4_ref(0, 16), "TUBE", idtmed[274], dtra4, 3);
3907 r2 = atheta + 90.;
3908 r3 = atheta + 90.;
3909 AliMatrix(idrotm[5225], 90., atheta, aphi2 + 90., r2, aphi2, r3);
3910 geant->Gspos(natra4_ref(0, 16), 1, "ITSV", xpos, ypos, zpos, idrotm[5225], "ONLY");
3911
3912
3913 } else if (fMinorVersion == 5) {
3914
3915
3916 rzero = 33.5;
3917 dtra2[0] = .92;
3918 dtra2[1] = 1.;
3919 dtra2[2] = dpcb[2] * 2. + 50. - 10.5;
3920 dtra3[0] = .92;
3921 dtra3[1] = 1.;
3922 dtra3[2] = 16.75;
3923 dtra4[0] = .92;
3924 dtra4[1] = 1.;
3925 dtra4[2] = TMath::Sqrt(dtra2[2] * dtra2[2] + 1088.0099999999998) / 2.;
3926 angle = 60.;
3927 offset = angle / 2.;
3928 for (i = 1; i <= 6; ++i) {
3929 xtra1[i - 1] = rzero * TMath::Cos(((i-1) * angle + offset) *degrad);
3930 ytra1[i - 1] = rzero * TMath::Sin(((i-1) * angle + offset) *degrad);
3931 ztra1[i - 1] = 0.;
3932 geant->Gsvolu(natra2_ref(0, i), "TUBE", idtmed[274], dtra2, 3);
3933 geant->Gspos(natra2_ref(0, i), 1, "ITSV", xtra1[i - 1], ytra1[i - 1], ztra1[i - 1], 0, "ONLY");
3934 }
3935
3936 atheta = 60.;
3937 aphi = 90.;
3938 xpos = (xtra1[0] + xtra1[1]) / 2.;
3939 ypos = (ytra1[0] + ytra1[1]) / 2.;
3940 zpos = 0.;
3941 geant->Gsvolu(natra3_ref(0, 1), "TUBE", idtmed[274], dtra3, 3);
3942 r2 = atheta + 90.;
3943 r3 = atheta + 90.;
3944 AliMatrix(idrotm[5200], 90., atheta, aphi + 90., r2, aphi, r3);
3945 geant->Gspos(natra3_ref(0, 1), 1, "ITSV", xpos, ypos, zpos, idrotm[5200], "ONLY");
3946
3947 atheta = 120.;
3948 aphi = 90.;
3949 xpos = (xtra1[1] + xtra1[2]) / 2.;
3950 ypos = (ytra1[1] + ytra1[2]) / 2.;
3951 zpos = 0.;
3952 geant->Gsvolu(natra3_ref(0, 2), "TUBE", idtmed[274], dtra3, 3);
3953 r2 = atheta + 90.;
3954 r3 = atheta + 90.;
3955 AliMatrix(idrotm[5201], 90., atheta, aphi + 90., r2, aphi, r3);
3956 geant->Gspos(natra3_ref(0, 2), 1, "ITSV", xpos, ypos, zpos, idrotm[5201], "ONLY");
3957
3958 atheta = 180.;
3959 aphi = 90.;
3960 xpos = (xtra1[2] + xtra1[3]) / 2.;
3961 ypos = (ytra1[2] + ytra1[3]) / 2.;
3962 zpos = 0.;
3963 geant->Gsvolu(natra3_ref(0, 3), "TUBE", idtmed[274], dtra3, 3);
3964 r2 = atheta + 90.;
3965 r3 = atheta + 90.;
3966 AliMatrix(idrotm[5202], 90., atheta, aphi + 90., r2, aphi, r3);
3967 geant->Gspos(natra3_ref(0, 3), 1, "ITSV", xpos, ypos, zpos, idrotm[5202], "ONLY");
3968
3969 atheta = 60.;
3970 aphi = 90.;
3971 xpos = (xtra1[3] + xtra1[4]) / 2.;
3972 ypos = (ytra1[3] + ytra1[4]) / 2.;
3973 zpos = 0.;
3974 geant->Gsvolu(natra3_ref(0, 4), "TUBE", idtmed[274], dtra3, 3);
3975 r2 = atheta + 90.;
3976 r3 = atheta + 90.;
3977 AliMatrix(idrotm[5203], 90., atheta, aphi + 90., r2, aphi, r3);
3978 geant->Gspos(natra3_ref(0, 4), 1, "ITSV", xpos, ypos, zpos, idrotm[5203], "ONLY");
3979
3980 atheta = 120.;
3981 aphi = 90.;
3982 xpos = (xtra1[4] + xtra1[5]) / 2.;
3983 ypos = (ytra1[4] + ytra1[5]) / 2.;
3984 zpos = 0.;
3985 geant->Gsvolu(natra3_ref(0, 5), "TUBE", idtmed[274], dtra3, 3);
3986 r2 = atheta + 90.;
3987 r3 = atheta + 90.;
3988 AliMatrix(idrotm[5204], 90., atheta, aphi + 90., r2, aphi, r3);
3989 geant->Gspos(natra3_ref(0, 5), 1, "ITSV", xpos, ypos, zpos, idrotm[5204], "ONLY");
3990
3991 atheta = 180.;
3992 aphi = 90.;
3993 xpos = (xtra1[5] + xtra1[0]) / 2.;
3994 ypos = (ytra1[5] + ytra1[0]) / 2.;
3995 zpos = 0.;
3996 geant->Gsvolu(natra3_ref(0, 6), "TUBE", idtmed[274], dtra3, 3);
3997 r2 = atheta + 90.;
3998 r3 = atheta + 90.;
3999 AliMatrix(idrotm[5205], 90., atheta, aphi + 90., r2, aphi, r3);
4000 geant->Gspos(natra3_ref(0, 6), 1, "ITSV", xpos, ypos, zpos, idrotm[5205], "ONLY");
4001
4002 atheta = 60.;
4003 aphi2 = TMath::ACos(dtra2[2] / TMath::Sqrt(dtra2[2] * dtra2[2] + (50. / cos(.59341194567807209) * (50. / cos(.59341194567807209))
4004 - 2500.))) * raddeg;
4005 aphi1 = 180. - aphi2;
4006 xpos = (xtra1[0] + xtra1[1]) / 2.;
4007 ypos = (ytra1[0] + ytra1[1]) / 2.;
4008 zpos = dtra2[2] / 2.;
4009 geant->Gsvolu(natra4_ref(0, 1), "TUBE", idtmed[274], dtra4, 3);
4010 r2 = atheta + 90.;
4011 r3 = atheta + 90.;
4012 AliMatrix(idrotm[5210], 90., atheta, aphi1 + 90., r2, aphi1, r3);
4013 geant->Gspos(natra4_ref(0, 1), 1, "ITSV", xpos, ypos, zpos, idrotm[5210], "ONLY");
4014 zpos = -dtra2[2] / 2.;
4015 geant->Gsvolu(natra4_ref(0, 2), "TUBE", idtmed[274], dtra4, 3);
4016 r2 = atheta + 90.;
4017 r3 = atheta + 90.;
4018 AliMatrix(idrotm[5211], 90., atheta, aphi2 + 90., r2, aphi2, r3);
4019 geant->Gspos(natra4_ref(0, 2), 1, "ITSV", xpos, ypos, zpos, idrotm[5211], "ONLY");
4020
4021 atheta = 120.;
4022 aphi1 = TMath::ACos(dtra2[2] / TMath::Sqrt(dtra2[2] * dtra2[2] + (50. / cos(.59341194567807209) * (50. / cos(.59341194567807209))
4023 - 2500.))) * raddeg;
4024 aphi2 = 180. - aphi1;
4025 xpos = (xtra1[1] + xtra1[2]) / 2.;
4026 ypos = (ytra1[1] + ytra1[2]) / 2.;
4027 zpos = dtra2[2] / 2.;
4028 geant->Gsvolu(natra4_ref(0, 3), "TUBE", idtmed[274], dtra4, 3);
4029 r2 = atheta + 90.;
4030 r3 = atheta + 90.;
4031 AliMatrix(idrotm[5212], 90., atheta, aphi1 + 90., r2, aphi1, r3);
4032 geant->Gspos(natra4_ref(0, 3), 1, "ITSV", xpos, ypos, zpos, idrotm[5212], "ONLY");
4033 zpos = -dtra2[2] / 2.;
4034 geant->Gsvolu(natra4_ref(0, 4), "TUBE", idtmed[274], dtra4, 3);
4035 r2 = atheta + 90.;
4036 r3 = atheta + 90.;
4037 AliMatrix(idrotm[5213], 90., atheta, aphi2 + 90., r2, aphi2, r3);
4038 geant->Gspos(natra4_ref(0, 4), 1, "ITSV", xpos, ypos, zpos, idrotm[5213], "ONLY");
4039
4040 atheta = 180.;
4041 aphi2 = TMath::ACos(dtra2[2] / TMath::Sqrt(dtra2[2] * dtra2[2] + (50. / cos(.59341194567807209) * (50. / cos(.59341194567807209))
4042 - 2500.))) * raddeg;
4043 aphi1 = 180. - aphi2;
4044 xpos = (xtra1[2] + xtra1[3]) / 2.;
4045 ypos = (ytra1[2] + ytra1[3]) / 2.;
4046 zpos = dtra2[2] / 2.;
4047 geant->Gsvolu(natra4_ref(0, 5), "TUBE", idtmed[274], dtra4, 3);
4048 r2 = atheta + 90.;
4049 r3 = atheta + 90.;
4050 AliMatrix(idrotm[5214], 90., atheta, aphi1 + 90., r2, aphi1, r3);
4051 geant->Gspos(natra4_ref(0, 5), 1, "ITSV", xpos, ypos, zpos, idrotm[5214], "ONLY");
4052 zpos = -dtra2[2] / 2.;
4053 geant->Gsvolu(natra4_ref(0, 6), "TUBE", idtmed[274], dtra4, 3);
4054 r2 = atheta + 90.;
4055 r3 = atheta + 90.;
4056 AliMatrix(idrotm[5215], 90., atheta, aphi2 + 90., r2, aphi2, r3);
4057 geant->Gspos(natra4_ref(0, 6), 1, "ITSV", xpos, ypos, zpos, idrotm[5215], "ONLY");
4058 atheta = 180.;
4059 aphi1 = TMath::ACos(dtra2[2] / TMath::Sqrt(dtra2[2] * dtra2[2] + (50. / cos(.59341194567807209) * (50. / cos(.59341194567807209))
4060 - 2500.))) * raddeg;
4061 aphi2 = 180. - aphi1;
4062 xpos = (xtra1[2] + xtra1[3]) / 2.;
4063 ypos = (ytra1[2] + ytra1[3]) / 2.;
4064 zpos = dtra2[2] / 2.;
4065 geant->Gsvolu(natra4_ref(0, 7), "TUBE", idtmed[274], dtra4, 3);
4066 r2 = atheta + 90.;
4067 r3 = atheta + 90.;
4068 AliMatrix(idrotm[5216], 90., atheta, aphi1 + 90., r2, aphi1, r3);
4069 geant->Gspos(natra4_ref(0, 7), 1, "ITSV", xpos, ypos, zpos, idrotm[5216], "ONLY");
4070 zpos = -dtra2[2] / 2.;
4071 geant->Gsvolu(natra4_ref(0, 8), "TUBE", idtmed[274], dtra4, 3);
4072 r2 = atheta + 90.;
4073 r3 = atheta + 90.;
4074 AliMatrix(idrotm[5217], 90., atheta, aphi2 + 90., r2, aphi2, r3);
4075 geant->Gspos(natra4_ref(0, 8), 1, "ITSV", xpos, ypos, zpos, idrotm[5217], "ONLY");
4076
4077 atheta = 60.;
4078 aphi2 = TMath::ACos(dtra2[2] / TMath::Sqrt(dtra2[2] * dtra2[2] + (50. / cos(.59341194567807209) * (50. / cos(.59341194567807209))
4079 - 2500.))) * raddeg;
4080 aphi1 = 180. - aphi2;
4081 xpos = (xtra1[3] + xtra1[4]) / 2.;
4082 ypos = (ytra1[3] + ytra1[4]) / 2.;
4083 zpos = dtra2[2] / 2.;
4084 geant->Gsvolu(natra4_ref(0, 9), "TUBE", idtmed[274], dtra4, 3);
4085 r2 = atheta + 90.;
4086 r3 = atheta + 90.;
4087 AliMatrix(idrotm[5218], 90., atheta, aphi1 + 90., r2, aphi1, r3);
4088 geant->Gspos(natra4_ref(0, 9), 1, "ITSV", xpos, ypos, zpos, idrotm[5218], "ONLY");
4089 zpos = -dtra2[2] / 2.;
4090 geant->Gsvolu(natra4_ref(0, 10), "TUBE", idtmed[274], dtra4, 3);
4091 r2 = atheta + 90.;
4092 r3 = atheta + 90.;
4093 AliMatrix(idrotm[5219], 90., atheta, aphi2 + 90., r2, aphi2, r3);
4094 geant->Gspos(natra4_ref(0, 10), 1, "ITSV", xpos, ypos, zpos, idrotm[5219], "ONLY");
4095
4096 atheta = 120.;
4097 aphi1 = TMath::ACos(dtra2[2] / TMath::Sqrt(dtra2[2] * dtra2[2] + (50. / cos(.59341194567807209) * (50. / cos(.59341194567807209))
4098 - 2500.))) * raddeg;
4099 aphi2 = 180. - aphi1;
4100 xpos = (xtra1[4] + xtra1[5]) / 2.;
4101 ypos = (ytra1[4] + ytra1[5]) / 2.;
4102 zpos = dtra2[2] / 2.;
4103 geant->Gsvolu(natra4_ref(0, 11), "TUBE", idtmed[274], dtra4, 3);
4104 r2 = atheta + 90.;
4105 r3 = atheta + 90.;
4106 AliMatrix(idrotm[5220], 90., atheta, aphi1 + 90., r2, aphi1, r3);
4107 geant->Gspos(natra4_ref(0, 11), 1, "ITSV", xpos, ypos, zpos, idrotm[5220], "ONLY");
4108 zpos = -dtra2[2] / 2.;
4109 geant->Gsvolu(natra4_ref(0, 12), "TUBE", idtmed[274], dtra4, 3);
4110 r2 = atheta + 90.;
4111 r3 = atheta + 90.;
4112 AliMatrix(idrotm[5221], 90., atheta, aphi2 + 90., r2, aphi2, r3);
4113 geant->Gspos(natra4_ref(0, 12), 1, "ITSV", xpos, ypos, zpos, idrotm[5221], "ONLY");
4114
4115 atheta = 180.;
4116 aphi2 = TMath::ACos(dtra2[2] / TMath::Sqrt(dtra2[2] * dtra2[2] + (50. / cos(.59341194567807209) * (50. / cos(.59341194567807209))
4117 - 2500.))) * raddeg;
4118 aphi1 = 180. - aphi2;
4119 xpos = (xtra1[5] + xtra1[0]) / 2.;
4120 ypos = (ytra1[5] + ytra1[0]) / 2.;
4121 zpos = dtra2[2] / 2.;
4122 geant->Gsvolu(natra4_ref(0, 13), "TUBE", idtmed[274], dtra4, 3);
4123 r2 = atheta + 90.;
4124 r3 = atheta + 90.;
4125 AliMatrix(idrotm[5222], 90., atheta, aphi1 + 90., r2, aphi1, r3);
4126 geant->Gspos(natra4_ref(0, 13), 1, "ITSV", xpos, ypos, zpos, idrotm[5222], "ONLY");
4127 zpos = -dtra2[2] / 2.;
4128 geant->Gsvolu(natra4_ref(0, 14), "TUBE", idtmed[274], dtra4, 3);
4129 r2 = atheta + 90.;
4130 r3 = atheta + 90.;
4131 AliMatrix(idrotm[5223], 90., atheta, aphi2 + 90., r2, aphi2, r3);
4132 geant->Gspos(natra4_ref(0, 14), 1, "ITSV", xpos, ypos, zpos, idrotm[5223], "ONLY");
4133 atheta = 180.;
4134 aphi1 = TMath::ACos(dtra2[2] / TMath::Sqrt(dtra2[2] * dtra2[2] + (50. / cos(.59341194567807209) * (50. / cos(.59341194567807209))
4135 - 2500.))) * raddeg;
4136 aphi2 = 180. - aphi1;
4137 xpos = (xtra1[5] + xtra1[0]) / 2.;
4138 ypos = (ytra1[5] + ytra1[0]) / 2.;
4139 zpos = dtra2[2] / 2.;
4140 geant->Gsvolu(natra4_ref(0, 15), "TUBE", idtmed[274], dtra4, 3);
4141 r2 = atheta + 90.;
4142 r3 = atheta + 90.;
4143 AliMatrix(idrotm[5224], 90., atheta, aphi1 + 90., r2, aphi1, r3);
4144 geant->Gspos(natra4_ref(0, 15), 1, "ITSV", xpos, ypos, zpos, idrotm[5224], "ONLY");
4145 zpos = -dtra2[2] / 2.;
4146 geant->Gsvolu(natra4_ref(0, 16), "TUBE", idtmed[274], dtra4, 3);
4147 r2 = atheta + 90.;
4148 r3 = atheta + 90.;
4149 AliMatrix(idrotm[5225], 90., atheta, aphi2 + 90., r2, aphi2, r3);
4150 geant->Gspos(natra4_ref(0, 16), 1, "ITSV", xpos, ypos, zpos, idrotm[5225], "ONLY");
4151 } else {
4152 goto L9123;
4153 }
4154
4155L9123:
4156// --- Define the end-caps
4157
4158// GOTO 9234 ! skip both end-caps
4159
4160// --- Define the Z>0 end-cap
4161
4162// GOTO 9345 ! skip the Z>0 end-cap
4163
4164 dcone[0] = 16.75;
4165 dcone[1] = 12.;
4166 dcone[2] = 12.02;
4167 dcone[3] = 46.899999999999999;
4168 dcone[4] = .02 / TMath::Cos(.78539816339744828) + 46.899999999999999;
4169 xpos = 0.;
4170 ypos = 0.;
4171 zpos = dpcb[2] * 2. + 45.899999999999999 - 10.5;
4172// end-ladder electro
4173 geant->Gsvolu("RCON", "CONE", idtmed[274], dcone, 5);
4174 ;
4175 geant->Gspos("RCON", 1, "ITSV", xpos, ypos, zpos, 0, "ONLY");
4176
4177 dtube[0] = .02 / TMath::Cos(.78539816339744828) +
4178 46.899999999999999;
4179 dtube[1] = 49.9;
4180// In the Simonetti's drawings 52. In the TP 50.
4181 dtube[2] = .15;
4182 xpos = 0.;
4183 ypos = 0.;
4184 zpos = dpcb[2] * 2. + 62.799999999999997 - 10.5;
4185// end-ladder electro
4186 geant->Gsvolu("RTB1", "TUBE", idtmed[274], dtube, 3);
4187 ;
4188 geant->Gspos("RTB1", 1, "ITSV", xpos, ypos, zpos, 0, "ONLY");
4189
4190 dtube[0] = 10.5;
4191 dtube[1] = 12.;
4192 dtube[2] = 1.3400000000000001;
4193 xpos = 0.;
4194 ypos = 0.;
4195 zpos = dpcb[2] * 2. + 21.59 - 10.5;
4196// end-ladder elect
4197 geant->Gsvolu("RTB2", "TUBE", idtmed[274], dtube, 3);
4198
4199 geant->Gspos("RTB2", 1, "ITSV", xpos, ypos, zpos, 0, "ONLY");
4200
4201 dpgon[0] = 15.;
4202 dpgon[1] = 360.;
4203 dpgon[2] = 12.;
4204 dpgon[3] = 2.;
4205 dpgon[4] = dpcb[2] * 2. + 22.93 - 10.5;
4206// end-ladder electron
4207 dpgon[5] = 12.;
4208 dpgon[6] = 13.5;
4209 dpgon[7] = dpcb[2] * 2. + 29.149999999999999 - 10.5;
4210// end-ladder electronics
4211 dpgon[8] = 12.;
4212 dpgon[9] = 13.5;
4213 xpos = 0.;
4214 ypos = 0.;
4215 zpos = 0.;
4216 geant->Gsvolu("RP03", "PGON", idtmed[274], dpgon, 10);
4217 geant->Gspos("RP03", 1, "ITSV", xpos, ypos, zpos, 0, "ONLY");
4218
4219 dpgon[0] = 7.5;
4220 dpgon[1] = 360.;
4221 dpgon[2] = 24.;
4222 dpgon[3] = 2.;
4223 dpgon[4] = dpcb[2] * 2. + 35.649999999999999 - 10.5;
4224// end-ladder e
4225 dpgon[5] = 21.;
4226 dpgon[6] = 23.;
4227 dpgon[7] = dpcb[2] * 2. + 37.149999999999999 - 10.5;
4228// end-ladde
4229 dpgon[8] = 21.;
4230 dpgon[9] = 23.;
4231 xpos = 0.;
4232 ypos = 0.;
4233 zpos = 0.;
4234 geant->Gsvolu("RP04", "PGON", idtmed[274], dpgon, 10);
4235 geant->Gspos("RP04", 1, "ITSV", xpos, ypos, zpos, 0, "ONLY");
4236
4237 if (fMinorVersion == 0 || fMinorVersion == 1 || fMinorVersion ==2) {
4238 offset2 = 5.2;
4239 dpgon[0] = offset2 + 5.1428571428571432;
4240 dpgon[1] = 360.;
4241 dpgon[2] = 35.;
4242 dpgon[3] = 2.;
4243 dpgon[4] = dpcb[2] * 2. + 52.350000000000001 - 10.5;
4244// end-ladde
4245 dpgon[5] = 37.7;
4246 dpgon[6] = 40.;
4247 dpgon[7] = dpcb[2] * 2. + 53.850000000000001 - 10.5;
4248// end-la
4249 dpgon[8] = 37.7;
4250 dpgon[9] = 40.;
4251 xpos = 0.;
4252 ypos = 0.;
4253 zpos = 0.;
4254 geant->Gsvolu("RP05", "PGON", idtmed[274], dpgon, 10);
4255 geant->Gspos("RP05", 1, "ITSV", xpos, ypos, zpos, 0, "ONLY");
4256
4257 dpgon[0] = offset2 + 4.615384615384615;
4258 dpgon[1] = 360.;
4259 dpgon[2] = 39.;
4260 dpgon[3] = 2.;
4261 dpgon[4] = dpcb[2] * 2. + 57.350000000000001 - 10.5;
4262// end-ladder
4263 dpgon[5] = 42.7;
4264 dpgon[6] = 45.;
4265 dpgon[7] = dpcb[2] * 2. + 58.850000000000001 - 10.5;
4266// end-la
4267 dpgon[8] = 42.7;
4268 dpgon[9] = 45.;
4269 xpos = 0.;
4270 ypos = 0.;
4271 zpos = 0.;
4272 geant->Gsvolu("RP06", "PGON", idtmed[274], dpgon, 10);
4273 geant->Gspos("RP06", 1, "ITSV", xpos, ypos, zpos, 0, "ONLY");
4274 }
4275 if (fMinorVersion == 3 || fMinorVersion == 4 || fMinorVersion == 5) {
4276 offset2 = 5.2;
4277 dpgon[0] = offset2 + 5.625;
4278 dpgon[1] = 360.;
4279 dpgon[2] = 32.;
4280 dpgon[3] = 2.;
4281 dpgon[4] = 52.350000000000001 - 3.3999999999999986 / TMath::Tan(.78539816339744828) + dpcb[2] * 2. - 10.5;
4282// end-ladder electronics
4283 dpgon[5] = 34.3;
4284 dpgon[6] = 36.6;
4285 dpgon[7] = 53.850000000000001 - 3.3999999999999986 / TMath::Tan(.78539816339744828) + dpcb[2] * 2. - 10.5;
4286// end-ladder electr
4287 dpgon[8] = 34.3;
4288 dpgon[9] = 36.6;
4289 xpos = 0.;
4290 ypos = 0.;
4291 zpos = 0.;
4292 geant->Gsvolu("RP05", "PGON", idtmed[274], dpgon, 10);
4293 geant->Gspos("RP05", 1, "ITSV", xpos, ypos, zpos, 0, "ONLY");
4294
4295 dpgon[0] = offset2 + 5.;
4296 dpgon[1] = 360.;
4297 dpgon[2] = 36.;
4298 dpgon[3] = 2.;
4299 dpgon[4] = 57.350000000000001 - 3.7999999999999972 / TMath::Tan(.78539816339744828) + dpcb[2] * 2. - 10.5;
4300// end-ladder electronics
4301 dpgon[5] = 38.9;
4302 dpgon[6] = 41.2;
4303 dpgon[7] = 58.850000000000001 - 3.7999999999999972 / TMath::Tan(.78539816339744828) + dpcb[2] * 2. - 10.5;
4304// end-ladder electr
4305 dpgon[8] = 38.9;
4306 dpgon[9] = 41.2;
4307 xpos = 0.;
4308 ypos = 0.;
4309 zpos = 0.;
4310 geant->Gsvolu("RP06", "PGON", idtmed[274], dpgon, 10);
4311 geant->Gspos("RP06", 1, "ITSV", xpos, ypos, zpos, 0, "ONLY");
4312 }
4313
4314// 9345 CONTINUE
4315
4316// --- Define the Z<0 end-cap
4317
4318// GOTO 9456 ! skip the Z<0 end-cap
4319
4320 dcone[0] = 16.75;
4321 dcone[1] = 46.899999999999999;
4322 dcone[2] = .02 / TMath::Cos(.78539816339744828) + 46.899999999999999;
4323 dcone[3] = 12.;
4324 dcone[4] = 12.02;
4325 xpos = 0.;
4326 ypos = 0.;
4327 zpos = -45.899999999999999 - dpcb[2] * 2. + 10.5;
4328// end-ladder electr
4329 geant->Gsvolu("LCON", "CONE", idtmed[274], dcone, 5);
4330 ;
4331 geant->Gspos("LCON", 1, "ITSV", xpos, ypos, zpos, 0, "ONLY");
4332
4333 dtube[0] = .02 / TMath::Cos(.78539816339744828) +
4334 46.899999999999999;
4335 dtube[1] = 49.9;
4336// In the Simonetti's drawings 52. In the TP 50.
4337 dtube[2] = .15;
4338 xpos = 0.;
4339 ypos = 0.;
4340 zpos = -62.799999999999997 - dpcb[2] * 2. + 10.5;
4341// end-ladder electr
4342 geant->Gsvolu("LTB1", "TUBE", idtmed[274], dtube, 3);
4343 ;
4344 geant->Gspos("LTB1", 1, "ITSV", xpos, ypos, zpos, 0, "ONLY");
4345
4346 dtube[0] = 10.5;
4347 dtube[1] = 12.;
4348 dtube[2] = 1.3400000000000001;
4349 xpos = 0.;
4350 ypos = 0.;
4351 zpos = -21.59 - dpcb[2] * 2. + 10.5;
4352// end-ladder elec
4353 geant->Gsvolu("LTB2", "TUBE", idtmed[274], dtube, 3);
4354
4355 geant->Gspos("LTB2", 1, "ITSV", xpos, ypos, zpos, 0, "ONLY");
4356
4357 dpgon[0] = 15.;
4358 dpgon[1] = 360.;
4359 dpgon[2] = 12.;
4360 dpgon[3] = 2.;
4361 dpgon[4] = -29.149999999999999 - dpcb[2] * 2. + 10.5;
4362// end-ladder electronics
4363 dpgon[5] = 12.;
4364 dpgon[6] = 13.5;
4365 dpgon[7] = -22.93 - dpcb[2] * 2. + 10.5;
4366// end-ladder electro
4367 dpgon[8] = 12.;
4368 dpgon[9] = 13.5;
4369 xpos = 0.;
4370 ypos = 0.;
4371 zpos = 0.;
4372 geant->Gsvolu("LP03", "PGON", idtmed[274], dpgon, 10);
4373 geant->Gspos("LP03", 1, "ITSV", xpos, ypos, zpos, 0, "ONLY");
4374
4375 dpgon[0] = 7.5;
4376 dpgon[1] = 360.;
4377 dpgon[2] = 24.;
4378 dpgon[3] = 2.;
4379 dpgon[4] = -37.149999999999999 - dpcb[2] * 2. + 10.5;
4380// end-ladd
4381 dpgon[5] = 21.;
4382 dpgon[6] = 23.;
4383 dpgon[7] = -35.649999999999999 - dpcb[2] * 2. + 10.5;
4384// end-ladder
4385 dpgon[8] = 21.;
4386 dpgon[9] = 23.;
4387 xpos = 0.;
4388 ypos = 0.;
4389 zpos = 0.;
4390 geant->Gsvolu("LP04", "PGON", idtmed[274], dpgon, 10);
4391 geant->Gspos("LP04", 1, "ITSV", xpos, ypos, zpos, 0, "ONLY");
4392
4393 if (fMinorVersion == 0 || fMinorVersion == 1 || fMinorVersion ==2) {
4394 offset2 = 5.2;
4395 dpgon[0] = offset2 + 5.1428571428571432;
4396 dpgon[1] = 360.;
4397 dpgon[2] = 35.;
4398 dpgon[3] = 2.;
4399 dpgon[4] = -52.350000000000001 - dpcb[2] * 2. + 10.5;
4400// end-ladd
4401 dpgon[5] = 37.7;
4402 dpgon[6] = 40.;
4403 dpgon[7] = -53.850000000000001 - dpcb[2] * 2. + 10.5;
4404// end-l
4405 dpgon[8] = 37.7;
4406 dpgon[9] = 40.;
4407 xpos = 0.;
4408 ypos = 0.;
4409 zpos = 0.;
4410 geant->Gsvolu("LP05", "PGON", idtmed[274], dpgon, 10);
4411 geant->Gspos("LP05", 1, "ITSV", xpos, ypos, zpos, 0, "ONLY");
4412
4413 dpgon[0] = offset2 + 4.615384615384615;
4414 dpgon[1] = 360.;
4415 dpgon[2] = 39.;
4416 dpgon[3] = 2.;
4417 dpgon[4] = -57.350000000000001 - dpcb[2] * 2. + 10.5;
4418// end-ladde
4419 dpgon[5] = 42.7;
4420 dpgon[6] = 45.;
4421 dpgon[7] = -58.850000000000001 - dpcb[2] * 2. + 10.5;
4422// end-l
4423 dpgon[8] = 42.7;
4424 dpgon[9] = 45.;
4425 xpos = 0.;
4426 ypos = 0.;
4427 zpos = 0.;
4428 geant->Gsvolu("LP06", "PGON", idtmed[274], dpgon, 10);
4429 geant->Gspos("LP06", 1, "ITSV", xpos, ypos, zpos, 0, "ONLY");
4430 }
4431 if (fMinorVersion == 3 || fMinorVersion == 4 || fMinorVersion == 5) {
4432 offset2 = 5.2;
4433 dpgon[0] = offset2 + 5.625;
4434 dpgon[1] = 360.;
4435 dpgon[2] = 32.;
4436 dpgon[3] = 2.;
4437 dpgon[4] = 3.4 / TMath::Tan(.78539816339744828) - 52.35 - dpcb[2] * 2. + 10.5;
4438// end-ladder electronics
4439 dpgon[5] = 34.3;
4440 dpgon[6] = 36.6;
4441 dpgon[7] = 3.4 / TMath::Tan(.78539816339744828) - 53.85 - dpcb[2] * 2. + 10.5;
4442// end-ladder electr
4443 dpgon[8] = 34.3;
4444 dpgon[9] = 36.6;
4445 xpos = 0.;
4446 ypos = 0.;
4447 zpos = 0.;
4448 geant->Gsvolu("LP05", "PGON", idtmed[274], dpgon, 10);
4449 geant->Gspos("LP05", 1, "ITSV", xpos, ypos, zpos, 0, "ONLY");
4450
4451 dpgon[0] = offset2 + 5.;
4452 dpgon[1] = 360.;
4453 dpgon[2] = 36.;
4454 dpgon[3] = 2.;
4455 dpgon[4] = 3.8 / TMath::Tan(.78539816339744828) - 57.35 - dpcb[2] * 2. + 10.5;
4456// end-ladder electronics
4457 dpgon[5] = 38.9;
4458 dpgon[6] = 41.2;
4459 dpgon[7] = 3.8 / TMath::Tan(.78539816339744828) - 58.85 - dpcb[2] * 2. + 10.5;
4460// end-ladder electr
4461 dpgon[8] = 38.9;
4462 dpgon[9] = 41.2;
4463 xpos = 0.;
4464 ypos = 0.;
4465 zpos = 0.;
4466 geant->Gsvolu("LP06", "PGON", idtmed[274], dpgon, 10);
4467 geant->Gspos("LP06", 1, "ITSV", xpos, ypos, zpos, 0, "ONLY");
4468 }
4469
4470// 9456 CONTINUE
4471
4472
4473// --- Outputs the geometry tree in the EUCLID/CAD format
4474
4475 if (fEuclidOut) {
4476 geant->WriteEuclid("ITSgeometry", "ITSV", 1, 5);
4477 }
4478#endif
4479// --- Outputs the geometry tree in the EUCLID/CAD format
4480
4481 if (fEuclidOut) {
4482 pMC->WriteEuclid("ITSgeometry", "ITSV", 1, 5);
4483 }
4484}
4485
4486//_____________________________________________________________________________
4487void AliITSv4::CreateMaterials()
4488{
4489 //
4490 // Create Materials for ITS
4491 //
4492 AliITS::CreateMaterials();
4493}
4494
4495//_____________________________________________________________________________
4496void AliITSv4::Init()
4497{
4498 //
4499 // Initialise ITS after it has been created
4500 //
4501 AliMC* pMC = AliMC::GetMC();
4502 fIdSens1=pMC->VolId("ITS1");
4503 fIdSens2=pMC->VolId("ITS2");
4504 fIdSens3=pMC->VolId("ITS3");
4505 fIdSens4=pMC->VolId("ITS4");
4506 fIdSens5=pMC->VolId("ITS5");
4507 fIdSens6=pMC->VolId("ITS6");
4508}
4509
4510//_____________________________________________________________________________
4511void AliITSv4::StepManager()
4512{
4513 //
4514 // Called at every step in the ITS
4515 //
4516 Int_t copy, id;
4517 Float_t hits[7];
4518 Int_t vol[3];
4519 Float_t position[3];
4520 Float_t momentum[4];
4521 TClonesArray &lhits = *fHits;
4522 AliMC* pMC = AliMC::GetMC();
4523 //
4524 if(pMC->TrackCharge() && pMC->Edep()) {
4525 //
4526 // Only entering charged tracks
4527 if((id=pMC->CurrentVol(0,copy))==fIdSens1) {
4528 vol[0]=1;
4529 id=pMC->CurrentVolOff(1,0,copy);
4530 vol[1]=copy;
4531 id=pMC->CurrentVolOff(2,0,copy);
4532 vol[2]=copy;
4533 } else if(id==fIdSens2) {
4534 vol[0]=2;
4535 id=pMC->CurrentVolOff(1,0,copy);
4536 vol[1]=copy;
4537 id=pMC->CurrentVolOff(2,0,copy);
4538 vol[2]=copy;
4539 } else if(id==fIdSens3) {
4540 vol[0]=3;
4541 vol[1]=copy;
4542 id=pMC->CurrentVolOff(1,0,copy);
4543 vol[2]=copy;
4544 } else if(id==fIdSens4) {
4545 vol[0]=4;
4546 vol[1]=copy;
4547 id=pMC->CurrentVolOff(1,0,copy);
4548 vol[2]=copy;
4549 } else if(id==fIdSens5) {
4550 vol[0]=5;
4551 vol[1]=copy;
4552 id=pMC->CurrentVolOff(1,0,copy);
4553 vol[2]=copy;
4554 } else if(id==fIdSens6) {
4555 vol[0]=6;
4556 vol[1]=copy;
4557 id=pMC->CurrentVolOff(1,0,copy);
4558 vol[2]=copy;
4559 } else return;
4560 pMC->TrackPosition(position);
4561 pMC->TrackMomentum(momentum);
4562 hits[0]=position[0];
4563 hits[1]=position[1];
4564 hits[2]=position[2];
4565 hits[3]=momentum[0]*momentum[3];
4566 hits[4]=momentum[1]*momentum[3];
4567 hits[5]=momentum[2]*momentum[3];
4568 hits[6]=pMC->Edep();
4569 new(lhits[fNhits++]) AliITShit(fIshunt,gAlice->CurrentTrack(),vol,hits);
4570 }
4571}