ATLAS Offline Software
Loading...
Searching...
No Matches
CrackDMConstruction.cxx
Go to the documentation of this file.
1/*
2 Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
3*/
4
5#include "CrackRegionGeoModel/CrackDMConstruction.h"
6
7#include "GeoModelKernel/GeoElement.h"
8#include "GeoModelKernel/GeoMaterial.h"
9#include "GeoModelKernel/GeoFullPhysVol.h"
10#include "GeoModelKernel/GeoPhysVol.h"
11#include "GeoModelKernel/GeoVPhysVol.h"
12#include "GeoModelKernel/GeoLogVol.h"
13#include "GeoModelKernel/GeoTransform.h"
14#include "GeoModelKernel/GeoAlignableTransform.h"
15#include "GeoModelKernel/GeoIdentifierTag.h"
16#include "GeoModelKernel/GeoNameTag.h"
17#include "GeoModelKernel/GeoSerialTransformer.h"
18#include "GeoModelKernel/GeoXF.h"
19#include "GeoModelKernel/GeoTube.h"
20#include "GeoModelKernel/GeoPcon.h"
21#include "GeoModelKernel/GeoTubs.h"
22#include "GeoModelKernel/GeoCons.h"
23#include "GeoModelKernel/GeoBox.h"
24#include "GeoModelKernel/GeoTrap.h"
25#include "GeoModelKernel/GeoTrd.h"
26#include "GeoModelKernel/GeoShape.h"
27#include "GeoModelKernel/GeoShapeUnion.h"
28#include "GeoModelKernel/GeoShapeShift.h"
29#include "GeoModelKernel/GeoShapeSubtraction.h"
30#include "GeoModelKernel/Units.h"
31
32#include "StoreGate/StoreGateSvc.h"
33#include "GeoModelInterfaces/StoredMaterialManager.h"
34#include "GeoModelInterfaces/IGeoModelSvc.h"
35#include "GeoModelUtilities/DecodeVersionKey.h"
36#include "RDBAccessSvc/IRDBAccessSvc.h"
37#include "RDBAccessSvc/IRDBRecordset.h"
38#include "RDBAccessSvc/IRDBRecord.h"
39
40// For units:
41#include "GaudiKernel/PhysicalConstants.h"
42#include "GaudiKernel/MsgStream.h"
43#include "GaudiKernel/Bootstrap.h"
44
45// For functions:
46#include "GeoGenericFunctions/Variable.h"
47
48using namespace GeoGenfun;
49using namespace GeoXF;
50
51static const unsigned int NCrates=16;
52static const double Alfa=360*Gaudi::Units::deg/NCrates;
53static const double Enda=1155;
54static const double Endb=1695.2;
55static const double Endc=2771.6;
56static const double DYa=1155*tan(Alfa/2);
57static const double DYb=1695.2*tan(Alfa/2);
58static const double DYc=2771.6*tan(Alfa/2);
59
60
61void
62createSectorEnvelopes2FromDB (GeoFullPhysVol* envelope,
63 StoredMaterialManager* materialManager,
64 std::map<std::string, unsigned int>& trdMap,
65 IRDBRecordset& BarrelDMTrds,
66 std::map<std::string, unsigned int>& trapMap,
67 IRDBRecordset& BarrelDMTraps,
68 std::map<std::string, unsigned int>& boxMap,
69 IRDBRecordset& BarrelDMBoxes,
70 GENFUNCTION& f,
71 GeoBox* Box)
72{
73 unsigned int recordIndex = trdMap["SecE2"];
74 double SecE2xhlen1 = BarrelDMTrds[recordIndex]->getDouble("XHLEN1");
75 double SecE2xhlen2 = BarrelDMTrds[recordIndex]->getDouble("XHLEN2");
76 double SecE2ztr = BarrelDMTrds[recordIndex]->getDouble("ZTR");
77
78 recordIndex = trdMap["Spb0"];
79 double Spb0xhlen1 = BarrelDMTrds[recordIndex]->getDouble("XHLEN1");
80 double Spb0xhlen2 = BarrelDMTrds[recordIndex]->getDouble("XHLEN2");
81 double Spb0yhlen1 = BarrelDMTrds[recordIndex]->getDouble("YHLEN1");
82 double Spb0yhlen2 = BarrelDMTrds[recordIndex]->getDouble("YHLEN2");
83 double Spb0zhlen = BarrelDMTrds[recordIndex]->getDouble("ZHLEN");
84 double Spb0ytr = BarrelDMTrds[recordIndex]->getDouble("YTR");
85 double Spb0xrot = BarrelDMTrds[recordIndex]->getDouble("XROT");
86
87 recordIndex = boxMap["Box"];
88 //double Boxhlen = BarrelDMBoxes[recordIndex]->getDouble("HLEN");
89 //double Boxhwdt = BarrelDMBoxes[recordIndex]->getDouble("HWDT");
90 //double Boxhhgt = BarrelDMBoxes[recordIndex]->getDouble("HHGT");
91 double Boxxtr = BarrelDMBoxes[recordIndex]->getDouble("XTR");
92 double Boxytr = BarrelDMBoxes[recordIndex]->getDouble("YTR");
93 double Boxztr = BarrelDMBoxes[recordIndex]->getDouble("ZTR");
94 //double Boxxrot = BarrelDMBoxes[recordIndex]->getDouble("XROT");
95
96 recordIndex = boxMap["SplBox"];
97 double SplBoxhlen = BarrelDMBoxes[recordIndex]->getDouble("HLEN");
98 double SplBoxhwdt = BarrelDMBoxes[recordIndex]->getDouble("HWDT");
99 double SplBoxhhgt = BarrelDMBoxes[recordIndex]->getDouble("HHGT");
100 double SplBoxytr = BarrelDMBoxes[recordIndex]->getDouble("YTR");
101 double SplBoxztr = BarrelDMBoxes[recordIndex]->getDouble("ZTR");
102 double SplBoxxrot = BarrelDMBoxes[recordIndex]->getDouble("XROT");
103
104 recordIndex = trapMap["Spb1"];
105 double Spb1zhlen = BarrelDMTraps[recordIndex]->getDouble("ZHLEN");
106 double Spb1theta = BarrelDMTraps[recordIndex]->getDouble("THETA");
107 double Spb1phi = BarrelDMTraps[recordIndex]->getDouble("PHI");
108 double Spb1yzn = BarrelDMTraps[recordIndex]->getDouble("YZN");
109 double Spb1xynzn = BarrelDMTraps[recordIndex]->getDouble("XYNZN");
110 double Spb1xypzn = BarrelDMTraps[recordIndex]->getDouble("XYPZN");
111 double Spb1angn = BarrelDMTraps[recordIndex]->getDouble("ANGN");
112 double Spb1yzp = BarrelDMTraps[recordIndex]->getDouble("YZP");
113 double Spb1xynzp = BarrelDMTraps[recordIndex]->getDouble("XYNZP");
114 double Spb1xypzp = BarrelDMTraps[recordIndex]->getDouble("XYPZP");
115 double Spb1angp = BarrelDMTraps[recordIndex]->getDouble("ANGP");
116 double Spb1xtr = BarrelDMTraps[recordIndex]->getDouble("XTR");
117 double Spb1ytr = BarrelDMTraps[recordIndex]->getDouble("YTR");
118 double Spb1ztr = BarrelDMTraps[recordIndex]->getDouble("ZTR");
119 double Spb1xrot = BarrelDMTraps[recordIndex]->getDouble("XROT");
120
121 recordIndex = trapMap["Spb3"];
122 double Spb3zhlen = BarrelDMTraps[recordIndex]->getDouble("ZHLEN");
123 double Spb3theta = BarrelDMTraps[recordIndex]->getDouble("THETA");
124 double Spb3phi = BarrelDMTraps[recordIndex]->getDouble("PHI");
125 double Spb3yzn = BarrelDMTraps[recordIndex]->getDouble("YZN");
126 double Spb3xynzn = BarrelDMTraps[recordIndex]->getDouble("XYNZN");
127 double Spb3xypzn = BarrelDMTraps[recordIndex]->getDouble("XYPZN");
128 double Spb3angn = BarrelDMTraps[recordIndex]->getDouble("ANGN");
129 double Spb3yzp = BarrelDMTraps[recordIndex]->getDouble("YZP");
130 double Spb3xynzp = BarrelDMTraps[recordIndex]->getDouble("XYNZP");
131 double Spb3xypzp = BarrelDMTraps[recordIndex]->getDouble("XYPZP");
132 double Spb3angp = BarrelDMTraps[recordIndex]->getDouble("ANGP");
133 double Spb3ztr = BarrelDMTraps[recordIndex]->getDouble("ZTR");
134
135 recordIndex = trapMap["Spb2"];
136 double Spb2zhlen = BarrelDMTraps[recordIndex]->getDouble("ZHLEN");
137 double Spb2theta = BarrelDMTraps[recordIndex]->getDouble("THETA");
138 double Spb2phi = BarrelDMTraps[recordIndex]->getDouble("PHI");
139 double Spb2yzn = BarrelDMTraps[recordIndex]->getDouble("YZN");
140 double Spb2xynzn = BarrelDMTraps[recordIndex]->getDouble("XYNZN");
141 double Spb2xypzn = BarrelDMTraps[recordIndex]->getDouble("XYPZN");
142 double Spb2angn = BarrelDMTraps[recordIndex]->getDouble("ANGN");
143 double Spb2yzp = BarrelDMTraps[recordIndex]->getDouble("YZP");
144 double Spb2xynzp = BarrelDMTraps[recordIndex]->getDouble("XYNZP");
145 double Spb2xypzp = BarrelDMTraps[recordIndex]->getDouble("XYPZP");
146 double Spb2angp = BarrelDMTraps[recordIndex]->getDouble("ANGP");
147 double Spb2ytr = BarrelDMTraps[recordIndex]->getDouble("YTR");
148 double Spb2ztr = BarrelDMTraps[recordIndex]->getDouble("ZTR");
149
150 const GeoMaterial* matLArServices17 = materialManager->getMaterial("LAr::LArServices17");// 0.035*gram/cm3
151 const GeoMaterial* matLArServices18 = materialManager->getMaterial("LAr::LArServices18");// 0.240*gram/cm3
152 const GeoMaterial* matLArServices19 = materialManager->getMaterial("LAr::LArServices19");// 0.469*gram/cm3
153 const GeoMaterial* matLArServices20 = materialManager->getMaterial("LAr::LArServices20");// 0.353*gram/cm3
154 const GeoMaterial *alu = materialManager->getMaterial("std::Aluminium"); //2.7 g/cm3
155 const GeoMaterial *air = materialManager->getMaterial("std::Air"); //0.001214 g/cm3
156
157 GeoTrf::Transform3D Cut3Boxe = GeoTrf::Translate3D(Boxxtr, Boxytr, Boxztr)*GeoTrf::RotateX3D(-20*Gaudi::Units::deg)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
158 GeoTrf::Transform3D Cut4Boxe = GeoTrf::Translate3D(Boxxtr, -Boxytr,Boxztr)*GeoTrf::RotateX3D(20*Gaudi::Units::deg)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
159
160 // build 5 instances of SectorEnvelopes1 with 3 different materials!
161 GeoTrd *Trdair2 = new GeoTrd(SecE2xhlen1, SecE2xhlen2, DYb, DYc, (Endc-Endb)/2);
162 const GeoShape & SectorEnvelope= ((*Trdair2).
163 subtract((*Box) <<GeoTrf::Transform3D(Cut3Boxe)).
164 subtract((*Box) <<GeoTrf::Transform3D(Cut4Boxe)));
165
166 const GeoShape & SectorEnvelopes= ((SectorEnvelope).
167 add(SectorEnvelope << GeoTrf::TranslateY3D(-(DYb+DYc)*cos(Alfa/2)*cos(Alfa/2))*GeoTrf::TranslateZ3D(-(DYb+DYc)*0.5*sin(Alfa))*GeoTrf::RotateX3D(Alfa)));
168
169 GeoLogVol *lvse2r = new GeoLogVol("LAr::DM::SectorEnvelopes2r",&SectorEnvelopes,matLArServices20);
170 GeoPhysVol *sectorenvelopes2r = new GeoPhysVol(lvse2r); // for right-handed splice boxes
171
172 GeoLogVol *lvse2l = new GeoLogVol("LAr::DM::SectorEnvelopes2l",&SectorEnvelopes,matLArServices20);
173 GeoPhysVol *sectorenvelopes2l = new GeoPhysVol(lvse2l); // for left-handed splice boxes
174
175 GeoLogVol *lvse2h = new GeoLogVol("LAr::DM::SectorEnvelopes2h",&SectorEnvelopes,matLArServices19);
176 GeoPhysVol *sectorenvelopes2h = new GeoPhysVol(lvse2h); // no splice boxes horizontal at 0 & 180 Gaudi::Units::deg.
177
178 GeoLogVol *lvse2vup = new GeoLogVol("LAr::DM::SectorEnvelopes2vup",&SectorEnvelopes,matLArServices17);
179 GeoPhysVol *sectorenvelopes2vup = new GeoPhysVol(lvse2vup); // no splice boxes vertical up at 90 Gaudi::Units::deg
180
181 GeoLogVol *lvse2vd = new GeoLogVol("LAr::DM::SectorEnvelopes2Vd",&SectorEnvelopes,matLArServices18);
182 GeoPhysVol *sectorenvelopes2vd = new GeoPhysVol(lvse2vd); // no splice boxes vertical down at 270 Gaudi::Units::deg
183
184 //---------- Build Splice boxes for InDet optical fibers--------
185
186 GeoTrap *GeoTrap1 = new GeoTrap(Spb1zhlen, Spb1theta, Spb1phi, Spb1yzn, Spb1xynzn, Spb1xypzn, Spb1angn, Spb1yzp, Spb1xynzp, Spb1xypzp, Spb1angp);
187 GeoBox *Box1 = new GeoBox(SplBoxhlen, SplBoxhwdt, SplBoxhhgt);
188 const GeoShape & SpliceBox = ((*GeoTrap1).
189 subtract(*Box1 << GeoTrf::TranslateZ3D(SplBoxztr)*GeoTrf::TranslateY3D(-SplBoxytr)*GeoTrf::RotateX3D(SplBoxxrot*Gaudi::Units::deg)));
190
191 GeoTransform *xtr = new GeoTransform (GeoTrf::TranslateZ3D(Spb1ztr)*GeoTrf::TranslateY3D(-Spb1ytr)*GeoTrf::TranslateX3D(Spb1xtr)*GeoTrf::RotateX3D(Spb1xrot*Gaudi::Units::deg));
192 sectorenvelopes2r->add(xtr);
193 GeoLogVol *lvspbr = new GeoLogVol("LAr::DM::SPliceBoxr",&SpliceBox,alu);
194 GeoPhysVol *spliceboxr = new GeoPhysVol(lvspbr);
195 sectorenvelopes2r->add(spliceboxr);
196
197 GeoTransform *xtl = new GeoTransform (GeoTrf::TranslateZ3D(Spb1ztr)*GeoTrf::TranslateY3D(-Spb1ytr)*GeoTrf::TranslateX3D(Spb1xtr)*GeoTrf::RotateY3D(-180*Gaudi::Units::deg)*GeoTrf::RotateX3D(-(Alfa/2)));
198 sectorenvelopes2l->add(xtl);
199 GeoLogVol *lvspbl = new GeoLogVol("LAr::DM::SpliceBoxl",&SpliceBox,alu);
200 GeoPhysVol *spliceboxl = new GeoPhysVol(lvspbl);
201 sectorenvelopes2l->add(spliceboxl);
202
204 GeoTrd *Trd1 = new GeoTrd(Spb0xhlen1, Spb0xhlen2, Spb0yhlen1, Spb0yhlen2, Spb0zhlen);
205 GeoTrap *GeoTrap2 = new GeoTrap(Spb2zhlen, Spb2theta, Spb2phi, Spb2yzn, Spb2xynzn, Spb2xypzn, Spb2angn, Spb2yzp, Spb2xynzp, Spb2xypzp, Spb2angp);
206 GeoTrap *GeoTrap3 = new GeoTrap(Spb3zhlen, Spb3theta, Spb3phi, Spb3yzn, Spb3xynzn, Spb3xypzn, Spb3angn, Spb3yzp, Spb3xynzp, Spb3xypzp, Spb3angp);
207
208 GeoTransform *xt1 = new GeoTransform (GeoTrf::TranslateY3D(-Spb0ytr)*GeoTrf::RotateX3D(Spb0xrot*Gaudi::Units::deg));
209 spliceboxr->add(xt1);
210 spliceboxl->add(xt1);
211 GeoLogVol *lt1 = new GeoLogVol("LAr::DM::TBox1",Trd1,air);
212 GeoPhysVol *tbox1 = new GeoPhysVol(lt1);
213 spliceboxr->add(tbox1);
214 spliceboxl->add(tbox1);
215
216 GeoTransform *xt2 = new GeoTransform (GeoTrf::TranslateZ3D(Spb2ztr)*GeoTrf::TranslateY3D(Spb2ytr));
217 spliceboxr->add(xt2);
218 spliceboxl->add(xt2);
219 GeoLogVol *lt2 = new GeoLogVol("LAr::DM::TBox2",GeoTrap2,air);
220 GeoPhysVol *tbox2 = new GeoPhysVol(lt2);
221 spliceboxr->add(tbox2);
222 spliceboxl->add(tbox2);
223
224 GeoTransform *xt3 = new GeoTransform (GeoTrf::TranslateZ3D(-Spb3ztr));
225 spliceboxr->add(xt3);
226 spliceboxl->add(xt3);
227 GeoLogVol *lt3 = new GeoLogVol("LAr::DM::TBox3",GeoTrap3,air);
228 GeoPhysVol *tbox3 = new GeoPhysVol(lt3);
229 spliceboxr->add(tbox3);
230 spliceboxl->add(tbox3);
231
232 //-------------- Place volumes in LAr Envelope -------------------
233
234 TRANSFUNCTION seA2r = Pow(GeoTrf::RotateZ3D(1.0),8*f-(3*Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
235 TRANSFUNCTION seA2l = Pow(GeoTrf::RotateZ3D(1.0),8*f+(5*Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
236 TRANSFUNCTION seC2r = Pow(GeoTrf::RotateZ3D(1.0),8*f-(3*Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(-SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
237 TRANSFUNCTION seC2l = Pow(GeoTrf::RotateZ3D(1.0),8*f+(5*Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(-SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
238 TRANSFUNCTION seA2Vup = Pow(GeoTrf::RotateZ3D(1.0),f+(9*Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
239 TRANSFUNCTION seA2Vd = Pow(GeoTrf::RotateZ3D(1.0),f-(7*Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
240 TRANSFUNCTION seA2H = Pow(GeoTrf::RotateZ3D(1.0),8*f+(Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
241 TRANSFUNCTION seC2Vup = Pow(GeoTrf::RotateZ3D(1.0),f+(9*Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(-SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
242 TRANSFUNCTION seC2Vd = Pow(GeoTrf::RotateZ3D(1.0),f-(7*Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(-SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
243 TRANSFUNCTION seC2H = Pow(GeoTrf::RotateZ3D(1.0),8*f+(Alfa/2))*GeoTrf::TranslateX3D((Endb+Endc)/2)*GeoTrf::TranslateZ3D(-SecE2ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
244
245 GeoSerialTransformer *setA2r = new GeoSerialTransformer(sectorenvelopes2r,&seA2r, 2);
246 GeoSerialTransformer *setA2l = new GeoSerialTransformer(sectorenvelopes2l,&seA2l, 2);
247 GeoSerialTransformer *setC2r = new GeoSerialTransformer(sectorenvelopes2r,&seC2r, 2);
248 GeoSerialTransformer *setC2l = new GeoSerialTransformer(sectorenvelopes2l,&seC2l, 2);
249 GeoSerialTransformer *setA2Vup = new GeoSerialTransformer(sectorenvelopes2vup,&seA2Vup, 1);
250 GeoSerialTransformer *setA2Vd = new GeoSerialTransformer(sectorenvelopes2vd,&seA2Vd, 1);
251 GeoSerialTransformer *setA2H = new GeoSerialTransformer(sectorenvelopes2h,&seA2H, 2);
252 GeoSerialTransformer *setC2Vup = new GeoSerialTransformer(sectorenvelopes2vup,&seC2Vup, 1);
253 GeoSerialTransformer *setC2Vd = new GeoSerialTransformer(sectorenvelopes2vd,&seC2Vd, 1);
254 GeoSerialTransformer *setC2H = new GeoSerialTransformer(sectorenvelopes2h,&seC2H, 2);
255
256 envelope->add(setA2r);
257 envelope->add(setA2l);
258 envelope->add(setC2r);
259 envelope->add(setC2l);
260 envelope->add(setA2Vup);
261 envelope->add(setA2Vd);
262 envelope->add(setA2H);
263 envelope->add(setC2Vup);
264 envelope->add(setC2Vd);
265 envelope->add(setC2H);
266}
267
268
269void
270createBridgeEnvelopesFromDB (GeoFullPhysVol* envelope,
271 std::map<std::string, unsigned int>& trapMap,
272 IRDBRecordset& BarrelDMTraps,
273 const GeoMaterial* matLArServices8,
274 GENFUNCTION& f)
275{
276 unsigned int recordIndex = trapMap["BridgeE"];
277 const IRDBRecord* r = BarrelDMTraps[recordIndex];
278 double BridgeEzhlen = r->getDouble("ZHLEN");
279 double BridgeEtheta = r->getDouble("THETA");
280 double BridgeEphi = r->getDouble("PHI");
281 double BridgeEyzn = r->getDouble("YZN");
282 double BridgeExynzn = r->getDouble("XYNZN");
283 double BridgeExypzn = r->getDouble("XYPZN");
284 double BridgeEangn = r->getDouble("ANGN");
285 double BridgeEyzp = r->getDouble("YZP");
286 double BridgeExynzp = r->getDouble("XYNZP");
287 double BridgeExypzp = r->getDouble("XYPZP");
288 double BridgeEangp = r->getDouble("ANGP");
289 double BridgeExtr = r->getDouble("XTR");
290 double BridgeEztr = r->getDouble("ZTR");
291
292 GeoTrap *Trapair = new GeoTrap(BridgeEzhlen, BridgeEtheta*Gaudi::Units::deg, BridgeEphi, BridgeEyzn, BridgeExynzn, BridgeExypzn, BridgeEangn, BridgeEyzp, BridgeExynzp, BridgeExypzp, BridgeEangp);
293 GeoLogVol *lvbre = new GeoLogVol("LAr::DM::BridgeEnvelopes",Trapair,matLArServices8);//In the end Density at least >= than SE1 because of Cryo Pipes
294 GeoPhysVol *bridgeenvelopes = new GeoPhysVol(lvbre);
295
296 TRANSFUNCTION breA = Pow(GeoTrf::RotateZ3D(1.0),f-(Alfa/2))*GeoTrf::TranslateX3D(BridgeExtr)*GeoTrf::TranslateZ3D(BridgeEztr)*GeoTrf::RotateZ3D(90*Gaudi::Units::deg)*GeoTrf::RotateY3D(90*Gaudi::Units::deg)*GeoTrf::RotateX3D(90*Gaudi::Units::deg);
297 TRANSFUNCTION breC = Pow(GeoTrf::RotateZ3D(1.0),f-(Alfa/2))*GeoTrf::TranslateX3D(BridgeExtr)*GeoTrf::TranslateZ3D(-BridgeEztr)*GeoTrf::RotateZ3D(-90*Gaudi::Units::deg)*GeoTrf::RotateY3D(-90*Gaudi::Units::deg)*GeoTrf::RotateX3D(-90*Gaudi::Units::deg);
298 GeoSerialTransformer *bretA = new GeoSerialTransformer(bridgeenvelopes,&breA, NCrates);
299 GeoSerialTransformer *bretC = new GeoSerialTransformer(bridgeenvelopes,&breC, NCrates);
300 envelope->add(bretA);
301 envelope->add(bretC);
302}
303
304
305void
306createBaseEnvelopesFromDB (GeoFullPhysVol* envelope,
307 std::map<std::string, unsigned int>& trdMap,
308 IRDBRecordset& BarrelDMTrds,
309 const GeoMaterial* matLArServices8,
310 GENFUNCTION& f)
311{
312 unsigned int recordIndex = trdMap["BaseE"];
313 const IRDBRecord* r = BarrelDMTrds[recordIndex];
314 double BaseExhlen1 = r->getDouble("XHLEN1");
315 double BaseExhlen2 = r->getDouble("XHLEN2");
316 double BaseEyhlen1 = r->getDouble("YHLEN1");
317 double BaseEyhlen2 = r->getDouble("YHLEN2");
318 double BaseEzhlen = r->getDouble("ZHLEN");//
319 double BaseExtr = r->getDouble("XTR");
320 double BaseEztr = r->getDouble("ZTR");
321
322 GeoTrd *Trd1air = new GeoTrd(BaseExhlen1, BaseExhlen2, BaseEyhlen1, BaseEyhlen2, BaseEzhlen);
323 GeoLogVol *lvbe = new GeoLogVol("LAr::DM::BaseEnvelopes",Trd1air,matLArServices8); //In the end Density at least >= than SE1 because of Cryo Pipes
324 GeoPhysVol *baseenvelopes = new GeoPhysVol(lvbe);
325
326 TRANSFUNCTION beA = Pow(GeoTrf::RotateZ3D(1.0),f-(Alfa/2))*GeoTrf::TranslateX3D(BaseExtr)*GeoTrf::TranslateZ3D(BaseEztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
327 TRANSFUNCTION beC = Pow(GeoTrf::RotateZ3D(1.0),f+(Alfa/2))*GeoTrf::TranslateX3D(BaseExtr)*GeoTrf::TranslateZ3D(-BaseEztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
328 GeoSerialTransformer *betA = new GeoSerialTransformer(baseenvelopes,&beA, NCrates);
329 GeoSerialTransformer *betC = new GeoSerialTransformer(baseenvelopes,&beC, NCrates);
330 envelope->add(betA);
331 envelope->add(betC);
332}
333
334void createFromDB (GeoFullPhysVol* envelope,
335 IRDBAccessSvc* rdbAccess,
336 IGeoModelSvc* geoModel,
337 StoredMaterialManager* materialManager)
338{
339 // Use Geometry Database
340 DecodeVersionKey keyLAr(geoModel,"LAr");
341 IRDBRecordset_ptr BarrelDMTraps = rdbAccess->getRecordsetPtr("BarrelDMTraps",keyLAr.tag(),keyLAr.node());
342 IRDBRecordset_ptr BarrelDMTrds = rdbAccess->getRecordsetPtr("BarrelDMTrds",keyLAr.tag(),keyLAr.node());
343 IRDBRecordset_ptr BarrelDMTubes = rdbAccess->getRecordsetPtr("BarrelDMTubes",keyLAr.tag(),keyLAr.node());
344 IRDBRecordset_ptr BarrelDMBoxes = rdbAccess->getRecordsetPtr("BarrelDMBoxes",keyLAr.tag(),keyLAr.node());
345
346 std::map<std::string, unsigned int> tubeMap;
347 for (unsigned int i=0; i<BarrelDMTubes->size(); i++)
348 {
349 const std::string& key = (*BarrelDMTubes)[i]->getString("TUBENAME");
350 tubeMap[key] = i;
351 }
352 std::map<std::string, unsigned int> boxMap;
353 for (unsigned int j=0; j<BarrelDMBoxes->size(); j++)
354 {
355 const std::string& key = (*BarrelDMBoxes)[j]->getString("BOXNAME");
356 boxMap[key] = j;
357 }
358 std::map<std::string, unsigned int> trdMap;
359 for (unsigned int k=0; k<BarrelDMTrds->size(); k++)
360 {
361 const std::string& key = (*BarrelDMTrds)[k]->getString("TRDNAME");
362 trdMap[key] = k;
363 }
364 std::map<std::string, unsigned int> trapMap;
365 for (unsigned int l=0; l<BarrelDMTraps->size(); l++)
366 {
367 const std::string& key = (*BarrelDMTraps)[l]->getString("TRAPNAME");
368 trapMap[key] = l;
369 }
370
371 unsigned int recordIndex;
372
373 // Get materials
374 const GeoMaterial *alu = materialManager->getMaterial("std::Aluminium"); //2.7 g/cm3
375 const GeoMaterial* matBoardsEnvelope = materialManager->getMaterial("LAr::BoardsEnvelope");// 0.932*gram/cm3);
376 const GeoMaterial* matLArServices1 = materialManager->getMaterial("LAr::LArServices1");// 1.020*gram/cm3
377 const GeoMaterial* matLArServices2 = materialManager->getMaterial("LAr::LArServices2");// 0.955*gram/cm3
378 const GeoMaterial* matLArServices3 = materialManager->getMaterial("LAr::LArServices3");// 1.005*gram/cm3
379 const GeoMaterial* matLArServices4 = materialManager->getMaterial("LAr::LArServices4");// 0.460*gram/cm3
380 const GeoMaterial* matLArServices5 = materialManager->getMaterial("LAr::LArServices5");// 0.480*gram/cm3
381 const GeoMaterial* matLArServices6 = materialManager->getMaterial("LAr::LArServices6");// 1.000*gram/cm3
382 const GeoMaterial* matLArServices7 = materialManager->getMaterial("LAr::LArServices7");// 0.935*gram/cm3
383 const GeoMaterial* matLArServices8 = materialManager->getMaterial("LAr::LArServices8");// 1.070*gram/cm3
384 const GeoMaterial* matLArServices9 = materialManager->getMaterial("LAr::LArServices9");// 1.020*gram/cm3
385 const GeoMaterial* matLArServices10 = materialManager->getMaterial("LAr::LArServices10");// 0.995*gram/cm3
386 const GeoMaterial* matLArServices11 = materialManager->getMaterial("LAr::LArServices11");// 0.835*gram/cm3
387 const GeoMaterial* matLArServices12 = materialManager->getMaterial("LAr::LArServices12");// 0.640*gram/cm3
388 const GeoMaterial* matLArServices13 = materialManager->getMaterial("LAr::LArServices13");// 0.690*gram/cm3
389 const GeoMaterial* matLArServices14 = materialManager->getMaterial("LAr::LArServices14");// 0.825*gram/cm3
390 const GeoMaterial* matLArServices15 = materialManager->getMaterial("LAr::LArServices15");// 0.875*gram/cm3
391 const GeoMaterial* matLArServices16 = materialManager->getMaterial("LAr::LArServices16");// 1.035*gram/cm3
392
393 const double inv_Endab = 1. / (Endb - Enda);
394 Variable i;
395 GENFUNCTION f = Alfa*i;
396
398
399 recordIndex = tubeMap["Ped2"];
400 double ped2zhlen = (*BarrelDMTubes)[recordIndex]->getDouble("ZHLEN");
401 double ped2minr = (*BarrelDMTubes)[recordIndex]->getDouble("MINR");
402 double ped2maxr = (*BarrelDMTubes)[recordIndex]->getDouble("MAXR");
403 double ped2ytr = (*BarrelDMTubes)[recordIndex]->getDouble("YTR");
404 recordIndex = tubeMap["Ped3"];
405 double ped3zhlen = (*BarrelDMTubes)[recordIndex]->getDouble("ZHLEN");
406 double ped3minr = (*BarrelDMTubes)[recordIndex]->getDouble("MINR");
407 double ped3maxr = (*BarrelDMTubes)[recordIndex]->getDouble("MAXR");
408 double ped3xtr = (*BarrelDMTubes)[recordIndex]->getDouble("XTR");
409
410 recordIndex = boxMap["Pedest"];
411 double pedesthlen = (*BarrelDMBoxes)[recordIndex]->getDouble("HLEN");
412 double pedesthwdt = (*BarrelDMBoxes)[recordIndex]->getDouble("HWDT");
413 double pedesthhgt = (*BarrelDMBoxes)[recordIndex]->getDouble("HHGT");
414 double pedestxtr = (*BarrelDMBoxes)[recordIndex]->getDouble("XTR");
415 double pedestztr = (*BarrelDMBoxes)[recordIndex]->getDouble("ZTR");
416 recordIndex = boxMap["Ped1"];
417 double ped1hlen = (*BarrelDMBoxes)[recordIndex]->getDouble("HLEN");
418 double ped1hwdt = (*BarrelDMBoxes)[recordIndex]->getDouble("HWDT");
419 double ped1hhgt = (*BarrelDMBoxes)[recordIndex]->getDouble("HHGT");
420 recordIndex = boxMap["Crate1"];
421 double crate1hlen = (*BarrelDMBoxes)[recordIndex]->getDouble("HLEN");
422 double crate1hwdt = (*BarrelDMBoxes)[recordIndex]->getDouble("HWDT");
423 double crate1hhgt = (*BarrelDMBoxes)[recordIndex]->getDouble("HHGT");
424 double crate1xtr = (*BarrelDMBoxes)[recordIndex]->getDouble("XTR");
425 double crate1ztr = (*BarrelDMBoxes)[recordIndex]->getDouble("ZTR");
426 recordIndex = boxMap["Crate2"];
427 double crate2hlen = (*BarrelDMBoxes)[recordIndex]->getDouble("HLEN");
428 double crate2hwdt = (*BarrelDMBoxes)[recordIndex]->getDouble("HWDT");
429 double crate2hhgt = (*BarrelDMBoxes)[recordIndex]->getDouble("HHGT");
430 recordIndex = boxMap["Crate3"];
431 double crate3hlen = (*BarrelDMBoxes)[recordIndex]->getDouble("HLEN");
432 double crate3hwdt = (*BarrelDMBoxes)[recordIndex]->getDouble("HWDT");
433 double crate3hhgt = (*BarrelDMBoxes)[recordIndex]->getDouble("HHGT");
434 double crate3xtr = (*BarrelDMBoxes)[recordIndex]->getDouble("XTR");
435 recordIndex = boxMap["BoardE"];
436 double BoardEhlen = (*BarrelDMBoxes)[recordIndex]->getDouble("HLEN");
437 double BoardEhwdt = (*BarrelDMBoxes)[recordIndex]->getDouble("HWDT");
438 double BoardEhhgt = (*BarrelDMBoxes)[recordIndex]->getDouble("HHGT");
439 double BoardExtr = (*BarrelDMBoxes)[recordIndex]->getDouble("XTR");
440 double BoardEytr = (*BarrelDMBoxes)[recordIndex]->getDouble("YTR");
441 double BoardEztr = (*BarrelDMBoxes)[recordIndex]->getDouble("ZTR");
442 recordIndex = boxMap["Box"];
443 double Boxhlen = (*BarrelDMBoxes)[recordIndex]->getDouble("HLEN");
444 double Boxhwdt = (*BarrelDMBoxes)[recordIndex]->getDouble("HWDT");
445 double Boxhhgt = (*BarrelDMBoxes)[recordIndex]->getDouble("HHGT");
446 double Boxxtr = (*BarrelDMBoxes)[recordIndex]->getDouble("XTR");
447 double Boxytr = (*BarrelDMBoxes)[recordIndex]->getDouble("YTR");
448 double Boxxrot = (*BarrelDMBoxes)[recordIndex]->getDouble("XROT");
449
450 // ----- build pedestal -----
451 GeoBox *Pedestal = new GeoBox(pedesthlen, pedesthwdt, pedesthhgt);
452 GeoBox *Ped1 = new GeoBox(ped1hlen, ped1hwdt, ped1hhgt);
453 GeoTube *Ped2 = new GeoTube(ped2minr, ped2maxr, ped2zhlen);
454 GeoTube *Ped3 = new GeoTube(ped3minr,ped3maxr , ped3zhlen);
455 const GeoShape & CratePed=((*Pedestal).subtract(*Ped1).
456 subtract((*Ped2) <<GeoTrf::TranslateY3D(-ped2ytr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg)).
457 subtract((*Ped3) <<GeoTrf::TranslateX3D(-ped3xtr)).
458 subtract((*Ped2) <<GeoTrf::TranslateY3D(ped2ytr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg)));
459
460 GeoLogVol *lvped = new GeoLogVol("LAr::DM::Ped",&CratePed,alu);
461 GeoPhysVol *pedestal = new GeoPhysVol(lvped);
462
463 // ----- build crates -----
464 GeoBox *Crate1 = new GeoBox(crate1hlen, crate1hwdt, crate1hhgt);
465 GeoBox *Crate2 = new GeoBox(crate2hlen, crate2hwdt, crate2hhgt);
466 GeoBox *Crate3 = new GeoBox(crate3hlen, crate3hwdt, crate3hhgt);
467 const GeoShape & FEBCrate=(*Crate1).subtract(*Crate2).add((*Crate3) <<GeoTrf::TranslateX3D(-crate3xtr));
468
469 GeoLogVol *lvcrate = new GeoLogVol("LAr::DM::Crate",&FEBCrate,alu);
470 GeoPhysVol *crate = new GeoPhysVol(lvcrate);
471
472 // ----- build boardenvelopes -----
473 GeoBox *BoardEnvelope = new GeoBox(BoardEhlen, BoardEhwdt, BoardEhhgt);
474 GeoLogVol *lvbenv = new GeoLogVol("LAr::DM::FEBoard",BoardEnvelope,matBoardsEnvelope);
475 GeoPhysVol *boardenvelope = new GeoPhysVol(lvbenv);
476
477 //-------------- Place volumes in envelope ----------------------------
478
479 //Crates
480 TRANSFUNCTION crA = Pow(GeoTrf::RotateZ3D(1.0),f)*GeoTrf::TranslateX3D(crate1xtr)*GeoTrf::TranslateZ3D(crate1ztr);
481 TRANSFUNCTION crC = Pow(GeoTrf::RotateZ3D(1.0),f)*GeoTrf::TranslateX3D(crate1xtr)*GeoTrf::TranslateZ3D(-crate1ztr);
482 GeoSerialTransformer *crtA = new GeoSerialTransformer(crate,&crA, NCrates);
483 GeoSerialTransformer *crtC = new GeoSerialTransformer(crate,&crC, NCrates);
484 envelope->add(crtA);
485 envelope->add(crtC);
486
487 //Pedestals
488 TRANSFUNCTION pedA = Pow(GeoTrf::RotateZ3D(1.0),f)*GeoTrf::TranslateX3D(pedestxtr)*GeoTrf::TranslateZ3D(pedestztr);
489 TRANSFUNCTION pedC = Pow(GeoTrf::RotateZ3D(1.0),f)*GeoTrf::TranslateX3D(pedestxtr)*GeoTrf::TranslateZ3D(-pedestztr);
490 GeoSerialTransformer *pedtA = new GeoSerialTransformer(pedestal,&pedA, NCrates);
491 GeoSerialTransformer *pedtC = new GeoSerialTransformer(pedestal,&pedC, NCrates);
492 envelope->add(pedtA);
493 envelope->add(pedtC);
494
495 //FEBoards
496 TRANSFUNCTION feb1A = Pow(GeoTrf::RotateZ3D(1.0),f)*GeoTrf::TranslateY3D(BoardEytr)*GeoTrf::TranslateX3D(BoardExtr)*GeoTrf::TranslateZ3D(BoardEztr);
497 TRANSFUNCTION feb2A = Pow(GeoTrf::RotateZ3D(1.0),f)*GeoTrf::TranslateY3D(-BoardEytr)*GeoTrf::TranslateX3D(BoardExtr)*GeoTrf::TranslateZ3D(BoardEztr);
498 TRANSFUNCTION feb1C = Pow(GeoTrf::RotateZ3D(1.0),f)*GeoTrf::TranslateY3D(BoardEytr)*GeoTrf::TranslateX3D(BoardExtr)*GeoTrf::TranslateZ3D(-BoardEztr);
499 TRANSFUNCTION feb2C = Pow(GeoTrf::RotateZ3D(1.0),f)*GeoTrf::TranslateY3D(-BoardEytr)*GeoTrf::TranslateX3D(BoardExtr)*GeoTrf::TranslateZ3D(-BoardEztr);
500 GeoSerialTransformer *febt1A = new GeoSerialTransformer(boardenvelope,&feb1A, NCrates);
501 GeoSerialTransformer *febt1C = new GeoSerialTransformer(boardenvelope,&feb1C, NCrates);
502 GeoSerialTransformer *febt2A = new GeoSerialTransformer(boardenvelope,&feb2A, NCrates);
503 GeoSerialTransformer *febt2C = new GeoSerialTransformer(boardenvelope,&feb2C, NCrates);
504 envelope->add(febt1A);
505 envelope->add(febt1C);
506 envelope->add(febt2A);
507 envelope->add(febt2C);
508
510 recordIndex = trdMap["SecP"];
511 double SecPxhlen1 = (*BarrelDMTrds)[recordIndex]->getDouble("XHLEN1");
512 double SecPxhlen2 = (*BarrelDMTrds)[recordIndex]->getDouble("XHLEN2");
513 double SecPyhlen1 = (*BarrelDMTrds)[recordIndex]->getDouble("YHLEN1");//
514 double SecPyhlen2 = (*BarrelDMTrds)[recordIndex]->getDouble("YHLEN2");//
515 double SecPzhlen = (*BarrelDMTrds)[recordIndex]->getDouble("ZHLEN");//
516 double SecPxtr = (*BarrelDMTrds)[recordIndex]->getDouble("XTR");
517 double SecPztr = (*BarrelDMTrds)[recordIndex]->getDouble("ZTR");
518 recordIndex = trdMap["BaseP"];
519 double BasePxhlen1 = (*BarrelDMTrds)[recordIndex]->getDouble("XHLEN1");
520 double BasePxhlen2 = (*BarrelDMTrds)[recordIndex]->getDouble("XHLEN2");
521 double BasePyhlen1 = (*BarrelDMTrds)[recordIndex]->getDouble("YHLEN1");
522 double BasePyhlen2 = (*BarrelDMTrds)[recordIndex]->getDouble("YHLEN2");
523 double BasePzhlen = (*BarrelDMTrds)[recordIndex]->getDouble("ZHLEN");
524 double BasePxtr = (*BarrelDMTrds)[recordIndex]->getDouble("XTR");
525 double BasePztr = (*BarrelDMTrds)[recordIndex]->getDouble("ZTR");
526 recordIndex = trdMap["SecE1"];
527 double SecE1xhlen1 = (*BarrelDMTrds)[recordIndex]->getDouble("XHLEN1");
528 double SecE1xhlen2 = (*BarrelDMTrds)[recordIndex]->getDouble("XHLEN2");
529 double SecE1ztr = (*BarrelDMTrds)[recordIndex]->getDouble("ZTR");
530 recordIndex = trapMap["BridgeP"];
531 double BridgePzhlen = (*BarrelDMTraps)[recordIndex]->getDouble("ZHLEN");
532 double BridgePtheta = (*BarrelDMTraps)[recordIndex]->getDouble("THETA");
533 double BridgePphi = (*BarrelDMTraps)[recordIndex]->getDouble("PHI");
534 double BridgePyzn = (*BarrelDMTraps)[recordIndex]->getDouble("YZN");
535 double BridgePxynzn = (*BarrelDMTraps)[recordIndex]->getDouble("XYNZN");
536 double BridgePxypzn = (*BarrelDMTraps)[recordIndex]->getDouble("XYPZN");
537 double BridgePangn = (*BarrelDMTraps)[recordIndex]->getDouble("ANGN");
538 double BridgePyzp = (*BarrelDMTraps)[recordIndex]->getDouble("YZP");
539 double BridgePxynzp = (*BarrelDMTraps)[recordIndex]->getDouble("XYNZP");
540 double BridgePxypzp = (*BarrelDMTraps)[recordIndex]->getDouble("XYPZP");
541 double BridgePangp = (*BarrelDMTraps)[recordIndex]->getDouble("ANGP");
542 double BridgePxtr = (*BarrelDMTraps)[recordIndex]->getDouble("XTR");
543 double BridgePztr = (*BarrelDMTraps)[recordIndex]->getDouble("ZTR");
544
545 // transforms
546 GeoBox *Box = new GeoBox(Boxhlen, Boxhwdt, Boxhhgt);
547
548 GeoTrf::Transform3D Cut3Boxp = GeoTrf::Translate3D(Boxxtr, Boxytr, Boxxrot)*GeoTrf::RotateX3D(-20*Gaudi::Units::deg)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
549 GeoTrf::Transform3D Cut4Boxp = GeoTrf::Translate3D(Boxxtr, -Boxytr,Boxxrot)*GeoTrf::RotateX3D(20*Gaudi::Units::deg)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
550
551 // ----- build sector envelopes -----
552 // build 16 instances of SectorEnvelopes1 each with its own material!
553 GeoTrd *Trdair1 = new GeoTrd(SecE1xhlen1, SecE1xhlen2, DYa, DYb, (Endb-Enda)/2);
554
555 GeoLogVol *lvse1g1 = new GeoLogVol("LAr::DM::SectorEnvelopes1g1",Trdair1,matLArServices1);
556 GeoPhysVol *sectorenvelopes1g1 = new GeoPhysVol(lvse1g1);
557
558 GeoLogVol *lvse1g2 = new GeoLogVol("LAr::DM::SectorEnvelopes1g3",Trdair1,matLArServices2);
559 GeoPhysVol *sectorenvelopes1g2 = new GeoPhysVol(lvse1g2);
560
561 GeoLogVol *lvse1g3 = new GeoLogVol("LAr::DM::SectorEnvelopes1g3",Trdair1,matLArServices3);
562 GeoPhysVol *sectorenvelopes1g3 = new GeoPhysVol(lvse1g3);
563
564 GeoLogVol *lvse1g4 = new GeoLogVol("LAr::DM::SectorEnvelopes1g4",Trdair1,matLArServices4);
565 GeoPhysVol *sectorenvelopes1g4 = new GeoPhysVol(lvse1g4);
566
567 GeoLogVol *lvse1g5 = new GeoLogVol("LAr::DM::SectorEnvelopes1g5",Trdair1,matLArServices5);
568 GeoPhysVol *sectorenvelopes1g5 = new GeoPhysVol(lvse1g5);
569
570 GeoLogVol *lvse1g6 = new GeoLogVol("LAr::DM::SectorEnvelopes1g6",Trdair1,matLArServices6);
571 GeoPhysVol *sectorenvelopes1g6 = new GeoPhysVol(lvse1g6);
572
573 GeoLogVol *lvse1g7 = new GeoLogVol("LAr::DM::SectorEnvelopes1g7",Trdair1,matLArServices7);
574 GeoPhysVol *sectorenvelopes1g7 = new GeoPhysVol(lvse1g7);
575
576 GeoLogVol *lvse1g8 = new GeoLogVol("LAr::DM::SectorEnvelopes1g8",Trdair1,matLArServices8);
577 GeoPhysVol *sectorenvelopes1g8 = new GeoPhysVol(lvse1g8);
578
579 GeoLogVol *lvse1g9 = new GeoLogVol("LAr::DM::SectorEnvelopes1g9",Trdair1,matLArServices9);
580 GeoPhysVol *sectorenvelopes1g9 = new GeoPhysVol(lvse1g9);
581
582 GeoLogVol *lvse1g10 = new GeoLogVol("LAr::DM::SectorEnvelopes1g10",Trdair1,matLArServices10);
583 GeoPhysVol *sectorenvelopes1g10 = new GeoPhysVol(lvse1g10);
584
585 GeoLogVol *lvse1g11 = new GeoLogVol("LAr::DM::SectorEnvelopes1g11",Trdair1,matLArServices11);
586 GeoPhysVol *sectorenvelopes1g11 = new GeoPhysVol(lvse1g11);
587
588 GeoLogVol *lvse1g12 = new GeoLogVol("LAr::DM::SectorEnvelopes1g12",Trdair1,matLArServices12);
589 GeoPhysVol *sectorenvelopes1g12 = new GeoPhysVol(lvse1g12);
590
591 GeoLogVol *lvse1g13 = new GeoLogVol("LAr::DM::SectorEnvelopes1g13",Trdair1,matLArServices13);
592 GeoPhysVol *sectorenvelopes1g13 = new GeoPhysVol(lvse1g13);
593
594 GeoLogVol *lvse1g14 = new GeoLogVol("LAr::DM::SectorEnvelopes1g14",Trdair1,matLArServices14);
595 GeoPhysVol *sectorenvelopes1g14 = new GeoPhysVol(lvse1g14);
596
597 GeoLogVol *lvse1g15 = new GeoLogVol("LAr::DM::SectorEnvelopes1g15",Trdair1,matLArServices15);
598 GeoPhysVol *sectorenvelopes1g15 = new GeoPhysVol(lvse1g15);
599
600 GeoLogVol *lvse1g16 = new GeoLogVol("LAr::DM::SectorEnvelopes1g16",Trdair1,matLArServices16);
601 GeoPhysVol *sectorenvelopes1g16 = new GeoPhysVol(lvse1g16);
602
603
604 std::vector<GeoPhysVol*> se1List;
605 se1List.push_back(sectorenvelopes1g1);
606 se1List.push_back(sectorenvelopes1g2);
607 se1List.push_back(sectorenvelopes1g3);
608 se1List.push_back(sectorenvelopes1g4);
609 se1List.push_back(sectorenvelopes1g5);
610 se1List.push_back(sectorenvelopes1g6);
611 se1List.push_back(sectorenvelopes1g7);
612 se1List.push_back(sectorenvelopes1g8);
613 se1List.push_back(sectorenvelopes1g9);
614 se1List.push_back(sectorenvelopes1g10);
615 se1List.push_back(sectorenvelopes1g11);
616 se1List.push_back(sectorenvelopes1g12);
617 se1List.push_back(sectorenvelopes1g13);
618 se1List.push_back(sectorenvelopes1g14);
619 se1List.push_back(sectorenvelopes1g15);
620 se1List.push_back(sectorenvelopes1g16);
621
622 // dedicated volumes to add in SectorEnvelopes1
623 // (a) volume in Trd shape to add slice of extra material in a given eta bin
624 IRDBRecordset_ptr BarrelDMRing = rdbAccess->getRecordsetPtr("LArDMEnv1Ring",keyLAr.tag(),keyLAr.node());
625 for (unsigned int i=0;i<BarrelDMRing->size();i++) {
626 double etaMin=(*BarrelDMRing)[i]->getDouble("ETAMIN");
627 double etaMax=(*BarrelDMRing)[i]->getDouble("ETAMAX");
628 double thicknessExtra=(*BarrelDMRing)[i]->getDouble("THICKNESS");
629 std::string ringName = "LAr::DM::SectorEnvelopes1::"+(*BarrelDMRing)[i]->getString("RINGNAME");
630
631 double radiusMin=SecE1ztr/sinh(etaMax);
632 double radiusMax=SecE1ztr/sinh(etaMin);
633 double dy1 = DYa + (DYb-DYa)*(radiusMin-Enda)*inv_Endab;
634 double dy2 = DYa + (DYb-DYa)*(radiusMax-Enda)*inv_Endab;
635 double zpos=0.5*(radiusMax+radiusMin - (Endb+Enda));
636 const GeoMaterial *matExtraTdr = materialManager->getMaterial((*BarrelDMRing)[i]->getString("MATERIAL"));
637
638 GeoTrd *extraMatTdr = new GeoTrd(thicknessExtra/2., thicknessExtra/2., dy1, dy2, (radiusMax-radiusMin)/2);
639 GeoLogVol *extraMatLog = new GeoLogVol(ringName,extraMatTdr,matExtraTdr);
640 GeoPhysVol *extraMatPhys = new GeoPhysVol(extraMatLog);
641 for (unsigned int isect=0;isect<se1List.size();isect++) {
642 se1List[isect]->add(new GeoTransform(GeoTrf::TranslateZ3D(zpos)));
643 se1List[isect]->add(extraMatPhys);
644 }
645 }
646
647 // (b) extra material at fixed phi locations in the PPF1 area
648 IRDBRecordset_ptr BarrelDMPhiBox = rdbAccess->getRecordsetPtr("LArDMEnv1PhiBox",keyLAr.tag(),keyLAr.node());
649 for (unsigned int i=0;i<BarrelDMPhiBox->size();i++) {
650 double eta=(*BarrelDMPhiBox)[i]->getDouble("ETA");
651 double phi0=(*BarrelDMPhiBox)[i]->getDouble("PHI0");
652 double deltaR=(*BarrelDMPhiBox)[i]->getDouble("DELTAR");
653 double deltaRphi=(*BarrelDMPhiBox)[i]->getDouble("DELTARPHI");
654 double thickness=(*BarrelDMPhiBox)[i]->getDouble("THICKNESS");
655 int nphi=(*BarrelDMPhiBox)[i]->getInt("NPHI");
656 int noHorizontal = (*BarrelDMPhiBox)[i]->getInt("NOHORIZ");
657
658 const GeoMaterial* matExtraPPF1 = materialManager->getMaterial((*BarrelDMPhiBox)[i]->getString("MATERIAL"));
659 std::string boxName = "LAr::DM::SectorEnvelopes1::"+(*BarrelDMPhiBox)[i]->getString("BOXNAME");
660
661 GeoBox *ppf1Box = new GeoBox(thickness/2.,deltaRphi/2.,deltaR/2.);
662 GeoLogVol *ppf1Log = new GeoLogVol(boxName,ppf1Box,matExtraPPF1);
663 GeoPhysVol *ppf1Phys = new GeoPhysVol(ppf1Log);
664
665 int nPerEnv1 = nphi/NCrates;
666 double dphi=2.*M_PI/((float)(nphi));
667 double radius=SecE1ztr/sinh(eta);
668 for (int iphi=0;iphi<nPerEnv1;iphi++) {
669 double xpos=0.;
670 double ypos=radius*sin(phi0+((float)(iphi))*dphi);
671 double zpos=radius*cos(phi0+((float)(iphi))*dphi) - (Endb+Enda)/2.;
672 for (unsigned int isect=0;isect<se1List.size();isect++) {
673 // no PPF1 box around phi=0 and phi=pi
674 if (noHorizontal>0 && ((isect==7 && iphi==1) || (isect==8 && iphi==0) || (isect==15 && iphi==1) || (isect==0 && iphi==0) ) ) continue;
675 se1List[isect]->add(new GeoTransform(GeoTrf::Translate3D(xpos,ypos,zpos)));
676 se1List[isect]->add(ppf1Phys);
677 }
678 }
679 }
680
681 // ----- build base plates -----
682 GeoTrd *Trd1alu = new GeoTrd(BasePxhlen1, BasePxhlen2, BasePyhlen1, BasePyhlen2, BasePzhlen);
683 GeoLogVol *lvbp = new GeoLogVol("LAr::DM::BasePlates",Trd1alu,alu);
684 GeoPhysVol *baseplates = new GeoPhysVol(lvbp);
685
686 // ----- build bridge plates -----
687 GeoTrap *Trapalu = new GeoTrap(BridgePzhlen, BridgePtheta*Gaudi::Units::deg, BridgePphi, BridgePyzn, BridgePxynzn, BridgePxypzn, BridgePangn, BridgePyzp, BridgePxynzp, BridgePxypzp, BridgePangp);
688 GeoLogVol *lvbrp = new GeoLogVol("LAr::DM::BridgePlates",Trapalu,alu);
689 GeoPhysVol *bridgeplates = new GeoPhysVol(lvbrp);
690
691
692 // ----- build sector plates -----
693 GeoTrd *Trd2alu = new GeoTrd(SecPxhlen1, SecPxhlen2, SecPyhlen1, SecPyhlen2, SecPzhlen );
694 const GeoShape & SectorPlates= ((*Trd2alu).
695 subtract((*Box) <<GeoTrf::Transform3D(Cut3Boxp)).
696 subtract((*Box) <<GeoTrf::Transform3D(Cut4Boxp)));
697 GeoLogVol *lvsp = new GeoLogVol("LAr::DM::SectorPlates",&SectorPlates,alu);
698 GeoPhysVol *sectorplates = new GeoPhysVol(lvsp);
699
700
701 //-------------- Place volumes in LAr Envelope -------------------
702
703 //sectorPlates
704 TRANSFUNCTION spA = Pow(GeoTrf::RotateZ3D(1.0),f-(Alfa/2))*GeoTrf::TranslateX3D(SecPxtr)*GeoTrf::TranslateZ3D(SecPztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
705 TRANSFUNCTION spC = Pow(GeoTrf::RotateZ3D(1.0),f+(Alfa/2))*GeoTrf::TranslateX3D(SecPxtr)*GeoTrf::TranslateZ3D(-SecPztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
706 GeoSerialTransformer *sptA = new GeoSerialTransformer(sectorplates,&spA, NCrates);
707 GeoSerialTransformer *sptC = new GeoSerialTransformer(sectorplates,&spC, NCrates);
708 envelope->add(sptA);
709 envelope->add(sptC);
710
711 //bridgePlates
712 TRANSFUNCTION brpA1 = Pow(GeoTrf::RotateZ3D(1.0),f-(5*Alfa/2))*GeoTrf::TranslateX3D(BridgePxtr)*GeoTrf::TranslateZ3D(BridgePztr)*GeoTrf::RotateZ3D(90*Gaudi::Units::deg)*GeoTrf::RotateY3D(90*Gaudi::Units::deg)*GeoTrf::RotateX3D(90*Gaudi::Units::deg);
713 TRANSFUNCTION brpA2 = Pow(GeoTrf::RotateZ3D(1.0),f+(13*Alfa/2))*GeoTrf::TranslateX3D(BridgePxtr)*GeoTrf::TranslateZ3D(BridgePztr)*GeoTrf::RotateZ3D(90*Gaudi::Units::deg)*GeoTrf::RotateY3D(90*Gaudi::Units::deg)*GeoTrf::RotateX3D(90*Gaudi::Units::deg);
714 TRANSFUNCTION brpC1 = Pow(GeoTrf::RotateZ3D(1.0),f-(5*Alfa/2))*GeoTrf::TranslateX3D(BridgePxtr)*GeoTrf::TranslateZ3D(-BridgePztr)*GeoTrf::RotateZ3D(-90*Gaudi::Units::deg)*GeoTrf::RotateY3D(-90*Gaudi::Units::deg)*GeoTrf::RotateX3D(-90*Gaudi::Units::deg);
715 TRANSFUNCTION brpC2 = Pow(GeoTrf::RotateZ3D(1.0),f+(13*Alfa/2))*GeoTrf::TranslateX3D(BridgePxtr)*GeoTrf::TranslateZ3D(-BridgePztr)*GeoTrf::RotateZ3D(-90*Gaudi::Units::deg)*GeoTrf::RotateY3D(-90*Gaudi::Units::deg)*GeoTrf::RotateX3D(-90*Gaudi::Units::deg); GeoSerialTransformer *brptA1 = new GeoSerialTransformer(bridgeplates,&brpA1, 5);
716 GeoSerialTransformer *brptA2 = new GeoSerialTransformer(bridgeplates,&brpA2, 5);
717 GeoSerialTransformer *brptC1 = new GeoSerialTransformer(bridgeplates,&brpC1, 5);
718 GeoSerialTransformer *brptC2 = new GeoSerialTransformer(bridgeplates,&brpC2, 5);
719 envelope->add(brptA1);
720 envelope->add(brptA2);
721 envelope->add(brptC1);
722 envelope->add(brptC2);
723
724 //basePlates
725 TRANSFUNCTION bpA = Pow(GeoTrf::RotateZ3D(1.0),f-(Alfa/2))*GeoTrf::TranslateX3D(BasePxtr)*GeoTrf::TranslateZ3D(BasePztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
726 TRANSFUNCTION bpC = Pow(GeoTrf::RotateZ3D(1.0),f+(Alfa/2))*GeoTrf::TranslateX3D(BasePxtr)*GeoTrf::TranslateZ3D(-BasePztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
727 GeoSerialTransformer *bptA = new GeoSerialTransformer(baseplates,&bpA, NCrates);
728 GeoSerialTransformer *bptC = new GeoSerialTransformer(baseplates,&bpC, NCrates);
729 envelope->add(bptA);
730 envelope->add(bptC);
731
732 //sectorEnvelopes1
733 //counter-clockwise from top if taking sideA for reference (clockwise for sideC)
734 TRANSFUNCTION seA1G5 = Pow(GeoTrf::RotateZ3D(1.0),f+(9*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
735 TRANSFUNCTION seC1G5 = Pow(GeoTrf::RotateZ3D(1.0),f+(9*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
736 TRANSFUNCTION seA1G6 = Pow(GeoTrf::RotateZ3D(1.0),f+(11*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
737 TRANSFUNCTION seC1G6 = Pow(GeoTrf::RotateZ3D(1.0),f+(11*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
738 TRANSFUNCTION seA1G7 = Pow(GeoTrf::RotateZ3D(1.0),f+(13*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
739 TRANSFUNCTION seC1G7 = Pow(GeoTrf::RotateZ3D(1.0),f+(13*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
740 TRANSFUNCTION seA1G8 = Pow(GeoTrf::RotateZ3D(1.0),f+(15*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
741 TRANSFUNCTION seC1G8 = Pow(GeoTrf::RotateZ3D(1.0),f+(15*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
742 TRANSFUNCTION seA1G9 = Pow(GeoTrf::RotateZ3D(1.0),f+(17*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
743 TRANSFUNCTION seC1G9 = Pow(GeoTrf::RotateZ3D(1.0),f+(17*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
744 TRANSFUNCTION seA1G10 = Pow(GeoTrf::RotateZ3D(1.0),f+(19*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
745 TRANSFUNCTION seC1G10 = Pow(GeoTrf::RotateZ3D(1.0),f+(19*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
746 TRANSFUNCTION seA1G11 = Pow(GeoTrf::RotateZ3D(1.0),f+(21*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
747 TRANSFUNCTION seC1G11 = Pow(GeoTrf::RotateZ3D(1.0),f+(21*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
748 TRANSFUNCTION seA1G12 = Pow(GeoTrf::RotateZ3D(1.0),f+(23*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
749 TRANSFUNCTION seC1G12 = Pow(GeoTrf::RotateZ3D(1.0),f+(23*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
750 //clockwise from top if taking sideA for reference (counter-clockwise for sideC)
751 TRANSFUNCTION seA1G4 = Pow(GeoTrf::RotateZ3D(1.0),f+(7*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
752 TRANSFUNCTION seC1G4 = Pow(GeoTrf::RotateZ3D(1.0),f+(7*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
753 TRANSFUNCTION seA1G3 = Pow(GeoTrf::RotateZ3D(1.0),f+(5*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
754 TRANSFUNCTION seC1G3 = Pow(GeoTrf::RotateZ3D(1.0),f+(5*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
755 TRANSFUNCTION seA1G2 = Pow(GeoTrf::RotateZ3D(1.0),f+(3*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
756 TRANSFUNCTION seC1G2 = Pow(GeoTrf::RotateZ3D(1.0),f+(3*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
757 TRANSFUNCTION seA1G1 = Pow(GeoTrf::RotateZ3D(1.0),f+(1*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
758 TRANSFUNCTION seC1G1 = Pow(GeoTrf::RotateZ3D(1.0),f+(1*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
759 TRANSFUNCTION seA1G16 = Pow(GeoTrf::RotateZ3D(1.0),f-(1*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
760 TRANSFUNCTION seC1G16 = Pow(GeoTrf::RotateZ3D(1.0),f-(1*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
761 TRANSFUNCTION seA1G15 = Pow(GeoTrf::RotateZ3D(1.0),f-(3*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
762 TRANSFUNCTION seC1G15 = Pow(GeoTrf::RotateZ3D(1.0),f-(3*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
763 TRANSFUNCTION seA1G14 = Pow(GeoTrf::RotateZ3D(1.0),f-(5*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
764 TRANSFUNCTION seC1G14 = Pow(GeoTrf::RotateZ3D(1.0),f-(5*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
765 TRANSFUNCTION seA1G13 = Pow(GeoTrf::RotateZ3D(1.0),f-(7*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
766 TRANSFUNCTION seC1G13 = Pow(GeoTrf::RotateZ3D(1.0),f-(7*Alfa/2))*GeoTrf::TranslateX3D((Endb+Enda)/2)*GeoTrf::TranslateZ3D(-SecE1ztr)*GeoTrf::RotateY3D(90*Gaudi::Units::deg);
767
768 GeoSerialTransformer *setA1G5 = new GeoSerialTransformer(sectorenvelopes1g5,&seA1G5, 1);
769 GeoSerialTransformer *setC1G5 = new GeoSerialTransformer(sectorenvelopes1g5,&seC1G5, 1);
770 GeoSerialTransformer *setA1G6 = new GeoSerialTransformer(sectorenvelopes1g6,&seA1G6, 1);
771 GeoSerialTransformer *setC1G6 = new GeoSerialTransformer(sectorenvelopes1g6,&seC1G6, 1);
772 GeoSerialTransformer *setA1G7 = new GeoSerialTransformer(sectorenvelopes1g7,&seA1G7, 1);
773 GeoSerialTransformer *setC1G7 = new GeoSerialTransformer(sectorenvelopes1g7,&seC1G7, 1);
774 GeoSerialTransformer *setA1G8 = new GeoSerialTransformer(sectorenvelopes1g8,&seA1G8, 1);
775 GeoSerialTransformer *setC1G8 = new GeoSerialTransformer(sectorenvelopes1g8,&seC1G8, 1);
776 GeoSerialTransformer *setA1G9 = new GeoSerialTransformer(sectorenvelopes1g9,&seA1G9, 1);
777 GeoSerialTransformer *setC1G9 = new GeoSerialTransformer(sectorenvelopes1g9,&seC1G9, 1);
778 GeoSerialTransformer *setA1G10 = new GeoSerialTransformer(sectorenvelopes1g10,&seA1G10, 1);
779 GeoSerialTransformer *setC1G10 = new GeoSerialTransformer(sectorenvelopes1g10,&seC1G10, 1);
780 GeoSerialTransformer *setA1G11 = new GeoSerialTransformer(sectorenvelopes1g11,&seA1G11, 1);
781 GeoSerialTransformer *setC1G11 = new GeoSerialTransformer(sectorenvelopes1g11,&seC1G11, 1);
782 GeoSerialTransformer *setA1G12 = new GeoSerialTransformer(sectorenvelopes1g12,&seA1G12, 1);
783 GeoSerialTransformer *setC1G12 = new GeoSerialTransformer(sectorenvelopes1g12,&seC1G12, 1);
784
785 GeoSerialTransformer *setA1G4 = new GeoSerialTransformer(sectorenvelopes1g4,&seA1G4, 1);
786 GeoSerialTransformer *setC1G4 = new GeoSerialTransformer(sectorenvelopes1g4,&seC1G4, 1);
787 GeoSerialTransformer *setA1G3 = new GeoSerialTransformer(sectorenvelopes1g3,&seA1G3, 1);
788 GeoSerialTransformer *setC1G3 = new GeoSerialTransformer(sectorenvelopes1g3,&seC1G3, 1);
789 GeoSerialTransformer *setA1G2 = new GeoSerialTransformer(sectorenvelopes1g2,&seA1G2, 1);
790 GeoSerialTransformer *setC1G2 = new GeoSerialTransformer(sectorenvelopes1g2,&seC1G2, 1);
791 GeoSerialTransformer *setA1G1 = new GeoSerialTransformer(sectorenvelopes1g1,&seA1G1, 1);
792 GeoSerialTransformer *setC1G1 = new GeoSerialTransformer(sectorenvelopes1g1,&seC1G1, 1);
793 GeoSerialTransformer *setA1G16 = new GeoSerialTransformer(sectorenvelopes1g16,&seA1G16, 1);
794 GeoSerialTransformer *setC1G16 = new GeoSerialTransformer(sectorenvelopes1g16,&seC1G16, 1);
795 GeoSerialTransformer *setA1G15 = new GeoSerialTransformer(sectorenvelopes1g15,&seA1G15, 1);
796 GeoSerialTransformer *setC1G15 = new GeoSerialTransformer(sectorenvelopes1g15,&seC1G15, 1);
797 GeoSerialTransformer *setA1G14 = new GeoSerialTransformer(sectorenvelopes1g14,&seA1G14, 1);
798 GeoSerialTransformer *setC1G14 = new GeoSerialTransformer(sectorenvelopes1g14,&seC1G14, 1);
799 GeoSerialTransformer *setA1G13 = new GeoSerialTransformer(sectorenvelopes1g13,&seA1G13, 1);
800 GeoSerialTransformer *setC1G13 = new GeoSerialTransformer(sectorenvelopes1g13,&seC1G13, 1);
801
802 envelope->add(setA1G5);
803 envelope->add(setC1G5);
804 envelope->add(setA1G6);
805 envelope->add(setC1G6);
806 envelope->add(setA1G7);
807 envelope->add(setC1G7);
808 envelope->add(setA1G8);
809 envelope->add(setC1G8);
810 envelope->add(setA1G9);
811 envelope->add(setC1G9);
812 envelope->add(setA1G10);
813 envelope->add(setC1G10);
814 envelope->add(setA1G11);
815 envelope->add(setC1G11);
816 envelope->add(setA1G12);
817 envelope->add(setC1G12);
818
819 envelope->add(setA1G4);
820 envelope->add(setC1G4);
821 envelope->add(setA1G3);
822 envelope->add(setC1G3);
823 envelope->add(setA1G2);
824 envelope->add(setC1G2);
825 envelope->add(setA1G1);
826 envelope->add(setC1G1);
827 envelope->add(setA1G16);
828 envelope->add(setC1G16);
829 envelope->add(setA1G15);
830 envelope->add(setC1G15);
831 envelope->add(setA1G14);
832 envelope->add(setC1G14);
833 envelope->add(setA1G13);
834 envelope->add(setC1G13);
835
836 createSectorEnvelopes2FromDB (envelope, materialManager,
837 trdMap, *BarrelDMTrds,
838 trapMap, *BarrelDMTraps,
839 boxMap, *BarrelDMBoxes,
840 f, Box);
841 createBridgeEnvelopesFromDB (envelope, trapMap, *BarrelDMTraps,
842 matLArServices8, f);
843 createBaseEnvelopesFromDB (envelope, trdMap, *BarrelDMTrds,
844 matLArServices8, f);
845}
846
847
848CrackDMConstruction::CrackDMConstruction(IRDBAccessSvc* rdbAccess
849 , IGeoModelSvc* geoModel
850 , StoredMaterialManager* materialManager
851 , bool activateFT)
852 : AthMessaging("CrackDMConstruction")
853 , m_rdbAccess(rdbAccess)
854 , m_geoModel(geoModel)
855 , m_materialManager(materialManager)
856 , m_activateFT(activateFT)
857{
858}
859
860void CrackDMConstruction::create(GeoFullPhysVol* envelope)
861{
862 createFromDB(envelope
863 , m_rdbAccess
864 , m_geoModel
865 , m_materialManager);
866
867 if(m_activateFT){
868 std::string name = "LAr::Barrel::SignalFT::";
869 ATH_MSG_DEBUG("creating " << name << " volumes");
870
871 const GeoMaterial* iron = m_materialManager->getMaterial("std::Iron");
872
873 const double wflange_height = 37.*Gaudi::Units::mm;
874 const double wflange_R = 0.5*360.*Gaudi::Units::mm;
875 const GeoMaterial* wflange_mat = m_materialManager->getMaterial("LAr::FT::WarmFlange");
876
877 GeoShape* wflange = new GeoTube(0., wflange_R, wflange_height/2);
878 GeoLogVol* wflangeLV = new GeoLogVol(name + "WarmFlange", wflange, wflange_mat);
879 GeoPhysVol* wflangePV = new GeoPhysVol(wflangeLV);
880
881 const double bellow_height = 225.*Gaudi::Units::mm;
882 const double bellow_Router = 0.5*299.*Gaudi::Units::mm; // this also to be cut in warm wall
883 const double bellow_wall = 15.*Gaudi::Units::mm;
884 const GeoMaterial* bellow_mat = m_materialManager->getMaterial("LAr::FT::Bellow");
885
886 const double bellow_Rinner = bellow_Router - bellow_wall;
887 GeoShape* bellow = new GeoTube(bellow_Rinner, bellow_Router, bellow_height/2);
888 GeoLogVol* bellowLV = new GeoLogVol(name + "Bellow", bellow, bellow_mat);
889 GeoPhysVol* bellowPV = new GeoPhysVol(bellowLV);
890
891 const GeoMaterial* vcables_mat = m_materialManager->getMaterial("LAr::FT::VacuumCables");
892
893 GeoShape* vcables = new GeoTube(0., bellow_Rinner, bellow_height/2);
894 GeoLogVol* vcablesLV = new GeoLogVol(name + "VacuumCables", vcables, vcables_mat);
895 GeoPhysVol* vcablesPV = new GeoPhysVol(vcablesLV);
896
897 const double cflange_height = 35.*Gaudi::Units::mm;
898 const double cflange_Router = 0.5*283.*Gaudi::Units::mm;
899 const GeoMaterial* cflange_mat = m_materialManager->getMaterial("LAr::FT::ColdFlange");
900 GeoShape* cflange = new GeoTube(0., cflange_Router, cflange_height/2);
901 GeoLogVol* cflangeLV = new GeoLogVol(name + "ColdFlange", cflange, cflange_mat);
902 GeoPhysVol* cflangePV = new GeoPhysVol(cflangeLV);
903
904 const double coldbox1_Router = cflange_Router;
905 const double coldbox1_wall = 0.134*2.54*Gaudi::Units::cm;
906 const double coldbox1_height = 90.*Gaudi::Units::mm;
907 const double coldbox2_height = 16.*Gaudi::Units::mm;
908 const double hole_r = 0.5*133.*Gaudi::Units::mm;
909 const double hole_shift = -31.*Gaudi::Units::mm;
910 const double coldbox3_Router = 0.5*140.*Gaudi::Units::mm; // this also to be cut in cold wall
911 const double barrel_dist_from_outer_warm_wall_to_inner_cold_wall =
912 500.*Gaudi::Units::mm;
913 const double coldbox3_height = // adjust to match dist between cryo walls
914 barrel_dist_from_outer_warm_wall_to_inner_cold_wall
915 - bellow_height - coldbox1_height - coldbox2_height;
916 ATH_MSG_DEBUG("funnel tube len = " << coldbox3_height / Gaudi::Units::mm << " mm ");
917 const GeoMaterial* coldbox_mat = iron;
918 GeoShape* coldbox1 = new GeoTube(coldbox1_Router - coldbox1_wall, coldbox1_Router, coldbox1_height/2); // wide part
919 GeoShape* coldbox11 = new GeoTube(0., coldbox1_Router, coldbox1_height/2); // wide part for FTenvelope
920 GeoLogVol* coldbox1LV = new GeoLogVol(name + "ColdBox1", coldbox1, coldbox_mat);
921 GeoPhysVol* coldbox1PV = new GeoPhysVol(coldbox1LV);
922 GeoShape* coldbox21 = new GeoTube(0., coldbox1_Router, coldbox2_height/2); // plate
923 GeoShape* coldbox22 = new GeoTube(0., hole_r, coldbox2_height); // hole in the plate
924 const GeoShape& coldbox2 = coldbox21->subtract((*coldbox22) << GeoTrf::TranslateY3D(hole_shift));
925 GeoLogVol* coldbox2LV = new GeoLogVol(name + "ColdBox2", &coldbox2, coldbox_mat);
926 GeoPhysVol* coldbox2PV = new GeoPhysVol(coldbox2LV);
927 GeoShape* coldbox3 = new GeoTube(hole_r, coldbox3_Router, coldbox3_height/2); // narrow part
928 GeoShape* coldbox31 = new GeoTube(0., coldbox3_Router, coldbox3_height/2); // narrow part for FTenvelope
929 GeoLogVol* coldbox3LV = new GeoLogVol(name + "ColdBox3", coldbox3, coldbox_mat);
930 GeoPhysVol* coldbox3PV = new GeoPhysVol(coldbox3LV);
931
932 GeoTrf::TranslateZ3D bellow_pos(-wflange_height/2 - bellow_height/2);
933 GeoTrf::TranslateZ3D cflange_pos(-wflange_height/2 - bellow_height - cflange_height/2);
934 GeoTrf::TranslateZ3D coldbox1_pos(-wflange_height/2 - bellow_height - cflange_height - coldbox1_height/2);
935 GeoTrf::TranslateZ3D coldbox2_pos(-wflange_height/2 - bellow_height - cflange_height - coldbox1_height - coldbox2_height/2);
936 GeoTrf::Translate3D coldbox2hole_pos(
937 0., hole_shift,
938 -wflange_height/2 - bellow_height - cflange_height - coldbox1_height - coldbox2_height/2
939 );
940 GeoTrf::Translate3D coldbox3_pos(
941 0., hole_shift,
942 -wflange_height/2 - bellow_height - cflange_height - coldbox1_height - coldbox2_height - coldbox3_height/2
943 );
944
945 const GeoShape& FTenvelope = wflange->add(
946 (*bellow) << bellow_pos
947 ).add(
948 (*cflange) << cflange_pos
949 ).add(
950 (*coldbox11) << coldbox1_pos
951 ).add(
952 (*coldbox21) << coldbox2_pos
953 ).add(
954 (*coldbox31) << coldbox3_pos
955 );
956
957 GeoLogVol* FTLV = new GeoLogVol(name + "Envelope",
958 &FTenvelope,
959 m_materialManager->getMaterial("std::Air")
960 );
961 GeoPhysVol* FTPV = new GeoPhysVol(FTLV);
962
963 FTPV->add(wflangePV);
964 GeoTransform *bellow_trf = new GeoTransform(bellow_pos);
965 FTPV->add(bellow_trf);
966 FTPV->add(bellowPV);
967 FTPV->add(bellow_trf);
968 FTPV->add(vcablesPV);
969 GeoTransform *cflange_trf = new GeoTransform(cflange_pos);
970 FTPV->add(cflange_trf);
971 FTPV->add(cflangePV);
972 GeoTransform *coldbox1_trf = new GeoTransform(coldbox1_pos);
973 FTPV->add(coldbox1_trf);
974 FTPV->add(coldbox1PV);
975 GeoTransform *coldbox2_trf = new GeoTransform(coldbox2_pos);
976 FTPV->add(coldbox2_trf);
977 FTPV->add(coldbox2PV);
978 GeoTransform *coldbox3_trf = new GeoTransform(coldbox3_pos);
979 FTPV->add(coldbox3_trf);
980 FTPV->add(coldbox3PV);
981
982 const GeoMaterial* lar_mat = m_materialManager->getMaterial("std::LiquidArgon");
983 GeoShape* lar1 = new GeoTube(0., coldbox1_Router - coldbox1_wall, coldbox1_height/2);
984 GeoShape* lar2 = new GeoTube(0., hole_r, coldbox1_height);
985 const GeoShape& lar = lar1->subtract((*lar2) << GeoTrf::TranslateY3D(hole_shift));
986 GeoLogVol* larLV = new GeoLogVol(name + "LAr", &lar, lar_mat);
987 GeoPhysVol* larPV = new GeoPhysVol(larLV);
988 FTPV->add(coldbox1_trf);
989 FTPV->add(larPV);
990
991 const GeoMaterial *pigtail_mat = m_materialManager->getMaterial("LAr::FT::Pigtail");
992
993 const double pth = (coldbox1_height + coldbox2_height + coldbox3_height) / 2;
994 GeoTransform *pigtail_trf = new GeoTransform(
995 GeoTrf::Translate3D(0, hole_shift,
996 -wflange_height/2 - bellow_height - cflange_height - pth
997 ));
998 GeoShape* pigtail = new GeoTube(0., hole_r, pth);
999 GeoLogVol* pigtailLV = new GeoLogVol(name + "Pigtail", pigtail, pigtail_mat);
1000 GeoPhysVol* pigtailPV = new GeoPhysVol(pigtailLV);
1001 FTPV->add(pigtail_trf);
1002 FTPV->add(pigtailPV);
1003
1004 // total lenght should be 28.5 cm
1005 const double ocable_len = 10.*Gaudi::Units::cm;
1006 const double ocable_R = (1.1/2)*sqrt(1920*2.85)*Gaudi::Units::mm;
1007 const GeoMaterial* ocable_mat = m_materialManager->getMaterial("LAr::FT::Cable");
1008 GeoShape* ocable = new GeoTube(0., ocable_R, ocable_len / 2);
1009 GeoLogVol* ocableLV = new GeoLogVol("LAr::Barrel::FTCables", ocable, ocable_mat);
1010 GeoPhysVol* ocablePV = new GeoPhysVol(ocableLV);
1011
1012 // todo: take cryostat parameters from DB
1013 const double r0 = 277.5*Gaudi::Units::cm // cryo warm wall outer radius
1014 + wflange_height/2; // warm flange is abobe the wall
1015 const double r2 = r0 + wflange_height/2 + ocable_len/2;
1016 const double z_pos = 3170.*Gaudi::Units::mm;
1017 const int NCrates = 16;
1018 const double dphi = 4.*Gaudi::Units::deg;
1019
1020 auto put1 = [&envelope](GeoPhysVol *object, double r, double phi, double z)
1021 {
1022 envelope->add(new GeoTransform(
1023 GeoTrf::Translate3D(r*cos(phi), r*sin(phi), z) *
1024 GeoTrf::RotateX3D(90.*Gaudi::Units::deg) *
1025 GeoTrf::RotateY3D(phi + 90.*Gaudi::Units::deg)
1026 ));
1027 if(z < 0){
1028 envelope->add(new GeoTransform(GeoTrf::RotateZ3D(180.*Gaudi::Units::deg)));
1029 }
1030 envelope->add(object);
1031 };
1032
1033 auto put = [&put1, &FTPV, &ocablePV, r0, r2](double phi, double z)
1034 {
1035 put1(FTPV, r0, phi, z);
1036 put1(ocablePV, r2, phi, z);
1037 };
1038
1039 for(int i = 0; i < 16; ++ i){
1040 const double phi = 360.*Gaudi::Units::deg / NCrates * i;
1041 put(phi - dphi, z_pos);
1042 put(phi + dphi, z_pos);
1043 put(phi - dphi, -z_pos);
1044 put(phi + dphi, -z_pos);
1045 }
1046 }
1047
1048}
1049
M_PI
#define M_PI
Definition ActiveFraction.h:11
Box
Acts::Volume::BoundingBox Box
Definition ActsCaloTrackingVolumeBuilder.cxx:28
eta
Scalar eta() const
pseudorapidity method
Definition AmgMatrixBasePlugin.h:83
deltaR
Scalar deltaR(const MatrixBase< Derived > &vec) const
Definition AmgMatrixBasePlugin.h:105
phi
Scalar phi() const
phi method
Definition AmgMatrixBasePlugin.h:67
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition AthMsgStreamMacros.h:29
Endc
static const double Endc
Definition CrackDMConstruction.cxx:55
Enda
static const double Enda
Definition CrackDMConstruction.cxx:53
DYa
static const double DYa
Definition CrackDMConstruction.cxx:56
createBridgeEnvelopesFromDB
void createBridgeEnvelopesFromDB(GeoFullPhysVol *envelope, std::map< std::string, unsigned int > &trapMap, IRDBRecordset &BarrelDMTraps, const GeoMaterial *matLArServices8, GENFUNCTION &f)
Definition CrackDMConstruction.cxx:270
createBaseEnvelopesFromDB
void createBaseEnvelopesFromDB(GeoFullPhysVol *envelope, std::map< std::string, unsigned int > &trdMap, IRDBRecordset &BarrelDMTrds, const GeoMaterial *matLArServices8, GENFUNCTION &f)
Definition CrackDMConstruction.cxx:306
NCrates
static const unsigned int NCrates
Definition CrackDMConstruction.cxx:51
DYc
static const double DYc
Definition CrackDMConstruction.cxx:58
createFromDB
void createFromDB(GeoFullPhysVol *envelope, IRDBAccessSvc *rdbAccess, IGeoModelSvc *geoModel, StoredMaterialManager *materialManager)
Definition CrackDMConstruction.cxx:334
Endb
static const double Endb
Definition CrackDMConstruction.cxx:54
DYb
static const double DYb
Definition CrackDMConstruction.cxx:57
createSectorEnvelopes2FromDB
void createSectorEnvelopes2FromDB(GeoFullPhysVol *envelope, StoredMaterialManager *materialManager, std::map< std::string, unsigned int > &trdMap, IRDBRecordset &BarrelDMTrds, std::map< std::string, unsigned int > &trapMap, IRDBRecordset &BarrelDMTraps, std::map< std::string, unsigned int > &boxMap, IRDBRecordset &BarrelDMBoxes, GENFUNCTION &f, GeoBox *Box)
Definition CrackDMConstruction.cxx:62
Alfa
static const double Alfa
Definition CrackDMConstruction.cxx:52
IRDBAccessSvc.h
Definition of the abstract IRDBAccessSvc interface.
IRDBRecordset_ptr
std::shared_ptr< IRDBRecordset > IRDBRecordset_ptr
Definition IRDBAccessSvc.h:27
IRDBRecord.h
Definition of the abstract IRDBRecord interface.
IRDBRecordset.h
Definition of the abstract IRDBRecordset interface.
z
#define z
AthMessaging::AthMessaging
AthMessaging(IMessageSvc *msgSvc, const std::string &name)
Constructor.
Definition AthMessaging.cxx:13
CrackDMConstruction::m_rdbAccess
IRDBAccessSvc * m_rdbAccess
Definition CrackDMConstruction.h:27
CrackDMConstruction::CrackDMConstruction
CrackDMConstruction()=delete
CrackDMConstruction::create
void create(GeoFullPhysVol *envelope)
Definition CrackDMConstruction.cxx:860
CrackDMConstruction::m_materialManager
StoredMaterialManager * m_materialManager
Definition CrackDMConstruction.h:29
DecodeVersionKey
This is a helper class to query the version tags from GeoModelSvc and determine the appropriate tag a...
Definition DecodeVersionKey.h:18
DecodeVersionKey::tag
const std::string & tag() const
Return version tag.
Definition DecodeVersionKey.cxx:91
DecodeVersionKey::node
const std::string & node() const
Return the version node.
Definition DecodeVersionKey.cxx:97
IRDBAccessSvc
IRDBAccessSvc is an abstract interface to the athena service that provides the following functionalit...
Definition IRDBAccessSvc.h:42
IRDBAccessSvc::getRecordsetPtr
virtual IRDBRecordset_ptr getRecordsetPtr(const std::string &node, const std::string &tag, const std::string &tag2node="", const std::string &connName="ATLASDD")=0
Provides access to the Recordset object containing HVS-tagged data.
IRDBRecord
IRDBRecord is one record in the IRDBRecordset object.
Definition IRDBRecord.h:27
IRDBRecordset
IRDBRecordset is an abstract interface to the RDB Recordsets, which represent a snapshot of HVS-tagge...
Definition IRDBRecordset.h:35
IRDBRecordset::size
virtual unsigned int size() const =0
StoredMaterialManager
This class holds one or more material managers and makes them storeable, under StoreGate.
Definition StoredMaterialManager.h:28
StoredMaterialManager::getMaterial
virtual const GeoMaterial * getMaterial(const std::string &name)=0
Variable
Wrapper around a histogram which allows for some additional filling patterns and data manipulation.
Definition Trigger/TrigCost/TrigCostAnalysis/src/Variable.h:39
add
bool add(const std::string &hname, TKey *tobj)
Definition fastadd.cxx:55
r
int r
Definition globals.cxx:22
object
Generated on for ATLAS Offline Software by doxygen 1.14.0