FTFinder.h

Go to the documentation of this file.

#pragma once
 
#include "CLHEP/Vector/ThreeVector.h"
#include "GaudiKernel/Algorithm.h"
#include "HepPDT/ParticleDataTable.hh"
 
#include "MagneticFieldSvc/IBesMagFieldSvc.h"
#include "RawDataProviderSvc/IRawDataProviderSvc.h"
 
#include "FTGeom.h"
#include "FTLayer.h"
#include "FTList.h"
#include "FTSegment.h"
#include "FTSuperLayer.h"
#include "FTTrack.h"
#include "FTWire.h"
#include "MdcParameter.h"
 
#ifndef ENABLE_BACKWARDS_COMPATIBILITY
typedef HepGeom::Point3D<double> HepPoint3D;
#endif
 
using namespace CLHEP;
 
class FTFinder {
public:
  /// returns FTFinder pointer
  // static FTFinder * GetPointer(void);
 
  void setAlgorithmPointer( Algorithm* );
 
public:
  /// Constructors and destructor
  FTFinder();
 
  /// initializer(creates geometry)
  void init();
 
  /// terminator
  void term();
 
  /// begin run function(reads constants)
  void begin_run();
 
  /// track finder core
  void event();
 
public: // Selectors
  /// returns superlayer
  FTSuperLayer* superLayer( int id ) const;
 
  /// returns track list
  FTList<FTTrack*>& tracks();
 
  /// returns event primary vertex
  CLHEP::Hep3Vector vertex() const;
 
  /// returns wire ID for given FTWire object
  int getWireId( FTWire* ) const;
 
  /// convert t to x
  float t2x( const FTLayer* l, const float t ) const;
 
  /// convert x to t
  float x2t( const FTLayer* l, const float x ) const;
 
  /// set bunchtime for MC events
  void setBunchtime( double bunchtime ) { _bunchtime = bunchtime; };
  /// set crude mdc time calibconst (~200 ns) for t0 calculation
  void setT0cal( double t0cal ) { _t0cal = t0cal; };
 
private: // private member functions
  Algorithm* m_algorithm;
 
  /// clear object
  void clear();
 
  /// unpack RAWCDC and create wire-hit
  int updateMdc();
 
  /// get event start time from segment linear fit
  void getTestime();
 
  /// create track list
  void mkTrackList();
 
  /// create 3D track list
  void mkTrack3D();
 
  /// finds event primary vertex
  int VertexFit( int z_flag );
 
  /// finds event primary vertex from 2D tracks
  int VertexFit2D();
 
  /// find vertex closest to IP
  int findBestVertex();
 
  /// corrects event timing after 2nd r-phi fit and returns event timing
  int CorrectEvtTiming();
 
  /// link axial segments to make track
  FTTrack* linkAxialSegments( FTSegment* inner, FTSegment* outer );
 
  // link 2,3,4 axial segments
  void linkAxialSuperLayer234( FTList<FTSegment*>& inner_segments );
 
  // link 9,10 axial segments
  void linkAxialSuperLayer910( FTList<FTSegment*>& outer_segments );
 
  void linkAxialSegments_step( FTList<FTSegment*>& InnerSegments,
                               FTList<FTSegment*>& OuterSegments,
                               FTList<FTSegment*>& Connected, float maxDphi, float maxChi2 );
 
#ifndef OnlineMode
  /// make Mdst_charged/trk/trk_fit table from reconstruted tracks
  void makeMdst();
#endif
 
  /// output tracking results into TDS
  //   by wangdy
  StatusCode makeTds();
 
  // register T0 into TDS by max
  // StatusCode registT0();
 
public: // public data menbers
  int                   eventNo;
  int                   runNo;
  int                   expflag;
  float                 t0Estime;
  float                 tOffSet;
  float                 evtTiming;
  const HepPoint3D      pivot;
  float                 Testime;
  int                   i_rPhiFit;
  int                   tEstFlag;
  FTList<FTList<float>> tEstime;
 
private: // private data members
  // static FTFinder * _tFinder;
  FTGeom* m_geom;
 
  FTList<FTTrack*>   _tracks;
  FTList<FTSegment*> _linkedSegments;
  FTList<FTSegment*> _axialSegCollect;
 
  double _vx;
  double _vy;
  double _vz;
  int    _ExpNo;
  int    _RunNo;
  double _bunchtime;
  double _t0cal;
  int    m_total_trk;
 
  HepPDT::ParticleDataTable* m_particleTable;
  IRawDataProviderSvc*       m_rawDataProviderSvc;
  IBesMagFieldSvc*           m_pmgnIMF;
 
  static MdcParameter* param;
};