NeutralDReconstruction.h

Go to the documentation of this file.

#ifndef Reconstruction_NeutralDReconstruction_H
#define Reconstruction_NeutralDReconstruction_H
 
#include "GaudiKernel/Algorithm.h"
 
#include "EvtRecEvent/EvtRecDTag.h"
#include "ReadBeamInfFromDb.h"
 
#include "BesDChain/CDChargedKaonList.h"
#include "BesDChain/CDChargedPionList.h"
#include "BesDChain/CDDecayList.h"
#include "BesDChain/CDEta.h"
#include "BesDChain/CDKs.h"
#include "BesDChain/CDPhoton.h"
#include "BesDChain/CDPi0.h"
 
#include "VertexDbSvc/IVertexDbSvc.h"
 
#include <string>
 
#include "ITools.h"
 
class NeutralDReconstruction : public Algorithm {
 
public:
  NeutralDReconstruction( const std::string& name, ISvcLocator* pSvcLocator );
  StatusCode initialize();
  StatusCode execute();
  StatusCode finalize();
 
  void saveD0Info( CDDecayList::iterator, double, int, int, EvtRecDTag* );
  void updateKsInfo( CDDecayList::iterator, double, int, int, EvtRecDTag*, vector<int>,
                     IVertexDbSvc*, bool );
  void savetrack( vector<int>, vector<int>, EvtRecTrackIterator, EvtRecTrackIterator,
                  EvtRecTrackIterator, EvtRecTrackIterator, EvtRecDTag* );
  void pidtag( vector<int>, vector<int>, CDChargedKaonList&, CDChargedPionList&, EvtRecDTag* );
  // HepLorentzVector vfit(string channel, vector<int> kaonid, vector<int> pionid, HepPoint3D
  // vx, EvtRecTrackIterator charged_begin);
  vector<string> getlist( string& filename );
 
private:
  StatusCode registerEvtRecDTagCol( EvtRecDTagCol* dtagCol, MsgStream& log );
 
private:
  IPiKPSelector<CDChargedPion>* m_pionSelector{ nullptr };
  IPiKPSelector<CDChargedKaon>* m_kaonSelector{ nullptr };
  ISelectorTool<CDPhoton>*      m_photonSelector{ nullptr };
  ISelectorTool<CDKs>*          m_ksSelector{ nullptr };
  ISelectorTool<CDPi0>*         m_pi0Selector{ nullptr };
  ISelectorTool<CDEta>*         m_etatoGGSelector{ nullptr };
  IEbeamBetaSelector*           m_neutralDSelector{ nullptr };
  ISelectorTool<CDDecay>*       m_eptoPiPiEtaSelector{ nullptr };
  ISelectorTool<CDDecay>*       m_rhotoPiPiSelector{ nullptr };
  ISelectorTool<CDDecay>*       m_eptoRhoGamSelector{ nullptr };
 
  bool m_debug;
 
  int m_irun;
  int m_ievt;
  int m_nChrg;
  int m_nNeu;
  int m_nPion;
  int m_nKaon;
  int m_nPi0;
  int m_nKs;
 
  string         m_decaylist;
  vector<string> chanlist;
 
  double            m_beamE;
  Hep3Vector        m_beta;
  bool              m_ReadBeamEFromDB;
  bool              m_usecalibBeamE;
  bool              m_usevertexfit;
  ReadBeamInfFromDb m_readDb;
 
  bool m_useVFrefine;
  bool m_useBFC;
};
#endif