PrimaryVertex.h

Go to the documentation of this file.

#ifndef Analysis_PRIMARYVTX_H
#define Analysis_PRIMARYVTX_H
 
#include "CLHEP/Matrix/Vector.h"
#include "EvtRecEvent/EvtRecPrimaryVertex.h"
#include "EvtRecEvent/EvtRecTrack.h"
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/NTuple.h"
#include "TH1D.h"
 
using CLHEP::HepVector;
 
class EvtRecEvent;
 
class PrimaryVertex : public Algorithm {
public:
  PrimaryVertex( const std::string& name, ISvcLocator* pSvcLocator );
  StatusCode initialize();
  StatusCode execute();
  StatusCode finalize();
 
private:
  StatusCode RegisterEvtRecPrimaryVertex( EvtRecPrimaryVertex* aNewEvtRecPrimaryVertex,
                                          MsgStream&           log );
 
  void SelectGoodChargedTracks( SmartDataPtr<EvtRecEvent>&    recEvent,
                                SmartDataPtr<EvtRecTrackCol>& recTrackCol,
                                std::vector<int>& icp, std::vector<int>& icm,
                                std::vector<int>& iGood );
 
private:
  int m_sel_number[15];
 
  int    m_trackNumberCut;
  double m_vz0Cut;
  double m_cosThetaCut;
  int    m_fitMethod;
  double m_globalChisqCut;
  double m_chisqCut;
  int    m_trackIteration;
  int    m_vertexIteration;
  double m_chi2CutforTrkIter;
  int    m_freedomCut;
  double m_chi2CutforSmooth;
  /*
  //particle ID cut
  bool m_useDedx;
  bool m_useTof1;
  bool m_useTof2;
  bool m_useTofE;
  bool m_useTofQ;
  bool m_useEmc;
  bool m_useMuc;
  double m_PidProbCut;*/
 
  int m_output; // switch of output PrimaryVertex.root
 
  // define NTuples
  NTuple::Tuple*       m_tuple1; // Global vertex fit results
  NTuple::Item<double> m_chig;
  NTuple::Item<long>   m_ndofg;
  NTuple::Item<double> m_probg;
  NTuple::Item<double> m_gvx;
  NTuple::Item<double> m_gvy;
  NTuple::Item<double> m_gvz;
 
  NTuple::Tuple*       m_tuple2; // chi-square of smooth
  NTuple::Item<double> m_chis;
  NTuple::Item<double> m_chif;
  NTuple::Item<double> m_probs;
  NTuple::Item<double> m_probf;
 
  NTuple::Tuple*       m_tuple3;
  NTuple::Item<double> m_pull_drho;
  NTuple::Item<double> m_pull_phi;
  NTuple::Item<double> m_pull_kapha;
  NTuple::Item<double> m_pull_dz;
  NTuple::Item<double> m_pull_lamb;
  NTuple::Item<double> m_pull_momentum;
 
  NTuple::Tuple*       m_tuple4; // Kalman vertex fit results
  NTuple::Item<double> m_chik;
  NTuple::Item<long>   m_ndofk;
  NTuple::Item<double> m_probk;
  NTuple::Item<double> m_kvx;
  NTuple::Item<double> m_kvy;
  NTuple::Item<double> m_kvz;
};
#endif