67 Fatal(
"ComputeBBox",
"Assembly shape %s without volume",
GetName());
81 for (
Int_t i=0; i<nd; i++) {
86 box->SetBoxPoints(vert);
87 for (
Int_t ipt=0; ipt<8; ipt++) {
93 if (
pt[2]<zmin) zmin=
pt[2];
94 if (
pt[2]>zmax) zmax=
pt[2];
101 fDZ = 0.5*(zmax-zmin);
134 box->SetBoxPoints(vert);
135 for (
Int_t ipt=0; ipt<8; ipt++) {
141 if (
pt[2]<zmin) zmin=
pt[2];
142 if (
pt[2]>zmax) zmax=
pt[2];
148 fDZ = 0.5*(zmax-zmin);
186 Int_t *check_list = 0;
198 for (
id=0;
id<ncheck;
id++) {
213 for (
id=0;
id<nd;
id++) {
239 Info(
"DistFromInside",
"Cannot compute distance from inside the assembly (but from a component)");
258 if (iact<3 && safe) {
293 for (i=0; i<3; i++)
pt[i] +=
snext*dir[i];
310 if (nd<5 || !voxels) {
311 for (i=0; i<nd; i++) {
317 if (idebug>4) printf(
"%s[%d] distance to %s ...\n", sindent.
Data(),
indent, node->
GetName());
323 printf(
"%s[%d] %s -> from local=(%19.16f, %19.16f, %19.16f, %19.16f, %19.16f, %19.16f)\n",
325 printf(
"%s[%d] -> (l)to: %s shape %s snext=%g\n", sindent.
Data(),
indent, node->
GetName(),
357 for (i=0; i<ncheck; i++) {
362 if (idebug>4) printf(
"%s[%d] distance to %s ...\n", sindent.
Data(),
indent, node->
GetName());
368 printf(
"%s[%d] %s -> from local=(%19.16f, %19.16f, %19.16f, %19.16f, %19.16f, %19.16f)\n",
370 printf(
"%s[%d] -> to: %s shape %s snext=%g\n", sindent.
Data(),
indent, node->
GetName(),
402 Error(
"Divide",
"Assemblies cannot be divided. Division volume %s not created", divname);
412 Error(
"GetMakeRuntimeShape",
"Assemblies cannot be parametrized.");
421 printf(
"*** Shape %s: TGeoShapeAssembly ***\n",
GetName());
423 printf(
" Bounding box:\n");
433 Error(
"SetSegsAndPols",
"Drawing functions should not be called for assemblies, but rather for their content");
449 memcpy(loc, point, 3*
sizeof(
Double_t));
467 if (voxels) boxes = voxels->
GetBoxes();
469 for (
Int_t id=0;
id<nd;
id++) {
474 if (dxyz0 > safety)
continue;
476 if (dxyz1 > safety)
continue;
478 if (dxyz2 > safety)
continue;
479 if (dxyz0>0) dxyz+=dxyz0*dxyz0;
480 if (dxyz1>0) dxyz+=dxyz1*dxyz1;
481 if (dxyz2>0) dxyz+=dxyz2*dxyz2;
482 if (dxyz >= safety*safety)
continue;
486 if (safe<=0.0)
return 0.0;
487 if (safe<safety) safety = safe;
504 Error(
"SetPoints",
"Drawing functions should not be called for assemblies, but rather for their content");
512 Error(
"SetPoints",
"Drawing functions should not be called for assemblies, but rather for their content");
static void indent(ostringstream &buf, int indent_level)
R__EXTERN TGeoManager * gGeoManager
Generic 3D primitive description class.
virtual Double_t DistFromOutside(const Double_t *point, const Double_t *dir, Int_t iact=1, Double_t step=TGeoShape::Big(), Double_t *safe=0) const
Compute distance from outside point to surface of the box.
virtual void InspectShape() const
Prints shape parameters.
virtual Bool_t Contains(const Double_t *point) const
Test if point is inside this shape.
TGeoNavigator * GetCurrentNavigator() const
Returns current navigator for the calling thread.
static Int_t GetVerboseLevel()
Set verbosity level (static function).
Geometrical transformation package.
virtual void MasterToLocal(const Double_t *master, Double_t *local) const
convert a point by multiplying its column vector (x, y, z, 1) to matrix
Class providing navigation API for TGeo geometries.
TGeoNodeCache * GetCache() const
TGeoStateInfo * GetInfo()
Get next state info pointer.
void ReleaseInfo()
Release last used state info pointer.
A node represent a volume positioned inside another.They store links to both volumes and to the TGeoM...
TGeoVolume * GetVolume() const
virtual TGeoMatrix * GetMatrix() const =0
virtual void LocalToMasterVect(const Double_t *local, Double_t *master) const
Convert a vector from local reference system to mother reference.
virtual void LocalToMaster(const Double_t *local, Double_t *master) const
Convert the point coordinates from local reference system to mother reference.
virtual void MasterToLocal(const Double_t *master, Double_t *local) const
Convert the point coordinates from mother reference to local reference system.
virtual void MasterToLocalVect(const Double_t *master, Double_t *local) const
Convert a vector from mother reference to local reference system.
Double_t Safety(const Double_t *point, Bool_t in=kTRUE) const
computes the closest distance from given point to this shape
The shape encapsulating an assembly (union) of volumes.
virtual void Contains_v(const Double_t *points, Bool_t *inside, Int_t vecsize) const
Check the inside status for each of the points in the array.
virtual void GetMeshNumbers(Int_t &nvert, Int_t &nsegs, Int_t &npols) const
Returns numbers of vertices, segments and polygons composing the shape mesh.
virtual void ComputeNormal_v(const Double_t *points, const Double_t *dirs, Double_t *norms, Int_t vecsize)
Compute the normal for an array o points so that norm.dot.dir is positive Input: Arrays of point coor...
virtual void ComputeBBox()
Compute bounding box of the assembly.
virtual void Safety_v(const Double_t *points, const Bool_t *inside, Double_t *safe, Int_t vecsize) const
Compute safe distance from each of the points in the input array.
TGeoShapeAssembly()
Default constructor.
virtual Double_t DistFromOutside(const Double_t *point, const Double_t *dir, Int_t iact=1, Double_t step=TGeoShape::Big(), Double_t *safe=0) const
compute distance from outside point to surface of the hyperboloid.
virtual void SavePrimitive(std::ostream &out, Option_t *option="")
Save a primitive as a C++ statement(s) on output stream "out".
virtual TGeoShape * GetMakeRuntimeShape(TGeoShape *mother, TGeoMatrix *mat) const
in case shape has some negative parameters, these has to be computed in order to fit the mother
TGeoVolumeAssembly * fVolume
virtual void DistFromInside_v(const Double_t *points, const Double_t *dirs, Double_t *dists, Int_t vecsize, Double_t *step) const
Compute distance from array of input points having directions specified by dirs. Store output in dist...
virtual Bool_t Contains(const Double_t *point) const
Test if point is inside the assembly.
virtual Int_t DistancetoPrimitive(Int_t px, Int_t py)
compute closest distance from point px,py to each vertex. Should not be called.
virtual TGeoVolume * Divide(TGeoVolume *voldiv, const char *divname, Int_t iaxis, Int_t ndiv, Double_t start, Double_t step)
Cannot divide assemblies.
virtual ~TGeoShapeAssembly()
destructor
virtual void SetPoints(Double_t *points) const
No mesh for assemblies.
virtual Double_t DistFromInside(const Double_t *point, const Double_t *dir, Int_t iact=1, Double_t step=TGeoShape::Big(), Double_t *safe=0) const
Compute distance from inside point to surface of the hyperboloid.
void RecomputeBoxLast()
Recompute bounding box of the assembly after adding a node.
virtual void SetSegsAndPols(TBuffer3D &buff) const
Fill TBuffer3D structure for segments and polygons.
virtual void ComputeNormal(const Double_t *point, const Double_t *dir, Double_t *norm)
Compute normal to closest surface from POINT. Should not be called.
virtual void DistFromOutside_v(const Double_t *points, const Double_t *dirs, Double_t *dists, Int_t vecsize, Double_t *step) const
Compute distance from array of input points having directions specified by dirs. Store output in dist...
virtual void InspectShape() const
print shape parameters
virtual Double_t Safety(const Double_t *point, Bool_t in=kTRUE) const
computes the closest distance from given point to this shape, according to option.
Base abstract class for all shapes.
virtual void ComputeNormal(const Double_t *point, const Double_t *dir, Double_t *norm)=0
static Bool_t IsSameWithinTolerance(Double_t a, Double_t b)
Check if two numbers differ with less than a tolerance.
virtual Double_t Safety(const Double_t *point, Bool_t in=kTRUE) const =0
virtual const char * GetName() const
Get the shape name.
virtual Double_t DistFromOutside(const Double_t *point, const Double_t *dir, Int_t iact=1, Double_t step=TGeoShape::Big(), Double_t *safe=0) const =0
virtual void ComputeBBox()=0
virtual Bool_t Contains(const Double_t *point) const =0
virtual Int_t GetCurrentNodeIndex() const
virtual Int_t GetNextNodeIndex() const
void SetNextNodeIndex(Int_t index)
void SetCurrentNodeIndex(Int_t index)
TGeoVolume, TGeoVolumeMulti, TGeoVolumeAssembly are the volume classes.
Int_t GetNdaughters() const
TGeoNode * GetNode(const char *name) const
get the pointer to a daughter node
TGeoVoxelFinder * GetVoxels() const
Getter for optimization structure.
TGeoShape * GetShape() const
virtual Int_t GetCurrentNodeIndex() const
virtual Bool_t IsAssembly() const
Returns true if the volume is an assembly or a scaled assembly.
Finder class handling voxels.
Double_t * GetBoxes() const
virtual Int_t * GetCheckList(const Double_t *point, Int_t &nelem, TGeoStateInfo &td)
get the list of daughter indices for which point is inside their bbox
virtual Int_t * GetNextVoxel(const Double_t *point, const Double_t *dir, Int_t &ncheck, TGeoStateInfo &td)
get the list of new candidates for the next voxel crossed by current ray printf("### GetNextVoxel\n")...
Bool_t IsSafeVoxel(const Double_t *point, Int_t inode, Double_t minsafe) const
Computes squared distance from POINT to the voxel(s) containing node INODE.
virtual void SortCrossedVoxels(const Double_t *point, const Double_t *dir, TGeoStateInfo &td)
get the list in the next voxel crossed by a ray
virtual const char * GetName() const
Returns name of object.
virtual const char * ClassName() const
Returns name of class to which the object belongs.
virtual void Warning(const char *method, const char *msgfmt,...) const
Issue warning message.
virtual void Error(const char *method, const char *msgfmt,...) const
Issue error message.
virtual void Fatal(const char *method, const char *msgfmt,...) const
Issue fatal error message.
virtual void Info(const char *method, const char *msgfmt,...) const
Issue info message.
const char * Data() const
void box(Int_t pat, Double_t x1, Double_t y1, Double_t x2, Double_t y2)
double dist(Rotation3D const &r1, Rotation3D const &r2)
Short_t Min(Short_t a, Short_t b)
Statefull info for the current geometry level.
#define snext(osub1, osub2)