#include <EvtVector4R.hh>

Public Member Functions
	EvtVector4R ()

	EvtVector4R (double e, double px, double py, double pz)

void	set (int i, double d)

void	set (double e, double px, double py, double pz)

EvtVector4R &	operator*= (double c)

EvtVector4R &	operator/= (double c)

EvtVector4R &	operator= (const EvtVector4R &v2)

EvtVector4R &	operator+= (const EvtVector4R &v2)

EvtVector4R &	operator-= (const EvtVector4R &v2)

double	get (int i) const

double	cont (const EvtVector4R &v4) const

double	mass2 () const

double	mass () const

void	applyRotateEuler (double alpha, double beta, double gamma)

void	applyBoostTo (const EvtVector4R &p4)

void	applyBoostTo (const EvtVector3R &boost)

EvtVector4R	cross (const EvtVector4R &v2)

double	dot (const EvtVector4R &v2) const

double	d3mag () const

double	scalartripler3 (const EvtVector4R &p1, const EvtVector4R &p2, const EvtVector4R &p3) const

double	dotr3 (const EvtVector4R &p1, const EvtVector4R &p2) const

double	mag2r3 (const EvtVector4R &p1) const

double	magr3 (const EvtVector4R &p1) const

double	theta ()

double	phi ()

Friends
EvtVector4R	rotateEuler (const EvtVector4R &rs, double alpha, double beta, double gamma)

EvtVector4R	boostTo (const EvtVector4R &rs, const EvtVector4R &p4)

EvtVector4R	boostTo (const EvtVector4R &rs, const EvtVector3R &boost)

EvtVector4R	operator* (double d, const EvtVector4R &v2)

EvtVector4R	operator* (const EvtVector4R &v2, double d)

EvtVector4R	operator/ (const EvtVector4R &v2, double d)

double	operator* (const EvtVector4R &v1, const EvtVector4R &v2)

EvtVector4R	operator+ (const EvtVector4R &v1, const EvtVector4R &v2)

EvtVector4R	operator- (const EvtVector4R &v1, const EvtVector4R &v2)

std::ostream &	operator<< (std::ostream &s, const EvtVector4R &v)

Detailed Description

Definition at line 29 of file EvtVector4R.hh.

Constructor & Destructor Documentation

◆ EvtVector4R() [1/2]

EvtVector4R::EvtVector4R ( )

inline

Definition at line 47 of file EvtVector4R.hh.

47{}

◆ EvtVector4R() [2/2]

EvtVector4R::EvtVector4R	(	double	e,
		double	px,
		double	py,
		double	pz )

Definition at line 34 of file EvtVector4R.cc.

                                                               {
  
  v[0]=e; v[1]=p1; v[2]=p2; v[3]=p3;
}

Member Function Documentation

◆ applyBoostTo() [1/2]

void EvtVector4R::applyBoostTo ( const EvtVector3R & boost )

Definition at line 120 of file EvtVector4R.cc.

                                                      {
 
  double bx,by,bz,gamma,b2;
 
  bx=boost.get(0);
  by=boost.get(1);
  bz=boost.get(2);
 
  double bxx=bx*bx;
  double byy=by*by;
  double bzz=bz*bz;
 
  b2=bxx+byy+bzz;
 
 
  if (b2==0.0){
    return;
  }
 
  assert(b2<1.0);
 
  gamma=1.0/sqrt(1-b2);
 
 
  double gb2=(gamma-1.0)/b2;
 
  double gb2xy=gb2*bx*by;
  double gb2xz=gb2*bx*bz;
  double gb2yz=gb2*by*bz;
 
  double gbx=gamma*bx;
  double gby=gamma*by;
  double gbz=gamma*bz;
 
  double e2=v[0];
  double px2=v[1];
  double py2=v[2];
  double pz2=v[3];
 
  v[0]=gamma*e2+gbx*px2+gby*py2+gbz*pz2;
 
  v[1]=gbx*e2+gb2*bxx*px2+px2+gb2xy*py2+gb2xz*pz2;
 
  v[2]=gby*e2+gb2*byy*py2+py2+gb2xy*px2+gb2yz*pz2;
 
  v[3]=gbz*e2+gb2*bzz*pz2+pz2+gb2yz*py2+gb2xz*px2;
 
  return;
 
}

◆ applyBoostTo() [2/2]

void EvtVector4R::applyBoostTo ( const EvtVector4R & p4 )

Definition at line 108 of file EvtVector4R.cc.

                                                   {
 
  double e=p4.get(0);
 
  EvtVector3R boost(p4.get(1)/e,p4.get(2)/e,p4.get(3)/e);
 
  applyBoostTo(boost);
 
  return;
 
}

Referenced by applyBoostTo(), and EvtPto3P::initDaughters().

◆ applyRotateEuler()

void EvtVector4R::applyRotateEuler	(	double	alpha,
		double	beta,
		double	gamma )

Definition at line 81 of file EvtVector4R.cc.

                                                                    {
 
  double sp=sin(phi);
  double st=sin(theta);
  double sk=sin(ksi);
  double cp=cos(phi);
  double ct=cos(theta);
  double ck=cos(ksi);
 
  double x=( ck*ct*cp-sk*sp)*v[1]+( -sk*ct*cp-ck*sp)*v[2]+st*cp*v[3];
  double y=( ck*ct*sp+sk*cp)*v[1]+(-sk*ct*sp+ck*cp)*v[2]+st*sp*v[3];
  double z=-ck*st*v[1]+sk*st*v[2]+ct*v[3];
 
  v[1]=x;
  v[2]=y;
  v[3]=z;
  
}

Referenced by EvtPto3P::initDaughters().

◆ cont()

double EvtVector4R::cont ( const EvtVector4R & v4 ) const

inline

Definition at line 163 of file EvtVector4R.hh.

                                                           {
  
  return v[0]*v4.v[0]-v[1]*v4.v[1]-
         v[2]*v4.v[2]-v[3]*v4.v[3];
}

◆ cross()

EvtVector4R EvtVector4R::cross ( const EvtVector4R & v2 )

Definition at line 171 of file EvtVector4R.cc.

                                                     {
 
  //Calcs the cross product.  Added by djl on July 27, 1995.
  //Modified for real vectros by ryd Aug 28-96
 
  EvtVector4R temp;
  
  temp.v[0] = 0.0; 
  temp.v[1] = v[2]*p2.v[3] - v[3]*p2.v[2];
  temp.v[2] = v[3]*p2.v[1] - v[1]*p2.v[3];
  temp.v[3] = v[1]*p2.v[2] - v[2]*p2.v[1];
 
  return temp;
}

Referenced by EvtAngSam3::decay(), EvtHAngSam3::decay(), EvtJTO3P::decay(), EvtDecayAngleChi(), EvtmyEulerAngles::EvtmyEulerAngles(), EvtmyEulerAngles::EvtmyEulerAngles(), EvtmyEulerAngles::getphi(), and EvtmyEulerAngles::getpsi().

◆ d3mag()

double EvtVector4R::d3mag ( ) const

Definition at line 186 of file EvtVector4R.cc.

{
  double temp;
 
  temp = v[1]*v[1]+v[2]*v[2]+v[3]*v[3];
 
  temp = sqrt( temp );
  
  return temp;
} // r3mag

Referenced by EvtMHelAmp::amplitude(), AngularSam::amps(), rhopifull::amps1(), VVS::amps1(), EvtHelSys::Angles(), EvtAngH2::decay(), EvtAngSam3::decay(), EvtBody3::decay(), EvtCalHelAmp::decay(), EvtConExc::decay(), EvtDeBD::decay(), EvtHAngSam3::decay(), EvtHypWK::decay(), EvtJ2BB1::decay(), EvtJPE::decay(), EvtMBody3::decay(), EvtSingleParticle2::decay(), EvtSingleParticle::decay(), EvtSSDCP::decay(), EvtSTS::decay(), EvtSTSCP::decay(), EvtSVPCP::decay(), EvtSVPHelAmp::decay(), EvtSVS::decay(), EvtSVSCP::decay(), EvtSVSNONCPEIGEN::decay(), EvtTSS::decay(), EvtTVSPwave::decay(), EvtVSPPwave::decay(), EvtVSS::decay(), EvtVSSMix::decay(), EvtVVSPwave::decay(), EvtConExc::difgamXs(), EvtAmpA1(), EvtAmpA2(), EvtDecayAngleChi(), EvtConExc::findMaxXS(), EvtConExc::gamHXSection(), EvtHelSys::getHelAng(), EvtmyEulerAngles::getphi(), EvtmyEulerAngles::getpsi(), EvtmyEulerAngles::gettheta(), and EvtConExc::Rad2difXs().

◆ dot()

double EvtVector4R::dot ( const EvtVector4R & v2 ) const

Definition at line 199 of file EvtVector4R.cc.

                                                    {
 
  //Returns the dot product of the 3 momentum.  Added by
  //djl on July 27, 1995.  for real!!!
 
  double temp;
 
  temp = v[1]*p2.v[1];
  temp += v[2]*p2.v[2];
  temp += v[3]*p2.v[3];
 
  return temp;
 
} //dot

Referenced by EvtHAngSam3::decay(), EvtDecayAngleChi(), EvtmyEulerAngles::getphi(), EvtmyEulerAngles::getpsi(), and EvtmyEulerAngles::gettheta().

◆ dotr3()

double EvtVector4R::dotr3	(	const EvtVector4R &	p1,
		const EvtVector4R &	p2 ) const

Definition at line 229 of file EvtVector4R.cc.

{
    return 1/mass2() * ((*this) * p1) * ((*this) * p2) - p1 * p2;
}

Referenced by EvtDecayAnglePhi().

◆ get()

double EvtVector4R::get ( int i ) const

inline

Definition at line 179 of file EvtVector4R.hh.

                                          {
  return v[i];
}

Referenced by EvtMHelAmp::amplitude(), AngularSam::amps(), EvtHelSys::Angles(), EvtTensor4C::cont1(), EvtTensor4C::cont2(), EvtAngH2::decay(), EvtAngSam3::decay(), EvtBody3::decay(), EvtBsquark::decay(), EvtBtoXsll::decay(), EvtConExc::decay(), EvtD0ToKKpipi::decay(), EvtD0ToKpiEtap::decay(), EvtD0ToKpipi0::decay(), EvtD0ToKpipi0pi0::decay(), EvtD0ToKpipipi::decay(), EvtD0ToKSKppimpi0::decay(), EvtD0ToKSLKK::decay(), EvtD0ToKSpi0eta::decay(), EvtD0ToKSpi0pi0::decay(), EvtD0ToKSpi0pi0pi0::decay(), EvtD0ToKSpipipi0pi0::decay(), EvtD0TopipiEta::decay(), EvtDsToEta3pi::decay(), EvtDsToEtapi2pi0::decay(), EvtDsToEtapipi0::decay(), EvtDsToEtappipi0::decay(), EvtDsToKKpi::decay(), EvtDsToKKpipi0::decay(), EvtDsToKKpipipi::decay(), EvtDsToKpipi::decay(), EvtDsToKpPipPimPi0::decay(), EvtDsToKSKmPipPip::decay(), EvtDsToKSKpi0::decay(), EvtDsToKSKpPipPim::decay(), EvtDsToKSKSpi::decay(), EvtDsToKSpipi0::decay(), EvtDsTopipi0pi0::decay(), EvtDsTopipipipi0pi0::decay(), EvtDsToPipPipPimPi0::decay(), EvtDToKKpi::decay(), EvtDToKmPipPipPi0::decay(), EvtDToKppipi::decay(), EvtDToKppipipi0::decay(), EvtDToKSKmpippip::decay(), EvtDToKSKSK::decay(), EvtDToKSKSpi::decay(), EvtDToKSpieta::decay(), EvtDToKSpipi0pi0::decay(), EvtDToKSpipipi::decay(), EvtDTopiEtaEta::decay(), EvtDTopipi0Eta::decay(), EvtDToPiPi0Etap::decay(), EvtDTopipi0pi0::decay(), EvtEtaDalitz::decay(), EvtFDC::decay(), EvtFlatQ2::decay(), EvtHAngSam3::decay(), EvtHelPPJ::decay(), EvtJ2BB1::decay(), EvtJetSet::decay(), EvtJTO3P::decay(), EvtKstarstargamma::decay(), EvtLambdaP_BarGamma::decay(), EvtLNuGamma::decay(), EvtLunda::decay(), EvtLundCharm::decay(), EvtMBody3::decay(), EvtOmegaDalitz::decay(), EvtOpenCharm::decay(), EvtSingleParticle2::decay(), EvtSingleParticle::decay(), EvtTauola::decay(), EvtTVSPwave::decay(), EvtVectorIsr::decay(), EvtVVSPwave::decay(), EvtConExc::difgamXs(), EvtPHOTOS::doRadCorr(), EvtPhotonParticle::epsParentPhoton(), EvtAmpA1(), EvtAmpA2(), EvtEulerAngles::EvtEulerAngles(), EvtmyEulerAngles::EvtmyEulerAngles(), EvtVector4C::EvtVector4C(), EvtConExc::findMaxXS(), EvtConExc::gamHXSection(), EvtGen::generateEvent(), EvtConExc::hadron_angle_sampling(), EvtRexc::hadron_angle_sampling(), EvtPto3P::initDaughters(), operator<<(), EvtFDC::setEvtMomentum(), and EvtSVVHelAmp::SVVHel().

◆ mag2r3()

double EvtVector4R::mag2r3 ( const EvtVector4R & p1 ) const

Definition at line 236 of file EvtVector4R.cc.

{
    return Square((*this) * p1)/mass2() - p1.mass2();
}

Referenced by EvtDecayAnglePhi(), and magr3().

◆ magr3()

double EvtVector4R::magr3 ( const EvtVector4R & p1 ) const

Definition at line 242 of file EvtVector4R.cc.

{
    return sqrt(mag2r3(p1));
}

◆ mass()

double EvtVector4R::mass ( ) const

Definition at line 39 of file EvtVector4R.cc.

                              {
 
  double m2=v[0]*v[0]-v[1]*v[1]-v[2]*v[2]-v[3]*v[3];
 
  if (m2>0.0) {
    return sqrt(m2);
  }
  else{
    return 0.0;
  }
}

Referenced by EvtSemiLeptonicTensorAmp::CalcAmp(), EvtPsi3Sdecay::choseDecay(), EvtBtoXsll::decay(), EvtConExc::decay(), EvtKstarstargamma::decay(), EvtmH2::decay(), EvtmPhsp::decay(), EvtPhiDalitz::decay(), EvtPhsp::decay(), EvtPhspStepf::decay(), EvtSSDCP::decay(), EvtSTS::decay(), EvtSVPCP::decay(), EvtSVPHelAmp::decay(), EvtSVSCP::decay(), EvtSVSCPLH::decay(), EvtSVSNONCPEIGEN::decay(), EvtVPHOtoVISR::decay(), EvtVPHOtoVISRHi::decay(), EvtConExc::difgamXs(), EvtDecayAnglePhi(), EvtDecayPlaneNormalAngle(), EvtConExc::findMaxXS(), EvtConExc::gamHXSection(), EvtParticle::generateMassTree(), EvtParticle::mass(), EvtConExc::Rad1difXs(), EvtConExc::Rad2difXs(), EvtResonance::relBrWig(), EvtResonance2::resAmpl(), scalartripler3(), and EvtMBreitWigner::shape().

◆ mass2()

double EvtVector4R::mass2 ( ) const

inline

Definition at line 116 of file EvtVector4R.hh.

                                      {
 
  return v[0]*v[0]-v[1]*v[1]-v[2]*v[2]-v[3]*v[3];
}

Referenced by rhopifull::amps1(), VVS::amps1(), EvtSemiLeptonicTensorAmp::CalcAmp(), EvtBtoXsll::decay(), EvtKstarnunu::decay(), EvtOmegaDalitz::decay(), EvtPhiDalitz::decay(), EvtPi0Dalitz::decay(), EvtTauHadnu::decay(), dotr3(), EvtDecayAngle(), EvtDecayPlaneNormalAngle(), EvtConExc::gamHXSection(), mag2r3(), and EvtResonance::relBrWig().

◆ operator*=()

EvtVector4R & EvtVector4R::operator*= ( double c )

inline

Definition at line 136 of file EvtVector4R.hh.

                                                   {
 
  v[0]*=c;  
  v[1]*=c;  
  v[2]*=c;  
  v[3]*=c;  
 
  return *this;
}

◆ operator+=()

EvtVector4R & EvtVector4R::operator+= ( const EvtVector4R & v2 )

inline

Definition at line 96 of file EvtVector4R.hh.

                                                                {
 
  v[0]+=v2.v[0];
  v[1]+=v2.v[1];
  v[2]+=v2.v[2];
  v[3]+=v2.v[3];
  
  return *this; 
}

◆ operator-=()

EvtVector4R & EvtVector4R::operator-= ( const EvtVector4R & v2 )

inline

Definition at line 106 of file EvtVector4R.hh.

                                                                {
 
  v[0]-=v2.v[0];
  v[1]-=v2.v[1];
  v[2]-=v2.v[2];
  v[3]-=v2.v[3];
  
  return *this; 
}

◆ operator/=()

EvtVector4R & EvtVector4R::operator/= ( double c )

inline

Definition at line 146 of file EvtVector4R.hh.

                                                   {
 
  double cinv=1.0/c;  
  v[0]*=cinv;  
  v[1]*=cinv;  
  v[2]*=cinv;  
  v[3]*=cinv;  
 
  return *this;
}

◆ operator=()

EvtVector4R & EvtVector4R::operator= ( const EvtVector4R & v2 )

inline

Definition at line 86 of file EvtVector4R.hh.

                                                               {
 
  v[0]=v2.v[0];
  v[1]=v2.v[1];
  v[2]=v2.v[2];
  v[3]=v2.v[3];
  
  return *this; 
}

◆ phi()

double EvtVector4R::phi ( )

Definition at line 254 of file EvtVector4R.cc.

                       {
  return (v[1]==0 && v[2]==0)?0: atan2(v[2],v[1]);
}

Referenced by applyRotateEuler().

◆ scalartripler3()

double EvtVector4R::scalartripler3	(	const EvtVector4R &	p1,
		const EvtVector4R &	p2,
		const EvtVector4R &	p3 ) const

Definition at line 217 of file EvtVector4R.cc.

{
    EvtVector4C lc=dual(directProd(*this, p1)).cont2(p2);
    EvtVector4R  l(real(lc.get(0)), real(lc.get(1)), real(lc.get(2)),
            real(lc.get(3)));
 
    return -1.0/mass() * (l * p3);
}

Referenced by EvtDecayAnglePhi().

◆ set() [1/2]

void EvtVector4R::set	(	double	e,
		double	px,
		double	py,
		double	pz )

inline

Definition at line 188 of file EvtVector4R.hh.

                                                                   {
 
  v[0]=e;
  v[1]=p1;
  v[2]=p2;
  v[3]=p3;
}

◆ set() [2/2]

void EvtVector4R::set	(	int	i,
		double	d )

inline

Definition at line 183 of file EvtVector4R.hh.

Referenced by EvtbTosllScalarAmp::CalcAmp(), EvtbTosllVectorAmp::CalcAmp(), EvtSemiLeptonicBaryonAmp::CalcAmp(), EvtSemiLeptonicTensorAmp::CalcAmp(), EvtSemiLeptonicVectorAmp::CalcAmp(), EvtbTosllAmp::CalcMaxProb(), EvtSemiLeptonicAmp::CalcMaxProb(), EvtBHadronic::decay(), EvtDDalitz::decay(), EvtHelPPJ::decay(), EvtKstarnunu::decay(), EvtPhokhara_Lambda::decay(), EvtSll::decay(), EvtSLN::decay(), EvtTVSPwave::decay(), EvtVSPPwave::decay(), EvtPHOTOS::doRadCorr(), EvtGen::generateEvent(), EvtParticle::get4Pos(), EvtConExc::hadron_angle_sampling(), EvtRexc::hadron_angle_sampling(), EvtHelSys::Helrotate(), EvtPto3P::initDaughters(), EvtParticle::setp(), EvtConExc::SetP4(), and EvtConExc::SetP4Rvalue().

◆ theta()

double EvtVector4R::theta ( )

Definition at line 249 of file EvtVector4R.cc.

                         {
  return (v[1]==0 && v[2]==0 && v[3]==0 )?0: atan2(sqrt(v[1]*v[1]+v[2]*v[2]),v[3]);
}

Referenced by EvtConExc::angularSampling(), EvtRexc::angularSampling(), applyRotateEuler(), EvtAngSamLab::decay(), and EvtConExc::photonSampling().

Friends And Related Symbol Documentation

◆ boostTo [1/2]

EvtVector4R boostTo	(	const EvtVector4R &	rs,
		const EvtVector3R &	boost )

friend

Definition at line 70 of file EvtVector4R.cc.

                                     {
 
  EvtVector4R tmp(rs);
  tmp.applyBoostTo(boost);
  return tmp;
 
}

◆ boostTo [2/2]

EvtVector4R boostTo	(	const EvtVector4R &	rs,
		const EvtVector4R &	p4 )

friend

Definition at line 61 of file EvtVector4R.cc.

                                  {
 
  EvtVector4R tmp(rs);
  tmp.applyBoostTo(p4);
  return tmp;
 
}

◆ operator* [1/3]

double operator*	(	const EvtVector4R &	v1,
		const EvtVector4R &	v2 )

friend

Definition at line 157 of file EvtVector4R.hh.

                                                                    {
 
  return v1.v[0]*v2.v[0]-v1.v[1]*v2.v[1]-
         v1.v[2]*v2.v[2]-v1.v[3]*v2.v[3];
}

◆ operator* [2/3]

EvtVector4R operator*	(	const EvtVector4R &	v2,
		double	d )

friend

Definition at line 126 of file EvtVector4R.hh.

                                                            {
  
  return EvtVector4R(v2)*=c;
}

◆ operator* [3/3]

EvtVector4R operator*	(	double	d,
		const EvtVector4R &	v2 )

friend

Definition at line 121 of file EvtVector4R.hh.

                                                            {
  
  return EvtVector4R(v2)*=c;
}

◆ operator+

EvtVector4R operator+	(	const EvtVector4R &	v1,
		const EvtVector4R &	v2 )

friend

Definition at line 174 of file EvtVector4R.hh.

                                                                         {
  
  return EvtVector4R(v1)+=v2;
}

◆ operator-

EvtVector4R operator-	(	const EvtVector4R &	v1,
		const EvtVector4R &	v2 )

friend

Definition at line 169 of file EvtVector4R.hh.

                                                                         {
  
  return EvtVector4R(v1)-=v2;
}

◆ operator/

EvtVector4R operator/	(	const EvtVector4R &	v2,
		double	d )

friend

Definition at line 131 of file EvtVector4R.hh.

                                                            {
  
  return EvtVector4R(v2)/=c;
}

◆ operator<<

std::ostream & operator<<	(	std::ostream &	s,
		const EvtVector4R &	v )

friend

◆ rotateEuler

EvtVector4R rotateEuler	(	const EvtVector4R &	rs,
		double	alpha,
		double	beta,
		double	gamma )

friend

Definition at line 52 of file EvtVector4R.cc.

                                                  {
 
  EvtVector4R tmp(rs);
  tmp.applyRotateEuler(alpha,beta,gamma);
  return tmp;
 
}

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

Public Member Functions

Friends

Detailed Description

Constructor & Destructor Documentation

◆ EvtVector4R() [1/2]

◆ EvtVector4R() [2/2]

Member Function Documentation

◆ applyBoostTo() [1/2]

◆ applyBoostTo() [2/2]

◆ applyRotateEuler()

◆ cont()

◆ cross()

◆ d3mag()

◆ dot()

◆ dotr3()

◆ get()

◆ mag2r3()

◆ magr3()

◆ mass()

◆ mass2()

◆ operator*=()

◆ operator+=()

◆ operator-=()

◆ operator/=()

◆ operator=()

◆ phi()

◆ scalartripler3()

◆ set() [1/2]

◆ set() [2/2]

◆ theta()

Friends And Related Symbol Documentation

◆ boostTo [1/2]

◆ boostTo [2/2]

◆ operator* [1/3]

◆ operator* [2/3]

◆ operator* [3/3]

◆ operator+

◆ operator-

◆ operator/

◆ operator<<

◆ rotateEuler