Changeset View
Changeset View
Standalone View
Standalone View
EvtGenBase/EvtDalitzPlot.hh
| Show First 20 Lines • Show All 46 Lines • ▼ Show 20 Lines | public: | ||||
| // Absolute limits for Zemach coordinate qres and qhel (approximate) | // Absolute limits for Zemach coordinate qres and qhel (approximate) | ||||
| // qHelAbsMin(BC,CA) means absolute minimum for (qCA-qAB)/2. | // qHelAbsMin(BC,CA) means absolute minimum for (qCA-qAB)/2. | ||||
| double qResAbsMin( EvtCyclic3::Pair i ) const; | double qResAbsMin( EvtCyclic3::Pair i ) const; | ||||
| double qResAbsMax( EvtCyclic3::Pair i ) const; | double qResAbsMax( EvtCyclic3::Pair i ) const; | ||||
| double qHelAbsMin( EvtCyclic3::Pair i ) const; | double qHelAbsMin( EvtCyclic3::Pair i ) const; | ||||
| double qHelAbsMax( EvtCyclic3::Pair i ) const; | double qHelAbsMax( EvtCyclic3::Pair i ) const; | ||||
| inline double qSumMin() const { return sum() + _ldel; } | inline double qSumMin() const { return sum() + m_ldel; } | ||||
| inline double qSumMax() const { return sum() + _rdel; } | inline double qSumMax() const { return sum() + m_rdel; } | ||||
| inline bool fuzzy() const { return ( _rdel - _ldel != 0. ); } | inline bool fuzzy() const { return ( m_rdel - m_ldel != 0. ); } | ||||
| // Find the area of the Dalitz plot by numeric integration. (N bins for variable q(i) are used). | // Find the area of the Dalitz plot by numeric integration. (N bins for variable q(i) are used). | ||||
| // Very large numbers of N can result in a very long calculation. It should not | // Very large numbers of N can result in a very long calculation. It should not | ||||
| // matter which two pairs f variables are used. The integral should eventually | // matter which two pairs f variables are used. The integral should eventually | ||||
| // converge to the same number | // converge to the same number | ||||
| double getArea( int N = 1000, EvtCyclic3::Pair i = EvtCyclic3::AB, | double getArea( int N = 1000, EvtCyclic3::Pair i = EvtCyclic3::AB, | ||||
| EvtCyclic3::Pair j = EvtCyclic3::BC ) const; | EvtCyclic3::Pair j = EvtCyclic3::BC ) const; | ||||
| Show All 21 Lines | public: | ||||
| // and birth vertices | // and birth vertices | ||||
| EvtTwoBodyVertex vD( EvtCyclic3::Pair iRes, double m0, int L ) const; | EvtTwoBodyVertex vD( EvtCyclic3::Pair iRes, double m0, int L ) const; | ||||
| EvtTwoBodyVertex vB( EvtCyclic3::Pair iRes, double m0, int L ) const; | EvtTwoBodyVertex vB( EvtCyclic3::Pair iRes, double m0, int L ) const; | ||||
| // Accessors | // Accessors | ||||
| double sum() const; | double sum() const; | ||||
| inline double bigM() const { return _bigM; } | inline double bigM() const { return m_bigM; } | ||||
| inline double mA() const { return _mA; } | inline double mA() const { return m_mA; } | ||||
| inline double mB() const { return _mB; } | inline double mB() const { return m_mB; } | ||||
| inline double mC() const { return _mC; } | inline double mC() const { return m_mC; } | ||||
| double m( EvtCyclic3::Index i ) const; | double m( EvtCyclic3::Index i ) const; | ||||
| void print() const; | void print() const; | ||||
| void sanityCheck() const; | void sanityCheck() const; | ||||
| protected: | protected: | ||||
| // Defines two dimensional dalitz plot | // Defines two dimensional dalitz plot | ||||
| double _mA; | double m_mA; | ||||
| double _mB; | double m_mB; | ||||
| double _mC; | double m_mC; | ||||
| double _bigM; | double m_bigM; | ||||
| // Defines third dimension, or fuzziness. M^2 + ldel < M^2 < M^2 + rdel | // Defines third dimension, or fuzziness. M^2 + ldel < M^2 < M^2 + rdel | ||||
| double _ldel; | double m_ldel; | ||||
| double _rdel; | double m_rdel; | ||||
| }; | }; | ||||
| #endif | #endif | ||||