]> git.uio.no Git - u/mrichter/AliRoot.git/blame - MUON/AliMUONSt1GeometryBuilderV2.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
1. Added detElemId to all constructors; it is kept as the UniqueID of TObject 2.Added...
[u/mrichter/AliRoot.git] / MUON / AliMUONSt1GeometryBuilderV2.cxx
Content-type: text/html ]> git.uio.no Git - u/mrichter/AliRoot.git/blame - MUON/AliMUONSt1GeometryBuilderV2.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:


500 - Internal Server Error
Malformed UTF-8 character (fatal) at (eval 5) line 1, <$fd> line 1389.
CommitLineData
ba030c0e 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
d1cd2474 16// $Id$
17//
5f1df83a 18// Authors: David Guez, Ivana Hrivnacova, Marion MacCormick; IPN Orsay
19//
d1cd2474 20// Class AliMUONSt1GeometryBuilderV2
21// ---------------------------------
22// MUON Station1 detailed geometry construction class.
5f1df83a 23// (Originally defined in AliMUONv2.cxx - now removed.)
24// Included in AliRoot 2004/01/23
d1cd2474 25
26#ifdef ST1_WITH_STL
27 #include <vector>
28#endif
ba030c0e 29
d1cd2474 30#ifdef ST1_WITH_ROOT
31 #include "TArrayI.h"
32#endif
ba030c0e 33
34#include <TVector2.h>
d1cd2474 35#include <TVector3.h>
36#include <TGeoMatrix.h>
ba030c0e 37#include <TClonesArray.h>
5f91c9e8 38#include <Riostream.h>
39#include <TSystem.h>
5d12ce38 40#include <TVirtualMC.h>
ba030c0e 41
5f91c9e8 42#include "AliMpFiles.h"
43#include "AliMpReader.h"
44#include "AliMpSector.h"
45#include "AliMpRow.h"
46#include "AliMpVRowSegment.h"
47#include "AliMpMotifMap.h"
48#include "AliMpMotifPosition.h"
ba030c0e 49
d1cd2474 50#include "AliMUONSt1GeometryBuilderV2.h"
51#include "AliMUONSt1SpecialMotif.h"
52#include "AliMUON.h"
53#include "AliMUONChamber.h"
54#include "AliMUONChamberGeometry.h"
ba030c0e 55#include "AliRun.h"
56#include "AliMagF.h"
ba030c0e 57
d1cd2474 58ClassImp(AliMUONSt1GeometryBuilderV2)
ba030c0e 59
5f91c9e8 60// Thickness Constants
d1cd2474 61const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzPadPlane=0.0148/2.; //Pad plane
62const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzFoam = 2.083/2.; //Foam of mechanicalplane
63const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzFR4 = 0.0031/2.; //FR4 of mechanical plane
64const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzSnPb = 0.0091/2.; //Pad/Kapton connection (66 pt)
65const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzKapton = 0.0122/2.; //Kapton
66const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzBergPlastic = 0.3062/2.;//Berg connector
67const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzBergCopper = 0.1882/2.; //Berg connector
68const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzDaughter = 0.0156/2.; //Daughter board
69const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzGas = 0.2/2.; //Gas thickness
5f91c9e8 70
71// Quadrant Mother volume - TUBS1 - Middle layer of model
d1cd2474 72const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherIR1 = 18.3;
73const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherOR1 = 105.673;
74const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherThick1 = 6.5/2;
75const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherPhiL1 = 0.;
76const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherPhiU1 = 90.;
5f91c9e8 77
78// Quadrant Mother volume - TUBS2 - near and far layers of model
d1cd2474 79const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherIR2 = 20.7;
80const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherOR2 = 100.073;
81const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherThick2 = 3.0/2;
82const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherPhiL2 = 0.;
83const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherPhiU2 = 90.;
5f91c9e8 84
85// Sensitive copper pads, foam layer, PCB and electronics model parameters
d1cd2474 86const GReal_t AliMUONSt1GeometryBuilderV2::fgkHxHole=1.5/2.;
87const GReal_t AliMUONSt1GeometryBuilderV2::fgkHyHole=6./2.;
88const GReal_t AliMUONSt1GeometryBuilderV2::fgkHxBergPlastic=0.74/2.;
89const GReal_t AliMUONSt1GeometryBuilderV2::fgkHyBergPlastic=5.09/2.;
90const GReal_t AliMUONSt1GeometryBuilderV2::fgkHxBergCopper=0.25/2.;
91const GReal_t AliMUONSt1GeometryBuilderV2::fgkHyBergCopper=3.6/2.;
92const GReal_t AliMUONSt1GeometryBuilderV2::fgkHxKapton=0.8/2.;
93const GReal_t AliMUONSt1GeometryBuilderV2::fgkHyKapton=5.7/2.;
94const GReal_t AliMUONSt1GeometryBuilderV2::fgkHxDaughter=2.3/2.;
95const GReal_t AliMUONSt1GeometryBuilderV2::fgkHyDaughter=6.3/2.;
96const GReal_t AliMUONSt1GeometryBuilderV2::fgkOffsetX=1.46;
97const GReal_t AliMUONSt1GeometryBuilderV2::fgkOffsetY=0.71;
98const GReal_t AliMUONSt1GeometryBuilderV2::fgkDeltaFilleEtamX=1.46;
99const GReal_t AliMUONSt1GeometryBuilderV2::fgkDeltaFilleEtamY=0.051;
100
101const GReal_t AliMUONSt1GeometryBuilderV2::fgkDeltaQuadLHC=2.6; // LHC Origin wrt Quadrant Origin
102const GReal_t AliMUONSt1GeometryBuilderV2::fgkFrameOffset=5.0;
103
104const char* AliMUONSt1GeometryBuilderV2::fgkHoleName="MCHL";
105const char* AliMUONSt1GeometryBuilderV2::fgkDaughterName="MCDB";
106const char AliMUONSt1GeometryBuilderV2::fgkFoamLayerSuffix='F'; // prefix for automatic volume naming
107const char* AliMUONSt1GeometryBuilderV2::fgkQuadrantMLayerName="SQM";
108const char* AliMUONSt1GeometryBuilderV2::fgkQuadrantNLayerName="SQN";
109const char* AliMUONSt1GeometryBuilderV2::fgkQuadrantFLayerName="SQF";
5f91c9e8 110
111//______________________________________________________________________________
d1cd2474 112AliMUONSt1GeometryBuilderV2::AliMUONSt1GeometryBuilderV2(AliMUON* muon)
113 : AliMUONVGeometryBuilder(&muon->Chamber(0), &muon->Chamber(1)),
114 fMUON(muon)
ba030c0e 115{
5f91c9e8 116 // set path to mapping data files
117 if (! gSystem->Getenv("MINSTALL")) {
118 TString dirPath = gSystem->Getenv("ALICE_ROOT");
119 dirPath += "/MUON/mapping";
120 AliMpFiles::Instance()->SetTopPath(dirPath);
121 gSystem->Setenv("MINSTALL", dirPath.Data());
122 //cout << "AliMpFiles top path set to " << dirPath << endl;
123 }
124 //else
d1cd2474 125 // cout << gSystem->Getenv("MINSTALL") << endl;
ba030c0e 126}
127
5f91c9e8 128//______________________________________________________________________________
d1cd2474 129AliMUONSt1GeometryBuilderV2::AliMUONSt1GeometryBuilderV2()
130 : AliMUONVGeometryBuilder(),
131 fMUON(0)
ba030c0e 132{
d1cd2474 133// Default Constructor
134// --
ba030c0e 135}
136
5f91c9e8 137//______________________________________________________________________________
d1cd2474 138AliMUONSt1GeometryBuilderV2::AliMUONSt1GeometryBuilderV2(const AliMUONSt1GeometryBuilderV2& rhs)
139 : AliMUONVGeometryBuilder(rhs)
ba030c0e 140{
141// Dummy copy constructor
d1cd2474 142
143 Fatal("Copy constructor",
144 "Copy constructor is not implemented.");
ba030c0e 145}
146
5f91c9e8 147//______________________________________________________________________________
d1cd2474 148AliMUONSt1GeometryBuilderV2::~AliMUONSt1GeometryBuilderV2()
ba030c0e 149{
150// Destructor
ba030c0e 151}
152
d1cd2474 153
154//______________________________________________________________________________
155AliMUONSt1GeometryBuilderV2&
156AliMUONSt1GeometryBuilderV2::operator = (const AliMUONSt1GeometryBuilderV2& rhs)
157{
158 // check assignement to self
159 if (this == &rhs) return *this;
160
161 Fatal("operator=",
162 "Assignment operator is not implemented.");
163
164 return *this;
165}
166
5f91c9e8 167//
168// Private methods
169//
170
171//______________________________________________________________________________
d1cd2474 172void AliMUONSt1GeometryBuilderV2::CreateHole()
ba030c0e 173{
5f91c9e8 174// Create all the elements found inside a foam hole
ba030c0e 175// --
d1cd2474 176 Int_t* idtmed = fMUON->GetIdtmed()->GetArray()-1099;
177 Int_t idAir = idtmed[1100]; // medium 1
178 //Int_t idCopper = idtmed[1109]; // medium 10 = copper
179 Int_t idCopper = idtmed[1121]; // medium 22 = copper
ba030c0e 180
5f91c9e8 181 GReal_t par[3];
182 GReal_t posX,posY,posZ;
183
184 par[0] = fgkHxHole;
185 par[1] = fgkHyHole;
186 par[2] = fgkHzFoam;
187 gMC->Gsvolu(fgkHoleName,"BOX",idAir,par,3);
ba030c0e 188
5f91c9e8 189 par[0] = fgkHxKapton;
190 par[1] = fgkHyKapton;
191 par[2] = fgkHzSnPb;
192 gMC->Gsvolu("SNPB", "BOX", idCopper, par, 3);
193 posX = 0.;
194 posY = 0.;
195 posZ = -fgkHzFoam+fgkHzSnPb;
196 gMC->Gspos("SNPB",1,fgkHoleName, posX, posY, posZ, 0,"ONLY");
ba030c0e 197
5f91c9e8 198 par[0] = fgkHxHole;
199 par[1] = fgkHyBergPlastic;
200 par[2] = fgkHzKapton;
201 gMC->Gsvolu("KAPT", "BOX", idCopper, par, 3);
202 posX = 0.;
203 posY = 0.;
204 posZ = 0.;
205 gMC->Gspos("KAPT",1,fgkHoleName, posX, posY, posZ, 0,"ONLY");
206}
ba030c0e 207
5f91c9e8 208//______________________________________________________________________________
d1cd2474 209void AliMUONSt1GeometryBuilderV2::CreateDaughterBoard()
5f91c9e8 210{
211// Create all the elements in a daughter board
212// --
d1cd2474 213 Int_t* idtmed = fMUON->GetIdtmed()->GetArray()-1099;
5f91c9e8 214 Int_t idAir = idtmed[1100]; // medium 1
d1cd2474 215 //Int_t idCopper = idtmed[1109]; // medium 10 = copper
216 //Int_t idPlastic =idtmed[1116]; // medium 17 = Plastic
217 Int_t idCopper = idtmed[1121]; // medium 22 = copper
218 Int_t idPlastic =idtmed[1127]; // medium 28 = Plastic
ba030c0e 219
5f91c9e8 220 GReal_t par[3];
221 GReal_t posX,posY,posZ;
222
223 par[0]=fgkHxDaughter;
224 par[1]=fgkHyDaughter;
225 par[2]=TotalHzDaughter();
226 gMC->Gsvolu(fgkDaughterName,"BOX",idAir,par,3);
ba030c0e 227
5f91c9e8 228 par[0]=fgkHxBergPlastic;
229 par[1]=fgkHyBergPlastic;
230 par[2]=fgkHzBergPlastic;
231 gMC->Gsvolu("BRGP","BOX",idPlastic,par,3);
232 posX=0.;
233 posY=0.;
234 posZ = -TotalHzDaughter() + fgkHzBergPlastic;
235 gMC->Gspos("BRGP",1,fgkDaughterName,posX,posY,posZ,0,"ONLY");
236
237 par[0]=fgkHxBergCopper;
238 par[1]=fgkHyBergCopper;
239 par[2]=fgkHzBergCopper;
240 gMC->Gsvolu("BRGC","BOX",idCopper,par,3);
241 posX=0.;
242 posY=0.;
243 posZ=0.;
244 gMC->Gspos("BRGC",1,"BRGP",posX,posY,posZ,0,"ONLY");
245
246 par[0]=fgkHxDaughter;
247 par[1]=fgkHyDaughter;
248 par[2]=fgkHzDaughter;
249 gMC->Gsvolu("DGHT","BOX",idCopper,par,3);
250 posX=0.;
251 posY=0.;
252 posZ = -TotalHzDaughter() + 2.*fgkHzBergPlastic + fgkHzDaughter;
253 gMC->Gspos("DGHT",1,fgkDaughterName,posX,posY,posZ,0,"ONLY");
254}
255
256//______________________________________________________________________________
d1cd2474 257void AliMUONSt1GeometryBuilderV2::CreateInnerLayers()
5f91c9e8 258{
259// Create the layer of sensitive volumes with gas
260// and the copper layer.
261// --
262
263// Gas Medium
d1cd2474 264 Int_t* idtmed = fMUON->GetIdtmed()->GetArray()-1099;
265 //Int_t idArCO2 = idtmed[1108]; // medium 9 (ArCO2 80%)
266 //Int_t idCopper = idtmed[1109]; // medium 10 = copper
267 Int_t idArCO2 = idtmed[1124]; // medium 25 (ArCO2 80%)
268 Int_t idCopper = idtmed[1121]; // medium 22 = copper
5f91c9e8 269
270 Float_t par[11];
271
272//Make gas volume - composed of 11 trapezoids
273// section 1 of 11
274 par[0] = fgkHzGas;
275 par[1] = 0.;
276 par[2] = 0.;
277 par[3] = 71.33/2.;
278 par[4] = 9.76/2.;
279 par[5] = 48.77/2.;
280 par[6] = 15.3;
281 par[7] = 71.33/2.;
282 par[8] = 9.76/2.;
283 par[9] = 48.77/2.;
284 par[10] = 15.3;
285
286 gMC->Gsvolu("SA1G", "TRAP", idArCO2, par, 11);
287 gMC->Gsvolu("SA2G", "TRAP", idArCO2, par, 11);
ba030c0e 288
5f91c9e8 289 par[0] = fgkHzPadPlane;
290 gMC->Gsvolu("SA1C", "TRAP", idCopper,par, 11);
291
292// section 2 of 11
293 par[0] = fgkHzGas;
294 par[1] = 0.;
295 par[2] = 0.;
296 par[3] = 79.68/2.;
297 par[4] = 10.4/2.;
298 par[5] = 57.0/2.;
299 par[6] = 0.;
300 par[7] = 79.68/2.;
301 par[8] = 10.4/2.;
302 par[9] = 57.0/2.;
303 par[10] = 0.;
304 gMC->Gsvolu("SB1G", "TRAP", idArCO2, par, 11);
305 gMC->Gsvolu("SB2G", "TRAP", idArCO2, par, 11);
306
307 par[0] = fgkHzPadPlane;
308 gMC->Gsvolu("SB1C", "TRAP", idCopper,par, 11);
309
310// section 3 of 11
311 par[0] = fgkHzGas;
312 par[1] = 0.;
313 par[2] = 0.;
314 par[3] = 71.33/2.;
315 par[4] = 48.77/2.;
316 par[5] = 9.73/2.;
317 par[6] = -15.3;
318 par[7] = 71.33/2.;
319 par[8] = 48.77/2.;
320 par[9] = 9.73/2.;
321 par[10] = -15.3;
322
323 gMC->Gsvolu("SC1G", "TRAP", idArCO2, par, 11);
324 gMC->Gsvolu("SC2G", "TRAP", idArCO2, par, 11);
325
326 par[0] = fgkHzPadPlane;
327 gMC->Gsvolu("SC1C", "TRAP", idCopper,par, 11);
328
329// section 4 of 11
330 par[0] = fgkHzGas;
331 par[1] = 0.;
332 par[2] = 0.;
333 par[3] = 6.00/2.;
334 par[4] = 0.;
335 par[5] = 1.56/2.;
336 par[6] = 7.41;
337 par[7] = 6.00/2.;
338 par[8] = 0.;
339 par[9] = 1.56/2.;
340 par[10] = 7.41;
341 gMC->Gsvolu("SD1G", "TRAP", idArCO2, par, 11);
342 gMC->Gsvolu("SD2G", "TRAP", idArCO2, par, 11);
343
344 par[0] = fgkHzPadPlane;
345 gMC->Gsvolu("SD1C", "TRAP", idCopper,par, 11);
346
347// section 5 of 11
348 par[0] = fgkHzGas;
349 par[1] = 0.;
350 par[2] = 0.;
351 par[3] = 1.516/2.;
352 par[4] = 0.;
353 par[5] = 0.829/2.;
354 par[6] = 15.3;
355 par[7] = 1.516/2.;
356 par[8] = 0.;
357 par[9] = 0.829/2.;
358 par[10] = 15.3;
359 gMC->Gsvolu("SE1G", "TRAP", idArCO2, par, 11);
360 gMC->Gsvolu("SE2G", "TRAP", idArCO2, par, 11);
361
362 par[0] = fgkHzPadPlane;
363 gMC->Gsvolu("SE1C", "TRAP", idCopper,par, 11);
364
365// section 6 of 11
366 par[0] = fgkHzGas;
367 par[1] = 0.;
368 par[2] = 0.;
369 par[3] = 3.92/2.;
370 par[4] = 0.;
371 par[5] = 0.562/2.;
372 par[6] = -4.1;
373 par[7] = 3.92/2.;
374 par[8] = 0.;
375 par[9] = 0.562/2.;
376 par[10] = -4.1;
377 gMC->Gsvolu("SF1G", "TRAP", idArCO2, par, 11);
378 gMC->Gsvolu("SF2G", "TRAP", idArCO2, par, 11);
379
380 par[0] = fgkHzPadPlane;
381 gMC->Gsvolu("SF1C", "TRAP", idCopper,par, 11);
382
383// section 7 of 11
384 par[0] = fgkHzGas;
385 par[1] = 0.;
386 par[2] = 0.;
387 par[3] = 0.941/2.;
388 par[4] = 0.562/2.;
389 par[5] = 0.;
390 par[6] = -16.6;
391 par[7] = 0.941/2.;
392 par[8] = 0.562/2.;
393 par[9] = 0.;
394 par[10] =-16.6;
395 gMC->Gsvolu("SG1G", "TRAP", idArCO2, par, 11);
396 gMC->Gsvolu("SG2G", "TRAP", idArCO2, par, 11);
397
398 par[0] = fgkHzPadPlane;
399 gMC->Gsvolu("SG1C", "TRAP", idCopper,par, 11);
400
401// section 8 of 11
402 par[0] = fgkHzGas;
403 par[1] = 0.;
404 par[2] = 0.;
405 par[3] = 3.94/2.;
406 par[4] = 0.57/2.;
407 par[5] = 0.;
408 par[6] = 4.14;
409 par[7] = 3.94/2.;
410 par[8] = 0.57/2.;
411 par[9] = 0.;
412 par[10] = 4.14;
413 gMC->Gsvolu("SH1G", "TRAP", idArCO2, par, 11);
414 gMC->Gsvolu("SH2G", "TRAP", idArCO2, par, 11);
415
416 par[0] = fgkHzPadPlane;
417 gMC->Gsvolu("SH1C", "TRAP", idCopper,par, 11);
418
419// section 9 of 11
420 par[0] = fgkHzGas;
421 par[1] = 0.;
422 par[2] = 0.;
423 par[3] = 0.95/2.;
424 par[4] = 0.;
425 par[5] = 0.57/2;
426 par[6] = 16.7;
427 par[7] = 0.95/2.;
428 par[8] = 0.;
429 par[9] = 0.57/2;
430 par[10] = 16.7;
431 gMC->Gsvolu("SI1G", "TRAP", idArCO2, par, 11);
432 gMC->Gsvolu("SI2G", "TRAP", idArCO2, par, 11);
433
434 par[0] = fgkHzPadPlane;
435 gMC->Gsvolu("SI1C", "TRAP", idCopper,par, 11);
436
437// section 10 of 11
438 par[0] = fgkHzGas;
439 par[1] = 0.;
440 par[2] = 0.;
441 par[3] = 1.49/2.;
442 par[4] = 0.;
443 par[5] = 0.817/2.;
444 par[6] = -15.4;
445 par[7] = 1.49/2.;
446 par[8] = 0.;
447 par[9] = 0.817/2.;
448 par[10] = -15.4;
449 gMC->Gsvolu("SJ1G", "TRAP", idArCO2, par, 11);
450 gMC->Gsvolu("SJ2G", "TRAP", idArCO2, par, 11);
451
452 par[0] = fgkHzPadPlane;
453 gMC->Gsvolu("SJ1C", "TRAP", idCopper,par, 11);
454
455// section 11 of 11
456 par[0] = fgkHzGas;
457 par[1] = 0.;
458 par[2] = 0.;
459 par[3] = 5.93/2.;
460 par[4] = 0.;
461 par[5] = 1.49/2.;
462 par[6] = -7.16;
463 par[7] = 5.93/2.;
464 par[8] = 0.;
465 par[9] = 1.49/2.;
466 par[10] = -7.16;
467 gMC->Gsvolu("SK1G", "TRAP", idArCO2, par, 11);
468 gMC->Gsvolu("SK2G", "TRAP", idArCO2, par, 11);
469
470 par[0] = fgkHzPadPlane;
471 gMC->Gsvolu("SK1C", "TRAP", idCopper,par, 11);
ba030c0e 472}
5f91c9e8 473
474//______________________________________________________________________________
d1cd2474 475void AliMUONSt1GeometryBuilderV2::CreateQuadrant(Int_t chamber)
ba030c0e 476{
477// create the quadrant (bending and non-bending planes)
478// for the given chamber
479// --
5f91c9e8 480
ba030c0e 481 CreateFrame(chamber);
482
d1cd2474 483#ifdef ST1_WITH_STL
484 SpecialMap specialMap;
5f91c9e8 485 specialMap[1001] = AliMUONSt1SpecialMotif(TVector2( 0.1, 0.84), 90.);
486 specialMap[1002] = AliMUONSt1SpecialMotif(TVector2( 0.5, 0.36));
487 specialMap[1003] = AliMUONSt1SpecialMotif(TVector2(1.01, 0.36));
d1cd2474 488#endif
489
490#ifdef ST1_WITH_ROOT
491 SpecialMap specialMap;
492 specialMap.Add(1001, (Long_t) new AliMUONSt1SpecialMotif(TVector2( 0.1, 0.84), 90.));
493 specialMap.Add(1002, (Long_t) new AliMUONSt1SpecialMotif(TVector2( 0.5, 0.36)));
494 specialMap.Add(1003, (Long_t) new AliMUONSt1SpecialMotif(TVector2(1.01, 0.36)));
495#endif
496
4a147321 497 AliMpReader reader1(kStation1, kBendingPlane);
5f91c9e8 498 AliMpSector* sector1 = reader1.BuildSector();
ba030c0e 499
5f91c9e8 500 Bool_t reflectZ = true;
501 TVector3 where = TVector3(2.5+0.1+0.56+0.001, 2.5+0.1+0.001, 0.);
502 PlaceSector(sector1, specialMap, where, reflectZ, chamber);
ba030c0e 503
d1cd2474 504#ifdef ST1_WITH_STL
ba030c0e 505 specialMap.clear();
506 specialMap[4001] = AliMUONSt1SpecialMotif(TVector2(1.01,0.59),90.);
507 specialMap[4002] = AliMUONSt1SpecialMotif(TVector2(1.96, 0.17));
508 specialMap[4003] = AliMUONSt1SpecialMotif(TVector2(1.61,-1.18));
509 specialMap[4004] = AliMUONSt1SpecialMotif(TVector2(0.2 ,-0.08));
510 specialMap[4005] = AliMUONSt1SpecialMotif(TVector2(0.2 , 0.25));
511 specialMap[4006] = AliMUONSt1SpecialMotif(TVector2(0.28, 0.21));
d1cd2474 512#endif
513
514#ifdef ST1_WITH_ROOT
515 specialMap.Delete();
516 specialMap.Add(4001,(Long_t) new AliMUONSt1SpecialMotif(TVector2(1.01,0.59),90.));
517 specialMap.Add(4002,(Long_t) new AliMUONSt1SpecialMotif(TVector2(1.96, 0.17)));
518 specialMap.Add(4003,(Long_t) new AliMUONSt1SpecialMotif(TVector2(1.61,-1.18)));
519 specialMap.Add(4004,(Long_t) new AliMUONSt1SpecialMotif(TVector2(0.2 ,-0.08)));
520 specialMap.Add(4005,(Long_t) new AliMUONSt1SpecialMotif(TVector2(0.2 , 0.25)));
521 specialMap.Add(4006,(Long_t) new AliMUONSt1SpecialMotif(TVector2(0.28, 0.21)));
522#endif
523
4a147321 524 AliMpReader reader2(kStation1, kNonBendingPlane);
5f91c9e8 525 AliMpSector* sector2 = reader2.BuildSector();
ba030c0e 526
5f91c9e8 527 reflectZ = false;
528 where = TVector3(where.X()+0.63/2.,where.Y()+0.42/2., 0.); //add a half pad shift
529 PlaceSector(sector2, specialMap, where, reflectZ, chamber);
d1cd2474 530
531#ifdef ST1_WITH_ROOT
532 specialMap.Delete();
533#endif
ba030c0e 534}
535
5f91c9e8 536//______________________________________________________________________________
d1cd2474 537void AliMUONSt1GeometryBuilderV2::CreateFoamBox(const char* name,const TVector2& dimensions)
ba030c0e 538{
5f91c9e8 539// create all the elements in the copper plane
540// --
ba030c0e 541
d1cd2474 542 Int_t* idtmed = fMUON->GetIdtmed()->GetArray()-1099;
5f91c9e8 543 Int_t idAir = idtmed[1100]; // medium 1
d1cd2474 544 //Int_t idFoam = idtmed[1115]; // medium 16 = Foam
545 //Int_t idFR4 = idtmed[1114]; // medium 15 = FR4
546 Int_t idFoam = idtmed[1125]; // medium 26 = Foam
547 Int_t idFR4 = idtmed[1122]; // medium 23 = FR4
ba030c0e 548
5f91c9e8 549 // mother volume
550 GReal_t par[3];
551 par[0] = dimensions.X();
552 par[1] = dimensions.Y();
553 par[2] = TotalHzPlane();
554 gMC->Gsvolu(name,"BOX",idAir,par,3);
555
556 // foam layer
557 GReal_t posX,posY,posZ;
558 char eName[5];
559 strcpy(eName,name);
560 eName[3]=fgkFoamLayerSuffix;
561 par[0] = dimensions.X();
562 par[1] = dimensions.Y();
563 par[2] = fgkHzFoam;
564 gMC->Gsvolu(eName,"BOX",idFoam,par,3);
565 posX=0.;
566 posY=0.;
567 posZ = -TotalHzPlane() + fgkHzFoam;
568 gMC->Gspos(eName,1,name,posX,posY,posZ,0,"ONLY");
ba030c0e 569
5f91c9e8 570 // mechanical plane FR4 layer
571 eName[3]='R';
572 par[0] = dimensions.X();
573 par[1] = dimensions.Y();
574 par[2] = fgkHzFR4;
575 gMC->Gsvolu(eName,"BOX",idFR4,par,3);
576 posX=0.;
577 posY=0.;
578 posZ = -TotalHzPlane()+ 2.*fgkHzFoam + fgkHzFR4;
579 gMC->Gspos(eName,1,name,posX,posY,posZ,0,"ONLY");
580}
ba030c0e 581
5f91c9e8 582//______________________________________________________________________________
d1cd2474 583void AliMUONSt1GeometryBuilderV2::CreatePlaneSegment(const char* name,const TVector2& dimensions,
5f91c9e8 584 Int_t nofHoles)
585{
586// Create a segment of a plane (this includes a foam layer,
587// holes in the foam to feed the kaptons through, kapton connectors
588// and the mother board.)
589// --
ba030c0e 590
5f91c9e8 591 CreateFoamBox(name,dimensions);
ba030c0e 592
5f91c9e8 593 char eName[5];
594 strcpy(eName,name);
595 eName[3]=fgkFoamLayerSuffix;
ba030c0e 596
5f91c9e8 597 for (Int_t holeNum=0;holeNum<nofHoles;holeNum++) {
598 GReal_t posX = ((2.*holeNum+1.)/nofHoles-1.)*dimensions.X();
599 GReal_t posY = 0.;
600 GReal_t posZ = 0.;
ba030c0e 601
5f91c9e8 602 gMC->Gspos(fgkHoleName,holeNum+1,eName,posX,posY,posZ,0,"ONLY");
603 }
604}
ba030c0e 605
5f91c9e8 606//______________________________________________________________________________
d1cd2474 607void AliMUONSt1GeometryBuilderV2::CreateFrame(Int_t chamber)
5f91c9e8 608{
609// Create the non-sensitive elements of the frame for the <chamber>
610//
611//
612// Model and notation:
613//
614// The Quadrant volume name starts with SQ
615// The volume segments are numbered 00 to XX.
616//
617// OutTopFrame
618// (SQ02-16)
619// ------------
620// OutEdgeFrame / |
621// (SQ17-24) / | InVFrame (SQ00-01)
622// / |
623// | |
624// OutVFrame | _- -
625// (SQ25-39) | | InArcFrame (SQ42-45)
626// | |
627// -------------
628// InHFrame (SQ40-41)
629//
630//
631// 06 February 2003 - Overlapping volumes resolved.
632// One quarter chamber is comprised of three TUBS volumes: SQMx, SQNx, and SQFx,
633// where SQMx is the Quadrant Middle layer for chamber <x> ( posZ in [-3.25,3.25]),
634// SQNx is the Quadrant Near side layer for chamber <x> ( posZ in [-6.25,3-.25) ), and
635// SQFx is the Quadrant Far side layer for chamber <x> ( posZ in (3.25,6.25] ).
636//---
ba030c0e 637
d1cd2474 638 const Float_t kNearFarLHC=2.4; // Near and Far TUBS Origin wrt LHC Origin
ba030c0e 639
ba030c0e 640 // tracking medias
d1cd2474 641 Int_t* idtmed = fMUON->GetIdtmed()->GetArray()-1099;
5f91c9e8 642
ba030c0e 643 Int_t idAir = idtmed[1100]; // medium 1
d1cd2474 644 //Int_t idFrameEpoxy = idtmed[1115]; // medium 16 = Frame Epoxy ME730
645 //Int_t idInox = idtmed[1116]; // medium 17 Stainless Steel (18%Cr,9%Ni,Fe)
646 //Int_t idFR4 = idtmed[1110]; // medium 11 FR4
647 //Int_t idCopper = idtmed[1109]; // medium 10 Copper
648 //Int_t idAlu = idtmed[1103]; // medium 4 Aluminium
649 Int_t idFrameEpoxy = idtmed[1123]; // medium 24 = Frame Epoxy ME730 // was 20 not 16
650 Int_t idInox = idtmed[1128]; // medium 29 Stainless Steel (18%Cr,9%Ni,Fe) // was 21 not 17
651 Int_t idFR4 = idtmed[1122]; // medium 23 FR4 // was 15 not 11
652 Int_t idCopper = idtmed[1121]; // medium 22 Copper
653 Int_t idAlu = idtmed[1120]; // medium 21 Aluminium
5f91c9e8 654
655
656// Rotation Matrices
657 Int_t rot1, rot2, rot3;
658
659// Rotation matrices
d1cd2474 660 fMUON->AliMatrix(rot1, 90., 90., 90., 180., 0., 0.); // +90 deg in x-y plane
661 fMUON->AliMatrix(rot2, 90., 45., 90., 135., 0., 0.); // +45 deg in x-y plane
662