// @(#)root/geom:$Id: TGeoCone.h 24870 2008-07-17 16:30:58Z brun $ // 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_TGeoCone #define ROOT_TGeoCone #ifndef ROOT_TGeoBBox #include "TGeoBBox.h" #endif //////////////////////////////////////////////////////////////////////////// // // // TGeoCone - conical tube class. It has 5 parameters : // // dz - half length in z // // Rmin1, Rmax1 - inside and outside radii at -dz // // Rmin2, Rmax2 - inside and outside radii at +dz // // // //////////////////////////////////////////////////////////////////////////// class TGeoCone : public TGeoBBox { protected : // data members Double_t fDz; // half length Double_t fRmin1; // inner radius at -dz Double_t fRmax1; // outer radius at -dz Double_t fRmin2; // inner radius at +dz Double_t fRmax2; // outer radius at +dz // methods public: // constructors TGeoCone(); TGeoCone(Double_t dz, Double_t rmin1, Double_t rmax1, Double_t rmin2, Double_t rmax2); TGeoCone(const char *name, Double_t dz, Double_t rmin1, Double_t rmax1, Double_t rmin2, Double_t rmax2); TGeoCone(Double_t *params); // destructor virtual ~TGeoCone(); // methods virtual Double_t Capacity() const; static Double_t Capacity(Double_t dz, Double_t rmin1, Double_t rmax1, Double_t rmin2, Double_t rmax2); virtual void ComputeBBox(); virtual void ComputeNormal(Double_t *point, Double_t *dir, Double_t *norm); static void ComputeNormalS(Double_t *point, Double_t *dir, Double_t *norm, Double_t dz, Double_t rmin1, Double_t rmax1, Double_t rmin2, Double_t rmax2); virtual Bool_t Contains(Double_t *point) const; virtual Int_t DistancetoPrimitive(Int_t px, Int_t py); static void DistToCone(Double_t *point, Double_t *dir, Double_t dz, Double_t r1, Double_t r2, Double_t &b, Double_t &delta); static Double_t DistFromInsideS(Double_t *point, Double_t *dir, Double_t dz, Double_t rmin1, Double_t rmax1, Double_t rmin2, Double_t rmax2); virtual Double_t DistFromInside(Double_t *point, Double_t *dir, Int_t iact=1, Double_t step=TGeoShape::Big(), Double_t *safe=0) const; static Double_t DistFromOutsideS(Double_t *point, Double_t *dir, Double_t dz, Double_t rmin1, Double_t rmax1, Double_t rmin2, Double_t rmax2); virtual Double_t DistFromOutside(Double_t *point, Double_t *dir, Int_t iact=1, Double_t step=TGeoShape::Big(), Double_t *safe=0) const; 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 Int_t GetByteCount() const {return 56;} virtual const TBuffer3D &GetBuffer3D(Int_t reqSections, Bool_t localFrame) const; virtual Double_t GetDz() const {return fDz;} 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; virtual Double_t GetRmin1() const {return fRmin1;} virtual Double_t GetRmax1() const {return fRmax1;} virtual Double_t GetRmin2() const {return fRmin2;} virtual Double_t GetRmax2() const {return fRmax2;} virtual void InspectShape() const; virtual Bool_t IsCylType() const {return kTRUE;} virtual TBuffer3D *MakeBuffer3D() const; virtual Double_t Safety(Double_t *point, Bool_t in=kTRUE) const; static Double_t SafetyS(Double_t *point, Bool_t in, Double_t dz, Double_t rmin1, Double_t rmax1, Double_t rmin2, Double_t rmax2, Int_t skipz=0); virtual void SavePrimitive(ostream &out, Option_t *option = ""); void SetConeDimensions(Double_t dz, Double_t rmin1, Double_t rmax1, Double_t rmin2, Double_t rmax2); virtual void SetDimensions(Double_t *param); virtual void SetPoints(Double_t *points) const; virtual void SetPoints(Float_t *points) const; virtual void SetSegsAndPols(TBuffer3D &buffer) const; virtual void Sizeof3D() const; ClassDef(TGeoCone, 1) // conical tube class }; //////////////////////////////////////////////////////////////////////////// // // // TGeoConeSeg - a phi segment of a conical tube. Has 7 parameters : // // - the same 5 as a cone; // // - first phi limit (in degrees) // // - second phi limit // // // //////////////////////////////////////////////////////////////////////////// class TGeoConeSeg : public TGeoCone { protected: // data members Double_t fPhi1; // first phi limit Double_t fPhi2; // second phi limit public: // constructors TGeoConeSeg(); TGeoConeSeg(Double_t dz, Double_t rmin1, Double_t rmax1, Double_t rmin2, Double_t rmax2, Double_t phi1, Double_t phi2); TGeoConeSeg(const char *name, Double_t dz, Double_t rmin1, Double_t rmax1, Double_t rmin2, Double_t rmax2, Double_t phi1, Double_t phi2); TGeoConeSeg(Double_t *params); // destructor virtual ~TGeoConeSeg(); // methods virtual Double_t Capacity() const; static Double_t Capacity(Double_t dz, Double_t rmin1, Double_t rmax1, Double_t rmin2, Double_t rmax2, Double_t phi1, Double_t phi2); virtual void ComputeBBox(); virtual void ComputeNormal(Double_t *point, Double_t *dir, Double_t *norm); static void ComputeNormalS(Double_t *point, Double_t *dir, Double_t *norm, Double_t dz, Double_t rmin1, Double_t rmax1, Double_t rmin2, Double_t rmax2, Double_t c1, Double_t s1, Double_t c2, Double_t s2); virtual Bool_t Contains(Double_t *point) const; virtual Int_t DistancetoPrimitive(Int_t px, Int_t py); static Double_t DistToCons(Double_t *point, Double_t *dir, Double_t r1, Double_t z1, Double_t r2, Double_t z2, Double_t phi1, Double_t phi2); static Double_t DistFromInsideS(Double_t *point, Double_t *dir, Double_t dz, Double_t rmin1, Double_t rmax1, Double_t rmin2, Double_t rmax2, Double_t c1, Double_t s1, Double_t c2, Double_t s2, Double_t cm, Double_t sm, Double_t cdfi); virtual Double_t DistFromInside(Double_t *point, Double_t *dir, Int_t iact=1, Double_t step=TGeoShape::Big(), Double_t *safe=0) const; static Double_t DistFromOutsideS(Double_t *point, Double_t *dir, Double_t dz, Double_t rmin1, Double_t rmax1, Double_t rmin2, Double_t rmax2, Double_t c1, Double_t s1, Double_t c2, Double_t s2, Double_t cm, Double_t sm, Double_t cdfi); virtual Double_t DistFromOutside(Double_t *point, Double_t *dir, Int_t iact=1, Double_t step=TGeoShape::Big(), Double_t *safe=0) 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 const TBuffer3D &GetBuffer3D(Int_t reqSections, Bool_t localFrame) const; virtual Int_t GetByteCount() const {return 64;} 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; Double_t GetPhi1() const {return fPhi1;} Double_t GetPhi2() const {return fPhi2;} virtual void InspectShape() const; virtual TBuffer3D *MakeBuffer3D() const; virtual Double_t Safety(Double_t *point, Bool_t in=kTRUE) const; static Double_t SafetyS(Double_t *point, Bool_t in, Double_t dz, Double_t rmin1, Double_t rmax1, Double_t rmin2, Double_t rmax2, Double_t phi1, Double_t phi2, Int_t skipz=0); virtual void SavePrimitive(ostream &out, Option_t *option = ""); void SetConsDimensions(Double_t dz, Double_t rmin1, Double_t rmax1, Double_t rmin2, Double_t rmax2, Double_t phi1, Double_t phi2); virtual void SetDimensions(Double_t *param); virtual void SetPoints(Double_t *points) const; virtual void SetPoints(Float_t *points) const; virtual void SetSegsAndPols(TBuffer3D &buffer) const; virtual void Sizeof3D() const; ClassDef(TGeoConeSeg, 1) // conical tube segment class }; #endif