MdcTrkRecon.h

Go to the documentation of this file.

#ifndef MDCTRKRECON_H
#define MDCTRKRECON_H
 
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/NTuple.h"
#include "MagneticFieldSvc/IBesMagFieldSvc.h"
#include "MdcPrintSvc/IMdcPrintSvc.h"
#include "MdcRecEvent/RecMdcTrack.h"
#include "MdcTrkRecon/MdcFlagHold.h"
#include "RawDataProviderSvc/IRawDataProviderSvc.h"
#include <string>
 
class MdcDetector;
class MdcSegData;
class MdcSegList;
class MdcTrackListBase;
class MdcSegFinder;
namespace MdcPatRec {
  class BField;
}
class IHistogram1D;
class McParticle;
#ifdef MDCPATREC_TIMETEST
class IBesTimerSvc;
class BesTimer;
#endif
 
#ifdef UseCalib
class MdcCalibFunSvc;
#endif
 
class MdcTrkRecon : public Algorithm {
public:
  MdcTrkRecon( const std::string& name, ISvcLocator* pSvcLocator );
  ~MdcTrkRecon();
  StatusCode initialize();
  StatusCode execute();
  StatusCode finalize();
  StatusCode beginRun();
 
protected:
  bool       ignoringUsedHits() const { return m_onlyUnusedHits; }
  bool       poisoningHits() const { return m_poisonHits; }
  void       fillSegList();
  void       fillTrackList();
  void       dumpDigi();
  void       dumpTdsTrack( RecMdcTrackCol* trackList );
  void       fillMcTruth();
  void       fillEvent();
  StatusCode bookNTuple();
 
private:
  bool m_beginRun{ false };
 
  int         t_iExecute;
  int         m_fittingMethod;
  int         m_allHit; // if all hit (active and inactive) include in the number of nHit
  std::string m_configFile;
  std::string m_paramFile;
  std::string m_pdtFile;
  const MdcDetector*              _gm;
  MdcFlagHold                     m_flags; // run-time-settable flags and params
  std::auto_ptr<MdcSegData>       m_hitData;
  std::auto_ptr<MdcSegList>       m_segs;
  std::auto_ptr<MdcTrackListBase> m_tracks;
  std::auto_ptr<MdcSegFinder>     m_segFinder;
  IRawDataProviderSvc*            m_rawDataProviderSvc;
  IMdcPrintSvc*                   m_mdcPrintSvc;
 
  MdcPatRec::BField* m_bfield;
 
  std::vector<int>   m_selEvtNo;
  bool               m_onlyUnusedHits;
  bool               m_poisonHits;
  bool               m_doLineFit;
  bool               m_tryBunch;
  bool               m_recForEsTime;
  std::vector<float> m_helixHitsSigma;
 
  uint32_t m_getDigiFlag;
  int      m_maxMdcDigi;
  bool     m_keepBadTdc;
  bool     m_dropHot;
  bool     m_keepUnmatch;
  int      m_minMdcDigi;
 
  bool   m_fieldCosmic;
  double m_d0Cut;
  double m_z0Cut;
  double m_dropTrkPt;
  bool   m_mcHist;
  int    m_hist;
  int    m_debug;
  bool   m_doSagFlag;
  bool   m_arbitrateHits;
  bool   m_combineTracking;
 
  double mcDrift[43][288]; // FIXME assume max 43 track/event, 288 hit/track
  double mcLayer[43][288];
  double mcWire[43][288];
  double mcLR[43][288];
  double mcX[43][288];
  double mcY[43][288];
  double mcZ[43][288];
  double hitOnSeg[43];
  double hitOnMcTk[43];
  int    hitPoisoned[43][288];
  int    hitInSegList[43][288];
  int    t_nHitInTk[100];
  int    t_t0Stat;
  double t_t0;
  double t_t0Truth;
  double t_mcTkNum;
  int    t_nDigi;
  int    t_eventNo;
  int    t_nRecTk;
  int    t_iExexute;
  int    isPrimaryOfMcTk[100];
  int    pdgOfMcTk[100];
  bool   m_noInner;
#ifdef MDCPATREC_TIMETEST
  IBesTimerSvc* m_timersvc;
  BesTimer*     m_timer[2];
  double        ti_nHit;
  double        ti_p;
  double        ti_pt;
#endif
  IBesMagFieldSvc* m_pIMF;
};
 
#endif