#include "QCMCFilterAlg/QCMCFilter.h"
#include "QCMCFilterAlg/Dalitz.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/AlgFactory.h"
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/PropertyMgr.h"
#include "GaudiKernel/Bootstrap.h"
#include "GaudiKernel/RegistryEntry.h"
#include "TMath.h"
#include <cmath>
#include "HepPDT/ParticleDataTable.hh"
#include "HepPDT/ParticleData.hh"
#include "GaudiKernel/IPartPropSvc.h"
#include "McTruth/McParticle.h"
#include "McTruth/MdcMcHit.h"
#include "EventModel/EventModel.h"
#include "EventModel/Event.h"
#include "EventModel/EventHeader.h"
#include "EvtRecEvent/EvtRecEvent.h"
#include "CLHEP/Random/RandFlat.h"
#include "CLHEP/Matrix/Vector.h"
#include "CLHEP/Matrix/Matrix.h"
#include "CLHEP/Matrix/SymMatrix.h"
#include "CLHEP/Vector/ThreeVector.h"
#include "CLHEP/Vector/LorentzVector.h"
#include "CLHEP/Vector/TwoVector.h"
#include <vector>
Go to the source code of this file.
Functions | |
HepSymMatrix | m_weights (10, 0) |
HepMatrix | m_modeCounter (10, 10, 0) |
HepMatrix | m_keptModeCounter (10, 10, 0) |
Variables | |
const double | xmpi0 = 0.134976 |
const double | xmeta = 0.54784 |
const double | xmkaon = 0.49368 |
const double | xmpion = 0.13957 |
const double | xmk0 = 0.49761 |
const double | xmrho = 0.77549 |
const double | xmd0 = 1.86484 |
const double | PI = 3.1415926 |
static const int | kPsi3770ID = 30443 |
static const int | kD0ID = 421 |
static const int | kD0BarID = -421 |
static const int | kDpID = 411 |
static const int | kDmID = -411 |
static const int | kGammaID = 22 |
static const int | kGammaFSRID = -22 |
static const int | kPiPlusID = 211 |
static const int | kPiMinusID = -211 |
static const int | kPi0ID = 111 |
static const int | kEtaID = 221 |
static const int | kEtaPrimeID = 331 |
static const int | kF0ID = 9010221 |
static const int | kFPrime0ID = 10221 |
static const int | kF0m1500ID = 50221 |
static const int | kF2ID = 225 |
static const int | kA00ID = 10111 |
static const int | kA0PlusID = 10211 |
static const int | kA0MinusID = -10211 |
static const int | kRhoPlusID = 213 |
static const int | kRhoMinusID = -213 |
static const int | kRho0ID = 113 |
static const int | kRho2SPlusID = 30213 |
static const int | kRho2SMinusID = -30213 |
static const int | kRho2S0ID = 30113 |
static const int | kA1PlusID = 20213 |
static const int | kA1MinusID = -20213 |
static const int | kA10ID = 20113 |
static const int | kOmegaID = 223 |
static const int | kPhiID = 333 |
static const int | kKPlusID = 321 |
static const int | kKMinusID = -321 |
static const int | kK0SID = 310 |
static const int | kK0LID = 130 |
static const int | kK0ID = 311 |
static const int | kK0BarID = -311 |
static const int | kKStarPlusID = 323 |
static const int | kKStarMinusID = -323 |
static const int | kKStar0ID = 313 |
static const int | kKStar0BarID = -313 |
static const int | kK0Star0ID = 10311 |
static const int | kK0Star0BarID = -10311 |
static const int | kK0StarPlusID = 10321 |
static const int | kK0StarMinusID = -10321 |
static const int | kK1PlusID = 10323 |
static const int | kK1MinusID = -10323 |
static const int | kK10ID = 10313 |
static const int | kK10BarID = -10313 |
static const int | kK1PrimePlusID = 20323 |
static const int | kK1PrimeMinusID = -20323 |
static const int | kK1Prime0ID = 20313 |
static const int | kK1Prime0BarID = -20313 |
static const int | kK2StarPlusID = 325 |
static const int | kK2StarMinusID = -325 |
static const int | kK2Star0ID = 315 |
static const int | kK2Star0BarID = -315 |
static const int | kEMinusID = 11 |
static const int | kEPlusID = -11 |
static const int | kMuMinusID = 13 |
static const int | kMuPlusID = -13 |
static const int | kNuEID = 12 |
static const int | kNuEBarID = -12 |
static const int | kNuMuID = 14 |
static const int | kNuMuBarID = -14 |
static const int | kFlavoredCF = 0 |
static const int | kFlavoredCFBar = 1 |
static const int | kFlavoredCS = 2 |
static const int | kFlavoredCSBar = 3 |
static const int | kSLPlus = 4 |
static const int | kSLMinus = 5 |
static const int | kCPPlus = 6 |
static const int | kCPMinus = 7 |
static const int | kDalitz = 8 |
static const int | kNDecayTypes = 9 |
int | m_nUnknownEvents = 0 |
int | m_nUnknownDecays = 0 |
int | m_nD0D0barEvents = 0 |
int | m_nD0bar = 0 |
int | m_nDpDmEvents = 0 |
int | m_nD0D0barDiscarded = 0 |
int | m_nDpDmDiscarded = 0 |
int | m_nCPPlus = 0 |
int | m_nCPMinus = 0 |
int | m_nFlavoredCFD0 = 0 |
int | m_nFlavoredCSD0 = 0 |
int | m_nFlavoredDCSD0 = 0 |
int | m_nSL = 0 |
int | m_nDalitz = 0 |
double | m_dalitzNumer1 = 0 |
double | m_dalitzNumer2 = 0 |
double | m_dalitzDenom = 0 |
double | dalitzNumer1_fil = 0 |
double | dalitzNumer2_fil = 0 |
double | dalitzDenom_fil = 0 |
double | m_rwsCF = 0. |
double | m_rwsCS = 0. |
double | m_deltaCF = 0. |
double | m_deltaCS = 0. |
int | m_rho_flag |
HepMatrix m_keptModeCounter | ( | 10 | , | |
10 | , | |||
0 | ||||
) |
Referenced by QCMCFilter::execute(), and QCMCFilter::finalize().
HepMatrix m_modeCounter | ( | 10 | , | |
10 | , | |||
0 | ||||
) |
Referenced by QCMCFilter::execute(), and QCMCFilter::finalize().
HepSymMatrix m_weights | ( | 10 | , | |
0 | ||||
) |
double dalitzDenom_fil = 0 |
Definition at line 164 of file QCMCFilter.cxx.
Referenced by QCMCFilter::execute(), and QCMCFilter::finalize().
double dalitzNumer1_fil = 0 |
Definition at line 162 of file QCMCFilter.cxx.
Referenced by QCMCFilter::execute(), and QCMCFilter::finalize().
double dalitzNumer2_fil = 0 |
Definition at line 163 of file QCMCFilter.cxx.
Referenced by QCMCFilter::execute(), and QCMCFilter::finalize().
const int kA00ID = 10111 [static] |
const int kA0MinusID = -10211 [static] |
const int kA0PlusID = 10211 [static] |
const int kA10ID = 20113 [static] |
Definition at line 92 of file QCMCFilter.cxx.
const int kA1MinusID = -20213 [static] |
const int kA1PlusID = 20213 [static] |
const int kCPMinus = 7 [static] |
Definition at line 138 of file QCMCFilter.cxx.
Referenced by QCMCFilter::findD0Decay(), and QCMCFilter::initialize().
const int kCPPlus = 6 [static] |
Definition at line 137 of file QCMCFilter.cxx.
Referenced by QCMCFilter::findD0Decay(), and QCMCFilter::initialize().
const int kD0BarID = -421 [static] |
Definition at line 67 of file QCMCFilter.cxx.
Referenced by QCMCFilter::execute(), and QCMCFilter::findD0Decay().
const int kD0ID = 421 [static] |
Definition at line 66 of file QCMCFilter.cxx.
Referenced by QCMCFilter::execute(), and QCMCFilter::findD0Decay().
const int kDalitz = 8 [static] |
Definition at line 139 of file QCMCFilter.cxx.
Referenced by QCMCFilter::execute(), QCMCFilter::findD0Decay(), and QCMCFilter::initialize().
const int kDmID = -411 [static] |
Definition at line 69 of file QCMCFilter.cxx.
const int kDpID = 411 [static] |
const int kEMinusID = 11 [static] |
const int kEPlusID = -11 [static] |
const int kEtaID = 221 [static] |
const int kEtaPrimeID = 331 [static] |
const int kF0ID = 9010221 [static] |
const int kF0m1500ID = 50221 [static] |
const int kF2ID = 225 [static] |
const int kFlavoredCF = 0 [static] |
Definition at line 131 of file QCMCFilter.cxx.
Referenced by QCMCFilter::findD0Decay(), and QCMCFilter::initialize().
const int kFlavoredCFBar = 1 [static] |
Definition at line 132 of file QCMCFilter.cxx.
Referenced by QCMCFilter::findD0Decay(), and QCMCFilter::initialize().
const int kFlavoredCS = 2 [static] |
Definition at line 133 of file QCMCFilter.cxx.
Referenced by QCMCFilter::findD0Decay(), and QCMCFilter::initialize().
const int kFlavoredCSBar = 3 [static] |
Definition at line 134 of file QCMCFilter.cxx.
Referenced by QCMCFilter::findD0Decay(), and QCMCFilter::initialize().
const int kFPrime0ID = 10221 [static] |
const int kGammaFSRID = -22 [static] |
const int kGammaID = 22 [static] |
const int kK0BarID = -311 [static] |
const int kK0ID = 311 [static] |
const int kK0LID = 130 [static] |
const int kK0SID = 310 [static] |
const int kK0Star0BarID = -10311 [static] |
const int kK0Star0ID = 10311 [static] |
const int kK0StarMinusID = -10321 [static] |
const int kK0StarPlusID = 10321 [static] |
const int kK10BarID = -10313 [static] |
const int kK10ID = 10313 [static] |
const int kK1MinusID = -10323 [static] |
const int kK1PlusID = 10323 [static] |
const int kK1Prime0BarID = -20313 [static] |
const int kK1Prime0ID = 20313 [static] |
const int kK1PrimeMinusID = -20323 [static] |
const int kK1PrimePlusID = 20323 [static] |
const int kK2Star0BarID = -315 [static] |
Definition at line 120 of file QCMCFilter.cxx.
const int kK2Star0ID = 315 [static] |
Definition at line 119 of file QCMCFilter.cxx.
const int kK2StarMinusID = -325 [static] |
const int kK2StarPlusID = 325 [static] |
const int kKMinusID = -321 [static] |
const int kKPlusID = 321 [static] |
const int kKStar0BarID = -313 [static] |
const int kKStar0ID = 313 [static] |
const int kKStarMinusID = -323 [static] |
const int kKStarPlusID = 323 [static] |
const int kMuMinusID = 13 [static] |
const int kMuPlusID = -13 [static] |
const int kNDecayTypes = 9 [static] |
const int kNuEBarID = -12 [static] |
const int kNuEID = 12 [static] |
const int kNuMuBarID = -14 [static] |
const int kNuMuID = 14 [static] |
const int kOmegaID = 223 [static] |
const int kPhiID = 333 [static] |
const int kPi0ID = 111 [static] |
const int kPiMinusID = -211 [static] |
const int kPiPlusID = 211 [static] |
const int kPsi3770ID = 30443 [static] |
const int kRho0ID = 113 [static] |
const int kRho2S0ID = 30113 [static] |
const int kRho2SMinusID = -30213 [static] |
const int kRho2SPlusID = 30213 [static] |
const int kRhoMinusID = -213 [static] |
const int kRhoPlusID = 213 [static] |
const int kSLMinus = 5 [static] |
Definition at line 136 of file QCMCFilter.cxx.
Referenced by QCMCFilter::findD0Decay(), and QCMCFilter::initialize().
const int kSLPlus = 4 [static] |
Definition at line 135 of file QCMCFilter.cxx.
Referenced by QCMCFilter::findD0Decay(), and QCMCFilter::initialize().
double m_dalitzDenom = 0 |
Definition at line 161 of file QCMCFilter.cxx.
double m_dalitzNumer1 = 0 |
Definition at line 159 of file QCMCFilter.cxx.
double m_dalitzNumer2 = 0 |
Definition at line 160 of file QCMCFilter.cxx.
double m_deltaCF = 0. |
Definition at line 168 of file QCMCFilter.cxx.
double m_deltaCS = 0. |
Definition at line 169 of file QCMCFilter.cxx.
int m_nCPMinus = 0 |
Definition at line 153 of file QCMCFilter.cxx.
int m_nCPPlus = 0 |
Definition at line 152 of file QCMCFilter.cxx.
int m_nD0bar = 0 |
Definition at line 148 of file QCMCFilter.cxx.
int m_nD0D0barDiscarded = 0 |
Definition at line 150 of file QCMCFilter.cxx.
int m_nD0D0barEvents = 0 |
Definition at line 147 of file QCMCFilter.cxx.
int m_nDalitz = 0 |
Definition at line 158 of file QCMCFilter.cxx.
int m_nDpDmDiscarded = 0 |
Definition at line 151 of file QCMCFilter.cxx.
int m_nDpDmEvents = 0 |
Definition at line 149 of file QCMCFilter.cxx.
int m_nFlavoredCFD0 = 0 |
Definition at line 154 of file QCMCFilter.cxx.
int m_nFlavoredCSD0 = 0 |
Definition at line 155 of file QCMCFilter.cxx.
int m_nFlavoredDCSD0 = 0 |
Definition at line 156 of file QCMCFilter.cxx.
int m_nSL = 0 |
Definition at line 157 of file QCMCFilter.cxx.
int m_nUnknownDecays = 0 |
Definition at line 146 of file QCMCFilter.cxx.
int m_nUnknownEvents = 0 |
Definition at line 145 of file QCMCFilter.cxx.
int m_rho_flag |
double m_rwsCF = 0. |
Definition at line 166 of file QCMCFilter.cxx.
double m_rwsCS = 0. |
Definition at line 167 of file QCMCFilter.cxx.
const double PI = 3.1415926 |
Definition at line 61 of file QCMCFilter.cxx.
const double xmd0 = 1.86484 |
Definition at line 59 of file QCMCFilter.cxx.
const double xmeta = 0.54784 |
Definition at line 54 of file QCMCFilter.cxx.
const double xmk0 = 0.49761 |
const double xmkaon = 0.49368 |
Definition at line 55 of file QCMCFilter.cxx.
const double xmpi0 = 0.134976 |
Definition at line 53 of file QCMCFilter.cxx.
const double xmpion = 0.13957 |
const double xmrho = 0.77549 |
Definition at line 58 of file QCMCFilter.cxx.