/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Muc/MucRecEvent/MucRecEvent-00-02-52/MucRecEvent/MucRec3DRoad.h

Go to the documentation of this file.

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

/*
 *    2003/12/25   Zhengyun You      Peking University
 * 
 *    2005/02/27   Zhengyun You      Peking University
 *                 transplanted to Gaudi framework
 */

#ifndef MUC_REC_3DROAD_H
#define MUC_REC_3DROAD_H

#include "GaudiKernel/ContainedObject.h"
#include "GaudiKernel/SmartRef.h" 
#include "GaudiKernel/ObjectVector.h"

#include "MucRecEvent/MucRec2DRoad.h"

class MucRec3DRoad : virtual public ContainedObject 
{
public:

  MucRec3DRoad(MucRec2DRoad *road0, MucRec2DRoad *road1);

  MucRec3DRoad &operator=(const MucRec3DRoad &orig);
  
  MucRec3DRoad(const MucRec3DRoad &source);
  ~MucRec3DRoad();

  void SetIndex(const int &index);

  void SetGroup(const int &Group);

  void SetSimpleFitParams(const float &m_SlopeZX, 
                          const float &m_InterceptZX,
                          const float &m_SlopeZY, 
                          const float &m_InterceptZY);
  int RefitNoCurrentGap(const int &gap,
                         float &slopeZX,
                         float &interceptZX,
                         float &slopeZY,
                         float &interceptZY);

  vector<Identifier> ProjectToStrip(const int& part, const int& gap, const HepPoint3D& gPoint, const Hep3Vector&gDirection);
  
  vector<Identifier> ProjectToStrip(const int& part, const int& gap, const HepPoint3D& gPoint, const Hep3Vector&gDirection, vector<int>&padID, vector<float>&intersect_x, vector<float>&intersect_y, vector<float>&intersect_z);
  
  void ProjectWithSigma(const int &gap, float &x, float &y, float &z, float &sigmaX, float &sigmaY, float &sigmaZ);

  void ProjectNoCurrentGap(const int &gap, float &x, float &y, float &z,  float &sigmaX, float &sigmaY, float &sigmaZ,
                           float& quad_x1, float& quad_y1, float& quad_z1, float& quad_x2, float& quad_y2, float& quad_z2);

  void Project(const int &gap, float &x1, float &y1, float &z1, float &x2, float &y2, float &z2);

  void Project(const int &gap, const float zx, const float x0,const float zy,const float y0 , float &x1, float &y1, float &z1);

  //void ProjectToReference(float &x, float &y, float &z);

  //void ProjectToZ(float &x, float &y, float &zplane);

  void TransformPhiRToXY(const float &vk, const float &vr,
                         const float &k0, const float &r0,
                         float &vx, float &vy, 
                         float &x0, float &y0) const;

  float GetVxSign(const float vk) const;

  int GetIndex() const;
  
  int GetPart() const;

  int GetSeg() const;

  int GetGroup() const;
  
  int GetLastGap() const;
  
  int GetNGapsWithHits() const;

  int GetTotalHits() const;

  int GetHitsPerGap(const int &gap) const;

  int GetMaxHitsPerGap() const;

  bool HasHitInGap(const int &gap) const;
  
  int GetNSharedHits(const MucRec3DRoad *road) const;

  // Are the two X and Y clusters in each gap in the same panel?
  // bool ClusterPanelsMatch() const;

  int GetLastGapDelta() const;

  int GetTotalHitsDelta() const;
  
  int GetDegreesOfFreedom() const;

  float GetReducedChiSquare() const;

  float GetSlopeZX() const;

  float GetInterceptZX() const;

  float GetSlopeZY() const;

  float GetInterceptZY() const;

  MucRecHit *GetHit(const int &gap) const;
  
  //bool IsIntersectionOK(const TMuiCluster *hcluster,
  //            const TMuiCluster *vcluster) const;

  MucRec2DRoad *Get2DRoad(const int &orient = 0) const;
  
  vector<Identifier> GetHitsID() const;

  void PrintHitsInfo();

private:

  Hep3Vector m_VertexPos;    // Event vertex position.
  Hep3Vector m_VertexSigma;  // Uncertainty of event vertex position.

  MucRec2DRoad *m_Road0;  // object describing the 2-D road orient 0
  MucRec2DRoad *m_Road1;  // object describing the 2-D road orient 1

  int   m_Index;            // uniquely identifies this road in the event
  int   m_Part;            // the part in which this road was found
  int   m_Seg;             // the segment in which this road was found
  int   m_LastGap;         // Last gap with attached hits.
  float m_Chi2;            // Chi-square parameter of trajectory fit.
  int   m_DOF;             // Number of degrees of freedom of fit.

  float m_SlopeZX;         // X = m_InterceptZX + m_SlopeZX * Z.
  float m_InterceptZX;     // X = m_InterceptZX + m_SlopeZX * Z.
  float m_SlopeZY;         // Y = m_InterceptZY + m_SlopeZY * Z.
  float m_InterceptZY;     // Y = m_InterceptZY + m_SlopeZY * Z.

  // == true if X and Y clusters are in the same panels in all gaps.
  //bool  m_ClusterPanelsMatch;
  int m_Group;            // Which group does this road belong to?
};

typedef ObjectVector<MucRec3DRoad> MucRec3DRoadContainer;

#endif  /* MUC_REC_3DROAD_H */

---