// @(#)root/base:$Id$
// Author: Andrei Gheata   28/07/03

/*************************************************************************
 * 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_TGeoTorus
#define ROOT_TGeoTorus

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

////////////////////////////////////////////////////////////////////////////
//                                                                        //
// TGeoTorus - Torus segment class. A torus has 5 parameters :            //
//            R    - axial radius                                         //
//            Rmin - inner radius                                         //
//            Rmax - outer radius                                         //
//            Phi1 - starting phi                                         //
//            Dphi - phi extent                                           //
//                                                                        //
////////////////////////////////////////////////////////////////////////////

class TGeoTorus : public TGeoBBox
{
protected :
// data members
   Double_t              fR;    // axial radius
   Double_t              fRmin; // inner radius
   Double_t              fRmax; // outer radius
   Double_t              fPhi1; // starting phi
   Double_t              fDphi; // phi extent
// methods

public:
   virtual Double_t      Capacity() const;
   Double_t              Daxis(const Double_t *pt, const Double_t *dir, Double_t t) const;
   Double_t              DDaxis(const Double_t *pt, const Double_t *dir, Double_t t) const;
   Double_t              DDDaxis(const Double_t *pt, const Double_t *dir, Double_t t) const;
   Double_t              ToBoundary(const Double_t *pt, const Double_t *dir, Double_t r, Bool_t in) const;
   Int_t                 SolveCubic(Double_t a, Double_t b, Double_t c, Double_t *x) const;
   Int_t                 SolveQuartic(Double_t a, Double_t b, Double_t c, Double_t d, Double_t *x) const;
public:
   // constructors
   TGeoTorus();
   TGeoTorus(Double_t r, Double_t rmin, Double_t rmax, Double_t phi1=0, Double_t dphi=360);
   TGeoTorus(const char * name, Double_t r, Double_t rmin, Double_t rmax, Double_t phi1=0, Double_t dphi=360);
   TGeoTorus(Double_t *params);
   // destructor
   virtual ~TGeoTorus() {}
   // methods

   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 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);
   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 56;}
   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;}
   Double_t              GetR() const    {return fR;}
   Double_t              GetRmin() const {return fRmin;}
   Double_t              GetRmax() const {return fRmax;}
   Double_t              GetPhi1() const {return fPhi1;}
   Double_t              GetDphi() const {return fDphi;}
   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                  SetTorusDimensions(Double_t r, Double_t rmin, Double_t rmax, Double_t phi1, Double_t dphi);
   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(TGeoTorus, 1)         // torus class

};

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