/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Muc/MucGeomSvc/MucGeomSvc-00-02-25/MucGeomSvc/MucGeoGap.h

Go to the documentation of this file.

//$id$
//
//$log$

/*
 *    2003/08/30   Zhengyun You     Peking University 
 * 
 *    2004/09/09   Zhengyun You     Peking University
 *                 transplanted to Gaudi framework
 */

#ifndef MUC_GEO_GAP_H
#define MUC_GEO_GAP_H

#include <iostream>
#include <vector>
#include <CLHEP/Vector/ThreeVector.h>
#include <CLHEP/Geometry/Point3D.h>
#include <CLHEP/Geometry/Plane3D.h>
#include <CLHEP/Vector/Rotation.h>
#include <CLHEP/Matrix/Matrix.h>
#include "TGeoNode.h"
#include "TGeoMatrix.h"
#include "TGeoPhysicalNode.h"

#ifndef ENABLE_BACKWARDS_COMPATIBILITY
typedef HepGeom::Point3D<double> HepPoint3D;
#endif
#ifndef ENABLE_BACKWARDS_COMPATIBILITY
typedef HepGeom::Vector3D<double> HepVector3D;
#endif

using namespace std;
using namespace CLHEP;

class MucGeoStrip;

class MucGeoGap
{
 public:

  MucGeoGap();

  MucGeoGap(const int part,
            const int seg,
            const int gap,
            const int orient,
            const int stripNum,
            const float xSize,
            const float ySize,
            const float zSize,
            const float xTarget1Global,
            const float yTarget1Global,
            const float zTarget1Global,
            const float xTarget2Global,
            const float yTarget2Global,
            const float zTarget2Global,
            const float xTarget3Global,
            const float yTarget3Global,
            const float zTarget3Global,
            const float dzHighEdge,
            const float dzFarFrontGas,
            const float dzNearFrontGas,
            const float dzNearBackGas,
            const float dzFarBackGas,
            const float dxTarget1ToFiducial,
            const float dyTarget1ToFiducial,
            const float dxFiducialToCenter,
            const float dyFiducialToCenter);

  MucGeoGap(const int part,
            const int seg,
            const int gap,
            const int orient,
            const int stripNum,
            const TGeoPhysicalNode *gapPhyscialNode,
            const float ironThickness);
 
  MucGeoGap& operator=(const MucGeoGap& orig);

  MucGeoGap(const MucGeoGap& orig);

  ~MucGeoGap();

  int Part() const { return m_Part; }

  int Seg() const { return m_Seg; }
  
  int Gap() const { return m_Gap; }

  int Orient() const { return m_Orient; }

  // MucSoftID GetSoftID() const { return m_GapID; }

  int GetStripNum() const { return m_StripNum; }

  float GetIronThickness() const { return m_IronThickness; }

  bool GetHitStatus() const { return m_HitStatus; }

  float GetdzFarFrontGas()  const { return m_dzFarFrontGas; }  

  float GetdzNearFrontGas() const { return m_dzNearFrontGas; }  

  float GetdzFarBackGas()   const { return m_dzFarBackGas; }  

  float GetdzNearBackGas()  const { return m_dzNearBackGas; }  

  MucGeoStrip *GetStrip(const int strip) const;

  HepPoint3D GetCenter() const { return m_Center; }

  void GetSize(float &xSize, float &ySize, float &zSize) const;
 
  void GetRotationMatrix(float &thetaX, float &phiX,
                         float &thetaY, float &phiY,
                         float &thetaZ, float &phiZ) const;
  
  int GuessStrip(const float x, const float y,const float z) const;
  
  HepPoint3D ProjectToGap(const HepPoint3D gPoint, 
                          const Hep3Vector gVect) const;

  HepPoint3D ProjectToGapWithSigma(const HepPoint3D gPoint,
              const Hep3Vector gVect,
              const HepPoint3D gPointSigma,
              const Hep3Vector gVectSigma,
              HepPoint3D &gCross,
              HepPoint3D &gCrossSigma) const;

  void ProjectToGapSurface(const HepPoint3D gPoint, 
                                 const Hep3Vector gVect,
                                 HepPoint3D& cross1,
                                 HepPoint3D& cross2) const;
  
  HepPoint3D ProjectToGapQuadLocal(const int part,
                          const int orient,   const float a,
                          const float b,  const float c, const int whichhalf,
                          HepPoint3D& cross1,
                          HepPoint3D& cross2) const;



  HepPoint3D ProjectToGap(const HepPoint3D gPoint, 
                          const float vy,const float y0, const float a,
                          const float b,  const float c, const int whichhalf,
                          HepPoint3D& cross1,
                          HepPoint3D& cross2) const;

  void ProjectToGapSurface(const HepPoint3D gPoint, 
                           const float vy,const float y0, const float a,
                           const float b,  const float c, const int whichhalf,
                           HepPoint3D& cross1,
                           HepPoint3D& cross2) const;
  

  HepPoint3D CompareIntersection(const int whichhalf, const HepPoint3D gCross1, const HepPoint3D gCross2,
                                 const float a, const float b, const float c )const;

  
  Hep3Vector RotateToGlobal(const Hep3Vector pVect) const; 

  Hep3Vector RotateToGap(const Hep3Vector gVect) const;
  
  HepPoint3D TransformToGlobal(const HepPoint3D pPoint) const; 
  
  HepPoint3D TransformToGap(const HepPoint3D gPoint) const;
  
  bool IsInGap(const float x, const float y,const float z) const;

  MucGeoStrip *AddStrip(const int strip);

  //void SetID(const int part, const int seg, const int gap)
  //{ m_GapID.SetID(part, seg, gap); }

  void SetHitStatus(const bool hit) { m_HitStatus = hit; }
  
 private:
   
  int   m_Part;             // Part No of this gap.
  int   m_Seg;              // Segment No of this gap.
  int   m_Gap;              // Gap No of this gap.
  int   m_StripNum;         // Number of strips in this gap.
  int   m_Orient;           // Orientation of this gap.
  bool  m_HitStatus;        // Whether this gap is hit.

  float m_XSize;            // X size of this gap.
  float m_YSize;            // Y size of this gap.
  float m_ZSize;            // Z size of this gap.
  float m_dzHighEdge;       // Z distance from higer Z edge to strip center.
  
  float m_dzFarFrontGas;    // Z distance from low  edge of front gas chamber to strip center. 
  float m_dzNearFrontGas;   // Z distance from high edge of front gas chamber to strip center.
  float m_dzNearBackGas;    // Z distance from low  edge of back  gas chamber to strip center.
  float m_dzFarBackGas;     // Z distance from high edge of back  gas chamber to strip center.

  float m_IronThickness;    // Thickness of iron before this gap 
  
  HepPoint3D  m_Center;     // Position of gap center in global coordinate.
  HepRotation m_Rotation;   // Rotaion matrix from global to gap coordinate.
  HepRotation m_RotationT;  // Rotaion matrix from gap to global coordinate.
  Hep3Vector  m_Translation;// Transform matrix from global to gap coordinate.

  HepPoint3D  m_SurfaceInner;     // Position of the inner surface center in this gap in global coordinate.
  HepPoint3D  m_SurfaceOuter;     // Position of the outer surface center in this gap in global coordinate.

  TGeoPhysicalNode *m_PhyNode;

  vector<MucGeoStrip*> m_pMucGeoStrip; 
};
 
ostream& operator << (ostream& s, const MucGeoGap& gap);

#endif    /* MUC_GEO_GAP_H */

---