TGeoNode.h

Go to the documentation of this file.

// @(#)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_TGeoNode
#define ROOT_TGeoNode
 
#include "TGeoAtt.h"
 
#include "TNamed.h"
 
#include "TGeoVolume.h"
 
#include "TGeoPatternFinder.h"
 
// forward declarations
class TString;
class TGeoVolume;
class TGeoShape;
class TGeoMedium;
class TGeoMatrix;
class TGeoHMatrix;
class TGeoExtension;
 
////////////////////////////////////////////////////////////////////////////
//                                                                        //
// TGeoNode - base class for logical nodes. They represent volumes        //
//   positioned inside a mother volume                                    //
//                                                                        //
////////////////////////////////////////////////////////////////////////////
 
class TGeoNode : public TNamed,
                 public TGeoAtt
{
protected:
   TGeoVolume       *fVolume = nullptr;        // volume associated with this
   TGeoVolume       *fMother = nullptr;        // mother volume
   Int_t             fNumber = 0;              // copy number
   Int_t             fNovlp = 0;               // number of overlaps
   Int_t            *fOverlaps = nullptr;      //[fNovlp] list of indices for overlapping brothers
   TGeoExtension    *fUserExtension = nullptr; //! Transient user-defined extension to volumes
   TGeoExtension    *fFWExtension = nullptr;   //! Transient framework-defined extension to volumes
 
   void              CopyOverlaps(Int_t *ovlp, Int_t novlp);
 
private:
   TGeoNode(const TGeoNode&) = delete;
   TGeoNode& operator=(const TGeoNode&) = delete;
 
public:
   enum {
      kGeoNodeMatrix = BIT(14),
      kGeoNodeOffset = BIT(15),
      kGeoNodeVC     = BIT(16),
      kGeoNodeOverlap = BIT(17),
      kGeoNodeCloned = BIT(18)
   };
 
   // constructors
   TGeoNode();
   TGeoNode(const TGeoVolume *vol);
   // destructor
   ~TGeoNode() override;
 
   void              Browse(TBrowser *b) override;
   virtual void      cd() const {}
   void              CheckOverlaps(Double_t ovlp=0.1, Option_t *option=""); // *MENU*
   void              CheckShapes();
   Int_t             CountDaughters(Bool_t unique_volumes=kFALSE);
   Int_t     DistancetoPrimitive(Int_t px, Int_t py) override;
   void              Draw(Option_t *option="") override;
   void              DrawOnly(Option_t *option="");
   void              DrawOverlaps();
   void      ExecuteEvent(Int_t event, Int_t px, Int_t py) override;
   void              FillIdArray(Int_t &ifree, Int_t &nodeid, Int_t *array) const;
   Int_t             FindNode(const TGeoNode *node, Int_t level);
   virtual Int_t     GetByteCount() const {return 44;}
   TGeoNode         *GetDaughter(Int_t ind) const {return fVolume->GetNode(ind);}
   virtual TGeoMatrix *GetMatrix() const = 0;
 
   Int_t             GetColour() const {return fVolume->GetLineColor();}
   virtual Int_t     GetIndex() const                    {return 0;}
   virtual TGeoPatternFinder *GetFinder() const          {return nullptr;}
   TGeoMedium       *GetMedium() const                   {return fVolume->GetMedium();}
   TGeoVolume       *GetMotherVolume() const             {return fMother;}
   Int_t             GetNdaughters() const {return fVolume->GetNdaughters();}
   TObjArray        *GetNodes() const {return fVolume->GetNodes();}
   Int_t             GetNumber() const {return fNumber;}
   Int_t            *GetOverlaps(Int_t &novlp) const {novlp=fNovlp; return fOverlaps;}
   TGeoVolume       *GetVolume() const                   {return fVolume;}
   char     *GetObjectInfo(Int_t px, Int_t py) const override;
   virtual Int_t     GetOptimalVoxels() const {return 0;}
   void              InspectNode() const; // *MENU*
   Bool_t            IsCloned() const {return TObject::TestBit(kGeoNodeCloned);}
   Bool_t    IsFolder() const override {return (GetNdaughters()?kTRUE:kFALSE);}
   Bool_t            IsOffset() const {return TObject::TestBit(kGeoNodeOffset);}
   Bool_t            IsOnScreen() const; // *MENU*
   Bool_t            IsOverlapping() const {return TObject::TestBit(kGeoNodeOverlap);}
   Bool_t            IsVirtual() const {return TObject::TestBit(kGeoNodeVC);}
   Bool_t            IsVisible() const {return (TGeoAtt::IsVisible() && fVolume->IsVisible());}
   Bool_t            IsVisDaughters() const {return (TGeoAtt::IsVisDaughters() && fVolume->IsVisDaughters());}
   Bool_t            MayOverlap(Int_t iother) const;
 
   virtual TGeoNode *MakeCopyNode() const {return nullptr;}
   Double_t          Safety(const Double_t *point, Bool_t in=kTRUE) const;
   void              SaveAttributes(std::ostream &out);
   void              SetCurrentPoint(Double_t x, Double_t y, Double_t z) {fVolume->SetCurrentPoint(x,y,z);}// *MENU*
   void              SetVolume(TGeoVolume *volume)       {fVolume = volume;}
   void              SetNumber(Int_t number)             {fNumber=number;}
   void              SetCloned(Bool_t flag=kTRUE)        {TObject::SetBit(kGeoNodeCloned, flag);}
   void              SetOverlapping(Bool_t flag=kTRUE)   {TObject::SetBit(kGeoNodeOverlap, flag);}
   void              SetVirtual()                        {TObject::SetBit(kGeoNodeVC, kTRUE);}
   void              SetVisibility(Bool_t vis=kTRUE) override; // *MENU*
   void              SetInvisible()                      {SetVisibility(kFALSE);} // *MENU*
   void              SetAllInvisible()                   {VisibleDaughters(kFALSE);} // *MENU*
   void              SetMotherVolume(TGeoVolume *mother) {fMother = mother;}
   void              SetOverlaps(Int_t *ovlp, Int_t novlp);
   void              SetUserExtension(TGeoExtension *ext);
   void              SetFWExtension(TGeoExtension *ext);
   TGeoExtension    *GetUserExtension() const {return fUserExtension;}
   TGeoExtension    *GetFWExtension() const   {return fFWExtension;}
   TGeoExtension    *GrabUserExtension() const;
   TGeoExtension    *GrabFWExtension() const;
 
   virtual void      MasterToLocal(const Double_t *master, Double_t *local) const;
   virtual void      MasterToLocalVect(const Double_t *master, Double_t *local) const;
   virtual void      LocalToMaster(const Double_t *local, Double_t *master) const;
   virtual void      LocalToMasterVect(const Double_t *local, Double_t *master) const;
 
   void      ls(Option_t *option = "") const override;
   void      Paint(Option_t *option = "") override;
   void              PrintCandidates() const; // *MENU*
   void              PrintOverlaps() const; // *MENU*
   void              VisibleDaughters(Bool_t vis=kTRUE); // *MENU*
 
   ClassDefOverride(TGeoNode, 2)               // base class for all geometry nodes
};
 
////////////////////////////////////////////////////////////////////////////
//                                                                        //
// TGeoNodeMatrix - node containing a general transformation              //
//                                                                        //
////////////////////////////////////////////////////////////////////////////
 
class TGeoNodeMatrix : public TGeoNode
{
private:
   TGeoMatrix       *fMatrix = nullptr; // transf. matrix of fNode in fMother system
 
   TGeoNodeMatrix(const TGeoNodeMatrix& gnm) = delete;
   TGeoNodeMatrix& operator=(const TGeoNodeMatrix& gnm) = delete;
 
public:
   // constructors
   TGeoNodeMatrix();
   TGeoNodeMatrix(const TGeoVolume *vol, const TGeoMatrix *matrix);
   // destructor
   ~TGeoNodeMatrix() override;
 
   Int_t     GetByteCount() const override;
   Int_t     GetOptimalVoxels() const override;
   Bool_t    IsFolder() const override {return kTRUE;}
   TGeoMatrix *GetMatrix() const override   {return fMatrix;}
   TGeoNode *MakeCopyNode() const override;
   void              SetMatrix(const TGeoMatrix *matrix);
 
   ClassDefOverride(TGeoNodeMatrix, 1)               // a geometry node in the general case
};
 
////////////////////////////////////////////////////////////////////////////
//                                                                        //
// TGeoNodeOffset - node containing only an translation offset            //
//                                                                        //
////////////////////////////////////////////////////////////////////////////
 
class TGeoNodeOffset : public TGeoNode
{
private:
   Double_t          fOffset = 0.;       // X offset for this node with respect to its mother
   Int_t             fIndex = 0;         // index of this node in the division
   TGeoPatternFinder *fFinder = nullptr; // finder for this node
 
   TGeoNodeOffset(const TGeoNodeOffset&) = delete;
   TGeoNodeOffset& operator=(const TGeoNodeOffset&) = delete;
 
public:
   // constructors
   TGeoNodeOffset();
   TGeoNodeOffset(const TGeoVolume *vol, Int_t index, Double_t offset);
   // destructor
   ~TGeoNodeOffset() override;
 
   void      cd() const override           {fFinder->cd(fIndex);}
   Double_t          GetOffset() const {return fOffset;}
   Int_t     GetIndex() const override;
   TGeoPatternFinder *GetFinder() const override {return fFinder;}
   TGeoMatrix *GetMatrix() const override {cd(); return fFinder->GetMatrix();}
   TGeoNode *MakeCopyNode() const override;
   void              SetFinder(TGeoPatternFinder *finder) {fFinder = finder;}
 
   ClassDefOverride(TGeoNodeOffset, 1)      // a geometry node with just an offset
};
 
////////////////////////////////////////////////////////////////////////////
//                                                                        //
// TGeoIteratorPlugin - Plugin for a TGeoIterator providing the method    //
//                      ProcessNode each time Next is called.             //
//                                                                        //
////////////////////////////////////////////////////////////////////////////
 
class TGeoIterator;
 
class TGeoIteratorPlugin : public TObject
{
protected:
   const TGeoIterator *fIterator = nullptr; // Caller iterator
private:
   // No copy
   TGeoIteratorPlugin(const TGeoIteratorPlugin &);
   TGeoIteratorPlugin &operator=(const TGeoIteratorPlugin &);
public:
   TGeoIteratorPlugin() : TObject(),fIterator(nullptr) {}
   ~TGeoIteratorPlugin() override {}
 
   virtual void      ProcessNode() = 0;
   void              SetIterator(const TGeoIterator *iter) {fIterator = iter;}
 
   ClassDefOverride(TGeoIteratorPlugin, 0)  // ABC for user plugins connecter to a geometry iterator.
};
 
////////////////////////////////////////////////////////////////////////////
//                                                                        //
// TGeoIterator - iterator for the node tree                              //
//                                                                        //
////////////////////////////////////////////////////////////////////////////
 
class TGeoIterator
{
private:
   TGeoVolume       *fTop = nullptr;        // Top volume of the iterated branch
   Bool_t            fMustResume = kFALSE;  // Private flag to resume from current node.
   Bool_t            fMustStop = kFALSE;    // Private flag to signal that the iterator has finished.
   Int_t             fLevel = 0;            // Current level in the tree
   Int_t             fType = 0;             // Type of iteration
   Int_t            *fArray = nullptr;      // Array of node indices for the current path
   TGeoHMatrix      *fMatrix = nullptr;     // Current global matrix
   TString           fTopName;              // User name for top
   TGeoIteratorPlugin
                    *fPlugin = nullptr;     // User iterator plugin
   Bool_t            fPluginAutoexec = kFALSE; // Plugin automatically executed during next()
 
   void            IncreaseArray();
protected:
   TGeoIterator() : fTop(nullptr), fMustResume(kFALSE), fMustStop(kFALSE), fLevel(0), fType(0),
                    fArray(nullptr), fMatrix(nullptr), fTopName(), fPlugin(nullptr), fPluginAutoexec(kFALSE) { }
 
public:
   TGeoIterator(TGeoVolume *top);
   TGeoIterator(const TGeoIterator &iter);
   virtual          ~TGeoIterator();
 
   TGeoIterator   &operator=(const TGeoIterator &iter);
   TGeoNode       *operator()();
   TGeoNode       *Next();
   void            Up() { if (fLevel > 0) fLevel--; }
 
   const TGeoMatrix *GetCurrentMatrix() const;
   Int_t           GetIndex(Int_t i) const {return ((i<=fLevel)?fArray[i]:-1);}
   Int_t           GetLevel() const {return fLevel;}
   TGeoNode       *GetNode(Int_t level) const;
   void            GetPath(TString &path) const;
   TGeoIteratorPlugin
                  *GetUserPlugin() const {return fPlugin;}
 
   TGeoVolume     *GetTopVolume() const {return fTop;}
   Int_t           GetType() const {return fType;}
   void            Reset(TGeoVolume *top=nullptr);
   void            SetUserPlugin(TGeoIteratorPlugin *plugin);
   void            SetPluginAutoexec(Bool_t mode) {fPluginAutoexec = mode;}
   void            SetType(Int_t type) {fType = type;}
   void            SetTopName(const char* name);
   void            Skip();
 
   ClassDef(TGeoIterator,0)  //Iterator for geometry.
};
 
#endif