Corrected setting of magnetic field
[u/mrichter/AliRoot.git] / ITS / AliITSvPPRsymm.cxx
CommitLineData
5d5bdf03 1 /**************************************************************************
db915fea 2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
88cb7938 16/* $Id$ */
db915fea 17
bae7e562 18//////////////////////////////////////////////////////////////////////////////
19// //
20// Inner Traking System version PPR symmetric //
21// This class contains the base procedures for the Inner Tracking System //
22// //
23// Authors: R. Barbera //
24// version 6. //
25// Created 2000. //
26// //
27// NOTE: THIS IS THE SYMMETRIC PPR geometry of the ITS. //
28// THIS WILL NOT WORK //
29// with the geometry or module classes or any analysis classes. You are //
30// strongly encouraged to uses AliITSv5. //
31// //
32//////////////////////////////////////////////////////////////////////////////
9d2da7b3 33// See AliITSvPPRsymm::StepManager().
4ae5bbc4 34#include <Riostream.h>
9d2da7b3 35#include <stdio.h>
36#include <stdlib.h>
88cb7938 37
38#include <TBRIK.h>
39#include <TCanvas.h>
40#include <TClonesArray.h>
41#include <TFile.h> // only required for Tracking function?
db915fea 42#include <TGeometry.h>
88cb7938 43#include <TLorentzVector.h>
44#include <TMath.h>
db915fea 45#include <TNode.h>
db915fea 46#include <TObjArray.h>
9d2da7b3 47#include <TObjString.h>
88cb7938 48#include <TPCON.h>
9d2da7b3 49#include <TSystem.h>
88cb7938 50#include <TTUBE.h>
51#include <TTUBS.h>
52#include <TVirtualMC.h>
db915fea 53
9d2da7b3 54#include "AliRun.h"
db915fea 55#include "AliMagF.h"
56#include "AliConst.h"
bae7e562 57#include "AliITSGeant3Geometry.h"
bae7e562 58#include "AliITSgeom.h"
bae7e562 59#include "AliITSgeomSDD.h"
88cb7938 60#include "AliITSgeomSPD.h"
bae7e562 61#include "AliITSgeomSSD.h"
88cb7938 62#include "AliITShit.h"
2c198b76 63#include "AliITSresponseSDD.h"
88cb7938 64#include "AliITSresponseSPD.h"
2c198b76 65#include "AliITSresponseSSD.h"
2c198b76 66#include "AliITSsegmentationSDD.h"
88cb7938 67#include "AliITSsegmentationSPD.h"
2c198b76 68#include "AliITSsegmentationSSD.h"
82873bc1 69#include "AliITSsimulationSDD.h"
88cb7938 70#include "AliITSsimulationSPD.h"
82873bc1 71#include "AliITSsimulationSSD.h"
88cb7938 72#include "AliITSvPPRsymm.h"
73#include "AliMagF.h"
74#include "AliRun.h"
db915fea 75
76ClassImp(AliITSvPPRsymm)
77
78//_____________________________________________________________________________
6fd3393a 79AliITSvPPRsymm::AliITSvPPRsymm() : AliITSvPPRasymm() {
db915fea 80////////////////////////////////////////////////////////////////////////
81// Standard default constructor for the ITS version 9.
82////////////////////////////////////////////////////////////////////////
db915fea 83}
84//_____________________________________________________________________________
6fd3393a 85AliITSvPPRsymm::AliITSvPPRsymm(const char *name, const char *title) : AliITSvPPRasymm(name, title){
db915fea 86////////////////////////////////////////////////////////////////////////
87// Standard constructor for the ITS version 9.
88////////////////////////////////////////////////////////////////////////
db915fea 89}
90//_____________________________________________________________________________
91AliITSvPPRsymm::~AliITSvPPRsymm() {
92////////////////////////////////////////////////////////////////////////
9d2da7b3 93// Standard destructor for the ITS version 9.
db915fea 94////////////////////////////////////////////////////////////////////////
db915fea 95}
96//_____________________________________________________________________________
97void AliITSvPPRsymm::CreateGeometry(){
98////////////////////////////////////////////////////////////////////////
99// This routine defines and Creates the geometry for version 9 of the ITS.
100////////////////////////////////////////////////////////////////////////
101
102 //INNER RADII OF THE SILICON LAYERS
e4f08ac3 103 // Float_t rl[6] = { 3.8095,7.,15.,24.,38.1,43.5765 };
db915fea 104 //THICKNESSES OF LAYERS (in % radiation length)
105 Float_t drl[6] = { 1.03,1.03,0.94,0.95,0.91,0.87 };
106 //HALF LENGTHS OF LAYERS
e4f08ac3 107 // Float_t dzl[6] = { 14.35,14.35,25.1,32.1,49.405,55.27 };
db915fea 108 //LENGTHS OF END-LADDER BOXES (ALL INCLUDED)
e4f08ac3 109 // Float_t dzb[6] = { 12.4,12.4,13.5,15.,7.5,7.5 };
db915fea 110 //THICKNESSES OF END-LADDER BOXES (ALL INCLUDED)
e4f08ac3 111 // Float_t drb[6] = { rl[1]-rl[0],0.2,5.,5.,4.,4. };
9d2da7b3 112
db915fea 113
e4f08ac3 114 Float_t dits[100], rlim, zmax;
9d2da7b3 115 // Float_t zpos;
116 // Float_t pcits[50]
e4f08ac3 117 Float_t ztpc;
9d2da7b3 118 Int_t idrotm[1999], i;
e4f08ac3 119 Float_t dgh[100];
9d2da7b3 120
121
122 // Define some variables for SPD
123
124 Float_t dits1[3], di101[3], di107[3], di10b[3], di106[3]; // for layer 1
125 Float_t di103[3], di10a[3], di102[3]; // for layer 1
126 Float_t dits2[3], di1d1[3], di1d7[3], di20b[3], di1d6[3]; // for layer 2
127 Float_t di1d3[3], di20a[3], di1d2[3]; // for layer 2
128 Float_t di108[3], di104[3]; // for both layers
129
a8c83dcb 130 Float_t ddet1=200.; // total detector thickness on layer 1 (micron)
131 Float_t dchip1=200.; // total chip thickness on layer 1 (micron)
9d2da7b3 132
a8c83dcb 133 Float_t ddet2=200.; // total detector thickness on layer 2 (micron)
134 Float_t dchip2=200.; // total chip thickness on layer 2 (micron)
9d2da7b3 135
214030ed 136 Float_t dbus=300.; // total bus thickness on both layers (micron)
9d2da7b3 137
138 ddet1 = GetThicknessDet1();
139 ddet2 = GetThicknessDet2();
140 dchip1 = GetThicknessChip1();
141 dchip2 = GetThicknessChip2();
142
9d2da7b3 143 if(ddet1 < 100. || ddet1 > 300.) {
1e9111db 144 cout << "ITS - WARNING: the detector thickness for layer 1 is outside the range of [100,300] microns."
a8c83dcb 145 " The default value of 200 microns will be used." << endl;
146 ddet1=200.;
9d2da7b3 147 }
148
149 if(ddet2 < 100. || ddet2 > 300.) {
1e9111db 150 cout << "ITS - WARNING: the detector thickness for layer 2 is outside the range of [100,300] microns."
a8c83dcb 151 " The default value of 200 microns will be used." << endl;
152 ddet2=200.;
9d2da7b3 153 }
154
214030ed 155 if(dchip1 < 100. || dchip1 > 300.) {
156 cout << "ITS - WARNING: the chip thickness for layer 1 is outside the range of [100,300] microns."
a8c83dcb 157 " The default value of 200 microns will be used." << endl;
158 dchip1=200.;
9d2da7b3 159 }
160
214030ed 161 if(dchip2 < 100. || dchip2 > 300.) {
162 cout << "ITS - WARNING: the chip thickness for layer 2 is outside the range of [100,300] microns."
a8c83dcb 163 " The default value of 200 microns will be used." << endl;
164 dchip2=200.;
9d2da7b3 165 }
e4f08ac3 166
167 Int_t rails = 1; // flag for rails (1 --> rails in; 0 --> rails out)
9d2da7b3 168
e4f08ac3 169 Int_t fluid = 1; // flag for the cooling fluid (1 --> water; 0 --> freon)
170
171 rails = GetRails();
172
173 fluid = GetCoolingFluid();
174
175 if(rails != 0 && rails != 1) {
1e9111db 176 cout << "ITS - WARNING: the switch for rails is not set neither to 0 (rails out) nor to 1 (rails in)."
e4f08ac3 177 " The default value of 1 (rails in) will be used." << endl;
178
179 }
180
181 if(fluid != 0 && fluid != 1) {
1e9111db 182 cout << "ITS - WARNING: the switch for cooling fluid is not set neither to 0 (freon) nor to 1 (water)."
e4f08ac3 183 " The default value of 1 (water) will be used." << endl;
184 }
185
1e9111db 186 cout << "ITS: Detector thickness on layer 1 is set to " << ddet1 << " microns." << endl;
187 cout << "ITS: Chip thickness on layer 1 is set to " << dchip1 << " microns." << endl;
188 cout << "ITS: Detector thickness on layer 2 is set to " << ddet2 << " microns." << endl;
189 cout << "ITS: Chip thickness on layer 2 is set to " << dchip2 << " microns." << endl;
e4f08ac3 190 if(rails == 0 ) {
1e9111db 191 cout << "ITS: Rails are out." << endl;
e4f08ac3 192 } else {
1e9111db 193 cout << "ITS: Rails are in." << endl;
e4f08ac3 194 }
195 if(fluid == 0 ) {
1e9111db 196 cout << "ITS: The cooling fluid is freon." << endl;
e4f08ac3 197 } else {
1e9111db 198 cout << "ITS: The cooling fluid is water." << endl;
e4f08ac3 199 }
200
9d2da7b3 201 ddet1 = ddet1*0.0001/2.; // conversion from tot length in um to half in cm
202 ddet2 = ddet2*0.0001/2.; // conversion from tot length in um to half in cm
203 dchip1 = dchip1*0.0001/2.;// conversion from tot length in um to half in cm
204 dchip2 = dchip2*0.0001/2.;// conversion from tot length in um to half in cm
e4f08ac3 205 dbus = dbus*0.0001/2.; // conversion from tot length in um to half in cm
9d2da7b3 206
207 Float_t deltax, deltay;
208
6fd3393a 209 Int_t thickness = GetMinorVersion()/10;
210 Int_t option = GetMinorVersion() - 10*thickness;
9d2da7b3 211
212
213 // Define some variables for SDD
214
215
216 Float_t sin30, cos30;
217
218 // SDD electronics+services main volumes
219 Float_t I018dits[3], I024dits[3], I047dits[3], I048dits[3];
220
221 // SDD detector ladder
222
223 Float_t I302dits[3], I402dits[3], I004dits[3], I005dits[3];
224 Float_t Y_SDD_sep = 0.20;
225 Float_t ySDD;
226 Int_t iSDD;
227 Float_t Z_SDD_lay3[6] = {18.55, 10.95, 3.70, -3.70, -11.20, -18.35};
228 Float_t Z_SDD_lay4[8] = {25.75, 18.60, 11.00, 3.70, -3.70, -11.20, -18.45, -26.05};
229
230 // ladder foot and end-ladder (frame and cooling)
231 Float_t I028dits[3], I420dits[3], I421dits[3], I422dits[6], I423dits[3];
232 Float_t I424dits[3], xI424, yI424;
233 Float_t I425dits[3];
234 Int_t indI425;
235 Float_t I029dits[4], I030dits[4], I031dits[3], I032dits[3];
236
237 // SDD ladder frame and cooling
238 Float_t SDD_CoolPipe[3] = {1.7000, -0.5500, 0.0000};
239 Float_t I035dits[3], I037dits[3], I038dits[3];
240 Float_t I039dits[3], xI039, yI039;
241 Float_t I041dits[5];
242
243 // SDD hybrid, chips and capacitors
244 Float_t I050dits[3], xI050, yI050;
245 Float_t I052dits[3], xI052, yI052;
246 Float_t I042dits[3], xI042, yI042;
247 Float_t xI042space = 0.17;
248 Float_t I043dits[3], xI043, yI043;
249 Float_t xI043space = 0.17;
250 Float_t zchip, zChipSpace;
251 Float_t I051dits[3], xI051, yI051, zI051, yI051space, xcap;
252 Int_t ichip, icap;
253
254 // SDD microcables
255 Float_t I044dits[4], xI044, yI044, volI044;
256 Float_t xHV, yHV, zHV, xLV, yLV, zLV;
257 Char_t HVname[5], LVname[5];
258
259
260 // Define media off-set
db915fea 261
262 Int_t *idtmed = fIdtmed->GetArray()-199;
9d2da7b3 263
db915fea 264
265 // Rotation matrices
266
9d2da7b3 267 // SPD - option 'a' (this is NOT the default so leave commented)
268
269
270 if (option == 1) {
271
272 AliMatrix(idrotm[201],90.0,90.0,90.0,180.0,0.0,0.0);
273 AliMatrix(idrotm[202],90.0,90.0,90.0,0.0,0.0,0.0);
274 AliMatrix(idrotm[203],90.0,350.0,90.0,260.0,0.0,0.0);
275 AliMatrix(idrotm[204],90.0,170.0,90.0,80.0,0.0,0.0);
276 AliMatrix(idrotm[205],90.0,10.0,90.0,100.0,0.0,0.0);
277 AliMatrix(idrotm[206],90.0,190.0,90.0,280.0,0.0,0.0);
278 AliMatrix(idrotm[207],90.0,342.0,90.0,72.0,0.0,0.0);
279 AliMatrix(idrotm[208],90.0,156.999893,90.0,246.999893,0.0,0.0);
280 AliMatrix(idrotm[209],90.0,147.999802,90.0,237.999893,0.0,0.0);
281 AliMatrix(idrotm[210],90.0,138.999802,90.0,228.999802,0.0,0.0);
282 AliMatrix(idrotm[211],90.0,129.999802,90.0,219.999802,0.0,0.0);
283 AliMatrix(idrotm[212],90.0,36.7896,90.0,126.789597,0.0,0.0);
284 AliMatrix(idrotm[213],90.0,343.579712,90.0,73.579697,0.0,0.0);
285 AliMatrix(idrotm[214],90.0,95.413696,90.0,185.413696,0.0,0.0);
286 AliMatrix(idrotm[215],90.0,5.4141,90.0,95.414101,0.0,0.0);
287 AliMatrix(idrotm[216],90.0,318.296906,90.0,48.296902,0.0,0.0);
288 AliMatrix(idrotm[217],90.0,67.000099,90.0,157.000107,0.0,0.0);
289 AliMatrix(idrotm[218],90.0,337.003998,90.0,67.003998,0.0,0.0);
290 AliMatrix(idrotm[219],90.0,247.000305,90.0,337.000305,0.0,0.0);
291 AliMatrix(idrotm[220],90.0,305.633514,90.0,35.633499,0.0,0.0);
292 AliMatrix(idrotm[221],90.0,58.000198,90.0,148.000198,0.0,0.0);
293 AliMatrix(idrotm[222],90.0,327.997101,90.0,57.997101,0.0,0.0 );
294 AliMatrix(idrotm[223],90.0,237.994202,90.0,327.994202,0.0,0.0);
295 AliMatrix(idrotm[224],90.0,296.627502,90.0,26.627399,0.0,0.0);
296 AliMatrix(idrotm[225],90.0,48.994099,90.0,138.994095,0.0,0.0);
297 AliMatrix(idrotm[226],90.0,318.990997,90.0,48.991001,0.0,0.0);
298 AliMatrix(idrotm[227],90.0,228.988205,90.0,318.98819,0.0,0.0);
299 AliMatrix(idrotm[228],90.0,287.621399,90.0,17.621401,0.0,0.0);
300 AliMatrix(idrotm[229],90.0,39.988098,90.0,129.988098,0.0,0.0);
301 AliMatrix(idrotm[230],90.0,309.984985,90.0,39.985001,0.0,0.0);
302 AliMatrix(idrotm[231],90.0,327.2612,90.0,57.2612,0.0,0.0);
303 AliMatrix(idrotm[232],90.0,237.261398,90.0,327.261414,0.0,0.0);
304 AliMatrix(idrotm[233],90.0,252.000504,90.0,342.000488,0.0,0.0 );
305 AliMatrix(idrotm[234],90.0,71.9991,90.0,161.9991,0.0,0.0);
306 AliMatrix(idrotm[235],90.0,270.0,90.0,0.0,0.0,0.0);
307 AliMatrix(idrotm[236],90.0,180.013702,90.0,270.013702,0.0,0.0);
9bf5e8ff 308 AliMatrix(idrotm[237],90.0,180.0,90.0,90.0,0.0,0.0);
9d2da7b3 309 AliMatrix(idrotm[238],90.0,144.0,90.0,234.0,0.0,0.0);
310 AliMatrix(idrotm[239],90.0,216.0,90.0,306.0,0.0,0.0);
311 AliMatrix(idrotm[240],90.0,288.0,90.0,18.0,0.0,0.0);
312 AliMatrix(idrotm[241],90.0,324.0,90.0,54.0,0.0,0.0);
313 AliMatrix(idrotm[242],90.0,36.0,90.0,126.0,0.0,0.0);
314 AliMatrix(idrotm[243],90.0,108.0,90.0,198.0,0.0,0.0);
9bf5e8ff 315 AliMatrix(idrotm[244],90.0,180.0,90.0,270.0,0.0,0.0);
316 AliMatrix(idrotm[245],90.0,162.0,90.0,252.0,0.0,0.0);
317 AliMatrix(idrotm[246],90.0,310.0,90.0,40.0,0.0,0.0);
318 AliMatrix(idrotm[247],90.0,319.0,90.0,49.0,0.0,0.0);
319 AliMatrix(idrotm[248],90.0,328.0,90.0,58.0,0.0,0.0);
320 AliMatrix(idrotm[249],90.0,337.0,90.0,67.0,0.0,0.0);
9d2da7b3 321 AliMatrix(idrotm[1003],90.0,73.5,90.0,163.5,0.0,0.0);
322 AliMatrix(idrotm[1011],90.0,342.0,90.0,72.0,0.0,0.0);
323 AliMatrix(idrotm[1039],90.0,72.0,90.0,162.0,0.0,0.0);
324 AliMatrix(idrotm[1043],90.0,66.91,90.0,156.91,0.0,0.0);
325 AliMatrix(idrotm[1065],90.0,144.0,90.0,234.0,0.0,0.0);
326 AliMatrix(idrotm[1078],90.0,180.0,90.0,270.0,0.0,0.0);
327 AliMatrix(idrotm[1088],90.0,57.41,90.0,147.41,0.0,0.0);
328 AliMatrix(idrotm[1089],90.0,333.0,90.0,63.0,0.0,0.0);
329 AliMatrix(idrotm[1090],90.0,351.0,90.0,81.0,0.0,0.0);
330 AliMatrix(idrotm[1091],90.0,216.0,90.0,306.0,0.0,0.0);
331 AliMatrix(idrotm[1092],90.0,27.0,90.0,117.0,0.0,0.0);
332 AliMatrix(idrotm[1093],90.0,18.0,90.0,108.0,0.0,0.0);
333 AliMatrix(idrotm[1094],90.0,9.0,90.0,99.0,0.0,0.0);
334 AliMatrix(idrotm[1104],90.0,252.0,90.0,342.0,0.0,0.0);
335 AliMatrix(idrotm[1106],90.0,36.0,90.0,126.0,0.0,0.0);
336 AliMatrix(idrotm[1107],90.0,108.0,90.0,198.0,0.0,0.0);
337 AliMatrix(idrotm[1108],90.0,324.0,90.0,54.0,180.0,0.0);
338 AliMatrix(idrotm[1109],90.0,0.0,90.0,90.0,180.0,0.0);
339 AliMatrix(idrotm[1110],90.0,36.0,90.0,126.0,180.0,0.0);
340 AliMatrix(idrotm[1111],90.0,72.0,90.0,162.0,180.0,0.0);
341 AliMatrix(idrotm[1112],90.0,108.0,90.0,198.0,180.0,0.0);
342 AliMatrix(idrotm[1113],90.0,144.0,90.0,234.0,180.0,0.0);
343 AliMatrix(idrotm[1114],90.0,180.0,90.0,270.0,180.0,0.0);
344 AliMatrix(idrotm[1115],90.0,216.0,90.0,306.0,180.0,0.0);
345 AliMatrix(idrotm[1116],90.0,252.0,90.0,342.0,180.0,0.0);
346 AliMatrix(idrotm[1117],90.0,288.0,90.0,18.0,0.0,0.0);
347 AliMatrix(idrotm[1118],90.0,288.0,90.0,18.0,180.0,0.0);
348 AliMatrix(idrotm[1130],90.0,324.0,90.0,54.0,0.0,0.0);
349
350 }
db915fea 351
9d2da7b3 352 // SPD - option 'b' (this is the default)
353
354 if (option == 2) {
355
9bf5e8ff 356 AliMatrix(idrotm[201],90.0,0.0,90.0,90.0,0.0,0.0);
9d2da7b3 357 AliMatrix(idrotm[202],90.0,90.0,90.0,0.0,0.0,0.0);
358 AliMatrix(idrotm[203],90.0,350.0,90.0,260.0,0.0,0.0);
359 AliMatrix(idrotm[204],90.0,170.0,90.0,80.0,0.0,0.0);
360 AliMatrix(idrotm[205],90.0,10.0,90.0,100.0,0.0,0.0);
361 AliMatrix(idrotm[206],90.0,190.0,90.0,280.0,0.0,0.0);
362 AliMatrix(idrotm[207],90.0,342.0,90.0,72.0,0.0,0.0);
363 AliMatrix(idrotm[208],90.0,156.999893,90.0,246.999893,0.0,0.0);
364 AliMatrix(idrotm[209],90.0,147.999802,90.0,237.999893,0.0,0.0);
365 AliMatrix(idrotm[210],90.0,138.999802,90.0,228.999802,0.0,0.0);
366 AliMatrix(idrotm[211],90.0,129.999802,90.0,219.999802,0.0,0.0);
367 AliMatrix(idrotm[212],90.0,36.7896,90.0,126.789597,0.0,0.0);
368 AliMatrix(idrotm[213],90.0,343.579712,90.0,73.579697,0.0,0.0);
369 AliMatrix(idrotm[214],90.0,95.413696,90.0,185.413696,0.0,0.0);
370 AliMatrix(idrotm[215],90.0,5.4141,90.0,95.414101,0.0,0.0);
371 AliMatrix(idrotm[216],90.0,318.296906,90.0,48.296902,0.0,0.0);
372 AliMatrix(idrotm[217],90.0,67.000099,90.0,157.000107,0.0,0.0);
373 AliMatrix(idrotm[218],90.0,337.003998,90.0,67.003998,0.0,0.0);
374 AliMatrix(idrotm[219],90.0,247.000305,90.0,337.000305,0.0,0.0);
375 AliMatrix(idrotm[220],90.0,305.633514,90.0,35.633499,0.0,0.0);
376 AliMatrix(idrotm[221],90.0,58.000198,90.0,148.000198,0.0,0.0);
377 AliMatrix(idrotm[222],90.0,327.997101,90.0,57.997101,0.0,0.0);
378 AliMatrix(idrotm[223],90.0,237.994202,90.0,327.994202,0.0,0.0);
379 AliMatrix(idrotm[224],90.0,296.627502,90.0,26.627399,0.0,0.0);
380 AliMatrix(idrotm[225],90.0,48.994099,90.0,138.994095,0.0,0.0);
381 AliMatrix(idrotm[226],90.0,318.990997,90.0,48.991001,0.0,0.0);
382 AliMatrix(idrotm[227],90.0,228.988205,90.0,318.98819,0.0,0.0);
383 AliMatrix(idrotm[228],90.0,287.621399,90.0,17.621401,0.0,0.0);
384 AliMatrix(idrotm[229],90.0,39.988098,90.0,129.988098,0.0,0.0);
385 AliMatrix(idrotm[230],90.0,309.984985,90.0,39.985001,0.0,0.0);
386 AliMatrix(idrotm[231],90.0,327.2612,90.0,57.2612,0.0,0.0);
387 AliMatrix(idrotm[232],90.0,237.261398,90.0,327.261414,0.0,0.0);
388 AliMatrix(idrotm[233],90.0,252.000504,90.0,342.000488,0.0,0.0);
389 AliMatrix(idrotm[234],90.0,71.9991,90.0,161.9991,0.0,0.0);
390 AliMatrix(idrotm[235],90.0,270.0,90.0,0.0,0.0,0.0);
391 AliMatrix(idrotm[236],90.0,180.013702,90.0,270.013702,0.0,0.0);
392 AliMatrix(idrotm[237],90.0,90.0,90.0,180.0,0.0,0.0);
9bf5e8ff 393 AliMatrix(idrotm[238],90.0,180.0,90.0,270.0,0.0,0.0);
394 AliMatrix(idrotm[239],90.0,162.0,90.0,252.0,0.0,0.0);
395 AliMatrix(idrotm[240],90.0,310.0,90.0,40.0,0.0,0.0);
396 AliMatrix(idrotm[241],90.0,319.0,90.0,49.0,0.0,0.0);
397 AliMatrix(idrotm[242],90.0,328.0,90.0,58.0,0.0,0.0);
398 AliMatrix(idrotm[243],90.0,337.0,90.0,67.0,0.0,0.0);
9d2da7b3 399 AliMatrix(idrotm[244],90.0,216.0,90.0,306.0,0.0,0.0);
400 AliMatrix(idrotm[245],90.0,36.0,90.0,126.0,0.0,0.0);
401 AliMatrix(idrotm[246],90.0,108.0,90.0,198.0,0.0,0.0);
402 AliMatrix(idrotm[247],90.0,144.0,90.0,234.0,0.0,0.0);
403 AliMatrix(idrotm[248],90.0,288.0,90.0,18.0,0.0,0.0);
404 AliMatrix(idrotm[249],90.0,324.0,90.0,54.0,0.0,0.0);
405 AliMatrix(idrotm[1003],90.0,73.5,90.0,163.5,0.0,0.0);
406 AliMatrix(idrotm[1011],90.0,342.0,90.0,72.0,0.0,0.0);
407 AliMatrix(idrotm[1039],90.0,72.0,90.0,162.0,0.0,0.0);
408 AliMatrix(idrotm[1043],90.0,66.91,90.0,156.91,0.0,0.0);
409 AliMatrix(idrotm[1065],90.0,144.0,90.0,234.0,0.0,0.0);
410 AliMatrix(idrotm[1078],90.0,180.0,90.0,270.0,0.0,0.0);
411 AliMatrix(idrotm[1088],90.0,57.41,90.0,147.41,0.0,0.0);
412 AliMatrix(idrotm[1089],90.0,333.0,90.0,63.0,0.0,0.0);
413 AliMatrix(idrotm[1090],90.0,351.0,90.0,81.0,0.0,0.0);
414 AliMatrix(idrotm[1091],90.0,216.0,90.0,306.0,0.0,0.0);
415 AliMatrix(idrotm[1092],90.0,27.0,90.0,117.0,0.0,0.0);
416 AliMatrix(idrotm[1093],90.0,18.0,90.0,108.0,0.0,0.0);
417 AliMatrix(idrotm[1094],90.0,9.0,90.0,99.0,0.0,0.0);
418 AliMatrix(idrotm[1104],90.0,252.0,90.0,342.0,0.0,0.0);
419 AliMatrix(idrotm[1106],90.0,36.0,90.0,126.0,0.0,0.0);
420 AliMatrix(idrotm[1107],90.0,108.0,90.0,198.0,0.0,0.0);
421 AliMatrix(idrotm[1108],90.0,324.0,90.0,54.0,180.0,0.0);
422 AliMatrix(idrotm[1109],90.0,0.0,90.0,90.0,180.0,0.0);
423 AliMatrix(idrotm[1110],90.0,36.0,90.0,126.0,180.0,0.0);
424 AliMatrix(idrotm[1111],90.0,72.0,90.0,162.0,180.0,0.0);
425 AliMatrix(idrotm[1112],90.0,108.0,90.0,198.0,180.0,0.0);
426 AliMatrix(idrotm[1113],90.0,144.0,90.0,234.0,180.0,0.0);
427 AliMatrix(idrotm[1114],90.0,180.0,90.0,270.0,180.0,0.0);
428 AliMatrix(idrotm[1115],90.0,216.0,90.0,306.0,180.0,0.0);
429 AliMatrix(idrotm[1116],90.0,252.0,90.0,342.0,180.0,0.0);
430 AliMatrix(idrotm[1117],90.0,288.0,90.0,18.0,0.0,0.0);
431 AliMatrix(idrotm[1118],90.0,288.0,90.0,18.0,180.0,0.0);
432 AliMatrix(idrotm[1130],90.0,324.0,90.0,54.0,0.0,0.0);
433
434 }
435
db915fea 436 // SDD
437
438 AliMatrix(idrotm[301],0.0,0.0,90.0,90.0,90.0,180.0);
439 AliMatrix(idrotm[302],0.0,0.0,90.0,90.0,90.0,0.0);
440 AliMatrix(idrotm[303],180.0,0.0,90.0,90.0,90.0,0.0);
441 AliMatrix(idrotm[304],180.0,0.0,90.0,90.0,90.0,180.0);
442 AliMatrix(idrotm[305],90.0,347.14,90.0,77.14,0.0,0.0);
443 AliMatrix(idrotm[306],90.0,321.43,90.0,51.43,0.0,0.0);
444 AliMatrix(idrotm[307],90.0,295.71,90.0,25.71,0.0,0.0);
445 AliMatrix(idrotm[308],90.0,244.29,90.0,334.29,0.0,0.0);
446 AliMatrix(idrotm[309],90.0,218.57,90.0,308.57,0.0,0.0);
447 AliMatrix(idrotm[310],90.0,167.14,90.0,257.14,0.0,0.0);
448 AliMatrix(idrotm[311],90.0,141.43,90.0,231.43,0.0,0.0);
449 AliMatrix(idrotm[312],90.0,0.0,0.0,0.0,90.0,270.0);
450 AliMatrix(idrotm[313],90.0,115.71,90.0,205.71,0.0,0.0);
451 AliMatrix(idrotm[314],90.0,335.45,90.0,65.45,0.0,0.0);
452 AliMatrix(idrotm[315],90.0,319.09,90.0,49.09,0.0,0.0);
453 AliMatrix(idrotm[316],90.0,302.73,90.0,32.73,0.0,0.0);
454 AliMatrix(idrotm[317],90.0,286.36,90.0,16.36,0.0,0.0);
455 AliMatrix(idrotm[318],90.0,270.0,90.0,360.0,0.0,0.0);
456 AliMatrix(idrotm[319],90.0,253.64,90.0,343.64,0.0,0.0);
457 AliMatrix(idrotm[320],90.0,237.27,90.0,327.27,0.0,0.0);
458 AliMatrix(idrotm[321],90.0,12.86,90.0,102.86,0.0,0.0);
459 AliMatrix(idrotm[322],90.0,220.91,90.0,310.91,0.0,0.0);
460 AliMatrix(idrotm[323],90.0,204.55,90.0,294.55,0.0,0.0);
461 AliMatrix(idrotm[324],90.0,188.18,90.0,278.18,0.0,0.0);
462 AliMatrix(idrotm[325],90.0,171.82,90.0,261.82,0.0,0.0);
463 AliMatrix(idrotm[326],90.0,155.45,90.0,245.45,0.0,0.0);
464 AliMatrix(idrotm[327],90.0,139.09,90.0,229.09,0.0,0.0);
465 AliMatrix(idrotm[328],90.0,122.73,90.0,212.73,0.0,0.0);
466 AliMatrix(idrotm[329],90.0,106.36,90.0,196.36,0.0,0.0);
467 AliMatrix(idrotm[330],90.0,73.64,90.0,163.64,0.0,0.0);
468 AliMatrix(idrotm[331],90.0,40.91,90.0,130.91,0.0,0.0);
469 AliMatrix(idrotm[332],90.0,24.55,90.0,114.55,0.0,0.0);
470 AliMatrix(idrotm[333],90.0,38.57,90.0,128.57,0.0,0.0);
471 AliMatrix(idrotm[334],90.0,351.82,90.0,81.82,0.0,0.0);
472 AliMatrix(idrotm[335],90.0,8.18,90.0,98.18,0.0,0.0);
473 AliMatrix(idrotm[336],90.0,64.29,90.0,154.29,0.0,0.0);
474 AliMatrix(idrotm[337],111.0,300.0,21.0,300.0,90.0,30.0);
475 AliMatrix(idrotm[338],69.0,240.0,159.0,240.0,90.0,150.0);
476 AliMatrix(idrotm[339],111.0,240.0,21.0,240.0,90.0,150.0);
477 AliMatrix(idrotm[340],69.0,300.0,159.0,300.0,90.0,30.0);
478 AliMatrix(idrotm[341],128.0,0.0,38.0,0.0,90.0,270.0);
479 AliMatrix(idrotm[342],90.0,240.0,180.0,0.0,90.0,330.);
480 AliMatrix(idrotm[343],90.0,120.0,180.0,0.0,90.0,210.0);
481 AliMatrix(idrotm[344],90.0,0.0,180.0,0.0,90.0,90.0);
482 AliMatrix(idrotm[345],90.0,180.0,90.0,90.0,0.0,0.0);
483 AliMatrix(idrotm[346],90.0,300.0,90.0,30.0,0.0,0.0);
484 AliMatrix(idrotm[347],90.0,240.0,90.0,150.0,0.0,0.0);
485 AliMatrix(idrotm[348],90.0,180.0,0.0,0.0,90.0,270.0);
486 AliMatrix(idrotm[349],90.0,235.0,90.0,145.0,0.0,0.0);
487 AliMatrix(idrotm[350],90.0,90.0,90.0,180.0,0.0,0.0);
488 AliMatrix(idrotm[351],90.0,305.0,90.0,35.0,0.0,0.0);
489 AliMatrix(idrotm[352],0.0,0.0,90.0,0.0,90.0,90.0);
490 AliMatrix(idrotm[353],90.0,60.0,90.0,150.0,0.0,0.0);
491 AliMatrix(idrotm[354],90.0,120.0,90.0,30.0,0.0,0.0);
492 AliMatrix(idrotm[355],90.0,180.0,90.0,90.0,180.0,0.0);
493 AliMatrix(idrotm[356],90.0,270.0,90.0,0.0,0.0,0.0);
494 AliMatrix(idrotm[366],90.0,57.27,90.0,147.27,0.0,0.0);
495 AliMatrix(idrotm[386],90.0,192.86,90.0,282.86,0.0,0.0);
496
497 // SSD
498
499 AliMatrix(idrotm[501],90.0,148.24,90.0,238.24,0.0,0.0);
500 AliMatrix(idrotm[503],90.0,137.65,90.0,227.65,0.0,0.0);
501 AliMatrix(idrotm[504],90.0,127.06,90.0,217.06,0.0,0.0);
502 AliMatrix(idrotm[505],90.0,116.47,90.0,206.47,0.0,0.0);
503 AliMatrix(idrotm[506],90.0,105.88,90.0,195.88,0.0,0.0);
504 AliMatrix(idrotm[507],90.0,95.29,90.0,185.29,0.0,0.0);
505 AliMatrix(idrotm[508],90.0,84.71,90.0,174.71,0.0,0.0);
506 AliMatrix(idrotm[509],90.0,74.12,90.0,164.12,0.0,0.0);
507 AliMatrix(idrotm[510],90.0,63.53,90.0,153.53,0.0,0.0);
508 AliMatrix(idrotm[511],90.0,52.94,90.0,142.94,0.0,0.0);
509 AliMatrix(idrotm[512],90.0,42.35,90.0,132.35,0.0,0.0);
510 AliMatrix(idrotm[513],90.0,31.76,90.0,121.76,0.0,0.0);
511 AliMatrix(idrotm[514],90.0,10.59,90.0,100.59,0.0,0.0);
512 AliMatrix(idrotm[515],90.0,349.41,90.0,79.41,0.0,0.0);
513 AliMatrix(idrotm[516],90.0,338.82,90.0,68.82,0.0,0.0);
514 AliMatrix(idrotm[517],90.0,328.24,90.0,58.24,0.0,0.0);
515 AliMatrix(idrotm[518],90.0,317.65,90.0,47.65,0.0,0.0);
516 AliMatrix(idrotm[519],90.0,307.06,90.0,37.06,0.0,0.0);
517 AliMatrix(idrotm[520],90.0,296.47,90.0,26.47,0.0,0.0);
518 AliMatrix(idrotm[521],90.0,285.88,90.0,15.88,0.0,0.0);
519 AliMatrix(idrotm[522],90.0,275.29,90.0,5.29,0.0,0.0);
520 AliMatrix(idrotm[523],90.0,264.71,90.0,354.71,0.0,0.0);
521 AliMatrix(idrotm[524],90.0,254.12,90.0,344.12,0.0,0.0);
522 AliMatrix(idrotm[525],90.0,243.53,90.0,333.53,0.0,0.0);
523 AliMatrix(idrotm[526],90.0,232.94,90.0,322.94,0.0,0.0);
524 AliMatrix(idrotm[527],90.0,222.35,90.0,312.35,0.0,0.0);
525 AliMatrix(idrotm[528],90.0,211.76,90.0,301.76,0.0,0.0);
526 AliMatrix(idrotm[529],90.0,190.59,90.0,280.59,0.0,0.0);
527 AliMatrix(idrotm[530],90.0,169.41,90.0,259.41,0.0,0.0);
528 AliMatrix(idrotm[531],90.0,158.82,90.0,248.82,0.0,0.0);
529 AliMatrix(idrotm[532],90.0,360.0,90.0,90.0,0.0,0.0);
530 AliMatrix(idrotm[533],90.0,180.0,90.0,270.0,0.0,0.0);
531 AliMatrix(idrotm[534],90.0,189.47,90.0,279.47,0.0,0.0);
9d2da7b3 532 AliMatrix(idrotm[535],90.0,198.95,90.0,288.95,0.0,0.0);
db915fea 533 AliMatrix(idrotm[537],90.0,217.89,90.0,307.89,0.0,0.0);
534 AliMatrix(idrotm[538],90.0,227.37,90.0,317.37,0.0,0.0);
535 AliMatrix(idrotm[539],90.0,236.84,90.0,326.84,0.0,0.0);
536 AliMatrix(idrotm[540],90.0,246.32,90.0,336.32,0.0,0.0);
537 AliMatrix(idrotm[541],90.0,255.79,90.0,345.79,0.0,0.0);
538 AliMatrix(idrotm[542],90.0,265.26,90.0,355.26,0.0,0.0);
539 AliMatrix(idrotm[543],90.0,274.74,90.0,4.74,0.0,0.0);
540 AliMatrix(idrotm[544],90.0,284.21,90.0,14.21,0.0,0.0);
541 AliMatrix(idrotm[545],90.0,293.68,90.0,23.68,0.0,0.0);
542 AliMatrix(idrotm[546],90.0,303.16,90.0,33.16,0.0,0.0);
543 AliMatrix(idrotm[547],90.0,312.63,90.0,42.63,0.0,0.0);
544 AliMatrix(idrotm[548],90.0,322.11,90.0,52.11,0.0,0.0);
545 AliMatrix(idrotm[549],90.0,331.58,90.0,61.58,0.0,0.0);
546 AliMatrix(idrotm[550],90.0,341.05,90.0,71.05,0.0,0.0);
547 AliMatrix(idrotm[551],90.0,350.53,90.0,80.53,0.0,0.0);
548 AliMatrix(idrotm[552],90.0,9.47,90.0,99.47,0.0,0.0);
9d2da7b3 549 AliMatrix(idrotm[553],90.0,18.95,90.0,108.95,0.0,0.0);
db915fea 550 AliMatrix(idrotm[555],90.0,37.89,90.0,127.89,0.0,0.0);
551 AliMatrix(idrotm[556],90.0,47.37,90.0,137.37,0.0,0.0);
552 AliMatrix(idrotm[557],90.0,56.84,90.0,146.84,0.0,0.0);
553 AliMatrix(idrotm[558],90.0,66.32,90.0,156.32,0.0,0.0);
554 AliMatrix(idrotm[559],90.0,75.79,90.0,165.79,0.0,0.0);
555 AliMatrix(idrotm[560],90.0,85.26,90.0,175.26,0.0,0.0);
556 AliMatrix(idrotm[561],90.0,94.74,90.0,184.74,0.0,0.0);
557 AliMatrix(idrotm[562],90.0,104.21,90.0,194.21,0.0,0.0);
558 AliMatrix(idrotm[563],90.0,113.68,90.0,203.68,0.0,0.0);
559 AliMatrix(idrotm[564],90.0,123.16,90.0,213.16,0.0,0.0);
560 AliMatrix(idrotm[565],90.0,132.63,90.0,222.63,0.0,0.0);
561 AliMatrix(idrotm[566],90.0,142.11,90.0,232.11,0.0,0.0);
562 AliMatrix(idrotm[567],90.0,151.58,90.0,241.58,0.0,0.0);
563 AliMatrix(idrotm[568],90.0,161.05,90.0,251.05,0.0,0.0);
564 AliMatrix(idrotm[569],90.0,170.53,90.0,260.53,0.0,0.0);
565 AliMatrix(idrotm[570],90.0,180.0,90.0,90.0,180.0,0.0);
566 AliMatrix(idrotm[571],90.0,0.0,0.0,0.0,90.0,270.0);
567 AliMatrix(idrotm[572],90.0,180.0,0.0,0.0,90.0,270.0);
568 AliMatrix(idrotm[573],90.0,180.0,90.0,90.0,0.0,0.0);
569 AliMatrix(idrotm[575],90.0,120.0,180.0,0.0,90.0,210.0);
570 AliMatrix(idrotm[576],65.71,300.0,90.0,30.0,24.29,120.0);
571 AliMatrix(idrotm[577],114.29,300.0,90.0,30.0,155.71,120.0);
572 AliMatrix(idrotm[579],65.71,240.0,90.0,150.0,24.29,60.0);
573 AliMatrix(idrotm[580],114.29,240.0,90.0,150.0,155.71,60.0);
574 AliMatrix(idrotm[581],90.0,240.0,180.0,0.0,90.0,330.0);
575 AliMatrix(idrotm[583],90.0,0.0,180.0,0.0,90.0,90.0);
576 AliMatrix(idrotm[584],90.0,180.0,180.0,0.0,90.0,90.0);
577 AliMatrix(idrotm[586],180.0,0.0,90.0,90.0,90.0,0.0);
578 AliMatrix(idrotm[618],90.0,201.18,90.0,291.18,0.0,0.0);
579 AliMatrix(idrotm[620],90.0,28.42,90.0,118.42,0.0,0.0);
580 AliMatrix(idrotm[623],90.0,208.42,90.0,298.42,0.0,0.0);
581 AliMatrix(idrotm[633],132.46,0.0,90.0,90.0,42.46,360.0);
582 AliMatrix(idrotm[653],90.0,21.18,90.0,111.18,0.0,0.0);
583
9d2da7b3 584
585 // SDD cone
586
587 AliMatrix(idrotm[846],90.0,300.0,90.0,30.0,0.0,0.0);
588 AliMatrix(idrotm[851],90.0,305.0,90.0,35.0,0.0,0.0);
589 AliMatrix(idrotm[853],90.0,60.0,90.0,150.0,0.0,0.0);
590 AliMatrix(idrotm[856],90.0,0.0,90.0,90.0,180.0,0.0);
591 AliMatrix(idrotm[857],90.0,5.0,90.0,95.0,180.0,0.0);
592 AliMatrix(idrotm[858],90.0,65.0,90.0,155.0,180.0,0.0);
593 AliMatrix(idrotm[859],90.0,305.0,90.0,35.0,180.0,0.0);
594 AliMatrix(idrotm[860],90.0,245.0,90.0,335.0,180.0,0.0);
595 AliMatrix(idrotm[861],90.0,185.0,90.0,275.0,180.0,0.0);
596 AliMatrix(idrotm[862],90.0,125.0,90.0,215.0,180.0,0.0);
597 AliMatrix(idrotm[863],90.0,257.5,90.0,347.5,180.0,0.0);
598 AliMatrix(idrotm[864],90.0,227.5,90.0,317.5,180.0,0.0);
599 AliMatrix(idrotm[865],90.0,197.5,90.0,287.5,180.0,0.0);
600 AliMatrix(idrotm[867],90.0,167.5,90.0,257.5,180.0,0.0);
601 AliMatrix(idrotm[868],90.0,287.5,90.0,17.5,0.0,0.0);
602 AliMatrix(idrotm[869],90.0,137.5,90.0,227.5,180.0,0.0);
603 AliMatrix(idrotm[870],90.0,107.5,90.0,197.5,180.0,0.0);
604 AliMatrix(idrotm[871],90.0,77.5,90.0,167.5,180.0,0.0);
605 AliMatrix(idrotm[872],90.0,47.5,90.0,137.5,180.0,0.0);
606 AliMatrix(idrotm[873],90.0,17.5,90.0,107.5,180.0,0.0);
607 AliMatrix(idrotm[874],90.0,347.5,90.0,77.5,180.0,0.0);
608 AliMatrix(idrotm[875],90.0,317.5,90.0,47.5,180.0,0.0);
609 AliMatrix(idrotm[876],90.0,287.5,90.0,17.5,180.0,0.0);
610 AliMatrix(idrotm[877],90.0,185.0,90.0,275.0,0.0,0.0);
611 AliMatrix(idrotm[878],90.0,180.0,90.0,270.0,0.0,0.0);
612 AliMatrix(idrotm[879],90.0,125.0,90.0,215.0,0.0,0.0);
613 AliMatrix(idrotm[880],90.0,65.0,90.0,155.0,0.0,0.0);
614 AliMatrix(idrotm[881],90.0,5.0,90.0,95.0,0.0,0.0);
615 AliMatrix(idrotm[882],90.0,245.0,90.0,335.0,0.0,0.0);
616 AliMatrix(idrotm[883],90.0,47.5,90.0,137.5,0.0,0.0);
617 AliMatrix(idrotm[884],90.0,77.5,90.0,167.5,0.0,0.0);
618 AliMatrix(idrotm[885],90.0,107.5,90.0,197.5,0.0,0.0);
619 AliMatrix(idrotm[887],90.0,137.5,90.0,227.5,0.0,0.0);
620 AliMatrix(idrotm[888],90.0,167.5,90.0,257.5,0.0,0.0);
621 AliMatrix(idrotm[889],90.0,197.5,90.0,287.5,0.0,0.0);
622 AliMatrix(idrotm[890],90.0,227.5,90.0,317.5,0.0,0.0);
623 AliMatrix(idrotm[891],90.0,347.5,90.0,77.5,0.0,0.0);
624 AliMatrix(idrotm[892],90.0,317.5,90.0,47.5,0.0,0.0);
625 AliMatrix(idrotm[893],90.0,257.5,90.0,347.5,0.0,0.0);
626 AliMatrix(idrotm[894],90.0,270.0,0.0,0.0,90.0,180.0);
627 AliMatrix(idrotm[895],90.0,286.36,0.0,0.0,90.0,196.36);
628 AliMatrix(idrotm[896],90.0,302.73,0.0,0.0,90.0,212.73);
629 AliMatrix(idrotm[897],90.0,319.09,0.0,0.0,90.0,229.09);
630 AliMatrix(idrotm[898],90.0,17.5,90.0,107.5,0.0,0.0);
631 AliMatrix(idrotm[899],90.0,335.45,0.0,0.0,90.0,245.45);
632 AliMatrix(idrotm[900],90.0,351.82,0.0,0.0,90.0,261.82);
633 AliMatrix(idrotm[901],90.0,8.18,0.0,0.0,90.0,278.18);
634 AliMatrix(idrotm[902],90.0,24.55,0.0,0.0,90.0,294.55);
635 AliMatrix(idrotm[903],90.0,40.91,0.0,0.0,90.0,310.91);
636 AliMatrix(idrotm[904],90.0,57.27,0.0,0.0,90.0,327.27);
637 AliMatrix(idrotm[905],90.0,73.64,0.0,0.0,90.0,343.64);
638 AliMatrix(idrotm[906],90.0,90.0,0.0,0.0,90.0,360.0);
639 AliMatrix(idrotm[907],90.0,106.36,0.0,0.0,90.0,16.36);
640 AliMatrix(idrotm[908],90.0,122.73,0.0,0.0,90.0,32.73);
641 AliMatrix(idrotm[909],90.0,139.09,0.0,0.0,90.0,49.09);
642 AliMatrix(idrotm[910],90.0,155.45,0.0,0.0,90.0,65.45);
643 AliMatrix(idrotm[911],90.0,171.82,0.0,0.0,90.0,81.82);
644 AliMatrix(idrotm[912],90.0,188.18,0.0,0.0,90.0,98.18);
645 AliMatrix(idrotm[913],90.0,204.55,0.0,0.0,90.0,114.55);
646 AliMatrix(idrotm[914],90.0,220.91,0.0,0.0,90.0,130.91);
647 AliMatrix(idrotm[915],90.0,237.27,0.0,0.0,90.0,147.27);
648 AliMatrix(idrotm[916],90.0,253.64,0.0,0.0,90.0,163.64);
649 AliMatrix(idrotm[917],90.0,295.71,0.0,0.0,90.0,205.71);
650 AliMatrix(idrotm[918],90.0,321.43,0.0,0.0,90.0,231.43);
651 AliMatrix(idrotm[919],90.0,347.14,0.0,0.0,90.0,257.14);
652 AliMatrix(idrotm[920],90.0,12.86,0.0,0.0,90.0,282.86);
653 AliMatrix(idrotm[921],90.0,38.57,0.0,0.0,90.0,308.57);
654 AliMatrix(idrotm[922],90.0,64.29,0.0,0.0,90.0,334.29);
655 AliMatrix(idrotm[923],90.0,115.71,0.0,0.0,90.0,25.71);
656 AliMatrix(idrotm[924],90.0,141.43,0.0,0.0,90.0,51.43);
657 AliMatrix(idrotm[925],90.0,167.14,0.0,0.0,90.0,77.14);
658 AliMatrix(idrotm[926],90.0,192.86,0.0,0.0,90.0,102.86);
659 AliMatrix(idrotm[927],90.0,218.57,0.0,0.0,90.0,128.57);
660 AliMatrix(idrotm[928],90.0,244.29,0.0,0.0,90.0,154.29);
661 AliMatrix(idrotm[929],90.0,120.0,90.0,210.0,0.0,0.0);
662 AliMatrix(idrotm[930],90.0,240.0,90.0,330.0,0.0,0.0);
663 AliMatrix(idrotm[931],90.0,60.0,90.0,150.0,180.0,0.0);
664 AliMatrix(idrotm[932],90.0,120.0,90.0,210.0,180.0,0.0);
665 AliMatrix(idrotm[933],90.0,180.0,90.0,270.0,180.0,0.0);
666 AliMatrix(idrotm[934],90.0,240.0,90.0,330.0,180.0,0.0);
667 AliMatrix(idrotm[935],90.0,300.0,90.0,30.0,180.0,0.0);
668
669 // SSD cone
670
671 AliMatrix(idrotm[701],90.0,0.0,90.0,90.0,180.0,0.0);
672 AliMatrix(idrotm[702],90.0,347.5,90.0,77.5,180.0,0.0);
673 AliMatrix(idrotm[703],90.0,17.5,90.0,107.5,180.0,0.0);
674 AliMatrix(idrotm[704],90.0,47.5,90.0,137.5,180.0,0.0);
675 AliMatrix(idrotm[705],90.0,77.5,90.0,167.5,180.0,0.0);
676 AliMatrix(idrotm[706],90.0,107.5,90.0,197.5,180.0,0.0);
677 AliMatrix(idrotm[707],90.0,137.5,90.0,227.5,180.0,0.0);
678 AliMatrix(idrotm[708],90.0,167.5,90.0,257.5,180.0,0.0);
679 AliMatrix(idrotm[709],90.0,197.5,90.0,287.5,180.0,0.0);
680 AliMatrix(idrotm[710],90.0,227.5,90.0,317.5,180.0,0.0);
681 AliMatrix(idrotm[711],90.0,257.5,90.0,347.5,180.0,0.0);
682 AliMatrix(idrotm[712],90.0,287.5,90.0,17.5,180.0,0.0);
683 AliMatrix(idrotm[713],90.0,317.5,90.0,47.5,180.0,0.0);
684 AliMatrix(idrotm[714],90.0,328.4,90.0,58.4,180.0,0.0);
685 AliMatrix(idrotm[715],90.0,28.4,90.0,118.4,180.0,0.0);
686 AliMatrix(idrotm[716],90.0,88.4,90.0,178.4,180.0,0.0);
687 AliMatrix(idrotm[717],90.0,148.4,90.0,238.4,180.0,0.0);
688 AliMatrix(idrotm[718],90.0,208.4,90.0,298.4,180.0,0.0);
689 AliMatrix(idrotm[719],90.0,268.4,90.0,358.4,180.0,0.0);
690 AliMatrix(idrotm[720],90.0,28.4,90.0,118.4,0.0,0.0);
691 AliMatrix(idrotm[721],90.0,88.4,90.0,178.4,0.0,0.0);
692 AliMatrix(idrotm[722],90.0,148.4,90.0,238.4,0.0,0.0);
693 AliMatrix(idrotm[723],90.0,208.4,90.0,298.4,0.0,0.0);
694 AliMatrix(idrotm[724],90.0,268.4,90.0,358.4,0.0,0.0);
695 AliMatrix(idrotm[725],90.0,328.4,90.0,58.4,0.0,0.0);
696 AliMatrix(idrotm[726],90.0,77.5,90.0,167.5,0.0,0.0);
697 AliMatrix(idrotm[727],90.0,107.5,90.0,197.5,0.0,0.0);
698 AliMatrix(idrotm[728],90.0,137.5,90.0,227.5,0.0,0.0);
699 AliMatrix(idrotm[729],90.0,167.5,90.0,257.5,0.0,0.0);
700 AliMatrix(idrotm[730],90.0,227.5,90.0,317.5,0.0,0.0);
701 AliMatrix(idrotm[731],90.0,257.5,90.0,347.5,0.0,0.0);
702 AliMatrix(idrotm[732],90.0,317.5,90.0,47.5,0.0,0.0);
703 AliMatrix(idrotm[733],90.0,197.5,90.0,287.5,0.0,0.0);
704 AliMatrix(idrotm[734],90.0,347.5,90.0,77.5,0.0,0.0);
705 AliMatrix(idrotm[735],90.0,47.5,90.0,137.5,0.0,0.0);
706 AliMatrix(idrotm[768],90.0,287.5,90.0,17.5,0.0,0.0);
707 AliMatrix(idrotm[798],90.0,17.5,90.0,107.5,0.0,0.0);
708
709 // Services
710
711 AliMatrix(idrotm[200], 90., 0., 90., 90., 180., 0.);
db915fea 712
713
714 // CONVERT INTO CM (RL(SI)=9.36 CM)
9d2da7b3 715
db915fea 716 for (i = 0; i < 6; ++i) {
717 drl[i] = drl[i] / 100. * 9.36;
718 }
719
720 // FIELD CAGE HALF LENGTH
721
1e9111db 722 rlim = 50.;
e4f08ac3 723 zmax = 74.;
724 ztpc = 284.;
9d2da7b3 725
db915fea 726 // --- Define ghost volume containing the whole ITS (including services)
727 // and fill it with air
728
729 dgh[0] = 0.;
730 dgh[1] = 360.;
6d62fed9 731 dgh[2] = 16.;
e4f08ac3 732 dgh[3] = -ztpc-5.-0.1;
733 dgh[4] = 62.4;
db915fea 734 dgh[5] = 85.;
e4f08ac3 735 dgh[6] = -ztpc;
736 dgh[7] = 62;
db915fea 737 dgh[8] = 85.;
e4f08ac3 738 dgh[9] = -ztpc;
739 dgh[10] = 62;
740 dgh[11] = 62+4.;
741 dgh[12] = -97.5;
742 dgh[13] = 46;
1e9111db 743 dgh[14] = rlim+6;
e4f08ac3 744 dgh[15] = -zmax;
745 dgh[16] = 46;
1e9111db 746 dgh[17] = rlim+6;
e4f08ac3 747 dgh[18] = -48;
6d62fed9 748 dgh[19] = 6;
1e9111db 749 dgh[20] = rlim+6;
e4f08ac3 750 dgh[21] = -28.6;
6d62fed9 751 dgh[22] = 6;
1e9111db 752 dgh[23] = rlim+6;
e4f08ac3 753 dgh[24] = -27.6;
6d62fed9 754 dgh[25] = 3.295;
1e9111db 755 dgh[26] = rlim+6;
e4f08ac3 756 dgh[27] = 27.6;
6d62fed9 757 dgh[28] = 3.295;
1e9111db 758 dgh[29] = rlim+6;
e4f08ac3 759 dgh[30] = 28.6;
6d62fed9 760 dgh[31] = 6;
1e9111db 761 dgh[32] = rlim+6;
e4f08ac3 762 dgh[33] = 48;
6d62fed9 763 dgh[34] = 6;
1e9111db 764 dgh[35] = rlim+6;
e4f08ac3 765 dgh[36] = zmax;
6d62fed9 766 dgh[37] = 46;
1e9111db 767 dgh[38] = rlim+6;
e4f08ac3 768 dgh[39] = 97.5;
769 dgh[40] = 46;
1e9111db 770 dgh[41] = rlim+6;
e4f08ac3 771 dgh[42] = ztpc;
772 dgh[43] = 62;
773 dgh[44] = 62+4.;
774 dgh[45] = ztpc;
775 dgh[46] = 62;
6d62fed9 776 dgh[47] = 85.;
e4f08ac3 777 dgh[48] = ztpc+4.+0.1;
6d62fed9 778 dgh[49] = 62.4;
779 dgh[50] = 85.;
780 gMC->Gsvolu("ITSV", "PCON", idtmed[205], dgh, 51);
45725ec9 781
6d62fed9 782
db915fea 783
784 // --- Place the ghost volume in its mother volume (ALIC) and make it
785 // invisible
786
787 gMC->Gspos("ITSV", 1, "ALIC", 0., 0., 0., 0, "ONLY");
788 //gMC->Gsatt("ITSV", "SEEN", 0);
9d2da7b3 789
790
db915fea 791 // --- Define ghost volume containing the six layers and fill it with air
792
793 dgh[0] = 0.;
794 dgh[1] = 360.;
6d62fed9 795 dgh[2] = 8.;
e4f08ac3 796 dgh[3] = -zmax;
45725ec9 797 dgh[4] = 46.;
e4f08ac3 798 dgh[5] = rlim;
799 dgh[6] = -47.5;
6d62fed9 800 dgh[7] = 6.005;
e4f08ac3 801 dgh[8] = rlim;
802 dgh[9] = -28.5;
6d62fed9 803 dgh[10] = 6.005;
e4f08ac3 804 dgh[11] = rlim;
805 dgh[12] = -27.5;
6d62fed9 806 dgh[13] = 3.3;
e4f08ac3 807 dgh[14] = rlim;
808 dgh[15] = 27.5;
6d62fed9 809 dgh[16] = 3.3;
e4f08ac3 810 dgh[17] = rlim;
811 dgh[18] = 28.5;
6d62fed9 812 dgh[19] = 6.005;
e4f08ac3 813 dgh[20] = rlim;
814 dgh[21] = 47.5;
6d62fed9 815 dgh[22] = 6.005;
e4f08ac3 816 dgh[23] = rlim;
817 dgh[24] = zmax;
6d62fed9 818 dgh[25] = 46.;
e4f08ac3 819 dgh[26] = rlim;
6d62fed9 820 gMC->Gsvolu("ITSD", "PCON", idtmed[205], dgh, 27);
db915fea 821
9d2da7b3 822 // --- Place the ghost volume in its mother volume (ITSV) and make it
db915fea 823 // invisible
824
825 gMC->Gspos("ITSD", 1, "ITSV", 0., 0., 0., 0, "ONLY");
826 //gMC->Gsatt("ITSD", "SEEN", 0);
db915fea 827
9d2da7b3 828 // --- Define SPD (option 'a') volumes ----------------------------
829
830 // SPD - option 'a'
831 // (this is NOT the default)
832
833 if (option == 1) {
834
835 dits[0] = 3.7;
836 dits[1] = 7.75;
45725ec9 837 dits[2] = 26.1;
9d2da7b3 838 gMC->Gsvolu("IT12", "TUBE", idtmed[254], dits, 3);
839
840 dits[0] = 3.7;
841 dits[1] = 7.7;
842 dits[2] = 24;
843 dits[3] = 57;
844 dits[4] = 100;
845 gMC->Gsvolu("I12A", "TUBS", idtmed[254], dits, 5); // sector
846
847 di10a[0] = 0.843;
848 di10a[1] = ddet1+dchip1+dbus+0.0025;
849 di10a[2] = 19.344;
850 gMC->Gsvolu("I10A", "BOX ", idtmed[254], di10a, 3); // mother volume
851 // on layer 1
852 di20a[0] = 0.843;
853 di20a[1] = ddet2+dchip2+dbus+0.0025;
854 di20a[2] = 19.344;
855 gMC->Gsvolu("I20A", "BOX ", idtmed[254], di20a, 3); // mother volume
856 // on layer 2
857 dits[0] = 1.3673;
858 dits[1] = 0.01;
859 dits[2] = 24;
860 gMC->Gsvolu("I123", "BOX ", idtmed[253], dits, 3);
861
862 dits[0] = 0.06;
863 dits[1] = 0.08;
864 dits[2] = 24;
865 dits[3] = -36.79;
866 dits[4] = 21.834;
867 gMC->Gsvolu("I121", "TUBS", idtmed[253], dits, 5);
868
869 dits[0] = 0.1253;
870 dits[1] = 0.01;
871 dits[2] = 24;
872 gMC->Gsvolu("I122", "BOX ", idtmed[253], dits, 3);
873
874 dits[0] = 0.04;
875 dits[1] = 0.06 ;
876 dits[2] = 24;
877 dits[3] = 126.79;
878 dits[4] = 270;
879 gMC->Gsvolu("I120", "TUBS", idtmed[253], dits, 5);
880
881 dits[0] = 0.1134;
882 dits[1] = 0.01;
883 dits[2] = 24;
884 gMC->Gsvolu("I144", "BOX ", idtmed[253], dits, 3);
885
886 dits[0] = 0.25;
887 dits[1] = 0.06;
888 dits[2] = 24;
889 gMC->Gsvolu("I113", "BOX ", idtmed[254], dits, 3);
890
891 dits[0] = 0.077;
892 dits[1] = 0.01;
893 dits[2] = 24;
894 gMC->Gsvolu("I143", "BOX ", idtmed[253], dits, 3);
895
896 dits[0] = 0.04;
897 dits[1] = 0.06;
898 dits[2] = 24;
899 dits[3] = 0;
900 dits[4] = 90;
901 gMC->Gsvolu("I142", "TUBS", idtmed[253], dits, 5);
902
903 dits[0] = 0.0695;
904 dits[1] = 0.01;
905 dits[2] = 24;
906 gMC->Gsvolu("I141", "BOX ", idtmed[253], dits, 3);
907
908 dits[0] = 0.06;
909 dits[1] = 0.08;
910 dits[2] = 24;
911 dits[3] = 0;
912 dits[4] = 108;
913 gMC->Gsvolu("I140", "TUBS", idtmed[253], dits, 5);
914
915 dits[0] = 0.1835;
916 dits[1] = 0.01;
917 dits[2] = 24;
918 gMC->Gsvolu("I139", "BOX ", idtmed[253], dits, 3);
919
920 dits[0] = 0.1894 ;
921 dits[1] = 0.01;
922 dits[2] = 24;
923 gMC->Gsvolu("I138", "BOX ", idtmed[253], dits, 3);
924
925 dits[0] = 0.04;
926 dits[1] = 0.06;
927 dits[2] = 24;
928 dits[3] = 0;
929 dits[4] = 75.261;
930 gMC->Gsvolu("I137", "TUBS", idtmed[253], dits, 5);
931
932 dits[0] = 1.3401;
933 dits[1] = 0.01;
934 dits[2] = 24;
935 gMC->Gsvolu("I136", "BOX ", idtmed[253], dits, 3);
936
937 dits[0] = 0.05;
938 dits[1] = 0.07;
939 dits[2] = 24;
940 dits[3] = 0;
941 dits[4] = 72.739;
942 gMC->Gsvolu("I135", "TUBS", idtmed[253], dits, 5);
943
944 dits[0] = 0.1193;
945 dits[1] = 0.01;
946 dits[2] = 24;
947 gMC->Gsvolu("I134", "BOX ", idtmed[253], dits, 3);
948
949 dits[0] = 0.163;
950 dits[1] = 0.01;
951 dits[2] = 24;
952 gMC->Gsvolu("I133", "BOX ", idtmed[253], dits, 3);
953
954 dits[0] = 0.04;
955 dits[1] = 0.06;
956 dits[2] = 24;
957 dits[3] = 0;
958 dits[4] = 157.633;
959 gMC->Gsvolu("I132", "TUBS", idtmed[253], dits, 5);
960
961 dits[0] = 0.2497;
962 dits[1] = 0.01;
963 dits[2] = 24;
964 gMC->Gsvolu("I131", "BOX ", idtmed[253], dits, 3);
965
966 dits[0] = 0.06;
967 dits[1] = 0.08;
968 dits[2] = 24;
969 dits[3] = 0;
970 dits[4] = 148.633;
971 gMC->Gsvolu("I130", "TUBS", idtmed[253], dits, 5);
972
973 dits[0] = 0.292;
974 dits[1] = 0.01;
975 dits[2] = 24;
976 gMC->Gsvolu("I129", "BOX ", idtmed[253], dits, 3);
977
978 dits[0] = 0.163;
979 dits[1] = 0.01;
980 dits[2] = 24;
981 gMC->Gsvolu("I128", "BOX ", idtmed[253], dits, 3);
982
983 dits[0] = 0.04;
984 dits[1] = 0.06;
985 dits[2] = 24;
986 dits[3] = 0;
987 dits[4] = 161.297;
988 gMC->Gsvolu("I126", "TUBS", idtmed[253], dits, 5);
989
990 dits[0] = 0.2433;
991 dits[1] = 0.01;
992 dits[2] = 24;
993 gMC->Gsvolu("I125", "BOX ", idtmed[253], dits, 3);
994
995 dits[0] = 0.06;
996 dits[1] = 0.08;
997 dits[2] = 24;
998 dits[3] = 0;
999 dits[4] = 42.883;
1000 gMC->Gsvolu("I124", "TUBS", idtmed[253], dits, 5);
1001
1002 di103[0] = 0.793;
1003 di103[1] = ddet1+dchip1;
1004 di103[2] = 3.536;
1005 gMC->Gsvolu("I103", "BOX ", idtmed[254], di103, 3); // contains det and chip
1006 // layer 1
1007 dits[0] = 0.793;
fc44ab90 1008 dits[1] = 0.475; //0.685; 0.015
9d2da7b3 1009 dits[2] = 2.5;
fc44ab90 1010 gMC->Gsvolu("I105", "BOX ", idtmed[290], dits, 3);
9d2da7b3 1011
1012 di104[0] = 0.843;
1013 di104[1] = dbus;
1014 di104[2] = 14.344;
1015 gMC->Gsvolu("I104", "BOX ", idtmed[275], di104, 3); // bus for both layers
1016
1017 di1d3[0] = 0.793;
1018 di1d3[1] = ddet2+dchip2;
1019 di1d3[2] = 3.536;
1020 gMC->Gsvolu("I1D3", "BOX ", idtmed[254], di1d3, 3); // contains det and chip
1021 // layer 2
1022 dits[0] = 0.793;
1023 dits[0] = 0.06;
1024 dits[1] = 0.08;
1025 dits[2] = 24;
1026 dits[3] = 0;
1027 dits[4] = 80;
1028 gMC->Gsvolu("I112", "TUBS", idtmed[253], dits, 5);
1029
1030 dits[0] = 0.04;
1031 dits[1] = 0.06;
1032 dits[2] = 24;
1033 dits[3] = 0;
1034 dits[4] = 80;
1035 gMC->Gsvolu("I111", "TUBS", idtmed[253], dits, 5);
1036
1037 dits[0] = 0.15;
1038 dits[1] = 0.0146;
1039 dits[2] = 24;
1040 gMC->Gsvolu("I118", "BOX ", idtmed[273], dits, 3);
1041
1042 dits[0] = 0.1315;
1043 dits[1] = 0.01;
1044 dits[2] = 24;
1045 gMC->Gsvolu("I110", "BOX ", idtmed[253], dits, 3);
1046
1047 dits[0] = 0.025;
1048 dits[1] = 0.035;
1049 dits[2] = 24;
1050 dits[3] = 0;
1051 dits[4] = 180;
1052 gMC->Gsvolu("I114", "TUBS", idtmed[264], dits, 5);
1053
e4f08ac3 1054 if (fluid == 1) {
1055 dits[0] = 0;
1056 dits[1] = 0.025;
1057 dits[2] = 24;
1058 dits[3] = 0;
1059 dits[4] = 180;
1060 gMC->Gsvolu("I115", "TUBS", idtmed[211], dits, 5); // set water as cooling fluid
1061 } else {
1062 dits[0] = 0;
1063 dits[1] = 0.025;
1064 dits[2] = 24;
1065 dits[3] = 0;
1066 dits[4] = 180;
1067 gMC->Gsvolu("I115", "TUBS", idtmed[212], dits, 5); // set freon as cooling fluid
1068 }
1069
9d2da7b3 1070 dits[0] = 0.063;
1071 dits[1] = 0.035;
1072 dits[2] = 24;
1073 gMC->Gsvolu("I116", "BOX ", idtmed[264], dits, 3);
1074
1075 di102[0] = 0.793;
1076 di102[1] = dchip1;
1077 di102[2] = 0.68;
1078 gMC->Gsvolu("I102", "BOX ", idtmed[201], di102, 3); // chip layer 1
1079
1080 di1d2[0] = 0.793;
1081 di1d2[1] = dchip2;
1082 di1d2[2] = 0.68;
1083 gMC->Gsvolu("I1D2", "BOX ", idtmed[201], di1d2, 3); // chip layer 2
1084
1085 di101[0] = 0.705;
1086 di101[1] = ddet1;
1087 di101[2] = 3.536;
1088 gMC->Gsvolu("I101", "BOX ", idtmed[250], di101, 3); // contains detector
1089 // layer 1
1090 di1d1[0] = 0.705;
1091 di1d1[1] = ddet2;
1092 di1d1[2] = 3.536;
1093 gMC->Gsvolu("I1D1", "BOX ", idtmed[250], di1d1, 3); // contains detector
1094 // layer 2
e4f08ac3 1095 if (fluid == 1) {
1096 dits[0] = 0.063;
1097 dits[1] = 0.025;
1098 dits[2] = 24;
1099 gMC->Gsvolu("I117", "BOX ", idtmed[211], dits, 3); // set water as cooling fuid
1100 } else {
1101 dits[0] = 0.063;
1102 dits[1] = 0.025;
1103 dits[2] = 24;
1104 gMC->Gsvolu("I117", "BOX ", idtmed[212], dits, 3); // set freon as cooling fluid
1105 }
9d2da7b3 1106
1107 dits1[0] = 0.64;
1108 dits1[1] = ddet1;
1109 dits1[2] = 3.48;
1110 gMC->Gsvolu("ITS1", "BOX ", idtmed[200], dits1, 3); // detector layer 1
1111
1112 dits2[0] = 0.64;
1113 dits2[1] = ddet2;
1114 dits2[2] = 3.48;
1115 gMC->Gsvolu("ITS2", "BOX ", idtmed[200], dits2, 3); // detector layer 2
1116
1117 dits[0] = 3.701;
1118 dits[1] = 7.699;
1119 dits[2] = 4;
1120 dits[3] = 57.1;
1121 dits[4] = 99.9;
1122 gMC->Gsvolu("I650", "TUBS", idtmed[254], dits, 5); // was I150 in old geom.
1123
1124 dits[0] = 0;
1125 dits[1] = 0.5;
1126 dits[2] = 1.5;
1127 gMC->Gsvolu("I676", "TUBE", idtmed[274], dits, 3); // was I176 in old geom.
1128
1129 dits[0] = 0;
1130 dits[1] = 0.18;
1131 dits[2] = 0.8;
1132 gMC->Gsvolu("I673", "TUBE", idtmed[274], dits, 3); // was I173 in old geom.
1133
1134 dits[0] = 0;
1135 dits[1] = 0.18;
1136 dits[2] = 3;
1137 gMC->Gsvolu("I671", "TUBE", idtmed[274], dits, 3); // was I171 in old geom.
1138
1139 dits[0] = 0;
1140 dits[1] = 0.075;
1141 dits[2] = 0.8;
1142 gMC->Gsvolu("I669", "TUBE", idtmed[264], dits, 3); // was I169 in old geom.
1143
1144 dits[0] = 3.5;
1145 dits[1] = 5.6;
1146 dits[2] = 0.55;
1147 dits[3] = 0;
1148 dits[4] = 38;
1149 gMC->Gsvolu("I667", "TUBS", idtmed[263], dits, 5); // was I167 in old geom.
1150
1151 dits[0] = 6.6;
1152 dits[1] = 7.6;
1153 dits[2] = 0.5;
1154 dits[3] = 0;
1155 dits[4] = 9;
1156 gMC->Gsvolu("I666", "TUBS", idtmed[263], dits, 5); // was I166 in old geom.
1157
1158 dits[0] = 0.26;
1159 dits[1] = 0.32;
1160 dits[2] = 0.55;
1161 gMC->Gsvolu("I678", "TUBE", idtmed[263], dits, 3); // was I178 in old geom.
1162
e4f08ac3 1163 if (fluid == 1) {
1164 dits[0] = 0;
1165 dits[1] = 0.3;
1166 dits[2] = 1.5;
1167 gMC->Gsvolu("I677", "TUBE", idtmed[211], dits, 3); // set water as cooling fluid
1168 // was I177 in old geom.
1169 } else {
1170 dits[0] = 0;
1171 dits[1] = 0.3;
1172 dits[2] = 1.5;
1173 gMC->Gsvolu("I677", "TUBE", idtmed[212], dits, 3); // set freon as cooling fluid
1174 // was I177 in old geom.
1175 }
1176
9d2da7b3 1177 dits[0] = 0.07;
1178 dits[1] = 0.125;
1179 dits[2] = 0.3;
1180 gMC->Gsvolu("I675", "TUBE", idtmed[263], dits, 3); // was I175 in old geom.
1181
e4f08ac3 1182 if (fluid == 1) {
1183 dits[0] = 0;
1184 dits[1] = 0.1;
1185 dits[2] = 0.8;
1186 gMC->Gsvolu("I674", "TUBE", idtmed[211], dits, 3); // set water as cooling fluid
1187 // was I174 in old geom.
1188 } else {
1189 dits[0] = 0;
1190 dits[1] = 0.1;
1191 dits[2] = 0.8;
1192 gMC->Gsvolu("I674", "TUBE", idtmed[212], dits, 3); // set freon as cooling fluid
1193 // was I174 in old geom.
1194 }
1195
1196 if (fluid == 1) {
1197 dits[0] = 0;
1198 dits[1] = 0.1;
1199 dits[2] = 3;
1200 gMC->Gsvolu("I672", "TUBE", idtmed[211], dits, 3); // set water as cooling fluid
1201 // was I172 in old geom.
1202 } else {
1203 dits[0] = 0;
1204 dits[1] = 0.1;
1205 dits[2] = 3;
1206 gMC->Gsvolu("I672", "TUBE", idtmed[212], dits, 3); // set freon as cooling fluid
1207 // was I172 in old geom.
1208 }
1209
1210 if (fluid == 1) {
1211 dits[0] = 0;
1212 dits[1] = 0.0746;
1213 dits[2] = 0.8;
1214 gMC->Gsvolu("I670", "TUBE", idtmed[211], dits, 3); // set water as cooling fluid
1215 // was I170 in old geom.
1216 } else {
1217 dits[0] = 0;
1218 dits[1] = 0.0746;
1219 dits[2] = 0.8;
1220 gMC->Gsvolu("I670", "TUBE", idtmed[212], dits, 3); // set freon as cooling fluid
1221 // was I170 in old geom.
1222 }
1223
1224 if (fluid == 1) {
1225 dits[0] = 3.7;
1226 dits[1] = 5.4;
1227 dits[2] = 0.35;
1228 dits[3] = 2;
1229 dits[4] = 36;
1230 gMC->Gsvolu("I668", "TUBS", idtmed[211], dits, 5); // set water as cooling fluid
1231 // was I168 in old geom.
1232 } else {
1233 dits[0] = 3.7;
1234 dits[1] = 5.4;
1235 dits[2] = 0.35;
1236 dits[3] = 2;
1237 dits[4] = 36;
1238 gMC->Gsvolu("I668", "TUBS", idtmed[212], dits, 5); // set freon as cooling fluid
1239 // was I168 in old geom.
1240 }
9d2da7b3 1241
db915fea 1242
9d2da7b3 1243 }
db915fea 1244
9d2da7b3 1245 // --- Define SPD (option 'b') volumes ----------------------------
1246
1247 // SPD - option 'b'
1248 // (this is the default)
1249
1250 if (option == 2) {
1251
1252 dits[0] = 3.7;
1253 dits[1] = 7.75;
45725ec9 1254 dits[2] = 26.1;
9d2da7b3 1255 gMC->Gsvolu("IT12", "TUBE", idtmed[254], dits, 3);
1256
1257 dits[0] = 3.7;
1258 dits[1] = 7.7;
1259 dits[2] = 24;
1260 dits[3] = 57;
1261 dits[4] = 100;
1262 gMC->Gsvolu("I12B", "TUBS", idtmed[254], dits, 5); // sector
1263
1264 di10b[0] = 0.843;
1265 di10b[1] = ddet1+dchip1+dbus+0.0025;
1266 di10b[2] = 19.344;
1267 gMC->Gsvolu("I10B", "BOX ", idtmed[254], di10b, 3); // mother volume
1268 // on layer 1
1269
1270 di20b[0] = 0.843;
1271 di20b[1] = ddet2+dchip2+dbus+0.0025;
1272 di20b[2] = 19.344;
1273 gMC->Gsvolu("I20B", "BOX ", idtmed[254], di20b, 3); // mother volume
1274 // layer 2
1275
1276 dits[0] = 1.3673;
1277 dits[1] = 0.01;
1278 dits[2] = 24;
1279 gMC->Gsvolu("I123", "BOX ", idtmed[253], dits, 3);
1280
1281 dits[0] = 0.06;
1282 dits[1] = 0.08;
1283 dits[2] = 24;
1284 dits[3] = -36.79;
1285 dits[4] = 21.834;
1286 gMC->Gsvolu("I121", "TUBS", idtmed[253], dits, 5);
1287
1288 dits[0] = 0.1253;
1289 dits[1] = 0.01;
1290 dits[2] = 24;
1291 gMC->Gsvolu("I122", "BOX ", idtmed[253], dits, 3);
1292
1293 dits[0] = 0.04;
1294 dits[1] = 0.06 ;
1295 dits[2] = 24;
1296 dits[3] = 126.79;
1297 dits[4] = 270;
1298 gMC->Gsvolu("I120", "TUBS", idtmed[253], dits, 5);
1299
1300 dits[0] = 0.1134;
1301 dits[1] = 0.01;
1302 dits[2] = 24;
1303 gMC->Gsvolu("I144", "BOX ", idtmed[253], dits, 3);
1304
1305 dits[0] = 0.25;
1306 dits[1] = 0.06;
1307 dits[2] = 24;
1308 gMC->Gsvolu("I113", "BOX ", idtmed[254], dits, 3);
1309
1310 dits[0] = 0.077;
1311 dits[1] = 0.01;
1312 dits[2] = 24;
1313 gMC->Gsvolu("I143", "BOX ", idtmed[253], dits, 3);
1314
1315 dits[0] = 0.04;
1316 dits[1] = 0.06;
1317 dits[2] = 24;
1318 dits[3] = 0;
1319 dits[4] = 90;
1320 gMC->Gsvolu("I142", "TUBS", idtmed[253], dits, 5);
1321
1322 dits[0] = 0.0695;
1323 dits[1] = 0.01;
1324 dits[2] = 24;
1325 gMC->Gsvolu("I141", "BOX ", idtmed[253], dits, 3);
1326
1327 dits[0] = 0.06;
1328 dits[1] = 0.08;
1329 dits[2] = 24;
1330 dits[3] = 0;
1331 dits[4] = 108;
1332 gMC->Gsvolu("I140", "TUBS", idtmed[253], dits, 5);
1333
1334 dits[0] = 0.1835;
1335 dits[1] = 0.01;
1336 dits[2] = 24;
1337 gMC->Gsvolu("I139", "BOX ", idtmed[253], dits, 3);
1338
1339 dits[0] = 0.1894 ;
1340 dits[1] = 0.01;
1341 dits[2] = 24;
1342 gMC->Gsvolu("I138", "BOX ", idtmed[253], dits, 3);
1343
1344 dits[0] = 0.04;
1345 dits[1] = 0.06;
1346 dits[2] = 24;
1347 dits[3] = 0;
1348 dits[4] = 75.261;
1349 gMC->Gsvolu("I137", "TUBS", idtmed[253], dits, 5);
1350
1351 dits[0] = 1.3401;
1352 dits[1] = 0.01;
1353 dits[2] = 24;
1354 gMC->Gsvolu("I136", "BOX ", idtmed[253], dits, 3);
1355
1356 dits[0] = 0.05;
1357 dits[1] = 0.07;
1358 dits[2] = 24;
1359 dits[3] = 0;
1360 dits[4] = 72.739;
1361 gMC->Gsvolu("I135", "TUBS", idtmed[253], dits, 5);
1362
1363 dits[0] = 0.1193;
1364 dits[1] = 0.01;
1365 dits[2] = 24;
1366 gMC->Gsvolu("I134", "BOX ", idtmed[253], dits, 3);
1367
1368 dits[0] = 0.163;
1369 dits[1] = 0.01;
1370 dits[2] = 24;
1371 gMC->Gsvolu("I133", "BOX ", idtmed[253], dits, 3);
1372
1373 dits[0] = 0.04;
1374 dits[1] = 0.06;
1375 dits[2] = 24;
1376 dits[3] = 0;
1377 dits[4] = 157.633;
1378 gMC->Gsvolu("I132", "TUBS", idtmed[253], dits, 5);
1379
1380 dits[0] = 0.2497;
1381 dits[1] = 0.01;
1382 dits[2] = 24;
1383 gMC->Gsvolu("I131", "BOX ", idtmed[253], dits, 3);
1384
1385 dits[0] = 0.06;
1386 dits[1] = 0.08;
1387 dits[2] = 24;
1388 dits[3] = 0;
1389 dits[4] = 148.633;
1390 gMC->Gsvolu("I130", "TUBS", idtmed[253], dits, 5);
1391
1392 dits[0] = 0.292;
1393 dits[1] = 0.01;
1394 dits[2] = 24;
1395 gMC->Gsvolu("I129", "BOX ", idtmed[253], dits, 3);
1396
1397 dits[0] = 0.163;
1398 dits[1] = 0.01;
1399 dits[2] = 24;
1400 gMC->Gsvolu("I128", "BOX ", idtmed[253], dits, 3);
1401
1402 dits[0] = 0.04;
1403 dits[1] = 0.06;
1404 dits[2] = 24;
1405 dits[3] = 0;
1406 dits[4] = 161.297;
1407 gMC->Gsvolu("I126", "TUBS", idtmed[253], dits, 5);
1408
1409 dits[0] = 0.2433;
1410 dits[1] = 0.01;
1411 dits[2] = 24;
1412 gMC->Gsvolu("I125", "BOX ", idtmed[253], dits, 3);
1413
1414 dits[0] = 0.06;
1415 dits[1] = 0.08;
1416 dits[2] = 24;
1417 dits[3] = 0;
1418 dits[4] = 42.883;
1419 gMC->Gsvolu("I124", "TUBS", idtmed[253], dits, 5);
1420
1421 dits[0] = 0.793;
fc44ab90 1422 dits[1] = 0.475; //0.685; 0.015
9d2da7b3 1423 dits[2] = 2.5;
fc44ab90 1424 gMC->Gsvolu("I105", "BOX ", idtmed[290], dits, 3);
9d2da7b3 1425
1426 di107[0] = 0.793;
1427 di107[1] = ddet1+dchip1;
1428 di107[2] = 3.536;
1429 gMC->Gsvolu("I107", "BOX ", idtmed[254], di107, 3); // contains det and chip
1430 // layer 1
1431 dits[0] = 0.705;
1432 dits[1] = 0.01;
1433 dits[2] = 2.5;
1434 gMC->Gsvolu("I109", "BOX ", idtmed[275], dits, 3);
1435
1436 di108[0] = 0.705;
1437 di108[1] = dbus;
1438 di108[2] = 14.344;
1439 gMC->Gsvolu("I108", "BOX ", idtmed[275], di108, 3); // bus for both layers
1440
1441 di1d7[0] = 0.7975;
1442 di1d7[1] = ddet2+dchip2;
1443 di1d7[2] = 3.536;
1444 gMC->Gsvolu("I1D7", "BOX ", idtmed[254], di1d7, 3); // contains det and chip
1445 // layer 2
1446 dits[0] = 0.06;
1447 dits[1] = 0.08;
1448 dits[2] = 24;
1449 dits[3] = 0;
1450 dits[4] = 80;
1451 gMC->Gsvolu("I112", "TUBS", idtmed[253], dits, 5);
1452
1453 dits[0] = 0.04;
1454 dits[1] = 0.06;
1455 dits[2] = 24;
1456 dits[3] = 0;
1457 dits[4] = 80;
1458 gMC->Gsvolu("I111", "TUBS", idtmed[253], dits, 5);
1459
1460 dits[0] = 0.15;
1461 dits[1] = 0.0146;
1462 dits[2] = 24;
1463 gMC->Gsvolu("I118", "BOX ", idtmed[273], dits, 3);
1464
1465 dits[0] = 0.1315;
1466 dits[1] = 0.01;
1467 dits[2] = 24;
1468 gMC->Gsvolu("I110", "BOX ", idtmed[253], dits, 3);
1469
1470 dits[0] = 0.025;
1471 dits[1] = 0.035;
1472 dits[2] = 24;
1473 dits[3] = 0;
1474 dits[4] = 180;
1475 gMC->Gsvolu("I114", "TUBS", idtmed[264], dits, 5);
1476
e4f08ac3 1477 if (fluid == 1) {
1478 dits[0] = 0;
1479 dits[1] = 0.025;
1480 dits[2] = 24;
1481 dits[3] = 0;
1482 dits[4] = 180;
1483 gMC->Gsvolu("I115", "TUBS", idtmed[211], dits, 5); // set water as cooling fluid
1484 } else {
1485 dits[0] = 0;
1486 dits[1] = 0.025;
1487 dits[2] = 24;
1488 dits[3] = 0;
1489 dits[4] = 180;
1490 gMC->Gsvolu("I115", "TUBS", idtmed[212], dits, 5); // set freon as cooling fluid
1491 }
1492
9d2da7b3 1493 dits[0] = 0.063;
1494 dits[1] = 0.035;
1495 dits[2] = 24;
1496 gMC->Gsvolu("I116", "BOX ", idtmed[264], dits, 3);
1497
1498 di106[0] = 0.7975;
1499 di106[1] = dchip1;
1500 di106[2] = 0.68;
1501 gMC->Gsvolu("I106", "BOX ", idtmed[201], di106, 3); // chip layer 1
1502
1503 di1d6[0] = 0.7975;
1504 di1d6[1] = dchip2;
1505 di1d6[2] = 0.68;
1506 gMC->Gsvolu("I1D6", "BOX ", idtmed[201], di1d6, 3); // chip layer 2
1507
1508 di101[0] = 0.705;
1509 di101[1] = ddet1;
1510 di101[2] = 3.536;
1511 gMC->Gsvolu("I101", "BOX ", idtmed[250], di101, 3); // contains detector
1512 // layer 1
1513 di1d1[0] = 0.705;
1514 di1d1[1] = ddet2;
1515 di1d1[2] = 3.536;
1516 gMC->Gsvolu("I1D1", "BOX ", idtmed[250], di1d1, 3); // contains detector
1517 // layer 2
e4f08ac3 1518
1519 if (fluid == 1) {
1520 dits[0] = 0.063;
1521 dits[1] = 0.025;
1522 dits[2] = 24;
1523 gMC->Gsvolu("I117", "BOX ", idtmed[211], dits, 3); // set water as cooling fluid
1524 } else {
1525 dits[0] = 0.063;
1526 dits[1] = 0.025;
1527 dits[2] = 24;
1528 gMC->Gsvolu("I117", "BOX ", idtmed[212], dits, 3); // set freon as cooling fluid
1529 }
9d2da7b3 1530
1531 dits1[0] = 0.64;
1532 dits1[1] = ddet1;
1533 dits1[2] = 3.48;
1534 gMC->Gsvolu("ITS1", "BOX ", idtmed[200], dits1, 3); // detector layer 1
1535
1536 dits2[0] = 0.64;
1537 dits2[1] = ddet2;
1538 dits2[2] = 3.48;
1539 gMC->Gsvolu("ITS2", "BOX ", idtmed[200], dits2, 3); // detector layer 2
1540
1541 dits[0] = 3.701;
1542 dits[1] = 7.699;
1543 dits[2] = 4;
1544 dits[3] = 57.1;
1545 dits[4] = 99.9;
1546 gMC->Gsvolu("I650", "TUBS", idtmed[254], dits, 5); // was I150 in old geom.
1547
1548 dits[0] = 0;
1549 dits[1] = 0.5;
1550 dits[2] = 1.5;
1551 gMC->Gsvolu("I676", "TUBE", idtmed[274], dits, 3); // was I176 in old geom.
1552
1553 dits[0] = 0;
1554 dits[1] = 0.18;
1555 dits[2] = 0.8;
1556 gMC->Gsvolu("I673", "TUBE", idtmed[274], dits, 3); // was I173 in old geom.
1557
1558 dits[0] = 0;
1559 dits[1] = 0.18;
1560 dits[2] = 3;
1561 gMC->Gsvolu("I671", "TUBE", idtmed[274], dits, 3); // was I171 in old geom.
1562
1563 dits[0] = 0;
1564 dits[1] = 0.075;
1565 dits[2] = 0.8;
1566 gMC->Gsvolu("I669", "TUBE", idtmed[264], dits, 3); // was I169 in old geom.
1567
1568 dits[0] = 3.5;
1569 dits[1] = 5.6;
1570 dits[2] = 0.55;
1571 dits[3] = 0;
1572 dits[4] = 38;
1573 gMC->Gsvolu("I667", "TUBS", idtmed[263], dits, 5); // was I167 in old geom.
1574
1575 dits[0] = 6.6;
1576 dits[1] = 7.6;
1577 dits[2] = 0.5;
1578 dits[3] = 0;
1579 dits[4] = 9;
1580 gMC->Gsvolu("I666", "TUBS", idtmed[263], dits, 5); // was I166 in old geom.
1581
1582 dits[0] = 0.26;
1583 dits[1] = 0.32;
1584 dits[2] = 0.55;
1585 gMC->Gsvolu("I678", "TUBE", idtmed[263], dits, 3); // was I178 in old geom.
1586
e4f08ac3 1587 if (fluid == 1) {
1588 dits[0] = 0;
1589 dits[1] = 0.3;
1590 dits[2] = 1.5;
1591 gMC->Gsvolu("I677", "TUBE", idtmed[211], dits, 3); //set water as cooling fluid
1592 // was I177 in old geom.
1593 } else {
1594 dits[0] = 0;
1595 dits[1] = 0.3;
1596 dits[2] = 1.5;
1597 gMC->Gsvolu("I677", "TUBE", idtmed[212], dits, 3); //set freon as cooling fluid
1598 // was I177 in old geom.
1599 }
9d2da7b3 1600
1601 dits[0] = 0.07;
1602 dits[1] = 0.125;
1603 dits[2] = 0.3;
1604 gMC->Gsvolu("I675", "TUBE", idtmed[263], dits, 3); // was I175 in old geom.
1605
e4f08ac3 1606 if (fluid == 1) {
1607 dits[0] = 0;
1608 dits[1] = 0.1;
1609 dits[2] = 0.8;
1610 gMC->Gsvolu("I674", "TUBE", idtmed[211], dits, 3); //set water as cooling fluid
1611 // was I174 in old geom.
1612 } else {
1613 dits[0] = 0;
1614 dits[1] = 0.1;
1615 dits[2] = 0.8;
1616 gMC->Gsvolu("I674", "TUBE", idtmed[212], dits, 3); //set freon as cooling fluid
1617 // was I174 in old geom.
1618 }
1619
1620 if (fluid == 1) {
1621 dits[0] = 0;
1622 dits[1] = 0.1;
1623 dits[2] = 3;
1624 gMC->Gsvolu("I672", "TUBE", idtmed[211], dits, 3); //set water as cooling fluid
1625 // was I172 in old geom.
1626 } else {
1627 dits[0] = 0;
1628 dits[1] = 0.1;
1629 dits[2] = 3;
1630 gMC->Gsvolu("I672", "TUBE", idtmed[212], dits, 3); //set freon as cooling fluid
1631 // was I172 in old geom.
1632 }
1633
1634 if (fluid == 1) {
1635 dits[0] = 0;
1636 dits[1] = 0.0746;
1637 dits[2] = 0.8;
1638 gMC->Gsvolu("I670", "TUBE", idtmed[211], dits, 3); //set water as cooling fluid
1639 // was I170 in old geom.
1640 } else {
1641 dits[0] = 0;
1642 dits[1] = 0.0746;
1643 dits[2] = 0.8;
1644 gMC->Gsvolu("I670", "TUBE", idtmed[212], dits, 3); //set freon as cooling fluid
1645 // was I170 in old geom.
1646 }
1647
1648 if (fluid == 1) {
1649 dits[0] = 3.7;
1650 dits[1] = 5.4;
1651 dits[2] = 0.35;
1652 dits[3] = 2;
1653 dits[4] = 36;
1654 gMC->Gsvolu("I668", "TUBS", idtmed[211], dits, 5); //set water as cooling fluid
1655 // was I168 in old geom.
1656 } else {
1657 dits[0] = 3.7;
1658 dits[1] = 5.4;
1659 dits[2] = 0.35;
1660 dits[3] = 2;
1661 dits[4] = 36;
1662 gMC->Gsvolu("I668", "TUBS", idtmed[212], dits, 5); //set freon as cooling fluid
1663 // was I168 in old geom.
1664 }
1665
db915fea 1666
9d2da7b3 1667 }
db915fea 1668
1669 // --- Define SDD volumes ------------------------------------------
1670
9d2da7b3 1671
1672 cos30 = cos(30.*3.14159/180.);
1673 sin30 = sin(30.*3.14159/180.);
1674
1675
db915fea 1676 dits[0] = 0;
1677 dits[1] = 360;
1678 dits[2] = 6;
e4f08ac3 1679 dits[3] = -34.6;
1680 dits[4] = 23.49;
1681 dits[5] = 28;
1682 dits[6] = -23.65;
1683 dits[7] = 23.49;
1684 dits[8] = 28;
1685 dits[9] = -23.65;
1686 dits[10] = 14.59;
1687 dits[11] = 28;
1688 dits[12] = 23.65;
1689 dits[13] = 14.59;
1690 dits[14] = 28;
1691 dits[15] = 23.65;
1692 dits[16] = 23.49;
1693 dits[17] = 28;
1694 dits[18] = 34.6;
1695 dits[19] = 23.49;
1696 dits[20] = 28;
9d2da7b3 1697 gMC->Gsvolu("IT34", "PCON", idtmed[209], dits, 21);
1698
1699 // block of the SDD electronics and related ladder frame
1700 I018dits[0] = 3.2;
1701 I018dits[1] = 2;
1702 I018dits[2] = 3.65;
1703 gMC->Gsvolu("I018", "BOX ", idtmed[209], I018dits, 3);
1704
1705 // block of the SDD end ladder
1706 I024dits[0] = 3.2;
1707 I024dits[1] = 2;
1708 I024dits[2] = 2.725;
1709 gMC->Gsvolu("I024", "BOX ", idtmed[209], I024dits, 3);
1710
1711 // ladder frame of layer 3 - F.T. March,7-2001
1712 I047dits[0] = I018dits[0];
1713 I047dits[1] = I018dits[1];
1714 I047dits[2] = 6*I018dits[2] + 2*I024dits[2];
1715 gMC->Gsvolu("I047", "BOX ", idtmed[209], I047dits, 3);
1716
1717 // ladder frame of layer 4 - F.T. March,7-2001
1718 I048dits[0] = I018dits[0];
1719 I048dits[1] = I018dits[1];
1720 I048dits[2] = 8*I018dits[2] + 2*I024dits[2];
1721 gMC->Gsvolu("I048", "BOX ", idtmed[209], I048dits, 3);
1722
1723
1724 // global SDD volume (sensitive + insensitive)
1725 I302dits[0] = 3.6250;
1726 I302dits[1] = 0.0150;
1727 I302dits[2] = 4.3794;
1728 gMC->Gsvolu("I302", "BOX ", idtmed[278], I302dits, 3);
1729
1730 // Like for I302 - F.T. March,7-2001
1731 I402dits[0] = 3.6250;
1732 I402dits[1] = 0.0150;
1733 I402dits[2] = 4.3794;
1734 gMC->Gsvolu("I402", "BOX ", idtmed[278], I402dits, 3);
1735
1736 // SDD ladder of layer 3 - F.T. March,7-2001
1737 I004dits[0] = I302dits[0]+0.005;
1738 I004dits[1] = 2*I302dits[1]+Y_SDD_sep/2.;
767187dd 1739 I004dits[2] = TMath::Abs(Z_SDD_lay3[0]);
1740 if (I004dits[2] < TMath::Abs(Z_SDD_lay3[5])) {
1741 I004dits[2] = TMath::Abs(Z_SDD_lay3[5]);
9d2da7b3 1742 }
1743 I004dits[2] = I004dits[2] + I302dits[2];
1744 gMC->Gsvolu("I004", "BOX ", idtmed[209], I004dits, 3);
1745
1746 // SDD ladder of layer 4 - F.T. March,7-2001
1747 I005dits[0] = I402dits[0]+0.005;
1748 I005dits[1] = 2*I402dits[1]+Y_SDD_sep/2.;
767187dd 1749 I005dits[2] = TMath::Abs(Z_SDD_lay4[0]);
1750 if (I005dits[2] < TMath::Abs(Z_SDD_lay4[7])) {
1751 I005dits[2] = TMath::Abs(Z_SDD_lay4[7]);
9d2da7b3 1752 }
1753 I005dits[2] = I005dits[2] + I402dits[2];
1754 gMC->Gsvolu("I005", "BOX ", idtmed[209], I005dits, 3);
1755
1756
1757 // -- block of the SDD ladder foot and end ladder
1758
1759 // ladder foot mother volume
1760 I028dits[0] = 3.0000;
1761 I028dits[1] = 0.4000;
1762 I028dits[2] = 0.9000;
1763 gMC->Gsvolu("I028", "BOX ", idtmed[224], I028dits, 3);
1764
1765 // positioning-box #1 at SDD end-ladder - F.T. March,7-2001
1766 I420dits[0] = 0.4500;
1767 I420dits[1] = 0.4000;
1768 I420dits[2] = 0.4500;
1769 gMC->Gsvolu("I420", "BOX ", idtmed[264], I420dits, 3);
1770
1771 // positioning-box #2 at SDD end-ladder - F.T. March,7-2001
e4f08ac3 1772 I421dits[0] = 0.;
1773 I421dits[1] = 0.25;
1774 I421dits[2] = I420dits[1];
1775 gMC->Gsvolu("I421", "TUBE", idtmed[209], I421dits, 3);
9d2da7b3 1776
1777 // reference ruby-sphere at SDD end-ladder - F.T. March,7-2001
1778 I422dits[0] = 0.0000;
1779 I422dits[1] = 0.2000;
1780 I422dits[2] = 0.0000;
1781 I422dits[3] = 180.00;
1782 I422dits[4] = 0.0000;
1783 I422dits[5] = 360.00;
1784 gMC->Gsvolu("I422", "SPHE", idtmed[277], I422dits, 6);
1785
1786 // support for ruby-sphere (I422) - F.T. March,7-2001
1787 I423dits[0] = 0.0000;
1788 I423dits[1] = 0.1000;
e4f08ac3 1789 I423dits[2] = (I420dits[1]-I422dits[1])/2.;
9d2da7b3 1790 gMC->Gsvolu("I423", "TUBE", idtmed[264], I423dits, 3);
1791
1792 // passage for HV microcables - F.T. March,7-2001
1793 I424dits[0] = 1.5000;
1794 I424dits[1] = 0.1500;
1795 I424dits[2] = I421dits[2];
1796 gMC->Gsvolu("I424", "BOX ", idtmed[209], I424dits, 3);
1797
1798 // HV microcables segment at the end ladder - F.T. March,7-2001
1799 I425dits[0] = 1.350000;
1800 I425dits[1] = 0.015250;
1801 I425dits[2] = I024dits[2];
1802 gMC->Gsvolu("I425", "BOX ", idtmed[279], I425dits, 3);
1803
1804 // lower edge of SDD ladder frame at end-ladder - part 1
db915fea 1805 dits[0] = 0.2;
1806 dits[1] = 0.1815;
9d2da7b3 1807 dits[2] = I024dits[2];
db915fea 1808 dits[3] = 0.015;
9d2da7b3 1809 gMC->Gsvolu("I025", "TRD1", idtmed[208], dits, 4);
db915fea 1810
9d2da7b3 1811 // lower edge of SDD ladder frame at end-ladder - part 2
db915fea 1812 dits[0] = 0.183;
1813 dits[1] = 0.165;
9d2da7b3 1814 dits[2] = I024dits[2];
1815 dits[3] = 0.015;
1816 gMC->Gsvolu("I026", "TRD1", idtmed[208], dits, 4);
1817
1818 // new: for the 1st top rod of the structure
1819 // at the end-ladder - F.T. March,7-2001
1820 I029dits[0] = 0.2;
1821 I029dits[1] = 0.1815;
1822 I029dits[2] = 1.0100;
1823 I029dits[3] = 0.015;
1824 gMC->Gsvolu("I029", "TRD1", idtmed[208], I029dits, 4);
1825
1826 // new: for the 2nd top rod of the structure
1827 // at the end-ladder - F.T. March,7-2001
1828 I030dits[0] = 0.1830;
1829 I030dits[1] = 0.1650;
1830 I030dits[2] = 1.0100;
1831 I030dits[3] = 0.0150;
1832 gMC->Gsvolu("I030", "TRD1", idtmed[208], I030dits, 4);
1833
1834 // inox cooling tubes for the end ladder - F.T. March,7-2001
1835 I031dits[0] = 0.093;
1836 I031dits[1] = 0.1;
1837 I031dits[2] = I024dits[2];
1838 gMC->Gsvolu("I031", "TUBE", idtmed[264], I031dits, 3);
1839
e4f08ac3 1840 if (fluid == 1) {
1841 // cooling water for the end ladder - F.T. March,7-2001
1842 I032dits[0] = 0;
1843 I032dits[1] = I031dits[0];
1844 I032dits[2] = I024dits[2];
1845 gMC->Gsvolu("I032", "TUBE", idtmed[211], I032dits, 3);
1846 } else {
1847 // cooling freon for the end ladder - R.B. March,21-2001
1848 I032dits[0] = 0;
1849 I032dits[1] = I031dits[0];
1850 I032dits[2] = I024dits[2];
1851 gMC->Gsvolu("I032", "TUBE", idtmed[212], I032dits, 3);
1852 }
1853
9d2da7b3 1854 // -- block of the SDD ladder frame holding the electronics
1855
1856 // edge of the ladder frame - part 1
1857 dits[0] = 0.2;
1858 dits[1] = 0.182;
1859 dits[2] = 3.65;
db915fea 1860 dits[3] = 0.015;
9d2da7b3 1861 gMC->Gsvolu("I019", "TRD1", idtmed[208], dits, 4);
db915fea 1862
9d2da7b3 1863 // edge of the ladder frame - part 2
1864 dits[0] = 0.183;
1865 dits[1] = 0.165;
1866 dits[2] = 3.65;
1867 dits[3] = 0.015;
1868 gMC->Gsvolu("I020", "TRD1", idtmed[208], dits, 4);
1869
1870 // inclined segments of the ladder frame
db915fea 1871 dits[0] = 2.23;
1872 dits[1] = 2.1;
1873 dits[2] = 0.05;
1874 dits[3] = 0.03;
9d2da7b3 1875 gMC->Gsvolu("I021", "TRD1", idtmed[208], dits, 4);
db915fea 1876
9d2da7b3 1877 // horiz.segments of the ladders, normal to ladder edges
db915fea 1878 dits[0] = 2.1;
1879 dits[1] = 2;
1880 dits[2] = 0.06;
1881 dits[3] = 0.04;
9d2da7b3 1882 gMC->Gsvolu("I022", "TRD1", idtmed[208], dits, 4);
db915fea 1883
9d2da7b3 1884 // horiz.segments of the ladders, at 45 deg. to ladder edges
1885 dits[0] = 2.615;
1886 dits[1] = 2.465;
1887 dits[2] = 0.06;
1888 dits[3] = 0.04;
1889 gMC->Gsvolu("I023", "TRD1", idtmed[208], dits, 4);
db915fea 1890
9d2da7b3 1891 // supports of the ceramic pins holding the detectors
db915fea 1892 dits[0] = 0.3;
1893 dits[1] = 0.05;
1894 dits[2] = 0.15;
9d2da7b3 1895 gMC->Gsvolu("I033", "BOX ", idtmed[208], dits, 3);
db915fea 1896
9d2da7b3 1897 // ceramic pins holding the detectors
db915fea 1898 dits[0] = 0;
1899 dits[1] = 0.05;
1900 dits[2] = 0.225;
9d2da7b3 1901 gMC->Gsvolu("I034", "TUBE", idtmed[277], dits, 3);
db915fea 1902
9d2da7b3 1903 // holders of cooling tubes
1904 I035dits[0] = 0.1;
1905 I035dits[1] = 0.15;
1906 I035dits[2] = 0.2;
1907 gMC->Gsvolu("I035", "TUBE", idtmed[208], I035dits, 3);
db915fea 1908
9d2da7b3 1909 // top holders of microcables
1910 dits[0] = 0.2;
db915fea 1911 dits[1] = 0.01;
9d2da7b3 1912 dits[2] = 0.05;
1913 gMC->Gsvolu("I036", "BOX ", idtmed[208], dits, 3);
db915fea 1914
9d2da7b3 1915 // inox cooling tubes - F.T. March,7-2001
1916 I037dits[0] = 0.093;
1917 I037dits[1] = 0.1;
1918 I037dits[2] = I018dits[2];
1919 gMC->Gsvolu("I037", "TUBE", idtmed[264], I037dits, 3);
db915fea 1920
e4f08ac3 1921 if (fluid == 1) {
1922 // cooling water - F.T. March,7-2001
1923 I038dits[0] = 0;
1924 I038dits[1] = I037dits[0];
1925 I038dits[2] = I018dits[2];
1926 gMC->Gsvolu("I038", "TUBE", idtmed[211], I038dits, 3);
1927 } else {
1928 // cooling freon - R.B. March,21-2001
1929 I038dits[0] = 0;
1930 I038dits[1] = I037dits[0];
1931 I038dits[2] = I018dits[2];
1932 gMC->Gsvolu("I038", "TUBE", idtmed[212], I038dits, 3);
1933 }
9d2da7b3 1934 // -- block of the SDD electronics (heat bridge, chips, hybrid, anode microcable)
db915fea 1935
9d2da7b3 1936 // SDD heat bridge - F.T. March,7-2001
1937 I039dits[0] = 1.1000;
1938 I039dits[1] = 0.0087;
1939 I039dits[2] = 3.2500;
1940 gMC->Gsvolu("I039", "BOX ", idtmed[268], I039dits, 3);
db915fea 1941
9d2da7b3 1942 // SDD clip part 1
1943 dits[0] = 0.25;
1944 dits[1] = 0.01;
1945 dits[2] = I039dits[2];
1946 gMC->Gsvolu("I040", "BOX ", idtmed[268], dits, 3);
1947
1948 // SDD clip part 2
1949 I041dits[0] = 0.1;
1950 I041dits[1] = 0.12;
1951 I041dits[2] = I039dits[2];
1952 I041dits[3] = 90;
1953 I041dits[4] = 320;
1954 gMC->Gsvolu("I041", "TUBS", idtmed[268], I041dits, 5);
1955
1956
1957 // SDD PASCAL - F.T. March,7-2001
1958 I042dits[0] = 0.5000;
1959 I042dits[1] = 0.0175;
1960 I042dits[2] = 0.5000;
1961 gMC->Gsvolu("I042", "BOX ", idtmed[206], I042dits, 3);
1962
1963 // SDD AMBRA - F.T. March,7-2001
1964 I043dits[0] = 0.3500;
1965 I043dits[1] = 0.0175;
1966 I043dits[2] = 0.5000;
1967 gMC->Gsvolu("I043", "BOX ", idtmed[206], I043dits, 3);
1968
1969 // SDD capacitors - F.T. March,7-2001
1970 I051dits[0] = 0.1400;
1971 I051dits[1] = 0.0350;
1972 I051dits[2] = 0.0625;
1973 gMC->Gsvolu("I051", "BOX ", idtmed[276], I051dits, 3);
1974
1975 // SDD hybrid circuit - F.T. March,7-2001
1976 I052dits[0] = 1.725000;
1977 I052dits[1] = 0.003743;
1978 I052dits[2] = I039dits[2];
e4f08ac3 1979 gMC->Gsvolu("I052", "BOX ", idtmed[281], I052dits, 3);
9d2da7b3 1980
1981 // SDD anode microcable : changed - F.T. March,7-2001
1982 I044dits[0] = I018dits[2];
1983 I044dits[1] = I039dits[2];
1984 I044dits[2] = 0.00084;
1985 I044dits[3] = (15.189149/(I044dits[0]+I044dits[1]))/2;
1986 gMC->Gsvolu("I044", "TRD1", idtmed[282], I044dits, 4);
1987 volI044 = ((2*I044dits[0] + 2*I044dits[1]) * 2*I044dits[2])/2 * 2*I044dits[3];
1988
1989 // SDD electronics box - F.T. March,7-2001
1990 I050dits[1] = I039dits[1]+I052dits[1]+I051dits[1]+I044dits[2];
1991 I050dits[0] = I018dits[1]/cos(30.*3.14159/180.)-I050dits[1]*sin(30.*3.14159/180.);
1992 I050dits[2] = I018dits[2];
1993 gMC->Gsvolu("I050", "BOX ", idtmed[209], I050dits, 3);
1994
1995 // SDD sensitive volume
1996 dits[0] = 3.50850;
1997 dits[1] = 0.01499; // not 0.015 because it is included into I302 which is 0.015
1998 dits[2] = 3.76320;
db915fea 1999 gMC->Gsvolu("ITS3", "BOX ", idtmed[200], dits, 3);
2000
9d2da7b3 2001 // Like for ITS3 - F.T. March,7-2001
2002 dits[0] = 3.50850;
2003 dits[1] = 0.01499; // not 0.015 because it is included into I402 which is 0.015
2004 dits[2] = 3.76320;
db915fea 2005 gMC->Gsvolu("ITS4", "BOX ", idtmed[200], dits, 3);
2006
2007
2008 // --- Define SSD volumes ------------------------------------------
2009
2010
2011 dits[0] = 0;
2012 dits[1] = 360;
2013 dits[2] = 6;
e4f08ac3 2014 dits[3] = -57.45;
2015 dits[4] = 43.6;
db915fea 2016 dits[5] = 48;
e4f08ac3 2017 dits[6] = -49.15;
2018 dits[7] = 43.6;
db915fea 2019 dits[8] = 48;
e4f08ac3 2020 dits[9] = -49.15;
2021 dits[10] = 36.9;
db915fea 2022 dits[11] = 48;
e4f08ac3 2023 dits[12] = 50.55;
2024 dits[13] = 36.9;
db915fea 2025 dits[14] = 48;
e4f08ac3 2026 dits[15] = 50.55;
2027 dits[16] = 43.6;
db915fea 2028 dits[17] = 48;
e4f08ac3 2029 dits[18] = 57.45;
2030 dits[19] = 43.6;
db915fea 2031 dits[20] = 48;
9d2da7b3 2032 gMC->Gsvolu("IT56", "PCON", idtmed[220], dits, 21);
db915fea 2033
2034 dits[0] = 3.4;
2035 dits[1] = 1.955;
9d2da7b3 2036 dits[2] = 56.5;
2037 gMC->Gsvolu("I570", "BOX ", idtmed[204], dits, 3);
db915fea 2038
2039 dits[0] = 3.75;
2040 dits[1] = 0.045;
2041 dits[2] = 50.975;
9d2da7b3 2042 gMC->Gsvolu("I569", "BOX ", idtmed[204], dits, 3);
db915fea 2043
2044 dits[0] = 3.4;
2045 dits[1] = 1.955;
9d2da7b3 2046 dits[2] = 47;
2047 gMC->Gsvolu("I571", "BOX ", idtmed[204], dits, 3);
db915fea 2048
2049 dits[0] = 3.75;
2050 dits[1] = 0.045;
9d2da7b3 2051 dits[2] = 43.3;
2052 gMC->Gsvolu("I565", "BOX ", idtmed[204], dits, 3);
db915fea 2053
2054 dits[0] = 3.4;
2055 dits[1] = 1.955;
2056 dits[2] = 3.15;
9d2da7b3 2057 gMC->Gsvolu("I553", "BOX ", idtmed[204], dits, 3);
db915fea 2058
2059 dits[0] = 3.405;
2060 dits[1] = 1.955;
2061 dits[2] = 1.955;
9d2da7b3 2062 gMC->Gsvolu("I523", "BOX ", idtmed[204], dits, 3);
db915fea 2063
2064 dits[0] = 3.75;
2065 dits[1] = 0.015;
2066 dits[2] = 2.1;
9d2da7b3 2067 gMC->Gsvolu("I566", "BOX ", idtmed[206], dits, 3);
db915fea 2068
2069 dits[0] = 3.4;
2070 dits[1] = 1.955;
2071 dits[2] = 3.15;
9d2da7b3 2072 gMC->Gsvolu("I544", "BOX ", idtmed[204], dits, 3);
db915fea 2073
2074 dits[0] = 3.41;
2075 dits[1] = 1.955;
2076 dits[2] = 1.955;
9d2da7b3 2077 gMC->Gsvolu("I516", "BOX ", idtmed[204], dits, 3);
db915fea 2078
2079 dits[0] = 3.75;
2080 dits[1] = 0.015;
2081 dits[2] = 2.1;
9d2da7b3 2082 gMC->Gsvolu("I562", "BOX ", idtmed[206], dits, 3);
db915fea 2083
e4f08ac3 2084 if (fluid == 1) {
2085 dits[0] = 0;
2086 dits[1] = 0.07;
2087 dits[2] = 3.15;
2088 gMC->Gsvolu("I559", "TUBE", idtmed[211], dits, 3); // set water as cooling fluid
2089 } else {
2090 dits[0] = 0;
2091 dits[1] = 0.07;
2092 dits[2] = 3.15;
2093 gMC->Gsvolu("I559", "TUBE", idtmed[212], dits, 3); // set freon as cooling fluid
2094 }
db915fea 2095
2096 dits[0] = 0.07;
2097 dits[1] = 0.1;
2098 dits[2] = 3.15;
9d2da7b3 2099 gMC->Gsvolu("I560", "TUBE", idtmed[210], dits, 3);
db915fea 2100
2101 dits[0] = 0.225;
2102 dits[1] = 0.195;
2103 dits[2] = 3.15;
2104 dits[3] = 0.025;
9d2da7b3 2105 gMC->Gsvolu("I558", "TRD1", idtmed[203], dits, 4);
db915fea 2106
2107 dits[0] = 0.25;
2108 dits[1] = 0.22;
2109 dits[2] = 3.15;
2110 dits[3] = 0.025;
9d2da7b3 2111 gMC->Gsvolu("I557", "TRD1", idtmed[203], dits, 4);
db915fea 2112
2113 dits[0] = 2.17;
2114 dits[1] = 0.035;
2115 dits[2] = 0.05;
9d2da7b3 2116 gMC->Gsvolu("I556", "BOX ", idtmed[203], dits, 3);
db915fea 2117
2118 dits[0] = 2 ;
2119 dits[1] = 0.035;
2120 dits[2] = 0.05;
9d2da7b3 2121 gMC->Gsvolu("I554", "BOX ", idtmed[203], dits, 3);
db915fea 2122
2123 dits[0] = 2.675;
2124 dits[1] = 0.035;
2125 dits[2] = 0.05;
9d2da7b3 2126 gMC->Gsvolu("I555", "BOX ", idtmed[203], dits, 3);
db915fea 2127
2128 dits[0] = 0.3;
2129 dits[1] = 0.15;
2130 dits[2] = 0.15;
9d2da7b3 2131 gMC->Gsvolu("I561", "BOX ", idtmed[203], dits, 3);
db915fea 2132
2133 dits[0] = 0.025;
2134 dits[1] = 0.025;
2135 dits[2] = 0.05;
9d2da7b3 2136 gMC->Gsvolu("I519", "BOX ", idtmed[214], dits, 3);
db915fea 2137
2138 dits[0] = 0.304;
2139 dits[1] = 0.0275;
2140 dits[2] = 0.432;
9d2da7b3 2141 gMC->Gsvolu("I521", "BOX ", idtmed[206], dits, 3);
db915fea 2142
2143 dits[0] = 0.16;
2144 dits[1] = 0.08;
2145 dits[2] = 0.08;
9d2da7b3 2146 gMC->Gsvolu("I520", "BOX ", idtmed[214], dits, 3);
db915fea 2147
2148 dits[0] = 3.4;
2149 dits[1] = 0.015;
2150 dits[2] = 0.525;
9d2da7b3 2151 gMC->Gsvolu("I518", "BOX ", idtmed[203], dits, 3);
db915fea 2152
2153 dits[0] = 0.15;
2154 dits[1] = 0.105;
2155 dits[2] = 0.29;
2156 dits[3] = 0.08;
9d2da7b3 2157 gMC->Gsvolu("I522", "TRD1", idtmed[203], dits, 4);
db915fea 2158
2159 dits[0] = 0.07;
2160 dits[1] = 0.1;
2161 dits[2] = 1.955;
9d2da7b3 2162 gMC->Gsvolu("I542", "TUBE", idtmed[210], dits, 3);
db915fea 2163
e4f08ac3 2164 if (fluid == 1) {
2165 dits[0] = 0;
2166 dits[1] = 0.07;
2167 dits[2] = 1.955;
2168 gMC->Gsvolu("I541", "TUBE", idtmed[211], dits, 3); // set water as cooling fluid
2169 } else {
2170 dits[0] = 0;
2171 dits[1] = 0.07;
2172 dits[2] = 1.955;
2173 gMC->Gsvolu("I541", "TUBE", idtmed[212], dits, 3); // set freon as cooling fluid
2174 }
db915fea 2175
2176 dits[0] = 0.3;
2177 dits[1] = 0.15;
2178 dits[2] = 0.15;
9d2da7b3 2179 gMC->Gsvolu("I543", "BOX ", idtmed[203], dits, 3);
db915fea 2180
2181 dits[0] = 0.25;
2182 dits[1] = 0.22;
2183 dits[2] = 1.955;
2184 dits[3] = 0.025;
9d2da7b3 2185 gMC->Gsvolu("I537", "TRD1", idtmed[203], dits, 4);
db915fea 2186
2187 dits[0] = 0.225;
2188 dits[1] = 0.195;
2189 dits[2] = 1.955;
2190 dits[4] = 0.025;
9d2da7b3 2191 gMC->Gsvolu("I538", "TRD1", idtmed[203], dits, 4);
db915fea 2192
2193 dits[0] = 2.17;
2194 dits[1] = 0.035;
2195 dits[2] = 0.05;
9d2da7b3 2196 gMC->Gsvolu("I536", "BOX ", idtmed[203], dits, 3);
db915fea 2197
2198 dits[0] = 2.675;
2199 dits[1] = 0.035;
2200 dits[2] = 0.05;
9d2da7b3 2201 gMC->Gsvolu("I535", "BOX ", idtmed[203], dits, 3);
db915fea 2202
2203 dits[0] = 2;
2204 dits[1] = 0.035;
2205 dits[2] = 0.05;
9d2da7b3 2206 gMC->Gsvolu("I534", "BOX ", idtmed[203], dits, 3);
db915fea 2207
2208 dits[0] = 0;
2209 dits[1] = 0.05;
2210 dits[2] = 0.17;
9d2da7b3 2211 gMC->Gsvolu("I540", "TUBE", idtmed[203], dits, 3);
db915fea 2212
2213 dits[0] = 0;
2214 dits[1] = 0.05;
2215 dits[2] = 0.205;
9d2da7b3 2216 gMC->Gsvolu("I539", "TUBE", idtmed[203], dits, 3);
db915fea 2217
2218 dits[0] = 3.65;
2219 dits[1] = 0.015;
2220 dits[2] = 2;
2221 gMC->Gsvolu("ITS6", "BOX ", idtmed[200], dits, 3);
2222
e4f08ac3 2223 if (fluid == 1) {
2224 dits[0] = 0;
2225 dits[1] = 0.07;
2226 dits[2] = 3.15;
2227 gMC->Gsvolu("I550", "TUBE", idtmed[211], dits, 3); // set water as cooling fluid
2228 } else {
2229 dits[0] = 0;
2230 dits[1] = 0.07;
2231 dits[2] = 3.15;
2232 gMC->Gsvolu("I550", "TUBE", idtmed[212], dits, 3); // set freon as cooling fluid
2233 }
db915fea 2234
2235 dits[0] = 0.07;
2236 dits[1] = 0.1;
2237 dits[2] = 3.15;
9d2da7b3 2238 gMC->Gsvolu("I551", "TUBE", idtmed[210], dits, 3);
db915fea 2239
2240 dits[0] = 0.225;
2241 dits[1] = 0.195;
2242 dits[2] = 3.15;
2243 dits[3] = 0.025;
9d2da7b3 2244 gMC->Gsvolu("I549", "TRD1", idtmed[203], dits, 4);
db915fea 2245
2246 dits[0] = 0.25;
2247 dits[1] = 0.22;
2248 dits[2] = 3.15;
2249 dits[3] = 0.025;
9d2da7b3 2250 gMC->Gsvolu("I548", "TRD1", idtmed[203], dits, 4);
db915fea 2251
2252 dits[0] = 2.17;
2253 dits[1] = 0.035;
2254 dits[2] = 0.05;
9d2da7b3 2255 gMC->Gsvolu("I547", "BOX ", idtmed[203], dits, 3);
db915fea 2256
2257 dits[0] = 2;
2258 dits[1] = 0.035;
2259 dits[2] = 0.05;
9d2da7b3 2260 gMC->Gsvolu("I545", "BOX ", idtmed[203], dits, 3);
db915fea 2261
2262 dits[0] = 2.675;
2263 dits[1] = 0.035;
2264 dits[2] = 0.05;
9d2da7b3 2265 gMC->Gsvolu("I546", "BOX ", idtmed[203], dits, 3);
db915fea 2266
2267 dits[0] = 0.3;
2268 dits[1] = 0.15;
2269 dits[2] = 0.15;
9d2da7b3 2270 gMC->Gsvolu("I552", "BOX ", idtmed[203], dits, 3);
db915fea 2271
2272 dits[0] = 0.304;
2273 dits[1] = 0.0275;
2274 dits[2] = 0.4322;
9d2da7b3 2275 gMC->Gsvolu("I515", "BOX ", idtmed[206], dits, 3);
db915fea 2276
2277 dits[0] = 0.025;
2278 dits[1] = 0.025;
2279 dits[2] = 0.05;
9d2da7b3 2280 gMC->Gsvolu("I513", "BOX ", idtmed[214], dits, 3);
db915fea 2281
2282 dits[0] = 0.16;
2283 dits[1] = 0.08;
2284 dits[2] = 0.08;
9d2da7b3 2285 gMC->Gsvolu("I514", "BOX ", idtmed[214], dits, 3);
db915fea 2286
2287 dits[0] = 3.4;
2288 dits[1] = 0.015;
2289 dits[2] = 0.525;
9d2da7b3 2290 gMC->Gsvolu("I512", "BOX ", idtmed[203], dits, 3);
db915fea 2291
2292 dits[0] = 0.225;
2293 dits[1] = 0.195;
2294 dits[2] = 1.955;
2295 dits[3] = 0.025;
9d2da7b3 2296 gMC->Gsvolu("I528", "TRD1", idtmed[203], dits, 4);
db915fea 2297
2298 dits[0] = 0.25;
2299 dits[1] = 0.22;
2300 dits[2] = 1.955;
2301 dits[3] = 0.025;
9d2da7b3 2302 gMC->Gsvolu("I527", "TRD1", idtmed[203], dits, 4);
db915fea 2303
2304 dits[0] = 2.17;
2305 dits[1] = 0.035;
2306 dits[2] = 0.05;
9d2da7b3 2307 gMC->Gsvolu("I526", "BOX ", idtmed[203], dits, 3);
db915fea 2308
2309 dits[0] = 2.675;
2310 dits[1] = 0.035;
2311 dits[2] = 0.05;
9d2da7b3 2312 gMC->Gsvolu("I525", "BOX ", idtmed[203], dits, 3);
db915fea 2313
2314 dits[0] = 2;
2315 dits[1] = 0.035;
2316 dits[2] = 0.05;
9d2da7b3 2317 gMC->Gsvolu("I524", "BOX ", idtmed[203], dits, 3);
db915fea 2318
2319 dits[0] = 0;
2320 dits[1] = 0.05;
2321 dits[2] = 0.205;
9d2da7b3 2322 gMC->Gsvolu("I529", "TUBE", idtmed[203], dits, 3);
db915fea 2323
2324 dits[0] = 0;
2325 dits[1] = 0.05;
2326 dits[2] = 0.17;
9d2da7b3 2327 gMC->Gsvolu("I530", "TUBE", idtmed[203], dits, 3);
db915fea 2328
2329 dits[0] = 0.15;
2330 dits[1] = 0.105;
2331 dits[2] = 0.29;
2332 dits[3] = 0.08;
9d2da7b3 2333 gMC->Gsvolu("I517", "TRD1", idtmed[203], dits, 4);
db915fea 2334
e4f08ac3 2335 if (fluid == 1) {
2336 dits[0] = 0;
2337 dits[1] = 0.07;
2338 dits[2] = 1.955;
2339 gMC->Gsvolu("I531", "TUBE", idtmed[211], dits, 3); // set water as cooling fluid
2340 } else {
2341 dits[0] = 0;
2342 dits[1] = 0.07;
2343 dits[2] = 1.955;
2344 gMC->Gsvolu("I531", "TUBE", idtmed[212], dits, 3); // set freon as cooling fluid
2345 }
db915fea 2346
2347 dits[0] = 0.07;
2348 dits[1] = 0.1;
2349 dits[2] = 1.955;
9d2da7b3 2350 gMC->Gsvolu("I532", "TUBE", idtmed[210], dits, 3);
db915fea 2351
2352 dits[0] = 0.3;
2353 dits[1] = 0.15;
2354 dits[2] = 0.15;
9d2da7b3 2355 gMC->Gsvolu("I533", "BOX ", idtmed[203], dits, 3);
db915fea 2356
2357 dits[0] = 3.65;
2358 dits[1] = 0.015;
2359 dits[2] = 2;
2360 gMC->Gsvolu("ITS5", "BOX ", idtmed[200], dits, 3);
2361
2362
e4f08ac3 2363
2364 // --- Define volumes of shield of SPD ----------------
9d2da7b3 2365
2366
2367 dits[0] = 8.37;
2368 dits[1] = 9.93;
2369 dits[2] = 25;
2370 gMC->Gsvolu("IC01", "TUBE", idtmed[289], dits, 3);
2371
2372 dits[0] = 8.3;
e4f08ac3 2373 dits[1] = 9.995;
9d2da7b3 2374 dits[2] = 17.5/2.;
2375 gMC->Gsvolu("IC02", "TUBE", idtmed[289], dits, 3);
2376
2377
e4f08ac3 2378 // --- Define volume of first cylinder between SPD and SDD --------------
9d2da7b3 2379
e4f08ac3 2380 dits[0] = (21.-0.128)/2.;
2381 dits[1] = 21./2.;
2382 dits[2] = 39.4;
2383 gMC->Gsvolu("ICY1", "TUBE", idtmed[208], dits, 3);
9d2da7b3 2384
e4f08ac3 2385 // --- Define volume of second cylinder between SDD and SSD --------------
9d2da7b3 2386
e4f08ac3 2387 dits[0] = (59.5-0.128)/2.;
2388 dits[1] = 59.5/2.;
2389 dits[2] = 56.2; // was 57
2390 gMC->Gsvolu("ICY2", "TUBE", idtmed[208], dits, 3);
9d2da7b3 2391
2392 // --- Define volumes of SDD cone ----------------------------------
2393
2394 dits[0] = 0;
2395 dits[1] = 360;
2396 dits[2] = 12;
2397 dits[3] = -59.7;
2398 dits[4] = 27;
e4f08ac3 2399 dits[5] = 28.6;
9d2da7b3 2400 dits[6] = -42.7;
2401 dits[7] = 10;
e4f08ac3 2402 dits[8] = 28.6;
9d2da7b3 2403 dits[9] = -34.65;
2404 dits[10] = 10;
e4f08ac3 2405 dits[11] = 28.6;
9d2da7b3 2406 dits[12] = -34.65;
2407 dits[13] = 10;
2408 dits[14] = 23.495;
2409 dits[15] = -23.7;
2410 dits[16] = 10;
2411 dits[17] = 23.495;
2412 dits[18] = -23.7;
2413 dits[19] = 10;
2414 dits[20] = 14.595;
2415 dits[21] = 23.7;
2416 dits[22] = 10;
2417 dits[23] = 14.595;
2418 dits[24] = 23.7;
2419 dits[25] = 10;
2420 dits[26] = 23.495;
2421 dits[27] = 34.65;
2422 dits[28] = 10;
2423 dits[29] = 23.495;
2424 dits[30] = 34.65;
2425 dits[31] = 10;
e4f08ac3 2426 dits[32] = 28.6;
2427 dits[33] = 42.7;
9d2da7b3 2428 dits[34] = 10;
e4f08ac3 2429 dits[35] = 28.6;
9d2da7b3 2430 dits[36] = 59.7;
2431 dits[37] = 27.2637;
e4f08ac3 2432 dits[38] = 28.6;
9d2da7b3 2433 gMC->Gsvolu("IS02", "PCON", idtmed[204], dits, 39);
2434
2435 dits[0] = 0;
2436 dits[1] = 360;
2437 dits[2] = 6;
e4f08ac3 2438 dits[3] = 38.65;
2439 dits[4] = 10.75;
2440 dits[5] = 12.25;
2441 dits[6] = 40.15;
2442 dits[7] = 10.75;
2443 dits[8] = 13.96;
2444 dits[9] = 40.15;
2445 dits[10] = 12.46;
9d2da7b3 2446 dits[11] = 13.96;
e4f08ac3 2447 dits[12] = 55.75;
2448 dits[13] = 27;
2449 dits[14] = 28.5;
2450 dits[15] = 55.75;
2451 dits[16] = 27;
2452 dits[17] = 28.5;
2453 dits[18] = 57.25;
2454 dits[19] = 27;
2455 dits[20] = 28.5;
9d2da7b3 2456 gMC->Gsvolu("I093", "PCON", idtmed[272], dits, 21); // SDD cone
2457
2458 dits[0] = 0;
2459 dits[1] = 50;
2460 dits[2] = 3;
e4f08ac3 2461 dits[3] = 39;
9d2da7b3 2462 dits[4] = 14;
2463 dits[5] = 18.75;
2464 dits[6] = 46.7-3;
2465 dits[7] = 14;
2466 dits[8] = 18.75;
2467 dits[9] = 51.45-3;
2468 dits[10] = 18.75;
2469 dits[11] = 18.75;
e4f08ac3 2470 gMC->Gsvolu("I099", "PCON", idtmed[204], dits, 12); // SDD 3 cone hole
9d2da7b3 2471
2472 dits[0] = 0;
2473 dits[1] = 25;
2474 dits[2] = 3;
e4f08ac3 2475 dits[3] = 49;
9d2da7b3 2476 dits[4] = 23.4;
2477 dits[5] = 26.4;
2478 dits[6] = 56.1-3;
2479 dits[7] = 23.4;
2480 dits[8] = 26.4;
2481 dits[9] = 59.1-3;
2482 dits[10] = 26.4;
2483 dits[11] = 26.4;
e4f08ac3 2484 gMC->Gsvolu("I200", "PCON", idtmed[204], dits, 12); // SDD 4 cone hole
9d2da7b3 2485
e4f08ac3 2486 dits[0] = 10.0;
2487 dits[1] = 10.5;
9d2da7b3 2488 dits[2] = 0.25;
e4f08ac3 2489 gMC->Gsvolu("I090", "TUBE", idtmed[224], dits, 3); // SDD cylinder flange
9d2da7b3 2490
2491 dits[0] = 21.95;
e4f08ac3 2492 dits[1] = 22.95;
9d2da7b3 2493 dits[2] = 1;
e4f08ac3 2494 gMC->Gsvolu("I098", "TUBE", idtmed[283], dits, 3); // ladder support on layer 4
9d2da7b3 2495
e4f08ac3 2496 dits[0] = 13.1;
2497 dits[1] = 14.1;
9d2da7b3 2498 dits[2] = 1;
e4f08ac3 2499 gMC->Gsvolu("I097", "TUBE", idtmed[283], dits, 3); // ladder support on layer 3
9d2da7b3 2500
2501 dits[0] = 1;
2502 dits[1] = 1;
2503 dits[2] = 7.74;
2504 gMC->Gsvolu("I202", "BOX ", idtmed[272], dits, 3);
2505
2506 dits[0] = 1;
2507 dits[1] = 1;
2508 dits[2] = 9.14;
2509 gMC->Gsvolu("I203", "BOX ", idtmed[272], dits, 3);
2510
2511 dits[0] = 21.95;
2512 dits[1] = 22.95;
2513 dits[2] = 1;
2514 gMC->Gsvolu("I095", "TUBE", idtmed[224], dits, 3);
2515
2516 dits[0] = 3;
2517 dits[1] = 2.7;
2518 dits[2] = 1;
2519 dits[3] = 0.63;
2520 gMC->Gsvolu("I096", "TRD1", idtmed[264], dits, 4);
2521
2522 dits[0] = 13.1;
2523 dits[1] = 14.1;
2524 dits[2] = 1;
2525 gMC->Gsvolu("I094", "TUBE", idtmed[224], dits, 3);
2526
2527
2528 // --- Define volumes of SSD cone ----------------------------------
2529
2530
2531 dits[0] = 0;
2532 dits[1] = 360;
2533 dits[2] = 12;
e4f08ac3 2534 dits[3] = -zmax;
2535 dits[4] = 46;
2536 dits[5] = 49.25;
9d2da7b3 2537 dits[6] = -61.2;
e4f08ac3 2538 dits[7] = 28.7;
2539 dits[8] = 49.25;
9d2da7b3 2540 dits[9] = -57.5;
e4f08ac3 2541 dits[10] = 28.7;
2542 dits[11] = 49.25;
9d2da7b3 2543 dits[12] = -57.5;
e4f08ac3 2544 dits[13] = 28.7;
9d2da7b3 2545 dits[14] = 43.5;
e4f08ac3 2546 dits[15] = -49.2;
2547 dits[16] = 28.7;
9d2da7b3 2548 dits[17] = 43.5;
e4f08ac3 2549 dits[18] = -49.2;
2550 dits[19] = 28.7;
2551 dits[20] = 36.85;
2552 dits[21] = 50.6;
2553 dits[22] = 28.7;
2554 dits[23] = 36.85;
2555 dits[24] = 50.6;
2556 dits[25] = 28.7;
9d2da7b3 2557 dits[26] = 43.5;
e4f08ac3 2558 dits[27] = 57.5;
2559 dits[28] = 28.7;
9d2da7b3 2560 dits[29] = 43.5;
e4f08ac3 2561 dits[30] = 57.5;
2562 dits[31] = 28.7;
2563 dits[32] = 49.25;
2564 dits[33] = 61.2;
2565 dits[34] = 28.7;
2566 dits[35] = 49.25;
2567 dits[36] = zmax;
2568 dits[37] = 46;
2569 dits[38] = 49.25;
2570 gMC->Gsvolu("IS01", "PCON", idtmed[204], dits, 39); // SSD cone mother volume
9d2da7b3 2571
2572 dits[0] = 0;
2573 dits[1] = 360;
2574 dits[2] = 6;
e4f08ac3 2575 dits[3] = -zmax;
2576 dits[4] = 47.75;
2577 dits[5] = 49.25;
2578 dits[6] = -zmax+2.;
2579 dits[7] = 47.75;
2580 dits[8] = 49.25;
2581 dits[9] = -71.2819;
2582 dits[10] = 46.75;
2583 dits[11] = 49.0319;
2584 dits[12] = -57.25; // was 58.5
2585 dits[13] = 32.9681;
2586 dits[14] = 34.75;
2587 dits[15] = -57.25; // was 58.5
2588 dits[16] = 30;
2589 dits[17] = 34.75;
2590 dits[18] = -55.75; // was 57
2591 dits[19] = 30;
2592 dits[20] = 32.25; // was 31.5
9d2da7b3 2593 gMC->Gsvolu("I212", "PCON", idtmed[272], dits, 21); // SSD cone
2594
e4f08ac3 2595 dits[0] = 28.75;
2596 dits[1] = 29.75;
9d2da7b3 2597 dits[2] = 0.5;
e4f08ac3 2598 gMC->Gsvolu("I211", "TUBE", idtmed[224], dits, 3); // SSD cylinder flange
9d2da7b3 2599
e4f08ac3 2600 dits[0] = 35.8;
2601 dits[1] = 36.8;
9d2da7b3 2602 dits[2] = 1;
e4f08ac3 2603 gMC->Gsvolu("I217", "TUBE", idtmed[283], dits, 3); // ladder support on layer 5
9d2da7b3 2604
e4f08ac3 2605 dits[0] = 41.4;
2606 dits[1] = 42.4;
9d2da7b3 2607 dits[2] = 1;
e4f08ac3 2608 gMC->Gsvolu("I219", "TUBE", idtmed[283], dits, 3); // ladder support on layer 6
2609
2610 dits[0] = 42.05+5.;
2611 dits[1] = 42.55+5.;
9d2da7b3 2612 dits[2] = 1.25;
e4f08ac3 2613 gMC->Gsvolu("I214", "TUBE", idtmed[224], dits, 3); // layer 6 electronic support
2614 // this will change after PPR
2615 dits[0] = 37.05+5.;
2616 dits[1] = 37.55+5.;
9d2da7b3 2617 dits[2] = 1.25;
e4f08ac3 2618 gMC->Gsvolu("I213", "TUBE", idtmed[224], dits, 3); // layer 5 electronic support
2619 // this will change after PPR
9d2da7b3 2620 dits[0] = 0;
2621 dits[1] = 25;
2622 dits[2] = 5;
e4f08ac3 2623 dits[3] = -zmax+3;
2624 dits[4] = 45.5;
2625 dits[5] = 45.5;
9d2da7b3 2626 dits[6] = -69.7+3;
2627 dits[7] = 37;
e4f08ac3 2628 dits[8] = 45.5;
9d2da7b3 2629 dits[9] = -68.5+3;
2630 dits[10] = 37;
e4f08ac3 2631 dits[11] = 45.5;
9d2da7b3 2632 dits[12] = -68.5+4.8;
2633 dits[13] = 37;
e4f08ac3 2634 dits[14] = 45.5;
9d2da7b3 2635 dits[15] = -63.5+4.8;
2636 dits[16] = 37;
e4f08ac3 2637 dits[17] = 45.5;
2638 gMC->Gsvolu("I215", "PCON", idtmed[204], dits, 18); // SSD cone hole
9d2da7b3 2639
2640 dits[0] = 0;
2641 dits[1] = 3.2;
2642 dits[2] = 9;
e4f08ac3 2643 dits[3] = -14;
9d2da7b3 2644 dits[4] = 30.5;
2645 dits[5] = 33.5;
e4f08ac3 2646 dits[6] = -9.85;
9d2da7b3 2647 dits[7] = 30.5;
2648 dits[8] = 33.5;
e4f08ac3 2649 dits[9] = -9.85;
9d2da7b3 2650 dits[10] = 30.5;
e4f08ac3 2651 dits[11] = 43.45;
2652 dits[12] = -7.85;
9d2da7b3 2653 dits[13] = 30.5;
e4f08ac3 2654 dits[14] = 43.45;
2655 dits[15] = -7.85;
9d2da7b3 2656 dits[16] = 30.5;
2657 dits[17] = 36.5;
e4f08ac3 2658 dits[18] = -7;
9d2da7b3 2659 dits[19] = 30.5;
2660 dits[20] = 36.5;
e4f08ac3 2661 dits[21] = -4;
9d2da7b3 2662 dits[22] = 33.0173;
2663 dits[23] = 36.5;
e4f08ac3 2664 dits[24] = -4;
9d2da7b3 2665 dits[25] = 33.0173;
e4f08ac3 2666 dits[26] = 36.80;
2667 dits[27] = -2;
9d2da7b3 2668 dits[28] = 34.6955;
e4f08ac3 2669 dits[29] = 36.80;
2670 gMC->Gsvolu("I216", "PCON", idtmed[272], dits, 30); // supports (1-6) of the ladders
9d2da7b3 2671
2672
2673 // --- Place SPD (option 'a') volumes into their mother volume IT12
2674
2675 // SPD - option 'a'
2676 // (this is NOT the default)
2677
2678 if (option == 1) {
2679
2680 gMC->Gspos("I12A",5,"IT12",0.0,0.0,0.0,idrotm[238],"MANY");
2681 gMC->Gspos("I12A",6,"IT12",0.0,0.0,0.0,idrotm[236],"MANY");
2682 gMC->Gspos("I12A",7,"IT12",0.0,0.0,0.0,idrotm[239],"MANY");
2683 gMC->Gspos("I12A",8,"IT12",0.0,0.0,0.0,idrotm[233],"MANY");
2684 gMC->Gspos("I12A",9,"IT12",0.0,0.0,0.0,idrotm[240],"MANY");
2685 gMC->Gspos("I12A",10,"IT12",0.0,0.0,0.0,idrotm[241],"MANY");
2686 gMC->Gspos("I12A",2,"IT12",0.0,0.0,0.0,idrotm[242],"MANY");
2687 gMC->Gspos("I12A",3,"IT12",0.0,0.0,0.0,idrotm[234],"MANY");
2688 gMC->Gspos("I12A",4,"IT12",0.0,0.0,0.0,idrotm[243],"MANY");
2689 gMC->Gspos("I12A",1,"IT12",0.0,0.0,0.0,0,"MANY");
2690 deltax=((ddet1-0.01/2.)+(dchip1-0.015/2.))*TMath::Cos(270.*TMath::Pi()/180.); // see definition of idrotm[244]
2691 deltay=((ddet1-0.01/2.)+(dchip1-0.015/2.))*TMath::Sin(270.*TMath::Pi()/180.); // see definition of idrotm[244]
2692 gMC->Gspos("I10A",2,"I12A",0.203+deltax,3.8206+deltay,0.0,idrotm[244],"ONLY");
2693 deltax=((ddet1-0.01/2.)+(dchip1-0.015/2.))*TMath::Cos(252.*TMath::Pi()/180.); // see definition of idrotm[245]
2694 deltay=((ddet1-0.01/2.)+(dchip1-0.015/2.))*TMath::Sin(252.*TMath::Pi()/180.); // see definition of idrotm[245]
2695 gMC->Gspos("I10A",1,"I12A",1.4531+deltax,3.8152+deltay,0.0,idrotm[245],"ONLY");
2696 deltax=((ddet2-0.01/2.)+(dchip2-0.015/2.))*TMath::Cos(40.*TMath::Pi()/180.); // see definition of idrotm[246]
2697 deltay=((ddet2-0.01/2.)+(dchip2-0.015/2.))*TMath::Sin(40.*TMath::Pi()/180.); // see definition of idrotm[246]
2698 gMC->Gspos("I20A",1,"I12A",3.0174+deltax,6.5143+deltay,0.0,idrotm[246],"ONLY");
2699 deltax=((ddet2-0.01/2.)+(dchip2-0.015/2.))*TMath::Cos(49.*TMath::Pi()/180.); // see definition of idrotm[247]
2700 deltay=((ddet2-0.01/2.)+(dchip2-0.015/2.))*TMath::Sin(49.*TMath::Pi()/180.); // see definition of idrotm[247]
2701 gMC->Gspos("I20A",2,"I12A",1.9612+deltax,6.9062+deltay,0.0,idrotm[247],"ONLY");
2702 deltax=((ddet2-0.01/2.)+(dchip2-0.015/2.))*TMath::Cos(58.*TMath::Pi()/180.); // see definition of idrotm[248]
2703 deltay=((ddet2-0.01/2.)+(dchip2-0.015/2.))*TMath::Sin(58.*TMath::Pi()/180.); // see definition of idrotm[248]
2704 gMC->Gspos("I20A",3,"I12A",0.8567+deltax,7.1279+deltay,0.0,idrotm[248],"ONLY");
2705 deltax=((ddet2-0.01/2.)+(dchip2-0.015/2.))*TMath::Cos(67.*TMath::Pi()/180.); // see definition of idrotm[249]
2706 deltay=((ddet2-0.01/2.)+(dchip2-0.015/2.))*TMath::Sin(67.*TMath::Pi()/180.); // see definition of idrotm[249]
2707 gMC->Gspos("I20A",4,"I12A",-0.2689+deltax,7.1742+deltay,0.0,idrotm[249],"ONLY");
2708 gMC->Gspos("I123",2,"I12A",-0.2978,5.5196,0.0,idrotm[214],"ONLY");
2709 gMC->Gspos("I121",2,"I12A",-0.2385,4.1518,0.0,idrotm[213],"ONLY");
2710 gMC->Gspos("I122",2,"I12A",-0.2968,4.0207,0.0,idrotm[212],"ONLY");
2711 gMC->Gspos("I120",2,"I12A",-0.3672,3.9056,0.0,0,"ONLY");
2712 gMC->Gspos("I144",1,"I12A",-0.2538,3.8556,0.0,0,"ONLY");
2713 gMC->Gspos("I113",3,"I12A",0.1095,3.9056,0.0,0,"ONLY");
2714 gMC->Gspos("I143",1,"I12A",0.4365,3.8556,0.0,idrotm[236],"ONLY");
2715 gMC->Gspos("I142",1,"I12A",0.5136,3.9056,0.0,idrotm[235],"ONLY");
2716 gMC->Gspos("I141",1,"I12A",0.5636,3.9752,0.0,idrotm[201],"ONLY");
2717 gMC->Gspos("I140",1,"I12A",0.6336,4.0447,0.0,idrotm[234],"ONLY");
2718 gMC->Gspos("I139",1,"I12A",0.8297,4.0545,0.0,idrotm[207],"ONLY");
2719 gMC->Gspos("I113",5,"I12A",1.2575,3.9681,0.0,idrotm[207],"ONLY");
2720 gMC->Gspos("I138",1,"I12A",1.66,3.7848,0.0,idrotm[207],"ONLY");
2721 gMC->Gspos("I137",1,"I12A",1.8556,3.7738,0.0,idrotm[233],"ONLY");
2722 gMC->Gspos("I136",1,"I12A",2.6224,4.874,0.0,idrotm[232],"ONLY");
2723 gMC->Gspos("I135",1,"I12A",3.2967,6.0337,0.0,idrotm[231],"ONLY");
2724 gMC->Gspos("I134",1,"I12A",3.266,6.1636,0.0,idrotm[230],"ONLY");
2725 gMC->Gspos("I113",1,"I12A",2.9903,6.4144,0.0,idrotm[211],"ONLY");
2726 gMC->Gspos("I133",3,"I12A",2.7631,6.7627,0.0,idrotm[230],"ONLY");
2727 gMC->Gspos("I132",3,"I12A",2.62,6.8555,0.0,idrotm[229],"ONLY");
2728 gMC->Gspos("I131",3,"I12A",2.648,6.6023,0.0,idrotm[228],"ONLY");
2729 gMC->Gspos("I130",3,"I12A",2.6569,6.3431,0.0,idrotm[227],"ONLY");
2730 gMC->Gspos("I129",3,"I12A",2.3906,6.4819,0.0,idrotm[226],"ONLY");
2731 gMC->Gspos("I113",2,"I12A",1.9488,6.7998,0.0,idrotm[210],"ONLY");
2732 gMC->Gspos("I133",2,"I12A",1.6699,7.1085,0.0,idrotm[226],"ONLY");
2733 gMC->Gspos("I132",2,"I12A",1.5142,7.1777,0.0,idrotm[225],"ONLY");
2734 gMC->Gspos("I131",2,"I12A",1.5814,6.932,0.0,idrotm[224],"ONLY");
2735 gMC->Gspos("I130",2,"I12A",1.6308,6.6774,0.0,idrotm[223],"ONLY");
2736 gMC->Gspos("I129",2,"I12A",1.346,6.7728,0.0,idrotm[222],"ONLY");
2737 gMC->Gspos("I113",6,"I12A",0.8599,7.0176,0.0,idrotm[209],"ONLY");
2738 gMC->Gspos("I133",1,"I12A",0.5362,7.2789,0.0,idrotm[222],"ONLY");
2739 gMC->Gspos("I132",1,"I12A",0.3715,7.3228,0.0,idrotm[221],"ONLY");
2740 gMC->Gspos("I131",1,"I12A",0.4763,7.0907,0.0,idrotm[220],"ONLY");
2741 gMC->Gspos("I130",1,"I12A",0.5649,6.8469,0.0,idrotm[219],"ONLY");
2742 gMC->Gspos("I129",1,"I12A",0.2688,6.8966,0.0,idrotm[218],"ONLY");
2743 gMC->Gspos("I113",4,"I12A",-0.2497,7.0624,0.0,idrotm[208],"ONLY");
2744 gMC->Gspos("I128",1,"I12A",-0.6103,7.2698,0.0,idrotm[218],"ONLY");
2745 gMC->Gspos("I126",2,"I12A",-0.7799,7.2874,0.0,idrotm[217],"ONLY");
2746 gMC->Gspos("I125",2,"I12A",-0.6315,7.0883,0.0,idrotm[216],"ONLY");
2747 gMC->Gspos("I124",2,"I12A",-0.4965,6.8742,0.0,idrotm[215],"ONLY");
9bf5e8ff 2748 gMC->Gspos("I103",3,"I10A",-0.05,-di10a[1]+2.*di104[1]+di103[1],-3.536,0,"ONLY");
2749 gMC->Gspos("I103",4,"I10A",-0.05,-di10a[1]+2.*di104[1]+di103[1],-10.708,0,"ONLY");
9d2da7b3 2750 gMC->Gspos("I103",1,"I10A",-0.05,-di10a[1]+2.*di104[1]+di103[1],10.708,0,"ONLY");
2751 gMC->Gspos("I103",2,"I10A",-0.05,-di10a[1]+2.*di104[1]+di103[1],3.536,0,"ONLY");
2752 gMC->Gspos("I105",1,"I10A",-0.05,0.01,-16.844,idrotm[237],"ONLY");
2753 gMC->Gspos("I105",2,"I10A",-0.05,0.01,16.844,0,"ONLY");
2754 gMC->Gspos("I104",1,"I10A",0.0,-di10a[1]+di104[1],0.0,0,"ONLY");
9bf5e8ff 2755 gMC->Gspos("I1D3",3,"I20A",-0.05,-di20a[1]+2.*di104[1]+di1d3[1],-3.536,0,"ONLY");
2756 gMC->Gspos("I1D3",4,"I20A",-0.05,-di20a[1]+2.*di104[1]+di1d3[1],-10.708,0,"ONLY");
9d2da7b3 2757 gMC->Gspos("I1D3",1,"I20A",-0.05,-di20a[1]+2.*di104[1]+di1d3[1],10.708,0,"ONLY");
2758 gMC->Gspos("I1D3",2,"I20A",-0.05,-di20a[1]+2.*di104[1]+di1d3[1],3.536,0,"ONLY");
2759 gMC->Gspos("I105",3,"I20A",-0.05,0.01,-16.844,idrotm[237],"ONLY");
2760 gMC->Gspos("I105",4,"I20A",-0.05,0.01,16.844,0,"ONLY");
2761 gMC->Gspos("I104",2,"I20A",0.0,-di20a[1]+di104[1],0.0,0,"ONLY");
2762 gMC->Gspos("I112",2,"I113",0.25,0.02,0.0,idrotm[206],"ONLY");
2763 gMC->Gspos("I111",2,"I113",0.1318,-0.0008,0.0,idrotm[205],"ONLY");
2764 gMC->Gspos("I118",1,"I113",0.0,-0.0454,0.0,0,"ONLY");
2765 gMC->Gspos("I110",1,"I113",0.0,0.0492,0.0,0,"ONLY");
2766 gMC->Gspos("I114",1,"I113",0.063,0.0042,0.0,idrotm[202],"ONLY");
2767 gMC->Gspos("I115",1,"I113",0.063,0.0042,0.0,idrotm[202],"ONLY");
2768 gMC->Gspos("I115",2,"I113",-0.063,0.0042,0.0,idrotm[201],"ONLY");
2769 gMC->Gspos("I114",2,"I113",-0.063,0.0042,0.0,idrotm[201],"ONLY");
2770 gMC->Gspos("I116",1,"I113",0.0,0.0042,0.0,0,"ONLY");
2771 gMC->Gspos("I111",1,"I113",-0.1318,-0.0008,0.0,idrotm[204],"ONLY");
2772 gMC->Gspos("I112",1,"I113",-0.25,0.02,0.0,idrotm[203],"ONLY");
2773 gMC->Gspos("I101",1,"I103",-0.088,ddet1,0.0,0,"ONLY");
2774 gMC->Gspos("I102",1,"I103",0.0,-dchip1,-2.8,0,"ONLY");
2775 gMC->Gspos("I102",2,"I103",0.0,-dchip1,-1.4,0,"ONLY");
2776 gMC->Gspos("I102",3,"I103",0.0,-dchip1,0.0,0,"ONLY");
2777 gMC->Gspos("I102",4,"I103",0.0,-dchip1,1.4,0,"ONLY");
2778 gMC->Gspos("I102",5,"I103",0.0,-dchip1,2.8,0,"ONLY");
2779 gMC->Gspos("I1D1",1,"I1D3",-0.088,ddet2,0.0,0,"ONLY");
2780 gMC->Gspos("I1D2",1,"I1D3",0.0,-dchip2,-2.8,0,"ONLY");
2781 gMC->Gspos("I1D2",2,"I1D3",0.0,-dchip2,-1.4,0,"ONLY");
2782 gMC->Gspos("I1D2",3,"I1D3",0.0,-dchip2,0.0,0,"ONLY");
2783 gMC->Gspos("I1D2",4,"I1D3",0.0,-dchip2,1.4,0,"ONLY");
2784 gMC->Gspos("I1D2",5,"I1D3",0.0,-dchip2,2.8,0,"ONLY");
2785 gMC->Gspos("I117",1,"I116",0.0,0.0,0.0,0,"ONLY");
2786 gMC->Gspos("ITS1",1,"I101",0.0,0.0,0.0,0,"ONLY");
2787 gMC->Gspos("ITS2",1,"I1D1",0.0,0.0,0.0,0,"ONLY");
2788 gMC->Gspos("I650",16,"IT12",0.0,0.0,22.0,idrotm[1104],"MANY");
2789 gMC->Gspos("I650",20,"IT12",0.0,0.0,22.0,idrotm[1130],"MANY");
2790 gMC->Gspos("I650",18,"IT12",0.0,0.0,22.0,idrotm[1117],"MANY");
2791 gMC->Gspos("I650",1,"IT12",0.0,0.0,22.0,0,"MANY");
2792 gMC->Gspos("I650",4,"IT12",0.0,0.0,22.0,idrotm[1106],"MANY");
2793 gMC->Gspos("I650",6,"IT12",0.0,0.0,22.0,idrotm[1039],"MANY");
2794 gMC->Gspos("I650",8,"IT12",0.0,0.0,22.0,idrotm[1107],"MANY");
2795 gMC->Gspos("I650",10,"IT12",0.0,0.0,22.0,idrotm[1065],"MANY");
2796 gMC->Gspos("I650",12,"IT12",0.0,0.0,22.0,idrotm[1078],"MANY");
2797 gMC->Gspos("I650",14,"IT12",0.0,0.0,22.0,idrotm[1091],"MANY");
2798 gMC->Gspos("I650",19,"IT12",0.0,0.0,-22.0,idrotm[1108],"MANY");
2799 gMC->Gspos("I650",2,"IT12",0.0,0.0,-22.0,idrotm[1109],"MANY");
2800 gMC->Gspos("I650",3,"IT12",0.0,0.0,-22.0,idrotm[1110],"MANY");
2801 gMC->Gspos("I650",5,"IT12",0.0,0.0,-22.0,idrotm[1111],"MANY");
2802 gMC->Gspos("I650",7,"IT12",0.0,0.0,-22.0,idrotm[1112],"MANY");
2803 gMC->Gspos("I650",9,"IT12",0.0,0.0,-22.0,idrotm[1113],"MANY");
2804 gMC->Gspos("I650",11,"IT12",0.0,0.0,-22.0,idrotm[1114],"MANY");
2805 gMC->Gspos("I650",13,"IT12",0.0,0.0,-22.0,idrotm[1115],"MANY");
2806 gMC->Gspos("I650",15,"IT12",0.0,0.0,-22.0,idrotm[1116],"MANY");
2807 gMC->Gspos("I650",17,"IT12",0.0,0.0,-22.0,idrotm[1118],"MANY");
2808 gMC->Gspos("I666",1,"I650",0.0,0.0,0.25,idrotm[1003],"MANY");
2809 gMC->Gspos("I667",1,"I650",0.1102,0.9945,0.45,idrotm[1088],"ONLY");
2810 gMC->Gspos("I669",3,"I650",0.1883,4.0372,-3.2,0,"ONLY");
2811 gMC->Gspos("I671",3,"I650",0.1883,4.0372,0.6,0,"ONLY");
2812 gMC->Gspos("I669",2,"I650",1.3343,4.0609,-3.2,0,"ONLY");
2813 gMC->Gspos("I671",2,"I650",1.3343,4.0609,0.6,0,"ONLY");
2814 gMC->Gspos("I669",6,"I650",2.9567,6.1959,-3.2,idrotm[1089],"ONLY");
2815 gMC->Gspos("I671",6,"I650",2.9567,6.1959,0.6,idrotm[1089],"ONLY");
2816 gMC->Gspos("I669",5,"I650",1.9511,6.5822,-3.2,idrotm[1011],"ONLY");
2817 gMC->Gspos("I671",5,"I650",1.9511,6.5822,0.6,idrotm[1011],"ONLY");
2818 gMC->Gspos("I669",4,"I650",0.8974,6.8064,-3.2,idrotm[1090],"ONLY");
2819 gMC->Gspos("I671",4,"I650",0.8974,6.8064,0.6,idrotm[1090],"ONLY");
2820 gMC->Gspos("I669",1,"I650",-0.1784,6.863,-3.2,0,"ONLY");
2821 gMC->Gspos("I671",1,"I650",-0.1784,6.863,0.6,0,"ONLY");
2822 gMC->Gspos("I673",1,"I650",0.2173,4.8037,1.8,0,"ONLY");
2823 gMC->Gspos("I673",6,"I650",1.5093,4.5605,1.8,0,"ONLY");
2824 gMC->Gspos("I673",4,"I650",-0.173,6.2531,1.8,idrotm[1092],"ONLY");
2825 gMC->Gspos("I673",3,"I650",0.8073,6.2032,1.8,idrotm[1093],"ONLY");
2826 gMC->Gspos("I673",2,"I650",1.7678,6.0005,1.8,idrotm[1094],"ONLY");
2827 gMC->Gspos("I673",5,"I650",2.6847,5.6501,1.8,0,"ONLY");
2828 gMC->Gspos("I676",2,"I650",1.7618,5.2269,2.5,0,"ONLY");
2829 gMC->Gspos("I676",1,"I650",0.4018,5.5869,2.5,0,"ONLY");
2830 gMC->Gspos("I668",1,"I667",0.0,0.0,0.0,0,"ONLY");
2831 gMC->Gspos("I670",1,"I669",0.0,0.0,0.0,0,"ONLY");
2832 gMC->Gspos("I672",1,"I671",0.0,0.0,0.0,0,"ONLY");
2833 gMC->Gspos("I674",1,"I673",0.0,0.0,0.0,0,"MANY");
2834 gMC->Gspos("I675",1,"I673",0.0,0.0,-0.5,0,"ONLY");
2835 gMC->Gspos("I677",1,"I676",0.0,0.0,0.0,0,"MANY");
2836 gMC->Gspos("I678",1,"I676",0.0,0.0,-0.95,0,"ONLY");
2837
2838 }
2839
2840
2841 // --- Place SPD (option 'b') volumes into their mother volume IT12
2842
2843 // SPD - option 'b'
2844 // (this is the default)
2845
2846 if (option == 2) {
2847
2848 gMC->Gspos("I12B",1,"IT12",0.0,0.0,0.0,0,"MANY");
2849 gMC->Gspos("I12B",8,"IT12",0.0,0.0,0.0,idrotm[233],"MANY");
2850 gMC->Gspos("I12B",7,"IT12",0.0,0.0,0.0,idrotm[244],"MANY");
2851 gMC->Gspos("I12B",6,"IT12",0.0,0.0,0.0,idrotm[236],"MANY");
2852 gMC->Gspos("I12B",2,"IT12",0.0,0.0,0.0,idrotm[245],"MANY");
2853 gMC->Gspos("I12B",3,"IT12",0.0,0.0,0.0,idrotm[234],"MANY");
2854 gMC->Gspos("I12B",4,"IT12",0.0,0.0,0.0,idrotm[246],"MANY");
2855 gMC->Gspos("I12B",5,"IT12",0.0,0.0,0.0,idrotm[247],"MANY");
2856 gMC->Gspos("I12B",9,"IT12",0.0,0.0,0.0,idrotm[248],"MANY");
2857 gMC->Gspos("I12B",10,"IT12",0.0,0.0,0.0,idrotm[249],"MANY");
2858 deltax=((ddet1-0.01/2.)+(dchip1-0.015/2.))*TMath::Cos(270.*TMath::Pi()/180.); // see definition of idrotm[238]
2859 deltay=((ddet1-0.01/2.)+(dchip1-0.015/2.))*TMath::Sin(270.*TMath::Pi()/180.); // see definition of idrotm[238]
2860 gMC->Gspos("I10B",2,"I12B",0.203+deltax,3.8206+deltay,0.0,idrotm[238],"ONLY");
2861 deltax=((ddet1-0.01/2.)+(dchip1-0.015/2.))*TMath::Cos(252.*TMath::Pi()/180.); // see definition of idrotm[239]
2862 deltay=((ddet1-0.01/2.)+(dchip1-0.015/2.))*TMath::Sin(252.*TMath::Pi()/180.); // see definition of idrotm[239]
2863 gMC->Gspos("I10B",1,"I12B",1.4531+deltax,3.8152+deltay,0.0,idrotm[239],"ONLY");
2864 deltax=((ddet2-0.01/2.)+(dchip2-0.015/2.))*TMath::Cos(40.*TMath::Pi()/180.); // see definition of idrotm[240]
2865 deltay=((ddet2-0.01/2.)+(dchip2-0.015/2.))*TMath::Sin(40.*TMath::Pi()/180.); // see definition of idrotm[240]
2866 gMC->Gspos("I20B",1,"I12B",3.0174+deltax,6.5143+deltay,0.0,idrotm[240],"ONLY");
2867 deltax=((ddet2-0.01/2.)+(dchip2-0.015/2.))*TMath::Cos(49.*TMath::Pi()/180.); // see definition of idrotm[241]
2868 deltay=((ddet2-0.01/2.)+(dchip2-0.015/2.))*TMath::Sin(49.*TMath::Pi()/180.); // see definition of idrotm[241]
2869 gMC->Gspos("I20B",2,"I12B",1.9612+deltax,6.9062+deltay,0.0,idrotm[241],"ONLY");
2870 deltax=((ddet2-0.01/2.)+(dchip2-0.015/2.))*TMath::Cos(58.*TMath::Pi()/180.); // see definition of idrotm[242]
2871 deltay=((ddet2-0.01/2.)+(dchip2-0.015/2.))*TMath::Sin(58.*TMath::Pi()/180.); // see definition of idrotm[242]
2872 gMC->Gspos("I20B",3,"I12B",0.8567+deltax,7.1279+deltay,0.0,idrotm[242],"ONLY");
2873 deltax=((ddet2-0.01/2.)+(dchip2-0.015/2.))*TMath::Cos(67.*TMath::Pi()/180.); // see definition of idrotm[243]
2874 deltay=((ddet2-0.01/2.)+(dchip2-0.015/2.))*TMath::Sin(67.*TMath::Pi()/180.); // see definition of idrotm[243]
2875 gMC->Gspos("I20B",4,"I12B",-0.2689+deltax,7.1742+deltay,0.0,idrotm[243],"ONLY");
2876 gMC->Gspos("I123",1,"I12B",-0.2978,5.5196,0.0,idrotm[214],"ONLY");
2877 gMC->Gspos("I121",1,"I12B",-0.2385,4.1518,0.0,idrotm[213],"ONLY");
2878 gMC->Gspos("I122",1,"I12B",-0.2968,4.0207,0.0,idrotm[212],"ONLY");
2879 gMC->Gspos("I120",1,"I12B",-0.3672,3.9056,0.0,0,"ONLY");
2880 gMC->Gspos("I144",1,"I12B",-0.2538,3.8556,0.0,0,"ONLY");
2881 gMC->Gspos("I113",3,"I12B",0.1095,3.9056,0.0,0,"ONLY");
2882 gMC->Gspos("I143",1,"I12B",0.4365,3.8556,0.0,idrotm[236],"ONLY");
2883 gMC->Gspos("I142",1,"I12B",0.5136,3.9056,0.0,idrotm[235],"ONLY");
2884 gMC->Gspos("I141",1,"I12B",0.5636,3.9752,0.0,idrotm[237],"ONLY");
2885 gMC->Gspos("I140",1,"I12B",0.6336,4.0447,0.0,idrotm[234],"ONLY");
2886 gMC->Gspos("I139",1,"I12B",0.8297,4.0545,0.0,idrotm[207],"ONLY");
2887 gMC->Gspos("I113",5,"I12B",1.2575,3.9681,0.0,idrotm[207],"ONLY");
2888 gMC->Gspos("I138",1,"I12B",1.66,3.7848,0.0,idrotm[207],"ONLY");
2889 gMC->Gspos("I137",1,"I12B",1.8556,3.7738,0.0,idrotm[233],"ONLY");
2890 gMC->Gspos("I136",1,"I12B",2.6224,4.874,0.0,idrotm[232],"ONLY");
2891 gMC->Gspos("I135",1,"I12B",3.2967,6.0337,0.0,idrotm[231],"ONLY");
2892 gMC->Gspos("I134",1,"I12B",3.266,6.1636,0.0,idrotm[230],"ONLY");
2893 gMC->Gspos("I113",1,"I12B",2.9903,6.4144,0.0,idrotm[211],"ONLY");
2894 gMC->Gspos("I133",3,"I12B",2.7631,6.7627,0.0,idrotm[230],"ONLY");
2895 gMC->Gspos("I132",3,"I12B",2.62,6.8555,0.0,idrotm[229],"ONLY");
2896 gMC->Gspos("I131",3,"I12B",2.648,6.6023,0.0,idrotm[228],"ONLY");
2897 gMC->Gspos("I130",3,"I12B",2.6569,6.3431,0.0,idrotm[227],"ONLY");
2898 gMC->Gspos("I129",3,"I12B",2.3906,6.4819,0.0,idrotm[226],"ONLY");
2899 gMC->Gspos("I113",2,"I12B",1.9488,6.7998,0.0,idrotm[210],"ONLY");
2900 gMC->Gspos("I133",2,"I12B",1.6699,7.1085,0.0,idrotm[226],"ONLY");
2901 gMC->Gspos("I132",2,"I12B",1.5142,7.1777,0.0,idrotm[225],"ONLY");
2902 gMC->Gspos("I131",2,"I12B",1.5814,6.932,0.0,idrotm[224],"ONLY");
2903 gMC->Gspos("I130",2,"I12B",1.6308,6.6774,0.0,idrotm[223],"ONLY");
2904 gMC->Gspos("I129",2,"I12B",1.346,6.7728,0.0,idrotm[222],"ONLY");
2905 gMC->Gspos("I113",6,"I12B",0.8599,7.0176,0.0,idrotm[209],"ONLY");
2906 gMC->Gspos("I133",1,"I12B",0.5362,7.2789,0.0,idrotm[222],"ONLY");
2907 gMC->Gspos("I132",1,"I12B",0.3715,7.3228,0.0,idrotm[221],"ONLY");
2908 gMC->Gspos("I131",1,"I12B",0.4763,7.0907,0.0,idrotm[220],"ONLY");
2909 gMC->Gspos("I130",1,"I12B",0.5649,6.8469,0.0,idrotm[219],"ONLY");
2910 gMC->Gspos("I129",1,"I12B",0.2688,6.8966,0.0,idrotm[218],"ONLY");
2911 gMC->Gspos("I113",4,"I12B",-0.2497,7.0624,0.0,idrotm[208],"ONLY");
2912 gMC->Gspos("I128",1,"I12B",-0.6103,7.2698,0.0,idrotm[218],"ONLY");
2913 gMC->Gspos("I126",1,"I12B",-0.7799,7.2874,0.0,idrotm[217],"ONLY");
2914 gMC->Gspos("I125",1,"I12B",-0.6315,7.0883,0.0,idrotm[216],"ONLY");
2915 gMC->Gspos("I124",1,"I12B",-0.4965,6.8742,0.0,idrotm[215],"ONLY");
2916 gMC->Gspos("I105",3,"I10B",-0.05,-0.01,-16.844,idrotm[201],"ONLY");
2917 gMC->Gspos("I105",4,"I10B",-0.05,-0.01,16.844,0,"ONLY");
2918 gMC->Gspos("I107",2,"I10B",-0.0455,-di10b[1]+di107[1],3.536,0,"ONLY");
2919 gMC->Gspos("I107",1,"I10B",-0.0455,-di10b[1]+di107[1],10.708,0,"ONLY");
9bf5e8ff 2920 gMC->Gspos("I107",4,"I10B",-0.0455,-di10b[1]+di107[1],-10.708,0,"ONLY");
2921 gMC->Gspos("I107",3,"I10B",-0.0455,-di10b[1]+di107[1],-3.536,0,"ONLY");
9d2da7b3 2922 gMC->Gspos("I109",1,"I10B",-0.138,0.015,-16.844,idrotm[201],"ONLY");
2923 gMC->Gspos("I109",2,"I10B",-0.138,0.015,16.844,0,"ONLY");
2924 gMC->Gspos("I108",1,"I10B",-0.138,-di10b[1]+2.*di107[1]+di108[1],0.0,0,"ONLY");
2925 gMC->Gspos("I105",1,"I20B",-0.05,-0.01,-16.844,idrotm[201],"ONLY");
2926 gMC->Gspos("I105",2,"I20B",-0.05,-0.01,16.844,0,"ONLY");
2927 gMC->Gspos("I1D7",2,"I20B",-0.0455,-di20b[1]+di1d7[1],3.536,0,"ONLY");
2928 gMC->Gspos("I1D7",1,"I20B",-0.0455,-di20b[1]+di1d7[1],10.708,0,"ONLY");
2929 gMC->Gspos("I1D7",4,"I20B",-0.0455,-di20b[1]+di1d7[1],-10.708,idrotm[201],"ONLY");
2930 gMC->Gspos("I1D7",3,"I20B",-0.0455,-di20b[1]+di1d7[1],-3.536,idrotm[201],"ONLY");
2931 gMC->Gspos("I109",3,"I20B",-0.138,0.015,-16.844,idrotm[201],"ONLY");
2932 gMC->Gspos("I109",4,"I20B",-0.138,0.015,16.844,0,"ONLY");
2933 gMC->Gspos("I108",2,"I20B",-0.138,-di20b[1]+2.*di1d7[1]+di108[1],0.0,0,"ONLY");
2934 gMC->Gspos("I112",2,"I113",0.25,0.02,0.0,idrotm[206],"ONLY");
2935 gMC->Gspos("I111",2,"I113",0.1318,-0.0008,0.0,idrotm[205],"ONLY");
2936 gMC->Gspos("I118",1,"I113",0.0,-0.0454,0.0,0,"ONLY");
2937 gMC->Gspos("I110",1,"I113",0.0,0.0492,0.0,0,"ONLY");
2938 gMC->Gspos("I114",1,"I113",0.063,0.0042,0.0,idrotm[202],"ONLY");
2939 gMC->Gspos("I115",1,"I113",0.063,0.0042,0.0,idrotm[202],"ONLY");
2940 gMC->Gspos("I115",2,"I113",-0.063,0.0042,0.0,idrotm[237],"ONLY");
2941 gMC->Gspos("I114",2,"I113",-0.063,0.0042,0.0,idrotm[237],"ONLY");
2942 gMC->Gspos("I116",1,"I113",0.0,0.0042,0.0,0,"ONLY");
2943 gMC->Gspos("I111",1,"I113",-0.1318,-0.0008,0.0,idrotm[204],"ONLY");
2944 gMC->Gspos("I112",1,"I113",-0.25,0.02,0.0,idrotm[203],"ONLY");
2945 gMC->Gspos("I106",1,"I107",0.0,-dchip1,-1.4,0,"ONLY");
2946 gMC->Gspos("I106",2,"I107",0.0,-dchip1,0.0,0,"ONLY");
2947 gMC->Gspos("I106",3,"I107",0.0,-dchip1,1.4,0,"ONLY");
2948 gMC->Gspos("I106",4,"I107",0.0,-dchip1,2.8,0,"ONLY");
2949 gMC->Gspos("I106",5,"I107",0.0,-dchip1,-2.8,0,"ONLY");
2950 gMC->Gspos("I101",1,"I107",0.0,ddet1,0.0,0,"ONLY");
2951 gMC->Gspos("I1D6",1,"I1D7",0.0,-dchip2,-1.4,0,"ONLY");
2952 gMC->Gspos("I1D6",2,"I1D7",0.0,-dchip2,0.0,0,"ONLY");
2953 gMC->Gspos("I1D6",3,"I1D7",0.0,-dchip2,1.4,0,"ONLY");
2954 gMC->Gspos("I1D6",4,"I1D7",0.0,-dchip2,2.8,0,"ONLY");
2955 gMC->Gspos("I1D6",5,"I1D7",0.0,-dchip2,-2.8,0,"ONLY");
2956 gMC->Gspos("I1D1",1,"I1D7",0.0,ddet2,0.0,0,"ONLY");
2957 gMC->Gspos("I117",1,"I116",0.0,0.0,0.0,0,"ONLY");
2958 gMC->Gspos("ITS1",1,"I101",0.0,0.0,0.0,0,"ONLY");
2959 gMC->Gspos("ITS2",1,"I1D1",0.0,0.0,0.0,0,"ONLY");
2960 gMC->Gspos("I650",16,"IT12",0.0,0.0,22.0,idrotm[1104],"MANY");
2961 gMC->Gspos("I650",20,"IT12",0.0,0.0,22.0,idrotm[1130],"MANY");
2962 gMC->Gspos("I650",18,"IT12",0.0,0.0,22.0,idrotm[1117],"MANY");
2963 gMC->Gspos("I650",1,"IT12",0.0,0.0,22.0,0,"MANY");
2964 gMC->Gspos("I650",4,"IT12",0.0,0.0,22.0,idrotm[1106],"MANY");
2965 gMC->Gspos("I650",6,"IT12",0.0,0.0,22.0,idrotm[1039],"MANY");
2966 gMC->Gspos("I650",8,"IT12",0.0,0.0,22.0,idrotm[1107],"MANY");
2967 gMC->Gspos("I650",10,"IT12",0.0,0.0,22.0,idrotm[1065],"MANY");
2968 gMC->Gspos("I650",12,"IT12",0.0,0.0,22.0,idrotm[1078],"MANY");
2969 gMC->Gspos("I650",14,"IT12",0.0,0.0,22.0,idrotm[1091],"MANY");
2970 gMC->Gspos("I650",19,"IT12",0.0,0.0,-22.0,idrotm[1108],"MANY");
2971 gMC->Gspos("I650",2,"IT12",0.0,0.0,-22.0,idrotm[1109],"MANY");
2972 gMC->Gspos("I650",3,"IT12",0.0,0.0,-22.0,idrotm[1110],"MANY");
2973 gMC->Gspos("I650",5,"IT12",0.0,0.0,-22.0,idrotm[1111],"MANY");
2974 gMC->Gspos("I650",7,"IT12",0.0,0.0,-22.0,idrotm[1112],"MANY");
2975 gMC->Gspos("I650",9,"IT12",0.0,0.0,-22.0,idrotm[1113],"MANY");
2976 gMC->Gspos("I650",11,"IT12",0.0,0.0,-22.0,idrotm[1114],"MANY");
2977 gMC->Gspos("I650",13,"IT12",0.0,0.0,-22.0,idrotm[1115],"MANY");
2978 gMC->Gspos("I650",15,"IT12",0.0,0.0,-22.0,idrotm[1116],"MANY");
2979 gMC->Gspos("I650",17,"IT12",0.0,0.0,-22.0,idrotm[1118],"MANY");
2980 gMC->Gspos("I666",1,"I650",0.0,0.0,0.25,idrotm[1003],"MANY");
2981 gMC->Gspos("I667",1,"I650",0.1102,0.9945,0.45,idrotm[1088],"ONLY");
2982 gMC->Gspos("I669",3,"I650",0.1883,4.0372,-3.2,0,"ONLY");
2983 gMC->Gspos("I671",3,"I650",0.1883,4.0372,0.6,0,"ONLY");
2984 gMC->Gspos("I669",2,"I650",1.3343,4.0609,-3.2,0,"ONLY");
2985 gMC->Gspos("I671",2,"I650",1.3343,4.0609,0.6,0,"ONLY");
2986 gMC->Gspos("I669",6,"I650",2.9567,6.1959,-3.2,idrotm[1089],"ONLY");
2987 gMC->Gspos("I671",6,"I650",2.9567,6.1959,0.6,idrotm[1089],"ONLY");
2988 gMC->Gspos("I669",5,"I650",1.9511,6.5822,-3.2,idrotm[1011],"ONLY");
2989 gMC->Gspos("I671",5,"I650",1.9511,6.5822,0.6,idrotm[1011],"ONLY");
2990 gMC->Gspos("I669",4,"I650",0.8974,6.8064,-3.2,idrotm[1090],"ONLY");
2991 gMC->Gspos("I671",4,"I650",0.8974,6.8064,0.6,idrotm[1090],"ONLY");
2992 gMC->Gspos("I669",1,"I650",-0.1784,6.863,-3.2,0,"ONLY");
2993 gMC->Gspos("I671",1,"I650",-0.1784,6.863,0.6,0,"ONLY");
2994 gMC->Gspos("I673",1,"I650",0.2173,4.8037,1.8,0,"ONLY");
2995 gMC->Gspos("I673",6,"I650",1.5093,4.5605,1.8,0,"ONLY");
2996 gMC->Gspos("I673",4,"I650",-0.173,6.2531,1.8,idrotm[1092],"ONLY");
2997 gMC->Gspos("I673",3,"I650",0.8073,6.2032,1.8,idrotm[1093],"ONLY");
2998 gMC->Gspos("I673",2,"I650",1.7678,6.0005,1.8,idrotm[1094],"ONLY");
2999 gMC->Gspos("I673",5,"I650",2.6847,5.6501,1.8,0,"ONLY");
3000 gMC->Gspos("I676",2,"I650",1.7618,5.2269,2.5,0,"ONLY");
3001 gMC->Gspos("I676",1,"I650",0.4018,5.5869,2.5,0,"ONLY");
3002 gMC->Gspos("I668",1,"I667",0.0,0.0,0.0,0,"ONLY");
3003 gMC->Gspos("I670",1,"I669",0.0,0.0,0.0,0,"ONLY");
3004 gMC->Gspos("I672",1,"I671",0.0,0.0,0.0,0,"ONLY");
3005 gMC->Gspos("I674",1,"I673",0.0,0.0,0.0,0,"MANY");
3006 gMC->Gspos("I675",1,"I673",0.0,0.0,-0.5,0,"ONLY");
3007 gMC->Gspos("I677",1,"I676",0.0,0.0,0.0,0,"MANY");
3008 gMC->Gspos("I678",1,"I676",0.0,0.0,-0.95,0,"ONLY");
3009
3010 }
db915fea 3011
3012 // --- Place SDD volumes into their mother volume IT34
9d2da7b3 3013
db915fea 3014
9d2da7b3 3015 // -- position SDD detectors of ladder 3 / layer 3
3016
3017 gMC->Gspos("ITS3", 1,"I302", 0.0, 0.0, 0.0, 0, "ONLY");
3018 ySDD = Y_SDD_sep/2.+I302dits[1];
3019 for (iSDD=0; iSDD<6; iSDD++) {
3020 gMC->Gspos("I302", iSDD+1, "I004", 0.0, ySDD, Z_SDD_lay3[iSDD], 0, "ONLY");
3021 ySDD = -ySDD;
3022 }
3023
3024 gMC->Gspos("I004", 1,"IT34", -3.2777, 14.3607, 0.0, idrotm[321],"ONLY");
3025 gMC->Gspos("I004", 2,"IT34", -9.5581, 11.9855, 0.0, idrotm[333],"ONLY");
3026 gMC->Gspos("I004", 3,"IT34",-13.2713, 6.3911, 0.0, idrotm[336],"ONLY");
3027 gMC->Gspos("I004", 4,"IT34",-15.33, 0.0, 0.0, idrotm[350],"ONLY");
3028 gMC->Gspos("I004", 5,"IT34",-13.2713, -6.3911, 0.0, idrotm[313],"ONLY");
3029 gMC->Gspos("I004", 6,"IT34", -9.5581, -11.9855, 0.0, idrotm[311],"ONLY");
3030 gMC->Gspos("I004", 7,"IT34", -3.2777, -14.3607, 0.0, idrotm[310],"ONLY");
3031 gMC->Gspos("I004", 8,"IT34", 3.4112, -14.9456, 0.0, idrotm[386],"ONLY");
3032 gMC->Gspos("I004", 9,"IT34", 9.184, -11.5164, 0.0, idrotm[309],"ONLY");
3033 gMC->Gspos("I004",10,"IT34", 13.8119, -6.6514, 0.0, idrotm[308],"ONLY");
3034 gMC->Gspos("I004",11,"IT34", 14.73, 0.0, 0.0, idrotm[356],"ONLY");
3035 gMC->Gspos("I004",12,"IT34", 13.8119, 6.6514, 0.0, idrotm[307],"ONLY");
3036 gMC->Gspos("I004",13,"IT34", 9.184, 11.5164, 0.0, idrotm[306],"ONLY");
3037 gMC->Gspos("I004",14,"IT34", 3.4113, 14.9456, 0.0, idrotm[305],"ONLY");
3038
3039
3040 // -- position SDD detectors of ladder 4 / layer 4
3041
3042 gMC->Gspos("ITS4", 1,"I402", 0.0, 0.000, 0.0, 0,"ONLY");
3043 ySDD = -(Y_SDD_sep/2.+I402dits[1]);
3044 for (iSDD=0; iSDD<8; iSDD++) {
3045 gMC->Gspos("I402", iSDD+1, "I005", 0.0, ySDD, Z_SDD_lay4[iSDD], 0, "ONLY");
3046 ySDD = -ySDD;
3047 }
db915fea 3048
9d2da7b3 3049 gMC->Gspos("I005", 1,"IT34", -3.3629, 23.3895,-0.15, idrotm[335],"ONLY");
3050 gMC->Gspos("I005", 2,"IT34",-10.0447, 21.9949,-0.15, idrotm[332],"ONLY");
3051 gMC->Gspos("I005", 3,"IT34",-15.4744, 17.8584,-0.15, idrotm[331],"ONLY");
3052 gMC->Gspos("I005", 4,"IT34",-20.3415, 13.0727,-0.15, idrotm[366],"ONLY");
3053 gMC->Gspos("I005", 5,"IT34",-22.6728, 6.6573,-0.15, idrotm[330],"ONLY");
3054 gMC->Gspos("I005", 6,"IT34",-24.18, 0.0, -0.15, idrotm[350],"ONLY");
3055 gMC->Gspos("I005", 7,"IT34",-22.6728, -6.6573,-0.15, idrotm[329],"ONLY");
3056 gMC->Gspos("I005", 8,"IT34",-20.3415, -13.0727,-0.15, idrotm[328],"ONLY");
3057 gMC->Gspos("I005", 9,"IT34",-15.4744, -17.8584,-0.15, idrotm[327],"ONLY");
3058 gMC->Gspos("I005",10,"IT34",-10.0447, -21.9949,-0.15, idrotm[326],"ONLY");
3059 gMC->Gspos("I005",11,"IT34", -3.3629, -23.3895,-0.15, idrotm[325],"ONLY");
3060 gMC->Gspos("I005",12,"IT34", 3.4412, -23.9339,-0.15, idrotm[324],"ONLY");
3061 gMC->Gspos("I005",13,"IT34", 9.8163, -21.4946,-0.15, idrotm[323],"ONLY");
3062 gMC->Gspos("I005",14,"IT34", 15.8345, -18.274, -0.15, idrotm[322],"ONLY");
3063 gMC->Gspos("I005",15,"IT34", 19.8788, -12.7753,-0.15, idrotm[320],"ONLY");
3064 gMC->Gspos("I005",16,"IT34", 23.2005, -6.8123,-0.15, idrotm[319],"ONLY");
3065 gMC->Gspos("I005",17,"IT34", 23.63, 0.0, -0.15, idrotm[318],"ONLY");
3066 gMC->Gspos("I005",18,"IT34", 23.2005, 6.8123,-0.15, idrotm[317],"ONLY");
3067 gMC->Gspos("I005",19,"IT34", 19.8788, 12.7753,-0.15, idrotm[316],"ONLY");
3068 gMC->Gspos("I005",20,"IT34", 15.8345, 18.274, -0.15, idrotm[315],"ONLY");
3069 gMC->Gspos("I005",21,"IT34", 9.8163, 21.4946,-0.15, idrotm[314],"ONLY");
3070 gMC->Gspos("I005",22,"IT34", 3.4412, 23.9339,-0.15, idrotm[334],"ONLY");
3071
3072
3073 // -- build block of the SDD ladder frame holding the electronics
3074
3075 gMC->Gspos("I019", 1,"I018", -1.9, -1.735, 0.0, idrotm[344], "ONLY");
3076 gMC->Gspos("I019", 2,"I018", 1.987, -1.5843, 0.0, idrotm[343], "ONLY");
3077 gMC->Gspos("I019", 3,"I018", -0.087, 1.7066, 0.0, idrotm[342], "ONLY");
3078
3079 gMC->Gspos("I020", 1,"I018", -1.9782, -1.569, 0.0, idrotm[342], "ONLY");
3080 gMC->Gspos("I020", 2,"I018", 1.8824, -1.735, 0.0, idrotm[344], "ONLY");
3081 gMC->Gspos("I020", 3,"I018", 0.0958, 1.6913, 0.0, idrotm[343], "ONLY");
3082
3083 gMC->Gspos("I021", 1,"I018", 1.0761, 0.0835, 2.6008, idrotm[340], "ONLY");
3084 gMC->Gspos("I021", 2,"I018", -1.0761, 0.0835,-2.8008, idrotm[339], "ONLY");
3085 gMC->Gspos("I021", 3,"I018", -1.0761, 0.0835,-1.0492, idrotm[338], "ONLY");
3086 gMC->Gspos("I021", 4,"I018", 1.0761, 0.0835,-2.8008, idrotm[337], "ONLY");
3087 gMC->Gspos("I021", 5,"I018", 1.0761, 0.0835,-1.0492, idrotm[340], "ONLY");
3088 gMC->Gspos("I021", 6,"I018", -1.0761, 0.0835, 0.8492, idrotm[339], "ONLY");
3089 gMC->Gspos("I021", 7,"I018", -1.0761, 0.0835, 2.6008, idrotm[338], "ONLY");
3090 gMC->Gspos("I021", 8,"I018", 1.0761, 0.0835, 0.8492, idrotm[337], "ONLY");
3091
3092 gMC->Gspos("I022", 1,"I018", 0.0, -1.79, 3.55, idrotm[312], "ONLY");
3093 gMC->Gspos("I022", 2,"I018", 0.0, -1.79, -0.1, idrotm[312], "ONLY");
3094
3095 gMC->Gspos("I023", 1,"I018", 0.0, -1.79, 1.725, idrotm[341], "ONLY");
3096 gMC->Gspos("I023", 2,"I018", 0.0, -1.79, -1.925, idrotm[341], "ONLY");
3097
3098 gMC->Gspos("I033", 1,"I018", 1.8, -1.75, 1.35, 0, "MANY");
3099 gMC->Gspos("I033", 2,"I018", -1.8, -1.75, -2.65, idrotm[345], "MANY");
3100 gMC->Gspos("I033", 3,"I018", -1.8, -1.75, 1.35, idrotm[345], "MANY");
3101 gMC->Gspos("I033", 4,"I018", 1.8, -1.75, -2.65, 0, "MANY");
3102
3103 gMC->Gspos("I034", 1,"I018", 1.6, -1.775, 1.35, idrotm[312], "ONLY");
3104 gMC->Gspos("I034", 2,"I018", -1.6, -1.775, -2.65, idrotm[348], "ONLY");
3105 gMC->Gspos("I034", 3,"I018", -1.6, -1.775, 1.35, idrotm[348], "ONLY");
3106 gMC->Gspos("I034", 4,"I018", 1.6, -1.775, -2.65, idrotm[312], "ONLY");
3107
3108 gMC->Gspos("I035", 1,"I018", 1.7, -0.55, I018dits[2]-I035dits[2], 0, "MANY");
3109 gMC->Gspos("I035", 2,"I018", -1.7, -0.55, I018dits[2]-I035dits[2], 0, "MANY");
3110
3111 gMC->Gspos("I036", 1,"I018", 0.3087, 1.7191, 3.56, idrotm[346], "ONLY");
3112 gMC->Gspos("I036", 2,"I018", 0.3087, 1.7191,-0.11, idrotm[346], "ONLY");
3113 gMC->Gspos("I036", 3,"I018", -0.3087, 1.7191,-0.11, idrotm[347], "ONLY");
3114 gMC->Gspos("I036", 4,"I018", -0.3087, 1.7191, 3.56, idrotm[347], "ONLY");
3115
3116 gMC->Gspos("I037", 1,"I018", SDD_CoolPipe[0], SDD_CoolPipe[1], SDD_CoolPipe[2], 0 , "ONLY");
3117 gMC->Gspos("I037", 2,"I018", -SDD_CoolPipe[0], SDD_CoolPipe[1], SDD_CoolPipe[2], 0 , "ONLY");
3118
3119 gMC->Gspos("I038", 1,"I018", SDD_CoolPipe[0], SDD_CoolPipe[1], SDD_CoolPipe[2], 0 , "ONLY");
3120 gMC->Gspos("I038", 2,"I018", -SDD_CoolPipe[0], SDD_CoolPipe[1], SDD_CoolPipe[2], 0 , "ONLY");
3121
3122 gMC->Gspos("I040", 1,"I018", 1.9204, -0.7118, 0.0, idrotm[346],"ONLY");
3123 gMC->Gspos("I040", 2,"I018", -1.9204, -0.7118, 0.0, idrotm[347],"ONLY");
3124 gMC->Gspos("I041", 1,"I018", SDD_CoolPipe[0], SDD_CoolPipe[1], SDD_CoolPipe[2], idrotm[346], "ONLY");
3125 gMC->Gspos("I041", 2,"I018", -SDD_CoolPipe[0], SDD_CoolPipe[1], SDD_CoolPipe[2], idrotm[347], "ONLY");
3126
3127
3128 // -- build block of the SDD electronics (heat bridge, chips, hybrid, anode microcable)
3129
3130 xI050 = SDD_CoolPipe[0]+SDD_CoolPipe[1]*sin30+I050dits[1]/cos30+I041dits[1];
3131 yI050 = 0;
3132 xI039 = -SDD_CoolPipe[1]/cos30;
3133 yI039 = -I050dits[1]+I039dits[1];
3134 gMC->Gspos("I039", 1,"I050", xI039, yI039, 0.0, 0, "ONLY");
3135 xI042 = xI039+I039dits[0]-xI042space-I042dits[0];
3136 yI042 = yI039+I039dits[1]+I042dits[1];
3137 xI043 = xI039-I039dits[0]+xI043space+I043dits[0];
3138 yI043 = yI039+I039dits[1]+I043dits[1];
3139 zChipSpace = I042dits[2];
3140 if (zChipSpace < I043dits[2]) {
3141 zChipSpace = I043dits[2];
3142 }
3143 zChipSpace = zChipSpace * 2;
3144 yI051space = (2*I039dits[2] - 4*zChipSpace)/5;
3145 zchip = -I039dits[2] + yI051space + zChipSpace/2.;
3146 for (ichip=0; ichip<4; ichip++) {
3147 gMC->Gspos("I042", ichip+1, "I050", xI042, yI042, zchip, 0, "ONLY");
3148 gMC->Gspos("I043", ichip+1, "I050", xI043, yI043, zchip, 0, "ONLY");
3149 zchip += zChipSpace + yI051space;
3150 }
3151 xcap = 2*I039dits[0]/5.;
3152 yI051 = yI039+I039dits[1]+I051dits[1];
3153 zI051 = -I039dits[2] + yI051space/3.;
3154 icap = 1;
3155 for (ichip=0; ichip<5; ichip++) {
3156 xI051 = xI039-I039dits[0]+xcap;
3157 gMC->Gspos("I051", icap++,"I050", xI051, yI051, zI051, 0, "ONLY");
3158 zI051 += yI051space/3.;
3159 gMC->Gspos("I051", icap++,"I050", xI051, yI051, zI051, 0, "ONLY");
3160 xI051 += xcap;
3161 gMC->Gspos("I051", icap++,"I050", xI051, yI051, zI051, 0, "ONLY");
3162 xI051 += xcap;
3163 gMC->Gspos("I051", icap++,"I050", xI051, yI051, zI051, 0, "ONLY");
3164 xI051 += xcap;
3165 gMC->Gspos("I051", icap++,"I050", xI051, yI051, zI051, 0, "ONLY");
3166 zI051 -= yI051space/3.;
3167 if (ichip == 0) {
3168 gMC->Gspos("I051", icap++,"I050", xI051, yI051, zI051, 0, "ONLY");
3169 }
3170 zI051 += zChipSpace + yI051space;
3171 }
3172 xI052 = -I050dits[0]+I052dits[0];
3173 yI052 = yI051+I051dits[1]+I052dits[1];
3174 gMC->Gspos("I052", 1,"I050", xI052, yI052, 0.0, 0, "ONLY");
3175 xI044 = I050dits[0]-I044dits[3];
3176 yI044 = yI052+I052dits[1]+I044dits[2];
3177 gMC->Gspos("I044", 1,"I050", xI044, yI044, 0.0, idrotm[301], "ONLY");
3178 gMC->Gspos("I050", 1,"I018", xI050, yI050, 0.0, idrotm[346],"ONLY");
3179 gMC->Gspos("I050", 2,"I018", -xI050, yI050, 0.0, idrotm[347],"ONLY");
3180
3181
3182 // -- build block of the SDD ladder frame at the end ladders
3183
3184 gMC->Gspos("I021",12,"I024", 1.0761, 0.0836,-0.1242, idrotm[340], "ONLY");
3185 gMC->Gspos("I021",11,"I024", -1.0761, 0.0836,-0.1242, idrotm[338], "ONLY");
3186 gMC->Gspos("I021",13,"I024", -1.0761, 0.0836,-1.8758, idrotm[339], "ONLY");
3187 gMC->Gspos("I021",14,"I024", 1.0761, 0.0836,-1.8758, idrotm[337], "ONLY");
3188
3189 gMC->Gspos("I022", 3,"I024", 0.0, -1.7899, 0.825, idrotm[312], "ONLY");
3190
3191 gMC->Gspos("I023", 3,"I024", 0.0, -1.7899,-1.0, idrotm[341], "ONLY");
3192
3193 gMC->Gspos("I025", 1,"I024", -1.9, -1.7349, 0.0, idrotm[344], "ONLY");
3194 gMC->Gspos("I025", 2,"I024", 1.987, -1.5842, 0.0, idrotm[343], "ONLY");
3195
3196 gMC->Gspos("I026", 1,"I024", -1.9782, -1.5689, 0.0, idrotm[342], "ONLY");
3197 gMC->Gspos("I026", 2,"I024", 1.8824, -1.7349, 0.0, idrotm[344], "ONLY");
3198
3199 gMC->Gspos("I029", 1,"I024", -0.087, 1.7067, I029dits[2]-I024dits[2], idrotm[342], "ONLY");
3200
3201 gMC->Gspos("I030", 1,"I024", 0.0958, 1.6914, I030dits[2]-I024dits[2], idrotm[343], "ONLY");
3202
3203 gMC->Gspos("I031", 1,"I024", SDD_CoolPipe[0], SDD_CoolPipe[1], SDD_CoolPipe[2], 0 ,"ONLY");
3204 gMC->Gspos("I031", 2,"I024", -SDD_CoolPipe[0], SDD_CoolPipe[1], SDD_CoolPipe[2], 0 ,"ONLY");
3205
3206 gMC->Gspos("I032", 1,"I024", SDD_CoolPipe[0], SDD_CoolPipe[1], SDD_CoolPipe[2], 0 ,"ONLY");
3207 gMC->Gspos("I032", 2,"I024", -SDD_CoolPipe[0], SDD_CoolPipe[1], SDD_CoolPipe[2], 0 ,"ONLY");
3208
3209
3210 xI424 = I028dits[0]/3.;
3211 yI424 = -I028dits[1]+I424dits[1];
3212 gMC->Gspos("I422", 1,"I421", 0.0, 0.0, 0.0, 0, "ONLY");
e4f08ac3 3213 gMC->Gspos("I423", 1,"I421", 0.0, 0.0, I421dits[2]-I423dits[2], 0, "ONLY");
3214 gMC->Gspos("I421", 1,"I420", 0.0, 0.0, 0.0, idrotm[312], "ONLY");
9d2da7b3 3215 gMC->Gspos("I420", 1,"I028", -I028dits[0]/3., I028dits[1]-I420dits[1], 0.0, 0, "ONLY");
3216 gMC->Gspos("I424", 1,"I028", xI424, yI424, 0.0, 0, "ONLY");
3217 gMC->Gspos("I028", 1,"I024", 0.0, I028dits[1]-I024dits[1], I024dits[2]-I028dits[2], 0, "MANY");
3218
3219
3220 // -- build the SDD ladder 3
3221
3222 indI425 = 1;
3223 gMC->Gspos("I024", 1,"I047", 0.0, 0.0, 24.625, 0, "ONLY");
3224 gMC->Gspos("I018", 1,"I047", 0.0, 0.0, 3.65, 0, "ONLY");
3225 gMC->Gspos("I018", 2,"I047", 0.0, 0.0, 10.95, 0, "ONLY");
3226 gMC->Gspos("I018", 3,"I047", 0.0, 0.0, 18.25, 0, "ONLY");
3227 gMC->Gspos("I018", 4,"I047", 0.0, 0.0, -3.65, 0, "ONLY");
3228 gMC->Gspos("I018", 5,"I047", 0.0, 0.0, -10.95, 0, "ONLY");
3229 gMC->Gspos("I018", 6,"I047", 0.0, 0.0, -18.25, 0, "ONLY");
3230 gMC->Gspos("I024", 2,"I047", 0.0, 0.0, -24.625, idrotm[355], "ONLY");
3231 HVname[0] = 'I';
3232 HVname[1] = '3';
3233 HVname[2] = '1';
3234 HVname[4] = '\0';
3235 for (iSDD=0; iSDD<3; iSDD++) {
3236 HVname[3] = (Char_t)(48+iSDD+5);
3237 dits[0] = 1.350000;
3238 dits[1] = I425dits[1];
3239 dits[2] = (I047dits[2] - 2*I024dits[2] - Z_SDD_lay3[iSDD])/2.;
3240 gMC->Gsvolu(HVname, "BOX ", idtmed[279], dits, 3);
3241 xHV = 0.0;
3242 yHV = -I047dits[1] + (2*iSDD+1)*dits[1];
3243 zHV = I047dits[2] - 2*I024dits[2] - dits[2];
3244 gMC->Gspos(HVname, 1,"I047", xHV, yHV, zHV, 0, "ONLY");
3245 gMC->Gspos(HVname, 2,"I047", xHV, yHV, -zHV, 0, "ONLY");
3246 gMC->Gspos("I425", indI425++,"I047", xI424, yHV, 24.625, 0, "ONLY");
3247 gMC->Gspos("I425", indI425++,"I047", -xI424, yHV, -24.625, 0, "ONLY");
3248 }
3249 LVname[0] = 'I';
3250 LVname[1] = '3';
3251 LVname[2] = '1';
3252 LVname[4] = '\0';
3253 for (iSDD=0; iSDD<3; iSDD++) {
3254 LVname[3] = (Char_t)(48+iSDD+1);
3255 dits[0] = 1.350000;
3256 dits[1] = 0.004423;
3257 dits[2] = (I047dits[2] - (2*iSDD+1)*I018dits[2] - I039dits[2])/2.;
3258 gMC->Gsvolu(LVname, "BOX ", idtmed[280], dits, 3);
3259 yLV = I018dits[1] - dits[0]*cos30 - dits[1]*sin30;
3260 xLV = xI050 -
3261 fabs(yI050-yLV)*sin30/cos30 +
3262 (I050dits[1]+(2*iSDD+1)*dits[1])/cos30;
3263 zLV = I047dits[2] - dits[2];
3264 gMC->Gspos(LVname, 1,"I047", xLV, yLV, zLV, idrotm[346], "ONLY");
3265 gMC->Gspos(LVname, 2,"I047", xLV, yLV, -zLV, idrotm[346], "ONLY");
3266 gMC->Gspos(LVname, 3,"I047", -xLV, yLV, zLV, idrotm[347], "ONLY");
3267 gMC->Gspos(LVname, 4,"I047", -xLV, yLV, -zLV, idrotm[347], "ONLY");
3268 }
3269
3270
3271 // -- build the SDD ladder 4
3272
3273
3274 gMC->Gspos("I024", 3,"I048", -0.0001, 0.0, 31.925, 0, "ONLY");
3275 gMC->Gspos("I018", 7,"I048", -0.0001, 0.0, -3.65, 0, "ONLY");
3276 gMC->Gspos("I018", 8,"I048", -0.0001, 0.0, 3.65, 0, "ONLY");
3277 gMC->Gspos("I018", 9,"I048", -0.0001, 0.0, 10.95, 0, "ONLY");
3278 gMC->Gspos("I018",10,"I048", -0.0001, 0.0, 18.25, 0, "ONLY");
3279 gMC->Gspos("I018",11,"I048", -0.0001, 0.0, 25.55, 0, "ONLY");
3280 gMC->Gspos("I018",12,"I048", -0.0001, 0.0, -10.95, 0, "ONLY");
3281 gMC->Gspos("I018",13,"I048", -0.0001, 0.0, -18.25, 0, "ONLY");
3282 gMC->Gspos("I018",14,"I048", -0.0001, 0.0, -25.55, 0, "ONLY");
3283 gMC->Gspos("I024", 4,"I048", -0.0001, 0.0, -31.925, idrotm[355], "ONLY");
3284 HVname[0] = 'I';
3285 HVname[1] = '4';
3286 HVname[2] = '1';
3287 HVname[4] = '\0';
3288 for (iSDD=0; iSDD<4; iSDD++) {
3289 HVname[3] = (Char_t)(48+iSDD+5);
3290 dits[0] = 1.350000;
3291 dits[1] = I425dits[1];
3292 dits[2] = (I048dits[2] - 2*I024dits[2] - Z_SDD_lay4[iSDD])/2.;
3293 gMC->Gsvolu(HVname, "BOX ", idtmed[279], dits, 3);
3294 xHV = -0.0001;
3295 yHV = -I048dits[1] + (2*iSDD+1)*dits[1];
3296 zHV = I048dits[2] - 2*I024dits[2] - dits[2];
3297 gMC->Gspos(HVname, 1,"I048", xHV, yHV, zHV, 0, "ONLY");
3298 gMC->Gspos(HVname, 2,"I048", xHV, yHV, -zHV, 0, "ONLY");
3299 gMC->Gspos("I425", indI425++,"I048", xI424, yHV, 31.925, 0, "ONLY");
3300 gMC->Gspos("I425", indI425++,"I048", -xI424, yHV, -31.925, 0, "ONLY");
3301 }
3302 LVname[0] = 'I';
3303 LVname[1] = '4';
3304 LVname[2] = '1';
3305 LVname[4] = '\0';
3306 for (iSDD=0; iSDD<4; iSDD++) {
3307 LVname[3] = (Char_t)(48+iSDD+1);
3308 dits[0] = 1.350000;
3309 dits[1] = 0.004423;
3310 dits[2] = (I048dits[2] - (2*iSDD+1)*I018dits[2] - I039dits[2])/2.;
3311 gMC->Gsvolu(LVname, "BOX ", idtmed[280], dits, 3);
3312 yLV = I018dits[1] - dits[0]*cos30 - dits[1]*sin30;
3313 xLV = xI050 -
3314 fabs(yI050-yLV)*sin30/cos30 +
3315 (I050dits[1]+(2*iSDD+1)*dits[1])/cos30;
3316 zLV = I048dits[2] - dits[2];
3317 gMC->Gspos(LVname, 1,"I048", xLV, yLV, zLV, idrotm[346], "ONLY");
3318 gMC->Gspos(LVname, 2,"I048", xLV, yLV, -zLV, idrotm[346], "ONLY");
3319 gMC->Gspos(LVname, 3,"I048", -xLV, yLV, zLV, idrotm[347], "ONLY");
3320 gMC->Gspos(LVname, 4,"I048", -xLV, yLV, -zLV, idrotm[347], "ONLY");
3321 }
3322
db915fea 3323
9d2da7b3 3324 // -- build the SDD barrel (layers 3 and 4)
3325
3326 gMC->Gspos("I047", 1,"IT34", -3.7528, 16.4422, 0.0, idrotm[321], "ONLY");
3327 gMC->Gspos("I047", 2,"IT34",-10.8892, 13.6547, 0.0, idrotm[333], "ONLY");
3328 gMC->Gspos("I047", 3,"IT34",-15.1948, 7.3175, 0.0, idrotm[336], "ONLY");
3329 gMC->Gspos("I047", 4,"IT34",-17.465, 0.0, 0.0, idrotm[350], "ONLY");
3330 gMC->Gspos("I047", 5,"IT34",-15.1948, -7.3174, 0.0, idrotm[313], "ONLY");
3331 gMC->Gspos("I047", 6,"IT34",-10.8893, -13.6547, 0.0, idrotm[311], "ONLY");
3332 gMC->Gspos("I047", 7,"IT34", -3.7528, -16.4422, 0.0, idrotm[310], "ONLY");
3333 gMC->Gspos("I047", 8,"IT34", 3.8863, -17.0271, 0.0, idrotm[386], "ONLY");
3334 gMC->Gspos("I047", 9,"IT34", 10.5152, -13.1856, 0.0, idrotm[309], "ONLY");
3335 gMC->Gspos("I047",10,"IT34", 15.7354, -7.5778, 0.0, idrotm[308], "ONLY");
3336 gMC->Gspos("I047",11,"IT34", 16.865, 0.0, 0.0, idrotm[356], "ONLY");
3337 gMC->Gspos("I047",12,"IT34", 15.7354, 7.5778, 0.0, idrotm[307], "ONLY");
3338 gMC->Gspos("I047",13,"IT34", 10.5152, 13.1856, 0.0, idrotm[306], "ONLY");
3339 gMC->Gspos("I047",14,"IT34", 3.8863, 17.0271, 0.0, idrotm[305], "ONLY");
3340
3341 gMC->Gspos("I048", 1,"IT34", -3.6667, 25.5027, 0.0, idrotm[335], "ONLY");
3342 gMC->Gspos("I048", 2,"IT34",-10.9317, 23.937, 0.0, idrotm[332], "ONLY");
3343 gMC->Gspos("I048", 3,"IT34",-16.8725, 19.4719, 0.0, idrotm[331], "ONLY");
3344 gMC->Gspos("I048", 4,"IT34",-22.1376, 14.227, 0.0, idrotm[366], "ONLY");
3345 gMC->Gspos("I048", 5,"IT34",-24.7213, 7.2588, 0.0, idrotm[330], "ONLY");
3346 gMC->Gspos("I048", 6,"IT34",-26.315, 0.0, 0.0, idrotm[350], "ONLY");
3347 gMC->Gspos("I048", 7,"IT34",-24.7213, -7.2588, 0.0, idrotm[329], "ONLY");
3348 gMC->Gspos("I048", 8,"IT34",-22.1376, -14.227, 0.0, idrotm[328], "ONLY");
3349 gMC->Gspos("I048", 9,"IT34",-16.8725, -19.4719, 0.0, idrotm[327], "ONLY");
3350 gMC->Gspos("I048",10,"IT34",-10.9316, -23.937, 0.0, idrotm[326], "ONLY");
3351 gMC->Gspos("I048",11,"IT34", -3.6667, -25.5027, 0.0, idrotm[325], "ONLY");
3352 gMC->Gspos("I048",12,"IT34", 3.745, -26.0472, 0.0, idrotm[324], "ONLY");
3353 gMC->Gspos("I048",13,"IT34", 10.7032, -23.4367, 0.0, idrotm[323], "ONLY");
3354 gMC->Gspos("I048",14,"IT34", 17.2327, -19.8876, 0.0, idrotm[322], "ONLY");
3355 gMC->Gspos("I048",15,"IT34", 21.6749, -13.9296, 0.0, idrotm[320], "ONLY");
3356 gMC->Gspos("I048",16,"IT34", 25.2491, -7.4138, 0.0, idrotm[319], "ONLY");
3357 gMC->Gspos("I048",17,"IT34", 25.765, 0.0, 0.0, idrotm[318], "ONLY");
3358 gMC->Gspos("I048",18,"IT34", 25.2491, 7.4138, 0.0, idrotm[317], "ONLY");
3359 gMC->Gspos("I048",19,"IT34", 21.6749, 13.9296, 0.0, idrotm[316], "ONLY");
3360 gMC->Gspos("I048",20,"IT34", 17.2327, 19.8876, 0.0, idrotm[315], "ONLY");
3361 gMC->Gspos("I048",21,"IT34", 10.7032, 23.4367, 0.0, idrotm[314], "ONLY");
3362 gMC->Gspos("I048",22,"IT34", 3.745, 26.0472, 0.0, idrotm[334], "ONLY");
3363
3364
3365 // --- Place SSD volumes into their mother volume IT56
db915fea 3366
3367
9d2da7b3 3368 gMC->Gspos("I570",14,"IT56",-28.0681,-36.0619,-0.27,idrotm[566],"ONLY");
db915fea 3369 gMC->Gspos("I570",15,"IT56",-21.677,-40.0556,-0.27,idrotm[567],"ONLY");
3370 gMC->Gspos("I570",16,"IT56",-14.838,-43.2217,-0.27,idrotm[568],"ONLY");
3371 gMC->Gspos("I570",17,"IT56",-7.4965,-44.9238,-0.27,idrotm[569],"ONLY");
9d2da7b3 3372 gMC->Gspos("I570",18,"IT56",-0.27,-45.6977,-0.27,idrotm[533],"ONLY");
db915fea 3373 gMC->Gspos("I570",19,"IT56",7.4965,-44.9238,-0.27,idrotm[534],"ONLY");
3374 gMC->Gspos("I570",20,"IT56",14.838,-43.2217,-0.27,idrotm[535],"ONLY");
3375 gMC->Gspos("I570",21,"IT56",21.677,-40.0556,-0.27,idrotm[623],"ONLY");
3376 gMC->Gspos("I570",22,"IT56",28.0681,-36.0619,-0.27,idrotm[537],"ONLY");
3377 gMC->Gspos("I570",23,"IT56",33.5085,-30.8468,-0.27,idrotm[538],"ONLY");
3378 gMC->Gspos("I570",24,"IT56",38.2566,-24.9943,-0.27,idrotm[539],"ONLY");
3379 gMC->Gspos("I570",25,"IT56",41.7089,-18.2952,-0.27,idrotm[540],"ONLY");
3380 gMC->Gspos("I570",26,"IT56",44.2994,-11.2181,-0.27,idrotm[541],"ONLY");
3381 gMC->Gspos("I570",27,"IT56",45.3894,-3.7611,-0.27,idrotm[542],"ONLY");
3382 gMC->Gspos("I570",28,"IT56",45.5416,3.7737,-0.27,idrotm[543],"ONLY");
3383 gMC->Gspos("I570",29,"IT56",44.1513,11.1806,-0.27,idrotm[544],"ONLY");
3384 gMC->Gspos("I570",30,"IT56",41.8487,18.3566,-0.27,idrotm[545],"ONLY");
3385 gMC->Gspos("I570",31,"IT56",38.1287,24.9107,-0.27,idrotm[546],"ONLY");
3386 gMC->Gspos("I570",32,"IT56",33.6209,30.9502,-0.27,idrotm[547],"ONLY");
3387 gMC->Gspos("I570",33,"IT56",27.9743,35.9414,-0.27,idrotm[548],"ONLY");
3388 gMC->Gspos("I570",34,"IT56",21.7497,40.1899,-0.27,idrotm[549],"ONLY");
3389 gMC->Gspos("I570",35,"IT56",14.7884,43.0772,-0.27,idrotm[550],"ONLY");
3390 gMC->Gspos("I570",36,"IT56",7.5216,45.0744,-0.27,idrotm[551],"ONLY");
8a3132f2 3391 gMC->Gspos("I570",37,"IT56",0.,45.545,-0.27,0,"ONLY");
db915fea 3392 gMC->Gspos("I570",38,"IT56",-7.5216,45.0744,-0.27,idrotm[552],"ONLY");
3393 gMC->Gspos("I570",1,"IT56",-14.7884,43.0772,-0.27,idrotm[553],"ONLY");
3394 gMC->Gspos("I570",2,"IT56",-21.7497,40.1899,-0.27,idrotm[620],"ONLY");
3395 gMC->Gspos("I570",3,"IT56",-27.9743,35.9414,-0.27,idrotm[555],"ONLY");
3396 gMC->Gspos("I570",4,"IT56",-33.6209,30.9502,-0.27,idrotm[556],"ONLY");
3397 gMC->Gspos("I570",5,"IT56",-38.1287,24.9108,-0.27,idrotm[557],"ONLY");
3398 gMC->Gspos("I570",6,"IT56",-41.8487,18.3566,-0.27,idrotm[558],"ONLY");