Removed all geometry & segmentation info:
[u/mrichter/AliRoot.git] / MUON / AliMUONSt1GeometryBuilder.cxx
CommitLineData
e118b27e 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//
18// Class AliMUONSt1GeometryBuilder
19// -------------------------------
20// MUON Station1 coarse geometry construction class.
d1cd2474 21// Extracted from AliMUONv1
22// by Ivana Hrivnacova, IPN Orsay
5f1df83a 23// Included in AliRoot 2004/01/23
d1cd2474 24
25#include <TVirtualMC.h>
26#include <TGeoMatrix.h>
27
e118b27e 28#include "AliLog.h"
29
d1cd2474 30#include "AliMUONSt1GeometryBuilder.h"
31#include "AliMUON.h"
32#include "AliMUONChamber.h"
e118b27e 33#include "AliMUONGeometryModule.h"
a432117a 34#include "AliMUONGeometryEnvelopeStore.h"
d1cd2474 35
36ClassImp(AliMUONSt1GeometryBuilder)
37
38//______________________________________________________________________________
39AliMUONSt1GeometryBuilder::AliMUONSt1GeometryBuilder(AliMUON* muon)
a432117a 40 : AliMUONVGeometryBuilder("st1.dat",
e118b27e 41 muon->Chamber(0).GetGeometry(),
42 muon->Chamber(1).GetGeometry()),
d1cd2474 43 fMUON(muon)
44{
45// Standard constructor
46
47}
48
49//______________________________________________________________________________
50AliMUONSt1GeometryBuilder::AliMUONSt1GeometryBuilder()
51 : AliMUONVGeometryBuilder(),
52 fMUON(0)
53{
54// Default constructor
55}
56
57
58//______________________________________________________________________________
59AliMUONSt1GeometryBuilder::AliMUONSt1GeometryBuilder(const AliMUONSt1GeometryBuilder& rhs)
60 : AliMUONVGeometryBuilder(rhs)
61{
30178c30 62// Protected copy constructor
63
8c343c7c 64 AliFatal("Copy constructor is not implemented.");
d1cd2474 65}
66
67//______________________________________________________________________________
68AliMUONSt1GeometryBuilder::~AliMUONSt1GeometryBuilder() {
69//
70}
71
72//______________________________________________________________________________
73AliMUONSt1GeometryBuilder&
74AliMUONSt1GeometryBuilder::operator = (const AliMUONSt1GeometryBuilder& rhs)
75{
30178c30 76// Protected assignement operator
77
d1cd2474 78 // check assignement to self
79 if (this == &rhs) return *this;
80
8c343c7c 81 AliFatal("Assignment operator is not implemented.");
d1cd2474 82
83 return *this;
84}
85
86//
87// public methods
88//
89
90//______________________________________________________________________________
91void AliMUONSt1GeometryBuilder::CreateGeometry()
92{
93// From AliMUONv1::CreateGeometry()
94// ---
95
96//********************************************************************
97// Station 1 **
98//********************************************************************
99// CONCENTRIC
100 // indices 1 and 2 for first and second chambers in the station
101 // iChamber (first chamber) kept for other quanties than Z,
102 // assumed to be the same in both chambers
103
104 // Get tracking medias Ids
105 Int_t *idtmed = fMUON->GetIdtmed()->GetArray()-1099;
106 Int_t idAir= idtmed[1100]; // medium 1
107 Int_t idAlu1=idtmed[1103]; // medium 4
108 Int_t idAlu2=idtmed[1104]; // medium 5
109 Int_t idGas=idtmed[1108]; // medium 9 = Ar-CO2 gas (80%+20%)
110 Bool_t frameCrosses=kTRUE;
111
112 // Rotation matrices in the x-y plane
113 // phi= 0 deg
114 Int_t irot1;
115 fMUON->AliMatrix(irot1, 90., 0., 90., 90., 0., 0.);
116 // phi= 90 deg
117 Int_t irot2;
118 fMUON->AliMatrix(irot2, 90., 90., 90., 180., 0., 0.);
119
e118b27e 120 AliMUONChamber* iChamber1 = &fMUON->Chamber(0);
121 AliMUONChamber* iChamber2 = &fMUON->Chamber(1);
d1cd2474 122 AliMUONChamber* iChamber = iChamber1;
123
124 // DGas decreased from standard one (0.5)
125 iChamber->SetDGas(0.4);
126 iChamber2->SetDGas(0.4);
127
128 // DAlu increased from standard one (3% of X0),
129 // because more electronics with smaller pads
130 iChamber->SetDAlu(3.5 * 8.9 / 100.);
131 iChamber2->SetDAlu(3.5 * 8.9 / 100.);
132
133 // Half of the total thickness of frame crosses (including DAlu)
134 // for each chamber in stations 1 and 2:
135 // 3% of X0 of composite material,
136 // but taken as Aluminium here, with same thickness in number of X0
137 Float_t dframez = 3. * 8.9 / 100;
138 Float_t zfpos=-(iChamber->DGas()+dframez+iChamber->DAlu())/2;
139 // The same parameters are defined in builder for station 2
140
141 // Mother volume
142 // Outer excess and inner recess for mother volume radius
143 // with respect to ROuter and RInner
144 Float_t dframep=.001; // Value for station 3 should be 6 ...
145 // Width (RdPhi) of the frame crosses for stations 1 and 2 (cm)
146 // Float_t dframep1=.001;
147 Float_t dframep1 = 11.0;
148 Float_t phi=2*TMath::Pi()/12/2;
149 // The same parameters are defined in builder for station 2
150
151 Float_t tpar[3];
152 Float_t dstation = (-iChamber2->Z()) - (-iChamber1->Z());
153 tpar[0] = iChamber->RInner()-dframep;
154 tpar[1] = (iChamber->ROuter()+dframep)/TMath::Cos(phi);
155 tpar[2] = dstation/5;
156
157 gMC->Gsvolu("S01M", "TUBE", idAir, tpar, 3);
158 gMC->Gsvolu("S02M", "TUBE", idAir, tpar, 3);
159
160 // CHANGED
161 //gMC->Gspos("S01M", 1, "ALIC", 0., 0., zpos1 , 0, "ONLY");
a432117a 162 //gMC->Gspos("S02M", 1, "ALIC", 0., 0., zpos2 , 0, "ONLY");
163
164 GetEnvelopes(0)->AddEnvelope("S01M", 100, false);
165 GetEnvelopes(1)->AddEnvelope("S02M", 200, false);
166
d1cd2474 167
168// // Aluminium frames
169// // Outer frames
170// pgpar[0] = 360/12/2;
171// pgpar[1] = 360.;
172// pgpar[2] = 12.;
173// pgpar[3] = 2;
174// pgpar[4] = -dframez/2;
175// pgpar[5] = iChamber->ROuter();
176// pgpar[6] = pgpar[5]+dframep1;
177// pgpar[7] = +dframez/2;
178// pgpar[8] = pgpar[5];
179// pgpar[9] = pgpar[6];
180// gMC->Gsvolu("S01O", "PGON", idAlu1, pgpar, 10);
181// gMC->Gsvolu("S02O", "PGON", idAlu1, pgpar, 10);
182// gMC->Gspos("S01O",1,"S01M", 0.,0.,-zfpos, 0,"ONLY");
183// gMC->Gspos("S01O",2,"S01M", 0.,0.,+zfpos, 0,"ONLY");
184// gMC->Gspos("S02O",1,"S02M", 0.,0.,-zfpos, 0,"ONLY");
185// gMC->Gspos("S02O",2,"S02M", 0.,0.,+zfpos, 0,"ONLY");
186// //
187// // Inner frame
188// tpar[0]= iChamber->RInner()-dframep1;
189// tpar[1]= iChamber->RInner();
190// tpar[2]= dframez/2;
191// gMC->Gsvolu("S01I", "TUBE", idAlu1, tpar, 3);
192// gMC->Gsvolu("S02I", "TUBE", idAlu1, tpar, 3);
193
194// gMC->Gspos("S01I",1,"S01M", 0.,0.,-zfpos, 0,"ONLY");
195// gMC->Gspos("S01I",2,"S01M", 0.,0.,+zfpos, 0,"ONLY");
196// gMC->Gspos("S02I",1,"S02M", 0.,0.,-zfpos, 0,"ONLY");
197// gMC->Gspos("S02I",2,"S02M", 0.,0.,+zfpos, 0,"ONLY");
198//
199// Frame Crosses
200 if (frameCrosses) {
201 // outside gas
202 // security for inside mother volume
203 Float_t bpar[3];
204 bpar[0] = (iChamber->ROuter() - iChamber->RInner())
205 * TMath::Cos(TMath::ASin(dframep1 /
206 (iChamber->ROuter() - iChamber->RInner())))
207 / 2.0;
208 bpar[1] = dframep1/2;
209 // total thickness will be (4 * bpar[2]) for each chamber,
210 // which has to be equal to (2 * dframez) - DAlu
211 bpar[2] = (2.0 * dframez - iChamber->DAlu()) / 4.0;
212 gMC->Gsvolu("S01B", "BOX", idAlu1, bpar, 3);
213 gMC->Gsvolu("S02B", "BOX", idAlu1, bpar, 3);
214
215 gMC->Gspos("S01B",1,"S01M", +iChamber->RInner()+bpar[0] , 0,-zfpos,
216 irot1,"ONLY");
217 gMC->Gspos("S01B",2,"S01M", -iChamber->RInner()-bpar[0] , 0,-zfpos,
218 irot1,"ONLY");
219 gMC->Gspos("S01B",3,"S01M", 0, +iChamber->RInner()+bpar[0] ,-zfpos,
220 irot2,"ONLY");
221 gMC->Gspos("S01B",4,"S01M", 0, -iChamber->RInner()-bpar[0] ,-zfpos,
222 irot2,"ONLY");
223 gMC->Gspos("S01B",5,"S01M", +iChamber->RInner()+bpar[0] , 0,+zfpos,
224 irot1,"ONLY");
225 gMC->Gspos("S01B",6,"S01M", -iChamber->RInner()-bpar[0] , 0,+zfpos,
226 irot1,"ONLY");
227 gMC->Gspos("S01B",7,"S01M", 0, +iChamber->RInner()+bpar[0] ,+zfpos,
228 irot2,"ONLY");
229 gMC->Gspos("S01B",8,"S01M", 0, -iChamber->RInner()-bpar[0] ,+zfpos,
230 irot2,"ONLY");
231
232 gMC->Gspos("S02B",1,"S02M", +iChamber->RInner()+bpar[0] , 0,-zfpos,
233 irot1,"ONLY");
234 gMC->Gspos("S02B",2,"S02M", -iChamber->RInner()-bpar[0] , 0,-zfpos,
235 irot1,"ONLY");
236 gMC->Gspos("S02B",3,"S02M", 0, +iChamber->RInner()+bpar[0] ,-zfpos,
237 irot2,"ONLY");
238 gMC->Gspos("S02B",4,"S02M", 0, -iChamber->RInner()-bpar[0] ,-zfpos,
239 irot2,"ONLY");
240 gMC->Gspos("S02B",5,"S02M", +iChamber->RInner()+bpar[0] , 0,+zfpos,
241 irot1,"ONLY");
242 gMC->Gspos("S02B",6,"S02M", -iChamber->RInner()-bpar[0] , 0,+zfpos,
243 irot1,"ONLY");
244 gMC->Gspos("S02B",7,"S02M", 0, +iChamber->RInner()+bpar[0] ,+zfpos,
245 irot2,"ONLY");
246 gMC->Gspos("S02B",8,"S02M", 0, -iChamber->RInner()-bpar[0] ,+zfpos,
247 irot2,"ONLY");
248 }
249//
250// Chamber Material represented by Alu sheet
251 tpar[0]= iChamber->RInner();
252 tpar[1]= iChamber->ROuter();
253 tpar[2] = (iChamber->DGas()+iChamber->DAlu())/2;
254 gMC->Gsvolu("S01A", "TUBE", idAlu2, tpar, 3);
255 gMC->Gsvolu("S02A", "TUBE",idAlu2, tpar, 3);
256 gMC->Gspos("S01A", 1, "S01M", 0., 0., 0., 0, "ONLY");
257 gMC->Gspos("S02A", 1, "S02M", 0., 0., 0., 0, "ONLY");
258//
259// Sensitive volumes
260 // tpar[2] = iChamber->DGas();
261 tpar[2] = iChamber->DGas()/2;
262 gMC->Gsvolu("S01G", "TUBE", idGas, tpar, 3);
263 gMC->Gsvolu("S02G", "TUBE", idGas, tpar, 3);
264 gMC->Gspos("S01G", 1, "S01A", 0., 0., 0., 0, "ONLY");
265 gMC->Gspos("S02G", 1, "S02A", 0., 0., 0., 0, "ONLY");
266//
267// Frame Crosses to be placed inside gas
268 // NONE: chambers are sensitive everywhere
269// if (frameCrosses) {
270
271// dr = (iChamber->ROuter() - iChamber->RInner());
272// bpar[0] = TMath::Sqrt(dr*dr-dframep1*dframep1/4)/2;
273// bpar[1] = dframep1/2;
274// bpar[2] = iChamber->DGas()/2;
275// gMC->Gsvolu("S01F", "BOX", idAlu1, bpar, 3);
276// gMC->Gsvolu("S02F", "BOX", idAlu1, bpar, 3);
277
278// gMC->Gspos("S01F",1,"S01G", +iChamber->RInner()+bpar[0] , 0, 0,
279// irot1,"ONLY");
280// gMC->Gspos("S01F",2,"S01G", -iChamber->RInner()-bpar[0] , 0, 0,
281// irot1,"ONLY");
282// gMC->Gspos("S01F",3,"S01G", 0, +iChamber->RInner()+bpar[0] , 0,
283// irot2,"ONLY");
284// gMC->Gspos("S01F",4,"S01G", 0, -iChamber->RInner()-bpar[0] , 0,
285// irot2,"ONLY");
286
287// gMC->Gspos("S02F",1,"S02G", +iChamber->RInner()+bpar[0] , 0, 0,
288// irot1,"ONLY");
289// gMC->Gspos("S02F",2,"S02G", -iChamber->RInner()-bpar[0] , 0, 0,
290// irot1,"ONLY");
291// gMC->Gspos("S02F",3,"S02G", 0, +iChamber->RInner()+bpar[0] , 0,
292// irot2,"ONLY");
293// gMC->Gspos("S02F",4,"S02G", 0, -iChamber->RInner()-bpar[0] , 0,
294// irot2,"ONLY");
295// }
296}
297
298//______________________________________________________________________________
299void AliMUONSt1GeometryBuilder::SetTransformations()
300{
301// Defines the transformations for the station2 chambers.
302// ---
303
e118b27e 304 AliMUONChamber* iChamber1 = &fMUON->Chamber(0);
d1cd2474 305 Double_t zpos1= - iChamber1->Z();
306 iChamber1->GetGeometry()
307 ->SetTranslation(TGeoTranslation(0., 0., zpos1));
308
e118b27e 309 AliMUONChamber* iChamber2 = &fMUON->Chamber(1);
d1cd2474 310 Double_t zpos2 = - iChamber2->Z();
311 iChamber2->GetGeometry()
312 ->SetTranslation(TGeoTranslation(0., 0., zpos2));
313}
314
315//______________________________________________________________________________
316void AliMUONSt1GeometryBuilder::SetSensitiveVolumes()
317{
318// Defines the sensitive volumes for station1 chambers.
319// ---
320
e118b27e 321 GetGeometry(0)->SetSensitiveVolume("S01G");
322 GetGeometry(1)->SetSensitiveVolume("S02G");
d1cd2474 323}