EdbTrackFitter.h

Go to the documentation of this file.

#ifndef ROOT_EdbTrackFitter
#define ROOT_EdbTrackFitter
 
//                                                                      //
// EdbTrackFitter                                                       //
//                                                                      //
// definitions and functions for track fit                              //
//                                                                      //
#include "EdbPattern.h"
 
class EdbScanCond;
 
//______________________________________________________________________________
class EdbTrackFitter : public TNamed {
 
 private:
  int   eNsegMax;        // max number of segments (for arrays allocation)
 
 public:
  float eX0;             // rad length of the media [microns]
  float eM;              // mass of the particle (if negative - use the mass setted in the track)
  float ePdef;           // default momentum
  float eTPb;            // ?
  float ePcut;           // minimal momentum?
  bool  eDE_correction;  // take into account the energy loss or not
 
 public:
  EdbTrackFitter();
  virtual ~EdbTrackFitter(){}
 
  void           SetNsegMax(int nseg) {eNsegMax=nseg;}
  void           SetDefaultBrick();
  int            FitTrackLine(EdbTrackP &tr);
  int            FitTrackLine(const EdbTrackP &tr, float &x,float &y,float &z,float &tx,float &ty,float &w);
  int            Fit3Pos(EdbTrackP &tr);
 
  static float   MaxKink(EdbTrackP &tr);
  static float   MeanKink(EdbTrackP &tr);
  static float   Theta( EdbSegP &s, EdbSegP &s1 );
  static float   Chi2Seg( EdbSegP *s1, EdbSegP *s2);
  static float   MaxChi2Seg(EdbTrackP &tr);
  static float   MeanChi2Seg(EdbTrackP &tr);
  bool           SplitTrack( EdbTrackP &t, EdbTrackP &t1, int isplit );
  int            SplitTrackByKink( EdbTrackP *t, TObjArray &tracks, float maxkink );
  float          PMS_KF(EdbTrackP &t, float p0=10., float probbest=0.5);
 
  static float   Chi2ACP( EdbSegP s1, EdbSegP s2, EdbScanCond &cond);
  static float   Chi2ASeg( EdbSegP &s1, EdbSegP &s2, EdbSegP &s, EdbScanCond &cond1, EdbScanCond &cond2);
  static float   Chi2PSeg( EdbSegP &s1, EdbSegP &s2, EdbSegP &seg, EdbScanCond &cond1, EdbScanCond &cond2);
  static float   Chi2ASegLL( EdbSegP &s1, EdbSegP &s2, EdbSegP &s, EdbScanCond &cond1, EdbScanCond &cond2);
  float          Chi2SegM( EdbSegP s1, EdbSegP s2, EdbSegP &s, EdbScanCond &cond1, EdbScanCond &cond2);
  float          Chi2SegMCS( const EdbSegP &s1, const EdbSegP &s2);
  double         ProbSegMCS(EdbSegP *s1, EdbSegP *s2);
  void           Print();
  
  static int RMSprojXY(EdbTrackP &tr, float &ex, float &ey );
 
  ClassDef(EdbTrackFitter,1)  // track fitter
};
 
#endif /* ROOT_EdbTrackFitter */