// @(#)root/geom:$Id$
// Author: Andrei Gheata   31/01/02

/*************************************************************************
 * 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_TGeoTrd2
#define ROOT_TGeoTrd2

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


////////////////////////////////////////////////////////////////////////////
//                                                                        //
// TGeoTrd2 - a trapezoid with both x and y lengths varying with z. It    //
//   has 5 parameters, the half lengths in x at -dz and +dz, the half     //
//  lengths in y at -dz and +dz, and the half length in z (dz).           //
//                                                                        //
////////////////////////////////////////////////////////////////////////////

class TGeoTrd2 : public TGeoBBox
{
protected:
   // data members
   Double_t              fDx1; // half length in X at lower Z surface (-dz)
   Double_t              fDx2; // half length in X at higher Z surface (+dz)
   Double_t              fDy1; // half length in Y at lower Z surface (-dz)
   Double_t              fDy2; // half length in Y at higher Z surface (+dz)
   Double_t              fDz;  // half length in Z

   // methods

public:
   // constructors
   TGeoTrd2();
   TGeoTrd2(Double_t dx1, Double_t dx2, Double_t dy1, Double_t dy2, Double_t dz);
   TGeoTrd2(const char *name, Double_t dx1, Double_t dx2, Double_t dy1, Double_t dy2, Double_t dz);
   TGeoTrd2(Double_t *params);
   // destructor
   virtual ~TGeoTrd2();
   // methods

   virtual Double_t      Capacity() const;
   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          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 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 Double_t      GetAxisRange(Int_t iaxis, Double_t &xlo, Double_t &xhi) const;
   virtual void          GetBoundingCylinder(Double_t *param) const;
   virtual Int_t         GetByteCount() const {return 56;}
   Double_t              GetDx1() const {return fDx1;}
   Double_t              GetDx2() const {return fDx2;}
   Double_t              GetDy1() const {return fDy1;}
   Double_t              GetDy2() const {return fDy2;}
   Double_t              GetDz() const  {return fDz;}
   virtual Int_t         GetFittingBox(const TGeoBBox *parambox, TGeoMatrix *mat, Double_t &dx, Double_t &dy, Double_t &dz) const;
   virtual TGeoShape    *GetMakeRuntimeShape(TGeoShape *mother, TGeoMatrix *mat) const;
   void                  GetVisibleCorner(const Double_t *point, Double_t *vertex, Double_t *normals) const;
   void                  GetOppositeCorner(const Double_t *point, Int_t inorm, Double_t *vertex, Double_t *normals) const;
   virtual void          InspectShape() const;
   virtual Bool_t        IsCylType() const {return kFALSE;}
   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 = "");
   virtual void          SetDimensions(Double_t *param);
   virtual void          SetPoints(Double_t *points) const;
   virtual void          SetPoints(Float_t *points) const;
   void                  SetVertex(Double_t *vertex) const;
   virtual void          Sizeof3D() const;

   ClassDef(TGeoTrd2, 1)         // TRD2 shape class
};

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