Reconstruction/MdcPatRec/MdcxReco/include/MdcxReco/MdcxHit.h

Go to the documentation of this file.

//--------------------------------------------------------------------------
// File and Version Information:
//  $Id: MdcxHit.h,v 1.10 2009/12/16 08:59:13 zhangy Exp $
//
// Description:
//  Class Header for |MdcxHit|: hit that can calculate derivatives
//      and plot hits self specialized to the drift chamber
//
// Environment:
//  Software developed for the BaBar Detector at the SLAC B-Factory.
//
// Author List:
//  A. Snyder
//  Zhang Yao(zhangyao@ihep.ac.cn) Migration for BESIII
//
// Copyright Information:
//  Copyright (C) 1995  BEPCII
//
// History:
//  Migration for BESIII MDC
//
//------------------------------------------------------------------------
 
#ifndef _MDCXHIT_
#define _MDCXHIT_
 
#include "MdcCalibFunSvc/IMdcCalibFunSvc.h"
#include "MdcData/MdcHit.h"
#include "MdcGeom/MdcDetector.h"
#include "MdcGeom/MdcLayer.h"
#include "MdcGeom/MdcSWire.h"
#include <fstream>
#include <iostream>
#include <math.h>
#include <vector>
 
class MdcxHel;
class MdcDigi;
 
using namespace EventModel;
 
class MdcxHit {
  // A MDC hit that can calculate derivatives
public:
  // constructors
  MdcxHit( const MdcDigi* pdcdatum, float c0 = 0, float cresol = .0180 );
  MdcxHit( const MdcHit* pdchhit, float c0 = 0, float cresol = .0180 );
 
  // destructor
  virtual ~MdcxHit(){};
  // private copy constructor and oprator
private:
  MdcxHit( const MdcxHit& );
  MdcxHit& operator=( const MdcxHit& );
  bool     operator==( const MdcxHit& ) const;
 
public:
  static void setMdcCalibFunSvc( const IMdcCalibFunSvc* calibSvc );
  static void setCountPropTime( bool countPropTime ); // yzhang for csmc
  static void setMdcDetector( const MdcDetector* gm );
 
  // accessors
  inline const MdcDigi* getDigi() const { return _mdcDigi; }
  inline const MdcHit*  getMdcHit() const { return _mdcHit; }
  inline int            SuperLayer() const { return _superlayer; } // SuperLayer#
  inline int            WireNo() const { return _wirenumber; }     // Wire#
  inline int            Layer() const { return _layernumber; }     // layer#
  inline float          t() const { return _t; }                   // drift time
  inline float          x() const { return _x; }                   // x of wire
  inline float          y() const { return _y; }                   // y of wire
  inline float          xpos() const { return _xpos; }
  inline float          ypos() const { return _ypos; }
  inline float          xneg() const { return _xneg; }
  inline float          yneg() const { return _yneg; }
  inline float          wx() const { return _wx; }
  inline float          wy() const { return _wy; }
  inline float          wz() const { return _wz; }
  inline float          pw() const { return _pw; }
  inline float          phiMid() const { return _pw; } // yzhang add
  inline float          stereo() const { return _s; }  // stereo angle of wire
  inline float          v() const { return _v; }       // drift velocity
  inline int            type() const { return _type; } // flags |MdcxHit|
 
  // workers
  inline void SetConstErr( int i ) { _consterr = i; }
  ///////////////////////////////////////////
  float tcor( float zh = 0.0, float tof = 0.0, float tzero = 0.0 ) const;
  float d( MdcxHel& hel ) const; // |drift dist| (changes hel's internal state)
  float d( float zh = 0.0, float tof = 0.0, float tzero = 0.0, int wamb = 0,
           float eang = 0.0 ) const; // |drift dist|
  // float tcor(float tof) const;
  // float d(int ambig = 0, double entranceAngle = 0) const;
  // float d(MdcxHel& hel, int ambig = 0, double entranceAngle = 0) const;
  ///////////////////////////
 
  float              e( float dd = 0.0 ) const;                 // drift error currently in use
  float              pull( MdcxHel& hel ) const;                // Chisq contribution to fit
  float              residual( MdcxHel& hel ) const;            // residual of this hit
  std::vector<float> derivatives( MdcxHel& hel ) const;         // Derivatives, etc.
  void               print( std::ostream& o, int i = 0 ) const; // print this hit
  void               printAll( std::ostream& o, int i = 0 ) const; // print this hit
  inline void        SetUsedOnHel( int i ) { usedonhel = i; }
  inline int         GetUsedOnHel() const { return usedonhel; }
 
protected:
  // functions
  void process();
 
  // data
  unsigned       _wirenumber;
  unsigned       _layernumber;
  unsigned       _superlayer;
  float          _t;
  float          _q;
  unsigned       _iTdc;
  unsigned       _iAdc;
  float          _T0Walk; // zoujh
  float          _x;
  float          _y;
  float          _r;
  float          _s;
  float          _L; // length of this layer
  float          _d;
  float          _v;
  int            _type; // yzhang add 2009-11-03
  float          _e;
  float          _xpos, _ypos, _xneg, _yneg;
  float          _p, _sp, _cp;
  double         _pw;
  double         _wx, _wy, _wz;
  const MdcHit*  _mdcHit;
  const MdcDigi* _mdcDigi; // pointer to |MdcDigi| defining this hit
  int            _consterr;
  int            usedonhel;
  float          _c0; // accumulated global time offset; changes in t0 go here
  float          _cresol;
 
private:
  double getSigma( float driftDist, int ambig = 0, double entrAngle = 0, double dipAngle = 0,
                   double z = 0 ) const;
  static const IMdcCalibFunSvc* m_mdcCalibFunSvc; // pointer to MdcCalibFunSvc
  static const MdcDetector*     m_gm;             // pointer to MdcDetector
  static bool                   m_countPropTime;
}; // endof MdcxHit
 
typedef ObjectVector<MdcxHit> MdcxHitCol;
 
#endif