#include <D0ToKSLKK.h>

Public Member Functions
	D0ToKSLKK ()

virtual	~D0ToKSLKK ()

void	init (int Daug0Id, int Uspin)

complex< double >	Amp (vector< double > k0l, vector< double > kp, vector< double > km, int Daug0Id, int Uspin)

Detailed Description

Definition at line 9 of file D0ToKSLKK.h.

Constructor & Destructor Documentation

◆ D0ToKSLKK()

D0ToKSLKK::D0ToKSLKK ( )

inline

Definition at line 12 of file D0ToKSLKK.h.

12{}

◆ ~D0ToKSLKK()

D0ToKSLKK::~D0ToKSLKK ( )

virtual

Definition at line 32 of file D0ToKSLKK.cxx.

32{}

Member Function Documentation

◆ Amp()

complex< double > D0ToKSLKK::Amp	(	vector< double >	k0l,
		vector< double >	kp,
		vector< double >	km,
		int	Daug0Id,
		int	Uspin )

Definition at line 118 of file D0ToKSLKK.cxx.

                                                                                                               {
    
    double P1[4], P2[4], P3[4];
    P1[0] = k0l[3]; P1[1] = k0l[0]; P1[2] = k0l[1]; P1[3] = k0l[2];
    P2[0] = km[3] ; P2[1] = km[0] ; P2[2] = km[1] ; P2[3] = km[2] ;
    P3[0] = kp[3] ; P3[1] = kp[0] ; P3[2] = kp[1] ; P3[3] = kp[2] ;
        
     if(Daug0Id==310) SorL = true; else SorL = false;
   //double value;
   int spin[5]={0,1,0,0,0};
     double PDFD0[2];
     if(SorL){
        int g0[5]={5,1,3,1,1};
            double r0[5] = {3,3,3,3,3};
            double r1[5] = {5,5,5,5,5};
        int nstates=5;
        //calEva(P1, P2, P3, mass, width, rho, phi, g0, spin, modetype, r0, r1, value, 0, nstates,charge,SorL);
            calPDF(P1, P2, P3, mass, width, rho, phi, g0, spin, modetype, r0, r1, 0, nstates, PDFD0);
     }else if((!SorL)&&Uspin==1){
      int g0[5]={5,1,3,1,1};
            double r0[5] = {6.82071036651904,2.89695915812383,3,3,-6.61721570909587};
            double r1[5] = {1.80820371593891,-2.15090641236998,5,5,-0.972853397735864}; 
            int nstates=5;
        //calEva(P1, P2, P3, mass, width, rho, phi, g0, spin, modetype, r0, r1, value, 0, nstates,charge,SorL);
            calPDF(P1, P2, P3, mass, width, rho, phi, g0, spin, modetype, r0, r1, 0, nstates, PDFD0);
     }else{
      int g0[5]={5,1,3,1,1};
            double r0[5] = {3,3,3,3,3};
            double r1[5] = {5,5,5,5,5};
            int nstates=5;
        //calEva(P1, P2, P3, mass, width, rho, phi, g0, spin, modetype, r0, r1, value, 0, nstates,charge,SorL);
            calPDF(P1, P2, P3, mass, width, rho, phi, g0, spin, modetype, r0, r1, 0, nstates, PDFD0);
     }
 
   return complex<double>(PDFD0[0],PDFD0[1]);
}

Referenced by QCMCFilter::findD0Decay().

◆ init()

void D0ToKSLKK::init	(	int	Daug0Id,
		int	Uspin )

Definition at line 34 of file D0ToKSLKK.cxx.

                                          {
 
  //std::cout << "D0ToKSLKK ==> Initialization !" << std::endl;
    
   phi[0] = 0;
   rho[0] = 1;
   phi[1] = 3.42566005625144;
   rho[1] = 0.481320242396153;
   phi[2] = 3.38836604528211;
   rho[2] = 0.692817392287617;
     phi[3] = 2.09747486815378;    
   rho[3] = 0.642052025575875;
   phi[4] = -1.19528927010597;            
   rho[4] = 0.363639226481389;
 
   modetype[0]= 23;
   modetype[1]= 23;
   modetype[2]= 13; 
   modetype[3]= 13;
   modetype[4]= 23;
 
     if(Uspin==1&&Daug0Id==130){
         //std::cout << "D0ToKLKK ==> Initialization !" << std::endl;
         modetype[0]= 232;
        modetype[1]= 232;
        modetype[2]= 13; 
        modetype[3]= 13;
        modetype[4]= 232;
     }
 
   /*for (int i=0; i<5; i++) {
     cout << i << " rho= " << rho[i] << " phi= " << phi[i] << endl;
   }*/
 
   width[0]  = 0.272;
   width[1]  = 0.004249;
   width[2]  = 0.272;
   width[3]  = 0.258;
   width[4]  = 0.258;
 
   mass[0]   = 0.919;
   mass[1]   = 1.019461;
   mass[2]   = 0.919;
   mass[3]   = 1.439;
   mass[4]   = 1.439;
 
     mDM = 1.86484;
     mK0  = 0.497614;
     mKa  = 0.49368;
     mPi  = 0.13957;
     mEta = 0.547862;
     mKa2 = 0.24371994;//0.49368^2;
     mPi2 = 0.01947978;//0.13957^2;
     mEta2 = 0.30015277;//0.547862^2;
     mass_EtaP = 0.95778;
     mass_Kaon = 0.49368;
     mass_KS = 0.4976;
    
     math_pi = 3.1415926;
     mass_Pion2 = 0.0194797849;
     mass_2Pion = 0.27914;
     math_2pi = 6.2831852;
     rD2 = 25.0;    // 5*5
     rRes2 = 9.0;   // 3*3
     g2  = 0.23;             //K*0(1430)
 
 
   GS1 = 0.636619783;
   GS2 = 0.01860182466;
   GS3 = 0.1591549458;   //  1/(2*math_2pi)
   GS4 = 0.00620060822;  //  mass_Pion2/math_pi
 
    rho_omega = 0.00294;
    phi_omega = -0.02; 
 
   int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
   for (int i=0; i<4; i++) {
      for (int j=0; j<4; j++) {
         G[i][j] = GG[i][j];
      }
   }
}