Removing warnings ...
[u/mrichter/AliRoot.git] / MUON / AliMUONv2.cxx
Content-type: text/html git.uio.no Git - u/mrichter/AliRoot.git/blame - MUON/AliMUONv2.cxx


500 - Internal Server Error

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