QCMCFilter.h

Go to the documentation of this file.

#ifndef QCMCFILTER_h
#define QCMCFILTER_h 
 
#include <string>
#include "HepPDT/ParticleDataTable.hh"
#include "GaudiKernel/Algorithm.h"
 
#include "EventModel/EventHeader.h"
#include "EvtRecEvent/EvtRecEvent.h"
 
#include "CLHEP/Matrix/SymMatrix.h"
#include <vector>
using namespace std;
 
class IDataProviderSvc;
 
class QCMCFilter : public Algorithm {
 public:
  QCMCFilter(const std::string& name, ISvcLocator* pSvcLocator);
 
  std::vector<double> findD0Decay(int charm);
 
  StatusCode initialize();
  StatusCode execute();
  StatusCode finalize();
 
  // ------------ constants, enums and typedefs --------------
 
 private:
  IDataProviderSvc* m_evtSvc;
 
  HepPDT::ParticleDataTable* m_particleTable;
    
    bool m_debug;
  // ------------ data members -------------------------------
  //int m_nUnknownEvents ;
  //int m_nUnknownDecays ;
  //int m_nD0D0barEvents ;
  //int m_nD0bar ;
  //int m_nDpDmEvents ;
  //int m_nD0D0barDiscarded ;
  //int m_nDpDmDiscarded ;
  //int m_nCPPlus ;
  //int m_nCPMinus ;
  //int m_nFlavoredCFD0_0 ;
  //int m_nFlavoredCFD0_1 ;
  //int m_nFlavoredCFD0_3 ;
  //int m_nFlavoredCSD0 ;
  //int m_nFlavoredDCSD0_0 ;
  //int m_nFlavoredDCSD0_1 ;
  //int m_nFlavoredDCSD0_3 ;
  //int m_nSL ;
  
  //int m_nDalitz ;
  //double m_dalitzNumer1 ;
  //double m_dalitzNumer2 ;
  //double m_dalitzDenom ;
  
  bool m_invertSelection ;
  double m_x ;
  double m_y ;
  double m_rCF ;
  double m_RCF ;
  double m_zCF ;
  bool m_wCFSign ; // true = +, false = -
  double m_rCF_0 ;
  double m_dCF_0 ;
  double m_RCF_0 ;
  double m_zCF_0 ;
  bool m_wCFSign_0 ; // true = +, false = -
  double m_rCF_1 ;
  double m_dCF_1 ;
  double m_RCF_1 ;
  double m_zCF_1 ;
  bool m_wCFSign_1 ; // true = +, false = -
  double m_rCF_3 ;
  double m_RCF_3 ;
  double m_dCF_3 ;
  double m_zCF_3 ;
  bool m_wCFSign_3 ; // true = +, false = -
  double m_rCS ;
  double m_zCS ;
  bool m_wCSSign ; // true = +, false = -
  double m_dplusDminusWeight ;
  int m_dalitzModel ; 
  bool m_useNewWeights ;
  
  //double m_rwsCF ;
  //double m_rwsCF_0 ;
  //double m_rwsCF_1 ;
  //double m_rwsCF_3 ;
  //double m_rwsCS ;
  //double m_deltaCF ;
  //double m_deltaCF_0 ;
  //double m_deltaCF_1 ;
  //double m_deltaCF_3 ;
  //double m_deltaCS ;
  
  double m_largestWeight ;
 
    //vector<HepLorentzVector> m_D0daughter_P4;
    //vector<int> m_D0daughter_PID;
 
};
 
#endif