3 // Class AliMUONTriggerGeometryBuilder
4 // -----------------------------------
5 // MUON Trigger stations geometry construction class.
7 // Author: Philippe Crochette, LPC Clermont-Ferrand
9 #include <TVirtualMC.h>
10 #include <TGeoMatrix.h>
12 #include "AliMUONTriggerGeometryBuilder.h"
14 #include "AliMUONChamber.h"
15 #include "AliMUONChamberGeometry.h"
17 ClassImp(AliMUONTriggerGeometryBuilder)
19 //______________________________________________________________________________
20 AliMUONTriggerGeometryBuilder::AliMUONTriggerGeometryBuilder(AliMUON* muon)
21 : AliMUONVGeometryBuilder(&muon->Chamber(10), &muon->Chamber(11),&muon->Chamber(12),&muon->Chamber(13)),
24 // Standard constructor
28 //______________________________________________________________________________
29 AliMUONTriggerGeometryBuilder::AliMUONTriggerGeometryBuilder()
30 : AliMUONVGeometryBuilder(),
33 // Default constructor
37 //______________________________________________________________________________
38 AliMUONTriggerGeometryBuilder::AliMUONTriggerGeometryBuilder(const AliMUONTriggerGeometryBuilder& rhs)
39 : AliMUONVGeometryBuilder(rhs)
41 Fatal("Copy constructor",
42 "Copy constructor is not implemented.");
45 //______________________________________________________________________________
46 AliMUONTriggerGeometryBuilder::~AliMUONTriggerGeometryBuilder() {
50 //______________________________________________________________________________
51 AliMUONTriggerGeometryBuilder&
52 AliMUONTriggerGeometryBuilder::operator = (const AliMUONTriggerGeometryBuilder& rhs)
54 // check assignement to self
55 if (this == &rhs) return *this;
58 "Assignment operator is not implemented.");
67 //______________________________________________________________________________
68 void AliMUONTriggerGeometryBuilder::CreateGeometry()
70 // From AliMUONv1::CreateGeometry()
73 zpos1 and zpos2 are the middle of the first and second
74 planes of station 1 (+1m for second station):
75 zpos1=(zpos1m+zpos1p)/2=(15999+16071)/2=16035 mm, thick/2=40 mm
76 zpos2=(zpos2m+zpos2p)/2=(16169+16241)/2=16205 mm, thick/2=40 mm
77 zposxm and zposxp= middles of gaz gaps within a detection plane
78 rem: the total thickness accounts for 1 mm of al on both
79 side of the RPCs (see zpos1 and zpos2)
82 Int_t *idtmed = fMUON->GetIdtmed()->GetArray()-1099;
83 Int_t idAir= idtmed[1100]; // medium 1
84 Int_t idAlu1=idtmed[1103]; // medium 4
88 // vertical gap between right and left chambers (kDXZERO*2=4cm)
89 const Float_t kDXZERO=2.;
90 // main distances for chamber definition in first plane/first station
91 const Float_t kXMIN=34.;
92 const Float_t kXMED=51.;
93 const Float_t kXMAX=272.;
94 // kXMAX will become 255. in real life. segmentation to be updated accordingly
95 // (see fig.2-4 & 2-5 of Local Trigger Board PRR)
96 const Float_t kYMIN=34.;
97 const Float_t kYMAX=51.;
98 // inner/outer radius of flange between beam shield. and chambers (1/station)
99 const Float_t kRMIN[2]={50.,50.};
100 const Float_t kRMAX[2]={64.,68.};
101 // z position of the middle of the gas gap in mother vol
102 const Float_t kZm=-3.6;
103 const Float_t kZp=+3.6;
105 AliMUONChamber *iChamber, *iChamber1;
106 iChamber1 = GetChamber(10);
107 Float_t zpos1=-iChamber1->Z();
109 // ratio of zpos1m/zpos1p and inverse for first plane
110 Float_t zmp=(zpos1-3.6)/(zpos1+3.6);
113 Int_t icount=0; // chamber counter (0 1 2 3)
115 for (Int_t istation=0; istation<2; istation++) { // loop on stations
116 for (Int_t iplane=0; iplane<2; iplane++) { // loop on detection planes
118 Int_t iVolNum=1; // counter Volume Number
119 icount = Int_t(iplane*TMath::Power(2,0))+
120 Int_t(istation*TMath::Power(2,1));
122 iChamber = GetChamber(10+icount);
123 Float_t zpos = - iChamber->Z();
125 // Flange between beam shielding and RPC
126 tpar[0]= kRMIN[istation];
127 tpar[1]= kRMAX[istation];
130 sprintf(volFlange,"SF%dA",icount+1);
131 gMC->Gsvolu(volFlange,"TUBE",idAlu1,tpar,3); // Al
133 //gMC->Gspos(volFlange,1,"ALIC",0.,0.,zpos,0,"MANY");
134 iChamber->GetGeometry()->AddEnvelope(volFlange, false, "MANY");
137 Float_t zRatio = zpos / zpos1;
139 // envelopes (same size except line 5, all virtual)
141 tpar[1] = kYMIN * zRatio;
143 Int_t i=0; // counter
144 for (Int_t icolumn=0; icolumn<2; icolumn++) {
145 for (Int_t iline=1; iline<10; iline++){
146 tpar[0] = (kXMAX/2.) * zRatio;
147 if (iline==5) tpar[0] = ((kXMAX-kXMED)/2.)*zRatio;
149 sprintf(volEnv[i],"S%dL%d",icount,iline);
151 sprintf(volEnv[i],"S%dR%d",icount,iline);
152 gMC->Gsvolu(volEnv[i],"BOX",idAir,tpar,0);
161 char volAlu[5]; // Alu
162 char volBak[5]; // Bakelite
163 char volGaz[5]; // Gas streamer
164 sprintf(volAlu,"SC%dA",icount+1);
165 sprintf(volBak,"SB%dA",icount+1);
166 sprintf(volGaz,"SG%dA",icount+1);
167 gMC->Gsvolu(volAlu,"BOX",idAlu1,tpar,0); // Al
168 gMC->Gsvolu(volBak,"BOX",idtmed[1107],tpar,0); // Bakelite
169 gMC->Gsvolu(volGaz,"BOX",idtmed[1106],tpar,0); // Gas streamer
173 gMC->Gsposp(volGaz,1,volBak,0.,0.,0.,0,"ONLY",tpar,3);
175 gMC->Gsposp(volBak,1,volAlu,0.,0.,0.,0,"ONLY",tpar,3);
178 Float_t xEnv = (kDXZERO+kXMED+(kXMAX-kXMED)/2.)*zRatio;
179 Float_t yEnvM = 0.; // y low position of envelope in chamber
180 Float_t yEnvP = 0.; // y up position of envelope in chamber
181 Float_t yEnvPsave = 0.; // tmp data
182 Float_t yEnvMsave = 0.; // tmp data
183 Float_t xpos = 0.; // x position of RPC in envelope
184 Float_t ypos = 0.; // y position of RPC in envelope
186 dpar[0] = ((kXMAX-kXMED)/2.)*zRatio;
187 dpar[1] = kYMIN * zRatio;
189 GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[4], true, TGeoTranslation(xEnv,yEnvM,kZm));
190 GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[13], true, TGeoTranslation(-xEnv,yEnvP,kZp));
192 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[4],iVolNum++,3, dpar);
193 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[13],iVolNum++,3, dpar);
195 // chamber type B (plus envelope chambers B & C)
196 xEnv = (kDXZERO+kXMAX/2.)*zRatio;
199 yEnvP = (yEnvMsave + kYMIN * zRatio ) * zpm + kYMIN * zRatio;
200 yEnvM = (yEnvPsave + kYMIN * zRatio ) * zmp + kYMIN * zRatio;
201 dpar[0] = ((kXMAX-kXMIN)/2.) * zRatio;
202 dpar[1] = ((kYMAX-kYMIN)/2.) * zRatio;
203 xpos = kXMIN/2. * zRatio;
204 ypos = (kYMIN - kYMIN/4.) * zRatio;
206 GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[3], true, TGeoTranslation( xEnv,-yEnvP,kZp));
207 GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[5], true, TGeoTranslation( xEnv, yEnvP,kZp));
208 GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[12], true, TGeoTranslation(-xEnv,-yEnvM,kZm));
209 GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[14], true, TGeoTranslation(-xEnv, yEnvM,kZm));
211 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[3],iVolNum++,TGeoTranslation(xpos, ypos,0.),3,dpar);
212 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[5],iVolNum++,TGeoTranslation(xpos,-ypos,0.),3,dpar);
213 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[12],iVolNum++,TGeoTranslation(-xpos, ypos,0.),3,dpar);
214 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[14],iVolNum++,TGeoTranslation(-xpos,-ypos,0.),3,dpar);
216 // chamber type C (note: same Z than type B)
217 dpar[0] = (kXMAX/2)*zRatio;
218 dpar[1] = (kYMAX/2)*zRatio;
220 ypos = ((kYMAX - kYMIN)/2.) * zRatio;
222 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[3],iVolNum++,TGeoTranslation(xpos,-ypos,0.),3,dpar);
223 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[5],iVolNum++,TGeoTranslation(xpos, ypos,0.),3,dpar);
224 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[12],iVolNum++,TGeoTranslation(-xpos,-ypos,0.),3,dpar);
225 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[14],iVolNum++,TGeoTranslation(-xpos, ypos,0.),3,dpar);
227 // chamber type D, E and F (same size)
231 yEnvP = (yEnvMsave + kYMIN * zRatio ) * zpm + kYMIN * zRatio;
232 yEnvM = (yEnvPsave + kYMIN * zRatio ) * zmp + kYMIN * zRatio;
233 dpar[0] = (kXMAX/2.)*zRatio;
234 dpar[1] = kYMIN*zRatio;
236 GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[2], true, TGeoTranslation(xEnv, yEnvM,kZm));
237 GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[6], true, TGeoTranslation(-xEnv, yEnvP,kZp));
238 GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[11], true, TGeoTranslation(xEnv,-yEnvM,kZm));
239 GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[15], true, TGeoTranslation(-xEnv,-yEnvP,kZp));
241 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[2],iVolNum++,3, dpar);
242 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[6],iVolNum++,3, dpar);
243 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[11],iVolNum++,3, dpar);
244 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[15],iVolNum++,3, dpar);
249 yEnvP = (yEnvMsave + kYMIN * zRatio ) * zpm + kYMIN * zRatio;
250 yEnvM = (yEnvPsave + kYMIN * zRatio ) * zmp + kYMIN * zRatio;
252 GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[1], true, TGeoTranslation(xEnv, yEnvP,kZp));
253 GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[7], true, TGeoTranslation(-xEnv, yEnvM,kZm));
254 GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[10], true, TGeoTranslation(xEnv,-yEnvP,kZp));
255 GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[16], true, TGeoTranslation(-xEnv,-yEnvM,kZm));
257 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[1],iVolNum++,3,dpar);
258 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[7],iVolNum++,3,dpar);
259 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[10],iVolNum++,3,dpar);
260 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[16],iVolNum++,3,dpar);
266 yEnvP = (yEnvMsave + kYMIN * zRatio ) * zpm + kYMIN * zRatio;
267 yEnvM = (yEnvPsave + kYMIN * zRatio ) * zmp + kYMIN * zRatio;
269 GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[0], true, TGeoTranslation(xEnv, yEnvM,kZm));
270 GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[8], true, TGeoTranslation(-xEnv, yEnvP,kZp));
271 GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[9], true, TGeoTranslation(xEnv,-yEnvM,kZm));
272 GetChamber(10+icount)->GetGeometry()->AddEnvelope(volEnv[17], true, TGeoTranslation(-xEnv,-yEnvP,kZp));
274 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[0],iVolNum++,3,dpar);
275 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[8],iVolNum++,3,dpar);
276 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[9],iVolNum++,3,dpar);
277 GetChamber(10+icount)->GetGeometry()->AddEnvelopeConstituentParam(volAlu,volEnv[17],iVolNum++,3,dpar);
279 } // end loop on detection planes
280 } // end loop on stations
283 //______________________________________________________________________________
284 void AliMUONTriggerGeometryBuilder::SetTransformations()
286 // Defines the transformations for the trigger chambers.
288 Double_t zpos1, zpos2;
289 AliMUONChamber *iChamber1, *iChamber2;
291 iChamber1 = GetChamber(10);
292 zpos1= - iChamber1->Z();
293 iChamber1->GetGeometry()
294 ->SetTranslation(TGeoTranslation(0., 0., zpos1));
296 iChamber2 = GetChamber(11);
297 zpos2 = - iChamber2->Z();
298 iChamber2->GetGeometry()
299 ->SetTranslation(TGeoTranslation(0., 0., zpos2));
301 iChamber1 = GetChamber(12);
302 zpos1 = - iChamber1->Z();
303 iChamber1->GetGeometry()
304 ->SetTranslation(TGeoTranslation(0., 0., zpos1));
306 iChamber2 = GetChamber(13);
307 zpos2 = - iChamber2->Z();
308 iChamber2->GetGeometry()
309 ->SetTranslation(TGeoTranslation(0., 0., zpos2));
312 //______________________________________________________________________________
313 void AliMUONTriggerGeometryBuilder::SetSensitiveVolumes()
315 // Defines the sensitive volumes for trigger station chambers.
318 GetChamber(10)->GetGeometry()->SetSensitiveVolume("SG1A");
319 GetChamber(11)->GetGeometry()->SetSensitiveVolume("SG2A");
320 GetChamber(12)->GetGeometry()->SetSensitiveVolume("SG3A");
321 GetChamber(13)->GetGeometry()->SetSensitiveVolume("SG4A");