Belle II Software  release-08-01-10
Ph1sustrCreator.cc
1 /**************************************************************************
2  * basf2 (Belle II Analysis Software Framework) *
3  * Author: The Belle II Collaboration *
4  * *
5  * See git log for contributors and copyright holders. *
6  * This file is licensed under LGPL-3.0, see LICENSE.md. *
7  **************************************************************************/
8 
9 #include <beast/ph1sustr/geometry/Ph1sustrCreator.h>
10 #include <beast/ph1sustr/simulation/SensitiveDetector.h>
11 
12 #include <geometry/Materials.h>
13 #include <geometry/CreatorFactory.h>
14 #include <framework/gearbox/GearDir.h>
15 #include <framework/logging/Logger.h>
16 
17 #include <cmath>
18 #include <boost/format.hpp>
19 #include <boost/foreach.hpp>
20 #include <boost/algorithm/string.hpp>
21 
22 #include <G4LogicalVolume.hh>
23 #include <G4PVPlacement.hh>
24 
25 //Shapes
26 #include <G4Box.hh>
27 #include "G4UnionSolid.hh"
28 #include "G4SubtractionSolid.hh"
29 
30 //Visualization
31 #include "G4Colour.hh"
32 #include <G4VisAttributes.hh>
33 
34 using namespace std;
35 using namespace boost;
36 
37 namespace Belle2 {
44  namespace ph1sustr {
45 
46  // Register the creator
48  geometry::CreatorFactory<Ph1sustrCreator> Ph1sustrFactory("PH1SUSTRCreator");
49 
50  Ph1sustrCreator::Ph1sustrCreator(): m_sensitive(0)
51  {
52  //m_sensitive = new SensitiveDetector();
53  }
54 
55  Ph1sustrCreator::~Ph1sustrCreator()
56  {
57  if (m_sensitive) delete m_sensitive;
58  }
59 
60  void Ph1sustrCreator::create(const GearDir& content, G4LogicalVolume& topVolume, geometry::GeometryTypes /* type */)
61  {
62 
63  m_sensitive = new SensitiveDetector();
64 
65  //lets get the stepsize parameter with a default value of 5 µm
66  //double stepSize = content.getLength("stepSize", 5 * CLHEP::um);
67 
68 
69  //no get the array. Notice that the default framework unit is cm, so the
70  //values will be automatically converted
71  vector<double> bar = content.getArray("bar");
72  B2INFO("Contents of bar: ");
73  BOOST_FOREACH(double value, bar) {
74  B2INFO("value: " << value);
75  }
76  /*double x_tpcbeamR = 0;
77  double y_tpcbeamR = 0;
78  double z_tpcbeamR = 0;
79  double x_tpcbeamL = 0;
80  double y_tpcbeamL = 0;
81  double z_tpcbeamL = 0;
82  double x_tpcbeamT = 0;
83  double y_tpcbeamT = 0;
84  double z_tpcbeamT = 0;
85  double x_tpcbeamB = 0;
86  double y_tpcbeamB = 0;
87  double z_tpcbeamB = 0;*/
88  //Lets loop over all the Active nodes
89  BOOST_FOREACH(const GearDir & activeParams, content.getNodes("Active")) {
90 
91  //plate positions
92  double x_tpcbeamR = activeParams.getLength("x_tpcbeamR") * CLHEP::cm;
93  double y_tpcbeamR = activeParams.getLength("y_tpcbeamR") * CLHEP::cm;
94  double z_tpcbeamR = activeParams.getLength("z_tpcbeamR") * CLHEP::cm;
95  double x_tpcbeamL = activeParams.getLength("x_tpcbeamL") * CLHEP::cm;
96  double y_tpcbeamL = activeParams.getLength("y_tpcbeamL") * CLHEP::cm;
97  double z_tpcbeamL = activeParams.getLength("z_tpcbeamL") * CLHEP::cm;
98  double x_tpcbeamT = activeParams.getLength("x_tpcbeamT") * CLHEP::cm;
99  double y_tpcbeamT = activeParams.getLength("y_tpcbeamT") * CLHEP::cm;
100  double z_tpcbeamT = activeParams.getLength("z_tpcbeamT") * CLHEP::cm;
101  double x_tpcbeamB = activeParams.getLength("x_tpcbeamB") * CLHEP::cm;
102  double y_tpcbeamB = activeParams.getLength("y_tpcbeamB") * CLHEP::cm;
103  double z_tpcbeamB = activeParams.getLength("z_tpcbeamB") * CLHEP::cm;
104 
105  //TPC vertical: 4x @ 1614/ea
106  //TPC horizontal + BGO base: 8x @ 1583/ea
107  //TPC railroad: 8x @ 2200/ea
108  //BGO vertical: 4x @ 928/ea
109  //BGO horizontal: 4x @ 318/ea
110  //G4double dz_20V2100bgov = 843.72 / 2.*CLHEP::mm;
111  //G4double dz_20V2100bgoh = 280.00 / 2.*CLHEP::mm;
112 
113  //Beam supporting the TPC-Tube-plate
114  //define tpc beam and plate dimensions
115  double betpcbeam = 190.8 / 2. * CLHEP::mm;
116  G4double dx_tpcbeam = 2.54 * 1.63 / 2.*CLHEP::cm;
117  G4double dy_tpcbeam = 2.54 * 1.63 / 2.*CLHEP::cm;
118  G4double dz_tpcbeam = 2200. / 2.*CLHEP::mm;
119  G4double dw_tpcbeam = 2.54 * 0.25 / 2.*CLHEP::cm;
120  //G4double dx_plate = 2.54 * 0.35 / 2.*CLHEP::cm;
121 
122  //G4double dx_plate_short = 0.4765 * CLHEP::cm;
123  G4double dx_plate_short = 0.5 / 2. * CLHEP::cm;
124  G4double dy_plate_short = 27.47788 / 2.*CLHEP::cm;
125  G4double dz_plate_short = 40. / 2.*CLHEP::cm;
126 
127  G4double dx_plate = 0.5 / 2. * CLHEP::cm;
128  G4double dy_plate = 32. / 2.*CLHEP::cm;
129  G4double dz_plate = 50. / 2.*CLHEP::cm;
130 
131  //Right from e^-~--~beam
132  //create plate volume
133  G4VSolid* s_plate = new G4Box("s_plate", dx_plate, dy_plate, dz_plate);
134  G4VSolid* s_plate_short = new G4Box("s_plate_short", dx_plate_short, dy_plate_short, dz_plate_short);
135 
136  //place plate volume
137  G4LogicalVolume* l_plate = new G4LogicalVolume(s_plate, geometry::Materials::get("Al"), "l_plate", 0, 0);
138  G4LogicalVolume* l_plate_short = new G4LogicalVolume(s_plate_short, geometry::Materials::get("Al"), "l_plate_short", 0, 0);
139  G4VisAttributes* white = new G4VisAttributes(G4Colour(1, 1, 1));
140  white->SetForceAuxEdgeVisible(true);
141  l_plate->SetVisAttributes(white);
142  l_plate_short->SetVisAttributes(white);
143 
144  G4VSolid* s_tpcbeam_a = new G4Box("s_tpcbeam_a", dx_tpcbeam, dy_tpcbeam, dz_tpcbeam);
145  G4VSolid* s_tpcbeam_b = new G4Box("s_tpcbeam_b", dx_tpcbeam - 2.*dw_tpcbeam, dy_tpcbeam - dw_tpcbeam, dz_tpcbeam);
146  G4VSolid* s_tpcbeampos = new G4SubtractionSolid("s_tpcbeampos", s_tpcbeam_a, s_tpcbeam_b, 0, G4ThreeVector(0, dw_tpcbeam, 0));
147  G4VSolid* s_tpcbeamneg = new G4SubtractionSolid("s_tpcbeamneg", s_tpcbeam_a, s_tpcbeam_b, 0, G4ThreeVector(0, -dw_tpcbeam, 0));
148  G4VSolid* s_tpcbeam = new G4UnionSolid("s_tpcbeam", s_tpcbeampos, s_tpcbeamneg, 0, G4ThreeVector(0, -2.*dy_tpcbeam, 0));
149 
150  //create tpc beam volumes
151  G4LogicalVolume* l_tpcbeam = new G4LogicalVolume(s_tpcbeam, geometry::Materials::get("FG_Epoxy"), "l_tpcbeam", 0, 0);
152 
153  //place plate volume
154  G4ThreeVector PH1SUSTRpos = G4ThreeVector(
155  x_tpcbeamR,
156  y_tpcbeamR,
157  z_tpcbeamR
158  );
159  new G4PVPlacement(0, PH1SUSTRpos, l_plate, "p_plateR", &topVolume, false, 1);
160 
161  //offsets hori
162  //G4double offset_h = /*110.*CLHEP::cm*/dz_tpcbeam + /*80.*CLHEP::cm*/ dz_plate - 2. * dz_tpcbeam;
163  G4double offset_h = /*110.*CLHEP::cm*/dz_tpcbeam + 80.*CLHEP::cm - 2. * dz_tpcbeam;
164 
165  //place 1st tpc beam volume
166  /*PH1SUSTRpos = G4ThreeVector(
167  x_tpcbeamR + dx_plate + dx_tpcbeam,
168  y_tpcbeamR + betpcbeam + 2. * dy_tpcbeam,
169  -offset_h / 2.
170  );
171  new G4PVPlacement(0, PH1SUSTRpos, l_tpcbeam, "p_tpcbeam", &topVolume, false, 1);*/
172  G4Transform3D TransForm = G4Translate3D(x_tpcbeamR + dx_plate + dx_tpcbeam,
173  y_tpcbeamR + betpcbeam + 2. * dy_tpcbeam,
174  -offset_h / 2.) * G4RotateZ3D(90.*CLHEP::deg);
175  new G4PVPlacement(TransForm, l_tpcbeam, "p_tpcbeamR1", &topVolume, false, 1);
176 
177  //place 2nd tpc beam volume
178  /*PH1SUSTRpos = G4ThreeVector(
179  x_tpcbeamR + dx_plate + dx_tpcbeam,
180  y_tpcbeamR - betpcbeam - 2. * dy_tpcbeam,
181  -offset_h / 2.
182  );
183  new G4PVPlacement(0, PH1SUSTRpos, l_tpcbeam, "p_tpcbeam", &topVolume, false, 1);*/
184  TransForm = G4Translate3D(x_tpcbeamR + dx_plate + dx_tpcbeam,
185  y_tpcbeamR - betpcbeam - 2. * dy_tpcbeam,
186  -offset_h / 2.) * G4RotateZ3D(90.*CLHEP::deg);
187  new G4PVPlacement(TransForm, l_tpcbeam, "p_tpcbeamR2", &topVolume, false, 1);
188 
189  //Left from e^-~--~beam
190  //place plate volume
191  PH1SUSTRpos = G4ThreeVector(
192  x_tpcbeamL,
193  y_tpcbeamL,
194  z_tpcbeamL
195  );
196  new G4PVPlacement(0, PH1SUSTRpos, l_plate, "p_plateL", &topVolume, false, 1);
197 
198  //place 1st tpc beam volume
199  /*PH1SUSTRpos = G4ThreeVector(
200  x_tpcbeamL - dx_plate - dx_tpcbeam,
201  y_tpcbeamL + betpcbeam + 2. * dy_tpcbeam,
202  -offset_h / 2.
203  );
204  new G4PVPlacement(0, PH1SUSTRpos, l_tpcbeam, "p_tpcbeam", &topVolume, false, 1);*/
205  TransForm = G4Translate3D(x_tpcbeamL - dx_plate - 3. * dx_tpcbeam,
206  y_tpcbeamL + betpcbeam + 2. * dy_tpcbeam,
207  -offset_h / 2.) * G4RotateZ3D(90.*CLHEP::deg);
208  new G4PVPlacement(TransForm, l_tpcbeam, "p_tpcbeamL1", &topVolume, false, 1);
209 
210  //place 2nd tpc beam volume
211  /*PH1SUSTRpos = G4ThreeVector(
212  x_tpcbeamL - dx_plate - dx_tpcbeam,
213  y_tpcbeamL - betpcbeam - 2. * dy_tpcbeam,
214  -offset_h / 2.
215  );
216  new G4PVPlacement(0, PH1SUSTRpos, l_tpcbeam, "p_tpcbeam", &topVolume, false, 1);*/
217  TransForm = G4Translate3D(x_tpcbeamL - dx_plate - 3. * dx_tpcbeam,
218  y_tpcbeamL - betpcbeam - 2. * dy_tpcbeam,
219  -offset_h / 2.) * G4RotateZ3D(90.*CLHEP::deg);
220  new G4PVPlacement(TransForm, l_tpcbeam, "p_tpcbeamL2", &topVolume, false, 1);
221 
222  //Bottom
223  G4RotationMatrix* rotXx = new G4RotationMatrix();
224  double Angle = 90. * CLHEP::deg;
225  rotXx->rotateZ(Angle);
226  //place bottom plate
227  PH1SUSTRpos = G4ThreeVector(
228  x_tpcbeamB,
229  y_tpcbeamB,
230  z_tpcbeamB
231  );
232  new G4PVPlacement(rotXx, PH1SUSTRpos, l_plate_short, "p_plateB", &topVolume, false, 1);
233 
234  //place 1st tpc beam volume
235  /*PH1SUSTRpos = G4ThreeVector(
236  x_tpcbeamB + betpcbeam + 2. * dx_tpcbeam,
237  y_tpcbeamB - dx_plate - dy_tpcbeam,
238  -offset_h / 2.
239  );
240  new G4PVPlacement(rotXx, PH1SUSTRpos, l_tpcbeam, "p_tpcbeam", &topVolume, false, 1);*/
241  TransForm = G4Translate3D(x_tpcbeamB + betpcbeam + 1.5 * dx_tpcbeam,
242  y_tpcbeamB - dx_plate - 1.5 * dy_tpcbeam,
243  -offset_h / 2.) /* G4RotateZ3D(90.*CLHEP::deg)*/;
244  new G4PVPlacement(TransForm, l_tpcbeam, "p_tpcbeamB1", &topVolume, false, 1);
245 
246  //place 2nd tpc beam volume
247  /*PH1SUSTRpos = G4ThreeVector(
248  x_tpcbeamB - betpcbeam - 2. * dx_tpcbeam,
249  y_tpcbeamB - dx_plate - dy_tpcbeam,
250  -offset_h / 2.
251  );
252  new G4PVPlacement(rotXx, PH1SUSTRpos, l_tpcbeam, "p_tpcbeam", &topVolume, false, 1);*/
253  TransForm = G4Translate3D(x_tpcbeamB - betpcbeam - 1.5 * dx_tpcbeam,
254  y_tpcbeamB - dx_plate - 1.5 * dy_tpcbeam,
255  -offset_h / 2.) /* G4RotateZ3D(90.*CLHEP::deg)*/;
256  new G4PVPlacement(TransForm, l_tpcbeam, "p_tpcbeamB2", &topVolume, false, 1);
257 
258  //Top
259  //rotXx = new G4RotationMatrix();
260  //rotXx->rotateZ( 90. );
261  //place top plate
262  PH1SUSTRpos = G4ThreeVector(
263  x_tpcbeamT,
264  y_tpcbeamT,
265  z_tpcbeamT
266  );
267  new G4PVPlacement(rotXx, PH1SUSTRpos, l_plate, "p_plateT", &topVolume, false, 1);
268 
269  //place 1st tpc beam volume
270  /*PH1SUSTRpos = G4ThreeVector(
271  x_tpcbeamT + betpcbeam + 2. * dx_tpcbeam,
272  y_tpcbeamT + dx_plate + dy_tpcbeam,
273  -offset_h / 2.
274  );
275  new G4PVPlacement(rotXx, PH1SUSTRpos, l_tpcbeam, "p_tpcbeam", &topVolume, false, 1);*/
276  TransForm = G4Translate3D(x_tpcbeamT + betpcbeam + 1.5 * dx_tpcbeam,
277  y_tpcbeamT + dx_plate + 3. * dy_tpcbeam,
278  -offset_h / 2.) /* G4RotateZ3D(90.*CLHEP::deg)*/;
279  new G4PVPlacement(TransForm, l_tpcbeam, "p_tpcbeamT1", &topVolume, false, 1);
280 
281  //place 2nd tpc beam volume
282  /*PH1SUSTRpos = G4ThreeVector(
283  x_tpcbeamT - betpcbeam - 2. * dx_tpcbeam,
284  y_tpcbeamT + dx_plate + dy_tpcbeam,
285  -offset_h / 2.
286  );
287  new G4PVPlacement(rotXx, PH1SUSTRpos, l_tpcbeam, "p_tpcbeam", &topVolume, false, 1);*/
288  TransForm = G4Translate3D(x_tpcbeamT - betpcbeam - 1.5 * dx_tpcbeam,
289  y_tpcbeamT + dx_plate + 3 * dy_tpcbeam,
290  -offset_h / 2.) /* G4RotateZ3D(90.*CLHEP::deg)*/;
291  new G4PVPlacement(TransForm, l_tpcbeam, "p_tpcbeamT2", &topVolume, false, 1);
292 
293  //vertical beams
294  //G4double dz_tpcbeamv = 1537.37 / 2.*CLHEP::mm;
295  G4double dz_tpcbeamv = 1583. / 2.*CLHEP::mm;
296  G4VSolid* s_tpcbeamv_a = new G4Box("s_tpcbeamv_a", dx_tpcbeam, dy_tpcbeam, dz_tpcbeamv);
297  G4VSolid* s_tpcbeamv_b = new G4Box("s_tpcbeamv_b", dx_tpcbeam - 2.*dw_tpcbeam, dy_tpcbeam - dw_tpcbeam, dz_tpcbeamv);
298  G4VSolid* s_tpcbeamvpos = new G4SubtractionSolid("s_tpcbeamvpos", s_tpcbeamv_a, s_tpcbeamv_b, 0, G4ThreeVector(0, dw_tpcbeam, 0));
299  G4VSolid* s_tpcbeamvneg = new G4SubtractionSolid("s_tpcbeamvneg", s_tpcbeamv_a, s_tpcbeamv_b, 0, G4ThreeVector(0, -dw_tpcbeam, 0));
300  G4VSolid* s_tpcbeamv = new G4UnionSolid("s_tpcbeamv", s_tpcbeamvpos, s_tpcbeamvneg, 0, G4ThreeVector(0, -2.*dy_tpcbeam, 0));
301  G4LogicalVolume* l_tpcbeamv = new G4LogicalVolume(s_tpcbeamv, geometry::Materials::get("FG_Epoxy"), "l_tpcbeamv", 0, 0);
302 
303  //offset verti
304  //G4double offset_v = fabs(76.*CLHEP::cm - 2. * dz_tpcbeamv) / 2.;
305  G4double offset_v = fabs(110.*CLHEP::cm - 2. * dz_tpcbeamv) / 2.;
306 
307  //place 1st vertical TPC beam
308  G4RotationMatrix* rotX = new G4RotationMatrix();
309  rotX->rotateX(90.*CLHEP::deg);
310  PH1SUSTRpos = G4ThreeVector(
311  x_tpcbeamL - dx_plate - 3 * dx_tpcbeam - 2.*dx_tpcbeam,
312  offset_v / 2.,
313  -800.*CLHEP::mm
314  );
315  new G4PVPlacement(rotX, PH1SUSTRpos, l_tpcbeamv, "p_tpcbeamv1", &topVolume, false, 0);
316  //place 2nd vertical TPC beam
317  PH1SUSTRpos = G4ThreeVector(
318  x_tpcbeamR + dx_plate + 3.*dx_tpcbeam + 2.*dx_tpcbeam,
319  offset_v / 2.,
320  -800.*CLHEP::mm
321  );
322  new G4PVPlacement(rotX, PH1SUSTRpos, l_tpcbeamv, "p_tpcbeamv2", &topVolume, false, 0);
323  //place 3rd vertical TPC beam
324  PH1SUSTRpos = G4ThreeVector(
325  x_tpcbeamL - dx_plate - 3.*dx_tpcbeam - 2.*dx_tpcbeam,
326  offset_v / 2.,
327  1100.*CLHEP::mm
328  );
329  new G4PVPlacement(rotX, PH1SUSTRpos, l_tpcbeamv, "p_tpcbeamv3", &topVolume, false, 0);
330  //place 4th vertical TPC beam
331  PH1SUSTRpos = G4ThreeVector(
332  x_tpcbeamR + dx_plate + 3.*dx_tpcbeam + 2.*dx_tpcbeam,
333  offset_v / 2.,
334  1100.*CLHEP::mm
335  );
336  new G4PVPlacement(rotX, PH1SUSTRpos, l_tpcbeamv, "p_tpcbeamv4", &topVolume, false, 0);
337 
338  //horizontal beams
339  //G4double dz_tpcbeamh = 1792. / 2.*CLHEP::mm;
340  G4double dz_tpcbeamh = 1614. / 2.*CLHEP::mm;
341  G4VSolid* s_tpcbeamh_a = new G4Box("s_tpcbeamh_a", dx_tpcbeam, dy_tpcbeam, dz_tpcbeamh);
342  G4VSolid* s_tpcbeamh_b = new G4Box("s_tpcbeamh_b", dx_tpcbeam - 2.*dw_tpcbeam, dy_tpcbeam - dw_tpcbeam, dz_tpcbeamh);
343  G4VSolid* s_tpcbeamhpos = new G4SubtractionSolid("s_tpcbeamhpos", s_tpcbeamh_a, s_tpcbeamh_b, 0, G4ThreeVector(0, dw_tpcbeam, 0));
344  G4VSolid* s_tpcbeamhneg = new G4SubtractionSolid("s_tpcbeanhneg", s_tpcbeamh_a, s_tpcbeamh_b, 0, G4ThreeVector(0, -dw_tpcbeam, 0));
345  G4VSolid* s_tpcbeamh = new G4UnionSolid("s_tpcbeamh", s_tpcbeamhpos, s_tpcbeamhneg, 0, G4ThreeVector(0, -2.*dy_tpcbeam, 0));
346 
347  G4LogicalVolume* l_tpcbeamh = new G4LogicalVolume(s_tpcbeamh, geometry::Materials::get("FG_Epoxy"), "l_tpcbeamh", 0, 0);
348 
349  //place 1st horizontal TPC beam
350  G4RotationMatrix* rotY = new G4RotationMatrix();
351  rotY->rotateY(90.*CLHEP::deg);
352  PH1SUSTRpos = G4ThreeVector(
353  0 * CLHEP::mm,
354  y_tpcbeamB - dx_plate - 4. * dy_tpcbeam - 2. * dy_tpcbeam,
355  1100.*CLHEP::mm - 2. * dy_tpcbeam
356  );
357  new G4PVPlacement(rotY, PH1SUSTRpos, l_tpcbeamh, "p_tpcbeamh1", &topVolume, false, 0);
358 
359  //place 2nd horizontal TPC beam
360  PH1SUSTRpos = G4ThreeVector(
361  0 * CLHEP::mm,
362  y_tpcbeamT + dx_plate + 5. * dy_tpcbeam + 2. * dy_tpcbeam,
363  1100.*CLHEP::mm - 2. * dy_tpcbeam
364  );
365  new G4PVPlacement(rotY, PH1SUSTRpos, l_tpcbeamh, "p_tpcbeamh2", &topVolume, false, 0);
366 
367  //place 3rd horizontal TPC beam
368  PH1SUSTRpos = G4ThreeVector(
369  0 * CLHEP::mm,
370  y_tpcbeamB - dx_plate - 4. * dy_tpcbeam - 2. * dy_tpcbeam,
371  -800.*CLHEP::mm + 4. * dy_tpcbeam
372  );
373  new G4PVPlacement(rotY, PH1SUSTRpos, l_tpcbeamh, "p_tpcbeamh3", &topVolume, false, 0);
374 
375  //place 4th horizontal TPC beam
376  PH1SUSTRpos = G4ThreeVector(
377  0 * CLHEP::mm,
378  y_tpcbeamT + dx_plate + 5. * dy_tpcbeam + 2. * dy_tpcbeam,
379  -800.*CLHEP::mm + 4. * dy_tpcbeam
380  );
381  new G4PVPlacement(rotY, PH1SUSTRpos, l_tpcbeamh, "p_tpcbeamh4", &topVolume, false, 0);
382 
383  G4VisAttributes* brown = new G4VisAttributes(G4Colour(.5, .5, 0));
384  brown->SetForceAuxEdgeVisible(true);
385  l_tpcbeam->SetVisAttributes(brown);
386  l_tpcbeamv->SetVisAttributes(brown);
387  l_tpcbeamh->SetVisAttributes(brown);
388 
389  //CsI box beams
390  G4double dz_csibeamh = activeParams.getLength("lcsiBeamh") * CLHEP::cm / 2.;
391  G4VSolid* s_csibeamh_a = new G4Box("s_csibeamh_a", dx_tpcbeam, dy_tpcbeam, dz_csibeamh);
392  G4VSolid* s_csibeamh_b = new G4Box("s_csibeamh_b", dx_tpcbeam - 2.*dw_tpcbeam, dy_tpcbeam - dw_tpcbeam, dz_csibeamh);
393  G4VSolid* s_csibeamhpos = new G4SubtractionSolid("s_csibeamhpos", s_csibeamh_a, s_csibeamh_b, 0, G4ThreeVector(0, dw_tpcbeam, 0));
394  G4VSolid* s_csibeamhneg = new G4SubtractionSolid("s_csibeamhneg", s_csibeamh_a, s_csibeamh_b, 0, G4ThreeVector(0, -dw_tpcbeam, 0));
395  //G4VSolid* s_csibeamhpos = new G4SubtractionSolid("s_csibeamhpos", s_csibeamh_a, s_csibeamh_b, 0, G4ThreeVector(dw_tpcbeam, 0, 0));
396  //G4VSolid* s_csibeamhneg = new G4SubtractionSolid("s_csibeamhneg", s_csibeamh_a, s_csibeamh_b, 0, G4ThreeVector(-dw_tpcbeam, 0, 0));
397  G4VSolid* s_csibeamh = new G4UnionSolid("s_csibeamh", s_csibeamhpos, s_csibeamhneg, 0, G4ThreeVector(0, -2.*dy_tpcbeam, 0));
398  G4LogicalVolume* l_csibeamh = new G4LogicalVolume(s_csibeamh, geometry::Materials::get("FG_Epoxy"), "l_csibeamh", 0, 0);
399  int xdimbase = 0;
400  int ydimbase = 0;
401  int zdimbase = 0;
402  int xdimcsiBeamh = 0;
403  int ydimcsiBeamh = 0;
404  int zdimcsiBeamh = 0;
405  int xdimcsiBeamvF = 0;
406  int ydimcsiBeamvF = 0;
407  int zdimcsiBeamvF = 0;
408  int xdimcsiBeamvB = 0;
409  int ydimcsiBeamvB = 0;
410  int zdimcsiBeamvB = 0;
411  int xdimbgobeamv = 0;
412  int ydimbgobeamv = 0;
413  int zdimbgobeamv = 0;
414  int xdimbgobeamh = 0;
415  int ydimbgobeamh = 0;
416  int zdimbgobeamh = 0;
417  int xdimbgobeamt = 0;
418  int ydimbgobeamt = 0;
419  int zdimbgobeamt = 0;
420  int xdimbgobeamb = 0;
421  int ydimbgobeamb = 0;
422  int zdimbgobeamb = 0;
423  int xdimtpcbeamb = 0;
424  int ydimtpcbeamb = 0;
425  int zdimtpcbeamb = 0;
426  /*int xdimtpcbeamhh = 0;
427  int ydimtpcbeamhh = 0;
428  int zdimtpcbeamhh = 0;
429  int xdimtpcbeamvv = 0;
430  int ydimtpcbeamvv = 0;
431  int zdimtpcbeamvv = 0;*/
432  double xbase[100];
433  double ybase[100];
434  double zbase[100];
435  double xcsibeamh[100];
436  double ycsibeamh[100];
437  double zcsibeamh[100];
438  double xcsibeamvF[100];
439  double ycsibeamvF[100];
440  double zcsibeamvF[100];
441  double xcsibeamvB[100];
442  double ycsibeamvB[100];
443  double zcsibeamvB[100];
444  double xbgobeamv[100];
445  double ybgobeamv[100];
446  double zbgobeamv[100];
447  double xbgobeamh[100];
448  double ybgobeamh[100];
449  double zbgobeamh[100];
450  double xbgobeamt[100];
451  double ybgobeamt[100];
452  double zbgobeamt[100];
453  double xbgobeamb[100];
454  double ybgobeamb[100];
455  double zbgobeamb[100];
456  double xtpcbeamb[100];
457  double ytpcbeamb[100];
458  double ztpcbeamb[100];
459  /*double xtpcbeamhh[100];
460  double ytpcbeamhh[100];
461  double ztpcbeamhh[100];
462  double xtpcbeamvv[100];
463  double ytpcbeamvv[100];
464  double ztpcbeamvv[100];*/
465 
466  double x_offset = activeParams.getLength("x_offset") * CLHEP::cm;
467  double y_offset = activeParams.getLength("y_offset") * CLHEP::cm;
468 
469  for (double xcsiBeamh : activeParams.getArray("xcsiBeamh", {0})) {
470  xcsiBeamh *= CLHEP::cm;
471  xcsibeamh[xdimcsiBeamh] = xcsiBeamh - x_offset;
472  xdimcsiBeamh++;
473  }
474  for (double ycsiBeamh : activeParams.getArray("ycsiBeamh", {0})) {
475  ycsiBeamh *= CLHEP::cm;
476  ycsibeamh[ydimcsiBeamh] = ycsiBeamh - y_offset;
477  ydimcsiBeamh++;
478  }
479  for (double zcsiBeamh : activeParams.getArray("zcsiBeamh", {0})) {
480  zcsiBeamh *= CLHEP::cm;
481  zcsibeamh[zdimcsiBeamh] = zcsiBeamh;
482  zdimcsiBeamh++;
483  }
484  for (int i = 0; i < xdimcsiBeamh; i++) {
485  G4Transform3D transform = G4Translate3D(xcsibeamh[i], ycsibeamh[i] - dy_tpcbeam,
486  zcsibeamh[i] + 2. * dy_tpcbeam) * G4RotateY3D(90.*CLHEP::deg) * G4RotateZ3D(90.*CLHEP::deg);
487  new G4PVPlacement(transform, l_csibeamh, TString::Format("p_csibeamh1_%d", i).Data(), &topVolume, false, 1);
488  transform = G4Translate3D(xcsibeamh[i], ycsibeamh[i] + dy_tpcbeam,
489  zcsibeamh[i] + 2. * dy_tpcbeam) * G4RotateY3D(90.*CLHEP::deg) * G4RotateZ3D(90.*CLHEP::deg);
490  new G4PVPlacement(transform, l_csibeamh, TString::Format("p_csibeamh2_%d", i).Data(), &topVolume, false, 1);
491  //PH1SUSTRpos = G4ThreeVector(xcsibeamh[i],ycsibeamh[i],zcsibeamh[i] - dy_tpcbeam);
492  //new G4PVPlacement(rotY, PH1SUSTRpos, l_csibeamh, "p_csibeamh", &topVolume, false, 0);
493  //PH1SUSTRpos = G4ThreeVector(xcsibeamh[i],ycsibeamh[i],zcsibeamh[i] + dy_tpcbeam);
494  //new G4PVPlacement(rotY+rotXx, PH1SUSTRpos, l_csibeamh, "p_csibeamh", &topVolume, false, 0);
495  }
496 
497  G4double dz_csibeamvF = activeParams.getLength("lcsiBeamvF") * CLHEP::cm / 2.;
498  G4VSolid* s_csibeamvF_a = new G4Box("s_csibeamvF_a", dx_tpcbeam, dy_tpcbeam, dz_csibeamvF);
499  G4VSolid* s_csibeamvF_b = new G4Box("s_csibeamvF_b", dx_tpcbeam - 2.*dw_tpcbeam, dy_tpcbeam - dw_tpcbeam, dz_csibeamvF);
500  G4VSolid* s_csibeamvFpos = new G4SubtractionSolid("s_csibeamvFpos", s_csibeamvF_a, s_csibeamvF_b, 0, G4ThreeVector(0, dw_tpcbeam,
501  0));
502  G4VSolid* s_csibeamvFneg = new G4SubtractionSolid("s_csibeamvFneg", s_csibeamvF_a, s_csibeamvF_b, 0, G4ThreeVector(0, -dw_tpcbeam,
503  0));
504  G4VSolid* s_csibeamvF = new G4UnionSolid("s_csibeamvF", s_csibeamvFpos, s_csibeamvFneg, 0, G4ThreeVector(0, -2.*dy_tpcbeam, 0));
505  G4LogicalVolume* l_csibeamvF = new G4LogicalVolume(s_csibeamvF, geometry::Materials::get("FG_Epoxy"), "l_csibeamvF", 0, 0);
506 
507  for (double xcsiBeamvF : activeParams.getArray("xcsiBeamvF", {0})) {
508  xcsiBeamvF *= CLHEP::cm;
509  xcsibeamvF[xdimcsiBeamvF] = xcsiBeamvF - x_offset;
510  xdimcsiBeamvF++;
511  }
512  for (double ycsiBeamvF : activeParams.getArray("ycsiBeamvF", {0})) {
513  ycsiBeamvF *= CLHEP::cm;
514  ycsibeamvF[ydimcsiBeamvF] = ycsiBeamvF - y_offset ;
515  ydimcsiBeamvF++;
516  }
517  for (double zcsiBeamvF : activeParams.getArray("zcsiBeamvF", {0})) {
518  zcsiBeamvF *= CLHEP::cm;
519  zcsibeamvF[zdimcsiBeamvF] = zcsiBeamvF ;
520  zdimcsiBeamvF++;
521  }
522  for (int i = 0; i < xdimcsiBeamvF; i++) {
523  //PH1SUSTRpos = G4ThreeVector(xcsibeamvF[i],0,zcsibeamh[0] - 2. * dx_tpcbeam);
524  PH1SUSTRpos = G4ThreeVector(xcsibeamvF[i], ycsibeamvF[i], zcsibeamvF[i]);
525  new G4PVPlacement(rotX, PH1SUSTRpos, l_csibeamvF, TString::Format("p_csibeamvF_%d", i).Data(), &topVolume, false, 0);
526  }
527 
528  G4double dz_csibeamvB = activeParams.getLength("lcsiBeamvB") * CLHEP::cm / 2.;
529  G4VSolid* s_csibeamvB_a = new G4Box("s_csibeamvB_a", dx_tpcbeam, dy_tpcbeam, dz_csibeamvB);
530  G4VSolid* s_csibeamvB_b = new G4Box("s_csibeamvB_b", dx_tpcbeam - 2.*dw_tpcbeam, dy_tpcbeam - dw_tpcbeam, dz_csibeamvB);
531  G4VSolid* s_csibeamvBpos = new G4SubtractionSolid("s_csibeamvBpos", s_csibeamvB_a, s_csibeamvB_b, 0, G4ThreeVector(0, dw_tpcbeam,
532  0));
533  G4VSolid* s_csibeamvBneg = new G4SubtractionSolid("s_csibeamvBneg", s_csibeamvB_a, s_csibeamvB_b, 0, G4ThreeVector(0, -dw_tpcbeam,
534  0));
535  G4VSolid* s_csibeamvB = new G4UnionSolid("s_csibeamvB", s_csibeamvBpos, s_csibeamvBneg, 0, G4ThreeVector(0, -2.*dy_tpcbeam, 0));
536  G4LogicalVolume* l_csibeamvB = new G4LogicalVolume(s_csibeamvB, geometry::Materials::get("FG_Epoxy"), "l_csibeamvB", 0, 0);
537 
538  for (double xcsiBeamvB : activeParams.getArray("xcsiBeamvB", {0})) {
539  xcsiBeamvB *= CLHEP::cm;
540  xcsibeamvB[xdimcsiBeamvB] = xcsiBeamvB - x_offset;
541  xdimcsiBeamvB++;
542  }
543  for (double ycsiBeamvB : activeParams.getArray("ycsiBeamvB", {0})) {
544  ycsiBeamvB *= CLHEP::cm;
545  ycsibeamvB[ydimcsiBeamvB] = ycsiBeamvB - y_offset ;
546  ydimcsiBeamvB++;
547  }
548  for (double zcsiBeamvB : activeParams.getArray("zcsiBeamvB", {0})) {
549  zcsiBeamvB *= CLHEP::cm;
550  zcsibeamvB[zdimcsiBeamvB] = zcsiBeamvB ;
551  zdimcsiBeamvB++;
552  }
553  for (int i = 0; i < xdimcsiBeamvB; i++) {
554  PH1SUSTRpos = G4ThreeVector(xcsibeamvB[i], ycsibeamvB[i], zcsibeamvB[i]);
555  new G4PVPlacement(rotX, PH1SUSTRpos, l_csibeamvB, TString::Format("p_csibeamvB_%d", i).Data(), &topVolume, false, 0);
556  }
557 
558  G4double dz_base = activeParams.getLength("lBase") * CLHEP::cm / 2.;
559  G4VSolid* s_base_a = new G4Box("s_base_a", dx_tpcbeam, dy_tpcbeam, dz_base);
560  G4VSolid* s_base_b = new G4Box("s_base_b", dx_tpcbeam - 2.*dw_tpcbeam, dy_tpcbeam - dw_tpcbeam, dz_base);
561  G4VSolid* s_base = new G4SubtractionSolid("s_base", s_base_a, s_base_b, 0, G4ThreeVector(0, dw_tpcbeam, 0));
562  G4LogicalVolume* l_base = new G4LogicalVolume(s_base, geometry::Materials::get("FG_Epoxy"), "l_base", 0, 0);
563 
564  for (double xBase : activeParams.getArray("xBase", {0})) {
565  xBase *= CLHEP::cm;
566  xbase[xdimbase] = xBase - x_offset;
567  xdimbase++;
568  }
569  for (double yBase : activeParams.getArray("yBase", {0})) {
570  yBase *= CLHEP::cm;
571  ybase[ydimbase] = yBase - y_offset;
572  ydimbase++;
573  }
574  for (double zBase : activeParams.getArray("zBase", {0})) {
575  zBase *= CLHEP::cm;
576  zbase[zdimbase] = zBase ;
577  zdimbase++;
578  }
579  for (int i = 0; i < xdimbase; i++) {
580  PH1SUSTRpos = G4ThreeVector(xbase[i], ybase[i], zbase[i]);
581  new G4PVPlacement(0, PH1SUSTRpos, l_base, TString::Format("p_base_%d", i).Data(), &topVolume, false, 0);
582  }
583 
584  G4double dz_bgobeamv = activeParams.getLength("lbgoBeamv") * CLHEP::cm / 2.;
585  G4VSolid* s_bgobeamv_a = new G4Box("s_bgobeamv_a", dx_tpcbeam, dy_tpcbeam, dz_bgobeamv);
586  G4VSolid* s_bgobeamv_b = new G4Box("s_bgobeamv_b", dx_tpcbeam - 2.*dw_tpcbeam, dy_tpcbeam - dw_tpcbeam, dz_bgobeamv);
587  G4VSolid* s_bgobeamvpos = new G4SubtractionSolid("s_bgobeamvpos", s_bgobeamv_a, s_bgobeamv_b, 0, G4ThreeVector(0, dw_tpcbeam, 0));
588  G4VSolid* s_bgobeamvneg = new G4SubtractionSolid("s_bgobeamvneg", s_bgobeamv_a, s_bgobeamv_b, 0, G4ThreeVector(0, -dw_tpcbeam, 0));
589  G4VSolid* s_bgobeamv = new G4UnionSolid("s_bgobeamv", s_bgobeamvpos, s_bgobeamvneg, 0, G4ThreeVector(0, -2.*dy_tpcbeam, 0));
590  G4LogicalVolume* l_bgobeamv = new G4LogicalVolume(s_bgobeamv, geometry::Materials::get("FG_Epoxy"), "l_bgobeamv", 0, 0);
591 
592  for (double xbgoBeamv : activeParams.getArray("xbgoBeamv", {0})) {
593  xbgoBeamv *= CLHEP::cm;
594  xbgobeamv[xdimbgobeamv] = xbgoBeamv;
595  xdimbgobeamv++;
596  }
597  for (double ybgoBeamv : activeParams.getArray("ybgoBeamv", {0})) {
598  ybgoBeamv *= CLHEP::cm;
599  ybgobeamv[ydimbgobeamv] = ybgoBeamv - y_offset;
600  ydimbgobeamv++;
601  }
602  for (double zbgoBeamv : activeParams.getArray("zbgoBeamv", {0})) {
603  zbgoBeamv *= CLHEP::cm;
604  zbgobeamv[zdimbgobeamv] = zbgoBeamv ;
605  zdimbgobeamv++;
606  }
607  for (int i = 0; i < xdimbgobeamv; i++) {
608  PH1SUSTRpos = G4ThreeVector(xbgobeamv[i], ybgobeamv[i], zbgobeamv[i]);
609  new G4PVPlacement(rotX, PH1SUSTRpos, l_bgobeamv, TString::Format("p_bgobeamv_%d", i).Data(), &topVolume, false, 0);
610  }
611 
612  G4double dz_bgobeamh = activeParams.getLength("lbgoBeamh") * CLHEP::cm / 2.;
613  G4VSolid* s_bgobeamh_a = new G4Box("s_bgobeamh_a", dx_tpcbeam, dy_tpcbeam, dz_bgobeamh);
614  G4VSolid* s_bgobeamh_b = new G4Box("s_bgobeamh_b", dx_tpcbeam - 2.*dw_tpcbeam, dy_tpcbeam - dw_tpcbeam, dz_bgobeamh);
615  G4VSolid* s_bgobeamhpos = new G4SubtractionSolid("s_bgobeamhpos", s_bgobeamh_a, s_bgobeamh_b, 0, G4ThreeVector(0, dw_tpcbeam, 0));
616  G4VSolid* s_bgobeamhneg = new G4SubtractionSolid("s_bgobeamhneg", s_bgobeamh_a, s_bgobeamh_b, 0, G4ThreeVector(0, -dw_tpcbeam, 0));
617  G4VSolid* s_bgobeamh = new G4UnionSolid("s_bgobeamh", s_bgobeamhpos, s_bgobeamhneg, 0, G4ThreeVector(0, -2.*dy_tpcbeam, 0));
618  G4LogicalVolume* l_bgobeamh = new G4LogicalVolume(s_bgobeamh, geometry::Materials::get("FG_Epoxy"), "l_bgobeamh", 0, 0);
619 
620  for (double xbgoBeamh : activeParams.getArray("xbgoBeamh", {0})) {
621  xbgoBeamh *= CLHEP::cm;
622  xbgobeamh[xdimbgobeamh] = xbgoBeamh;
623  xdimbgobeamh++;
624  }
625  for (double ybgoBeamh : activeParams.getArray("ybgoBeamh", {0})) {
626  ybgoBeamh *= CLHEP::cm;
627  ybgobeamh[ydimbgobeamh] = ybgoBeamh - y_offset;
628  ydimbgobeamh++;
629  }
630  for (double zbgoBeamh : activeParams.getArray("zbgoBeamh", {0})) {
631  zbgoBeamh *= CLHEP::cm;
632  zbgobeamh[zdimbgobeamh] = zbgoBeamh ;
633  zdimbgobeamh++;
634  }
635  for (int i = 0; i < xdimbgobeamh; i++) {
636  G4Transform3D transform = G4Translate3D(xbgobeamh[i], ybgobeamh[i],
637  zbgobeamh[i] + 2. * dy_tpcbeam) * G4RotateY3D(90.*CLHEP::deg) * G4RotateZ3D(90.*CLHEP::deg);
638  new G4PVPlacement(transform, l_bgobeamh, TString::Format("p_bgobeamh_%d", i).Data(), &topVolume, false, 0);
639  }
640 
641  G4double dz_bgobeamt = activeParams.getLength("lbgoBeamt") * CLHEP::cm / 2.;
642  G4VSolid* s_bgobeamt_a = new G4Box("s_bgobeamt_a", dx_tpcbeam, dy_tpcbeam, dz_bgobeamt);
643  G4VSolid* s_bgobeamt_b = new G4Box("s_bgobeamt_b", dx_tpcbeam - 2.*dw_tpcbeam, dy_tpcbeam - dw_tpcbeam, dz_bgobeamt);
644  G4VSolid* s_bgobeamt = new G4SubtractionSolid("s_bgobeamt", s_bgobeamt_a, s_bgobeamt_b, 0, G4ThreeVector(0, dw_tpcbeam, 0));
645  G4LogicalVolume* l_bgobeamt = new G4LogicalVolume(s_bgobeamt, geometry::Materials::get("FG_Epoxy"), "l_bgobeamt", 0, 0);
646 
647  for (double xbgoBeamt : activeParams.getArray("xbgoBeamt", {0})) {
648  xbgoBeamt *= CLHEP::cm;
649  xbgobeamt[xdimbgobeamt] = xbgoBeamt;
650  xdimbgobeamt++;
651  }
652  for (double ybgoBeamt : activeParams.getArray("ybgoBeamt", {0})) {
653  ybgoBeamt *= CLHEP::cm;
654  ybgobeamt[ydimbgobeamt] = ybgoBeamt - y_offset;
655  ydimbgobeamt++;
656  }
657  for (double zbgoBeamt : activeParams.getArray("zbgoBeamt", {0})) {
658  zbgoBeamt *= CLHEP::cm;
659  zbgobeamt[zdimbgobeamt] = zbgoBeamt ;
660  zdimbgobeamt++;
661  }
662  for (int i = 0; i < xdimbgobeamt; i++) {
663  PH1SUSTRpos = G4ThreeVector(xbgobeamt[i], ybgobeamt[i], zbgobeamt[i]);
664  new G4PVPlacement(0, PH1SUSTRpos, l_bgobeamt, TString::Format("p_bgobeamt_%d", i).Data(), &topVolume, false, 0);
665  }
666 
667  G4double dz_bgobeamb = activeParams.getLength("lbgoBeamb") * CLHEP::cm / 2.;
668  G4VSolid* s_bgobeamb_a = new G4Box("s_bgobeamb_a", dx_tpcbeam, dy_tpcbeam, dz_bgobeamb);
669  G4VSolid* s_bgobeamb_b = new G4Box("s_bgobeamb_b", dx_tpcbeam - 2.*dw_tpcbeam, dy_tpcbeam - dw_tpcbeam, dz_bgobeamb);
670  G4VSolid* s_bgobeambpos = new G4SubtractionSolid("s_bgobeambpos", s_bgobeamb_a, s_bgobeamb_b, 0, G4ThreeVector(0, dw_tpcbeam, 0));
671  G4VSolid* s_bgobeambneg = new G4SubtractionSolid("s_bgobeambneg", s_bgobeamb_a, s_bgobeamb_b, 0, G4ThreeVector(0, -dw_tpcbeam, 0));
672  G4VSolid* s_bgobeamb = new G4UnionSolid("s_bgobeamb", s_bgobeambpos, s_bgobeambneg, 0, G4ThreeVector(0, -2.*dy_tpcbeam, 0));
673  G4LogicalVolume* l_bgobeamb = new G4LogicalVolume(s_bgobeamb, geometry::Materials::get("FG_Epoxy"), "l_bgobeamb", 0, 0);
674  for (double xbgoBeamb : activeParams.getArray("xbgoBeamb", {0})) {
675  xbgoBeamb *= CLHEP::cm;
676  xbgobeamb[xdimbgobeamb] = xbgoBeamb - x_offset;
677  xdimbgobeamb++;
678  }
679  for (double ybgoBeamb : activeParams.getArray("ybgoBeamb", {0})) {
680  ybgoBeamb *= CLHEP::cm;
681  ybgobeamb[ydimbgobeamb] = ybgoBeamb - y_offset;
682  ydimbgobeamb++;
683  }
684  for (double zbgoBeamb : activeParams.getArray("zbgoBeamb", {0})) {
685  zbgoBeamb *= CLHEP::cm;
686  zbgobeamb[zdimbgobeamb] = zbgoBeamb ;
687  zdimbgobeamb++;
688  }
689  for (int i = 0; i < xdimbgobeamb; i++) {
690  G4Transform3D transform = G4Translate3D(xbgobeamb[i], ybgobeamb[i],
691  zbgobeamb[i] + 2. * dy_tpcbeam) * G4RotateY3D(90.*CLHEP::deg) * G4RotateZ3D(90.*CLHEP::deg);
692  new G4PVPlacement(transform, l_bgobeamb, TString::Format("p_bgobeamb_%d", i).Data(), &topVolume, false, 1);
693  }
694 
695  G4double dz_tpcbeamb = activeParams.getLength("ltpcBeamb") * CLHEP::cm / 2.;
696  G4VSolid* s_tpcbeamb_a = new G4Box("s_tpcbeamb_a", dx_tpcbeam, dy_tpcbeam, dz_tpcbeamb);
697  G4VSolid* s_tpcbeamb_b = new G4Box("s_tpcbeamb_b", dx_tpcbeam - 2.*dw_tpcbeam, dy_tpcbeam - dw_tpcbeam, dz_tpcbeamb);
698  G4VSolid* s_tpcbeambpos = new G4SubtractionSolid("s_tpcbeambpos", s_tpcbeamb_a, s_tpcbeamb_b, 0, G4ThreeVector(0, dw_tpcbeam, 0));
699  G4VSolid* s_tpcbeambneg = new G4SubtractionSolid("s_tpcbeambneg", s_tpcbeamb_a, s_tpcbeamb_b, 0, G4ThreeVector(0, -dw_tpcbeam, 0));
700  G4VSolid* s_tpcbeamb = new G4UnionSolid("s_tpcbeamb", s_tpcbeambpos, s_tpcbeambneg, 0, G4ThreeVector(0, -2.*dy_tpcbeam, 0));
701  G4LogicalVolume* l_tpcbeamb = new G4LogicalVolume(s_tpcbeamb, geometry::Materials::get("FG_Epoxy"), "l_tpcbeamb", 0, 0);
702  for (double xtpcBeamb : activeParams.getArray("xtpcBeamb", {0})) {
703  xtpcBeamb *= CLHEP::cm;
704  xtpcbeamb[xdimtpcbeamb] = xtpcBeamb - x_offset;
705  xdimtpcbeamb++;
706  }
707  for (double ytpcBeamb : activeParams.getArray("ytpcBeamb", {0})) {
708  ytpcBeamb *= CLHEP::cm;
709  ytpcbeamb[ydimtpcbeamb] = ytpcBeamb - y_offset;
710  ydimtpcbeamb++;
711  }
712  for (double ztpcBeamb : activeParams.getArray("ztpcBeamb", {0})) {
713  ztpcBeamb *= CLHEP::cm;
714  ztpcbeamb[zdimtpcbeamb] = ztpcBeamb ;
715  zdimtpcbeamb++;
716  }
717  for (int i = 0; i < xdimtpcbeamb; i++) {
718  G4Transform3D transform = G4Translate3D(xtpcbeamb[i], ytpcbeamb[i],
719  ztpcbeamb[i] + 2. * dy_tpcbeam) * G4RotateY3D(90.*CLHEP::deg) * G4RotateZ3D(90.*CLHEP::deg);
720  new G4PVPlacement(transform, l_tpcbeamb, TString::Format("p_tpcbeamb_%d", i).Data(), &topVolume, false, 1);
721  }
722 
723  /*
724  G4double dz_tpcbeamhh = activeParams.getLength("ltpcBeamh") * CLHEP::cm / 2.;
725  G4VSolid* s_tpcbeamhh_a = new G4Box("s_tpcbeamhh_a", dx_tpcbeam, dy_tpcbeam, dz_tpcbeamhh);
726  G4VSolid* s_tpcbeamhh_b = new G4Box("s_tpcbeamhh_b", dx_tpcbeam - 2.*dw_tpcbeam, dy_tpcbeam - dw_tpcbeam, dz_tpcbeamhh);
727  G4VSolid* s_tpcbeamhhpos = new G4SubtractionSolid("s_tpcbeamhhpos", s_tpcbeamhh_a, s_tpcbeamhh_b, 0, G4ThreeVector(0, dw_tpcbeam, 0));
728  G4VSolid* s_tpcbeamhhneg = new G4SubtractionSolid("s_tpcbeamhhneg", s_tpcbeamhh_a, s_tpcbeamhh_b, 0, G4ThreeVector(0, -dw_tpcbeam, 0));
729  G4VSolid* s_tpcbeamhh = new G4UnionSolid("s_tpcbeamhh", s_tpcbeamhhpos, s_tpcbeamhhneg, 0, G4ThreeVector(0, -2.*dy_tpcbeam, 0));
730  G4LogicalVolume* l_tpcbeamhh = new G4LogicalVolume(s_tpcbeamhh, geometry::Materials::get("FG_Epoxy") , "l_tpcbeamhh", 0, 0);
731 
732  for (double xtpcBeamh : activeParams.getArray("xtpcBeamh", {0})) {
733  xtpcBeamh *= CLHEP::cm;
734  xtpcbeamhh[xdimtpcbeamhh] = xtpcBeamh - x_offset;
735  xdimtpcbeamhh++;
736  }
737  for (double ytpcBeamh : activeParams.getArray("ytpcBeamh", {0})) {
738  ytpcBeamh *= CLHEP::cm;
739  ytpcbeamhh[ydimtpcbeamhh] = ytpcBeamh - y_offset;
740  ydimtpcbeamhh++;
741  }
742  for (double ztpcBeamh : activeParams.getArray("ztpcBeamh", {0})) {
743  ztpcBeamh *= CLHEP::cm;
744  ztpcbeamhh[zdimtpcbeamhh] = ztpcBeamh ;
745  zdimtpcbeamhh++;
746  }
747  for(int i = 0; i < xdimtpcbeamhh; i++) {
748  double z_tmp = 0;
749  if(i<1)z_tmp = ztpcbeamhh[i] + dy_tpcbeam/2.;
750  else z_tmp = ztpcbeamhh[i] - dy_tpcbeam *1.1;
751  G4Transform3D transform = G4Translate3D(xtpcbeamhh[i],ytpcbeamhh[i],z_tmp) * G4RotateY3D(90.*CLHEP::deg) ;
752  new G4PVPlacement(transform, l_tpcbeamhh, TString::Format("p_tpcbeamhh_%d",i).Data(), &topVolume, false, 1);
753  }
754 
755  G4double dz_tpcbeamvv = activeParams.getLength("ltpcBeamv") * CLHEP::cm / 2.;
756  G4VSolid* s_tpcbeamvv_a = new G4Box("s_tpcbeamvv_a", dx_tpcbeam, dy_tpcbeam, dz_tpcbeamvv);
757  G4VSolid* s_tpcbeamvv_b = new G4Box("s_tpcbeamvv_b", dx_tpcbeam - 2.*dw_tpcbeam, dy_tpcbeam - dw_tpcbeam, dz_tpcbeamvv);
758  G4VSolid* s_tpcbeamvvpos = new G4SubtractionSolid("s_tpcbeamvvpos", s_tpcbeamvv_a, s_tpcbeamvv_b, 0, G4ThreeVector(0, dw_tpcbeam, 0));
759  G4VSolid* s_tpcbeamvvneg = new G4SubtractionSolid("s_tpcbeamvvneg", s_tpcbeamvv_a, s_tpcbeamvv_b, 0, G4ThreeVector(0, -dw_tpcbeam, 0));
760  G4VSolid* s_tpcbeamvv = new G4UnionSolid("s_tpcbeamvv", s_tpcbeamvvpos, s_tpcbeamvvneg, 0, G4ThreeVector(0, -2.*dy_tpcbeam, 0));
761  G4LogicalVolume* l_tpcbeamvv = new G4LogicalVolume(s_tpcbeamvv, geometry::Materials::get("FG_Epoxy") , "l_tpcbeamvv", 0, 0);
762 
763  for (double xtpcBeamv : activeParams.getArray("xtpcBeamv", {0})) {
764  xtpcBeamv *= CLHEP::cm;
765  xtpcbeamvv[xdimtpcbeamvv] = xtpcBeamv - x_offset;
766  xdimtpcbeamvv++;
767  }
768  for (double ytpcBeamv : activeParams.getArray("ytpcBeamv", {0})) {
769  ytpcBeamv *= CLHEP::cm;
770  ytpcbeamvv[ydimtpcbeamvv] = ytpcBeamv - y_offset;
771  ydimtpcbeamvv++;
772  }
773  for (double ztpcBeamv : activeParams.getArray("ztpcBeamv", {0})) {
774  ztpcBeamv *= CLHEP::cm;
775  ztpcbeamvv[zdimtpcbeamvv] = ztpcBeamv ;
776  zdimtpcbeamvv++;
777  }
778  for(int i = 0; i < xdimtpcbeamvv; i++) {
779  //G4Transform3D transform = G4Translate3D(xtpcbeamvv[i],ytpcbeamvv[i],ztpcbeamvv[i]) * G4RotateY3D(90.*CLHEP::deg) ;
780  //new G4PVPlacement(transform, l_tpcbeamvv, TString::Format("p_tpcbeamvv_%d",i).Data(), &topVolume, false, 1);
781  PH1SUSTRpos = G4ThreeVector(xtpcbeamvv[i],ytpcbeamvv[i],ztpcbeamvv[i]);
782  new G4PVPlacement(rotX, PH1SUSTRpos, l_tpcbeamvv, TString::Format("p_tpcbeamvv_%d",i).Data(), &topVolume, false, 0);
783  }
784  */
785  G4double dz_fangsbeamhf = activeParams.getLength("lfangsBeamhf") * CLHEP::cm / 2.;
786  G4VSolid* s_fangsbeamhf_a = new G4Box("s_fangsbeamhf_a", dx_tpcbeam, dy_tpcbeam, dz_fangsbeamhf);
787  G4VSolid* s_fangsbeamhf_b = new G4Box("s_fangsbeamhf_b", dx_tpcbeam - 2.*dw_tpcbeam, dy_tpcbeam - dw_tpcbeam, dz_fangsbeamhf);
788  G4VSolid* s_fangsbeamhf = new G4SubtractionSolid("s_fangsbeamhf", s_fangsbeamhf_a, s_fangsbeamhf_b, 0, G4ThreeVector(0, dw_tpcbeam,
789  0));
790  G4LogicalVolume* l_fangsbeamhf = new G4LogicalVolume(s_fangsbeamhf, geometry::Materials::get("FG_Epoxy"), "l_fangsbeamhf", 0, 0);
791  double xfangsbeamhf = activeParams.getLength("xfangsBeamhf") * CLHEP::cm;
792  double yfangsbeamhf = activeParams.getLength("yfangsBeamhf") * CLHEP::cm;
793  double zfangsbeamhf = activeParams.getLength("zfangsBeamhf") * CLHEP::cm;
794  PH1SUSTRpos = G4ThreeVector(xfangsbeamhf, yfangsbeamhf, zfangsbeamhf);
795  new G4PVPlacement(rotY, PH1SUSTRpos, l_fangsbeamhf, "p_fangsbeamhf", &topVolume, false, 0);
796 
797  G4double dz_fangsbeamhb = activeParams.getLength("lfangsBeamhb") * CLHEP::cm / 2.;
798  G4VSolid* s_fangsbeamhb_a = new G4Box("s_fangsbeamhb_a", dx_tpcbeam, dy_tpcbeam, dz_fangsbeamhb);
799  G4VSolid* s_fangsbeamhb_b = new G4Box("s_fangsbeamhb_b", dx_tpcbeam - 2.*dw_tpcbeam, dy_tpcbeam - dw_tpcbeam, dz_fangsbeamhb);
800  G4VSolid* s_fangsbeamhbpos = new G4SubtractionSolid("s_fangsbeamhbpos", s_fangsbeamhb_a, s_fangsbeamhb_b, 0, G4ThreeVector(0,
801  dw_tpcbeam, 0));
802  G4VSolid* s_fangsbeamhbneg = new G4SubtractionSolid("s_fangsbeamhbneg", s_fangsbeamhb_a, s_fangsbeamhb_b, 0, G4ThreeVector(0,
803  -dw_tpcbeam, 0));
804  G4VSolid* s_fangsbeamhb = new G4UnionSolid("s_fangsbeamhb", s_fangsbeamhbpos, s_fangsbeamhbneg, 0, G4ThreeVector(0, -2.*dy_tpcbeam,
805  0));
806  G4LogicalVolume* l_fangsbeamhb = new G4LogicalVolume(s_fangsbeamhb, geometry::Materials::get("FG_Epoxy"), "l_fangsbeamhb", 0, 0);
807  double xfangsbeamhb = activeParams.getLength("xfangsBeamhb") * CLHEP::cm;
808  double yfangsbeamhb = activeParams.getLength("yfangsBeamhb") * CLHEP::cm;
809  double zfangsbeamhb = activeParams.getLength("zfangsBeamhb") * CLHEP::cm;
810  TransForm = G4Translate3D(xfangsbeamhb, yfangsbeamhb, zfangsbeamhb) /* G4RotateX3D(90.*CLHEP::deg)*/ * G4RotateY3D(90.*CLHEP::deg);
811  new G4PVPlacement(TransForm, l_fangsbeamhb, "p_fangsbeamhb", &topVolume, false, 1);
812  }
813  }
814  } // ph1sustr namespace
816 } // Belle2 namespace
Belle2::GearDir
GearDir is the basic class used for accessing the parameter store.
Definition: GearDir.h:31
Belle2::gearbox::Interface::getArray
std::vector< double > getArray(const std::string &path) const noexcept(false)
Get the parameter path as a list of double values converted to the standard unit.
Definition: Interface.cc:123
Belle2::gearbox::Interface::getLength
double getLength(const std::string &path="") const noexcept(false)
Get the parameter path as a double converted to the standard length unit.
Definition: Interface.h:259
Belle2::geometry::Materials::get
static G4Material * get(const std::string &name)
Find given material.
Definition: Materials.h:63
Belle2::ph1sustr::Ph1sustrCreator::create
virtual void create(const GearDir &content, G4LogicalVolume &topVolume, geometry::GeometryTypes type)
Creation of the detector geometry from Gearbox (XML).
Definition: Ph1sustrCreator.cc:60
Belle2::ph1sustr::Ph1sustrCreator::~Ph1sustrCreator
virtual ~Ph1sustrCreator()
Destructor.
Definition: Ph1sustrCreator.cc:55
Belle2::ph1sustr::Ph1sustrCreator::m_sensitive
SensitiveDetector * m_sensitive
SensitiveDetector phase 1 support structure.
Definition: Ph1sustrCreator.h:46
Belle2::ph1sustr::SensitiveDetector
Sensitive Detector implementation of the PH1SUSTR detector.
Definition: SensitiveDetector.h:22
Belle2::geometry::GeometryTypes
GeometryTypes
Flag indiciating the type of geometry to be used.
Definition: GeometryManager.h:36
Belle2::ph1sustr::Ph1sustrFactory
geometry::CreatorFactory< Ph1sustrCreator > Ph1sustrFactory("PH1SUSTRCreator")
Creator creates the phase 1 support structure geometry.
Belle2
Abstract base class for different kinds of events.
Definition: MillepedeAlgorithm.h:17
Belle2::geometry::CreatorFactory
Very simple class to provide an easy way to register creators with the CreatorManager.
Definition: CreatorFactory.h:31