#include "FullAmp.h" #include "Resonance.h" #include "Particle.h" #include #include #ifndef _INTENSITY_ #define _INTENSITY_ class Intensity{ public: ///Constructor & Dectructor Intensity(); /* Intensity(int J, std::vector Mj, std::vector lya1, int s1, std::vector lya2, int s2, */ /* std::vector l, int s, double mass1, double mass2, double massR, double widthR); */ Intensity(std::vector Rf1, std::vector Rf2, std::vector Rf3, Particle pi1, Particle pi2, Particle pi3); virtual ~Intensity(){}; /* ///Calculation intensity in helisity formalizm */ /* std::complex Calc(double phi, double theta, double invmass); */ ///Calculation intensity in helisity formalizm (as coherent sum for 12 and 23 decay) std::complex Calc(double phi13, double theta13, double invmass12, double phi21, double theta21, double invmass23, double phi12, double theta12, double invmass13); ///Debug printing void PrintPar(); private: /* double _J;//total moment of decay particle */ /* std::vector _Mj;// moment projection of decay particle */ double _s1; std::vector _lya1; // spin and helicity of 1st particle double _s2; std::vector _lya2;// spin and helicity of 2nd particle double _s3; std::vector _lya3;// spin and helicity of 2nd particle std::vector _lya;//helicity of decay particle // std::vector _l;//angular moment int _s;// spin double _mass1; // mass of 1st particle double _mass2; // mass of 2nd particle double _mass3; // mass of 3nd particle // double _massR; // mass of resonance // double _widthR;// width of resonance std::vector _Rf1;//resonances for (1,2) pair std::vector _Rf2;//resonances for (2,3) pair std::vector _Rf3;//resonances for (1,3) pair int _maxJr1n;//position of max spin for resonance(s) for (1,2) pair int _maxJr2n;//position of max spin for resonance(s) for (2,3) pair int _maxJr3n;//position max spin for resonance(s) for (1,3) pair Particle _pi1;// particle #1 Particle _pi2;// particle #2 Particle _pi3;// particle #3 FullAmp _a1; }; #endif