QCMCFilter.h

Go to the documentation of this file.

#ifndef QCMCFILTER_h
#define QCMCFILTER_h
 
#include "GaudiKernel/Algorithm.h"
#include "HepPDT/ParticleDataTable.hh"
#include <string>
 
#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_FCP_PiPiPi0;
 
  // 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