FTTrack.h

Go to the documentation of this file.

#pragma once
 
#include "TrackUtil/Helix.h"
#include "TrackUtil/Lpav.h"
#include "TrackUtil/zav.h"
 
#include "MagneticFieldSvc/IBesMagFieldSvc.h"
 
// #include "FTSegment.h"
#include "FTList.h"
#include "MdcParameter.h"
 
typedef HepGeom::Point3D<double> HepPoint3D;
 
class FTWire;
class FTSegment;
class FTFinder;
 
class FTTrack {
public:
  /// constructor
  FTTrack( FTList<FTSegment*>& axial_segments, float kappa, float chi2_kappa );
 
  /// destructor
  ~FTTrack();
 
public: // public member functions
  /// do r-phi circle fit
  int r_phiFit();
 
  /// do r-phi refit
  int r_phiReFit( float vx, float xy, int vtx_flag );
 
  // cut the largest residual hit when do r-phi fit
  int r_phi2Fit( float vx, float xy, int vtx_flag );
  int r_phi3Fit( int l, float vx, float xy, int vtx_flag );
  int r_phi4Fit( float vx, float xy, int vtx_flag );
 
  /// do s-z linear fit
  int s_zFit();
 
  /// update s and z information for linking
  void updateSZ();
 
  /// link stereo segments by tanLambda
  int linkStereoSegments();
 
public: // Selectors
  /// returns helix parameters
  Helix* helix() const;
 
  /// returns lpav
  const Lpav& lpav() const;
 
  /// returns zav
  const zav& Zav() const;
 
  /// returns axial segments
  const FTList<FTSegment*>& axial_segments();
 
  /// returns stereo_segments
  const FTList<FTSegment*>& stereo_segments();
 
  /// returns kappa at linking
  float kappa_tmp() const;
 
  /// returns sigmaKappa at linking
  float chi2_kappa_tmp() const;
 
  /// returns difference between z and estimated z
  /// at stereo segment linking if its valid
  float d_z( float s, float z ) const;
 
#ifndef OnlineMode
  /// printout _tracks infomation
  /// added by X.-R. Lu
  void printout();
#endif
 
public: // Modifiers
  /// add z for culculation of tanLambda
  float SigmaZ( float z );
 
  /// add s for culculation of tanLambda
  float SigmaS( float s );
 
  /// add s for culculation of dz, tanLambda
  float SigmaSS( float ss );
 
  /// add s for culculation of dz, tanLambda
  float SigmaSZ( float sz );
 
  /// append stereo segment to the stereo segment list
  void append_stereo( FTSegment*, float s, float z );
 
  /// append stereo segment to the cache
  void append_stereo_cache( FTSegment* );
 
  /// calculate the wire hits number
  int get_nhits();
 
public: // static data members
  float minPt;
  float minDr;
  float xtCoEff;
  float Testime;
 
public:
  void setFTFinder( FTFinder* );
 
private:
  FTFinder*        m_ftFinder;
  IBesMagFieldSvc* m_pmgnIMF;
 
public: // private data members
  Lpav*                      _la;
  zav*                       _za;
  Helix*                     _helix;
  FTList<FTSegment*>         _axial_segments;
  FTList<FTSegment*>         _stereo_segments;
  FTList<FTSegment*>         _stereo_segments_cache;
  FTList<FTList<FTSegment*>> _stereo_segments_by_superLayer;
  static MdcParameter*       param;
 
  union {
    float _kappa;
    float _SigmaS;
  };
 
  union {
    float _chi2_kappa;
    float _SigmaSS;
  };
 
  float _SigmaZ;
  float _SigmaSZ;
};