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

Go to the documentation of this file.

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

/*
 *    2003/12/21   Zhengyun You      Peking University
 * 
 *    2005/02/24   Zhengyun You      Peking University
 *                 transplanted to Gaudi framework
 */

#ifndef MUC_REC_2DROAD_H
#define MUC_REC_2DROAD_H

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

#include "MucRecEvent/MucRecHit.h"

using CLHEP::Hep3Vector;       // forward declaration

class MucRec2DRoad : virtual public ContainedObject  
{
public:

  MucRec2DRoad(const int &part = 0, 
               const int &seg  = 0, 
               const int &orient = 0,
               const float &xVertex = 0.0,
               const float &yVertex = 0.0,
               const float &zVertex = 0.0,
               const int &fittingMethod = 0);
  
  MucRec2DRoad &operator=(const MucRec2DRoad &orig);
  
  MucRec2DRoad(const MucRec2DRoad &source);

  ~MucRec2DRoad();
  
  void SetIndex(const int &index);

  void AttachHit(MucRecHit *hit);

  void AttachHitNoFit(MucRecHit *hit);

  void SetMaxNSkippedGaps(const int &nGaps);

  int SimpleFit(float &slope,      float &intercept,
                float &sigmaSlope, float &sigmaIntercept,
                float &chisq,      int &ndof);

  int SimpleFitNoCurrentGap(int currentgap,   float &slope,      float &intercept,
    float &sigmaSlope, float &sigmaIntercept,
    float &chisq,      int &ndof);

  //  This is a helper function for the Project methods.
  //  Give an estimated position (x,y,z) of the point to which we will
  //  project.
  int Fit(const float &x,
          const float &y,
          const float &z,
          float &slope,      float &intercept,
          float &sigmaSlope, float &sigmaIntercept,
          float &chisq,      int &ndof);

  void Project(const int &gap, float &x, float &y, float &z, float &x2, float &y2, float &z2);

  float WeightFunc(const float &chisq, const float &distance) const; 

  float WindowFunc(const float &chisq, const float &distance) const; 

  int GetIndex() const;

  int GetPart() const;

  int GetSeg() const;
  
  int GetOrient() const;

  void GetVertexPos(float &x, float &y, float &z) 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 MucRec2DRoad *road) const;

  float GetSlope() const;

  float GetIntercept() const;

  int GetDegreesOfFreedom() const;

  float GetReducedChiSquare() const;

  void GetSimpleFitParams(float &slope,      float &intercept,
                          float &sigmaSlope, float &sigmaIntercept,
                          float &chisq,      int &ndof) const;

  void GetSimpleFitParams(float &a,      float &b,  float &c, int &whichhalf,
                          float &sigmaa, float &sigmab, float &sigmac,
                          float &chisq,      int &ndof) const;
  bool GetQuadFitOk(){return m_QuadFitOK;}

  void GetPosSigma(float &possigma)const;

  MucRecHit *GetHit(const int &gap) const;
  
  float GetHitDistance(const MucRecHit *hit);
  
  float GetSearchWindowSize(const int &gap) const;

  bool HasHit(MucRecHit *hit) const;

  vector<Identifier> GetHitsID() const;

  vector<MucRecHit*> GetHits() const;

  void PrintHitsInfo() const;

private:

  void CountHits();

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

  int   m_Index;                 // uniquely identifies this road in the event
  int   m_Part;                  // Part in which the road was found.
  int   m_Seg;                   // Segment in which the road was found.
  int   m_Orient;                // orientation of hits attached to this road.

  float m_SimpleSlope;           // slope of fit with simple weights
  float m_SimpleIntercept;       // intercept of fit with simple weights
  float m_SimpleSlopeSigma;      // uncert. in slope of fit with simple weights
  float m_SimpleInterceptSigma;  // uncert. in intercept of fit with simple weights
  float m_SimplePosSigma;        // uncert. in pos of this 2D line with first gap of MUC

  
  float m_SimpleChi2;            // Chi-square of fit with simple weights

  float m_SimpleQuad_a;          //y = a*x*x + b*x +c
  float m_SimpleQuad_b;     
  float m_SimpleQuad_c;
  int   m_SimpleQuad_whichhalf;
  float m_SimpleQuad_aSigma;          //y = a*x*x + b*x +c
  float m_SimpleQuad_bSigma;     
  float m_SimpleQuad_cSigma;     
  int   m_fittingMethod;


  float m_Chi2;                  // Chi-square parameter of trajectory fit.

  int   m_SimpleDOF;             // degrees of freedom of fit with simple weights
  int   m_DOF;                   // Number of degrees of freedom of fit.

  int   m_MaxHitsPerGap;         // Max. number of hits attached in a single gap.
  int   m_LastGap;               // Last gap with attached clusters.
  int   m_TotalHits;             // Total number of hits attached.
  int   m_MaxNSkippedGaps;       // Max consecutive gaps allowed with no hits.

  bool  m_SimpleFitOK;           // If == true, we have a good "simple" fit.
  bool  m_FitOK;                 // If == true, we have a good fit.
  bool  m_QuadFitOK;

  vector<float> m_HitDistance;   // distance from best hit to road
  vector<MucRecHit*> m_pHits;    // Pointers to attached hits.
};

typedef ObjectVector<MucRec2DRoad> MucRec2DRoadContainer;

#endif  /* MUC_REC_2DROAD_H */

---