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

Go to the documentation of this file.

//$id$
/*
 *    2004/11/29   Zhengyun You      Peking University
 *                 Emc Geometry General for EventDisplay
 *
 *    2004/12/11   Zhengyun You      Peking University
 *                 named from EmcGeo to EmcROOTGeo
 *                 inherit from class SubDetectorROOTGeo
 *
 *    2005/04/04   Zhengyun You      Peking University
 *                 Emc End geometry added
 */

#ifndef Emc_ROOT_GEO_H
#define Emc_ROOT_GEO_H

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

#include "SubDetectorROOTGeo.h"
#include "Emc2DCrystal.h"

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

  ~EmcROOTGeo();

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

  void InitFromROOT( TGeoVolume *vol );

  void Init2DGeometry();

  void SetNode();

  void SetVolumeAppendInVis();

  void SetVolumeDefaultVis();

  void SetAllVisible();

  void SetQuarterVisible();

  void SetHalfVisible();

  void SetNoEndVisible();

  void SetPhysicalNode();
  
  //void SetPhysicalDefaultVis();

  void SetDetector();
  
  void SetHits();
  
  void SetVisEmcDetector();

  void SetVisEmcHits();

  int GetPartNb();

  int GetThetaNb(int part);

  int GetPhiNb(int part, int theta);

  void SetVolumeEmc(TGeoVolume *vol) { m_Emc = vol; }
  
  TGeoVolume *GetVolumeEmc() { return m_Emc; }
  
  TGeoVolume *GetVolumePart( int part );

  TGeoVolume *GetVolumePhi( int part, int phi );
  
  TGeoVolume *GetVolumeTheta( int part, int phi, int theta );

  TGeoVolume *GetVolumeCrystal( int part, int phi, int theta );

  TGeoNode *GetPart( int part );

  TGeoNode *GetPhi( int part, int phi );

  TGeoNode *GetTheta( int part, int phi, int theta );
  
  TGeoNode *GetTheta2( int part, int phi, int theta );
  
  TGeoNode *GetCrystal( int part, int phi, int theta );

  TGeoNode *GetCrystal2( int part, int phi, int theta );

  TGeoPhysicalNode *GetPhysicalCrystal( int part, int phi, int theta);
  
  TGeoPhysicalNode *GetPhysicalCrystal2( int part, int phi, int theta);
  
  bool HasTwoNodes( int part, int phi, int theta );

  void ComputeThetaPhi( int id, int sector, int nb, int &CryNumberTheta, int &CryNumberPhi);

  int  ComputeEndCopyNb(int num);

  Emc2DCrystal *Get2DCrystal( Int_t part, Int_t phi, Int_t theta ); 

  Int_t GetPart(TGeoPhysicalNode* phyNode);

  void Draw(Option_t *option);

  void DrawHits(Option_t *option);
 
  void Align(TGeoPhysicalNode *phyNode, TGeoMatrix *newmat=0, TGeoShape *newshape=0, Bool_t check=kFALSE);

  void RestorePhyNode(TGeoPhysicalNode *phyNode, TGeoNode *node);

 private:
  static const int m_kPart     = 3;
  static const int m_kPhiBr    = 120;
  static const int m_kThetaBr  = 44;

  int m_kPhiEc;
  int m_kThetaEc;
  int m_kSectorEc;
  int m_kNbEc;
  int m_kRealNbEc;
  
  int m_EmcColor;
  int m_partColor;
  int m_phiColor;
  int m_thetaColor;
  int m_brCrystalColor;
  int m_ecCrystalColor;

  static const int EMC_TIME_FACTOR   = 1;
  static const int EMC_CHARGE_FACTOR = 1000; 

  TGeoVolume       *m_Emc;
  TGeoNode         *m_NodePart[m_kPart]; 
  TGeoNode         *m_NodePhi[m_kPart][m_kPhiBr]; 
  TGeoNode         *m_NodeTheta[m_kPart][m_kPhiBr][m_kThetaBr]; // m_kPhiBr > m_kPhiEc, m_kThetaBr > m_kThetaEc
  TGeoNode         *m_NodeTheta2[m_kPart][m_kPhiBr][m_kThetaBr]; // contains theta phi has tow nodes 
  TGeoPhysicalNode *m_PhysicalCrystal[m_kPart][m_kPhiBr][m_kThetaBr];
  TGeoPhysicalNode *m_PhysicalCrystal2[m_kPart][m_kPhiBr][m_kThetaBr];

  TObjArray        *m_PhyNodeOrgArray;
  
  // Volume refer to logicalVolume,  only refer to the one exists.
  // Node   refer to physicalVolume, each real theta has a corresponding physical volume, althought they may be the same one.
  // PhysicalNode gives the path combined by physical volume, it's unique for each real theta.

  const TObjArray  *m_EmcDigiCol;

  // 2D
  Emc2DCrystal     *m_Emc2DCrystal[m_kPart][m_kPhiBr][m_kThetaBr];

};

#endif   /* Emc_ROOT_GEO_H */

---