include/EvtGenModels/EvtD0ToKspipi.hh

0001 #ifndef EVTD0TOKSPIPI_HH
0002 #define EVTD0TOKSPIPI_HH
0003
0004 #include "EvtGenBase/EvtComplex.hh"
0005 #include "EvtGenBase/EvtDalitzPoint.hh"
0006 #include "EvtGenBase/EvtDalitzReso.hh"
0007 #include "EvtGenBase/EvtDecayAmp.hh"
0008
0009 #include <string>
0010 #include <utility>
0011 #include <vector>
0012
0013 class EvtParticle;
0014
0015 class EvtD0ToKspipi : public EvtDecayAmp {
0016   public:
0017     std::string getName() override;
0018     EvtDecayBase* clone() override;
0019
0020     void init() override;
0021     void initProbMax() override;
0022     void decay( EvtParticle* parent ) override;
0023
0024   private:
0025     // Calculate the total amplitude given the Dalitz plot point
0026     EvtComplex calcTotAmp( const EvtDalitzPoint& point ) const;
0027
0028     // Set particle IDs and PDG masses
0029     void setPDGValues();
0030
0031     // Setup the Dalitz plot resonances and their amplitude coefficients
0032     void initResonances();
0033
0034     // Daughter IDs (updated according to decay file ordering)
0035     int m_d0 = 0;
0036     int m_d1 = 1;
0037     int m_d2 = 2;
0038
0039     // Resonance lineshape and complex amplitude coefficient pair
0040     typedef std::pair<EvtDalitzReso, EvtComplex> ResAmpPair;
0041
0042     // Vector of (resonance, coeff) pairs
0043     std::vector<ResAmpPair> m_resonances;
0044
0045     // IDs of the relevant particles
0046     EvtId m_BP;
0047     EvtId m_BM;
0048     EvtId m_B0;
0049     EvtId m_B0B;
0050     EvtId m_D0;
0051     EvtId m_D0B;
0052     EvtId m_KM;
0053     EvtId m_KP;
0054     EvtId m_K0;
0055     EvtId m_K0B;
0056     EvtId m_KL;
0057     EvtId m_KS;
0058     EvtId m_PIM;
0059     EvtId m_PIP;
0060
0061     // Flavor of the B parent and D
0062     EvtId m_bFlavor;
0063     EvtId m_dFlavor;
0064
0065     // Masses of the relevant particles
0066     double m_mD0;
0067     double m_mKs;
0068     double m_mPi;
0069     double m_mK;
0070 };
0071 #endif