82 SetTorusDimensions(
r, rmin, rmax, phi1, dphi);
83 if ((fRmin < 0) || (fRmax < 0))
84 SetShapeBit(kGeoRunTimeShape);
95 SetTorusDimensions(
r, rmin, rmax, phi1, dphi);
96 if ((fRmin < 0) || (fRmax < 0))
97 SetShapeBit(kGeoRunTimeShape);
112 SetDimensions(param);
113 if (fRmin < 0 || fRmax < 0)
114 SetShapeBit(kGeoRunTimeShape);
200 if ((
fRmax - daxis) > 1E-5) {
212 for (
Int_t i = 0; i < 3; i++) {
213 norm[i] = point[i] - r0[i];
214 normsq += norm[i] * norm[i];
221 if (dir[0] * norm[0] + dir[1] * norm[1] + dir[2] * norm[2] < 0) {
246 Double_t radsq = (rxy -
fR) * (rxy -
fR) + point[2] * point[2];
263 else if (
fDphi < 360)
274 for (
Int_t i = 0; i < 3; i++)
275 p[i] =
pt[i] + t * dir[i];
286 for (
Int_t i = 0; i < 3; i++)
287 p[i] =
pt[i] + t * dir[i];
290 return ((
p[2] * dir[2] -
fR *
TMath::Sqrt(dir[0] * dir[0] + dir[1] * dir[1])) /
295 Double_t dd = (
p[0] * dir[0] +
p[1] * dir[1] +
p[2] * dir[2] - (
p[0] * dir[0] +
p[1] * dir[1]) *
fR / rxy) /
d;
305 for (
Int_t i = 0; i < 3; i++)
306 p[i] =
pt[i] + t * dir[i];
314 (
p[0] * dir[0] +
p[1] * dir[1] +
p[2] * dir[2] - (
p[0] * dir[0] +
p[1] * dir[1]) *
fR / rxy) / daxis;
315 Double_t dddaxis = 1 - ddaxis * ddaxis - (1 - dir[2] * dir[2]) *
fR / rxy +
316 fR * (
p[0] * dir[0] +
p[1] * dir[1]) * (
p[0] * dir[0] +
p[1] * dir[1]) / (rxy * rxy * rxy);
327 if (iact < 3 && safe) {
331 if ((iact == 1) && (step <= *safe))
356 TGeoTubeSeg::DistFromInsideS(point, dir,
fR -
fRmax,
fR +
fRmax,
fRmax,
c1,
s1,
c2, s2, cm, sm, cdfi);
358 if (daxis >=
fRmin + 1.E-8 && daxis <=
fRmax - 1.E-8)
370 if (iact < 3 && safe) {
374 if ((iact == 1) && (step <= *safe))
384 Double_t c1 = 0,
s1 = 0,
c2 = 0, s2 = 0, cm = 0, sm = 0, cdfi = 0;
410 fio = 0.5 * (phi1 + phi2);
418 rxy2 = point[0] * point[0] + point[1] * point[1];
420 if (!hasphi || inphi)
443 daxis =
Daxis(point, dir, snext);
447 for (i = 0; i < 3; i++)
448 pt[i] = point[i] + snext * dir[i];
454 if (daxis <
fRmin + 1.E-8) {
459 for (i = 0; i < 3; i++)
460 pt[i] += 0.1 * eps * dir[i];
467 dring =
TGeoTubeSeg::DistFromInsideS(
pt, dir,
fR -
fRmin,
fR +
fRmin,
fRmin,
c1,
s1,
c2, s2, cm, sm, cdfi);
473 snext += dring + eps;
474 for (i = 0; i < 3; i++)
475 pt[i] = point[i] + snext * dir[i];
484 for (i = 0; i < 3; i++)
485 pt[i] += 0.1 * eps * dir[i];
491 s1,
c2, s2, cm, sm, cdfi);
499 snext += dring + eps;
500 for (i = 0; i < 3; i++)
501 pt[i] = point[i] + snext * dir[i];
523 case 2:
return "PHI";
525 default:
return "UNDEFINED";
548 case 3: dx = 0;
return dx;
572 Error(
"GetMakeRuntimeShape",
"parametrized toruses not supported");
581 printf(
"*** Shape %s: TGeoTorus ***\n",
GetName());
582 printf(
" R = %11.5f\n",
fR);
583 printf(
" Rmin = %11.5f\n",
fRmin);
584 printf(
" Rmax = %11.5f\n",
fRmax);
585 printf(
" Phi1 = %11.5f\n",
fPhi1);
586 printf(
" Dphi = %11.5f\n",
fDphi);
587 printf(
" Bounding box:\n");
606 Int_t nbSegs = (2 *
n - 1) * (
n - 1);
607 Int_t nbPols = (
n - 1) * (
n - 1);
609 nbSegs += (2 *
n - 1) * (
n - 1);
610 nbPols += (
n - 1) * (
n - 1);
613 nbSegs += 2 * (
n - 1);
614 nbPols += 2 * (
n - 1);
635 Int_t indx, indp, startcap = 0;
648 for (i = 0; i <
n; i++) {
649 for (j = 0; j <
n - 1; j++) {
650 buff.
fSegs[indx + (i * (
n - 1) + j) * 3] =
c;
651 buff.
fSegs[indx + (i * (
n - 1) + j) * 3 + 1] = i * (
n - 1) + j;
652 buff.
fSegs[indx + (i * (
n - 1) + j) * 3 + 2] = i * (
n - 1) + ((j + 1) % (
n - 1));
655 indx += 3 *
n * (
n - 1);
658 for (i = 0; i <
n - 1; i++) {
659 for (j = 0; j <
n - 1; j++) {
660 buff.
fSegs[indx + (i * (
n - 1) + j) * 3] =
c;
661 buff.
fSegs[indx + (i * (
n - 1) + j) * 3 + 1] = i * (
n - 1) + j;
662 buff.
fSegs[indx + (i * (
n - 1) + j) * 3 + 2] = (i + 1) * (
n - 1) + j;
665 indx += 3 * (
n - 1) * (
n - 1);
666 startcap = (2 *
n - 1) * (
n - 1);
671 for (i = 0; i <
n; i++) {
672 for (j = 0; j <
n - 1; j++) {
673 buff.
fSegs[indx + (i * (
n - 1) + j) * 3] =
c;
674 buff.
fSegs[indx + (i * (
n - 1) + j) * 3 + 1] = indp + i * (
n - 1) + j;
675 buff.
fSegs[indx + (i * (
n - 1) + j) * 3 + 2] = indp + i * (
n - 1) + ((j + 1) % (
n - 1));
678 indx += 3 *
n * (
n - 1);
681 for (i = 0; i <
n - 1; i++) {
682 for (j = 0; j <
n - 1; j++) {
683 buff.
fSegs[indx + (i * (
n - 1) + j) * 3] =
c;
684 buff.
fSegs[indx + (i * (
n - 1) + j) * 3 + 1] = indp + i * (
n - 1) + j;
685 buff.
fSegs[indx + (i * (
n - 1) + j) * 3 + 2] = indp + (i + 1) * (
n - 1) + j;
688 indx += 3 * (
n - 1) * (
n - 1);
689 startcap = (4 *
n - 2) * (
n - 1);
696 for (j = 0; j <
n - 1; j++) {
697 buff.
fSegs[indx + j * 3] =
c + 1;
698 buff.
fSegs[indx + j * 3 + 1] = (
n - 1) * i + j;
699 buff.
fSegs[indx + j * 3 + 2] = indp + (
n - 1) * i + j;
703 for (j = 0; j <
n - 1; j++) {
704 buff.
fSegs[indx + j * 3] =
c + 1;
705 buff.
fSegs[indx + j * 3 + 1] = (
n - 1) * i + j;
706 buff.
fSegs[indx + j * 3 + 2] = indp + (
n - 1) * i + j;
711 for (j = 0; j <
n - 1; j++) {
712 buff.
fSegs[indx + j * 3] =
c + 1;
713 buff.
fSegs[indx + j * 3 + 1] = (
n - 1) * i + j;
714 buff.
fSegs[indx + j * 3 + 2] =
n * (
n - 1);
718 for (j = 0; j <
n - 1; j++) {
719 buff.
fSegs[indx + j * 3] =
c + 1;
720 buff.
fSegs[indx + j * 3 + 1] = (
n - 1) * i + j;
721 buff.
fSegs[indx + j * 3 + 2] =
n * (
n - 1) + 1;
732 for (i = 0; i <
n - 1; i++) {
733 for (j = 0; j <
n - 1; j++) {
735 buff.
fPols[indx++] = 4;
736 buff.
fPols[indx++] =
n * (
n - 1) + (
n - 1) * i + ((j + 1) % (
n - 1));
737 buff.
fPols[indx++] = (
n - 1) * (i + 1) + j;
738 buff.
fPols[indx++] =
n * (
n - 1) + (
n - 1) * i + j;
739 buff.
fPols[indx++] = (
n - 1) * i + j;
743 indp = (2 *
n - 1) * (
n - 1);
746 for (i = 0; i <
n - 1; i++) {
747 for (j = 0; j <
n - 1; j++) {
749 buff.
fPols[indx++] = 4;
750 buff.
fPols[indx++] = indp +
n * (
n - 1) + (
n - 1) * i + j;
751 buff.
fPols[indx++] = indp + (
n - 1) * (i + 1) + j;
752 buff.
fPols[indx++] = indp +
n * (
n - 1) + (
n - 1) * i + ((j + 1) % (
n - 1));
753 buff.
fPols[indx++] = indp + (
n - 1) * i + j;
761 for (j = 0; j <
n - 1; j++) {
762 buff.
fPols[indx++] =
c + 1;
764 buff.
fPols[indx++] = j;
765 buff.
fPols[indx++] = startcap + j;
767 buff.
fPols[indx++] = indp + j;
768 buff.
fPols[indx++] = startcap + ((j + 1) % (
n - 1));
772 for (j = 0; j <
n - 1; j++) {
773 buff.
fPols[indx++] =
c + 1;
775 buff.
fPols[indx++] = (
n - 1) * i + j;
776 buff.
fPols[indx++] = startcap + (
n - 1) + ((j + 1) % (
n - 1));
778 buff.
fPols[indx++] = indp + (
n - 1) * i + j;
779 buff.
fPols[indx++] = startcap + (
n - 1) + j;
794 saf[0] = rad -
fRmin;
795 saf[1] =
fRmax - rad;
799 for (i = 0; i < 2; i++)
809 for (i = 0; i < 2; i++)
823 out <<
" r = " <<
fR <<
";" << std::endl;
824 out <<
" rmin = " <<
fRmin <<
";" << std::endl;
825 out <<
" rmax = " <<
fRmax <<
";" << std::endl;
826 out <<
" phi1 = " <<
fPhi1 <<
";" << std::endl;
827 out <<
" dphi = " <<
fDphi <<
";" << std::endl;
871 for (i = 0; i <
n; i++) {
875 for (j = 0; j <
n - 1; j++) {
887 for (i = 0; i <
n; i++) {
891 for (j = 0; j <
n - 1; j++) {
930 for (i = 0; i <
n; i++) {
934 for (j = 0; j <
n - 1; j++) {
947 for (i = 0; i <
n; i++) {
951 for (j = 0; j <
n - 1; j++) {
984 else if (
fDphi < 360.)
1012 t = (-3 *
q * sq3 + delta) / (6 * sq3);
1013 u = (3 *
q * sq3 + delta) / (6 * sq3);
1019 u = delta / (6 * sq3);
1024 t =
x[0] *
x[0] +
a *
x[0] +
b;
1026 delta = u * u - 4. * t;
1030 x[1] = 0.5 * (-u - delta);
1031 x[2] = 0.5 * (-u + delta);
1054 delta =
e *
e - 4. *
g;
1058 h = 0.5 * (-
e - delta);
1061 x[ireal++] = -
h - 0.25 *
a;
1062 x[ireal++] =
h - 0.25 *
a;
1064 h = 0.5 * (-
e + delta);
1067 x[ireal++] = -
h - 0.25 *
a;
1068 x[ireal++] =
h - 0.25 *
a;
1072 for (i = 0; i < ireal; i++)
1074 memcpy(
x, xx, ireal *
sizeof(
Double_t));
1080 x[ireal++] = -0.25 *
a;
1082 for (i = 0; i < ind[0]; i++)
1083 x[ireal++] = xx[i] - 0.25 *
a;
1086 for (i = 0; i < ireal; i++)
1088 memcpy(
x, xx, ireal *
sizeof(
Double_t));
1100 for (i = 0; i < 3; i++) {
1111 delta =
h *
h - 4. * j;
1114 x[ireal++] = 0.5 * (-
h - delta) - 0.25 *
a;
1115 x[ireal++] = 0.5 * (-
h + delta) - 0.25 *
a;
1117 delta =
h *
h - 4. *
g / j;
1120 x[ireal++] = 0.5 * (
h - delta) - 0.25 *
a;
1121 x[ireal++] = 0.5 * (
h + delta) - 0.25 *
a;
1125 for (i = 0; i < ireal; i++)
1127 memcpy(
x, xx, ireal *
sizeof(
Double_t));
1142 Double_t rdotn =
pt[0] * dir[0] +
pt[1] * dir[1] +
pt[2] * dir[2];
1145 Double_t b = 2. * (r0sq + 2. * rdotn * rdotn - rsumsq + 2. *
fR *
fR * dir[2] * dir[2]);
1146 Double_t c = 4. * (r0sq * rdotn - rsumsq * rdotn + 2. *
fR *
fR *
pt[2] * dir[2]);
1154 Double_t b0 = (
pt[0] * dir[0] +
pt[1] * dir[1]) / (dir[0] * dir[0] + dir[1] * dir[1]);
1155 Double_t c0 = (
pt[0] *
pt[0] + (
pt[1] - r0) * (
pt[1] + r0)) / (dir[0] * dir[0] + dir[1] * dir[1]);
1166 c0 = (
pt[0] *
pt[0] + (
pt[1] - r0) * (
pt[1] + r0)) / (dir[0] * dir[0] + dir[1] * dir[1]);
1167 delta = b0 * b0 -
c0;
1180 for (
Int_t j = 0; j < nsol; j++)
1192 for (
Int_t i = 0; i < nsol; i++) {
1199 for (
Int_t ipt = 0; ipt < 3; ipt++)
1200 norm[ipt] =
pt[ipt] +
x[i] * dir[ipt] - r0[ipt];
1201 ndotd = norm[0] * dir[0] + norm[1] * dir[1] + norm[2] * dir[2];
1211 Double_t delta = s * s * s * s +
a * s * s * s +
b * s * s +
c * s +
d;
1212 Double_t eps0 = -delta / (4. * s * s * s + 3. *
a * s * s + 2. *
b * s +
c);
1219 delta = s * s * s * s +
a * s * s * s +
b * s * s +
c * s +
d;
1220 eps = -delta / (4. * s * s * s + 3. *
a * s * s + 2. *
b * s +
c);
1238 nvert =
n * (
n - 1);
1245 nsegs = (2 *
n - 1) * (
n - 1);
1246 npols = (
n - 1) * (
n - 1);
1248 nsegs += (2 *
n - 1) * (
n - 1);
1249 npols += (
n - 1) * (
n - 1);
1252 nsegs += 2 * (
n - 1);
1253 npols += 2 * (
n - 1);
1276 Int_t nbSegs = (2 *
n - 1) * (
n - 1);
1277 Int_t nbPols = (
n - 1) * (
n - 1);
1279 nbSegs += (2 *
n - 1) * (
n - 1);
1280 nbPols += (
n - 1) * (
n - 1);
1283 nbSegs += 2 * (
n - 1);
1284 nbPols += 2 * (
n - 1);
1287 if (buffer.SetRawSizes(nbPnts, 3 * nbPnts, nbSegs, 3 * nbSegs, nbPols, 6 * nbPols)) {
1295 if (!buffer.fLocalFrame) {
1313 for (
Int_t i = 0; i < vecsize; i++)
1324 for (
Int_t i = 0; i < vecsize; i++)
1334 for (
Int_t i = 0; i < vecsize; i++)
1344 for (
Int_t i = 0; i < vecsize; i++)
1355 for (
Int_t i = 0; i < vecsize; i++)
winID h TVirtualViewer3D TVirtualGLPainter p
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t Float_t Float_t Int_t Int_t UInt_t UInt_t Rectangle_t Int_t Int_t Window_t TString Int_t GCValues_t GetPrimarySelectionOwner GetDisplay GetScreen GetColormap GetNativeEvent const char const char dpyName wid window const char font_name cursor keysym reg const char only_if_exist regb h Point_t np
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t r
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t points
R__EXTERN TGeoManager * gGeoManager
Generic 3D primitive description class.
void FillBuffer3D(TBuffer3D &buffer, Int_t reqSections, Bool_t localFrame) const override
Fill the supplied buffer, with sections in desired frame See TBuffer3D.h for explanation of sections,...
Double_t DistFromOutside(const Double_t *point, const Double_t *dir, Int_t iact=1, Double_t step=TGeoShape::Big(), Double_t *safe=nullptr) const override
void InspectShape() const override
Int_t GetNsegments() const
Get number of segments approximating circles.
Geometrical transformation package.
Base abstract class for all shapes.
Int_t GetBasicColor() const
Get the basic color (0-7).
void TransformPoints(Double_t *points, UInt_t NbPoints) const
Tranform a set of points (LocalToMaster)
void SetShapeBit(UInt_t f, Bool_t set)
Equivalent of TObject::SetBit.
static Double_t SafetyPhi(const Double_t *point, Bool_t in, Double_t phi1, Double_t phi2)
Static method to compute safety w.r.t a phi corner defined by cosines/sines of the angles phi1,...
static Bool_t IsSameWithinTolerance(Double_t a, Double_t b)
Check if two numbers differ with less than a tolerance.
const char * GetPointerName() const
Provide a pointer name containing uid.
Int_t ShapeDistancetoPrimitive(Int_t numpoints, Int_t px, Int_t py) const
Returns distance to shape primitive mesh.
static void NormalPhi(const Double_t *point, const Double_t *dir, Double_t *norm, Double_t c1, Double_t s1, Double_t c2, Double_t s2)
Static method to compute normal to phi planes.
const char * GetName() const override
Get the shape name.
static Double_t Tolerance()
Bool_t TestShapeBit(UInt_t f) const
Int_t SolveQuartic(Double_t a, Double_t b, Double_t c, Double_t d, Double_t *x) const
void Sizeof3D() const override
Double_t ToBoundary(const Double_t *pt, const Double_t *dir, Double_t r, Bool_t in) const
Double_t GetAxisRange(Int_t iaxis, Double_t &xlo, Double_t &xhi) const override
Double_t DDDaxis(const Double_t *pt, const Double_t *dir, Double_t t) const
Double_t DistFromInside(const Double_t *point, const Double_t *dir, Int_t iact=1, Double_t step=TGeoShape::Big(), Double_t *safe=nullptr) const override
const char * GetAxisName(Int_t iaxis) const override
void SavePrimitive(std::ostream &out, Option_t *option="") override
Save a primitive as a C++ statement(s) on output stream "out".
void DistFromInside_v(const Double_t *points, const Double_t *dirs, Double_t *dists, Int_t vecsize, Double_t *step) const override
Double_t DistFromOutside(const Double_t *point, const Double_t *dir, Int_t iact=1, Double_t step=TGeoShape::Big(), Double_t *safe=nullptr) const override
Double_t Daxis(const Double_t *pt, const Double_t *dir, Double_t t) const
void ComputeNormal(const Double_t *point, const Double_t *dir, Double_t *norm) override
void SetTorusDimensions(Double_t r, Double_t rmin, Double_t rmax, Double_t phi1, Double_t dphi)
Int_t DistancetoPrimitive(Int_t px, Int_t py) override
Computes distance from point (px,py) to the object.
TBuffer3D * MakeBuffer3D() const override
void SetPoints(Double_t *points) const override
Int_t GetNmeshVertices() const override
void SetSegsAndPols(TBuffer3D &buff) const override
void Safety_v(const Double_t *points, const Bool_t *inside, Double_t *safe, Int_t vecsize) const override
TGeoShape * GetMakeRuntimeShape(TGeoShape *mother, TGeoMatrix *mat) const override
void GetBoundingCylinder(Double_t *param) const override
void GetMeshNumbers(Int_t &nvert, Int_t &nsegs, Int_t &npols) const override
Double_t DDaxis(const Double_t *pt, const Double_t *dir, Double_t t) const
Bool_t Contains(const Double_t *point) const override
void Contains_v(const Double_t *points, Bool_t *inside, Int_t vecsize) const override
void InspectShape() const override
TGeoVolume * Divide(TGeoVolume *voldiv, const char *divname, Int_t iaxis, Int_t ndiv, Double_t start, Double_t step) override
const TBuffer3D & GetBuffer3D(Int_t reqSections, Bool_t localFrame) const override
Stub implementation to avoid forcing implementation at this stage.
void SetDimensions(Double_t *param) override
void ComputeBBox() override
void ComputeNormal_v(const Double_t *points, const Double_t *dirs, Double_t *norms, Int_t vecsize) override
Int_t SolveCubic(Double_t a, Double_t b, Double_t c, Double_t *x) const
void DistFromOutside_v(const Double_t *points, const Double_t *dirs, Double_t *dists, Int_t vecsize, Double_t *step) const override
Double_t Safety(const Double_t *point, Bool_t in=kTRUE) const override
Double_t Capacity() const override
static Double_t DistFromInsideS(const Double_t *point, const Double_t *dir, Double_t rmin, Double_t rmax, Double_t dz, Double_t c1, Double_t s1, Double_t c2, Double_t s2, Double_t cm, Double_t sm, Double_t cdfi)
static Double_t DistFromOutsideS(const Double_t *point, const Double_t *dir, Double_t rmin, Double_t rmax, Double_t dz, Double_t c1, Double_t s1, Double_t c2, Double_t s2, Double_t cm, Double_t sm, Double_t cdfi)
static Double_t DistFromInsideS(const Double_t *point, const Double_t *dir, Double_t rmin, Double_t rmax, Double_t dz)
static Double_t DistFromOutsideS(const Double_t *point, const Double_t *dir, Double_t rmin, Double_t rmax, Double_t dz)
TGeoVolume, TGeoVolumeMulti, TGeoVolumeAssembly are the volume classes.
R__ALWAYS_INLINE Bool_t TestBit(UInt_t f) const
virtual const char * ClassName() const
Returns name of class to which the object belongs.
void SetBit(UInt_t f, Bool_t set)
Set or unset the user status bits as specified in f.
virtual void Error(const char *method, const char *msgfmt,...) const
Issue error message.
Long64_t LocMin(Long64_t n, const T *a)
Returns index of array with the minimum element.
Short_t Max(Short_t a, Short_t b)
Returns the largest of a and b.
T1 Sign(T1 a, T2 b)
Returns a value with the magnitude of a and the sign of b.
Double_t ATan2(Double_t y, Double_t x)
Returns the principal value of the arc tangent of y/x, expressed in radians.
Long64_t LocMax(Long64_t n, const T *a)
Returns index of array with the maximum element.
constexpr Double_t DegToRad()
Conversion from degree to radian: .
Double_t Sqrt(Double_t x)
Returns the square root of x.
LongDouble_t Power(LongDouble_t x, LongDouble_t y)
Returns x raised to the power y.
Short_t Min(Short_t a, Short_t b)
Returns the smallest of a and b.
Double_t Cos(Double_t)
Returns the cosine of an angle of x radians.
Double_t Sin(Double_t)
Returns the sine of an angle of x radians.
void Sort(Index n, const Element *a, Index *index, Bool_t down=kTRUE)
Sort the n elements of the array a of generic templated type Element.
constexpr Double_t RadToDeg()
Conversion from radian to degree: .
Short_t Abs(Short_t d)
Returns the absolute value of parameter Short_t d.