1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
|
#include "VecGeom/volumes/LogicalVolume.h"
#include "VecGeom/volumes/Box.h"
#include "VecGeom/volumes/Trapezoid.h"
#include "VecGeomTest/Benchmarker.h"
#include "VecGeom/management/GeoManager.h"
#include "ArgParser.h"
using namespace vecgeom;
Precision pmax(Precision p1, Precision p2)
{
if (p1 > p2) return p1;
return p2;
}
int main(int argc, char *argv[])
{
OPTION_INT(npoints, 1024);
OPTION_INT(nrep, 1024);
OPTION_DOUBLE(dz, 15.);
OPTION_DOUBLE(p1x, -2);
OPTION_DOUBLE(p2x, 2);
OPTION_DOUBLE(p3x, -3);
OPTION_DOUBLE(p4x, 3);
OPTION_DOUBLE(p5x, -4);
OPTION_DOUBLE(p6x, 4);
OPTION_DOUBLE(p7x, -6);
OPTION_DOUBLE(p8x, 6);
OPTION_DOUBLE(p1y, -5);
OPTION_DOUBLE(p2y, -5);
OPTION_DOUBLE(p3y, 5);
OPTION_DOUBLE(p4y, 5);
OPTION_DOUBLE(p5y, -10);
OPTION_DOUBLE(p6y, -10);
OPTION_DOUBLE(p7y, 10);
OPTION_DOUBLE(p8y, 10);
// UnplacedBox worldUnplaced = UnplacedBox(pmax(p1x, pmax(p2x, pmax(p3x, p4x)))*4, pmax(p1y, pmax(p2y, pmax(p3y,
// p4y)))*4, dz*4);
UnplacedBox worldUnplaced = UnplacedBox(20., 20., 20.);
// validate construtor for input corner points -- add an xy-offset for non-zero theta,phi
TrapCorners xyz;
Precision xoffset = 9;
Precision yoffset = -6;
// define corner points
// convention: p0(---); p1(+--); p2(-+-); p3(++-); p4(--+); p5(+-+); p6(-++); p7(+++)
xyz[0] = Vector3D<Precision>(p1x + xoffset, p1y + yoffset, -dz);
xyz[1] = Vector3D<Precision>(p2x + xoffset, p2y + yoffset, -dz);
xyz[2] = Vector3D<Precision>(p3x + xoffset, p3y + yoffset, -dz);
xyz[3] = Vector3D<Precision>(p4x + xoffset, p4y + yoffset, -dz);
xyz[4] = Vector3D<Precision>(p5x - xoffset, p5y - yoffset, dz);
xyz[5] = Vector3D<Precision>(p6x - xoffset, p6y - yoffset, dz);
xyz[6] = Vector3D<Precision>(p7x - xoffset, p7y - yoffset, dz);
xyz[7] = Vector3D<Precision>(p8x - xoffset, p8y - yoffset, dz);
// xyz[0] = Vector3D<Precision>( -2+xoffset, -5+yoffset, -15 );
// xyz[1] = Vector3D<Precision>( 2+xoffset, -5+yoffset, -15 );
// xyz[2] = Vector3D<Precision>( -3+xoffset, 5+yoffset, -15 );
// xyz[3] = Vector3D<Precision>( 3+xoffset, 5+yoffset, -15 );
// xyz[4] = Vector3D<Precision>( -4-xoffset,-10-yoffset, 15 );
// xyz[5] = Vector3D<Precision>( 4-xoffset,-10-yoffset, 15 );
// xyz[6] = Vector3D<Precision>( -6-xoffset, 10-yoffset, 15 );
// xyz[7] = Vector3D<Precision>( 6-xoffset, 10-yoffset, 15 );
//
// create trapezoid
UnplacedTrapezoid trapUnplaced(xyz);
// UnplacedTrapezoid trapUnplaced2(1,0,0, 1,1,1,0, 1,1,1,0);
LogicalVolume world("world", &worldUnplaced);
LogicalVolume trap("trap", &trapUnplaced);
Transformation3D transf(5, 5, 5);
world.PlaceDaughter(&trap, &transf);
// world.PlaceDaughter(&trap, &Transformation3D::kIdentity);
VPlacedVolume *worldPlaced = world.Place();
GeoManager::Instance().SetWorldAndClose(worldPlaced);
Benchmarker tester(GeoManager::Instance().GetWorld());
tester.SetVerbosity(3);
tester.SetPointCount(npoints);
tester.SetRepetitions(nrep);
return tester.RunBenchmark();
// cleanup
// delete transf;
}
|
