ROOT logo
// @(#)root/geom:$Id$
// Author: Andrei Gheata   24/10/01

/*************************************************************************
 * Copyright (C) 1995-2000, Rene Brun and Fons Rademakers.               *
 * All rights reserved.                                                  *
 *                                                                       *
 * For the licensing terms see $ROOTSYS/LICENSE.                         *
 * For the list of contributors see $ROOTSYS/README/CREDITS.             *
 *************************************************************************/

#ifndef ROOT_TGeoPcon
#define ROOT_TGeoPcon

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

  
//////////////////////////////////////////////////////////////////////////////
//                                                                          //
// TGeoPcon - a composite polycone. It has at least 9 parameters :          //
//            - the lower phi limit;                                        //
//            - the range in phi;                                           //
//            - the number of z planes (at least two) where the inner/outer //
//              radii are changing;                                         //
//            - z coordinate, inner and outer radius for each z plane       //
//                                                                          //
//////////////////////////////////////////////////////////////////////////////

class TGeoPcon : public TGeoBBox
{
protected:
   // data members
   Int_t                 fNz;      // number of z planes (at least two)
   Double_t              fPhi1;    // lower phi limit (converted to [0,2*pi)
   Double_t              fDphi;    // phi range
   Double_t             *fRmin;    //[fNz] pointer to array of inner radii 
   Double_t             *fRmax;    //[fNz] pointer to array of outer radii 
   Double_t             *fZ;       //[fNz] pointer to array of Z planes positions 
   Bool_t                fFullPhi; //! Full phi range flag
   Double_t              fC1;      //! Cosine of phi1
   Double_t              fS1;      //! Sine of phi1
   Double_t              fC2;      //! Cosine of phi1+dphi
   Double_t              fS2;      //! Sine of phi1+dphi
   Double_t              fCm;      //! Cosine of (phi1+phi2)/2
   Double_t              fSm;      //! Sine of (phi1+phi2)/2
   Double_t              fCdphi;   //! Cosine of dphi
   
   // methods
   TGeoPcon(const TGeoPcon&); 
   TGeoPcon& operator=(const TGeoPcon&); 

public:
   // constructors
   TGeoPcon();
   TGeoPcon(Double_t phi, Double_t dphi, Int_t nz);
   TGeoPcon(const char *name, Double_t phi, Double_t dphi, Int_t nz);
   TGeoPcon(Double_t *params);
   // destructor
   virtual ~TGeoPcon();
   // 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;
   virtual void          DefineSection(Int_t snum, Double_t z, Double_t rmin, Double_t rmax);
   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;
   Double_t              DistToSegZ(const Double_t *point, const Double_t *dir, Int_t &iz) const;
   virtual Int_t         DistancetoPrimitive(Int_t px, Int_t py);
   virtual TGeoVolume   *Divide(TGeoVolume *voldiv, const char *divname, Int_t iaxis, Int_t ndiv, 
                                Double_t start, Double_t step);
   virtual const char   *GetAxisName(Int_t iaxis) const;
   virtual Double_t      GetAxisRange(Int_t iaxis, Double_t &xlo, Double_t &xhi) const;
   virtual void          GetBoundingCylinder(Double_t *param) const;
   virtual const TBuffer3D &GetBuffer3D(Int_t reqSections, Bool_t localFrame) const;
   virtual Int_t         GetByteCount() const {return 60+12*fNz;}
   Double_t              GetPhi1() const {return fPhi1;}
   Double_t              GetDphi() const {return fDphi;}
   Int_t                 GetNz() const   {return fNz;}
   virtual Int_t         GetNsegments() const;
   Double_t             *GetRmin() const {return fRmin;}
   Double_t              GetRmin(Int_t ipl) const;
   Double_t             *GetRmax() const {return fRmax;}
   Double_t              GetRmax(Int_t ipl) const;
   Double_t             *GetZ() const    {return fZ;}
   Double_t              GetZ(Int_t ipl) const;
   virtual TGeoShape    *GetMakeRuntimeShape(TGeoShape * /*mother*/, TGeoMatrix * /*mat*/) const {return 0;}
   virtual Int_t         GetNmeshVertices() const;
   virtual Bool_t        GetPointsOnSegments(Int_t /*npoints*/, Double_t * /*array*/) const {return kFALSE;}
   virtual void          GetMeshNumbers(Int_t &nvert, Int_t &nsegs, Int_t &npols) const;
   virtual void          InspectShape() const;
   virtual Bool_t        IsCylType() const {return kTRUE;}
   virtual TBuffer3D    *MakeBuffer3D() const;
   Double_t             &Phi1()          {return fPhi1;}
   Double_t             &Dphi()          {return fDphi;}
   Double_t             &Rmin(Int_t ipl) {return fRmin[ipl];}
   Double_t             &Rmax(Int_t ipl) {return fRmax[ipl];}
   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;
   Double_t              SafetyToSegment(const Double_t *point, Int_t ipl, Bool_t in=kTRUE, Double_t safmin=TGeoShape::Big()) const;
   virtual void          SavePrimitive(std::ostream &out, Option_t *option = "");
   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(TGeoPcon, 1)         // polycone class 
};

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