EvtVector4R Class Reference

#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 44 of file EvtVector4R.hh.

{}

Referenced by applyBoostTo(), boostTo, boostTo, cont(), cross(), dot(), dotr3(), mag2r3(), magr3(), operator*, operator*, operator*, operator*=(), operator+, operator+=(), operator-, operator-=(), operator/, operator/=(), operator<<, operator=(), rotateEuler, and scalartripler3().

EvtVector4R() [2/2]

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

Definition at line 32 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 107 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 96 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(), boostTo, boostTo, and EvtPto3P::initDaughters().

applyRotateEuler()

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

Definition at line 69 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(), EvtGenKine::PhaseSpace(), EvtGenKine::PhaseSpacePole(), and rotateEuler.

cont()

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

inline

Definition at line 156 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 153 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(), and EvtmyEulerAngles::EvtmyEulerAngles().

d3mag()

double EvtVector4R::d3mag

(

)

const

Definition at line 168 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(), rhopifull::amps1(), VVS::amps1(), EvtHelSys::Angles(), EvtAngH2::decay(), EvtAngSam3::decay(), EvtBody3::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(), EvtAmpA1(), EvtAmpA2(), EvtDecayAngleChi(), EvtConExc::findMaxXS(), EvtConExc::gamHXSection(), and EvtConExc::Rad2difXs().

dot()

double EvtVector4R::dot

(

const EvtVector4R &

v2

)

const

Definition at line 181 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(), and EvtDecayAngleChi().

dotr3()

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

Definition at line 210 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 171 of file EvtVector4R.hh.

{ return v[i]; }

Referenced by EvtMHelAmp::amplitude(), AngularSam::amps(), EvtHelSys::Angles(), EvtDiracSpinor::applyBoostTo(), EvtRaritaSchwinger::applyBoostTo(), EvtTensor4C::applyBoostTo(), EvtVector4C::applyBoostTo(), applyBoostTo(), 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(), EvtDToKKpipi0::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(), 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(), EvtTensor4C::directProd, EvtTensor4C::directProd, EvtPHOTOS::doRadCorr(), EvtPhotonParticle::epsParentPhoton(), EvtAmpA1(), EvtAmpA2(), EvtEulerAngles::EvtEulerAngles(), EvtmyEulerAngles::EvtmyEulerAngles(), EvtVector4C::EvtVector4C(), EvtConExc::findMaxXS(), EvtGen::generateEvent(), EvtConExc::hadron_angle_sampling(), EvtRexc::hadron_angle_sampling(), EvtNeutrinoParticle::init(), EvtPhotonParticle::init(), EvtTensorParticle::init(), EvtPto3P::initDaughters(), EvtVector4C::operator*, EvtVector4C::operator*, EvtVector4C::operator*, operator<<(), and EvtSVVHelAmp::SVVHel().

mag2r3()

double EvtVector4R::mag2r3

(

const EvtVector4R &

p1

)

const

Definition at line 216 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 221 of file EvtVector4R.cc.

{ return sqrt( mag2r3( p1 ) ); }

mass()

double EvtVector4R::mass

(

)

const

Definition at line 40 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(), 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(), EvtConExc::findMaxXS(), EvtConExc::gamHXSection(), EvtGen::generateEvent(), EvtParticle::generateMassTree(), EvtParticle::mass(), EvtConExc::Rad1difXs(), EvtConExc::Rad2difXs(), EvtResonance::relBrWig(), EvtResonance2::resAmpl(), and scalartripler3().

mass2()

double EvtVector4R::mass2 ( ) const

inline

Definition at line 110 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(), EvtOmegaDalitz::decay(), EvtPhiDalitz::decay(), EvtTauHadnu::decay(), dotr3(), EvtDecayAngle(), EvtDecayPlaneNormalAngle(), EvtConExc::gamHXSection(), mag2r3(), and EvtResonance::relBrWig().

operator*=()

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

inline

Definition at line 130 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 90 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 100 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 140 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 80 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 231 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 199 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 175 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 173 of file EvtVector4R.hh.

{ v[i] = d; }

Referenced by EvtbTosllScalarAmp::CalcAmp(), EvtbTosllVectorAmp::CalcAmp(), EvtSemiLeptonicBaryonAmp::CalcAmp(), EvtSemiLeptonicTensorAmp::CalcAmp(), EvtSemiLeptonicVectorAmp::CalcAmp(), EvtbTosllAmp::CalcMaxProb(), EvtSemiLeptonicAmp::CalcMaxProb(), EvtBHadronic::decay(), EvtBTo3piCP::decay(), EvtBToKpipiCP::decay(), EvtCBTo3piMPP::decay(), EvtCBTo3piP00::decay(), EvtDDalitz::decay(), EvtHelPPJ::decay(), EvtJetSet::decay(), EvtJscont::decay(), EvtKstarnunu::decay(), EvtLunda::decay(), EvtLundCharm::decay(), EvtPhokhara::decay(), EvtPhokhara_4pi::decay(), EvtPhokhara_etagam::decay(), EvtPhokhara_etapgam::decay(), EvtPhokhara_K0K0::decay(), EvtPhokhara_KK::decay(), EvtPhokhara_Lambda::decay(), EvtPhokhara_LLB::decay(), EvtPhokhara_nnbar::decay(), EvtPhokhara_pi0gam::decay(), EvtPhokhara_pi0pi0pipi::decay(), EvtPhokhara_pipi::decay(), EvtPhokhara_pipieta::decay(), EvtPhokhara_pipipi0::decay(), EvtPhokhara_ppbar::decay(), EvtPycont::decay(), EvtPyGaGa::decay(), EvtPythia::decay(), EvtSingleParticle2::decay(), EvtSingleParticle::decay(), EvtSinglePoint::decay(), EvtSll::decay(), EvtSLN::decay(), EvtTauola::decay(), EvtTVSPwave::decay(), EvtVSPPwave::decay(), EvtVub::decay(), EvtVubHybrid::decay(), EvtVubNLO::decay(), EvtPHOTOS::doRadCorr(), EvtGen::generateEvent(), EvtParticle::get4Pos(), EvtConExc::hadron_angle_sampling(), EvtRexc::hadron_angle_sampling(), EvtHelSys::Helrotate(), EvtLambdac2pKpi::init(), EvtTrackGen::init(), EvtPto3P::initDaughters(), EvtGenKine::PhaseSpace(), EvtGenKine::PhaseSpacePole(), EvtConExc::SetP4(), and EvtConExc::SetP4Rvalue().

theta()

double EvtVector4R::theta

(

)

Definition at line 224 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 62 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 55 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 151 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 120 of file EvtVector4R.hh.

                                                                {
 
  return EvtVector4R( v2 ) *= c;
}

operator* [3/3]

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

friend

Definition at line 115 of file EvtVector4R.hh.

                                                                {
 
  return EvtVector4R( v2 ) *= c;
}

operator+

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

friend

Definition at line 166 of file EvtVector4R.hh.

                                                                             {
 
  return EvtVector4R( v1 ) += v2;
}

operator-

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

friend

Definition at line 161 of file EvtVector4R.hh.

                                                                             {
 
  return EvtVector4R( v1 ) -= v2;
}

operator/

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

friend

Definition at line 125 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 48 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:

• EvtVector4R.hh
• EvtVector4R.cc