EvtDsToKSKpPipPim Class Reference

#include <EvtDsToKSKpPipPim.hh>

Inheritance diagram for EvtDsToKSKpPipPim:

Public Member Functions
	EvtDsToKSKpPipPim ()
virtual	~EvtDsToKSKpPipPim ()
void	getName (std::string &name)
EvtDecayBase *	clone ()
void	init ()
void	initProbMax ()
void	decay (EvtParticle *p)
Public Member Functions inherited from EvtDecayProb
void	makeDecay (EvtParticle *p)
void	setProb (double prob)
double	getProb ()
void	setWeight (double weight)
virtual	~EvtDecayProb ()
Public Member Functions inherited from EvtDecayBase
virtual std::string	commandName ()
virtual void	command (std::string cmd)
double	getProbMax (double prob)
double	resetProbMax (double prob)
	EvtDecayBase ()
virtual	~EvtDecayBase ()
virtual bool	matchingDecay (const EvtDecayBase &other) const
EvtId	getParentId ()
double	getBranchingFraction ()
void	disableCheckQ ()
void	checkQ ()
int	getNDaug ()
EvtId *	getDaugs ()
EvtId	getDaug (int i)
int	getNArg ()
int	getPHOTOS ()
void	setPHOTOS ()
void	setVerbose ()
void	setSummary ()
double *	getArgs ()
std::string *	getArgsStr ()
double	getArg (int j)
std::string	getArgStr (int j)
std::string	getModelName ()
int	getDSum ()
int	summary ()
int	verbose ()
void	saveDecayInfo (EvtId ipar, int ndaug, EvtId *daug, int narg, std::vector< std::string > &args, std::string name, double brfr)
void	printSummary ()
void	setProbMax (double prbmx)
void	noProbMax ()
void	checkNArg (int a1, int a2=-1, int a3=-1, int a4=-1)
void	checkNDaug (int d1, int d2=-1)
void	checkSpinParent (EvtSpinType::spintype sp)
void	checkSpinDaughter (int d1, EvtSpinType::spintype sp)
virtual int	nRealDaughters ()

Additional Inherited Members
Static Public Member Functions inherited from EvtDecayBase
static void	findMasses (EvtParticle *p, int ndaugs, EvtId daugs[10], double masses[10])
static void	findMass (EvtParticle *p)
static double	findMaxMass (EvtParticle *p)
Protected Member Functions inherited from EvtDecayBase
bool	daugsDecayedByParentModel ()
Protected Attributes inherited from EvtDecayBase
bool	_daugsDecayedByParentModel

Detailed Description

Definition at line 8 of file EvtDsToKSKpPipPim.hh.

Constructor & Destructor Documentation

EvtDsToKSKpPipPim()

EvtDsToKSKpPipPim::EvtDsToKSKpPipPim ( )

inline

Definition at line 11 of file EvtDsToKSKpPipPim.hh.

{}

Referenced by clone().

~EvtDsToKSKpPipPim()

EvtDsToKSKpPipPim::~EvtDsToKSKpPipPim ( )

virtual

Definition at line 33 of file EvtDsToKSKpPipPim.cc.

{}

Member Function Documentation

clone()

EvtDecayBase * EvtDsToKSKpPipPim::clone ( )

virtual

Implements EvtDecayBase.

Definition at line 37 of file EvtDsToKSKpPipPim.cc.

{ return new EvtDsToKSKpPipPim; }

decay()

void EvtDsToKSKpPipPim::decay ( EvtParticle * p )

virtual

Implements EvtDecayBase.

Definition at line 150 of file EvtDsToKSKpPipPim.cc.

                                              {
  /*
                   double maxprob=0;
                   for(int ir=0;ir<=60000000;ir++){
                   p->initializePhaseSpace(getNDaug(),getDaugs());
                   EvtVector4R _ks   = p->getDaug(0)->getP4();
                   EvtVector4R _kp   = p->getDaug(1)->getP4();
                   EvtVector4R _pip  = p->getDaug(2)->getP4();
                   EvtVector4R _pim = p->getDaug(3)->getP4();
 
                   double _Pip[4],_Pim[4],_Ks[4],_Kp[4];
                   _Ks[0] = _ks.get(0); _Kp[0] = _kp.get(0); _Pip[0] = _pip.get(0); _Pim[0] =
     _pim.get(0); _Ks[1] = _ks.get(1); _Kp[1] = _kp.get(1); _Pip[1] = _pip.get(1); _Pim[1] =
     _pim.get(1); _Ks[2] = _ks.get(2); _Kp[2] = _kp.get(2); _Pip[2] = _pip.get(2); _Pim[2] =
     _pim.get(2); _Ks[3] = _ks.get(3); _Kp[3] = _kp.get(3); _Pip[3] = _pip.get(3); _Pim[3] =
     _pim.get(3);
 
                   double _prob;
                   int g0[8]={1,1,1,1,1,1,1,1};
                   int nstates=10;
 
                   calEvaMy(_Ks, _Kp, _Pip, _Pim, mass, width, rho, phi, g0, modetype, nstates,
     _prob); if(_prob>maxprob) { maxprob=_prob; printf("ir = %d,maxprob = %.10f,prob =
     %.10f\n\n",ir,maxprob,_prob);
                   }
                   }
                   printf("maxprob = %.10f\n", maxprob);
  */
  p->initializePhaseSpace( getNDaug(), getDaugs() );
  EvtVector4R ks  = p->getDaug( 0 )->getP4();
  EvtVector4R kp  = p->getDaug( 1 )->getP4();
  EvtVector4R pip = p->getDaug( 2 )->getP4();
  EvtVector4R pim = p->getDaug( 3 )->getP4();
 
  double Pip[4], Pim[4], Kp[4], Ks[4];
  Ks[0]  = ks.get( 0 );
  Kp[0]  = kp.get( 0 );
  Pip[0] = pip.get( 0 );
  Pim[0] = pim.get( 0 );
  Ks[1]  = ks.get( 1 );
  Kp[1]  = kp.get( 1 );
  Pip[1] = pip.get( 1 );
  Pim[1] = pim.get( 1 );
  Ks[2]  = ks.get( 2 );
  Kp[2]  = kp.get( 2 );
  Pip[2] = pip.get( 2 );
  Pim[2] = pim.get( 2 );
  Ks[3]  = ks.get( 3 );
  Kp[3]  = kp.get( 3 );
  Pip[3] = pip.get( 3 );
  Pim[3] = pim.get( 3 );
 
  // Ks[0]   = 0.6096125554;   Ks[1]   =  0.2036145752;   Ks[2]   = -0.2256821642;   Ks[3]   =
  // 0.1884652565; Kp[0]   = 0.7106732604;   Kp[1]   = -0.2599714599;   Kp[2]   = 0.3643877399;
  // Kp[3]   = 0.2469330230; Pip[0]  = 0.4177584141;   Pip[1]  = -0.0562004660;   Pip[2]  =
  // 0.2821155804;    Pip[3]  =-0.0020089709; Pim[0] = 0.3945120988;   Pim[1] = -0.1225632042;
  // Pim[2]  = -0.3645928242;   Pim[3]  = 0.0521817527;
 
  // Ks[0] = 0.6439438863;  Ks[1] = 0.3368405824; Ks[2] =  0.2308441601; Ks[3] = 0.0171298406;
  // Kp[0] = 0.7146391007;  Kp[1] = 0.4061832461; Kp[2] = -0.1841669581; Kp[3] =  0.2609401579;
  // Pip[0] = 0.2817222681; Pip[1] = -0.1005161765;Pip[2] = -0.1172908044;Pip[3] =
  // -0.1898077097; Pim[0] = 0.3638598699;Pim[1] = -0.3106229764;Pim[2] = 0.0905020525;Pim[3] =
  // 0.0907574503;
 
  double value;
  int    g0[8]   = { 1, 1, 1, 1, 1, 1, 1, 1 };
  int    nstates = 8;
 
  calEvaMy( Ks, Kp, Pip, Pim, mass, width, rho, phi, g0, modetype, nstates, value );
  setProb( value );
 
  return;
}

getName()

void EvtDsToKSKpPipPim::getName ( std::string & name )

virtual

Implements EvtDecayBase.

Definition at line 35 of file EvtDsToKSKpPipPim.cc.

{ model_name = "DsToKSKpPipPim"; }

init()

void EvtDsToKSKpPipPim::init ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 39 of file EvtDsToKSKpPipPim.cc.

                             {
  // check that there are 0 arguments
  checkNArg( 0 );
  checkNDaug( 4 );
 
  checkSpinParent( EvtSpinType::SCALAR );
  checkSpinDaughter( 0, EvtSpinType::SCALAR );
  checkSpinDaughter( 1, EvtSpinType::SCALAR );
  checkSpinDaughter( 2, EvtSpinType::SCALAR );
  checkSpinDaughter( 3, EvtSpinType::SCALAR );
 
  int mode = 0;
  phi[1]   = -3.1173e-03;
  phi[2]   = 5.6119e+00;
  phi[3]   = -7.0591e-01;
  phi[4]   = -3.1220e+00;
  phi[5]   = 3.5914e+00;
  phi[6]   = 3.8510e+00;
  phi[7]   = -1.0830e+00;
 
  rho[1]      = 1.8516e+00;
  rho[2]      = 2.8585e+00;
  rho[3]      = 1.1065e+00;
  rho[4]      = 5.4084e-01;
  rho[5]      = 1.3816e+00;
  rho[6]      = 1.9177e+00;
  rho[7]      = 2.2119e+00;
  modetype[0] = 7;
  modetype[1] = 4;
  modetype[2] = 3;
  modetype[3] = 15;
  modetype[4] = 8;
  modetype[5] = 11;
  modetype[6] = 12;
  modetype[7] = 17;
  phi[0]      = 0;
  rho[0]      = 1;
 
  // cout << "DsToKSKpPipPim :" << endl;
  // for (int i=0; i<8; i++) {
  //    cout << i << " rho= " << rho[i] << " phi= " << phi[i] << endl;
  // }
 
  mD         = 1.86486;
  metap      = 0.95778;
  mk0        = 0.497614;
  mass_Kaon  = 0.49368;
  mass_Eta   = 0.547862;
  mass_Pion  = 0.13957;
  math_pi    = 3.1415926;
  mass_Pion2 = 0.0194797849;
  mass_2Pion = 0.27914;
  math_2pi   = 6.2831852;
 
  rD2      = 25.0;
  rRes2    = 9.0;
  GK1270   = 8.7000e-02;
  GK1400   = 1.7400e-01;
  GKst0    = 4.7300e-02;
  GKstp    = 5.0300e-02;
  GP1680   = 1.5000e-01;
  GRho     = 1.4780e-01;
  Ga0_980  = 5.0000e-02;
  Geta1475 = 9.0000e-02;
  mK1270   = 1.2720e+00;
  mK1400   = 1.4030e+00;
  mKst0    = 8.9555e-01;
  mKstp    = 8.9176e-01;
  mP1680   = 1.6800e+00;
  mRho     = 7.7526e-01;
  ma0_980  = 9.9000e-01;
  meta1475 = 1.4750e+00;
 
  GS1 = 0.636619783;
  GS2 = 0.01860182466;
  GS3 = 0.1591549458;  //  1/(2*math_2pi)
  GS4 = 0.00620060822; //  mass_Pion2/math_pi
 
  int GG[4][4]       = { { 1, 0, 0, 0 }, { 0, -1, 0, 0 }, { 0, 0, -1, 0 }, { 0, 0, 0, -1 } };
  int EE[4][4][4][4] = {
      { { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
        { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 1 }, { 0, 0, -1, 0 } },
        { { 0, 0, 0, 0 }, { 0, 0, 0, -1 }, { 0, 0, 0, 0 }, { 0, 1, 0, 0 } },
        { { 0, 0, 0, 0 }, { 0, 0, 1, 0 }, { 0, -1, 0, 0 }, { 0, 0, 0, 0 } } },
      { { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, -1 }, { 0, 0, 1, 0 } },
        { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
        { { 0, 0, 0, 1 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { -1, 0, 0, 0 } },
        { { 0, 0, -1, 0 }, { 0, 0, 0, 0 }, { 1, 0, 0, 0 }, { 0, 0, 0, 0 } } },
      { { { 0, 0, 0, 0 }, { 0, 0, 0, 1 }, { 0, 0, 0, 0 }, { 0, -1, 0, 0 } },
        { { 0, 0, 0, -1 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 1, 0, 0, 0 } },
        { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
        { { 0, 1, 0, 0 }, { -1, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } } },
      { { { 0, 0, 0, 0 }, { 0, 0, -1, 0 }, { 0, 1, 0, 0 }, { 0, 0, 0, 0 } },
        { { 0, 0, 1, 0 }, { 0, 0, 0, 0 }, { -1, 0, 0, 0 }, { 0, 0, 0, 0 } },
        { { 0, -1, 0, 0 }, { 1, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
        { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } } } };
  for ( int i = 0; i < 4; i++ )
  {
    for ( int j = 0; j < 4; j++ )
    {
      G[i][j] = GG[i][j];
      for ( int k = 0; k < 4; k++ )
      {
        for ( int l = 0; l < 4; l++ ) { E[i][j][k][l] = EE[i][j][k][l]; }
      }
    }
  }
}

initProbMax()

void EvtDsToKSKpPipPim::initProbMax ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 148 of file EvtDsToKSKpPipPim.cc.

{ setProbMax( 824000 ); }

---

The documentation for this class was generated from the following files:

• EvtDsToKSKpPipPim.hh
• EvtDsToKSKpPipPim.cc