/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Generator/BesEvtGen/BesEvtGen-00-03-58/src/EvtGen/EvtGenBase/EvtKine.hh File Reference

Go to the source code of this file.

Functions
double	EvtDecayAngle (const EvtVector4R &, const EvtVector4R &, const EvtVector4R &)
double	EvtDecayAngleChi (const EvtVector4R &, const EvtVector4R &, const EvtVector4R &, const EvtVector4R &, const EvtVector4R &)
double	EvtDecayPlaneNormalAngle (const EvtVector4R &p, const EvtVector4R &q, const EvtVector4R &d1, const EvtVector4R &d2)
double	EvtDecayAnglePhi (const EvtVector4R &g, const EvtVector4R &p, const EvtVector4R &q, const EvtVector4R &d)
EvtComplex	wignerD (int j, int m1, int m2, double phi, double theta, double gamma)

---

Function Documentation

double EvtDecayAngle	(	const EvtVector4R &	,
		const EvtVector4R &	,
		const EvtVector4R &
	)

Definition at line 33 of file EvtKine.cc.

References EvtVector4R::mass2(), and q.

Referenced by EvtMHelAmp::amplitude(), and EvtResonance::resAmpl().

                                           {

  double pd=p*d;
  double pq=p*q;
  double qd=q*d;
  double mp2=p.mass2();
  double mq2=q.mass2();
  double md2=d.mass2();

  double cost=(pd*mq2-pq*qd)/sqrt((pq*pq-mq2*mp2)*(qd*qd-mq2*md2));

  return cost;

}

double EvtDecayAngleChi	(	const EvtVector4R &	,
		const EvtVector4R &	,
		const EvtVector4R &	,
		const EvtVector4R &	,
		const EvtVector4R &
	)

Definition at line 49 of file EvtKine.cc.

References boostTo(), EvtVector4R::cross(), EvtVector4R::d3mag(), EvtVector4R::dot(), EvtConst::twoPi, and x.

                                                   {

  EvtVector4R p4_d1p,p4_h1p,p4_h2p,p4_d2p;
  

  // boost all vectors parent restframe

  p4_d1p=boostTo(p4_d1,p4_p);
  p4_d2p=boostTo(p4_d2,p4_p);
  p4_h1p=boostTo(p4_h1,p4_p);
  p4_h2p=boostTo(p4_h2,p4_p);
  

  EvtVector4R d1_perp,d1_prime,h1_perp;
  EvtVector4R D;
  
  D=p4_d1p+p4_d2p;

  d1_perp=p4_d1p-(D.dot(p4_d1p)/D.dot(D))*D;
  h1_perp=p4_h1p-(D.dot(p4_h1p)/D.dot(D))*D;
  
  // orthogonal to both D and d1_perp
  
  d1_prime=D.cross(d1_perp);
  
  d1_perp= d1_perp/d1_perp.d3mag();
  d1_prime= d1_prime/d1_prime.d3mag();
  
  double x,y;
  
  x=d1_perp.dot(h1_perp);
  y=d1_prime.dot(h1_perp);
  
  double chi=atan2(y,x);
  
  if (chi<0.0) chi+=EvtConst::twoPi;
  
  return chi;
  
}

double EvtDecayAnglePhi	(	const EvtVector4R &	g,
		const EvtVector4R &	p,
		const EvtVector4R &	q,
		const EvtVector4R &	d
	)

Definition at line 111 of file EvtKine.cc.

References alpha, EvtVector4R::dotr3(), EvtVector4R::mag2r3(), EvtVector4R::mass(), q, EvtVector4R::scalartripler3(), EvtConst::twoPi, and x.

Referenced by EvtMHelAmp::amplitude().

{
    double eq = (p * q) / p.mass();
    double ed = (p * d) / p.mass();
    double mq = q.mass();
    double q2 = p.mag2r3(q);
    double qd = p.dotr3(q,d);
    double zq = p.dotr3(z,q);
    double zd = p.dotr3(z,d);
    double alpha = (eq - mq)/(q2 * mq) * qd - ed/mq;
                                                                                                              
    double y = p.scalartripler3(z,q,d) + alpha * p.scalartripler3(z,q,q);
    double x = (zq * (qd + alpha * q2) - q2 * (zd + alpha * zq)) / sqrt(q2);
                                                                                                              
    double phi = atan2(y,x);
                                                                                                              
    return phi<0 ? (phi+EvtConst::twoPi) : phi;
}

double EvtDecayPlaneNormalAngle	(	const EvtVector4R &	p,
		const EvtVector4R &	q,
		const EvtVector4R &	d1,
		const EvtVector4R &	d2
	)

Definition at line 94 of file EvtKine.cc.

References EvtTensor4C::cont2(), directProd(), dual(), EvtVector4C::get(), EvtVector4R::mass2(), and q.

                                                                      {

  EvtVector4C lc=dual(directProd(d1,d2)).cont2(q);

  EvtVector4R l(real(lc.get(0)),real(lc.get(1)),
                real(lc.get(2)),real(lc.get(3)));

  double pq=p*q;

  return q.mass()*(p*l)/sqrt(-(pq*pq-p.mass2()*q.mass2())*l.mass2());


}

EvtComplex wignerD	(	int	j,
		int	m1,
		int	m2,
		double	phi,
		double	theta,
		double	gamma
	)

Definition at line 132 of file EvtKine.cc.

References EvtdFunction::d(), and exp().

Referenced by EvtMHelAmp::amplitude().

{
    EvtComplex gp(0.0, -phi*m1);  
    EvtComplex gm(0.0, -gamma*m2);
    //double d = EvtdFunction::d(j, m1, m2, theta);
                                                                                                              
    return exp( gp ) * EvtdFunction::d(2*j, 2*m1,2*m2, theta) * exp( gm ); //pingrg, 2010-11-29, d(j,m1,m2,theta)->d(2*j,2*m1,2*m2,theta)consistent with Djmn in EvtHelSys
}

---