ExtMucKal.h

Go to the documentation of this file.

#ifndef ExtMucKal_H_
#define ExtMucKal_H_
 
// #include "CLHEP/Alist/AList.h"
#include "CLHEP/Matrix/SymMatrix.h"
#include "CLHEP/Matrix/Vector.h"
// #ifndef CLHEP_POINT3D_H
#include "CLHEP/Geometry/Point3D.h"
#include "Ext_xp_err.h"
#include "G4LogicalVolume.hh"
#include "G4VPhysicalVolume.hh"
#include "MucRawEvent/MucDigi.h"
#include "MucRecEvent/MucRecHit.h"
// #endif
 
class ExtMucKal //: public HepAList<KLMKCluster>
{
 
public:
  ExtMucKal();
  ~ExtMucKal();
  bool ExtMucFilter();
  void SetPosMomErr( Hep3Vector pos, Hep3Vector mom, HepSymMatrix err ) {
    m_CurrentPosition = pos;
    m_CurrentMomentum = mom;
    m_CurrentXPErr    = err;
  }
  void               SetMucDigiCol( MucDigiCol* amucdigi ) { m_MucDigiCol = amucdigi; }
  void               SetMucWindow( int aMucWindow ) { n_sigm = aMucWindow; }
  void               SetGapID( Hep3Vector id ) { m_gapid = id; }
  bool               MucKalIniti();
  void               XPmod( Hep3Vector& pos, Hep3Vector& mom, HepSymMatrix& err );
  bool               JCB();
  double             Fit( Hep3Vector& pos, Hep3Vector& mom, HepSymMatrix& err );
  HepMatrix          GetRoationMatrix( MucGeoGap* box );
  bool               GetFilterStatus() { return FilterOK; }
  double             GetChi2() { return Chi2_sub; }
  double             GetDistance( const MucRecHit* hit );
  double             GetPull() { return m_pull; }
  int                GetOrient() { return m_orient; }
  Hep3Vector         GetHitGap();
  vector<MucRecHit*> GapHit();
  vector<MucRecHit*> TrackHit();
  bool               GetSameStrip() { return m_samestrip; }
 
protected:
private:
  Hep3Vector   m_gapid;
  Hep3Vector   m_hitgap;
  Hep3Vector   m_CurrentPosition;
  Hep3Vector   m_CurrentMomentum;
  HepSymMatrix m_CurrentXPErr;
  HepVector3D  m_CurrentInsct;
  HepSymMatrix m_CurrentInsctXPErr;
  HepMatrix    m_jcb;
  MucDigiCol*  m_MucDigiCol;
  bool         FilterOK;
  double       Chi2_sub;
  HepSymMatrix m_err_mod;
  Hep3Vector   m_pos_mod;
  Hep3Vector   m_mom_mod;
  MucRecHit*   m_nearesthit;
  bool         HitExist;
  double       m_pull;
  double       m_sigma;
  int          m_orient;
  int          n_sigm;
  HepVector    m_bm;
  HepSymMatrix m_Ebm;
  int          m_iStrip;
  bool         m_samestrip;
};
 
#endif