db486a6e |
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 | |
b7943f00 |
16 | |
db486a6e |
17 | //************************************************************************* |
18 | // SDD geometry, based on ROOT geometrical modeler |
19 | // |
b7943f00 |
20 | // Ludovic Gaudichet gaudichet@to.infn.it |
db486a6e |
21 | //************************************************************************* |
22 | |
23 | |
db486a6e |
24 | |
25 | // General Root includes |
b7943f00 |
26 | //#include <Riostream.h> |
db486a6e |
27 | #include <TMath.h> |
28 | |
29 | // Root Geometry includes |
30 | #include <TGeoManager.h> |
31 | #include <TGeoVolume.h> |
db486a6e |
32 | #include <TGeoCone.h> |
33 | #include <TGeoTube.h> |
34 | #include <TGeoArb8.h> |
35 | #include <TGeoCompositeShape.h> |
36 | #include <TGeoMatrix.h> |
c789ee28 |
37 | #include <TGeoNode.h> |
db486a6e |
38 | |
b7943f00 |
39 | #include "AliITSgeom.h" |
40 | #include "AliITSgeomSDD.h" |
db486a6e |
41 | #include "AliITSv11GeometrySDD.h" |
b7943f00 |
42 | #include "AliITSv11GeomCableFlat.h" |
43 | #include "AliITSv11GeomCableRound.h" |
531d6cdc |
44 | #include "TList.h" |
b7943f00 |
45 | |
46 | const char* AliITSv11GeometrySDD::fgSDDsensitiveVolName = "ITSsddSensitiv"; |
47 | const Double_t AliITSv11GeometrySDD::fgkSegmentLength = 37.2*2*fgkmm; |
48 | const Double_t AliITSv11GeometrySDD::fgkLadderWidth = 50.0*fgkmm; |
49 | const Double_t AliITSv11GeometrySDD::fgkLadderHeight = 30.0*fgkmm; |
50 | const Double_t AliITSv11GeometrySDD::fgkLadderSegBoxDW = 3.5*fgkmm; |
51 | const Double_t AliITSv11GeometrySDD::fgkLadderSegBoxDH = 3.*fgkmm; |
52 | |
53 | const Double_t AliITSv11GeometrySDD::fgkLadderBeamRadius = 0.6*fgkmm; |
54 | const Double_t AliITSv11GeometrySDD::fgkLadderLa = 3.*fgkmm; |
fa4639a3 |
55 | const Double_t AliITSv11GeometrySDD::fgkLadderHa = 0.721979*fgkmm; |
b7943f00 |
56 | const Double_t AliITSv11GeometrySDD::fgkLadderLb = 3.7*fgkmm; |
fa4639a3 |
57 | const Double_t AliITSv11GeometrySDD::fgkLadderHb = 0.890428*fgkmm; |
b7943f00 |
58 | const Double_t AliITSv11GeometrySDD::fgkLadderl = 0.25*fgkmm; |
59 | |
60 | const Double_t AliITSv11GeometrySDD::fgkBottomBeamAngle = 56.5; |
61 | const Double_t AliITSv11GeometrySDD::fgkBeamSidePhi = 65; |
62 | |
63 | const Double_t AliITSv11GeometrySDD::fgkLadWaferSep = 2*fgkmm; |
fa4639a3 |
64 | const Double_t AliITSv11GeometrySDD::fgkPinSuppWidth = 2.5*fgkmm; |
65 | const Double_t AliITSv11GeometrySDD::fgkPinSuppHeight = 2.*fgkmm; |
b7943f00 |
66 | const Double_t AliITSv11GeometrySDD::fgkPinSuppRmax = 2.5/2.*fgkmm; |
67 | const Double_t AliITSv11GeometrySDD::fgkPinR = 1.5/2.*fgkmm; |
68 | const Double_t AliITSv11GeometrySDD::fgkPinSuppLength = 5.*fgkmm; |
69 | const Double_t AliITSv11GeometrySDD::fgkPinSuppThickness = 0.5*fgkmm; |
70 | const Double_t AliITSv11GeometrySDD::fgkPinSuppConeAngle = 4; |
71 | const Double_t AliITSv11GeometrySDD::fgkPinDXminOnSensor = (39./2.)*fgkmm; |
72 | const Double_t AliITSv11GeometrySDD::fgkPinPinDDXOnSensor = 3*fgkmm; |
73 | const Double_t AliITSv11GeometrySDD::fgkPinDYOnSensor = (52.5/2.)*fgkmm; |
db486a6e |
74 | |
b7943f00 |
75 | // parameters from ALR-0752/3 |
76 | const Double_t AliITSv11GeometrySDD::fgkCoolPipeSuppHeight = 3.2*fgkmm; |
77 | const Double_t AliITSv11GeometrySDD::fgkCoolPipeSuppMaxLength = 14*fgkmm; |
78 | const Double_t AliITSv11GeometrySDD::fgkCoolPipeSuppWidthExt = 0.4*fgkmm; |
79 | const Double_t AliITSv11GeometrySDD::fgkCoolPipeSuppWidthIn = 0.65*fgkmm; |
80 | const Double_t AliITSv11GeometrySDD::fgkCoolPipeSuppHoleDiam = 2*fgkmm; |
81 | const Double_t AliITSv11GeometrySDD::fgkCoolPipeSuppFulWidth = 5.15*fgkmm; |
82 | const Double_t AliITSv11GeometrySDD::fgkCoolPipeSuppTongW = 0.8*fgkmm; |
83 | const Double_t AliITSv11GeometrySDD::fgkCoolPipeSuppAngle = 22.5; |
84 | const Double_t AliITSv11GeometrySDD::fgkCoolPipeSuppSlitL = 4.9*fgkmm; |
85 | const Double_t AliITSv11GeometrySDD::fgkCoolPipeSuppAxeDist = 3.05*fgkmm; |
86 | const Double_t AliITSv11GeometrySDD::fgkCoolPipeInnerDiam = 1.84*fgkmm; |
87 | const Double_t AliITSv11GeometrySDD::fgkCoolPipeOuterDiam = 2.*fgkmm; |
88 | |
89 | const Double_t AliITSv11GeometrySDD::fgkBTBthick = 0.25 *fgkmm; |
90 | const Double_t AliITSv11GeometrySDD::fgkBTBlength = 55. *fgkmm; |
91 | const Double_t AliITSv11GeometrySDD::fgkBTBwidth = 18*fgkmm; |
92 | const Double_t AliITSv11GeometrySDD::fgkBTBaxisAtoBottom = 4*fgkmm; |
93 | const Double_t AliITSv11GeometrySDD::fgkBTBaxisAtoBase = 1.2*fgkmm; |
94 | const Double_t AliITSv11GeometrySDD::fgkRadiusAminBTB = 1. *fgkmm; |
95 | const Double_t AliITSv11GeometrySDD::fgkRadiusBminBTB = 0.53 *fgkmm; |
96 | const Double_t AliITSv11GeometrySDD::fgkBTBHoleLength = 15 *fgkmm; |
97 | const Double_t AliITSv11GeometrySDD::fgkBTBHolewidth = 6 *fgkmm; |
98 | const Double_t AliITSv11GeometrySDD::fgkBTBHoleRefX = 10 *fgkmm; |
99 | const Double_t AliITSv11GeometrySDD::fgkBTBHoleRefY = 6.5 *fgkmm; |
100 | |
101 | const Double_t AliITSv11GeometrySDD::fgkLay3Rmin = 145.*fgkmm; // not min! Rmin virtual tube |
102 | const Double_t AliITSv11GeometrySDD::fgkLay3Rmax = 205.*fgkmm; // not min! Rmax virtual tube |
103 | const Double_t AliITSv11GeometrySDD::fgkLay3Length = (524.+0.)*fgkmm; // ladder+supporting rings (length of the virtual tube) |
104 | const Double_t AliITSv11GeometrySDD::fgkLay3LadderLength = 524.*fgkmm; |
105 | const Double_t AliITSv11GeometrySDD::fgkLay3DetShortRadius = 146.0*fgkmm; |
106 | const Double_t AliITSv11GeometrySDD::fgkLay3DetLongRadius = 152.0*fgkmm; |
107 | const Double_t AliITSv11GeometrySDD::fgkLay3LaddTopCornerEnd = 15.6*fgkmm; |
108 | const Int_t AliITSv11GeometrySDD::fgkLay3Ndet = 6; |
109 | const Int_t AliITSv11GeometrySDD::fgkLay3Nladd = 14; |
110 | const Double_t AliITSv11GeometrySDD::fgkLay3CoolPipeSuppH = 7.5*fgkmm; |
111 | |
112 | const Double_t AliITSv11GeometrySDD::fgkLay4Rmin = 220.*fgkmm; // not min! Rmin virtual tube |
113 | const Double_t AliITSv11GeometrySDD::fgkLay4Rmax = 290.*fgkmm; // not min! Rmax virtual tube |
114 | const Double_t AliITSv11GeometrySDD::fgkLay4Length = (671.+0.)*fgkmm; // ladder+supporting rings (length of the virtual tube) |
115 | const Double_t AliITSv11GeometrySDD::fgkLay4LadderLength = 671.*fgkmm; |
116 | const Double_t AliITSv11GeometrySDD::fgkLay4DetShortRadius = 235.0*fgkmm; |
117 | const Double_t AliITSv11GeometrySDD::fgkLay4DetLongRadius = 240.5*fgkmm; |
118 | const Double_t AliITSv11GeometrySDD::fgkLay4LaddTopCornerEnd = 15.6*fgkmm; |
119 | const Int_t AliITSv11GeometrySDD::fgkLay4Ndet = 8; |
120 | const Int_t AliITSv11GeometrySDD::fgkLay4Nladd = 22; |
121 | const Double_t AliITSv11GeometrySDD::fgkLay4CoolPipeSuppH = 7.5*fgkmm; |
122 | |
123 | //hybrid |
124 | const Double_t AliITSv11GeometrySDD::fgkHybridAngle = 46; // approx !!! |
125 | // Origine taken at the hybrid corner : |
126 | const Double_t AliITSv11GeometrySDD::fgkHybridLength = 65*fgkmm; |
127 | const Double_t AliITSv11GeometrySDD::fgkHybridWidth = 41*fgkmm; |
128 | const Double_t AliITSv11GeometrySDD::fgkHybRndHoleRad = 1.05*fgkmm; |
129 | const Double_t AliITSv11GeometrySDD::fgkHybRndHoleZ = 2.5*fgkmm; |
130 | const Double_t AliITSv11GeometrySDD::fgkHybRndHoleX = fgkHybridWidth-23.599*fgkmm; |
131 | |
132 | const Double_t AliITSv11GeometrySDD::fgkHybFLlowHoleDZ = 9.698*fgkmm; |
133 | const Double_t AliITSv11GeometrySDD::fgkHybFLlowHolePasDX = 10.754*fgkmm; |
134 | const Double_t AliITSv11GeometrySDD::fgkHybFLlowHoleAmbDX = 9.122*fgkmm; |
135 | // center of ships to the border |
136 | const Double_t AliITSv11GeometrySDD::fgkHybFLlowChipZ4 = fgkHybridLength-(4.654 )*fgkmm-fgkHybFLlowHoleDZ/2; |
137 | const Double_t AliITSv11GeometrySDD::fgkHybFLlowChipZ3 = fgkHybridLength-(4.654+15. )*fgkmm-fgkHybFLlowHoleDZ/2; |
138 | const Double_t AliITSv11GeometrySDD::fgkHybFLlowChipZ2 = fgkHybridLength-(4.654+15.*2)*fgkmm-fgkHybFLlowHoleDZ/2; |
139 | const Double_t AliITSv11GeometrySDD::fgkHybFLlowChipZ1 = fgkHybridLength-(4.654+15.*3)*fgkmm-fgkHybFLlowHoleDZ/2; |
140 | const Double_t AliITSv11GeometrySDD::fgkHybFLlowPasX = fgkHybridWidth-32.775*fgkmm; |
141 | const Double_t AliITSv11GeometrySDD::fgkHybFLlowAmbX = fgkHybridWidth-20.791*fgkmm; |
142 | const Double_t AliITSv11GeometrySDD::fgkHybChipsDZ = 9.221*fgkmm; |
143 | const Double_t AliITSv11GeometrySDD::fgkHybPascalDX = 10.245*fgkmm; |
144 | const Double_t AliITSv11GeometrySDD::fgkHybAmbraDX = 8.51*fgkmm; |
145 | const Double_t AliITSv11GeometrySDD::fgkHybFLUpperWidth = 15.012*fgkmm; |
146 | const Double_t AliITSv11GeometrySDD::fgkHybFLUpperLength = 59.878*fgkmm; |
147 | const Double_t AliITSv11GeometrySDD::fgkHybFLUpperAlDZ = 11.183*fgkmm; |
148 | const Double_t AliITSv11GeometrySDD::fgkHybFLUpperAldx = 2.307*fgkmm; |
531d6cdc |
149 | |
150 | const Double_t AliITSv11GeometrySDD::fgkmu = 1*fgkmicron; // 1*fgkmicron; // can be increase for checking thin objects |
b7943f00 |
151 | const Double_t AliITSv11GeometrySDD::fgkHybridThBridgeThick = 0.25*fgkmm; // ??? |
152 | const Double_t AliITSv11GeometrySDD::fgkHybAlThick = 30*fgkmu; |
153 | const Double_t AliITSv11GeometrySDD::fgkHybUpThick = 20*fgkmu; |
154 | const Double_t AliITSv11GeometrySDD::fgkHybGlueScrnThick = 50*fgkmu; // ??? ????? |
155 | const Double_t AliITSv11GeometrySDD::fgkHybGlueLowThick = 90*fgkmu; |
156 | const Double_t AliITSv11GeometrySDD::fgkHybGlueUpThick = 90*fgkmu; // sur ????? |
157 | const Double_t AliITSv11GeometrySDD::fgkHybAlCCThick = 12*fgkmu; |
158 | const Double_t AliITSv11GeometrySDD::fgkHybUpCCThick = 12*fgkmu; |
159 | const Double_t AliITSv11GeometrySDD::fgkHybChipThick = 150*fgkmu; |
160 | const Double_t AliITSv11GeometrySDD::fgkHybGlueAgThick = 50*fgkmu; // ??? ???? |
161 | const Double_t AliITSv11GeometrySDD::fgkHybUnderNiThick = 20*fgkmu; // ??? ???? |
162 | const Int_t AliITSv11GeometrySDD::fgkNHybSMD = 25; |
163 | const Double_t AliITSv11GeometrySDD::fgkHybSMDposX[fgkNHybSMD] = |
164 | {2.92*fgkmm,6.5*fgkmm,10.13*fgkmm,13.59*fgkmm,21.40*fgkmm, |
165 | 2.92*fgkmm,6.5*fgkmm,10.13*fgkmm,13.59*fgkmm,19.91*fgkmm, |
166 | 2.92*fgkmm,6.5*fgkmm,10.13*fgkmm,13.59*fgkmm,17.09*fgkmm,21.40*fgkmm, |
167 | 2.92*fgkmm,6.5*fgkmm,10.13*fgkmm,13.59*fgkmm,19.91*fgkmm, |
168 | 1.63*fgkmm,5.22*fgkmm,13.59*fgkmm,21.40*fgkmm}; |
169 | const Double_t AliITSv11GeometrySDD::fgkHybSMDposZ[fgkNHybSMD] = |
170 | { 2.3 *fgkmm, 2.3 *fgkmm, 2.3 *fgkmm, 2.3 *fgkmm, 2.3 *fgkmm, |
171 | 17.315*fgkmm,17.315*fgkmm,17.315*fgkmm,17.315*fgkmm,17.315*fgkmm, |
172 | 32.31*fgkmm,32.31*fgkmm,32.31*fgkmm,32.31*fgkmm,32.31*fgkmm,32.31*fgkmm, |
173 | 47.38*fgkmm,47.38*fgkmm,47.38*fgkmm,47.38*fgkmm,47.38*fgkmm, |
174 | 62.68*fgkmm,62.06*fgkmm,62.06*fgkmm,62.06*fgkmm}; |
175 | const Double_t AliITSv11GeometrySDD::fgkHybSMDmiddleW = 0.954*fgkmm; |
176 | const Double_t AliITSv11GeometrySDD::fgkHybSMDmiddleL = 0.47 *fgkmm; |
177 | const Double_t AliITSv11GeometrySDD::fgkHybSMDendW = 1.132*fgkmm; |
178 | const Double_t AliITSv11GeometrySDD::fgkHybSMDendL = 0.925*fgkmm; |
179 | const Double_t AliITSv11GeometrySDD::fgkHybSMDheight = 400.*fgkmu; // ??? ????!!!!!!! |
180 | |
181 | const Double_t AliITSv11GeometrySDD::fgkWaferThickness = 300.*fgkmu; |
182 | const Double_t AliITSv11GeometrySDD::fgkWaferWidth = 72.5 *fgkmm; |
183 | const Double_t AliITSv11GeometrySDD::fgkWaferLength = 87.6 *fgkmm; |
184 | const Double_t AliITSv11GeometrySDD::fgkWaferThickSens = 299.8*fgkmu; |
185 | const Double_t AliITSv11GeometrySDD::fgkWaferWidthSens = 70.17*fgkmm; |
186 | const Double_t AliITSv11GeometrySDD::fgkWaferLengthSens = 74.97*fgkmm; |
187 | |
188 | const Double_t AliITSv11GeometrySDD::fgkDigitCablWidth = 18.4*fgkmm; |
189 | const Double_t AliITSv11GeometrySDD::fgkDigitCablAlThick = (30+30*8./10.)*fgkmicron; // will probably change |
190 | const Double_t AliITSv11GeometrySDD::fgkDigitCablPolyThick = (20+12)*fgkmicron; // will probably change |
191 | |
192 | const Double_t AliITSv11GeometrySDD::fgkWaHVcableAlThick = 30*2./10.*fgkmu; // will probably change // Al ratio is random !!! |
193 | const Double_t AliITSv11GeometrySDD::fgkWaHVcablePolyThick = 175*fgkmu; // will probably change |
194 | const Double_t AliITSv11GeometrySDD::fgkWaHVcableLength = 67.08*fgkmm; |
195 | const Double_t AliITSv11GeometrySDD::fgkWaHVcableWitdh = 17.4 *fgkmm; // check !!! |
531d6cdc |
196 | const Double_t AliITSv11GeometrySDD::fgkWaHVcableDW = 13.*fgkmm; //5.24*fgkmm; // check !!! |
b7943f00 |
197 | |
198 | const Double_t AliITSv11GeometrySDD::fgkSensorGlassLX = 5. *fgkmm; |
199 | const Double_t AliITSv11GeometrySDD::fgkSensorGlassLZ = 5. *fgkmm; |
200 | const Double_t AliITSv11GeometrySDD::fgkSensorGlassLY = 150. *fgkmu; |
201 | const Double_t AliITSv11GeometrySDD::fgkGlassDXOnSensor = 26.28*fgkmm; // check !!! |
202 | const Double_t AliITSv11GeometrySDD::fgkGlassDZOnSensor = 22.50*fgkmm; // check !!! |
203 | |
204 | const Double_t AliITSv11GeometrySDD::fgkTransitHVAlThick = 30*2./10.*fgkmu; // check // will probably change //Al ratio is random |
205 | const Double_t AliITSv11GeometrySDD::fgkTransitHVPolyThick = 100*fgkmu; // check // will probably change |
206 | const Double_t AliITSv11GeometrySDD::fgkTransitHVHeadLX = 71.46*fgkmm; // check !!! |
207 | const Double_t AliITSv11GeometrySDD::fgkTransitHVHeadLZ = 21.3*fgkmm; |
208 | const Double_t AliITSv11GeometrySDD::fgkTransitHVBondingLZ = 3.6*fgkmm; |
209 | const Double_t AliITSv11GeometrySDD::fgkTransitHVtailLength = 27*fgkmm; // ???, not yet fixed ... |
210 | const Double_t AliITSv11GeometrySDD::fgkTransitHVtailWidth = 26*fgkmm; |
211 | const Double_t AliITSv11GeometrySDD::fgkTransitHVtailXpos = 8*fgkmm; // ???, a mesurer !!! |
212 | const Double_t AliITSv11GeometrySDD::fgkTransitHVsideLZ = 10.34*fgkmm; |
213 | const Double_t AliITSv11GeometrySDD::fgkTransitHVsideLeftZ = 4.11*fgkmm; |
214 | const Double_t AliITSv11GeometrySDD::fgkTransitHVsideRightZ = 3.5*fgkmm; // ???, a mesurer !!! |
215 | |
216 | const Double_t AliITSv11GeometrySDD::fgkLongHVcablePolyThick= (20+30+125+30+20+30+125+30+20)*fgkmu; // check // will probably change |
217 | const Double_t AliITSv11GeometrySDD::fgkLongHVcableAlThick = (30+30*2/10+30)*fgkmu; // check // will probably change |
218 | const Double_t AliITSv11GeometrySDD::fgkLongHVcableSeparation = 600*fgkmicron; |
c789ee28 |
219 | |
c789ee28 |
220 | |
b7943f00 |
221 | ClassImp(AliITSv11GeometrySDD) |
c789ee28 |
222 | |
b7943f00 |
223 | //________________________________________________________________________ |
33ddec7d |
224 | AliITSv11GeometrySDD::AliITSv11GeometrySDD(): |
225 | AliITSv11Geometry(), |
226 | fPinSupport(0), |
227 | fCoolPipeSupportL(0), |
228 | fCoolPipeSupportR(0), |
229 | fSDDsensor(0), |
230 | fBaseThermalBridge(0), |
231 | fHybrid(0), |
232 | fDigitCableLay3A(0), |
233 | fDigitCableLay3B(0), |
234 | fDigitCableLay4A(0), |
235 | fDigitCableLay4B(0), |
236 | fMotherVol(0), |
237 | fAddHybrids(kTRUE), |
238 | fAddSensors(kTRUE), |
239 | fAddHVcables(kTRUE), |
240 | fAddCables(kTRUE), |
241 | fAddCoolingSyst(kTRUE), |
242 | fCoolingOn(kTRUE), |
243 | fAddOnlyLadder3min(-1), |
244 | fAddOnlyLadder3max(-1), |
245 | fAddOnlyLadder4min(-1), |
246 | fAddOnlyLadder4max(-1), |
247 | fColorCarbonFiber(2), |
248 | fColorRyton(2), |
249 | fColorPhynox(2), |
250 | fColorSilicon(2), |
251 | fColorAl(2), |
252 | fColorPolyhamide(2), |
253 | fColorGlass(2), |
254 | fColorSMD(2), |
255 | fColorSMDweld(2), |
256 | fLay3LadderUnderSegDH(0), |
257 | fLay4LadderUnderSegDH(0), |
258 | fLay3LaddShortRadius(0), |
259 | fLay3LaddLongRadius(0), |
260 | fLay4LaddShortRadius(0), |
261 | fLay4LaddLongRadius(0) |
c789ee28 |
262 | { |
b7943f00 |
263 | // |
264 | // Standard constructor |
265 | // |
b7943f00 |
266 | SetParameters(); |
531d6cdc |
267 | } |
c789ee28 |
268 | |
b7943f00 |
269 | //________________________________________________________________________ |
270 | AliITSv11GeometrySDD::AliITSv11GeometrySDD(Int_t debug) : |
33ddec7d |
271 | AliITSv11Geometry(debug), |
272 | fPinSupport(0), |
273 | fCoolPipeSupportL(0), |
274 | fCoolPipeSupportR(0), |
275 | fSDDsensor(0), |
276 | fBaseThermalBridge(0), |
277 | fHybrid(0), |
278 | fDigitCableLay3A(0), |
279 | fDigitCableLay3B(0), |
280 | fDigitCableLay4A(0), |
281 | fDigitCableLay4B(0), |
282 | fMotherVol(0), |
283 | fAddHybrids(kTRUE), |
284 | fAddSensors(kTRUE), |
285 | fAddHVcables(kTRUE), |
286 | fAddCables(kTRUE), |
287 | fAddCoolingSyst(kTRUE), |
b7943f00 |
288 | fCoolingOn(kTRUE), |
33ddec7d |
289 | fAddOnlyLadder3min(-1), |
290 | fAddOnlyLadder3max(-1), |
291 | fAddOnlyLadder4min(-1), |
292 | fAddOnlyLadder4max(-1), |
293 | fColorCarbonFiber(2), |
294 | fColorRyton(2), |
295 | fColorPhynox(2), |
296 | fColorSilicon(2), |
297 | fColorAl(2), |
298 | fColorPolyhamide(2), |
299 | fColorGlass(2), |
300 | fColorSMD(2), |
301 | fColorSMDweld(2), |
302 | fLay3LadderUnderSegDH(0), |
303 | fLay4LadderUnderSegDH(0), |
304 | fLay3LaddShortRadius(0), |
305 | fLay3LaddLongRadius(0), |
306 | fLay4LaddShortRadius(0), |
307 | fLay4LaddLongRadius(0) |
c789ee28 |
308 | { |
b7943f00 |
309 | // |
310 | // Constructor setting debugging level |
311 | // |
b7943f00 |
312 | SetParameters(); |
531d6cdc |
313 | } |
c789ee28 |
314 | |
c789ee28 |
315 | //________________________________________________________________________ |
b7943f00 |
316 | AliITSv11GeometrySDD::AliITSv11GeometrySDD(const AliITSv11GeometrySDD &s) : |
33ddec7d |
317 | AliITSv11Geometry(s.GetDebug()), |
318 | fPinSupport(s.fPinSupport), |
319 | fCoolPipeSupportL(s.fCoolPipeSupportL), |
320 | fCoolPipeSupportR(s.fCoolPipeSupportR), |
321 | fSDDsensor(s.fSDDsensor), |
322 | fBaseThermalBridge(s.fBaseThermalBridge), |
323 | fHybrid(s.fHybrid), |
324 | fDigitCableLay3A(s.fDigitCableLay3A), |
325 | fDigitCableLay3B(s.fDigitCableLay3B), |
326 | fDigitCableLay4A(s.fDigitCableLay4A), |
327 | fDigitCableLay4B(s.fDigitCableLay4B), |
328 | fMotherVol(s.fMotherVol), |
329 | fAddHybrids(s.fAddHybrids), |
330 | fAddSensors(s.fAddSensors), |
331 | fAddHVcables(s.fAddHVcables), |
332 | fAddCables(s.fAddCables), |
333 | fAddCoolingSyst(s.fAddCoolingSyst), |
334 | fCoolingOn(s.fCoolingOn), |
335 | fAddOnlyLadder3min(s.fAddOnlyLadder3min), |
336 | fAddOnlyLadder3max(s.fAddOnlyLadder3max), |
337 | fAddOnlyLadder4min(s.fAddOnlyLadder4min), |
338 | fAddOnlyLadder4max(s.fAddOnlyLadder4max), |
339 | fColorCarbonFiber(s.fColorCarbonFiber), |
340 | fColorRyton(s.fColorRyton), |
341 | fColorPhynox(s.fColorPhynox), |
342 | fColorSilicon(s.fColorSilicon), |
343 | fColorAl(s.fColorAl), |
344 | fColorPolyhamide(s.fColorPolyhamide), |
345 | fColorGlass(s.fColorGlass), |
346 | fColorSMD(s.fColorSMD), |
347 | fColorSMDweld(s.fColorSMDweld), |
348 | fLay3LadderUnderSegDH(s.fLay3LadderUnderSegDH), |
349 | fLay4LadderUnderSegDH(s.fLay4LadderUnderSegDH), |
350 | fLay3LaddShortRadius(s.fLay3LaddShortRadius), |
351 | fLay3LaddLongRadius(s.fLay3LaddLongRadius), |
352 | fLay4LaddShortRadius(s.fLay4LaddShortRadius), |
353 | fLay4LaddLongRadius(s.fLay4LaddLongRadius) |
b7943f00 |
354 | { |
33ddec7d |
355 | // Copy Constructor |
356 | // do only a "shallow copy" ... |
b7943f00 |
357 | SetParameters(); |
358 | } |
c789ee28 |
359 | |
b7943f00 |
360 | //________________________________________________________________________ |
361 | AliITSv11GeometrySDD& AliITSv11GeometrySDD:: |
362 | operator=(const AliITSv11GeometrySDD &s) { |
363 | // Assignment operator |
364 | if(&s == this) return *this; |
365 | fMotherVol = s.fMotherVol; |
366 | fAddHybrids = s.fAddHybrids; |
367 | fAddSensors = s.fAddSensors; |
368 | fAddHVcables = s.fAddHVcables; |
369 | fAddCables = s.fAddCables; |
370 | fAddCoolingSyst = s.fAddCoolingSyst; |
371 | fCoolingOn = s.fCoolingOn; |
372 | fAddOnlyLadder3min = s.fAddOnlyLadder3min; |
373 | fAddOnlyLadder3max = s.fAddOnlyLadder3max; |
374 | fAddOnlyLadder4min = s.fAddOnlyLadder4min; |
375 | fAddOnlyLadder4max = s.fAddOnlyLadder4max; |
376 | return *this; |
377 | } |
db486a6e |
378 | |
c789ee28 |
379 | //________________________________________________________________________ |
b7943f00 |
380 | AliITSv11GeometrySDD::~AliITSv11GeometrySDD() { |
381 | // Look like a destructor |
382 | // Smell like a destructor |
383 | // And actually is the destructor |
384 | if (fDigitCableLay3A) delete [] fDigitCableLay3A; |
385 | if (fDigitCableLay3B) delete [] fDigitCableLay3B; |
386 | if (fDigitCableLay4A) delete [] fDigitCableLay4A; |
387 | if (fDigitCableLay4B) delete [] fDigitCableLay4B; |
531d6cdc |
388 | } |
c789ee28 |
389 | |
c789ee28 |
390 | //________________________________________________________________________ |
b7943f00 |
391 | void AliITSv11GeometrySDD::SetParameters() { |
392 | // |
393 | // Define display colors and the non constant geometry parameters |
394 | // |
395 | |
fa4639a3 |
396 | fColorCarbonFiber = 4; |
397 | fColorRyton = 5; |
398 | fColorPhynox = 7; |
399 | fColorSilicon = 3; |
400 | fColorAl = 7; |
401 | fColorPolyhamide = 5; |
402 | fColorGlass = 2; |
403 | fColorSMD = 12; |
404 | fColorSMDweld = 17; |
405 | |
406 | fPinSupport = 0; |
b7943f00 |
407 | fCoolPipeSupportL = 0; |
408 | fCoolPipeSupportR = 0; |
fa4639a3 |
409 | fSDDsensor = 0; |
410 | fBaseThermalBridge= 0; |
411 | fHybrid = 0; |
b7943f00 |
412 | |
413 | Double_t detLadderDist = 8*fgkmm; |
414 | |
415 | fLay3LadderUnderSegDH = detLadderDist - (fgkWaHVcableAlThick+fgkWaHVcablePolyThick); |
416 | fLay4LadderUnderSegDH = detLadderDist - (fgkWaHVcableAlThick+fgkWaHVcablePolyThick); |
417 | |
418 | // radius from the center to the CF ladder : |
419 | fLay3LaddShortRadius = (fgkLay3DetShortRadius |
420 | + fgkLadWaferSep+2*fgkWaferThickness |
421 | + detLadderDist); |
422 | fLay3LaddLongRadius = (fgkLay3DetLongRadius |
423 | + fgkLadWaferSep+2*fgkWaferThickness |
424 | + detLadderDist); |
425 | fLay4LaddShortRadius = (fgkLay4DetShortRadius |
426 | + fgkLadWaferSep+2*fgkWaferThickness |
427 | + detLadderDist); |
428 | fLay4LaddLongRadius = (fgkLay4DetLongRadius |
429 | + fgkLadWaferSep+2*fgkWaferThickness |
430 | + detLadderDist); |
c789ee28 |
431 | |
432 | fLay3sensorZPos[0]= ( 35.8+72.4+75.8 )*fgkmm; |
433 | fLay3sensorZPos[1]= ( 35.8+72.4 )*fgkmm; |
434 | fLay3sensorZPos[2]= ( 35.8 )*fgkmm; |
435 | fLay3sensorZPos[3]= ( -37.9 )*fgkmm; |
436 | fLay3sensorZPos[4]= ( -37.9-74.9 )*fgkmm; |
437 | fLay3sensorZPos[5]= ( -37.9-74.9-71.1 )*fgkmm; |
438 | |
439 | fLay4sensorZPos[0] = ( 38.5+73.2+75.4+71.6 )*fgkmm; |
440 | fLay4sensorZPos[1] = ( 38.5+73.2+75.4 )*fgkmm; |
441 | fLay4sensorZPos[2] = ( 38.5+73.2 )*fgkmm; |
442 | fLay4sensorZPos[3] = ( 38.5 )*fgkmm; |
443 | fLay4sensorZPos[4] = ( -35.6 )*fgkmm; |
444 | fLay4sensorZPos[5] = ( -35.6-74.8 )*fgkmm; |
445 | fLay4sensorZPos[6] = ( -35.6-74.8-72.4 )*fgkmm; |
446 | fLay4sensorZPos[7] = ( -35.6-74.8-72.4-76. )*fgkmm; |
531d6cdc |
447 | } |
448 | |
c789ee28 |
449 | |
c789ee28 |
450 | //________________________________________________________________________ |
b7943f00 |
451 | TGeoMedium* AliITSv11GeometrySDD::GetMedium(const char* mediumName) { |
452 | // |
453 | // Called to get a medium, checks that it exists. |
454 | // If not, prints an error and returns 0 |
455 | // |
c789ee28 |
456 | |
531d6cdc |
457 | char ch[30]; |
458 | sprintf(ch, "ITS_%s",mediumName); |
459 | TGeoMedium* medium = gGeoManager->GetMedium(ch); |
b7943f00 |
460 | if (! medium) |
461 | printf("Error(AliITSv11GeometrySDD)::medium %s not found !\n", mediumName); |
b7943f00 |
462 | return medium; |
531d6cdc |
463 | } |
464 | |
c789ee28 |
465 | |
fa4639a3 |
466 | //________________________________________________________________________ |
467 | Int_t AliITSv11GeometrySDD::GetLay3NLadders() const{ |
468 | // Get the actual number of ladder in layer 3 |
469 | if ( (fAddOnlyLadder3min<0) || |
470 | (fAddOnlyLadder3min >= fgkLay3Nladd) || |
471 | (fAddOnlyLadder3max<0) || |
472 | (fAddOnlyLadder3max >= fgkLay3Nladd) ) |
473 | return fgkLay3Nladd; |
474 | else return (fAddOnlyLadder3max-fAddOnlyLadder3min+1); |
475 | } |
476 | |
531d6cdc |
477 | |
fa4639a3 |
478 | //________________________________________________________________________ |
479 | Int_t AliITSv11GeometrySDD::GetLay4NLadders() const{ |
480 | // Get the actual number of ladder in layer 4 |
481 | if ( (fAddOnlyLadder4min<0) || |
482 | (fAddOnlyLadder4min >= fgkLay4Nladd) || |
483 | (fAddOnlyLadder4max<0) || |
484 | (fAddOnlyLadder4max >= fgkLay4Nladd) ) |
485 | return fgkLay4Nladd; |
486 | else return (fAddOnlyLadder4max-fAddOnlyLadder4min+1); |
487 | } |
488 | |
531d6cdc |
489 | |
b7943f00 |
490 | //________________________________________________________________________ |
491 | void AliITSv11GeometrySDD::CreateBasicObjects() { |
492 | // |
493 | // Create basics objets which will be assembled together |
494 | // in Layer3 and Layer4 functions |
495 | // |
496 | |
531d6cdc |
497 | |
fa4639a3 |
498 | fDigitCableLay3A = new AliITSv11GeomCableFlat[fgkLay3Ndet]; |
499 | fDigitCableLay3B = new AliITSv11GeomCableFlat[fgkLay3Ndet]; |
500 | fDigitCableLay4A = new AliITSv11GeomCableFlat[fgkLay4Ndet]; |
501 | fDigitCableLay4B = new AliITSv11GeomCableFlat[fgkLay4Ndet]; |
502 | |
b7943f00 |
503 | fPinSupport = CreatePinSupport(); |
504 | fCoolPipeSupportL = CreateCoolPipeSupportL(); |
505 | fCoolPipeSupportR = CreateCoolPipeSupportR(); |
506 | fSDDsensor = CreateSDDsensor(); |
507 | fBaseThermalBridge = CreateBaseThermalBridge(); |
508 | fHybrid = CreateHybrid(0); |
509 | |
531d6cdc |
510 | |
b7943f00 |
511 | TGeoMedium *carbonFiberLadderStruct = GetMedium("ITSsddCarbonM55J"); |
512 | TGeoMedium *polyhamideSDD = GetMedium("ITSsddKAPTON_POLYCH2"); |
513 | TGeoMedium *alSDD = GetMedium("ITSal"); |
514 | |
515 | //******************************************************************** |
516 | // pieces of the carbon fiber structure |
517 | //******************************************************************** |
518 | Double_t dy = fgkLadderSegBoxDH/2; |
519 | Double_t triangleHeight = fgkLadderHeight - fgkLadderBeamRadius; |
520 | Double_t halfTheta = TMath::ATan( 0.5*fgkLadderWidth/triangleHeight ); |
521 | Double_t alpha = TMath::Pi()*3./4. - halfTheta/2.; |
522 | Double_t beta = (TMath::Pi() - 2.*halfTheta)/4.; |
523 | Double_t dYTranslation = (fgkLadderHeight/2. |
524 | -0.5*fgkLadderWidth*TMath::Tan(beta) |
525 | -fgkLadderBeamRadius); |
526 | Double_t distCenterSideDown = 0.5*fgkLadderWidth/TMath::Cos(beta); |
527 | |
528 | //--- the top V of the Carbon Fiber Ladder (segment) |
529 | TGeoArb8 *cfLaddTop1 = CreateLadderSide( fgkSegmentLength/2., halfTheta, |
530 | -1, fgkLadderLa, fgkLadderHa, fgkLadderl); |
531 | TGeoVolume *cfLaddTopVol1 = new TGeoVolume("ITSsddCFladdTopCornerVol1", |
532 | cfLaddTop1,carbonFiberLadderStruct); |
533 | TGeoArb8 *cfLaddTop2 = CreateLadderSide( fgkSegmentLength/2., halfTheta, |
534 | 1, fgkLadderLa, fgkLadderHa, fgkLadderl); |
535 | TGeoVolume *cfLaddTopVol2 = new TGeoVolume("ITSsddCFladdTopCornerVol2", |
536 | cfLaddTop2, carbonFiberLadderStruct); |
537 | cfLaddTopVol1->SetLineColor(fColorCarbonFiber); |
538 | cfLaddTopVol2->SetLineColor(fColorCarbonFiber); |
539 | TGeoTranslation *trTop1 = new TGeoTranslation(0, fgkLadderHeight/2-dy, 0); |
540 | |
541 | //--- the 2 side V |
542 | TGeoArb8 *cfLaddSide1 = CreateLadderSide( fgkSegmentLength/2., beta, -1, |
543 | fgkLadderLb, fgkLadderHb, fgkLadderl); |
544 | TGeoVolume *cfLaddSideVol1 = new TGeoVolume( "ITSsddCFladdSideCornerVol1", |
545 | cfLaddSide1,carbonFiberLadderStruct); |
546 | TGeoArb8 *cfLaddSide2 = CreateLadderSide( fgkSegmentLength/2., beta, 1, |
547 | fgkLadderLb, fgkLadderHb, fgkLadderl); |
548 | TGeoVolume *cfLaddSideVol2 = new TGeoVolume( "ITSsddCFladdSideCornerVol2", |
549 | cfLaddSide2,carbonFiberLadderStruct); |
550 | cfLaddSideVol1->SetLineColor(fColorCarbonFiber); |
551 | cfLaddSideVol2->SetLineColor(fColorCarbonFiber); |
552 | TGeoCombiTrans *ctSideR = CreateCombiTrans("", distCenterSideDown, 0, |
553 | alpha*TMath::RadToDeg()); |
554 | AddTranslationToCombiTrans(ctSideR, 0, -dYTranslation-dy, 0); |
555 | TGeoCombiTrans *ctSideL = CreateCombiTrans("", distCenterSideDown,0, |
556 | -alpha*TMath::RadToDeg()); |
557 | AddTranslationToCombiTrans(ctSideL, 0, -dYTranslation-dy, 0); |
558 | |
559 | //--- The beams |
560 | // Beams on the sides |
561 | Double_t beamPhiPrime = TMath::ASin(1./TMath::Sqrt( (1+TMath::Sin(2*beta)* |
562 | TMath::Sin(2*beta)/(TanD(fgkBeamSidePhi)*TanD(fgkBeamSidePhi))) )); |
563 | //cout<<"Phi prime = "<<beamPhiPrime*TMath::RadToDeg()<<endl; |
564 | Double_t beamLength = TMath::Sqrt( fgkLadderHeight*fgkLadderHeight/ |
565 | ( TMath::Sin(beamPhiPrime)*TMath::Sin(beamPhiPrime)) |
566 | + fgkLadderWidth*fgkLadderWidth/4.)-fgkLadderLa/2-fgkLadderLb/2; |
567 | TGeoTubeSeg *sideBeamS = new TGeoTubeSeg(0, fgkLadderBeamRadius,beamLength/2., |
568 | 0, 180); |
569 | TGeoVolume *sideBeam = new TGeoVolume("ITSsddCFSideBeamVol", sideBeamS, |
570 | carbonFiberLadderStruct); |
571 | sideBeam->SetLineColor(fColorCarbonFiber); |
572 | |
573 | //Euler rotation : about Z, then new X, then new Z |
531d6cdc |
574 | TGeoRotation *beamRot1 = new TGeoRotation("", 90-2.*beta*TMath::RadToDeg(), |
575 | -beamPhiPrime*TMath::RadToDeg(),-90); |
576 | TGeoRotation *beamRot2 = new TGeoRotation("", 90-2.*beta*TMath::RadToDeg(), |
577 | beamPhiPrime*TMath::RadToDeg(), -90); |
578 | TGeoRotation *beamRot3 = new TGeoRotation("", 90+2.*beta*TMath::RadToDeg(), |
579 | beamPhiPrime*TMath::RadToDeg(), -90); |
580 | TGeoRotation *beamRot4 = new TGeoRotation("", 90+2.*beta*TMath::RadToDeg(), |
581 | -beamPhiPrime*TMath::RadToDeg(),-90); |
b7943f00 |
582 | |
583 | TGeoCombiTrans *beamTransf[8]; |
584 | beamTransf[0] = new TGeoCombiTrans( 0.5*triangleHeight* |
585 | TMath::Tan(halfTheta), |
586 | fgkLadderBeamRadius/2. - dy, |
531d6cdc |
587 | -3*fgkSegmentLength/8, beamRot1); |
b7943f00 |
588 | beamTransf[1] = new TGeoCombiTrans(*beamTransf[0]); |
589 | AddTranslationToCombiTrans(beamTransf[1], 0, 0, fgkSegmentLength/2); |
590 | |
591 | beamTransf[2] = new TGeoCombiTrans(0.5*triangleHeight* |
592 | TMath::Tan(halfTheta), |
593 | fgkLadderBeamRadius/2. - dy, |
531d6cdc |
594 | -fgkSegmentLength/8, beamRot2); |
b7943f00 |
595 | beamTransf[3] = new TGeoCombiTrans(*beamTransf[2]); |
596 | AddTranslationToCombiTrans(beamTransf[3], 0, 0, fgkSegmentLength/2); |
597 | |
598 | beamTransf[4] = new TGeoCombiTrans(-0.5*triangleHeight* |
599 | TMath::Tan(halfTheta), |
600 | fgkLadderBeamRadius/2. - dy, |
531d6cdc |
601 | -3*fgkSegmentLength/8, beamRot3); |
b7943f00 |
602 | beamTransf[5] = new TGeoCombiTrans(*beamTransf[4]); |
603 | AddTranslationToCombiTrans(beamTransf[5], 0, 0, fgkSegmentLength/2); |
604 | |
605 | beamTransf[6] = new TGeoCombiTrans(-0.5*triangleHeight* |
531d6cdc |
606 | TMath::Tan(halfTheta),fgkLadderBeamRadius/2.-dy, -fgkSegmentLength/8,beamRot4); |
b7943f00 |
607 | beamTransf[7] = new TGeoCombiTrans(-0.5*triangleHeight* |
531d6cdc |
608 | TMath::Tan(halfTheta),fgkLadderBeamRadius/2.-dy,3*fgkSegmentLength/8,beamRot4); |
b7943f00 |
609 | |
610 | //--- Beams of the bottom |
611 | TGeoTubeSeg *bottomBeam1 = new TGeoTubeSeg(0, fgkLadderBeamRadius, |
612 | fgkLadderWidth/2.-fgkLadderLb/3, 0, 180); |
613 | TGeoVolume *bottomBeam1Vol = new TGeoVolume("ITSsddBottomBeam1Vol", |
614 | bottomBeam1, carbonFiberLadderStruct); |
615 | bottomBeam1Vol->SetLineColor(fColorCarbonFiber); |
616 | TGeoTubeSeg *bottomBeam2 = new TGeoTubeSeg(0, fgkLadderBeamRadius, |
617 | fgkLadderWidth/2.-fgkLadderLb/3, 0, 90); |
618 | TGeoVolume *bottomBeam2Vol = new TGeoVolume("ITSsddBottomBeam2Vol", |
619 | bottomBeam2, carbonFiberLadderStruct); |
620 | bottomBeam2Vol->SetLineColor(fColorCarbonFiber); |
621 | TGeoTubeSeg *bottomBeam3 = new TGeoTubeSeg(0, fgkLadderBeamRadius, |
622 | 0.5*fgkLadderWidth/SinD(fgkBottomBeamAngle) |
623 | - fgkLadderLb/3, 0, 180); |
624 | TGeoVolume *bottomBeam3Vol = new TGeoVolume("ITSsddBottomBeam3Vol", |
625 | bottomBeam3, carbonFiberLadderStruct); |
626 | bottomBeam3Vol->SetLineColor(fColorCarbonFiber); |
627 | //bottomBeam3Vol->SetLineColor(2); |
531d6cdc |
628 | TGeoRotation *bottomBeamRot1 = new TGeoRotation("", 90, 90, 90); |
629 | TGeoRotation *bottomBeamRot2 = new TGeoRotation("",-90, 90, -90); |
b7943f00 |
630 | |
b7943f00 |
631 | TGeoCombiTrans *bottomBeamTransf1 = new TGeoCombiTrans |
531d6cdc |
632 | (0,-(fgkLadderHeight/2-fgkLadderBeamRadius)-dy,0, bottomBeamRot1); |
b7943f00 |
633 | TGeoCombiTrans *bottomBeamTransf2 = new TGeoCombiTrans(0, |
634 | -(fgkLadderHeight/2-fgkLadderBeamRadius)-dy, |
531d6cdc |
635 | -fgkSegmentLength/2, bottomBeamRot1); |
b7943f00 |
636 | TGeoCombiTrans *bottomBeamTransf3 = new TGeoCombiTrans(0, |
637 | -(fgkLadderHeight/2 - fgkLadderBeamRadius) |
531d6cdc |
638 | - dy, fgkSegmentLength/2, bottomBeamRot2); |
b7943f00 |
639 | // be careful for beams #3: when "reading" from -z to +z and |
640 | // from the bottom of the ladder, it should draw a Lambda, and not a V |
531d6cdc |
641 | TGeoRotation *bottomBeamRot4 = new TGeoRotation("", -90, fgkBottomBeamAngle, -90); |
642 | TGeoRotation *bottomBeamRot5 = new TGeoRotation("" ,-90,-fgkBottomBeamAngle, -90); |
b7943f00 |
643 | TGeoCombiTrans *bottomBeamTransf4 = new TGeoCombiTrans |
531d6cdc |
644 | (0,-(fgkLadderHeight/2-fgkLadderBeamRadius)-dy,-fgkSegmentLength/4,bottomBeamRot4); |
b7943f00 |
645 | TGeoCombiTrans *bottomBeamTransf5 = new TGeoCombiTrans |
531d6cdc |
646 | (0,-(fgkLadderHeight/2-fgkLadderBeamRadius)-dy,fgkSegmentLength/4, bottomBeamRot5); |
b7943f00 |
647 | |
648 | fLaddSegCommonVol[0] = cfLaddTopVol1; fLaddSegCommonTr[0] = trTop1; |
649 | fLaddSegCommonVol[1] = cfLaddTopVol2; fLaddSegCommonTr[1] = trTop1; |
650 | fLaddSegCommonVol[2] = cfLaddSideVol1; fLaddSegCommonTr[2] = ctSideR; |
651 | fLaddSegCommonVol[3] = cfLaddSideVol1; fLaddSegCommonTr[3] = ctSideL; |
652 | fLaddSegCommonVol[4] = cfLaddSideVol2; fLaddSegCommonTr[4] = ctSideR; |
653 | fLaddSegCommonVol[5] = cfLaddSideVol2; fLaddSegCommonTr[5] = ctSideL; |
654 | fLaddSegCommonVol[6] = sideBeam; fLaddSegCommonTr[6] = beamTransf[0]; |
655 | fLaddSegCommonVol[7] = sideBeam; fLaddSegCommonTr[7] = beamTransf[1]; |
656 | fLaddSegCommonVol[8] = sideBeam; fLaddSegCommonTr[8] = beamTransf[2]; |
657 | fLaddSegCommonVol[9] = sideBeam; fLaddSegCommonTr[9] = beamTransf[3]; |
658 | fLaddSegCommonVol[10]= sideBeam; fLaddSegCommonTr[10]= beamTransf[4]; |
659 | fLaddSegCommonVol[11]= sideBeam; fLaddSegCommonTr[11]= beamTransf[5]; |
660 | fLaddSegCommonVol[12]= sideBeam; fLaddSegCommonTr[12]= beamTransf[6]; |
661 | fLaddSegCommonVol[13]= sideBeam; fLaddSegCommonTr[13]= beamTransf[7]; |
662 | fLaddSegCommonVol[14]= bottomBeam1Vol; fLaddSegCommonTr[14]= bottomBeamTransf1; |
663 | fLaddSegCommonVol[15]= bottomBeam2Vol; fLaddSegCommonTr[15]= bottomBeamTransf2; |
664 | fLaddSegCommonVol[16]= bottomBeam2Vol; fLaddSegCommonTr[16]= bottomBeamTransf3; |
665 | fLaddSegCommonVol[17]= bottomBeam3Vol; fLaddSegCommonTr[17]= bottomBeamTransf4; |
666 | fLaddSegCommonVol[18]= bottomBeam3Vol; fLaddSegCommonTr[18]= bottomBeamTransf5; |
c789ee28 |
667 | |
b7943f00 |
668 | |
669 | //******************************************************************** |
670 | // cables |
671 | //******************************************************************** |
672 | char cableName[30]; |
673 | for (Int_t i=0; i<fgkLay3Ndet; i++) { |
674 | sprintf(cableName, "digitCableLay3A_%i",i); |
675 | fDigitCableLay3A[i].SetName(cableName); |
676 | fDigitCableLay3A[i].SetWidth(fgkDigitCablWidth); |
677 | fDigitCableLay3A[i].SetThickness(fgkDigitCablPolyThick+fgkDigitCablAlThick); |
678 | fDigitCableLay3A[i].SetNLayers(2); |
679 | fDigitCableLay3A[i].SetLayer( 0, fgkDigitCablPolyThick, polyhamideSDD, |
680 | fColorPolyhamide); |
681 | fDigitCableLay3A[i].SetLayer(1, fgkDigitCablAlThick, alSDD, fColorAl); |
682 | sprintf(cableName, "digitCableLay3B_%i",i); |
683 | fDigitCableLay3B[i].SetName(cableName); |
684 | fDigitCableLay3B[i].SetWidth(fgkDigitCablWidth); |
685 | fDigitCableLay3B[i].SetThickness(fgkDigitCablPolyThick+fgkDigitCablAlThick); |
686 | fDigitCableLay3B[i].SetNLayers(2); |
687 | fDigitCableLay3B[i].SetLayer( 0, fgkDigitCablPolyThick, polyhamideSDD, |
688 | fColorPolyhamide); |
689 | fDigitCableLay3B[i].SetLayer(1, fgkDigitCablAlThick, alSDD, fColorAl); |
690 | }; |
691 | for (Int_t i=0; i<fgkLay4Ndet; i++) { |
692 | sprintf(cableName, "digitCableLay4A_%i",i); |
693 | fDigitCableLay4A[i].SetName(cableName); |
694 | fDigitCableLay4A[i].SetWidth(fgkDigitCablWidth); |
695 | fDigitCableLay4A[i].SetThickness(fgkDigitCablPolyThick+fgkDigitCablAlThick); |
696 | fDigitCableLay4A[i].SetNLayers(2); |
697 | fDigitCableLay4A[i].SetLayer( 0, fgkDigitCablPolyThick, polyhamideSDD, |
698 | fColorPolyhamide); |
699 | fDigitCableLay4A[i].SetLayer(1, fgkDigitCablAlThick, alSDD, fColorAl); |
700 | sprintf(cableName, "digitCableLay4B_%i",i); |
701 | fDigitCableLay4B[i].SetName(cableName); |
702 | fDigitCableLay4B[i].SetWidth(fgkDigitCablWidth); |
703 | fDigitCableLay4B[i].SetThickness(fgkDigitCablPolyThick+fgkDigitCablAlThick); |
704 | fDigitCableLay4B[i].SetNLayers(2); |
705 | fDigitCableLay4B[i].SetLayer( 0, fgkDigitCablPolyThick, polyhamideSDD, |
706 | fColorPolyhamide); |
707 | fDigitCableLay4B[i].SetLayer(1, fgkDigitCablAlThick, alSDD, fColorAl); |
708 | }; |
709 | // Well, those digit cables could also include the analog cables |
710 | // which have the same width and the same path, at least in the ladder. |
711 | // It will gain some computing ressources (less volumes) and some |
712 | // coding efforts ... !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
713 | // The only thing to do is to change the names and put the correct total |
714 | // thicknesses |
531d6cdc |
715 | } |
c789ee28 |
716 | |
c789ee28 |
717 | |
b7943f00 |
718 | //________________________________________________________________________ |
719 | void AliITSv11GeometrySDD::CheckOverlaps(Double_t precision){ |
720 | // |
721 | // a debugging function for checking some possible overlaps |
722 | // |
723 | if (fSDDsensor) fSDDsensor->CheckOverlaps(precision); |
724 | if (fHybrid) fHybrid->CheckOverlaps(precision); |
531d6cdc |
725 | } |
c789ee28 |
726 | |
727 | |
db486a6e |
728 | //________________________________________________________________________ |
729 | TGeoCombiTrans *AliITSv11GeometrySDD:: |
730 | CreateCombiTrans(const char *name, Double_t dy, Double_t dz, Double_t dphi) { |
731 | // |
732 | // return the TGeoCombiTrans which make a translation in y and z |
733 | // and a rotation in phi in the global coord system |
734 | // |
735 | |
736 | TGeoTranslation t1(dy*CosD(90.+dphi),dy*SinD(90.+dphi), dz); |
737 | TGeoRotation r1("",0.,0.,dphi); |
738 | |
739 | TGeoCombiTrans *combiTrans1 = new TGeoCombiTrans(name); |
740 | combiTrans1->SetTranslation(t1); |
741 | combiTrans1->SetRotation(r1); |
742 | return combiTrans1; |
531d6cdc |
743 | } |
c789ee28 |
744 | |
745 | |
db486a6e |
746 | //________________________________________________________________________ |
c789ee28 |
747 | void AliITSv11GeometrySDD::AddTranslationToCombiTrans(TGeoCombiTrans* ct, |
db486a6e |
748 | Double_t dx, |
749 | Double_t dy, |
b7943f00 |
750 | Double_t dz) const{ |
751 | // Add a dx,dy,dz translation to the initial TGeoCombiTrans |
752 | const Double_t *vect = ct->GetTranslation(); |
753 | Double_t newVect[3] = {vect[0]+dx, vect[1]+dy, vect[2]+dz}; |
754 | ct->SetTranslation(newVect); |
531d6cdc |
755 | } |
c789ee28 |
756 | |
757 | |
758 | //________________________________________________________________________ |
759 | void AliITSv11GeometrySDD::ShowOnePiece(TGeoVolume *moth) { |
760 | // for code developpment and debugging purposes |
761 | |
531d6cdc |
762 | if (! fSDDsensor) CreateBasicObjects(); |
763 | |
764 | // moth->AddNode(fPinSupport, 1, 0); |
765 | // moth->AddNode(fCoolPipeSupportL, 1, 0); |
766 | moth->AddNode(fSDDsensor, 1, 0); |
767 | // moth->AddNode(fBaseThermalBridge, 1, 0); |
768 | // moth->AddNode(fHybrid,100,0); |
b7943f00 |
769 | |
770 | |
531d6cdc |
771 | // TGeoVolume* seg = CreateLadderSegment( 3, 0); |
772 | // moth->AddNode(seg, 1, 0); |
b7943f00 |
773 | |
b7943f00 |
774 | |
531d6cdc |
775 | // TGeoVolume *lay3Ladder = CreateLadder(3); |
776 | // moth->AddNode(lay3Ladder, 1, 0); |
777 | |
778 | // TGeoVolume *lay3Detectors = CreateDetectors(3); |
779 | // moth->AddNode(lay3Detectors, 1, 0); |
780 | |
781 | } |
c789ee28 |
782 | |
783 | |
db486a6e |
784 | //________________________________________________________________________ |
b7943f00 |
785 | void AliITSv11GeometrySDD::Layer3(TGeoVolume *moth) { |
786 | // |
787 | // Insert the layer 3 in the mother volume. This is a virtual volume |
788 | // containing ladders of layer 3 and the supporting rings |
789 | // |
790 | |
791 | if (! moth) { |
792 | printf("Error::AliITSv11GeometrySDD: Can't insert layer3, mother is null!\n"); |
793 | return; |
794 | }; |
795 | |
796 | TGeoMedium *airSDD = GetMedium("ITSair"); |
797 | |
798 | fMotherVol = moth; |
799 | if (! fSDDsensor) CreateBasicObjects(); |
800 | |
801 | TGeoVolume *lay3Ladder = CreateLadder(3); |
802 | TGeoVolume *lay3Detectors = CreateDetectors(3); |
803 | TGeoTube *virtualLayer3Shape = new TGeoTube("ITSsddLayer3Shape", |
804 | fgkLay3Rmin,fgkLay3Rmax,fgkLay3Length*0.5); |
805 | TGeoVolume *virtualLayer3 = new TGeoVolume("ITSsddLayer3", |
806 | virtualLayer3Shape, airSDD); |
807 | |
808 | Double_t dPhi = 360./fgkLay3Nladd; |
809 | Double_t detectorsThick = fgkLadWaferSep + 2*fgkWaferThickness; |
810 | // Placing virtual ladder and detectors volumes following |
811 | // ladder ordering convention |
812 | char rotName[30]; |
813 | Int_t iLaddMin = 0; |
814 | Int_t iLaddMax = fgkLay3Nladd; |
815 | if ((fAddOnlyLadder3min>=0)&&(fAddOnlyLadder3max<fgkLay3Nladd)) { |
816 | iLaddMin = fAddOnlyLadder3min; |
817 | iLaddMax = fAddOnlyLadder3max+1; |
818 | }; |
819 | |
820 | for (Int_t iLadd = iLaddMin; iLadd < iLaddMax; iLadd++) { |
821 | |
822 | sprintf(rotName, "ITSsddLay3Ladd%i",iLadd); |
823 | Double_t minRadiusLadBox = fLay3LaddShortRadius-fLay3LadderUnderSegDH; |
824 | if (iLadd%2 != 0) |
825 | minRadiusLadBox = fLay3LaddLongRadius-fLay3LadderUnderSegDH; |
826 | minRadiusLadBox += ((TGeoBBox*)lay3Ladder->GetShape())->GetDY(); |
827 | TGeoCombiTrans *ctLadd = CreateCombiTrans(rotName,minRadiusLadBox, |
828 | 0,-90+iLadd*dPhi); |
829 | virtualLayer3->AddNode(lay3Ladder, iLadd, ctLadd); |
830 | /////////////////////////////////////////////////// |
831 | sprintf(rotName, "ITSsddLay3DetBox%i",iLadd); |
832 | Double_t minRadiusDetBox = fgkLay3DetShortRadius; |
833 | if (iLadd%2 != 0) minRadiusDetBox = fgkLay3DetLongRadius; |
834 | minRadiusDetBox += detectorsThick/2; |
835 | TGeoCombiTrans *ctDet = CreateCombiTrans(rotName, minRadiusDetBox, |
836 | 0,-90+iLadd*dPhi); |
837 | virtualLayer3->AddNode(lay3Detectors, iLadd, ctDet); |
838 | /////////////////////////////////////////////////// |
839 | } |
840 | |
841 | if(GetDebug(1)) virtualLayer3->CheckOverlaps(0.01); |
842 | virtualLayer3->SetVisibility(kFALSE); |
843 | moth->AddNode(virtualLayer3, 1, 0); |
531d6cdc |
844 | } |
c789ee28 |
845 | |
846 | |
db486a6e |
847 | //________________________________________________________________________ |
b7943f00 |
848 | void AliITSv11GeometrySDD::Layer4(TGeoVolume *moth) { |
849 | // |
850 | // Insert the layer 4 in the mother volume. This is a virtual volume |
851 | // containing ladders of layer 4 and the supporting rings |
852 | // |
853 | |
854 | if (! moth) { |
855 | printf("Error::AliITSv11GeometrySDD: Can't insert layer4, mother is null!\n"); |
856 | return; |
857 | }; |
858 | |
859 | fMotherVol = moth; |
860 | |
861 | if (! fSDDsensor) CreateBasicObjects(); |
862 | |
863 | TGeoTube *virtualLayer4Shape =new TGeoTube("ITSsddLayer4Shape", |
864 | fgkLay4Rmin,fgkLay4Rmax,fgkLay4Length*0.5); |
865 | TGeoMedium *airSDD = GetMedium("ITSair"); |
866 | TGeoVolume *virtualLayer4 = new TGeoVolume("ITSsddLayer4", |
867 | virtualLayer4Shape, airSDD); |
868 | TGeoVolume *lay4Ladder = CreateLadder(4); |
869 | TGeoVolume *lay4Detectors = CreateDetectors(4); |
870 | Double_t dPhi = 360./fgkLay4Nladd; |
871 | Double_t detBoxThickness = fgkLadWaferSep + 2*fgkWaferThickness; |
872 | // placing virtual ladder and detectors volumes following ladder |
873 | // ordering convention |
874 | char rotName[20]; |
875 | Int_t iLaddMin = 0; |
876 | Int_t iLaddMax = fgkLay4Nladd; |
877 | if ((fAddOnlyLadder4min >= 0)&&(fAddOnlyLadder4max < fgkLay4Nladd)) { |
878 | iLaddMin = fAddOnlyLadder4min; |
879 | iLaddMax = fAddOnlyLadder4max+1; |
880 | } |
881 | for (Int_t iLadd = iLaddMin; iLadd < iLaddMax; iLadd++) { |
882 | sprintf(rotName, "ITSsddLay4Ladd%i",iLadd); |
883 | Double_t minRadiusLadBox = fLay4LaddShortRadius-fLay4LadderUnderSegDH; |
884 | if (iLadd%2 != 0) |
885 | minRadiusLadBox = fLay4LaddLongRadius-fLay4LadderUnderSegDH; |
886 | minRadiusLadBox += ((TGeoBBox*)lay4Ladder->GetShape())->GetDY(); |
887 | TGeoCombiTrans *ctLadd = CreateCombiTrans(rotName, minRadiusLadBox, |
888 | 0, -90+iLadd*dPhi); |
889 | virtualLayer4->AddNode(lay4Ladder, iLadd, ctLadd); |
890 | sprintf(rotName, "ITSsddLay4DetBox%i",iLadd); |
891 | Double_t minRadiusDetBox = fgkLay4DetShortRadius; |
892 | if (iLadd%2 != 0) |
893 | minRadiusDetBox = fgkLay4DetLongRadius; |
894 | minRadiusDetBox += detBoxThickness/2; |
895 | TGeoCombiTrans *ctDet = CreateCombiTrans(rotName, minRadiusDetBox, |
896 | 0, -90+iLadd*dPhi); |
897 | virtualLayer4->AddNode(lay4Detectors, iLadd, ctDet); |
898 | } |
899 | |
900 | if(GetDebug(1)) virtualLayer4->CheckOverlaps(0.01); |
901 | virtualLayer4->SetVisibility(kFALSE); |
902 | moth->AddNode(virtualLayer4,1,0); |
531d6cdc |
903 | } |
b7943f00 |
904 | |
905 | |
906 | //________________________________________________________________________ |
907 | TGeoVolume *AliITSv11GeometrySDD::CreateLadder(Int_t iLay) { |
908 | // |
909 | // return a box volume containing the CF ladder |
910 | // |
911 | |
912 | Int_t nDetectors = fgkLay3Ndet; |
913 | Double_t ladderLength = fgkLay3LadderLength; |
914 | Double_t underSegDH = fLay3LadderUnderSegDH; |
915 | Double_t *sensorZPos = fLay3sensorZPos; |
916 | AliITSv11GeomCableFlat *digitCableA = fDigitCableLay3A; |
917 | AliITSv11GeomCableFlat *digitCableB = fDigitCableLay3B; |
918 | |
919 | if (iLay==3) {} |
920 | else if (iLay==4) { |
921 | nDetectors = fgkLay4Ndet; |
922 | ladderLength = fgkLay4LadderLength; |
923 | digitCableA = fDigitCableLay4A; |
924 | digitCableB = fDigitCableLay4B; |
925 | underSegDH = fLay4LadderUnderSegDH; |
926 | sensorZPos = fLay4sensorZPos; |
927 | } |
928 | else { |
929 | printf("AliITSv11GeometrySDD::CreateLadder : error=wrong layer\n"); |
930 | }; |
931 | Double_t ladderBoxDH = fgkLadderHeight+fgkLadderSegBoxDH+underSegDH; |
932 | TGeoBBox *ladBox = new TGeoBBox("ITSsddLadBox", |
933 | fgkLadderWidth/2+fgkPinSuppWidth+fgkLadderSegBoxDW, |
934 | ladderBoxDH/2, ladderLength/2); |
935 | TGeoMedium *airSDD = GetMedium("ITSair"); |
936 | TGeoVolume *virtualLadder = new TGeoVolume("ITSsddLadder",ladBox, airSDD); |
937 | |
938 | // placing virtual ladder segment following detector ordering convention |
939 | //======================================================================= |
940 | char transName[30]; |
941 | |
942 | // adding segment this way to create cable points in the correct order ... |
943 | for (Int_t iSegment = nDetectors/2-1; iSegment >= 0; iSegment-- ) { |
944 | |
945 | TGeoVolume *laddSegment = CreateLadderSegment(iLay, iSegment); |
946 | sprintf(transName, "ITSsddLay%iLaddSeg%i", iLay, iSegment); |
947 | Double_t segmentPos = fgkSegmentLength*(nDetectors/2-1-iSegment) |
948 | + fgkSegmentLength/2; |
949 | TGeoTranslation *segTr = new TGeoTranslation(transName, 0, |
950 | underSegDH/2,segmentPos); |
951 | //// |
952 | virtualLadder->AddNode(laddSegment, iSegment, segTr); |
953 | }; |
954 | for (Int_t iSegment = nDetectors/2; iSegment < nDetectors; iSegment++ ) { |
955 | |
956 | TGeoVolume *laddSegment = CreateLadderSegment(iLay, iSegment); |
957 | sprintf(transName, "ITSsddLay%iLaddSeg%i", iLay, iSegment); |
958 | Double_t segmentPos = fgkSegmentLength*(nDetectors/2-1-iSegment) |
959 | + fgkSegmentLength/2; |
960 | TGeoTranslation *segTr = new TGeoTranslation(transName, 0, |
961 | underSegDH/2,segmentPos); |
962 | //// |
963 | virtualLadder->AddNode(laddSegment, iSegment, segTr); |
964 | }; |
965 | |
966 | // putting virtual volume corresponding to the end of ladder |
967 | //======================================================================= |
968 | TGeoVolume *endLadder = CreateEndLadder( iLay ); |
969 | Double_t endLength = (ladderLength - nDetectors*fgkSegmentLength)/2.; |
970 | TGeoTranslation *endTrZPos = new TGeoTranslation("ITSsddEndTrZPos",0,0, |
971 | fgkSegmentLength*(nDetectors/2)+endLength/2.); |
972 | // Euler rotation : about Z, then new X, then new Z |
973 | TGeoRotation *endZNegRot = new TGeoRotation("",90, 180, -90); |
974 | TGeoCombiTrans *endTrZNeg = new TGeoCombiTrans(0,0, |
975 | -fgkSegmentLength*(nDetectors/2)-endLength/2.,endZNegRot); |
976 | virtualLadder->AddNode(endLadder, 1, endTrZPos); |
977 | virtualLadder->AddNode(endLadder, 2, endTrZNeg); |
978 | |
979 | // creating and inserting cable segments |
980 | // (check points are placed while creating segments) |
981 | //======================================================================= |
982 | if (fAddCables) |
983 | for (Int_t iSegment = 0; iSegment < nDetectors; iSegment++ ) { |
984 | |
985 | digitCableA[iSegment].SetInitialNode(virtualLadder); |
986 | digitCableB[iSegment].SetInitialNode(virtualLadder); |
987 | |
988 | for (Int_t iPt=1; iPt<digitCableA[iSegment].GetNCheckPoints(); iPt++ ) { |
989 | Double_t rotation = 0; |
990 | if (iPt>1) rotation = 90-fgkHybridAngle; |
991 | digitCableA[iSegment].CreateAndInsertCableSegment(iPt, rotation); |
992 | }; |
993 | |
994 | for (Int_t iPt=1; iPt<digitCableB[iSegment].GetNCheckPoints(); iPt++ ) { |
995 | Double_t rotation = 0; |
996 | if (iPt>1) rotation = fgkHybridAngle-90; |
997 | digitCableB[iSegment].CreateAndInsertCableSegment(iPt, rotation); |
998 | }; |
999 | }; |
1000 | |
1001 | // HV cable |
1002 | //======================================================================= |
fa4639a3 |
1003 | if (fAddHVcables) { |
b7943f00 |
1004 | TGeoMedium *polyhamideSDD = GetMedium("ITSsddKAPTON_POLYCH2"); |
1005 | TGeoMedium *alSDD = GetMedium("ITSal"); |
1006 | |
1007 | AliITSv11GeomCableFlat cableHV[fgkLay4Ndet]; // temp !!! |
1008 | char cableHVname[30]; |
1009 | for (Int_t iSegment = 0; iSegment<nDetectors; iSegment++) { |
1010 | sprintf(cableHVname,"ITSsddHVcable%i", iSegment); |
1011 | cableHV[iSegment].SetName(cableHVname); |
1012 | cableHV[iSegment].SetThickness(fgkLongHVcablePolyThick+fgkLongHVcableAlThick); |
1013 | cableHV[iSegment].SetWidth(fgkTransitHVtailWidth); |
1014 | cableHV[iSegment].SetNLayers(2); |
1015 | cableHV[iSegment].SetLayer(0, fgkLongHVcablePolyThick, polyhamideSDD, |
1016 | fColorPolyhamide); |
1017 | cableHV[iSegment].SetLayer(1, fgkLongHVcableAlThick, alSDD, fColorAl); |
1018 | cableHV[iSegment].SetInitialNode(virtualLadder); |
1019 | }; |
1020 | Double_t x1[3], x2[3], x3[3], |
1021 | vY[3] = {0,1,0}, vZ[3] = {0,0,1}, vYZ[3]={0,1,1}; |
1022 | |
fa4639a3 |
1023 | x1[0] = -fgkTransitHVtailXpos; |
1024 | x2[0] = -fgkTransitHVtailXpos; |
1025 | x3[0] = -fgkTransitHVtailXpos; |
b7943f00 |
1026 | for (Int_t iSegment = nDetectors/2-1; iSegment >= 0; iSegment-- ) { |
1027 | Double_t cableSeparation = TMath::Abs(iSegment - (nDetectors/2-1)) |
1028 | *fgkLongHVcableSeparation; |
1029 | x1[1] = - ladderBoxDH/2; |
1030 | x2[1] = - ladderBoxDH/2 + underSegDH - cableSeparation |
1031 | - (fgkLongHVcablePolyThick+fgkLongHVcableAlThick)/2; |
1032 | x3[1] = x2[1]; |
1033 | x1[2] = sensorZPos[iSegment]+fgkTransitHVtailLength-5*fgkmm; |
1034 | x2[2] = x1[2]+5*fgkmm; |
1035 | x3[2] = ladderLength/2-endLength; |
1036 | cableHV[iSegment].AddCheckPoint( virtualLadder, 0, x1, vY ); |
1037 | cableHV[iSegment].AddCheckPoint( virtualLadder, 1, x2, vYZ ); |
1038 | cableHV[iSegment].AddCheckPoint( virtualLadder, 2, x3, vZ ); |
1039 | |
1040 | cableHV[iSegment].CreateAndInsertCableSegment(1,0); |
1041 | cableHV[iSegment].CreateAndInsertCableSegment(2,0); |
1042 | }; |
1043 | |
1044 | vYZ[2] = -1; |
fa4639a3 |
1045 | x1[0] = fgkTransitHVtailXpos; |
1046 | x2[0] = fgkTransitHVtailXpos; |
1047 | x3[0] = fgkTransitHVtailXpos; |
1048 | |
b7943f00 |
1049 | for (Int_t iSegment = nDetectors/2; iSegment < nDetectors; iSegment++ ) { |
1050 | Double_t cableSeparation = TMath::Abs(iSegment - (nDetectors/2-1)) |
1051 | *fgkLongHVcableSeparation; |
1052 | x1[1] = - ladderBoxDH/2; |
1053 | x2[1] = - ladderBoxDH/2 + underSegDH - cableSeparation |
1054 | - (fgkLongHVcablePolyThick+fgkLongHVcableAlThick)/2; |
1055 | x3[1] = x2[1]; |
1056 | x1[2] = sensorZPos[iSegment]-fgkTransitHVtailLength+5*fgkmm; |
1057 | x2[2] = x1[2]-5*fgkmm; |
1058 | x3[2] = -ladderLength/2+endLength; |
1059 | cableHV[iSegment].AddCheckPoint( virtualLadder, 0, x1, vY ); |
1060 | cableHV[iSegment].AddCheckPoint( virtualLadder, 1, x2, vYZ ); |
1061 | cableHV[iSegment].AddCheckPoint( virtualLadder, 2, x3, vZ ); |
1062 | |
1063 | cableHV[iSegment].CreateAndInsertCableSegment(1,0); |
1064 | cableHV[iSegment].CreateAndInsertCableSegment(2,0); |
1065 | }; |
fa4639a3 |
1066 | }; |
c789ee28 |
1067 | |
b7943f00 |
1068 | //********************************** |
1069 | if(GetDebug(1)) virtualLadder->CheckOverlaps(0.01); |
fa4639a3 |
1070 | virtualLadder->SetVisibility(kFALSE); |
b7943f00 |
1071 | return virtualLadder; |
531d6cdc |
1072 | } |
c789ee28 |
1073 | |
1074 | |
db486a6e |
1075 | //________________________________________________________________________ |
b7943f00 |
1076 | TGeoArb8 *AliITSv11GeometrySDD::CreateLadderSide(Double_t dz, Double_t angle, |
1077 | Double_t xSign, Double_t L, Double_t H, Double_t l) { |
db486a6e |
1078 | // Create one half of the V shape corner of CF ladder |
1079 | |
1080 | TGeoArb8 *cfLaddSide = new TGeoArb8(dz); |
c789ee28 |
1081 | cfLaddSide->SetVertex( 0, 0, 0); |
1082 | cfLaddSide->SetVertex( 1, 0, -H); |
1083 | cfLaddSide->SetVertex( 2, xSign*(L*TMath::Sin(angle)-l*TMath::Cos(angle)), |
1084 | -L*TMath::Cos(angle)-l*TMath::Sin(angle)); |
1085 | cfLaddSide->SetVertex( 3, xSign*L*TMath::Sin(angle), -L*TMath::Cos(angle)); |
1086 | cfLaddSide->SetVertex( 4, 0, 0); |
1087 | cfLaddSide->SetVertex( 5, 0, -H); |
1088 | cfLaddSide->SetVertex( 6, xSign*(L*TMath::Sin(angle)-l*TMath::Cos(angle)), |
1089 | -L*TMath::Cos(angle)-l*TMath::Sin(angle)); |
1090 | cfLaddSide->SetVertex(7, xSign*L*TMath::Sin(angle), -L*TMath::Cos(angle)); |
db486a6e |
1091 | return cfLaddSide; |
531d6cdc |
1092 | } |
c789ee28 |
1093 | |
1094 | |
db486a6e |
1095 | //________________________________________________________________________ |
b7943f00 |
1096 | TGeoVolume* AliITSv11GeometrySDD::CreateHybrid(Int_t iLRSide) { |
1097 | // |
1098 | // return a box containing the front-end hybrid |
1099 | // |
1100 | |
1101 | Double_t roundHoleX = -fgkHybridWidth/2+fgkHybRndHoleX; |
1102 | |
1103 | Double_t screenTotalThick = fgkHybGlueScrnThick+fgkHybUpThick+fgkHybAlThick; |
1104 | Double_t lowFLTotalThick = fgkHybGlueLowThick+fgkHybUpThick+fgkHybAlThick; |
1105 | // Double_t upFLTotalThick = fgkHybGlueUpThick +fgkHybUpThick+fgkHybAlThick; |
1106 | Double_t chipsCCTotThick = fgkHybUnderNiThick+fgkHybGlueAgThick |
1107 | +fgkHybChipThick+2*(fgkHybUpCCThick+fgkHybAlCCThick); |
1108 | Double_t ccUpLayerTotThick = fgkHybUpCCThick+fgkHybAlCCThick+fgkHybUpCCThick; |
1109 | // Double_t volumeThick = (fgkHybridThBridgeThick+screenTotalThick+lowFLTotalThick |
1110 | // + upFLTotalThick + ccUpLayerTotThick); |
1111 | Double_t volumeThick = (fgkHybridThBridgeThick+screenTotalThick+lowFLTotalThick |
1112 | +fgkHybSMDheight); |
1113 | Double_t lowLayerYmin = -volumeThick/2+fgkHybridThBridgeThick |
1114 | +screenTotalThick; |
1115 | Double_t flUpThick = fgkHybGlueUpThick+fgkHybUpThick; |
1116 | |
1117 | //**************************************************** media : |
1118 | TGeoMedium *airSDD = GetMedium("ITSair"); |
1119 | TGeoMedium *carbonFiberLadderStruct = GetMedium("ITSsddCarbonM55J"); |
1120 | TGeoMedium *alSDD = GetMedium("ITSal"); |
1121 | TGeoMedium *alSDD80p100 = GetMedium("ITSal"); // to code !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
1122 | TGeoMedium *alSDD50p100 = GetMedium("ITSal"); // to code !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
1123 | TGeoMedium *polyhamideSDD = GetMedium("ITSsddKAPTON_POLYCH2"); |
1124 | TGeoMedium *niSDD = GetMedium("COPPER"); // to code !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
1125 | TGeoMedium *glueAG = GetMedium("ITSsddKAPTON_POLYCH2"); // to code !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
1126 | TGeoMedium *siliconSDD = GetMedium("ITSsddSiChip"); |
1127 | TGeoMedium *medSMD = GetMedium("SDDX7Rcapacitors"); // TO CODE !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
1128 | TGeoMedium *medSMDweld = GetMedium("SDDX7Rcapacitors"); // TO CODE !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
1129 | |
1130 | //**************************************************** main volume : |
1131 | TGeoBBox *hybridBox = new TGeoBBox("",fgkHybridWidth/2, volumeThick/2, |
1132 | (fgkHybridLength)/2); |
1133 | TGeoVolume *hybrid = new TGeoVolume("ITSsddHybridVol", hybridBox, |
1134 | airSDD); |
1135 | |
1136 | TGeoBBox *sThermalBridge = new TGeoBBox( "", fgkHybridWidth/2, |
1137 | fgkHybridThBridgeThick/2, |
1138 | fgkHybridLength/2); |
c789ee28 |
1139 | |
b7943f00 |
1140 | //**************************************************** Thermal bridge : |
c789ee28 |
1141 | TGeoVolume *vThermalBridge = new TGeoVolume("ITSsddHybridThBridge", |
1142 | sThermalBridge, |
1143 | carbonFiberLadderStruct); |
1144 | vThermalBridge->SetLineColor(fColorCarbonFiber); |
b7943f00 |
1145 | TGeoTranslation *thBridgeTr = new TGeoTranslation(0, -volumeThick/2 |
1146 | +fgkHybridThBridgeThick/2, 0); |
1147 | hybrid->AddNode(vThermalBridge, 1, thBridgeTr); |
1148 | |
1149 | //**************************************************** Screen layer : |
1150 | TGeoBBox *sAlScreenLayer = new TGeoBBox("sAlScreenLayer", fgkHybridWidth/2, |
1151 | fgkHybAlThick/2, fgkHybridLength/2); |
1152 | //here the upedex and glue layers are both assumed to be polyimide |
1153 | TGeoBBox *sUpGlueScreenLayer = new TGeoBBox("sUpGlueScreenLayer", |
1154 | fgkHybridWidth/2, |
1155 | (fgkHybUpThick+fgkHybGlueScrnThick)/2, |
1156 | fgkHybridLength/2); |
1157 | TGeoTube *sRoundHole = new TGeoTube("sRoundHole", 0, fgkHybRndHoleRad, |
1158 | (screenTotalThick+lowFLTotalThick)/2); |
1159 | |
1160 | TGeoTranslation *upGlueScreenTr = new TGeoTranslation("upGlueScreenTr",0, |
1161 | -volumeThick/2+fgkHybridThBridgeThick+(fgkHybUpThick+fgkHybGlueScrnThick)/2,0); |
1162 | |
1163 | TGeoTranslation *alScreenTr = new TGeoTranslation("AlScreenTr", 0, |
1164 | -volumeThick/2+fgkHybridThBridgeThick+fgkHybUpThick+fgkHybGlueScrnThick |
1165 | +fgkHybAlThick/2, 0); |
1166 | |
1167 | TGeoTranslation hybHolePos1Tr(roundHoleX, |
1168 | -volumeThick/2+fgkHybridThBridgeThick+(screenTotalThick+lowFLTotalThick)/2, |
1169 | -fgkHybridLength/2+fgkHybRndHoleZ); |
1170 | TGeoTranslation hybHolePos2Tr(roundHoleX, |
1171 | -volumeThick/2+fgkHybridThBridgeThick+(screenTotalThick+lowFLTotalThick)/2, |
1172 | fgkHybridLength/2-fgkHybRndHoleZ); |
531d6cdc |
1173 | |
1174 | TGeoRotation *rotHole = new TGeoRotation("", 0, 90, 0); |
1175 | TGeoCombiTrans *hybHolePos1 = new TGeoCombiTrans(hybHolePos1Tr, *rotHole); |
b7943f00 |
1176 | hybHolePos1->SetName("hybHolePos1"); |
531d6cdc |
1177 | TGeoCombiTrans *hybHolePos2 = new TGeoCombiTrans(hybHolePos2Tr, *rotHole); |
b7943f00 |
1178 | hybHolePos2->SetName("hybHolePos2"); |
1179 | |
1180 | upGlueScreenTr->RegisterYourself(); |
1181 | alScreenTr->RegisterYourself(); |
1182 | hybHolePos1->RegisterYourself(); |
1183 | hybHolePos2->RegisterYourself(); |
531d6cdc |
1184 | delete rotHole; |
b7943f00 |
1185 | |
1186 | TGeoCompositeShape *sScreenAl = new TGeoCompositeShape( |
1187 | "sAlScreenLayer:AlScreenTr-(sRoundHole:hybHolePos1" |
1188 | "+sRoundHole:hybHolePos2)"); |
1189 | TGeoVolume *vScreenAl = new TGeoVolume("vScreenAl",sScreenAl, alSDD); |
1190 | vScreenAl->SetLineColor(fColorAl); |
1191 | TGeoCompositeShape *sScreenUpGlue = new TGeoCompositeShape( |
1192 | "sUpGlueScreenLayer:upGlueScreenTr-(sRoundHole:hybHolePos1" |
1193 | "+sRoundHole:hybHolePos2)"); |
1194 | TGeoVolume *vScreenUpGlue = new TGeoVolume("vScreenUpGlue", |
1195 | sScreenUpGlue,polyhamideSDD); |
1196 | vScreenUpGlue->SetLineColor(fColorPolyhamide); |
1197 | |
1198 | hybrid->AddNode(vScreenUpGlue, 1, 0); |
1199 | hybrid->AddNode(vScreenAl, 1, 0); |
1200 | |
1201 | //**************************************************** FL low layer : |
1202 | Double_t sideWidth1 = fgkHybFLlowChipZ1 - fgkHybFLlowHoleDZ/2; |
1203 | Double_t sideWidth2 = fgkHybridLength - fgkHybFLlowChipZ4 - fgkHybFLlowHoleDZ/2; |
1204 | |
1205 | //here the upedex and glue layers are both assumed to be polyimide |
1206 | TGeoBBox *sUpGlueBar1 = new TGeoBBox("sUpGlueBar1", fgkHybridWidth/2, |
1207 | (fgkHybGlueLowThick+fgkHybUpThick)/2, |
1208 | sideWidth1/2); |
1209 | TGeoBBox *sAlBar1 = new TGeoBBox("sAlBar1", fgkHybridWidth/2, |
1210 | fgkHybAlThick/2, sideWidth1/2); |
1211 | |
1212 | TGeoTranslation *upGlueBarTr1 = new TGeoTranslation("upGlueBarTr1", 0, |
1213 | lowLayerYmin+(fgkHybGlueLowThick+fgkHybUpThick)/2, |
1214 | -(fgkHybridLength-sideWidth1)/2); |
1215 | TGeoTranslation *alBarTr1 = new TGeoTranslation("alBarTr1", 0, |
1216 | lowLayerYmin+fgkHybGlueLowThick+fgkHybUpThick+fgkHybAlThick/2, |
1217 | -(fgkHybridLength-sideWidth1)/2); |
1218 | upGlueBarTr1->RegisterYourself(); |
1219 | alBarTr1->RegisterYourself(); |
1220 | |
1221 | TGeoCompositeShape *sLowUpGlueBar1 = new TGeoCompositeShape( |
1222 | "sUpGlueBar1:upGlueBarTr1-sRoundHole:hybHolePos1"); |
1223 | TGeoCompositeShape *sLowAlBar1 = new TGeoCompositeShape( |
1224 | "sAlBar1:alBarTr1-sRoundHole:hybHolePos1"); |
1225 | TGeoVolume *vLowUpGlueBar1 = new TGeoVolume("vLowUpGlueBar1", |
1226 | sLowUpGlueBar1, polyhamideSDD); |
1227 | TGeoVolume *vLowAlBar1 = new TGeoVolume("vLowAlBar1", |
1228 | sLowAlBar1, alSDD); |
1229 | vLowUpGlueBar1->SetLineColor(fColorPolyhamide); |
1230 | vLowAlBar1->SetLineColor(fColorAl); |
1231 | hybrid->AddNode(vLowUpGlueBar1,1,0); |
1232 | hybrid->AddNode(vLowAlBar1,1,0); |
1233 | |
1234 | //--- |
1235 | //here the upedex and glue layers are both assumed to be polyimide |
1236 | TGeoBBox *sUpGlueBar2 = new TGeoBBox("sUpGlueBar2", fgkHybridWidth/2, |
1237 | (fgkHybGlueLowThick+fgkHybUpThick)/2, |
1238 | sideWidth2/2); |
1239 | TGeoBBox *sAlBar2 = new TGeoBBox("sAlBar2", fgkHybridWidth/2, |
1240 | fgkHybAlThick/2, sideWidth2/2); |
1241 | |
1242 | TGeoTranslation *upGlueBarTr2 = new TGeoTranslation("upGlueBarTr2", 0, |
1243 | lowLayerYmin+(fgkHybGlueLowThick+fgkHybUpThick)/2, |
1244 | (fgkHybridLength-sideWidth2)/2); |
1245 | TGeoTranslation *alBarTr2 = new TGeoTranslation("alBarTr2", 0, |
1246 | lowLayerYmin+fgkHybGlueLowThick+fgkHybUpThick+fgkHybAlThick/2, |
1247 | (fgkHybridLength-sideWidth2)/2); |
1248 | upGlueBarTr2->RegisterYourself(); |
1249 | alBarTr2->RegisterYourself(); |
1250 | |
1251 | TGeoCompositeShape *sLowUpGlueBar2 = new TGeoCompositeShape( |
1252 | "sUpGlueBar2:upGlueBarTr2-sRoundHole:hybHolePos2"); |
1253 | TGeoCompositeShape *sLowAlBar2 = new TGeoCompositeShape( |
1254 | "sAlBar2:alBarTr2-sRoundHole:hybHolePos2"); |
1255 | TGeoVolume *vLowUpGlueBar2 = new TGeoVolume("vLowUpGlueBar2",sLowUpGlueBar2, |
1256 | polyhamideSDD); |
1257 | TGeoVolume *vLowAlBar2 = new TGeoVolume("vLowAlBar2",sLowAlBar2, |
1258 | alSDD); |
1259 | vLowUpGlueBar2->SetLineColor(fColorPolyhamide); |
1260 | vLowAlBar2->SetLineColor(fColorAl); |
1261 | hybrid->AddNode(vLowUpGlueBar2, 1, 0); |
1262 | hybrid->AddNode(vLowAlBar2, 1, 0); |
1263 | |
1264 | if(GetDebug(3)){ // Remove compiler warning. |
1265 | sAlScreenLayer->InspectShape(); |
1266 | sUpGlueScreenLayer->InspectShape(); |
1267 | sRoundHole->InspectShape(); |
1268 | sUpGlueBar1->InspectShape(); |
1269 | sUpGlueBar2->InspectShape(); |
1270 | sAlBar1->InspectShape(); |
1271 | sAlBar2->InspectShape(); |
1272 | }; |
1273 | //--- |
1274 | //using class AliITSv11GeomCableFlat to add 2-layer segments ... |
1275 | Double_t piece1width = fgkHybFLlowPasX-fgkHybFLlowHolePasDX/2; |
1276 | AliITSv11GeomCableFlat lowFLpiece("lowFLpiece1",piece1width, |
1277 | lowFLTotalThick); |
1278 | lowFLpiece.SetNLayers(2); |
1279 | lowFLpiece.SetLayer(0, fgkHybGlueLowThick+fgkHybUpThick, polyhamideSDD, |
1280 | fColorPolyhamide); |
1281 | lowFLpiece.SetLayer(1, fgkHybAlThick, alSDD80p100, fColorAl); |
1282 | // alSDD at 80% : mostly to take into account strips of piece 3 |
1283 | |
1284 | Double_t x1[3] = { -fgkHybridWidth/2 + piece1width/2, |
1285 | lowLayerYmin + lowFLTotalThick/2, |
1286 | -fgkHybridLength/2 + sideWidth1 }; |
1287 | Double_t x2[3] ={ x1[0], x1[1], fgkHybridLength/2 - sideWidth2 }; |
1288 | Double_t vZ[3] = {0,0,1}; |
1289 | lowFLpiece.AddCheckPoint( hybrid, 0, x2, vZ ); |
1290 | lowFLpiece.AddCheckPoint( hybrid, 1, x1, vZ ); |
1291 | lowFLpiece.SetInitialNode(hybrid); |
1292 | lowFLpiece.CreateAndInsertCableSegment(1); |
1293 | lowFLpiece.ResetPoints(); |
1294 | |
1295 | Double_t piece2width = fgkHybFLlowAmbX-fgkHybFLlowPasX |
1296 | -fgkHybFLlowHolePasDX/2-fgkHybFLlowHoleAmbDX/2; |
1297 | |
1298 | lowFLpiece.SetWidth(piece2width); |
1299 | lowFLpiece.SetName("lowFLpiece2"); |
1300 | x1[0] = piece2width/2+fgkHybFLlowPasX+fgkHybFLlowHolePasDX/2-fgkHybridWidth/2; |
1301 | x2[0] = x1[0]; |
1302 | lowFLpiece.AddCheckPoint( hybrid, 0, x2, vZ ); |
1303 | lowFLpiece.AddCheckPoint( hybrid, 1, x1, vZ ); |
1304 | lowFLpiece.CreateAndInsertCableSegment(1); |
1305 | lowFLpiece.ResetPoints(); |
1306 | |
1307 | Double_t piece3width = fgkHybridWidth - fgkHybFLlowAmbX |
1308 | - fgkHybFLlowHoleAmbDX/2; |
1309 | |
1310 | lowFLpiece.SetWidth(piece3width); |
1311 | lowFLpiece.SetName("lowFLpiece3"); |
1312 | x1[0] = fgkHybridWidth/2-piece3width/2; |
1313 | x2[0] = x1[0]; |
1314 | lowFLpiece.AddCheckPoint( hybrid, 0, x2, vZ ); |
1315 | lowFLpiece.AddCheckPoint( hybrid, 1, x1, vZ ); |
1316 | lowFLpiece.CreateAndInsertCableSegment(1); |
1317 | |
1318 | Double_t zChips[4] = {fgkHybFLlowChipZ1,fgkHybFLlowChipZ2, |
1319 | fgkHybFLlowChipZ3,fgkHybFLlowChipZ4}; |
1320 | Double_t vX[3] = {1,0,0}; |
1321 | for (Int_t i=0; i<3; i++) { |
1322 | char ch[20]; |
1323 | sprintf(ch, "lowFLpieceA%i", i+4); |
1324 | lowFLpiece.SetName(ch); |
1325 | lowFLpiece.SetWidth(zChips[i+1]-zChips[i]-fgkHybFLlowHoleDZ); |
1326 | |
1327 | lowFLpiece.SetLayer(1, fgkHybAlThick, alSDD, fColorAl); |
1328 | x1[0] = -fgkHybridWidth/2 + piece1width; |
1329 | x2[0] = x1[0] + fgkHybFLlowHolePasDX; |
1330 | Double_t zPiece = (zChips[i+1]+zChips[i])/2 - fgkHybridLength/2; |
1331 | x1[2] = zPiece; x2[2] = zPiece; |
1332 | lowFLpiece.AddCheckPoint( hybrid, 0, x2, vX ); |
1333 | lowFLpiece.AddCheckPoint( hybrid, 1, x1, vX ); |
1334 | lowFLpiece.CreateAndInsertCableSegment(1,90); |
1335 | lowFLpiece.ResetPoints(); |
1336 | |
1337 | sprintf(ch, "lowFLpieceB%i", i+4); |
1338 | lowFLpiece.SetName(ch); |
1339 | x1[0] = fgkHybridWidth/2 - piece3width; |
1340 | x2[0] = x1[0] - fgkHybFLlowHoleAmbDX; |
1341 | lowFLpiece.AddCheckPoint( hybrid, 0, x1, vX ); |
1342 | lowFLpiece.AddCheckPoint( hybrid, 1, x2, vX ); |
1343 | lowFLpiece.CreateAndInsertCableSegment(1,90); |
1344 | }; |
1345 | |
1346 | //**************************************************** chips+CC: |
1347 | AliITSv11GeomCableFlat chip("", fgkHybChipsDZ, chipsCCTotThick); |
1348 | chip.SetInitialNode(hybrid); |
1349 | chip.SetNLayers(5); |
1350 | chip.SetLayer(0, fgkHybUnderNiThick, niSDD, 2); |
1351 | chip.SetLayer(1, fgkHybGlueAgThick, glueAG, 4); |
1352 | chip.SetLayer(2, fgkHybChipThick, siliconSDD, fColorSilicon); |
1353 | chip.SetLayer(3, fgkHybUpCCThick+fgkHybUpCCThick, polyhamideSDD, |
1354 | fColorPolyhamide); |
1355 | chip.SetLayer(4, fgkHybAlCCThick+fgkHybAlCCThick, alSDD80p100, fColorAl); |
1356 | // Here the tho CC (low+up) are merged |
1357 | // In fact, the last layer has a smaller surface of Al -> I put 80% |
1358 | |
1359 | x1[1] = lowLayerYmin + chipsCCTotThick/2; |
1360 | x2[1] = x1[1]; |
1361 | char ch[20]; |
1362 | |
1363 | for (Int_t i=0; i<4; i++) { |
1364 | sprintf(ch, "pascalCC%i", i); |
1365 | chip.SetName(ch); |
1366 | x1[0] = fgkHybFLlowPasX - fgkHybridWidth/2 - fgkHybPascalDX/2; |
1367 | x2[0] = x1[0] + fgkHybPascalDX; |
1368 | x1[2] = zChips[i] - fgkHybridLength/2; |
1369 | x2[2] = x1[2]; |
1370 | chip.AddCheckPoint( hybrid, 0, x1, vX ); |
1371 | chip.AddCheckPoint( hybrid, 1, x2, vX ); |
1372 | chip.CreateAndInsertCableSegment(1,-90); |
1373 | chip.ResetPoints(); |
1374 | |
1375 | sprintf(ch, "ambraCC%i", i); |
1376 | chip.SetName(ch); |
1377 | x1[0] = fgkHybFLlowAmbX - fgkHybridWidth/2 - fgkHybAmbraDX/2; |
1378 | x2[0] = x1[0] + fgkHybAmbraDX; |
1379 | chip.AddCheckPoint( hybrid, 0, x1, vX ); |
1380 | chip.AddCheckPoint( hybrid, 1, x2, vX ); |
1381 | chip.CreateAndInsertCableSegment(1,-90); |
1382 | chip.ResetPoints(); |
1383 | }; |
c789ee28 |
1384 | |
b7943f00 |
1385 | //**************************************************** CC outside chips: |
1386 | // je crois qu'il n'y a pas de 2ieme couche d'alu ici ... |
1387 | for (Int_t i = 0; i<4; i++) { |
1388 | char ch[20]; |
1389 | sprintf(ch, "ccLayerA%i", i); |
1390 | |
1391 | AliITSv11GeomCableFlat ccLayer1(ch, 6.6*fgkmm, ccUpLayerTotThick); |
1392 | ccLayer1.SetInitialNode(hybrid); |
1393 | ccLayer1.SetNLayers(2); |
1394 | ccLayer1.SetLayer(0, 2*fgkHybUpCCThick, polyhamideSDD, fColorPolyhamide); |
1395 | ccLayer1.SetLayer(1, fgkHybAlCCThick, alSDD50p100, fColorAl); |
1396 | // Al at ~50% |
1397 | |
1398 | x1[0] = -fgkHybridWidth/2; |
1399 | x2[0] = fgkHybFLlowPasX - fgkHybridWidth/2 - fgkHybPascalDX/2; |
1400 | x1[1] = lowLayerYmin + fgkHybUnderNiThick + fgkHybGlueAgThick |
1401 | + fgkHybChipThick + ccUpLayerTotThick/2; |
1402 | x2[1] = x1[1]; |
1403 | x1[2] = zChips[i] - fgkHybridLength/2; |
1404 | x2[2] = x1[2]; |
1405 | ccLayer1.AddCheckPoint( hybrid, 0, x1, vX ); |
1406 | ccLayer1.AddCheckPoint( hybrid, 1, x2, vX ); |
1407 | ccLayer1.CreateAndInsertCableSegment(1,-90); |
1408 | |
1409 | sprintf(ch, "ccLayerB%i", i); |
1410 | AliITSv11GeomCableFlat ccLayer2(ch, fgkHybChipsDZ, ccUpLayerTotThick); |
1411 | ccLayer2.SetInitialNode(hybrid); |
1412 | ccLayer2.SetNLayers(2); |
1413 | ccLayer2.SetLayer(0, 2*fgkHybUpCCThick, polyhamideSDD, fColorPolyhamide); |
1414 | ccLayer2.SetLayer(1, fgkHybAlCCThick, alSDD50p100, fColorAl); |
1415 | // Al at ~50% |
1416 | |
1417 | x1[0] = -fgkHybridWidth/2 + fgkHybFLlowPasX + fgkHybPascalDX/2; |
1418 | x2[0] = -fgkHybridWidth/2 + fgkHybFLlowAmbX - fgkHybAmbraDX/2; |
1419 | ccLayer2.AddCheckPoint( hybrid, 0, x1, vX ); |
1420 | ccLayer2.AddCheckPoint( hybrid, 1, x2, vX ); |
1421 | ccLayer2.CreateAndInsertCableSegment(1,-90); |
1422 | ccLayer2.ResetPoints(); |
1423 | sprintf(ch, "ccLayerC%i", i); |
1424 | ccLayer2.SetName(ch); |
1425 | x1[0] = -fgkHybridWidth/2 + fgkHybFLlowAmbX + fgkHybAmbraDX/2; |
1426 | x2[0] = fgkHybridWidth/2 - fgkHybFLUpperWidth + 3*fgkmm; |
1427 | x1[1] = lowLayerYmin + lowFLTotalThick + flUpThick + fgkHybAlThick |
1428 | + ccUpLayerTotThick/2; |
1429 | x2[1] = x1[1]; |
1430 | |
1431 | ccLayer2.AddCheckPoint( hybrid, 0, x1, vX ); |
1432 | ccLayer2.AddCheckPoint( hybrid, 1, x2, vX ); |
1433 | ccLayer2.CreateAndInsertCableSegment(1,-90); |
1434 | }; |
c789ee28 |
1435 | |
b7943f00 |
1436 | //**************************************************** FL UP: |
1437 | // (last Al layer will be a special triangular shape) |
1438 | TGeoBBox *sFLupPolyhamide = new TGeoBBox("sFLupPolyhamide", |
1439 | fgkHybFLUpperWidth/2, flUpThick/2, |
1440 | fgkHybFLUpperLength/2); |
1441 | TGeoVolume *vFLupPolyhamide = new TGeoVolume("vFLupPolyhamide", |
1442 | sFLupPolyhamide, polyhamideSDD); |
1443 | vFLupPolyhamide->SetLineColor(fColorPolyhamide); |
1444 | TGeoTranslation *trFLupPolyhamide = |
1445 | new TGeoTranslation(fgkHybridWidth/2-fgkHybFLUpperWidth/2, |
1446 | lowLayerYmin+lowFLTotalThick+flUpThick/2,0); |
1447 | |
1448 | hybrid->AddNode(vFLupPolyhamide, 1, trFLupPolyhamide); |
1449 | |
1450 | TGeoArb8 *aluStrip = new TGeoArb8(fgkHybAlThick/2); |
1451 | aluStrip->SetVertex( 0,-fgkHybFLUpperAlDZ/2, fgkHybFLUpperWidth); |
1452 | aluStrip->SetVertex( 1, fgkHybFLUpperAlDZ/2, fgkHybFLUpperWidth); |
1453 | aluStrip->SetVertex( 2, fgkHybFLUpperAlDZ/2, fgkHybFLUpperWidth-fgkHybFLUpperAldx); |
1454 | aluStrip->SetVertex( 3,-fgkHybFLUpperAlDZ/2, 0); |
1455 | aluStrip->SetVertex( 4,-fgkHybFLUpperAlDZ/2, fgkHybFLUpperWidth); |
1456 | aluStrip->SetVertex( 5, fgkHybFLUpperAlDZ/2, fgkHybFLUpperWidth); |
1457 | aluStrip->SetVertex( 6, fgkHybFLUpperAlDZ/2, fgkHybFLUpperWidth-fgkHybFLUpperAldx); |
1458 | aluStrip->SetVertex( 7,-fgkHybFLUpperAlDZ/2, 0); |
1459 | TGeoVolume *vAluStrip = new TGeoVolume("vAluStrip",aluStrip, alSDD50p100); |
1460 | // Al at ~50% |
1461 | |
1462 | vAluStrip->SetLineColor(fColorAl); |
531d6cdc |
1463 | //TGeoRotation rotAluStrip("rotAluStrip",0, -90, 90); |
1464 | TGeoRotation *rotAluStrip = new TGeoRotation("rotAluStrip",0, -90, 90); |
1465 | |
b7943f00 |
1466 | Double_t yRotAluStrip = lowLayerYmin+lowFLTotalThick |
1467 | +flUpThick+fgkHybAlThick/2; |
1468 | TGeoCombiTrans *aluStripTr1 = new TGeoCombiTrans( |
1469 | fgkHybridWidth/2,yRotAluStrip, |
531d6cdc |
1470 | fgkHybridLength/2-fgkHybFLlowChipZ1+1*fgkmm, rotAluStrip); |
b7943f00 |
1471 | TGeoCombiTrans *aluStripTr2 = new TGeoCombiTrans(*aluStripTr1); |
1472 | AddTranslationToCombiTrans(aluStripTr2,0,0, |
1473 | fgkHybFLlowChipZ1-fgkHybFLlowChipZ2); |
1474 | TGeoCombiTrans *aluStripTr3 = new TGeoCombiTrans(*aluStripTr2); |
1475 | AddTranslationToCombiTrans(aluStripTr3,0,0, |
1476 | fgkHybFLlowChipZ2-fgkHybFLlowChipZ3); |
1477 | TGeoCombiTrans *aluStripTr4 = new TGeoCombiTrans(*aluStripTr3); |
1478 | AddTranslationToCombiTrans(aluStripTr4,0,0, |
1479 | fgkHybFLlowChipZ3-fgkHybFLlowChipZ4); |
1480 | |
1481 | hybrid->AddNode(vAluStrip, 1, aluStripTr1); |
1482 | hybrid->AddNode(vAluStrip, 2, aluStripTr2); |
1483 | hybrid->AddNode(vAluStrip, 3, aluStripTr3); |
1484 | hybrid->AddNode(vAluStrip, 4, aluStripTr4); |
b7943f00 |
1485 | //**************************************************** SMD: |
1486 | TGeoBBox *hybSMD = new TGeoBBox("ITSsddSMDshape", |
1487 | fgkHybSMDmiddleL/2+fgkHybSMDendL, |
1488 | fgkHybSMDheight/2,fgkHybSMDendW/2); |
1489 | TGeoVolume *vHybSMD = new TGeoVolume("ITSsddSMD",hybSMD,airSDD); |
1490 | |
1491 | TGeoBBox *hybSMDmiddle = new TGeoBBox("ITSsddSMDmiddleShape", |
1492 | fgkHybSMDmiddleL/2,fgkHybSMDheight/2, |
1493 | fgkHybSMDmiddleW/2); |
1494 | TGeoVolume *vHybSMDmiddle = new TGeoVolume("ITSsddSMDmiddle", |
1495 | hybSMDmiddle,medSMD); |
1496 | vHybSMDmiddle->SetLineColor(fColorSMD); |
1497 | TGeoBBox *hybSMDend = new TGeoBBox("ITSsddSMDendShape", |
1498 | fgkHybSMDendL/2,fgkHybSMDheight/2,fgkHybSMDendW/2); |
1499 | TGeoVolume *vHybSMDend = new TGeoVolume("ITSsddSMDend", |
1500 | hybSMDend,medSMDweld); |
1501 | vHybSMDend->SetLineColor(fColorSMDweld); |
1502 | TGeoTranslation *vHybSMDendTr1 = new TGeoTranslation("", |
1503 | (fgkHybSMDmiddleL+fgkHybSMDendL)/2,0,0); |
1504 | TGeoTranslation *vHybSMDendTr2 = new TGeoTranslation("", |
1505 | -(fgkHybSMDmiddleL+fgkHybSMDendL)/2,0,0); |
1506 | vHybSMD->AddNode(vHybSMDmiddle,1,0); |
1507 | vHybSMD->AddNode(vHybSMDend,1,vHybSMDendTr1); |
1508 | vHybSMD->AddNode(vHybSMDend,2,vHybSMDendTr2); |
1509 | for (Int_t i=0; i<fgkNHybSMD; i++) { |
1510 | TGeoTranslation *vHybSMDtr = new TGeoTranslation("", |
1511 | -fgkHybridWidth/2+fgkHybSMDposX[i], |
1512 | lowLayerYmin+lowFLTotalThick+fgkHybSMDheight/2, |
1513 | -fgkHybridLength/2+fgkHybSMDposZ[i]); |
1514 | hybrid->AddNode(vHybSMD, i+1, vHybSMDtr); |
1515 | }; |
1516 | |
1517 | |
1518 | if (iLRSide == 0) { |
1519 | }; |
1520 | |
1521 | if(GetDebug(1)) hybrid->CheckOverlaps(0.01); |
1522 | hybrid->SetVisibility(kFALSE); |
1523 | return hybrid; |
531d6cdc |
1524 | } |
c789ee28 |
1525 | |
1526 | |
db486a6e |
1527 | //________________________________________________________________________ |
c789ee28 |
1528 | TGeoVolume* AliITSv11GeometrySDD::CreateLadderSegment(Int_t iLay, Int_t iSeg) { |
b7943f00 |
1529 | // |
1530 | // Return a box volume containing a segment of a ladder. |
1531 | // |
1532 | |
1533 | TGeoMedium *airSDD = GetMedium("ITSair"); |
1534 | TGeoMedium *phynoxSDD = GetMedium("ITSal"); // phynoxSDD To code !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
1535 | TGeoMedium *coolerMediumSDD = GetMedium("WATER"); |
1536 | |
1537 | Double_t tDY = fgkLadderSegBoxDH/2; //space left on top of the ladder |
1538 | Double_t segmentLength = fgkSegmentLength; |
1539 | Double_t spaceBetweenCables = 500*fgkmicron; |
1540 | |
1541 | //***************************************** |
1542 | // Set parameters according to (iLay,iSeg): |
1543 | //***************************************** |
1544 | Int_t nDetectors = fgkLay3Ndet; |
1545 | Double_t coolPipeSuppH = fgkLay3CoolPipeSuppH; |
1546 | Double_t sensorCenterZPos = fLay3sensorZPos[iSeg]- |
1547 | (fgkSegmentLength*fgkLay3Ndet/2. - |
1548 | fgkSegmentLength/2-(iSeg)*fgkSegmentLength); |
1549 | // sensorCenterZPos = z in segment local coord syst. |
1550 | |
1551 | AliITSv11GeomCableFlat *digitCableA = fDigitCableLay3A; |
1552 | AliITSv11GeomCableFlat *digitCableB = fDigitCableLay3B; |
1553 | |
1554 | if (iLay==3) { |
1555 | } else if (iLay==4) { |
1556 | nDetectors = fgkLay4Ndet; |
1557 | coolPipeSuppH = fgkLay4CoolPipeSuppH; |
1558 | sensorCenterZPos = fLay4sensorZPos[iSeg]- |
1559 | (fgkSegmentLength*fgkLay4Ndet/2. - |
1560 | fgkSegmentLength/2-(iSeg)*fgkSegmentLength); |
1561 | digitCableA = fDigitCableLay4A; |
1562 | digitCableB = fDigitCableLay4B; |
1563 | } else |
1564 | printf("AliITSv11GeometrySDD::CreateLadderSegment Wrong layer index !"); |
1565 | |
1566 | |
1567 | Double_t cableSideSign = -1; |
1568 | if (iSeg<nDetectors/2) cableSideSign = 1; |
1569 | Double_t spaceForCables = spaceBetweenCables* |
1570 | (nDetectors-TMath::Abs(nDetectors-2*iSeg-1)-1)/2 |
1571 | +0.1*fgkmicron; |
1572 | // gives [0-1-2-2-1-0]*spaceBetweenCables |
1573 | // or [0-1-2-3-3-2-1-0]*spaceBetweenCables |
1574 | Int_t iUpdateCableMin; |
1575 | Int_t iUpdateCableMax; |
1576 | if (cableSideSign==-1) { |
1577 | iUpdateCableMin = nDetectors/2; |
1578 | iUpdateCableMax = iSeg-1; |
1579 | } else { |
1580 | iUpdateCableMin = iSeg+1; |
1581 | iUpdateCableMax = nDetectors/2-1; |
1582 | }; |
1583 | |
1584 | if(GetDebug(1)){ |
1585 | cout << "Segment ("<< iLay <<',' << iSeg |
1586 | << ") : sensor z shift in local segment coord.=" |
1587 | << sensorCenterZPos << endl; |
1588 | }; |
1589 | |
1590 | //**************************** |
1591 | // The segment volume |
1592 | //**************************** |
1593 | TGeoBBox *segBox = new TGeoBBox("ITSsddSegBox", |
1594 | fgkLadderWidth/2+fgkPinSuppWidth+fgkLadderSegBoxDW, |
1595 | fgkLadderHeight/2+fgkLadderSegBoxDH/2, |
1596 | segmentLength/2); |
1597 | |
1598 | TGeoVolume *virtualSeg = new TGeoVolume("ITSsddSegment", |
1599 | segBox, airSDD); |
1600 | |
1601 | //****************************** |
1602 | // Carbon fiber structure : |
1603 | //****************************** |
1604 | |
1605 | virtualSeg->AddNode(fLaddSegCommonVol[0], 1, fLaddSegCommonTr[0]); |
1606 | Int_t volumeIndex = 1; |
1607 | for (Int_t i = 1; i<fgkNladdSegCommonVol;i++ ) { |
1608 | if (fLaddSegCommonVol[i]==fLaddSegCommonVol[i-1]) |
1609 | volumeIndex++; |
1610 | else |
1611 | volumeIndex = 1; |
1612 | virtualSeg->AddNode(fLaddSegCommonVol[i], volumeIndex, |
1613 | fLaddSegCommonTr[i]); |
1614 | }; |
1615 | |
1616 | //********************************** |
1617 | // Pine support of the sensors : |
1618 | //********************************** |
531d6cdc |
1619 | TGeoRotation *rotPS1 = new TGeoRotation("",0,-90,90); |
1620 | TGeoRotation *rotPS2 = new TGeoRotation("",0,-90,-90); |
1621 | |
b7943f00 |
1622 | TGeoCombiTrans *transPS1 = new TGeoCombiTrans( fgkPinDYOnSensor, |
1623 | - fgkLadderHeight/2.-tDY |
1624 | + fgkPinSuppHeight/2., |
531d6cdc |
1625 | sensorCenterZPos+fgkPinDXminOnSensor,rotPS1); |
b7943f00 |
1626 | TGeoCombiTrans *transPS2 = new TGeoCombiTrans(*transPS1); |
1627 | AddTranslationToCombiTrans(transPS2, 0, 0, fgkPinPinDDXOnSensor); |
1628 | TGeoCombiTrans *transPS3 = new TGeoCombiTrans(*transPS1); |
1629 | AddTranslationToCombiTrans(transPS3, 0, 0, -2*fgkPinDXminOnSensor); |
1630 | TGeoCombiTrans *transPS4 = new TGeoCombiTrans(*transPS3); |
1631 | AddTranslationToCombiTrans(transPS4, 0, 0, -fgkPinPinDDXOnSensor); |
1632 | |
1633 | TGeoCombiTrans *transPS5 = new TGeoCombiTrans( -fgkPinDYOnSensor, |
1634 | - fgkLadderHeight/2. - tDY |
1635 | + fgkPinSuppHeight/2., |
531d6cdc |
1636 | sensorCenterZPos+fgkPinDXminOnSensor,rotPS2); |
b7943f00 |
1637 | TGeoCombiTrans *transPS6 = new TGeoCombiTrans(*transPS5); |
1638 | AddTranslationToCombiTrans(transPS6, 0, 0, fgkPinPinDDXOnSensor); |
1639 | TGeoCombiTrans *transPS7 = new TGeoCombiTrans(*transPS5); |
1640 | AddTranslationToCombiTrans(transPS7, 0, 0, -2*fgkPinDXminOnSensor); |
1641 | TGeoCombiTrans *transPS8 = new TGeoCombiTrans(*transPS7); |
1642 | AddTranslationToCombiTrans(transPS8, 0, 0, -fgkPinPinDDXOnSensor); |
1643 | |
1644 | virtualSeg->AddNode(fPinSupport, 1, transPS1); |
1645 | virtualSeg->AddNode(fPinSupport, 2, transPS2); |
1646 | virtualSeg->AddNode(fPinSupport, 3, transPS3); |
1647 | virtualSeg->AddNode(fPinSupport, 4, transPS4); |
1648 | virtualSeg->AddNode(fPinSupport, 5, transPS5); |
1649 | virtualSeg->AddNode(fPinSupport, 6, transPS6); |
1650 | virtualSeg->AddNode(fPinSupport, 7, transPS7); |
1651 | virtualSeg->AddNode(fPinSupport, 8, transPS8); |
1652 | |
1653 | //****************************** |
1654 | // Cooling pipe supports : |
1655 | //****************************** |
1656 | Double_t triangleHeight = fgkLadderHeight - fgkLadderBeamRadius; |
1657 | Double_t halfTheta = TMath::ATan( 0.5*fgkLadderWidth/triangleHeight ); |
1658 | Double_t triangleCPaxeDist = fgkCoolPipeSuppAxeDist-fgkCoolPipeSuppWidthExt- |
1659 | fgkCoolPipeSuppWidthIn+fgkLadderBeamRadius; |
1660 | |
1661 | Double_t coolPipeSuppL = TMath::Tan(halfTheta)* |
1662 | (triangleHeight+triangleCPaxeDist/ |
1663 | TMath::Sin(halfTheta)-coolPipeSuppH); |
1664 | if (fAddCoolingSyst) { |
531d6cdc |
1665 | TGeoRotation *rotCPS2 = new TGeoRotation("", -halfTheta*TMath::RadToDeg(), -90, 90); |
1666 | TGeoRotation *rotCPS1 = new TGeoRotation("", halfTheta*TMath::RadToDeg(), -90, -90); |
b7943f00 |
1667 | TGeoCombiTrans *transCPS1 = new TGeoCombiTrans(coolPipeSuppL, |
fa4639a3 |
1668 | -fgkLadderHeight/2. - tDY |
b7943f00 |
1669 | +coolPipeSuppH+fgkLadderBeamRadius, |
531d6cdc |
1670 | -segmentLength/2., rotCPS1); |
b7943f00 |
1671 | TGeoCombiTrans *transCPS3 = new TGeoCombiTrans(*transCPS1); |
1672 | AddTranslationToCombiTrans(transCPS3, 0, 0, segmentLength); |
1673 | |
1674 | TGeoCombiTrans *transCPS2 = new TGeoCombiTrans(-coolPipeSuppL, |
fa4639a3 |
1675 | -fgkLadderHeight/2.- tDY |
b7943f00 |
1676 | +coolPipeSuppH+fgkLadderBeamRadius, |
531d6cdc |
1677 | segmentLength/2., rotCPS2); |
b7943f00 |
1678 | TGeoCombiTrans *transCPS4 = new TGeoCombiTrans(*transCPS2); |
1679 | AddTranslationToCombiTrans(transCPS4, 0, 0, -segmentLength); |
1680 | |
1681 | virtualSeg->AddNode(fCoolPipeSupportL, 1, transCPS1); |
1682 | virtualSeg->AddNode(fCoolPipeSupportL, 2, transCPS2); |
1683 | virtualSeg->AddNode(fCoolPipeSupportR, 1, transCPS3); |
1684 | virtualSeg->AddNode(fCoolPipeSupportR, 2, transCPS4); |
1685 | }; |
1686 | |
1687 | //************************ |
1688 | // Cooling pipes : |
1689 | //************************ |
1690 | TGeoTranslation *pipeTr1 = new TGeoTranslation(coolPipeSuppL, |
1691 | -fgkLadderHeight/2. - tDY + |
1692 | fgkLadderBeamRadius+coolPipeSuppH, 0); |
1693 | TGeoTranslation *pipeTr2 = new TGeoTranslation(-coolPipeSuppL, |
1694 | -fgkLadderHeight/2.- tDY + |
1695 | fgkLadderBeamRadius+coolPipeSuppH, 0); |
1696 | |
1697 | if (fAddCoolingSyst) { |
1698 | TGeoTube *coolingPipeShape = new TGeoTube( fgkCoolPipeInnerDiam/2, |
1699 | fgkCoolPipeOuterDiam/2, |
1700 | segmentLength/2); |
1701 | TGeoTube *coolerShape = new TGeoTube( 0, fgkCoolPipeInnerDiam/2, |
1702 | segmentLength/2); |
1703 | |
1704 | TGeoVolume *coolingPipe = new TGeoVolume("ITSsddCoolingPipe", |
1705 | coolingPipeShape, phynoxSDD ); |
1706 | coolingPipe->SetLineColor(fColorPhynox); |
1707 | TGeoVolume *cooler = new TGeoVolume("ITSsddCoolingLiquid",coolerShape, |
1708 | coolerMediumSDD ); |
1709 | |
1710 | |
1711 | virtualSeg->AddNode(coolingPipe, 1, pipeTr1); |
1712 | virtualSeg->AddNode(coolingPipe, 2, pipeTr2); |
1713 | if (fCoolingOn) { |
1714 | virtualSeg->AddNode(cooler, 1, pipeTr1); |
1715 | virtualSeg->AddNode(cooler, 2, pipeTr2); |
1716 | }; |
1717 | }; |
c789ee28 |
1718 | |
b7943f00 |
1719 | //********************************** |
1720 | // Bases of hybrid thermal bridges |
1721 | //********************************** |
1722 | Double_t shiftHyb = 1.05; // shift between thermal Bridge base and thermal bridge |
1723 | // approx !!! not clear on 0752/14-A |
1724 | if (fAddCoolingSyst) { |
531d6cdc |
1725 | TGeoRotation *rotHybrid1 = new TGeoRotation("", 0, 0, -90 - fgkHybridAngle); |
1726 | TGeoRotation *rotHybrid2 = new TGeoRotation("", 0 ,180, 90 - fgkHybridAngle); |
1727 | TGeoCombiTrans *baseTr1 = new TGeoCombiTrans(*pipeTr2, *rotHybrid1); |
1728 | TGeoCombiTrans *baseTr2 = new TGeoCombiTrans(*pipeTr1, *rotHybrid2); |
b7943f00 |
1729 | |
1730 | virtualSeg->AddNode(fBaseThermalBridge, 1, baseTr1); |
1731 | virtualSeg->AddNode(fBaseThermalBridge, 2, baseTr2); |
1732 | }; |
1733 | |
1734 | //************************* |
1735 | // the 2 hybrids : |
1736 | //************************* |
1737 | Double_t hybDy = ((TGeoBBox*)fHybrid->GetShape())->GetDY(); |
1738 | Double_t distAxeToHybridCenter = fgkBTBaxisAtoBase+hybDy; |
1739 | |
1740 | Double_t hybrVolX = ( distAxeToHybridCenter*CosD(fgkHybridAngle) |
1741 | - shiftHyb*SinD(fgkHybridAngle) ); |
1742 | Double_t hybrVolY = ( distAxeToHybridCenter*SinD(fgkHybridAngle) |
1743 | + shiftHyb*CosD(fgkHybridAngle) ); |
1744 | if (fAddHybrids) { |
1745 | TGeoRotation rotHybrid3("", 0, 0, 90. - fgkHybridAngle); |
1746 | TGeoRotation rotHybrid4("", 0 ,180, -90. - fgkHybridAngle); |
1747 | TGeoCombiTrans *hybTr1 = new TGeoCombiTrans(*pipeTr2, rotHybrid3); |
1748 | TGeoCombiTrans *hybTr2 = new TGeoCombiTrans(*pipeTr1, rotHybrid4); |
1749 | AddTranslationToCombiTrans( hybTr1, -hybrVolX, hybrVolY, 0); |
1750 | AddTranslationToCombiTrans( hybTr2, hybrVolX, hybrVolY, 0); |
1751 | |
1752 | virtualSeg->AddNode(fHybrid, 1, hybTr1); |
1753 | virtualSeg->AddNode(fHybrid, 2, hybTr2); |
1754 | }; |
1755 | |
1756 | //*********** |
1757 | // cables |
1758 | //*********** |
1759 | if (fAddCables) { |
1760 | // Starting from this segment |
1761 | Double_t hybDz = ((TGeoBBox*)fHybrid->GetShape())->GetDZ(); |
1762 | Double_t hybDx = ((TGeoBBox*)fHybrid->GetShape())->GetDX(); |
1763 | Double_t posDigitCableAlongHyb = shiftHyb+ hybDx |
1764 | - digitCableA->GetWidth()/2; |
1765 | Double_t distAxeToDigitCableCenter = distAxeToHybridCenter+hybDy |
1766 | - digitCableA->GetThickness()/2; |
1767 | |
1768 | Double_t digitCableX = ( coolPipeSuppL |
1769 | + distAxeToDigitCableCenter*CosD(fgkHybridAngle) |
1770 | - posDigitCableAlongHyb*SinD(fgkHybridAngle) ); |
1771 | Double_t digitCableY = ( - fgkLadderHeight/2.-TMath::Abs(tDY) |
1772 | + fgkLadderBeamRadius+coolPipeSuppH |
1773 | + distAxeToDigitCableCenter*SinD(fgkHybridAngle) |
1774 | + posDigitCableAlongHyb*CosD(fgkHybridAngle) ); |
1775 | |
1776 | |
1777 | Double_t digitCableCenterA0[3]={ -cableSideSign*digitCableX, |
1778 | digitCableY, cableSideSign*hybDz }; |
1779 | Double_t digitCableCenterA1[3] = { |
1780 | -cableSideSign*(digitCableX+spaceForCables*CosD(fgkHybridAngle)), |
1781 | digitCableY+spaceForCables*SinD(fgkHybridAngle), |
1782 | cableSideSign*segmentLength/2 }; |
1783 | |
1784 | Double_t digitCableCenterB0[3]={ cableSideSign*digitCableX, |
1785 | digitCableY,cableSideSign*hybDz}; |
1786 | Double_t digitCableCenterB1[3]={ |
1787 | cableSideSign*(digitCableX+spaceForCables*CosD(fgkHybridAngle)), |
1788 | digitCableY+spaceForCables*SinD(fgkHybridAngle), |
1789 | cableSideSign*segmentLength/2 }; |
1790 | |
1791 | Double_t vZ[3] = {0,0,1}; |
1792 | digitCableA[iSeg].AddCheckPoint( virtualSeg, 0, digitCableCenterA0, vZ); |
1793 | digitCableA[iSeg].AddCheckPoint( virtualSeg, 1, digitCableCenterA1, vZ); |
1794 | digitCableB[iSeg].AddCheckPoint( virtualSeg, 0, digitCableCenterB0, vZ); |
1795 | digitCableB[iSeg].AddCheckPoint( virtualSeg, 1, digitCableCenterB1, vZ); |
1796 | |
1797 | // Updating the other cables |
1798 | for (Int_t iCable=iUpdateCableMin; iCable<=iUpdateCableMax; iCable++) { |
1799 | |
1800 | Int_t iPoint = TMath::Abs(iCable-iSeg)+1; |
1801 | Double_t coord[3]; |
1802 | digitCableA[iCable].GetPoint( 1, coord); |
1803 | digitCableA[iCable].AddCheckPoint( virtualSeg, iPoint, coord, vZ); |
1804 | digitCableB[iCable].GetPoint( 1, coord); |
1805 | digitCableB[iCable].AddCheckPoint( virtualSeg, iPoint, coord, vZ); |
1806 | }; |
1807 | }; |
1808 | |
1809 | //********************************** |
1810 | if(GetDebug(1)) virtualSeg->CheckOverlaps(0.01); |
fa4639a3 |
1811 | virtualSeg->SetVisibility(kFALSE); |
b7943f00 |
1812 | return virtualSeg; |
531d6cdc |
1813 | } |
db486a6e |
1814 | |
c789ee28 |
1815 | |
1816 | //________________________________________________________________________ |
1817 | TGeoVolume* AliITSv11GeometrySDD::CreatePinSupport() { |
1818 | // |
1819 | // Create a pine support |
1820 | // axis of rotation is the cone axis, center in its middle |
1821 | // |
b7943f00 |
1822 | TGeoCone *cone = new TGeoCone("ITSsddPinSuppCone",fgkPinSuppHeight/2., |
1823 | 0,fgkPinSuppRmax,0,fgkPinSuppRmax- |
1824 | fgkPinSuppHeight*TanD(fgkPinSuppConeAngle) ); |
1825 | TGeoBBox *tong = new TGeoBBox("ITSsddPinSuppTong",fgkPinSuppRmax, |
1826 | fgkPinSuppLength/2.,fgkPinSuppThickness/2.); |
1827 | TGeoTube *hole = new TGeoTube("ITSsddPinSuppHole",0,fgkPinR, |
1828 | fgkPinSuppHeight/2.); |
c789ee28 |
1829 | if(GetDebug(3)){// Remove compiler warning. |
db486a6e |
1830 | cone->InspectShape(); |
1831 | tong->InspectShape(); |
1832 | hole->InspectShape(); |
c789ee28 |
1833 | }; |
db486a6e |
1834 | |
1835 | TGeoTranslation *tongTrans = new TGeoTranslation("ITSsddPinSuppTongTr",0, |
b7943f00 |
1836 | fgkPinSuppLength/2.,-fgkPinSuppHeight/2.+fgkPinSuppThickness/2.); |
db486a6e |
1837 | tongTrans->RegisterYourself(); |
1838 | TGeoCompositeShape *pinSupportShape = new TGeoCompositeShape( |
b7943f00 |
1839 | "ITSsddPinSupportShape","(ITSsddPinSuppCone+" |
db486a6e |
1840 | "ITSsddPinSuppTong:ITSsddPinSuppTongTr)-ITSsddPinSuppHole"); |
1841 | |
b7943f00 |
1842 | |
1843 | TGeoMedium *rytonSDD = GetMedium("ITSsddCarbonM55J"); //medium = ryton ? To code !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
1844 | TGeoVolume *pinSupport = new TGeoVolume("ITSsddPinSupport",pinSupportShape, |
db486a6e |
1845 | rytonSDD); |
c789ee28 |
1846 | pinSupport->SetLineColor(fColorRyton); |
db486a6e |
1847 | return pinSupport; |
c789ee28 |
1848 | // include the pin itself !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
531d6cdc |
1849 | } |
c789ee28 |
1850 | |
b7943f00 |
1851 | |
db486a6e |
1852 | //________________________________________________________________________ |
c789ee28 |
1853 | TGeoVolume* AliITSv11GeometrySDD::CreateCoolPipeSupportL() { |
1854 | // |
1855 | // Create half of the cooling pipe support (ALR-0752/3) |
1856 | // |
1857 | |
b7943f00 |
1858 | Double_t diffX = fgkCoolPipeSuppHeight*TanD(fgkCoolPipeSuppAngle); |
c789ee28 |
1859 | |
b7943f00 |
1860 | TGeoArb8 *side1 = new TGeoArb8(fgkCoolPipeSuppHeight/2.); |
1861 | side1->SetVertex( 0, 0, -fgkCoolPipeSuppWidthExt/2.); |
1862 | side1->SetVertex( 1, fgkCoolPipeSuppMaxLength/2.-diffX, |
1863 | -fgkCoolPipeSuppWidthExt/2.); |
1864 | side1->SetVertex( 2, fgkCoolPipeSuppMaxLength/2.-diffX, |
1865 | fgkCoolPipeSuppWidthExt/2.); |
1866 | side1->SetVertex( 3, 0, fgkCoolPipeSuppWidthExt/2.); |
1867 | side1->SetVertex( 4, 0, -fgkCoolPipeSuppWidthExt/2.); |
1868 | side1->SetVertex( 5, fgkCoolPipeSuppMaxLength/2., |
1869 | -fgkCoolPipeSuppWidthExt/2.); |
1870 | side1->SetVertex( 6, fgkCoolPipeSuppMaxLength/2., |
1871 | fgkCoolPipeSuppWidthExt/2.); |
1872 | side1->SetVertex( 7, 0, fgkCoolPipeSuppWidthExt/2.); |
c789ee28 |
1873 | side1->SetName("ITSsddCPSside1"); |
1874 | |
1875 | TGeoTranslation *side1Tr = new TGeoTranslation("ITSsddCPStr1",0, |
b7943f00 |
1876 | - fgkCoolPipeSuppAxeDist |
1877 | + fgkCoolPipeSuppWidthExt/2., 0); |
c789ee28 |
1878 | side1Tr->RegisterYourself(); |
1879 | TGeoTranslation *side2Tr = new TGeoTranslation("ITSsddCPStr2",0, |
b7943f00 |
1880 | - fgkCoolPipeSuppAxeDist |
1881 | + fgkCoolPipeSuppWidthExt*3/2. |
1882 | + fgkCoolPipeSuppWidthIn,0); |
c789ee28 |
1883 | side2Tr->RegisterYourself(); |
1884 | |
1885 | TGeoBBox *middle = new TGeoBBox("ITSsddCPSmiddle", |
b7943f00 |
1886 | (fgkCoolPipeSuppMaxLength/2.-fgkCoolPipeSuppSlitL)/2., |
1887 | fgkCoolPipeSuppWidthIn/2., fgkCoolPipeSuppHeight/2.); |
c789ee28 |
1888 | TGeoTranslation *middleTr = |
1889 | new TGeoTranslation("ITSsddCPStr3", |
b7943f00 |
1890 | (fgkCoolPipeSuppMaxLength/2.-fgkCoolPipeSuppSlitL)/2., |
1891 | -fgkCoolPipeSuppAxeDist+fgkCoolPipeSuppWidthExt |
1892 | +fgkCoolPipeSuppWidthIn/2., 0); |
c789ee28 |
1893 | middleTr->RegisterYourself(); |
1894 | |
1895 | TGeoBBox *axeBox = new TGeoBBox("ITSsddCPSaxeBox", |
b7943f00 |
1896 | fgkCoolPipeSuppTongW/4., |
1897 | (fgkCoolPipeSuppFulWidth |
1898 | - 2*fgkCoolPipeSuppWidthExt |
1899 | - fgkCoolPipeSuppWidthIn)/2, |
1900 | fgkCoolPipeSuppHeight/2.); |
c789ee28 |
1901 | |
1902 | TGeoTranslation *axeBoxTr = new TGeoTranslation("ITSsddCPSAxBoxTr", |
b7943f00 |
1903 | fgkCoolPipeSuppTongW/4., |
1904 | - fgkCoolPipeSuppAxeDist |
1905 | + fgkCoolPipeSuppFulWidth |
c789ee28 |
1906 | - axeBox->GetDY(), 0); |
1907 | axeBoxTr->RegisterYourself(); |
1908 | |
b7943f00 |
1909 | TGeoTube *axe = new TGeoTube("ITSsddCPSaxe",0,fgkCoolPipeSuppHoleDiam/2., |
1910 | fgkCoolPipeSuppTongW/4.); |
c789ee28 |
1911 | |
531d6cdc |
1912 | TGeoRotation *axeRot = new TGeoRotation("ITSsddCPSaxeRot",90,90,0); |
c789ee28 |
1913 | TGeoCombiTrans *axeTrans = new TGeoCombiTrans("ITSsddCPSaxeTr", |
531d6cdc |
1914 | fgkCoolPipeSuppTongW/4.,0,0,axeRot); |
c789ee28 |
1915 | axeTrans->RegisterYourself(); |
531d6cdc |
1916 | //delete axeRot; // make the code crash, no idea of why !!! |
c789ee28 |
1917 | |
1918 | if(GetDebug(3)){ |
1919 | middle->InspectShape(); |
1920 | axe->InspectShape(); |
1921 | }; |
1922 | |
b7943f00 |
1923 | TGeoMedium *rytonSDD = GetMedium("ITSsddCarbonM55J"); //medium = ryton ? !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
c789ee28 |
1924 | |
1925 | TGeoCompositeShape *coolPipeSuppShape = new TGeoCompositeShape( |
1926 | "ITSsddCoolPipeSuppShapeL", |
1927 | "ITSsddCPSmiddle:ITSsddCPStr3" |
1928 | "+ITSsddCPSside1:ITSsddCPStr1" |
1929 | "+ITSsddCPSside1:ITSsddCPStr2" |
1930 | "+ITSsddCPSaxeBox:ITSsddCPSAxBoxTr" |
1931 | "-ITSsddCPSaxe:ITSsddCPSaxeTr"); |
1932 | TGeoVolume *coolPipeSupp = new TGeoVolume("ITSsddCoolPipeSupportL", |
1933 | coolPipeSuppShape, rytonSDD); |
1934 | |
1935 | coolPipeSupp->SetLineColor(fColorRyton); |
531d6cdc |
1936 | |
c789ee28 |
1937 | return coolPipeSupp; |
531d6cdc |
1938 | } |
1939 | |
c789ee28 |
1940 | |
1941 | //________________________________________________________________________ |
1942 | TGeoVolume* AliITSv11GeometrySDD::CreateCoolPipeSupportR() { |
1943 | // |
1944 | //Create half of the cooling pipe support (ALR-0752/3) |
1945 | // |
1946 | |
b7943f00 |
1947 | Double_t diffX = fgkCoolPipeSuppHeight*TanD(fgkCoolPipeSuppAngle); |
c789ee28 |
1948 | |
b7943f00 |
1949 | TGeoArb8 *side1 = new TGeoArb8(fgkCoolPipeSuppHeight/2.); |
1950 | side1->SetVertex( 0, 0, -fgkCoolPipeSuppWidthExt/2.); |
1951 | side1->SetVertex( 1, -(fgkCoolPipeSuppMaxLength/2.-diffX), |
1952 | -fgkCoolPipeSuppWidthExt/2.); |
1953 | side1->SetVertex( 2, -(fgkCoolPipeSuppMaxLength/2.-diffX), |
1954 | fgkCoolPipeSuppWidthExt/2.); |
1955 | side1->SetVertex( 3, 0, fgkCoolPipeSuppWidthExt/2.); |
1956 | side1->SetVertex( 4, 0, -fgkCoolPipeSuppWidthExt/2.); |
1957 | side1->SetVertex( 5, -fgkCoolPipeSuppMaxLength/2., |
1958 | -fgkCoolPipeSuppWidthExt/2.); |
1959 | side1->SetVertex( 6, -fgkCoolPipeSuppMaxLength/2., |
1960 | fgkCoolPipeSuppWidthExt/2.); |
1961 | side1->SetVertex( 7, 0, fgkCoolPipeSuppWidthExt/2.); |
c789ee28 |
1962 | side1->SetName("ITSsddCPSside1R"); |
1963 | |
1964 | TGeoTranslation *side1Tr = new TGeoTranslation("ITSsddCPStr1R",0, |
b7943f00 |
1965 | - fgkCoolPipeSuppAxeDist |
1966 | + fgkCoolPipeSuppWidthExt/2., 0); |
c789ee28 |
1967 | side1Tr->RegisterYourself(); |
1968 | TGeoTranslation *side2Tr = new TGeoTranslation("ITSsddCPStr2R",0, |
b7943f00 |
1969 | - fgkCoolPipeSuppAxeDist |
1970 | + fgkCoolPipeSuppWidthExt*3/2. |
1971 | + fgkCoolPipeSuppWidthIn, 0); |
c789ee28 |
1972 | side2Tr->RegisterYourself(); |
1973 | |
1974 | TGeoBBox *middle = new TGeoBBox("ITSsddCPSmiddleR", |
b7943f00 |
1975 | (fgkCoolPipeSuppMaxLength/2. |
1976 | - fgkCoolPipeSuppSlitL)/2., |
1977 | fgkCoolPipeSuppWidthIn/2., |
1978 | fgkCoolPipeSuppHeight/2.); |
c789ee28 |
1979 | TGeoTranslation *middleTr = |
1980 | new TGeoTranslation("ITSsddCPStr3R", |
b7943f00 |
1981 | -( fgkCoolPipeSuppMaxLength/2. |
1982 | -fgkCoolPipeSuppSlitL)/2., |
1983 | -fgkCoolPipeSuppAxeDist + fgkCoolPipeSuppWidthExt |
1984 | + fgkCoolPipeSuppWidthIn/2.,0); |
c789ee28 |
1985 | middleTr->RegisterYourself(); |
1986 | |
1987 | TGeoBBox *axeBox = new TGeoBBox("ITSsddCPSaxeBoxR", |
b7943f00 |
1988 | fgkCoolPipeSuppTongW/4., |
1989 | (fgkCoolPipeSuppFulWidth |
1990 | - 2*fgkCoolPipeSuppWidthExt |
1991 | - fgkCoolPipeSuppWidthIn)/2, |
1992 | fgkCoolPipeSuppHeight/2.); |
c789ee28 |
1993 | |
1994 | TGeoTranslation *axeBoxTr = new TGeoTranslation("ITSsddCPSAxBoxTrR", |
b7943f00 |
1995 | - fgkCoolPipeSuppTongW/4., |
1996 | - fgkCoolPipeSuppAxeDist |
1997 | + fgkCoolPipeSuppFulWidth |
c789ee28 |
1998 | - axeBox->GetDY(),0); |
1999 | axeBoxTr->RegisterYourself(); |
2000 | |
b7943f00 |
2001 | TGeoTube *axe = new TGeoTube("ITSsddCPSaxeR",0,fgkCoolPipeSuppHoleDiam/2., |
2002 | fgkCoolPipeSuppTongW/4.); |
531d6cdc |
2003 | |
2004 | TGeoRotation *axeRot = new TGeoRotation("ITSsddCPSaxeRotR",90,90,0); |
c789ee28 |
2005 | TGeoCombiTrans *axeTrans = new TGeoCombiTrans("ITSsddCPSaxeTrR", |
531d6cdc |
2006 | -fgkCoolPipeSuppTongW/4.,0,0,axeRot); |
c789ee28 |
2007 | axeTrans->RegisterYourself(); |
531d6cdc |
2008 | //delete axeRot; |
c789ee28 |
2009 | |
2010 | if(GetDebug(3)){ |
2011 | middle->InspectShape(); |
2012 | axe->InspectShape(); |
2013 | }; |
2014 | |
2015 | TGeoCompositeShape *coolPipeSuppShape = new TGeoCompositeShape( |
2016 | "ITSsddCoolPipeSuppShapeR", |
2017 | "ITSsddCPSmiddleR:ITSsddCPStr3R" |
2018 | "+ITSsddCPSside1R:ITSsddCPStr1R" |
2019 | "+ITSsddCPSside1R:ITSsddCPStr2R" |
2020 | "+ITSsddCPSaxeBoxR:ITSsddCPSAxBoxTrR" |
2021 | "-ITSsddCPSaxeR:ITSsddCPSaxeTrR"); |
b7943f00 |
2022 | |
2023 | TGeoMedium *rytonSDD = GetMedium("ITSsddCarbonM55J"); //medium = ryton ? To code !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
c789ee28 |
2024 | TGeoVolume *coolPipeSupp = new TGeoVolume( "ITSsddCoolPipeSupportR", |
2025 | coolPipeSuppShape, rytonSDD); |
2026 | coolPipeSupp->SetLineColor(fColorRyton); |
2027 | |
2028 | return coolPipeSupp; |
531d6cdc |
2029 | } |
c789ee28 |
2030 | |
2031 | //________________________________________________________________________ |
2032 | TGeoVolume* AliITSv11GeometrySDD::CreateBaseThermalBridge() { |
2033 | // ALR 0752/8 |
2034 | |
b7943f00 |
2035 | Double_t dy = fgkBTBaxisAtoBase - fgkRadiusBminBTB - fgkBTBthick; |
c789ee28 |
2036 | |
b7943f00 |
2037 | Double_t base1width = fgkBTBwidth - fgkBTBaxisAtoBottom - fgkRadiusBminBTB |
2038 | - (fgkRadiusAminBTB+fgkBTBthick); |
c789ee28 |
2039 | TGeoBBox *base1 = new TGeoBBox( "ITSsddBTBbase1", base1width/2., |
b7943f00 |
2040 | fgkBTBthick/2., fgkBTBlength/2.); |
c789ee28 |
2041 | TGeoTranslation *base1Tr = new TGeoTranslation("ITSsddBTBtr1", |
b7943f00 |
2042 | fgkBTBaxisAtoBottom-fgkBTBwidth+base1width/2., |
2043 | -(fgkBTBaxisAtoBase-fgkBTBthick/2.), 0); |
c789ee28 |
2044 | base1Tr->RegisterYourself(); |
2045 | |
b7943f00 |
2046 | Double_t base2width = fgkBTBaxisAtoBottom - fgkRadiusAminBTB - fgkBTBthick |
2047 | - fgkRadiusBminBTB; |
c789ee28 |
2048 | TGeoBBox *base2 = new TGeoBBox( "ITSsddBTBbase2", base2width/2., |
b7943f00 |
2049 | fgkBTBthick/2., fgkBTBlength/2.); |
c789ee28 |
2050 | TGeoTranslation *base2Tr = new TGeoTranslation("ITSsddBTBtr2", |
b7943f00 |
2051 | fgkBTBaxisAtoBottom - base2width/2., |
2052 | -(fgkBTBaxisAtoBase-fgkBTBthick/2.), 0); |
c789ee28 |
2053 | base2Tr->RegisterYourself(); |
2054 | |
2055 | TGeoBBox *side = new TGeoBBox( "ITSsddBTBside", |
b7943f00 |
2056 | fgkBTBthick/2., dy/2., fgkBTBlength/2.); |
c789ee28 |
2057 | TGeoTranslation *sideTr1 = new TGeoTranslation("ITSsddBTBsideTr1", |
b7943f00 |
2058 | -fgkRadiusAminBTB-fgkBTBthick/2., -dy/2., 0); |
c789ee28 |
2059 | TGeoTranslation *sideTr2 = new TGeoTranslation("ITSsddBTBsideTr2", |
b7943f00 |
2060 | fgkRadiusAminBTB+fgkBTBthick/2., -dy/2., 0); |
c789ee28 |
2061 | sideTr1->RegisterYourself(); |
2062 | sideTr2->RegisterYourself(); |
2063 | |
b7943f00 |
2064 | TGeoBBox *hole = new TGeoBBox( "ITSsddBTBhole", fgkBTBHolewidth/2., |
2065 | fgkBTBthick/2., fgkBTBHoleLength/2.); |
c789ee28 |
2066 | TGeoTranslation *holeTr1 = new TGeoTranslation("ITSsddBTBholeTr1", |
b7943f00 |
2067 | - fgkBTBHoleRefX + fgkBTBHolewidth/2., |
2068 | - (fgkBTBaxisAtoBase-fgkBTBthick/2.), |
2069 | fgkBTBHoleRefY+(fgkBTBHoleLength-fgkBTBlength)/2.); |
c789ee28 |
2070 | TGeoTranslation *holeTr2 = new TGeoTranslation("ITSsddBTBholeTr2", |
b7943f00 |
2071 | - fgkBTBHoleRefX + fgkBTBHolewidth/2., |
2072 | - (fgkBTBaxisAtoBase-fgkBTBthick/2.), |
2073 | - fgkBTBHoleRefY-(fgkBTBHoleLength-fgkBTBlength)/2.); |
c789ee28 |
2074 | holeTr1->RegisterYourself(); |
2075 | holeTr2->RegisterYourself(); |
2076 | |
b7943f00 |
2077 | Double_t radiusAmaxBTB = fgkRadiusAminBTB + fgkBTBthick; |
c789ee28 |
2078 | TGeoTubeSeg *mainAxis = new TGeoTubeSeg( "ITSsddBTBmainAxis", |
b7943f00 |
2079 | fgkRadiusAminBTB, radiusAmaxBTB, |
2080 | fgkBTBlength/2., 0., 180.); |
c789ee28 |
2081 | TGeoTubeSeg *round1 = new TGeoTubeSeg( "ITSsddBTBround1", |
b7943f00 |
2082 | fgkRadiusBminBTB, fgkRadiusBminBTB+fgkBTBthick, |
2083 | fgkBTBlength/2., 270., 360.); |
c789ee28 |
2084 | TGeoTranslation *roundTr1 = new TGeoTranslation("ITSsddBTBround1Tr", |
b7943f00 |
2085 | -(fgkRadiusAminBTB+fgkBTBthick+fgkRadiusBminBTB), |
c789ee28 |
2086 | -dy, 0); |
2087 | roundTr1->RegisterYourself(); |
2088 | |
2089 | TGeoTubeSeg *round2 = new TGeoTubeSeg( "ITSsddBTBround2", |
b7943f00 |
2090 | fgkRadiusBminBTB, fgkRadiusBminBTB+fgkBTBthick, |
2091 | fgkBTBlength/2., 180., 270.); |
c789ee28 |
2092 | TGeoTranslation *roundTr2 = new TGeoTranslation("ITSsddBTBround2Tr", |
b7943f00 |
2093 | (fgkRadiusAminBTB+fgkBTBthick+fgkRadiusBminBTB), |
c789ee28 |
2094 | -dy, 0); |
2095 | roundTr2->RegisterYourself(); |
2096 | |
2097 | TGeoCompositeShape *sBaseThermalBridge = new TGeoCompositeShape( |
2098 | "ITSsddBaseThermalBridgeShape", |
2099 | "ITSsddBTBbase1:ITSsddBTBtr1" |
2100 | "+ ITSsddBTBbase2:ITSsddBTBtr2" |
2101 | "+ ITSsddBTBround1:ITSsddBTBround1Tr" |
2102 | "+ ITSsddBTBround2:ITSsddBTBround2Tr" |
2103 | "+ ITSsddBTBside:ITSsddBTBsideTr1" |
2104 | "+ ITSsddBTBside:ITSsddBTBsideTr2" |
2105 | "- ITSsddBTBhole:ITSsddBTBholeTr1" |
2106 | "- ITSsddBTBhole:ITSsddBTBholeTr2" |
2107 | "+ ITSsddBTBmainAxis"); |
2108 | |
2109 | if(GetDebug(3)){// Remove compiler warning. |
2110 | base1->InspectShape(); |
2111 | base2->InspectShape(); |
2112 | side->InspectShape(); |
2113 | hole->InspectShape(); |
2114 | mainAxis->InspectShape(); |
2115 | round1->InspectShape(); |
2116 | round2->InspectShape(); |
2117 | }; |
2118 | |
b7943f00 |
2119 | TGeoMedium *carbonFiberLadderStruct = GetMedium("ITSsddCarbonM55J"); |
c789ee28 |
2120 | TGeoVolume *vBaseThermalBridge = new TGeoVolume( "ITSsddBaseThermalBridge", |
2121 | sBaseThermalBridge, |
2122 | carbonFiberLadderStruct); |
2123 | |
2124 | vBaseThermalBridge->SetLineColor(fColorCarbonFiber); |
2125 | return vBaseThermalBridge; |
531d6cdc |
2126 | } |
c789ee28 |
2127 | |
2128 | |
c789ee28 |
2129 | //________________________________________________________________________ |
b7943f00 |
2130 | TGeoVolume* AliITSv11GeometrySDD::CreateEndLadder(Int_t iLay) { |
2131 | // |
2132 | // Return a box volume containing a end of a CF ladder. |
2133 | // |
2134 | |
2135 | TGeoMedium *airSDD = GetMedium("ITSair"); |
2136 | TGeoMedium *carbonFiberLadderStruct = GetMedium("ITSsddCarbonM55J"); |
2137 | |
2138 | Double_t length = (fgkLay3LadderLength-fgkLay3Ndet*fgkSegmentLength)/2.; |
2139 | Double_t coolPipeSuppH = fgkLay3CoolPipeSuppH; |
2140 | Double_t underSegDH = fLay3LadderUnderSegDH; |
2141 | if (iLay==3) { |
2142 | } else if (iLay==4) { |
2143 | length = (fgkLay4LadderLength-fgkLay4Ndet*fgkSegmentLength)/2.; |
2144 | coolPipeSuppH = fgkLay4CoolPipeSuppH; |
2145 | underSegDH = fLay4LadderUnderSegDH; |
2146 | } else { |
2147 | printf("error in AliITSv11GeometrySDD::CreateEndLadder: Wrong layer"); |
2148 | return 0; |
2149 | }; |
c789ee28 |
2150 | |
b7943f00 |
2151 | Double_t tDY = (- fgkLadderSegBoxDH/2 //space left on top of the ladder |
2152 | + underSegDH/2); //space under ladder segment |
2153 | // here tDY is not the same as for the segment because the end ladder |
2154 | // does not have a space under it, inside the general ladder volume. |
2155 | Double_t segmentLength = fgkSegmentLength; |
2156 | Double_t topCornerLength = fgkSegmentLength/2.-fgkLay4LaddTopCornerEnd; |
2157 | |
2158 | TGeoBBox *endBox = new TGeoBBox("ITSsddEndLaddBox", |
2159 | (fgkLadderWidth)/2, |
2160 | fgkLadderHeight/2+fgkLadderSegBoxDH/2+underSegDH/2, |
2161 | length/2); |
2162 | TGeoVolume *virtualEnd = new TGeoVolume("ITSsddEnd",endBox, airSDD); |
2163 | |
2164 | //********************************** |
2165 | // coding real matter : |
2166 | //********************************** |
2167 | Double_t triangleHeight = fgkLadderHeight - fgkLadderBeamRadius; |
2168 | Double_t halfTheta = TMath::ATan( 0.5*fgkLadderWidth/triangleHeight ); |
2169 | Double_t beta = (TMath::Pi()-2.*halfTheta)/4.; |
2170 | Double_t alpha = TMath::Pi()*3./4. - halfTheta/2.; |
2171 | |
2172 | //--- The 3 V shape corners of the Carbon Fiber Ladder |
2173 | //--- the top V |
2174 | TGeoArb8 *cfLaddTop1 = CreateLadderSide(topCornerLength/2., halfTheta, -1, |
2175 | fgkLadderLa, fgkLadderHa, fgkLadderl); |
2176 | TGeoVolume *cfLaddTopVol1 = new TGeoVolume("ITSsddCFladdTopCornerVol1", |
2177 | cfLaddTop1,carbonFiberLadderStruct); |
2178 | cfLaddTopVol1->SetLineColor(fColorCarbonFiber); |
2179 | TGeoArb8 *cfLaddTop2 = CreateLadderSide( topCornerLength/2., halfTheta, 1, |
2180 | fgkLadderLa, fgkLadderHa, fgkLadderl); |
2181 | TGeoVolume *cfLaddTopVol2 = new TGeoVolume("ITSsddCFladdTopCornerV2", |
2182 | cfLaddTop2,carbonFiberLadderStruct); |
2183 | cfLaddTopVol2->SetLineColor(fColorCarbonFiber); |
2184 | TGeoTranslation *trTop1 = new TGeoTranslation(0, fgkLadderHeight/2+tDY, |
2185 | -(length-topCornerLength)/2.); |
2186 | virtualEnd->AddNode(cfLaddTopVol1, 1, trTop1); |
2187 | virtualEnd->AddNode(cfLaddTopVol2, 1, trTop1); |
2188 | |
2189 | //--- the 2 side V |
2190 | TGeoArb8 *cfLaddSide1 = CreateLadderSide( length/2., beta, -1, |
2191 | fgkLadderLb, fgkLadderHb, fgkLadderl); |
2192 | TGeoVolume *cfLaddSideVol1 = new TGeoVolume("ITSsddCFladdSideCornerV1", |
2193 | cfLaddSide1,carbonFiberLadderStruct); |
2194 | cfLaddSideVol1->SetLineColor(fColorCarbonFiber); |
2195 | TGeoArb8 *cfLaddSide2 = CreateLadderSide( length/2., beta, 1, |
2196 | fgkLadderLb, fgkLadderHb, fgkLadderl); |
2197 | TGeoVolume *cfLaddSideVol2 = new TGeoVolume("ITSsddCFladdSideCornerV2", |
2198 | cfLaddSide2,carbonFiberLadderStruct); |
2199 | cfLaddSideVol2->SetLineColor(fColorCarbonFiber); |
2200 | Double_t dYTranslation = ( fgkLadderHeight/2. - 0.5*fgkLadderWidth* |
2201 | TMath::Tan(beta) - fgkLadderBeamRadius ); |
2202 | |
2203 | // because center of the triangle doesn't correspond to virtual vol. center |
2204 | Double_t distCenterSideDown = 0.5*fgkLadderWidth/TMath::Cos(beta); |
2205 | TGeoCombiTrans *ctSideR = CreateCombiTrans("", distCenterSideDown, 0, |
2206 | alpha*TMath::RadToDeg()); |
2207 | AddTranslationToCombiTrans(ctSideR, 0, -dYTranslation+tDY, 0); |
2208 | TGeoCombiTrans *ctSideL = CreateCombiTrans("", distCenterSideDown, 0, |
2209 | -alpha*TMath::RadToDeg()); |
2210 | AddTranslationToCombiTrans(ctSideL, 0, -dYTranslation+tDY, 0); |
2211 | virtualEnd->AddNode(cfLaddSideVol1, 1, ctSideR); |
2212 | virtualEnd->AddNode(cfLaddSideVol2, 1, ctSideR); |
2213 | virtualEnd->AddNode(cfLaddSideVol1, 2, ctSideL); |
2214 | virtualEnd->AddNode(cfLaddSideVol2, 2, ctSideL); |
2215 | |
2216 | //--- The beams |
2217 | // Beams on the sides |
2218 | Double_t beamPhiPrime = TMath::ASin(1./TMath::Sqrt( (1+TMath::Sin(2*beta)* |
2219 | TMath::Sin(2*beta)/(TanD(fgkBeamSidePhi)*TanD(fgkBeamSidePhi))) )); |
2220 | |
2221 | //Euler rotation : about Z, then new X, then new Z |
531d6cdc |
2222 | TGeoRotation *beamRot1 = new TGeoRotation("", 90-2.*beta*TMath::RadToDeg(), |
b7943f00 |
2223 | -beamPhiPrime*TMath::RadToDeg(), -90); |
531d6cdc |
2224 | TGeoRotation *beamRot2 = new TGeoRotation("", 90-2.*beta*TMath::RadToDeg(), |
b7943f00 |
2225 | beamPhiPrime*TMath::RadToDeg(), -90); |
531d6cdc |
2226 | TGeoRotation *beamRot3 = new TGeoRotation("", 90+2.*beta*TMath::RadToDeg(), |
b7943f00 |
2227 | beamPhiPrime*TMath::RadToDeg(), -90); |
531d6cdc |
2228 | TGeoRotation *beamRot4 = new TGeoRotation("", 90+2.*beta*TMath::RadToDeg(), |
b7943f00 |
2229 | -beamPhiPrime*TMath::RadToDeg(), -90); |
2230 | TGeoCombiTrans *beamTransf1 = new TGeoCombiTrans(0.5*triangleHeight* |
2231 | TMath::Tan(halfTheta), |
2232 | fgkLadderBeamRadius/2. + tDY, |
531d6cdc |
2233 | -length/2 + segmentLength/8, beamRot1); |
b7943f00 |
2234 | TGeoCombiTrans *beamTransf3 = new TGeoCombiTrans( 0.5*triangleHeight* |
2235 | TMath::Tan(halfTheta), |
2236 | fgkLadderBeamRadius/2.+tDY, |
531d6cdc |
2237 | -length/2 + 3*segmentLength/8, beamRot2); |
b7943f00 |
2238 | TGeoCombiTrans *beamTransf5 = new TGeoCombiTrans(-0.5*triangleHeight* |
2239 | TMath::Tan(halfTheta), |
2240 | fgkLadderBeamRadius/2.+tDY, |
531d6cdc |
2241 | -length/2 + segmentLength/8, beamRot3); |
b7943f00 |
2242 | TGeoCombiTrans *beamTransf7 = new TGeoCombiTrans(-0.5*triangleHeight* |
2243 | TMath::Tan(halfTheta), |
2244 | fgkLadderBeamRadius/2. + tDY, |
531d6cdc |
2245 | -length/2+3*segmentLength/8, beamRot4); |
b7943f00 |
2246 | |
2247 | virtualEnd->AddNode(fLaddSegCommonVol[6], 1, beamTransf1); |
2248 | virtualEnd->AddNode(fLaddSegCommonVol[6], 2, beamTransf3); |
2249 | virtualEnd->AddNode(fLaddSegCommonVol[6], 3, beamTransf5); |
2250 | virtualEnd->AddNode(fLaddSegCommonVol[6], 4, beamTransf7); |
2251 | |
2252 | //--- Beams of the bottom |
2253 | TGeoTubeSeg *bottomBeam1 = new TGeoTubeSeg(0, fgkLadderBeamRadius, |
2254 | fgkLadderWidth/2.-fgkLadderLb/3, 0, 180); |
2255 | TGeoVolume *bottomBeam1Vol = new TGeoVolume("ITSsddBottomBeam1Vol", |
2256 | bottomBeam1, carbonFiberLadderStruct); |
2257 | bottomBeam1Vol->SetLineColor(fColorCarbonFiber); |
2258 | |
531d6cdc |
2259 | TGeoRotation *bottomBeamRot1 = new TGeoRotation("",90, 90, 90); |
b7943f00 |
2260 | TGeoCombiTrans *bottomBeamTransf1 = new TGeoCombiTrans(0, |
2261 | -(fgkLadderHeight/2-fgkLadderBeamRadius)+tDY, |
531d6cdc |
2262 | -length/2+fgkSegmentLength/2, bottomBeamRot1); |
b7943f00 |
2263 | virtualEnd->AddNode(bottomBeam1Vol, 1, bottomBeamTransf1); |
2264 | TGeoTubeSeg *bottomBeam2 = new TGeoTubeSeg(0, fgkLadderBeamRadius, |
2265 | fgkLadderWidth/2.-fgkLadderLb/3, 0, 90); |
2266 | TGeoVolume *bottomBeam2Vol = new TGeoVolume("ITSsddBottomBeam2Vol", |
2267 | bottomBeam2, carbonFiberLadderStruct); |
2268 | bottomBeam2Vol->SetLineColor(fColorCarbonFiber); |
2269 | TGeoCombiTrans *bottomBeamTransf2 = new TGeoCombiTrans(0, |
531d6cdc |
2270 | -(fgkLadderHeight/2-fgkLadderBeamRadius)+tDY,-length/2,bottomBeamRot1); |
b7943f00 |
2271 | virtualEnd->AddNode(bottomBeam2Vol, 1, bottomBeamTransf2); |
2272 | |
2273 | //********************************** |
2274 | //the cooling pipe supports |
2275 | Double_t triangleCPaxeDist = fgkCoolPipeSuppAxeDist-fgkCoolPipeSuppWidthExt- |
fa4639a3 |
2276 | fgkCoolPipeSuppWidthIn+fgkLadderBeamRadius; |
b7943f00 |
2277 | |
2278 | Double_t coolPipeSuppL = TMath::Tan(halfTheta)* |
fa4639a3 |
2279 | (triangleHeight+triangleCPaxeDist/ |
2280 | TMath::Sin(halfTheta)-coolPipeSuppH); |
b7943f00 |
2281 | |
2282 | if (fAddCoolingSyst) { |
531d6cdc |
2283 | TGeoRotation *rotCPS2 = new TGeoRotation("",-halfTheta*TMath::RadToDeg(),-90, 90); |
2284 | TGeoRotation *rotCPS1 = new TGeoRotation("", halfTheta*TMath::RadToDeg(),-90,-90); |
b7943f00 |
2285 | TGeoCombiTrans *transCPS1 = new TGeoCombiTrans(coolPipeSuppL, |
fa4639a3 |
2286 | -fgkLadderHeight/2.+ tDY + |
b7943f00 |
2287 | coolPipeSuppH+fgkLadderBeamRadius, |
531d6cdc |
2288 | -length/2., rotCPS1); |
b7943f00 |
2289 | TGeoCombiTrans *transCPS4 = new TGeoCombiTrans(-coolPipeSuppL, |
fa4639a3 |
2290 | -fgkLadderHeight/2.+ tDY + |
b7943f00 |
2291 | coolPipeSuppH+fgkLadderBeamRadius, |
531d6cdc |
2292 | -length/2., rotCPS2); |
b7943f00 |
2293 | |
2294 | virtualEnd->AddNode(fCoolPipeSupportL, 1, transCPS1); |
2295 | virtualEnd->AddNode(fCoolPipeSupportR, 1, transCPS4); |
2296 | }; |
2297 | |
2298 | //********************************** |
2299 | if(GetDebug(1)) virtualEnd->CheckOverlaps(0.01); |
fa4639a3 |
2300 | virtualEnd->SetVisibility(kFALSE); |
b7943f00 |
2301 | return virtualEnd; |
531d6cdc |
2302 | } |
c789ee28 |
2303 | |
2304 | |
db486a6e |
2305 | //________________________________________________________________________ |
2306 | TGeoVolume* AliITSv11GeometrySDD::CreateSDDsensor() { |
b7943f00 |
2307 | // |
2308 | // return a box containing the SDD sensor |
2309 | // |
2310 | |
2311 | TGeoMedium *airSDD = GetMedium("ITSair"); |
2312 | TGeoMedium *siliconSDD = GetMedium("ITSsddSi"); |
2313 | TGeoMedium *alSDD = GetMedium("ITSal"); |
2314 | TGeoMedium *polyhamideSDD = GetMedium("ITSsddKAPTON_POLYCH2"); |
2315 | TGeoMedium *glassSDD = GetMedium("ITSsddSi"); // To code !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
2316 | |
2317 | |
2318 | Double_t rWraping = fgkWaferThickness/2+fgkWaHVcableAlThick+fgkWaHVcablePolyThick; |
2319 | Double_t witdhCableBox = (fgkWaHVcableWitdh - TMath::Pi()*rWraping)/2; |
2320 | |
2321 | Double_t sensoxBoxLength = ( fgkWaferLength + |
2322 | 2*(rWraping+witdhCableBox-fgkWaHVcableDW) ); |
2323 | // Makes life easier to include the space for the WA HV cable on both sides |
2324 | Double_t sensoxBoxThick = fgkWaferThickness + |
2325 | 2*(fgkWaHVcableAlThick+fgkWaHVcablePolyThick); |
2326 | |
2327 | TGeoBBox *box = new TGeoBBox("ITSsddSensorBox", |
2328 | fgkWaferWidth/2, sensoxBoxThick/2, sensoxBoxLength/2); |
2329 | TGeoVolume *virtualSensor = new TGeoVolume("ITSsddSensor",box,airSDD); |
2330 | |
2331 | //**************************** |
2332 | // silicon wafer |
2333 | //**************************** |
2334 | if (fAddSensors) { |
c789ee28 |
2335 | TGeoBBox *waferShape = new TGeoBBox("ITSsddWaferShape", |
b7943f00 |
2336 | fgkWaferWidth/2, fgkWaferThickness/2, fgkWaferLength/2); |
2337 | TGeoVolume *wafer = new TGeoVolume("ITSsddWafer", waferShape, siliconSDD); |
2338 | wafer->SetLineColor(fColorSilicon); |
2339 | TGeoBBox *sensBox = new TGeoBBox("ITSsddSensorSensBox", |
2340 | fgkWaferWidthSens/2,fgkWaferThickSens/2,fgkWaferLengthSens/2); |
2341 | TGeoVolume *sensVol=new TGeoVolume(fgSDDsensitiveVolName,sensBox,siliconSDD); |
2342 | sensVol->SetLineColor(fColorSilicon); |
2343 | |
2344 | wafer->AddNode(sensVol, 1, 0); |
2345 | virtualSensor->AddNode(wafer, 1, 0); |
2346 | }; |
2347 | |
2348 | //**************************** |
2349 | // glass |
2350 | //**************************** |
2351 | TGeoBBox *glass = new TGeoBBox("ITSsddGlassBox", fgkSensorGlassLX/2, |
2352 | fgkSensorGlassLY/2, fgkSensorGlassLZ/2); |
2353 | TGeoVolume *vGlass = new TGeoVolume("ITSsddGlass",glass, glassSDD); |
2354 | vGlass->SetLineColor(fColorGlass); |
2355 | TGeoTranslation *glassTr1 = new TGeoTranslation("",fgkGlassDXOnSensor, |
2356 | fgkWaferThickness/2+fgkSensorGlassLY/2, |
2357 | fgkGlassDZOnSensor); |
2358 | TGeoTranslation *glassTr2 = new TGeoTranslation("",-fgkGlassDXOnSensor, |
2359 | fgkWaferThickness/2+fgkSensorGlassLY/2, |
2360 | fgkGlassDZOnSensor); |
2361 | TGeoTranslation *glassTr3 = new TGeoTranslation("",fgkGlassDXOnSensor, |
2362 | fgkWaferThickness/2+fgkSensorGlassLY/2, |
2363 | -fgkGlassDZOnSensor); |
2364 | TGeoTranslation *glassTr4 = new TGeoTranslation("",-fgkGlassDXOnSensor, |
2365 | fgkWaferThickness/2+fgkSensorGlassLY/2, |
2366 | -fgkGlassDZOnSensor); |
2367 | virtualSensor->AddNode(vGlass, 1, glassTr1); |
2368 | virtualSensor->AddNode(vGlass, 2, glassTr2); |
2369 | virtualSensor->AddNode(vGlass, 3, glassTr3); |
2370 | virtualSensor->AddNode(vGlass, 4, glassTr4); |
2371 | |
2372 | //**************************** |
2373 | // Wrap-around cable |
2374 | //**************************** |
2375 | if (fAddHVcables) { |
2376 | AliITSv11GeomCableFlat waHVCable("ITSsddWaHVCableU",witdhCableBox, |
2377 | fgkWaHVcableAlThick+fgkWaHVcablePolyThick); |
2378 | waHVCable.SetNLayers(2); |
2379 | waHVCable.SetLayer(0, fgkWaHVcablePolyThick,polyhamideSDD,fColorPolyhamide); |
2380 | waHVCable.SetLayer(1, fgkWaHVcableAlThick, alSDD, fColorAl); |
2381 | waHVCable.SetInitialNode(virtualSensor); |
2382 | |
2383 | Double_t x1[3], x2[3], vX[3] = {1,0,0}; |
2384 | x1[0] = -fgkWaHVcableLength/2; |
2385 | x2[0] = -x1[0]; |
2386 | x1[1] = (fgkWaferThickness + waHVCable.GetThickness())/2; |
2387 | x2[1] = x1[1]; |
2388 | x1[2] = fgkWaferLength/2+waHVCable.GetWidth()/2-fgkWaHVcableDW; |
2389 | x2[2] = x1[2]; |
2390 | |
2391 | waHVCable.AddCheckPoint(virtualSensor, 0, x1, vX); |
2392 | waHVCable.AddCheckPoint(virtualSensor, 1, x2, vX); |
2393 | waHVCable.CreateAndInsertCableSegment(1,-90); |
2394 | x1[1] = -x1[1]; |
2395 | x2[1] = x1[1]; |
2396 | waHVCable.SetName("ITSsddWaHVCableD"); |
2397 | waHVCable.ResetPoints(); |
2398 | waHVCable.AddCheckPoint(virtualSensor, 0, x1, vX); |
2399 | waHVCable.AddCheckPoint(virtualSensor, 1, x2, vX); |
2400 | waHVCable.CreateAndInsertCableSegment(1, 90); |
2401 | |
2402 | AliITSv11GeomCableRound waHVCableFold("ITSsddWaHVCableFold", |
2403 | rWraping); |
2404 | waHVCableFold.SetPhi(180,360); |
2405 | waHVCableFold.SetNLayers(2); |
2406 | waHVCableFold.SetLayer(0, fgkWaferThickness/2+fgkWaHVcablePolyThick, |
2407 | polyhamideSDD, fColorPolyhamide); |
2408 | waHVCableFold.SetLayer(1, fgkWaHVcableAlThick, alSDD, fColorAl); |
2409 | waHVCableFold.SetInitialNode(virtualSensor); |
2410 | x1[1] = 0; |
2411 | x2[1] = 0; |
2412 | x1[2] = fgkWaferLength/2-fgkWaHVcableDW+witdhCableBox; |
2413 | x2[2] = x1[2]; |
2414 | waHVCableFold.AddCheckPoint(virtualSensor, 0, x1, vX); |
2415 | waHVCableFold.AddCheckPoint(virtualSensor, 1, x2, vX); |
2416 | waHVCableFold.CreateAndInsertCableSegment(1); |
2417 | |
2418 | //**************************** |
2419 | // transition cable |
2420 | //**************************** |
2421 | Double_t headRadius = (fgkTransitHVHeadLX*fgkTransitHVHeadLX/4.+ |
2422 | fgkTransitHVHeadLZ*fgkTransitHVHeadLZ) |
2423 | /(2.*fgkTransitHVHeadLZ); |
2424 | Double_t theta = TMath::ATan2(fgkTransitHVHeadLX/2, |
2425 | headRadius-fgkTransitHVHeadLZ) |
2426 | *TMath::RadToDeg(); |
2427 | |
2428 | TGeoTubeSeg *headPoly = new TGeoTubeSeg(0,headRadius, |
2429 | fgkTransitHVPolyThick/2, |
2430 | 90-theta,90+theta); |
2431 | headPoly->SetName("headPoly"); |
2432 | TGeoTranslation *headPolyTr = new TGeoTranslation(0,0, |
2433 | -fgkTransitHVPolyThick/2); |
2434 | headPolyTr->SetName("headPolyTr"); |
2435 | headPolyTr->RegisterYourself(); |
2436 | |
2437 | TGeoTubeSeg *headAl = new TGeoTubeSeg(0,headRadius, |
2438 | fgkTransitHVAlThick/2, |
2439 | 90-theta,90+theta); |
2440 | headAl->SetName("headAl"); |
2441 | TGeoTranslation *headAlTr = new TGeoTranslation(0,0, |
2442 | -fgkTransitHVPolyThick |
2443 | -fgkTransitHVAlThick/2); |
2444 | headAlTr->SetName("headAlTr"); |
2445 | headAlTr->RegisterYourself(); |
2446 | |
2447 | TGeoBBox *cache = new TGeoBBox(fgkTransitHVHeadLX/2, |
2448 | (headRadius-fgkTransitHVHeadLZ)/2, |
2449 | (fgkTransitHVPolyThick+fgkTransitHVAlThick)/2); |
2450 | cache->SetName("cache"); |
2451 | |
2452 | TGeoTranslation *headCacheTr = new TGeoTranslation(0, |
2453 | (headRadius-fgkTransitHVHeadLZ)/2, |
2454 | -(fgkTransitHVPolyThick |
2455 | +fgkTransitHVAlThick)/2); |
2456 | headCacheTr->SetName("cacheTr"); |
2457 | headCacheTr->RegisterYourself(); |
2458 | |
2459 | TGeoCompositeShape *headPolyComp = new TGeoCompositeShape( |
2460 | "headPoly:headPolyTr-cache:cacheTr"); |
2461 | TGeoVolume *vHeadPolyComp = new TGeoVolume( |
2462 | "ITSsddHVtransitHeadPoly",headPolyComp, polyhamideSDD); |
2463 | vHeadPolyComp->SetLineColor(fColorPolyhamide); |
2464 | TGeoCompositeShape *headAlComp = new TGeoCompositeShape( |
2465 | "headAl:headAlTr-cache:cacheTr"); |
2466 | TGeoVolume *vHeadAlComp = new TGeoVolume( |
2467 | "ITSsddHVtransitHeadAl",headAlComp, alSDD); |
2468 | vHeadAlComp->SetLineColor(fColorAl); |
2469 | |
2470 | |
531d6cdc |
2471 | // TGeoRotation rotHead("",0,90,0); |
2472 | // TGeoCombiTrans *rotHeadTr = new TGeoCombiTrans(0,fgkWaferThickness/2, |
2473 | // -headRadius+fgkTransitHVHeadLZ+fgkTransitHVBondingLZ/2, |
2474 | // &rotHead); |
2475 | TGeoRotation *rotHead = new TGeoRotation("",0,90,0); |
b7943f00 |
2476 | TGeoCombiTrans *rotHeadTr = new TGeoCombiTrans(0,fgkWaferThickness/2, |
2477 | -headRadius+fgkTransitHVHeadLZ+fgkTransitHVBondingLZ/2, |
531d6cdc |
2478 | rotHead); |
2479 | |
b7943f00 |
2480 | virtualSensor->AddNode(vHeadPolyComp,1,rotHeadTr); |
2481 | virtualSensor->AddNode(vHeadAlComp,1,rotHeadTr); |
2482 | |
2483 | //--- |
2484 | AliITSv11GeomCableFlat transitHVCable("ITSsddHVtransitCenter", |
2485 | fgkTransitHVBondingLZ, |
2486 | fgkTransitHVPolyThick+fgkTransitHVAlThick); |
2487 | transitHVCable.SetNLayers(2); |
2488 | transitHVCable.SetLayer(0, fgkTransitHVPolyThick,polyhamideSDD, |
2489 | fColorPolyhamide); |
2490 | transitHVCable.SetLayer(1, fgkTransitHVAlThick, alSDD, fColorAl); |
2491 | transitHVCable.SetInitialNode(virtualSensor); |
2492 | |
2493 | x1[0] = -fgkTransitHVHeadLX/2; |
2494 | x2[0] = -x1[0]; |
2495 | x1[1] = (fgkWaferThickness+fgkTransitHVPolyThick+fgkTransitHVAlThick)/2; |
2496 | x2[1] = x1[1]; |
2497 | x1[2] = 0; |
2498 | x2[2] = 0; |
2499 | transitHVCable.AddCheckPoint(virtualSensor, 0, x1, vX); |
2500 | transitHVCable.AddCheckPoint(virtualSensor, 1, x2, vX); |
2501 | transitHVCable.CreateAndInsertCableSegment(1,-90); |
2502 | transitHVCable.ResetPoints(); |
2503 | transitHVCable.SetName("ITSsddHVtransitTail"); |
2504 | transitHVCable.SetWidth(fgkTransitHVtailWidth); |
2505 | x1[0] = fgkTransitHVtailXpos; |
2506 | x2[0] = fgkTransitHVtailXpos; |
2507 | x1[2] = -fgkTransitHVBondingLZ/2; |
2508 | x2[2] = -fgkTransitHVBondingLZ/2-fgkTransitHVtailLength; |
2509 | Double_t vZ[3] = {0,0,1}; |
2510 | transitHVCable.AddCheckPoint(virtualSensor, 0, x1, vZ); |
2511 | transitHVCable.AddCheckPoint(virtualSensor, 1, x2, vZ); |
2512 | transitHVCable.CreateAndInsertCableSegment(1,0); |
2513 | |
2514 | //--- |
2515 | TGeoArb8 *sideLeft = new TGeoArb8( fgkTransitHVPolyThick/2 ); |
2516 | sideLeft->SetVertex(0, fgkTransitHVtailXpos+fgkTransitHVtailWidth/2,0); |
2517 | sideLeft->SetVertex(1, fgkTransitHVtailXpos+fgkTransitHVtailWidth/2, |
2518 | fgkTransitHVsideLZ); |
2519 | sideLeft->SetVertex(2, fgkTransitHVHeadLX/2, fgkTransitHVsideLeftZ); |
2520 | sideLeft->SetVertex(3, fgkTransitHVHeadLX/2, 0); |
2521 | sideLeft->SetVertex(4, fgkTransitHVtailXpos+fgkTransitHVtailWidth/2,0); |
2522 | sideLeft->SetVertex(5, fgkTransitHVtailXpos+fgkTransitHVtailWidth/2, |
2523 | fgkTransitHVsideLZ); |
2524 | sideLeft->SetVertex(6, fgkTransitHVHeadLX/2, fgkTransitHVsideLeftZ); |
2525 | sideLeft->SetVertex(7, fgkTransitHVHeadLX/2, 0); |
2526 | |
2527 | TGeoArb8 *sideLeftAl = new TGeoArb8( fgkTransitHVAlThick/2 ); |
2528 | sideLeftAl->SetVertex(0, fgkTransitHVtailXpos+fgkTransitHVtailWidth/2,0); |
2529 | sideLeftAl->SetVertex(1, fgkTransitHVtailXpos+fgkTransitHVtailWidth/2, |
2530 | fgkTransitHVsideLZ); |
2531 | sideLeftAl->SetVertex(2, fgkTransitHVHeadLX/2, fgkTransitHVsideLeftZ); |
2532 | sideLeftAl->SetVertex(3, fgkTransitHVHeadLX/2, 0); |
2533 | sideLeftAl->SetVertex(4, fgkTransitHVtailXpos+fgkTransitHVtailWidth/2,0); |
2534 | sideLeftAl->SetVertex(5, fgkTransitHVtailXpos+fgkTransitHVtailWidth/2, |
2535 | fgkTransitHVsideLZ); |
2536 | sideLeftAl->SetVertex(6, fgkTransitHVHeadLX/2, fgkTransitHVsideLeftZ); |
2537 | sideLeftAl->SetVertex(7, fgkTransitHVHeadLX/2, 0); |
2538 | |
2539 | TGeoArb8 *sideRight = new TGeoArb8( fgkTransitHVPolyThick/2 ); |
2540 | sideRight->SetVertex(0, fgkTransitHVtailXpos-fgkTransitHVtailWidth/2,0); |
2541 | sideRight->SetVertex(1, fgkTransitHVtailXpos-fgkTransitHVtailWidth/2, |
2542 | fgkTransitHVsideLZ); |
2543 | sideRight->SetVertex(2, -fgkTransitHVHeadLX/2, fgkTransitHVsideRightZ); |
2544 | sideRight->SetVertex(3, -fgkTransitHVHeadLX/2, 0); |
2545 | sideRight->SetVertex(4, fgkTransitHVtailXpos-fgkTransitHVtailWidth/2,0); |
2546 | sideRight->SetVertex(5, fgkTransitHVtailXpos-fgkTransitHVtailWidth/2, |
2547 | fgkTransitHVsideLZ); |
2548 | sideRight->SetVertex(6, -fgkTransitHVHeadLX/2, fgkTransitHVsideRightZ); |
2549 | sideRight->SetVertex(7, -fgkTransitHVHeadLX/2, 0); |
2550 | |
531d6cdc |
2551 | // TGeoRotation rotSide("",0,-90,0); |
2552 | // TGeoCombiTrans *sideRightTr = new TGeoCombiTrans(0, |
2553 | // (fgkWaferThickness+fgkTransitHVPolyThick)/2, |
2554 | // -fgkTransitHVBondingLZ/2,&rotSide); |
2555 | // TGeoCombiTrans *sideLeftTr = new TGeoCombiTrans(0, |
2556 | // (fgkWaferThickness+fgkTransitHVPolyThick)/2, |
2557 | // -fgkTransitHVBondingLZ/2, &rotSide); |
2558 | // TGeoCombiTrans *sideLeftAlTr = new TGeoCombiTrans(0, |
2559 | // fgkTransitHVPolyThick+(fgkWaferThickness+fgkTransitHVAlThick)/2, |
2560 | // -fgkTransitHVBondingLZ/2, &rotSide); |
2561 | TGeoRotation *rotSide = new TGeoRotation("",0,-90,0); |
b7943f00 |
2562 | TGeoCombiTrans *sideRightTr = new TGeoCombiTrans(0, |
2563 | (fgkWaferThickness+fgkTransitHVPolyThick)/2, |
531d6cdc |
2564 | -fgkTransitHVBondingLZ/2,rotSide); |
b7943f00 |
2565 | TGeoCombiTrans *sideLeftTr = new TGeoCombiTrans(0, |
2566 | (fgkWaferThickness+fgkTransitHVPolyThick)/2, |
531d6cdc |
2567 | -fgkTransitHVBondingLZ/2, rotSide); |
b7943f00 |
2568 | TGeoCombiTrans *sideLeftAlTr = new TGeoCombiTrans(0, |
2569 | fgkTransitHVPolyThick+(fgkWaferThickness+fgkTransitHVAlThick)/2, |
531d6cdc |
2570 | -fgkTransitHVBondingLZ/2, rotSide); |
2571 | |
b7943f00 |
2572 | |
2573 | TGeoVolume *vSideLeft = new TGeoVolume("ITSsddHVtransitSideLeft", |
2574 | sideLeft,polyhamideSDD); |
2575 | vSideLeft->SetLineColor(fColorPolyhamide); |
2576 | TGeoVolume *vSideLeftAl = new TGeoVolume("ITSsddHVtransitSideLeftAl", |
2577 | sideLeftAl,alSDD); |
2578 | vSideLeftAl->SetLineColor(fColorAl); |
2579 | TGeoVolume *vSideRight = new TGeoVolume("ITSsddHVtransitSideRight", |
2580 | sideRight,polyhamideSDD); |
2581 | vSideRight->SetLineColor(fColorPolyhamide); |
2582 | |
2583 | virtualSensor->AddNode(vSideLeft, 1, sideLeftTr); |
2584 | virtualSensor->AddNode(vSideLeftAl, 1, sideLeftAlTr); |
2585 | virtualSensor->AddNode(vSideRight, 1, sideRightTr); |
2586 | }; |
db486a6e |
2587 | |
b7943f00 |
2588 | //**************************** |
2589 | if(GetDebug(1)) virtualSensor->CheckOverlaps(0.01); |
2590 | virtualSensor->SetVisibility(kFALSE); |
2591 | return virtualSensor; |
531d6cdc |
2592 | } |
c789ee28 |
2593 | |
2594 | |
db486a6e |
2595 | //________________________________________________________________________ |
b7943f00 |
2596 | TGeoVolume *AliITSv11GeometrySDD::CreateDetectors(Int_t iLay) { |
2597 | // |
2598 | // return a box volume containing the detectors |
2599 | // |
2600 | |
2601 | TGeoMedium *airSDD = GetMedium("ITSair"); |
c789ee28 |
2602 | |
b7943f00 |
2603 | Int_t nDetectors = fgkLay3Ndet; |
2604 | Double_t ladderLength = fgkLay3LadderLength; |
2605 | Double_t *sensorZPos = fLay3sensorZPos; |
2606 | |
2607 | if (iLay==3) {} |
2608 | else if (iLay==4) { |
2609 | nDetectors = fgkLay4Ndet; |
2610 | ladderLength = fgkLay4LadderLength; |
2611 | sensorZPos = fLay4sensorZPos; |
2612 | } else { |
fa4639a3 |
2613 | printf("AliITSv11GeometrySDD::CreateDetectors: Error : Wrong layer"); |
b7943f00 |
2614 | }; |
2615 | |
2616 | char name[30]; |
2617 | Double_t volThickness = ( fgkLadWaferSep + 2*fgkWaferThickness + |
2618 | 2*(fgkWaHVcableAlThick+fgkWaHVcablePolyThick)); |
2619 | |
2620 | sprintf(name,"ITSsddDetBox%i",iLay); |
2621 | TGeoBBox *detBox = new TGeoBBox(name, fgkWaferWidth/2, volThickness/2, |
2622 | ladderLength*((nDetectors-0.5)/nDetectors)/2); |
2623 | TGeoVolume *virtualDet = new TGeoVolume("ITSsddLadd",detBox, airSDD); |
2624 | |
2625 | for (Int_t i=0; i<nDetectors; i++) { |
2626 | Double_t localZ = sensorZPos[i]; |
2627 | Double_t localY = fgkLadWaferSep/2+fgkWaferThickness/2; |
2628 | if (iLay==3) if (i%2!=0) localY = -localY; |
2629 | if (iLay==4) if (i%2==0) localY = -localY; |
2630 | sprintf(name, "ITSsddLay%iSensorPos%i",iLay, i); |
fa4639a3 |
2631 | |
2632 | if (i >= nDetectors/2) { |
2633 | TGeoTranslation *sensorPos = new TGeoTranslation(0,localY,localZ); |
2634 | sensorPos->SetName(name); |
2635 | virtualDet->AddNode(fSDDsensor, i, sensorPos); |
2636 | } |
2637 | else { |
531d6cdc |
2638 | TGeoRotation *rotSensor = new TGeoRotation("",0, 180, 180); |
fa4639a3 |
2639 | TGeoCombiTrans *sensorPos = new TGeoCombiTrans(0,localY, |
531d6cdc |
2640 | localZ, rotSensor); |
fa4639a3 |
2641 | sensorPos->SetName(name); |
2642 | virtualDet->AddNode(fSDDsensor, i, sensorPos); |
2643 | }; |
db486a6e |
2644 | } |
b7943f00 |
2645 | |
2646 | if(GetDebug(1)) virtualDet->CheckOverlaps(0.01); |
fa4639a3 |
2647 | virtualDet->SetVisibility(kFALSE); |
db486a6e |
2648 | return virtualDet; |
531d6cdc |
2649 | } |
c789ee28 |
2650 | |
b7943f00 |
2651 | |
db486a6e |
2652 | //________________________________________________________________________ |
b7943f00 |
2653 | Int_t AliITSv11GeometrySDD::ExportSensorGeometry(AliITSgeom *geom, Int_t iLaySDD, |
2654 | Int_t startMod) { |
2655 | // |
2656 | // export the geometry in a AliITSgeom object |
2657 | // |
2658 | |
2659 | if (! geom) { |
2660 | printf("error:Try to fill null (AliITSgeom *) object"); |
2661 | return kFALSE; |
2662 | }; |
2663 | if (! fMotherVol) { |
2664 | printf("error:Try to set sensor geometry while geometry is not defined\n"); |
2665 | return kFALSE; |
2666 | }; |
c789ee28 |
2667 | |
fa4639a3 |
2668 | const Float_t kDxyz[3] = {fgkWaferWidthSens/2., fgkWaferThickSens/2., |
2669 | fgkWaferLengthSens/2.}; |
b7943f00 |
2670 | if(!(geom->IsShapeDefined(kSDD))) |
2671 | geom->ReSetShape(kSDD, new AliITSgeomSDD256(3, kDxyz)); |
2672 | |
2673 | char layerName[30]; |
2674 | char ladderName[30]; |
2675 | char sensorName[30]; |
2676 | char senstivName[30]; |
2677 | const Int_t kNLay = 2; |
fa4639a3 |
2678 | const Int_t kNLadd[kNLay] = {fgkLay3Nladd, fgkLay4Nladd}; |
2679 | const Int_t kNDet[kNLay] = {fgkLay3Ndet, fgkLay4Ndet}; |
b7943f00 |
2680 | |
2681 | if (GetDebug(1)) |
2682 | printf("AliITSv11GeometrySDD::SetSensorGeometry(), nodes found :\n"); |
2683 | |
fa4639a3 |
2684 | Int_t firstSDDmod = startMod; |
b7943f00 |
2685 | for (Int_t iLay=0; iLay<kNLay; iLay++) { |
2686 | ///////////////////////////////////////// |
2687 | sprintf(layerName, "ITSsddLayer%i_1",iLay+3); |
2688 | TGeoNode *layNode = fMotherVol->GetNode(layerName); |
2689 | if (layNode) { |
2690 | if (GetDebug(1)) printf("%s\n",layNode->GetName()); |
2691 | TGeoVolume *layVolume = layNode->GetVolume(); |
2692 | TGeoHMatrix layMatrix(*layNode->GetMatrix()); |
2693 | |
2694 | for (Int_t iLadd=0; iLadd<kNLadd[iLay]; iLadd++) { |
2695 | ///////////////////////////////////////// |
2696 | sprintf(ladderName, "ITSsddLadd_%i", iLadd); |
2697 | TGeoNode *laddNode = layVolume->GetNode(ladderName); |
2698 | if (laddNode) { |
2699 | if (GetDebug(1)) printf("| %s\n",laddNode->GetName()); |
2700 | TGeoVolume *laddVolume = laddNode->GetVolume(); |
2701 | TGeoHMatrix laddMatrix(layMatrix); |
2702 | laddMatrix.Multiply(laddNode->GetMatrix()); |
2703 | |
2704 | for (Int_t iDet=0; iDet<kNDet[iLay]; iDet++) { |
2705 | ///////////////////////////////////////// |
2706 | sprintf(sensorName, "ITSsddSensor_%i",iDet); |
2707 | TGeoNode *detNode = laddVolume->GetNode(sensorName); |
2708 | if (detNode) { |
2709 | if (GetDebug(1)) printf("| | %s\n",detNode->GetName()); |
2710 | TGeoVolume *detVolume = detNode->GetVolume(); |
2711 | TGeoHMatrix detMatrix(laddMatrix); |
2712 | detMatrix.Multiply(detNode->GetMatrix()); |
2713 | |
2714 | TGeoNode *wafNode = detVolume->GetNode("ITSsddWafer_1"); |
2715 | if (wafNode) { |
2716 | TGeoVolume *wafVolume = wafNode->GetVolume(); |
2717 | TGeoHMatrix wafMatrix(detMatrix); |
2718 | detMatrix.Multiply(wafNode->GetMatrix()); |
2719 | //-------------------------------------------------------- |
2720 | sprintf(senstivName, "%s%s", fgSDDsensitiveVolName,"_1"); |
2721 | TGeoNode *sensitivNode = wafVolume->GetNode(senstivName); |
2722 | if (sensitivNode) { |
2723 | TGeoHMatrix sensMatrix(wafMatrix); |
2724 | sensMatrix.Multiply(sensitivNode->GetMatrix()); |
2725 | |
fa4639a3 |
2726 | // Sticking to the convention for local wafer coordinate |
2727 | // in AliITSgeom : |
2728 | if (iDet >= kNDet[iLay]/2) { |
2729 | // TGeoRotation rotY("",0,180,0); |
2730 | TGeoRotation rotY("",-180,-180,0); |
2731 | sensMatrix.Multiply(&rotY); |
2732 | }; |
2733 | // Creating the matrix in AliITSgeom for |
2734 | // this sensitive volume : |
b7943f00 |
2735 | Double_t *trans = sensMatrix.GetTranslation(); |
2736 | Double_t *r = sensMatrix.GetRotationMatrix(); |
2737 | Double_t rot[10] = {r[0],r[1],r[2], |
2738 | r[3],r[4],r[5], |
fa4639a3 |
2739 | r[6],r[7],r[8], 1.0}; |
2740 | //rot[9]!=0.0 => not a unity matrix |
023ae34b |
2741 | geom->CreateMatrix(startMod,iLay+iLaySDD,iLadd+1,iDet+1, |
b7943f00 |
2742 | kSDD,trans,rot); |
2743 | // iLadd+1, iDet+1 because ladd. and det. start at +1 |
2744 | // elsewhere |
2745 | startMod++; |
fa4639a3 |
2746 | |
b7943f00 |
2747 | } else |
2748 | printf("Error (ExportSensorGeometry) %s not found !\n", |
2749 | senstivName); |
2750 | } else |
2751 | printf("Error (ExportSensorGeometry) %s not found !\n", |
2752 | "ITSsddWafer_1"); |
2753 | } else |
2754 | printf("Error (ExportSensorGeometry) %s not found !\n", |
2755 | sensorName); |
2756 | }; |
2757 | } else |
2758 | printf("Error (ExportSensorGeometry) %s not found !\n", |
2759 | ladderName); |
2760 | }; |
2761 | } else |
2762 | printf("Error (ExportSensorGeometry) %s not found !\n", |
2763 | layerName); |
2764 | }; |
c789ee28 |
2765 | |
b7943f00 |
2766 | return (startMod-firstSDDmod); |
531d6cdc |
2767 | } |
c789ee28 |
2768 | |
2769 | |
fa4639a3 |
2770 | //________________________________________________________________________ |
b7943f00 |
2771 | Int_t AliITSv11GeometrySDD:: |
2772 | GetCurrentLayLaddDet(Int_t &lay, Int_t &ladd, Int_t&det) const { |
2773 | // |
2774 | // Function which gives the layer, ladder and det. |
2775 | // index of the current volume. To be used in |
2776 | // AliITS::StepManager() |
2777 | // |
c789ee28 |
2778 | |
b7943f00 |
2779 | if (gGeoManager->GetLevel()<3) return kFALSE; |
2780 | // Get the det index : |
2781 | TGeoNode *node = gGeoManager->GetMother(2); |
2782 | if (!node) return kFALSE; |
2783 | det = node->GetNumber()+1; |
c789ee28 |
2784 | |
b7943f00 |
2785 | // Get the ladder index : |
2786 | node = gGeoManager->GetMother(3); |
2787 | if (!node) return kFALSE; |
2788 | ladd = node->GetNumber()+1; |
2789 | |
2790 | // Get the layer index : |
2791 | if (node->GetNdaughters()==fgkLay3Ndet) |
fa4639a3 |
2792 | lay = 3; // this has to be equal to the iLaySDD argument given to ExportSensorGeometry() !!! |
2793 | else lay = 4; |
b7943f00 |
2794 | |
2795 | return kTRUE; |
531d6cdc |
2796 | } |