]>
Commit | Line | Data |
---|---|---|
30178c30 | 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 AliMUONTriggerGeometryBuilder | |
19 | // ----------------------------------- | |
20 | // MUON Trigger stations geometry construction class. | |
21 | // | |
22 | // Author: Philippe Crochette, LPC Clermont-Ferrand | |
23 | ||
24 | #include <TVirtualMC.h> | |
25 | #include <TGeoMatrix.h> | |
26 | ||
27 | #include "AliMUONTriggerGeometryBuilder.h" | |
28 | #include "AliMUON.h" | |
29 | #include "AliMUONChamber.h" | |
30 | #include "AliMUONChamberGeometry.h" | |
31 | ||
32 | ClassImp(AliMUONTriggerGeometryBuilder) | |
33 | ||
34 | //______________________________________________________________________________ | |
35 | AliMUONTriggerGeometryBuilder::AliMUONTriggerGeometryBuilder(AliMUON* muon) | |
36 | : AliMUONVGeometryBuilder(&muon->Chamber(10), &muon->Chamber(11),&muon->Chamber(12),&muon->Chamber(13)), | |
37 | fMUON(muon) | |
38 | { | |
39 | // Standard constructor | |
40 | ||
41 | } | |
42 | ||
43 | //______________________________________________________________________________ | |
44 | AliMUONTriggerGeometryBuilder::AliMUONTriggerGeometryBuilder() | |
45 | : AliMUONVGeometryBuilder(), | |
46 | fMUON(0) | |
47 | { | |
48 | // Default constructor | |
49 | } | |
50 | ||
51 | ||
52 | //______________________________________________________________________________ | |
53 | AliMUONTriggerGeometryBuilder::AliMUONTriggerGeometryBuilder(const AliMUONTriggerGeometryBuilder& rhs) | |
54 | : AliMUONVGeometryBuilder(rhs) | |
55 | { | |
30178c30 | 56 | // Protected copy constructor |
57 | ||
d1cd2474 | 58 | Fatal("Copy constructor", |
59 | "Copy constructor is not implemented."); | |
60 | } | |
61 | ||
62 | //______________________________________________________________________________ | |
63 | AliMUONTriggerGeometryBuilder::~AliMUONTriggerGeometryBuilder() { | |
64 | // | |
65 | } | |
66 | ||
67 | //______________________________________________________________________________ | |
68 | AliMUONTriggerGeometryBuilder& | |
69 | AliMUONTriggerGeometryBuilder::operator = (const AliMUONTriggerGeometryBuilder& rhs) | |
70 | { | |
30178c30 | 71 | // Protected assignement operator |
72 | ||
d1cd2474 | 73 | // check assignement to self |
74 | if (this == &rhs) return *this; | |
75 | ||
76 | Fatal("operator=", | |
77 | "Assignment operator is not implemented."); | |
78 | ||
79 | return *this; | |
80 | } | |
81 | ||
82 | // | |
83 | // public methods | |
84 | // | |
85 | ||
86 | //______________________________________________________________________________ | |
87 | void AliMUONTriggerGeometryBuilder::CreateGeometry() | |
88 | { | |
89 | // From AliMUONv1::CreateGeometry() | |
90 | ||
91 | /* | |
92 | zpos1 and zpos2 are the middle of the first and second | |
93 | planes of station 1 (+1m for second station): | |
94 | zpos1=(zpos1m+zpos1p)/2=(15999+16071)/2=16035 mm, thick/2=40 mm | |
95 | zpos2=(zpos2m+zpos2p)/2=(16169+16241)/2=16205 mm, thick/2=40 mm | |
96 | zposxm and zposxp= middles of gaz gaps within a detection plane | |
97 | rem: the total thickness accounts for 1 mm of al on both | |
98 | side of the RPCs (see zpos1 and zpos2) | |
99 | */ | |
100 | ||
101 | Int_t *idtmed = fMUON->GetIdtmed()->GetArray()-1099; | |
102 | Int_t idAir= idtmed[1100]; // medium 1 | |
103 | Int_t idAlu1=idtmed[1103]; // medium 4 | |
104 | Float_t tpar[3]; | |
105 | Double_t dpar[3]; | |
106 | ||
107 | // vertical gap between right and left chambers (kDXZERO*2=4cm) | |
108 | const Float_t kDXZERO=2.; | |
109 | // main distances for chamber definition in first plane/first station | |
110 | const Float_t kXMIN=34.; | |
111 | const Float_t kXMED=51.; | |
112 | const Float_t kXMAX=272.; | |
113 | // kXMAX will become 255. in real life. segmentation to be updated accordingly | |
114 | // (see fig.2-4 & 2-5 of Local Trigger Board PRR) | |
115 | const Float_t kYMIN=34.; | |
116 | const Float_t kYMAX=51.; | |
117 | // inner/outer radius of flange between beam shield. and chambers (1/station) | |
118 | const Float_t kRMIN[2]={50.,50.}; | |
119 | const Float_t kRMAX[2]={64.,68.}; | |
120 | // z position of the middle of the gas gap in mother vol | |
121 | const Float_t kZm=-3.6; | |
122 | const Float_t kZp=+3.6; | |
123 | ||
124 | AliMUONChamber *iChamber, *iChamber1; | |
125 | iChamber1 = GetChamber(10); | |
126 | Float_t zpos1=-iChamber1->Z(); | |
127 | ||
128 | // ratio of zpos1m/zpos1p and inverse for first plane | |
129 | Float_t zmp=(zpos1-3.6)/(zpos1+3.6); | |
130 | Float_t zpm=1./zmp; | |
131 | ||
132 | Int_t icount=0; // chamber counter (0 1 2 3) | |
133 | ||
134 | for (Int_t istation=0; istation<2; istation++) { // loop on stations | |
135 | for (Int_t iplane=0; iplane<2; iplane++) { // loop on detection planes | |
136 | ||
137 | Int_t iVolNum=1; // counter Volume Number | |
138 | icount = Int_t(iplane*TMath::Power(2,0))+ | |
139 | Int_t(istation*TMath::Power(2,1)); | |
140 | ||
141 | iChamber = GetChamber(10+icount); | |
142 | Float_t zpos = - iChamber->Z(); | |
143 | ||
144 | // Flange between beam shielding and RPC | |
145 | tpar[0]= kRMIN[istation]; | |
146 | tpar[1]= kRMAX[istation]; | |
147 | tpar[2]= 4.0; | |
148 | char volFlange[5]; | |
149 | sprintf(volFlange,"SF%dA",icount+1); | |
150 | gMC->Gsvolu(volFlange,"TUBE",idAlu1,tpar,3); // Al | |
151 | // changed by ivana | |
152 | //gMC->Gspos(volFlange,1,"ALIC",0.,0.,zpos,0,"MANY"); | |
6b82c1f0 | 153 | iChamber->GetGeometry()->AddEnvelope(volFlange, false, "MANY"); |
d1cd2474 | 154 | |
155 | // scaling factor | |
156 | Float_t zRatio = zpos / zpos1; | |
157 | ||
158 | // envelopes (same size except line 5, all virtual) | |
159 | char volEnv[18][5]; | |
160 | tpar[1] = kYMIN * zRatio; | |
161 | tpar[2] = 0.4; | |
162 | Int_t i=0; // counter | |
163 | for (Int_t icolumn=0; icolumn<2; icolumn++) { | |
164 | for (Int_t iline=1; iline<10; iline++){ | |
165 | tpar[0] = (kXMAX/2.) * zRatio; | |
166 | if (iline==5) tpar[0] = ((kXMAX-kXMED)/2.)*zRatio; | |
167 | if (icolumn==0) | |
168 | sprintf(volEnv[i],"S%dL%d",icount,iline); | |
169 | else | |
170 | sprintf(volEnv[i],"S%dR%d",icount,iline); | |
171 | gMC->Gsvolu(volEnv[i],"BOX",idAir,tpar,0); | |
172 | i++; | |
173 | } | |
174 | } | |
175 | ||
176 | // chamber prototype | |
177 | tpar[0]= 0.; | |
178 | tpar[1]= 0.; | |
179 | tpar[2]= 0.; | |
180 | char volAlu[5]; // Alu | |
181 | char volBak[5]; // Bakelite | |
182 | char volGaz[5]; // Gas streamer | |
183 | sprintf(volAlu,"SC%dA",icount+1); | |
184 | sprintf(volBak,"SB%dA",icount+1); | |
185 | sprintf(volGaz,"SG%dA",icount+1); | |
186 | gMC->Gsvolu(volAlu,"BOX",idAlu1,tpar,0); // Al | |
187 | gMC->Gsvolu(volBak,"BOX",idtmed[1107],tpar,0); // Bakelite | |
188 | gMC->Gsvolu(volGaz,"BOX",idtmed[1106],tpar,0); // Gas streamer | |
189 | tpar[0] = -1.; | |
190 | tpar[1] = -1.; | |
191 | tpar[2] = 0.1; | |
192 | gMC->Gsposp(volGaz,1,volBak,0.,0.,0.,0,"ONLY",tpar,3); | |
193 | tpar[2] = 0.3; | |
194 | gMC->Gsposp(volBak,1,volAlu,0.,0.,0.,0,"ONLY",tpar,3); | |
195 | ||
196 | // chamber type A | |
197 | Float_t xEnv = (kDXZERO+kXMED+(kXMAX-kXMED)/2.)*zRatio; | |
198 | Float_t yEnvM = 0.; // y low position of envelope in chamber | |
199 | Float_t yEnvP = 0.; // y up position of envelope in chamber | |
200 | Float_t yEnvPsave = 0.; // tmp data | |
201 | Float_t yEnvMsave = 0.; // tmp data | |
202 | Float_t xpos = 0.; // x position of RPC in envelope | |
203 | Float_t ypos = 0.; // y position of RPC in envelope | |
204 | dpar[2] = 0.4; | |
205 | dpar[0] = ((kXMAX-kXMED)/2.)*zRatio; | |
206 | dpar[1] = kYMIN * zRatio; | |
207 | ||
208 | GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[4], true, TGeoTranslation(xEnv,yEnvM,kZm)); | |
209 | GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[13], true, TGeoTranslation(-xEnv,yEnvP,kZp)); | |
210 | ||
211 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[4],iVolNum++,3, dpar); | |
212 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[13],iVolNum++,3, dpar); | |
213 | ||
214 | // chamber type B (plus envelope chambers B & C) | |
215 | xEnv = (kDXZERO+kXMAX/2.)*zRatio; | |
216 | yEnvPsave = yEnvP; | |
217 | yEnvMsave = yEnvM; | |
218 | yEnvP = (yEnvMsave + kYMIN * zRatio ) * zpm + kYMIN * zRatio; | |
219 | yEnvM = (yEnvPsave + kYMIN * zRatio ) * zmp + kYMIN * zRatio; | |
220 | dpar[0] = ((kXMAX-kXMIN)/2.) * zRatio; | |
221 | dpar[1] = ((kYMAX-kYMIN)/2.) * zRatio; | |
222 | xpos = kXMIN/2. * zRatio; | |
223 | ypos = (kYMIN - kYMIN/4.) * zRatio; | |
224 | ||
225 | GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[3], true, TGeoTranslation( xEnv,-yEnvP,kZp)); | |
226 | GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[5], true, TGeoTranslation( xEnv, yEnvP,kZp)); | |
227 | GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[12], true, TGeoTranslation(-xEnv,-yEnvM,kZm)); | |
228 | GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[14], true, TGeoTranslation(-xEnv, yEnvM,kZm)); | |
229 | ||
230 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[3],iVolNum++,TGeoTranslation(xpos, ypos,0.),3,dpar); | |
231 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[5],iVolNum++,TGeoTranslation(xpos,-ypos,0.),3,dpar); | |
232 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[12],iVolNum++,TGeoTranslation(-xpos, ypos,0.),3,dpar); | |
233 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[14],iVolNum++,TGeoTranslation(-xpos,-ypos,0.),3,dpar); | |
234 | ||
235 | // chamber type C (note: same Z than type B) | |
236 | dpar[0] = (kXMAX/2)*zRatio; | |
237 | dpar[1] = (kYMAX/2)*zRatio; | |
238 | xpos = 0.; | |
239 | ypos = ((kYMAX - kYMIN)/2.) * zRatio; | |
240 | ||
241 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[3],iVolNum++,TGeoTranslation(xpos,-ypos,0.),3,dpar); | |
242 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[5],iVolNum++,TGeoTranslation(xpos, ypos,0.),3,dpar); | |
243 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[12],iVolNum++,TGeoTranslation(-xpos,-ypos,0.),3,dpar); | |
244 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[14],iVolNum++,TGeoTranslation(-xpos, ypos,0.),3,dpar); | |
245 | ||
246 | // chamber type D, E and F (same size) | |
247 | // D | |
248 | yEnvPsave = yEnvP; | |
249 | yEnvMsave = yEnvM; | |
250 | yEnvP = (yEnvMsave + kYMIN * zRatio ) * zpm + kYMIN * zRatio; | |
251 | yEnvM = (yEnvPsave + kYMIN * zRatio ) * zmp + kYMIN * zRatio; | |
252 | dpar[0] = (kXMAX/2.)*zRatio; | |
253 | dpar[1] = kYMIN*zRatio; | |
254 | ||
255 | GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[2], true, TGeoTranslation(xEnv, yEnvM,kZm)); | |
256 | GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[6], true, TGeoTranslation(-xEnv, yEnvP,kZp)); | |
257 | GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[11], true, TGeoTranslation(xEnv,-yEnvM,kZm)); | |
258 | GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[15], true, TGeoTranslation(-xEnv,-yEnvP,kZp)); | |
259 | ||
260 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[2],iVolNum++,3, dpar); | |
261 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[6],iVolNum++,3, dpar); | |
262 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[11],iVolNum++,3, dpar); | |
263 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[15],iVolNum++,3, dpar); | |
264 | ||
265 | // E | |
266 | yEnvPsave = yEnvP; | |
267 | yEnvMsave = yEnvM; | |
268 | yEnvP = (yEnvMsave + kYMIN * zRatio ) * zpm + kYMIN * zRatio; | |
269 | yEnvM = (yEnvPsave + kYMIN * zRatio ) * zmp + kYMIN * zRatio; | |
270 | ||
271 | GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[1], true, TGeoTranslation(xEnv, yEnvP,kZp)); | |
272 | GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[7], true, TGeoTranslation(-xEnv, yEnvM,kZm)); | |
273 | GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[10], true, TGeoTranslation(xEnv,-yEnvP,kZp)); | |
274 | GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[16], true, TGeoTranslation(-xEnv,-yEnvM,kZm)); | |
275 | ||
276 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[1],iVolNum++,3,dpar); | |
277 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[7],iVolNum++,3,dpar); | |
278 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[10],iVolNum++,3,dpar); | |
279 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[16],iVolNum++,3,dpar); | |
280 | ||
281 | ||
282 | // F | |
283 | yEnvPsave = yEnvP; | |
284 | yEnvMsave = yEnvM; | |
285 | yEnvP = (yEnvMsave + kYMIN * zRatio ) * zpm + kYMIN * zRatio; | |
286 | yEnvM = (yEnvPsave + kYMIN * zRatio ) * zmp + kYMIN * zRatio; | |
287 | ||
288 | GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[0], true, TGeoTranslation(xEnv, yEnvM,kZm)); | |
289 | GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[8], true, TGeoTranslation(-xEnv, yEnvP,kZp)); | |
290 | GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[9], true, TGeoTranslation(xEnv,-yEnvM,kZm)); | |
291 | GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[17], true, TGeoTranslation(-xEnv,-yEnvP,kZp)); | |
292 | ||
293 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[0],iVolNum++,3,dpar); | |
294 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[8],iVolNum++,3,dpar); | |
295 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[9],iVolNum++,3,dpar); | |
296 | GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[17],iVolNum++,3,dpar); | |
297 | ||
298 | } // end loop on detection planes | |
299 | } // end loop on stations | |
300 | } | |
301 | ||
302 | //______________________________________________________________________________ | |
303 | void AliMUONTriggerGeometryBuilder::SetTransformations() | |
304 | { | |
305 | // Defines the transformations for the trigger chambers. | |
306 | // --- | |
307 | Double_t zpos1, zpos2; | |
308 | AliMUONChamber *iChamber1, *iChamber2; | |
309 | ||
310 | iChamber1 = GetChamber(10); | |
311 | zpos1= - iChamber1->Z(); | |
312 | iChamber1->GetGeometry() | |
313 | ->SetTranslation(TGeoTranslation(0., 0., zpos1)); | |
314 | ||
315 | iChamber2 = GetChamber(11); | |
316 | zpos2 = - iChamber2->Z(); | |
317 | iChamber2->GetGeometry() | |
318 | ->SetTranslation(TGeoTranslation(0., 0., zpos2)); | |
319 | ||
320 | iChamber1 = GetChamber(12); | |
321 | zpos1 = - iChamber1->Z(); | |
322 | iChamber1->GetGeometry() | |
323 | ->SetTranslation(TGeoTranslation(0., 0., zpos1)); | |
324 | ||
325 | iChamber2 = GetChamber(13); | |
326 | zpos2 = - iChamber2->Z(); | |
327 | iChamber2->GetGeometry() | |
328 | ->SetTranslation(TGeoTranslation(0., 0., zpos2)); | |
329 | } | |
330 | ||
331 | //______________________________________________________________________________ | |
332 | void AliMUONTriggerGeometryBuilder::SetSensitiveVolumes() | |
333 | { | |
334 | // Defines the sensitive volumes for trigger station chambers. | |
335 | // --- | |
336 | ||
337 | GetChamber(10)->GetGeometry()->SetSensitiveVolume("SG1A"); | |
338 | GetChamber(11)->GetGeometry()->SetSensitiveVolume("SG2A"); | |
339 | GetChamber(12)->GetGeometry()->SetSensitiveVolume("SG3A"); | |
340 | GetChamber(13)->GetGeometry()->SetSensitiveVolume("SG4A"); | |
341 | } | |
342 |