EvtDToPiPi0Etap.hh

Go to the documentation of this file.

#ifndef EVTDTOPIPI0ETAP_HH 
#define EVTDTOPIPI0ETAP_HH
 
#include "EvtGenBase/EvtDecayProb.hh"
#include "TComplex.h"
 
class EvtParticle;
 
class EvtDToPiPi0Etap:public  EvtDecayProb  {
 
public:
    
    EvtDToPiPi0Etap() {}
    virtual ~EvtDToPiPi0Etap();
 
    void getName(std::string& name);
    EvtDecayBase* clone();
 
    void init();
    void initProbMax();
    void decay(EvtParticle *p);
 
private:
 
    void Com_Multi(double a1[2], double a2[2], double res[2]);
    void Com_Divide(double a1[2], double a2[2], double res[2]);
    void calt1(double daug1[4], double daug2[4], double t1[4]);
    void calt2(double daug1[4], double daug2[4], double t2[4][4]);
    void propagatorRBW(double mass2, double mass, double width, double sa, double sb, double sc, double r2, int l, double prop[2]);
    void propagatorGS(double mass2, double mass, double width, double sa, double sb, double sc, double r2, double prop[2]);
    void calEva(double* Ks01, double* Ks02, double* Pic , double *mass1, double *width1, double *amp, double *phase,int* g0, int* spin,int* modetype, int nstates, double & Result);
 
    double SCADot(double a1[4], double a2[4]);
    double barrier(int l, double sa, double sb, double sc, double r, double mass);
    double Barrier(int l, double sa, double sb, double sc, double r2);
    double wid(double mass2, double mass, double sa, double sb, double sc, double r2, int l);
    double widl1(double mass2, double mass, double sa, double sb, double sc, double r2);
    double DDalitz(double P1[4], double P2[4], double P3[4], int Ang, double mass);
 
    int modetype[6];
    double rho[6];
    double phi[6];
    double mass[6];
    double width[6];
    
    double massDp;
    double massD0;
    
    double massEtap;
    double massPip;
    double massPi0;
    
    double mathPi;
    
    double GS1;
    double GS2;
    double GS3;
    double GS4;
    
    int G[4][4];
};
 
#endif