/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/EventDisplay/BesVisLib/BesVisLib-00-04-04/BesVisLib/MucROOTGeo.h

Go to the documentation of this file.

//$id$
/*
 *    2004/9/20    Zhengyun You      Peking University
 *                 Muc Geometry General for EventDisplay
 *
 *    2004/12/11   Zhengyun You      Peking University
 *                 named from MucGeo to MucROOTGeo
 *                 inherit from class SubDetectorROOTGeo
 */

#ifndef Muc_ROOT_GEO_H
#define Muc_ROOT_GEO_H

#include <TGeoVolume.h>
#include <TGeoNode.h>
#include <TGeoPhysicalNode.h>

#include "SubDetectorROOTGeo.h"
#include "BesPolygon2D.h"
#include "Muc2DStrip.h"

class MucROOTGeo : public SubDetectorROOTGeo
{
 public:
  
  MucROOTGeo();

  ~MucROOTGeo();

  void InitFromGDML( const char *gdmlFile, const char *setupName );

  void InitFromROOT( TGeoVolume *vol );

  void Init2DGeometry();

  void SetNode();
  
  void SetVolumeDefaultVis();

  void SetAllVisible();

  void SetQuarterVisible();

  void SetHalfVisible();

  void SetNoEndVisible();

  void SetPhysicalNode();
  
  void SetDetector();
  
  void SetHits();
  
  void SetVisMucDetector();

  void SetVisMucHits();

  int GetPartNb();

  int GetSegNb(int part);
  
  int GetGapNb(int part);

  int GetStripNb(int part, int seg, int gap);

  int GetAbsorberNb(int part);
  
  void SetVolumeMuc(TGeoVolume *vol) { m_Muc = vol; }
  
  TGeoVolume *GetVolumeMuc() { return m_Muc; }
  
  TGeoVolume *GetVolumeAbsorber( int part, int seg, int absorber );

  TGeoVolume *GetVolumeAbsorberSmallBlock( int gap, int sb );

  TGeoVolume *GetVolumeAbsorberPanel( int part, int seg, int absorber, int panel );

  TGeoVolume *GetVolumeGap( int part, int seg, int gap );

  TGeoVolume *GetVolumeBox( int part, int seg, int gap);

  TGeoVolume *GetVolumeBoxSurface( int part, int seg, int gap , int up);

  TGeoVolume *GetVolumeStripPlane( int part, int seg, int gap );

  TGeoVolume *GetVolumeStrip( int part, int seg, int gap, int strip );
  
  TGeoVolume *GetVolumeGasChamber( int part, int seg, int gap, int panel, int gasChamber );

  TGeoVolume *GetVolumeGasBorder( int part, int seg, int gap, int panel, int gasChamber );

  TGeoVolume *GetVolumeBakelite( int part, int seg, int gap, int RpcUpDown, int panel, int bakelite );

  TGeoNode *GetAbsorber( int part, int seg, int absorber );

  TGeoNode *GetAbsorberPanel( int part, int seg, int absorber, int panel );

  TGeoNode *GetGap( int part, int seg, int gap );

  TGeoNode *GetStripPlane( int part, int seg, int gap);

  TGeoNode *GetStrip( int part, int seg, int gap, int strip );
  
  //TGeoNode *GetGasChamber( int part, int seg, int gap, int gasChamber );

  //TGeoNode *GetBakelite( int part, int seg, int gap, int bakelite );

  TGeoPhysicalNode *GetPhysicalAbsorber( int part, int seg, int gap, int panel );

  TGeoPhysicalNode *GetPhysicalGap( int part, int seg, int gap );

  TGeoPhysicalNode *GetPhysicalStrip( int part, int seg, int gap, int strip );

  Muc2DStrip *Get2DStrip( int part, int seg, int gap, int strip ); 

  Int_t GetPart(TGeoPhysicalNode* phyNode);

  Bool_t IsZRVisible( int part, int seg );

  Double_t Range360(Double_t input);

  void Draw(Option_t *option);

  void DrawHits(Option_t *option);
 
 private:

  static const int m_kPart        = 3;
  static const int m_kSegMax      = 8;
  static const int m_kAbsorberMax = 9;
  static const int m_kGapMax      = 9;
  static const int m_kPanelMax    = 4;
  static const int m_kStripMax    = 112;
  static const int m_kBakelite    = 2;
  static const int m_kGasChamber  = 1;
  static const int m_kSmallBlockMax = 2;
  static const int m_kUpDown      = 2;

  static const int m_kSeg[m_kPart];
  static const int m_kAbsorber[m_kPart];
  static const int m_kGap[m_kPart];
  static const int m_kPanel[m_kPart];
  static const int m_kBakelitePanel[m_kPart][m_kUpDown];

  int m_StripNum[m_kPart][m_kSegMax][m_kGapMax];

  int m_MucColor;
  int m_absorberColor;
  int m_gapColor;
  int m_gasChamberColor;
  int m_bakeliteColor;
  int m_stripColor;

  static const Int_t m_kStripZMuliple = 12;  // as dz is very thin, multiply so that it can be seen.

  TGeoVolume       *m_Muc;
  TGeoNode         *m_NodeAbsorber[m_kPart][m_kSegMax][m_kAbsorberMax];
  TGeoNode         *m_NodeAbsorberPanel[m_kPart][m_kSegMax][m_kAbsorberMax][m_kPanelMax];
  TGeoNode         *m_NodeGap[m_kPart][m_kSegMax][m_kGapMax];
  TGeoNode         *m_NodeStripPlane[m_kPart][m_kSegMax][m_kGapMax];
  TGeoNode         *m_NodeStrip[m_kPart][m_kSegMax][m_kGapMax][m_kStripMax];

  TGeoPhysicalNode *m_PhysicalAbsorber[m_kPart][m_kSegMax][m_kAbsorberMax][m_kPanelMax];
  TGeoPhysicalNode *m_PhysicalGap[m_kPart][m_kSegMax][m_kGapMax];
  TGeoPhysicalNode *m_PhysicalStrip[m_kPart][m_kSegMax][m_kGapMax][m_kStripMax];

  const TObjArray  *m_MucDigiCol;

  // 2D
  BesPolygon2D     *m_MucXY[m_kPart][m_kSegMax];
  BesPolygon2D     *m_MucZR[m_kPart][m_kSegMax];
  BesPolygon2D     *m_MucXYGap[m_kPart][m_kSegMax][m_kGapMax];
  BesPolygon2D     *m_MucZRGap[m_kPart][m_kSegMax][m_kGapMax];
  
  Muc2DStrip       *m_Muc2DStrip[m_kPart][m_kSegMax][m_kGapMax][m_kStripMax];
};

#endif   /* Muc_ROOT_GEO_H */

---