ROOT logo
// @(#)root/geom:$Id$
// Author: Mihaela Gheata   20/06/04

/*************************************************************************
 * 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_TGeoParaboloid
#define ROOT_TGeoParaboloid

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

////////////////////////////////////////////////////////////////////////////
//                                                                        
// TGeoParaboloid - Paraboloid  class. A paraboloid is the solid bounded by
//            the following surfaces:
//            - 2 planes parallel with XY cutting the Z axis at Z=-dz and Z=+dz
//            - the surface of revolution of a parabola described by:
//                 z = a*(x*x + y*y) + b
//       The parameters a and b are automatically computed from:
//            - rlo - the radius of the circle of intersection between the 
//              parabolic surface and the plane z = -dz
//            - rhi - the radius of the circle of intersection between the 
//              parabolic surface and the plane z = +dz
//         | -dz = a*rlo*rlo + b
//         |  dz = a*rhi*rhi + b      where: rlo != rhi, both >= 0
//                                                                        
////////////////////////////////////////////////////////////////////////////


class TGeoParaboloid : public TGeoBBox
{
private:
   Double_t              fRlo;                  // radius at z=-dz
   Double_t              fRhi;                  // radius at z=+dz
   Double_t              fDz;                   // range on Z axis [-dz, dz]
   Double_t              fA;                    // quadratic coeff.
   Double_t              fB;                    // Z value of parabola at x=y=0
public:
   // constructors
   TGeoParaboloid();
   TGeoParaboloid(Double_t rlo, Double_t rhi, Double_t dz);
   TGeoParaboloid(const char *name, Double_t rlo, Double_t rhi, Double_t dz);
   TGeoParaboloid(Double_t *params);
   // destructor
   virtual ~TGeoParaboloid();
   // 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 Int_t         DistancetoPrimitive(Int_t px, Int_t py);
   Double_t              DistToParaboloid(const Double_t *point, const Double_t *dir, Bool_t in) const;
   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 TGeoVolume   *Divide(TGeoVolume *voldiv, const char *divname, Int_t iaxis, Int_t ndiv, 
                                Double_t start, Double_t step);
   virtual const TBuffer3D &GetBuffer3D(Int_t reqSections, Bool_t localFrame) const;
   Double_t              GetRlo() const    {return fRlo;}
   Double_t              GetRhi() const    {return fRhi;}
   Double_t              GetDz() const     {return fDz;}
   
   virtual void          GetBoundingCylinder(Double_t *param) const;
   virtual TGeoShape    *GetMakeRuntimeShape(TGeoShape *mother, TGeoMatrix *mat) const;
   virtual void          GetMeshNumbers(Int_t &nvert, Int_t &nsegs, Int_t &npols) const;
   virtual Int_t         GetNmeshVertices() const;
   virtual Bool_t        GetPointsOnSegments(Int_t /*npoints*/, Double_t * /*array*/) const {return kFALSE;}
   virtual void          InspectShape() const;
   virtual Bool_t        IsCylType() const {return kTRUE;}
   virtual TBuffer3D    *MakeBuffer3D() const;
   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                  SetParaboloidDimensions(Double_t rlo, Double_t rhi, Double_t dz);
   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(TGeoParaboloid, 1)         // paraboloid class

};


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