// @(#)root/geom:\$Id\$
// Author: Mihaela Gheata   24/01/04

/*************************************************************************
*                                                                       *
* For the licensing terms see \$ROOTSYS/LICENSE.                         *
* For the list of contributors see \$ROOTSYS/README/CREDITS.             *
*************************************************************************/

#ifndef ROOT_TGeoXtru
#define ROOT_TGeoXtru

#ifndef ROOT_TGeoBBox
#include "TGeoBBox.h"
#endif

class TGeoPolygon;

////////////////////////////////////////////////////////////////////////////
//                                                                        //
// TGeoXtru - An extrusion with fixed outline shape in x-y and a sequence //
//   of z extents (segments).  The overall scale of the outline scales    //
//   linearly between z points and the center can have an x-y offset.     //
//                                                                        //
//   Based on the initial implementation of R. Hatcher                    //
//                                                                        //
////////////////////////////////////////////////////////////////////////////

class TGeoXtru : public TGeoBBox
{
public:
{
Int_t              fSeg;   // !current segment [0,fNvert-1]
Int_t              fIz;    // !current z plane [0,fNz-1]
Double_t          *fXc;    // ![fNvert] current X positions for polygon vertices
Double_t          *fYc;    // ![fNvert] current Y positions for polygon vertices
TGeoPolygon       *fPoly;  // !polygon defining section shape

};

protected:
// data members
Int_t                 fNvert; // number of vertices of the 2D polygon (at least 3)
Int_t                 fNz;    // number of z planes (at least two)
Double_t              fZcurrent; // current Z position
Double_t             *fX;     //[fNvert] X positions for polygon vertices
Double_t             *fY;     //[fNvert] Y positions for polygon vertices
Double_t             *fZ;     //[fNz] array of Z planes positions
Double_t             *fScale; //[fNz] array of scale factors (for each Z)
Double_t             *fX0;    //[fNz] array of X offsets (for each Z)
Double_t             *fY0;    //[fNz] array of Y offsets (for each Z)

TGeoXtru(const TGeoXtru&);
TGeoXtru& operator=(const TGeoXtru&);

// methods
Double_t              DistToPlane(const Double_t *point, const Double_t *dir, Int_t iz, Int_t ivert, Double_t stepmax, Bool_t in) const;
void                  GetPlaneVertices(Int_t iz, Int_t ivert, Double_t *vert) const;
void                  GetPlaneNormal(const Double_t *vert, Double_t *norm) const;
Bool_t                IsPointInsidePlane(const Double_t *point, Double_t *vert, Double_t *norm) const;
Double_t              SafetyToSector(const Double_t *point, Int_t iz, Double_t safmin, Bool_t in);
void                  SetIz(Int_t iz);
void                  SetSeg(Int_t iseg);

public:
// constructors
TGeoXtru();
TGeoXtru(Int_t nz);
TGeoXtru(Double_t *param);
// destructor
virtual ~TGeoXtru();
// methods
virtual Double_t      Capacity() const;
virtual void          ComputeBBox();
virtual void          ComputeNormal(const Double_t *point, const Double_t *dir, Double_t *norm);
virtual void          ComputeNormal_v(const Double_t *points, const Double_t *dirs, Double_t *norms, Int_t vecsize);
virtual Bool_t        Contains(const Double_t *point) const;
virtual void          Contains_v(const Double_t *points, Bool_t *inside, Int_t vecsize) const;
Bool_t                DefinePolygon(Int_t nvert, const Double_t *xv, const Double_t *yv);
virtual void          DefineSection(Int_t snum, Double_t z, Double_t x0=0., Double_t y0=0., Double_t scale=1.);
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;
virtual void          DistFromInside_v(const Double_t *points, const Double_t *dirs, Double_t *dists, Int_t vecsize, Double_t *step) const;
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;
virtual void          DistFromOutside_v(const Double_t *points, const Double_t *dirs, Double_t *dists, Int_t vecsize, Double_t *step) const;
virtual Int_t         DistancetoPrimitive(Int_t px, Int_t py);
void                  DrawPolygon(Option_t *option="");
virtual const TBuffer3D &GetBuffer3D(Int_t reqSections, Bool_t localFrame) const;
//   virtual Int_t         GetByteCount() const {return 60+12*fNz;}
Int_t                 GetNz() const    {return fNz;}
Int_t                 GetNvert() const {return fNvert;}
Double_t              GetX(Int_t i) const {return (i<fNvert&&i>-1 &&fX!=0) ? fX[i] : -1.0E10;}
Double_t              GetY(Int_t i) const {return (i<fNvert&&i>-1 &&fY!=0) ? fY[i] : -1.0E10;}
Double_t              GetXOffset(Int_t i) const {return (i<fNz&&i>-1 && fX0!=0) ? fX0[i] : 0.0;}
Double_t              GetYOffset(Int_t i) const {return (i<fNz&&i>-1 && fY0!=0) ? fY0[i] : 0.0;}
Double_t              GetScale(Int_t i) const {return (i<fNz&&i>-1 && fScale!=0) ? fScale[i] : 1.0;}
Double_t             *GetZ() const     {return fZ;}
Double_t              GetZ(Int_t ipl) const;
virtual TGeoShape    *GetMakeRuntimeShape(TGeoShape * /*mother*/, TGeoMatrix * /*mat*/) const {return 0;}
virtual void          GetMeshNumbers(Int_t &nvert, Int_t &nsegs, Int_t &npols) const;
virtual Int_t         GetNmeshVertices() const;
virtual void          InspectShape() const;
virtual TBuffer3D    *MakeBuffer3D() const;
Double_t             &Z(Int_t ipl) {return fZ[ipl];}
virtual Double_t      Safety(const Double_t *point, Bool_t in=kTRUE) const;
virtual void          Safety_v(const Double_t *points, const Bool_t *inside, Double_t *safe, Int_t vecsize) const;
virtual void          SavePrimitive(std::ostream &out, Option_t *option = "");
void                  SetCurrentZ(Double_t z, Int_t iz);
void                  SetCurrentVertices(Double_t x0, Double_t y0, Double_t scale);
virtual void          SetDimensions(Double_t *param);
virtual void          SetPoints(Double_t *points) const;
virtual void          SetPoints(Float_t *points) const;
virtual void          SetSegsAndPols(TBuffer3D &buff) const;
virtual void          Sizeof3D() const;

ClassDef(TGeoXtru, 3)         // extruded polygon class
};

#endif
TGeoXtru.h:1
TGeoXtru.h:2
TGeoXtru.h:3
TGeoXtru.h:4
TGeoXtru.h:5
TGeoXtru.h:6
TGeoXtru.h:7
TGeoXtru.h:8
TGeoXtru.h:9
TGeoXtru.h:10
TGeoXtru.h:11
TGeoXtru.h:12
TGeoXtru.h:13
TGeoXtru.h:14
TGeoXtru.h:15
TGeoXtru.h:16
TGeoXtru.h:17
TGeoXtru.h:18
TGeoXtru.h:19
TGeoXtru.h:20
TGeoXtru.h:21
TGeoXtru.h:22
TGeoXtru.h:23
TGeoXtru.h:24
TGeoXtru.h:25
TGeoXtru.h:26
TGeoXtru.h:27
TGeoXtru.h:28
TGeoXtru.h:29
TGeoXtru.h:30
TGeoXtru.h:31
TGeoXtru.h:32
TGeoXtru.h:33
TGeoXtru.h:34
TGeoXtru.h:35
TGeoXtru.h:36
TGeoXtru.h:37
TGeoXtru.h:38
TGeoXtru.h:39
TGeoXtru.h:40
TGeoXtru.h:41
TGeoXtru.h:42
TGeoXtru.h:43
TGeoXtru.h:44
TGeoXtru.h:45
TGeoXtru.h:46
TGeoXtru.h:47
TGeoXtru.h:48
TGeoXtru.h:49
TGeoXtru.h:50
TGeoXtru.h:51
TGeoXtru.h:52
TGeoXtru.h:53
TGeoXtru.h:54
TGeoXtru.h:55
TGeoXtru.h:56
TGeoXtru.h:57
TGeoXtru.h:58
TGeoXtru.h:59
TGeoXtru.h:60
TGeoXtru.h:61
TGeoXtru.h:62
TGeoXtru.h:63
TGeoXtru.h:64
TGeoXtru.h:65
TGeoXtru.h:66
TGeoXtru.h:67
TGeoXtru.h:68
TGeoXtru.h:69
TGeoXtru.h:70
TGeoXtru.h:71
TGeoXtru.h:72
TGeoXtru.h:73
TGeoXtru.h:74
TGeoXtru.h:75
TGeoXtru.h:76
TGeoXtru.h:77
TGeoXtru.h:78
TGeoXtru.h:79
TGeoXtru.h:80
TGeoXtru.h:81
TGeoXtru.h:82
TGeoXtru.h:83
TGeoXtru.h:84
TGeoXtru.h:85
TGeoXtru.h:86
TGeoXtru.h:87
TGeoXtru.h:88
TGeoXtru.h:89
TGeoXtru.h:90
TGeoXtru.h:91
TGeoXtru.h:92
TGeoXtru.h:93
TGeoXtru.h:94
TGeoXtru.h:95
TGeoXtru.h:96
TGeoXtru.h:97
TGeoXtru.h:98
TGeoXtru.h:99
TGeoXtru.h:100
TGeoXtru.h:101
TGeoXtru.h:102
TGeoXtru.h:103
TGeoXtru.h:104
TGeoXtru.h:105
TGeoXtru.h:106
TGeoXtru.h:107
TGeoXtru.h:108
TGeoXtru.h:109
TGeoXtru.h:110
TGeoXtru.h:111
TGeoXtru.h:112
TGeoXtru.h:113
TGeoXtru.h:114
TGeoXtru.h:115
TGeoXtru.h:116
TGeoXtru.h:117
TGeoXtru.h:118
TGeoXtru.h:119
TGeoXtru.h:120
TGeoXtru.h:121
TGeoXtru.h:122
TGeoXtru.h:123
TGeoXtru.h:124
TGeoXtru.h:125
TGeoXtru.h:126
TGeoXtru.h:127
TGeoXtru.h:128
TGeoXtru.h:129
TGeoXtru.h:130