#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 116 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]; }
  }
}