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
|
#include "../benchmark/ArgParser.h"
#include "ShapeTester.h"
#include "VecGeom/volumes/Polycone.h"
typedef vecgeom::SimplePolycone Poly_t;
using vecgeom::Precision;
int main(int argc, char *argv[])
{
OPTION_INT(npoints, 10000);
OPTION_BOOL(debug, false);
OPTION_BOOL(stat, false);
using namespace vecgeom;
Precision Z_ValP[15];
Z_ValP[0] = -1520.;
Z_ValP[1] = -804.;
Z_ValP[2] = -804.;
Z_ValP[3] = -515.345;
Z_ValP[4] = -515.345;
Z_ValP[5] = -177.;
Z_ValP[6] = -177.;
Z_ValP[7] = 149.561;
Z_ValP[8] = 149.561;
Z_ValP[9] = 575.;
Z_ValP[10] = 575.;
Z_ValP[11] = 982.812;
Z_ValP[12] = 982.812;
Z_ValP[13] = 1166.7;
Z_ValP[14] = 1524.;
Precision R_MinP[15];
R_MinP[0] = 1238.;
R_MinP[1] = 1238.;
R_MinP[2] = 1238.;
R_MinP[3] = 1238.;
R_MinP[4] = 1238.;
R_MinP[5] = 1238.;
R_MinP[6] = 1238.;
R_MinP[7] = 1238.;
R_MinP[8] = 1238.;
R_MinP[9] = 1238.;
R_MinP[10] = 1238.;
R_MinP[11] = 1238.;
R_MinP[12] = 1238.;
R_MinP[13] = 1238.;
R_MinP[14] = 1455.22;
Precision R_MaxP[15];
R_MaxP[0] = 1555.01;
R_MaxP[1] = 1555.01;
R_MaxP[2] = 1538.05;
R_MaxP[3] = 1538.05;
R_MaxP[4] = 1523.26;
R_MaxP[5] = 1523.26;
R_MaxP[6] = 1506.24;
R_MaxP[7] = 1506.24;
R_MaxP[8] = 1488.52;
R_MaxP[9] = 1488.52;
R_MaxP[10] = 1471.28;
R_MaxP[11] = 1471.28;
R_MaxP[12] = 1455.22;
R_MaxP[13] = 1455.22;
R_MaxP[14] = 1455.22;
int Nz = 15;
auto poly2 = new Poly_t("Test", 0., /* initial phi starting angle */
kTwoPi, /* total phi angle */
Nz, /* number corners in r,z space */
Z_ValP, R_MinP, R_MaxP); /* r coordinate of these corners */
ShapeTester<vecgeom::VPlacedVolume> tester;
tester.setDebug(debug);
tester.setStat(stat);
tester.SetMaxPoints(npoints);
tester.SetTestBoundaryErrors(false);
#ifdef VECGEOM_FLOAT_PRECISION
tester.SetSolidTolerance(1.e-4);
#else
tester.SetSolidTolerance(2.e-7);
#endif
int errCode = tester.Run(poly2);
std::cout << "Final Error count for Shape *** " << poly2->GetName() << "*** = " << errCode << "\n";
std::cout << "=========================================================" << std::endl;
if (poly2) delete poly2;
return 0;
}
|
