#include "EvtGenBase/EvtPatches.hh"
#include <math.h>
#include "EvtGenBase/EvtKine.hh"
#include "EvtGenBase/EvtConst.hh"
#include "EvtGenBase/EvtVector4R.hh"
#include "EvtGenBase/EvtVector4C.hh"
#include "EvtGenBase/EvtTensor4C.hh"
#include "EvtGenBase/EvtdFunction.hh"
#include "EvtGenBase/EvtReport.hh"

Functions
double	EvtDecayAngle (const EvtVector4R &p, const EvtVector4R &q, const EvtVector4R &d)

double	EvtDecayAngleChi (const EvtVector4R &p4_p, const EvtVector4R &p4_d1, const EvtVector4R &p4_d2, const EvtVector4R &p4_h1, const EvtVector4R &p4_h2)

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

double	EvtDecayAnglePhi (const EvtVector4R &z, 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

◆ EvtDecayAngle()

double EvtDecayAngle	(	const EvtVector4R &	p,
		const EvtVector4R &	q,
		const EvtVector4R &	d )

Definition at line 33 of file EvtKine.cc.

                                   {
 
  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;
 
}

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

◆ EvtDecayAngleChi()

double EvtDecayAngleChi	(	const EvtVector4R &	p4_p,
		const EvtVector4R &	p4_d1,
		const EvtVector4R &	p4_d2,
		const EvtVector4R &	p4_h1,
		const EvtVector4R &	p4_h2 )

Definition at line 49 of file EvtKine.cc.

                                       {
 
  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;
  
}

◆ EvtDecayAnglePhi()

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

Definition at line 111 of file EvtKine.cc.

{
    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;
}

Referenced by EvtMHelAmp::amplitude().

◆ EvtDecayPlaneNormalAngle()

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

Definition at line 94 of file EvtKine.cc.

                                                                      {
 
  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());
 
}

◆ wignerD()

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

Definition at line 132 of file EvtKine.cc.

{
    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
}