#include <EvtParticle.hh>

Inheritance diagram for EvtParticle:

Public Member Functions
	EvtParticle ()

virtual	~EvtParticle ()

virtual EvtVector4C	epsParent (int i) const

virtual EvtVector4C	eps (int i) const

virtual EvtVector4C	epsParentPhoton (int i)

virtual EvtVector4C	epsPhoton (int i)

virtual EvtDiracSpinor	spParent (int) const

virtual EvtDiracSpinor	sp (int) const

virtual EvtDiracSpinor	spParentNeutrino () const

virtual EvtDiracSpinor	spNeutrino () const

virtual EvtTensor4C	epsTensorParent (int i) const

virtual EvtTensor4C	epsTensor (int i) const

virtual void	init (EvtId part_n, const EvtVector4R &p4)=0

void	addDaug (EvtParticle *node)

void	decay ()

void	deleteTree ()

void	deleteDaughters (bool keepChannel=false)

void	setChannel (int i)

void	setGeneratorFlag (int flag)

int	getGeneratorFlag ()

void	setIntFlag (std::vector< int > vi)

std::vector< int >	getIntFlag ()

void	makeDaughters (int ndaug, EvtId *id)

double	initializePhaseSpace (int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)

EvtParticle *	getDaug (int i)

EvtParticle *	nextIter (EvtParticle *rootOfTree=0)

void	makeStdHep (EvtStdHep &stdhep, EvtSecondary &secondary, EvtId *stable_parent_ihep)

void	makeStdHep (EvtStdHep &stdhep)

EvtVector4R	getP4Lab ()

EvtVector4R	getP4Restframe ()

EvtVector4R	get4Pos ()

EvtParticle *	getParent ()

void	insertDaugPtr (int idaug, EvtParticle *partptr)

double	mass () const

int	firstornot () const

void	setFirstOrNot ()

void	resetFirstOrNot ()

EvtId	getId () const

EvtSpinType::spintype	getSpinType () const

int	getSpinStates () const

const EvtVector4R &	getP4 () const

void	setP4 (const EvtVector4R &p4)

int	getChannel () const

int	getNDaug () const

void	resetNDaug ()

void	printTree () const

void	printTreeRec (int level) const

std::string	writeTreeRec (std::string) const

void	dumpTree () const

void	dumpTreeRec (int level, int dj) const

std::string	treeStr () const

std::string	treeStrRec (int level) const

void	printParticle () const

void	setLifetime (double tau)

void	setLifetime ()

double	getLifetime ()

void	setDiagonalSpinDensity ()

void	setVectorSpinDensity ()

void	setPolarizedSpinDensity (double r00, double r11, double r22)

void	setSpinDensityForward (const EvtSpinDensity &rho)

void	setSpinDensityForwardHelicityBasis (const EvtSpinDensity &rho)

void	setSpinDensityForwardHelicityBasis (const EvtSpinDensity &rho, double alpha, double beta, double gamma)

virtual EvtSpinDensity	rotateToHelicityBasis () const =0

virtual EvtSpinDensity	rotateToHelicityBasis (double alpha, double beta, double gamma) const =0

EvtSpinDensity	getSpinDensityForward ()

void	setSpinDensityBackward (const EvtSpinDensity &rho)

EvtSpinDensity	getSpinDensityBackward ()

void	noLifeTime ()

void	setId (EvtId id)

void	initDecay (bool useMinMass=false)

void	generateMassTree ()

double	compMassProb ()

void	setMass (double m)

bool	isInitialized ()

bool	hasValidP4 ()

bool	isDecayed ()

double *	decayProb ()

void	setDecayProb (double p)

void	setInclusiveMode (int im)

int	getInclusiveMode ()

	EvtParticle ()

virtual	~EvtParticle ()

virtual EvtVector4C	epsParent (int i) const

virtual EvtVector4C	eps (int i) const

virtual EvtVector4C	epsParentPhoton (int i)

virtual EvtVector4C	epsPhoton (int i)

virtual EvtDiracSpinor	spParent (int) const

virtual EvtDiracSpinor	sp (int) const

virtual EvtDiracSpinor	spParentNeutrino () const

virtual EvtDiracSpinor	spNeutrino () const

virtual EvtTensor4C	epsTensorParent (int i) const

virtual EvtTensor4C	epsTensor (int i) const

virtual void	init (EvtId part_n, const EvtVector4R &p4)=0

void	addDaug (EvtParticle *node)

void	decay ()

void	deleteTree ()

void	deleteDaughters (bool keepChannel=false)

void	setChannel (int i)

void	setGeneratorFlag (int flag)

int	getGeneratorFlag ()

void	setIntFlag (std::vector< int > vi)

std::vector< int >	getIntFlag ()

void	makeDaughters (int ndaug, EvtId *id)

double	initializePhaseSpace (int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)

EvtParticle *	getDaug (int i)

EvtParticle *	nextIter (EvtParticle *rootOfTree=0)

void	makeStdHep (EvtStdHep &stdhep, EvtSecondary &secondary, EvtId *stable_parent_ihep)

void	makeStdHep (EvtStdHep &stdhep)

EvtVector4R	getP4Lab ()

EvtVector4R	getP4Restframe ()

EvtVector4R	get4Pos ()

EvtParticle *	getParent ()

void	insertDaugPtr (int idaug, EvtParticle *partptr)

double	mass () const

int	firstornot () const

void	setFirstOrNot ()

void	resetFirstOrNot ()

EvtId	getId () const

EvtSpinType::spintype	getSpinType () const

int	getSpinStates () const

const EvtVector4R &	getP4 () const

void	setP4 (const EvtVector4R &p4)

int	getChannel () const

int	getNDaug () const

void	resetNDaug ()

void	printTree () const

void	printTreeRec (int level) const

std::string	writeTreeRec (std::string) const

void	dumpTree () const

void	dumpTreeRec (int level, int dj) const

std::string	treeStr () const

std::string	treeStrRec (int level) const

void	printParticle () const

void	setLifetime (double tau)

void	setLifetime ()

double	getLifetime ()

void	setDiagonalSpinDensity ()

void	setVectorSpinDensity ()

void	setPolarizedSpinDensity (double r00, double r11, double r22)

void	setSpinDensityForward (const EvtSpinDensity &rho)

void	setSpinDensityForwardHelicityBasis (const EvtSpinDensity &rho)

void	setSpinDensityForwardHelicityBasis (const EvtSpinDensity &rho, double alpha, double beta, double gamma)

virtual EvtSpinDensity	rotateToHelicityBasis () const =0

virtual EvtSpinDensity	rotateToHelicityBasis (double alpha, double beta, double gamma) const =0

EvtSpinDensity	getSpinDensityForward ()

void	setSpinDensityBackward (const EvtSpinDensity &rho)

EvtSpinDensity	getSpinDensityBackward ()

void	noLifeTime ()

void	setId (EvtId id)

void	initDecay (bool useMinMass=false)

void	generateMassTree ()

double	compMassProb ()

void	setMass (double m)

bool	isInitialized ()

bool	hasValidP4 ()

bool	isDecayed ()

double *	decayProb ()

void	setDecayProb (double p)

void	setInclusiveMode (int im)

int	getInclusiveMode ()

Protected Member Functions
void	setp (double e, double px, double py, double pz)

void	setp (const EvtVector4R &p4)

void	setpart_num (EvtId particle_number)

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

void	setp (const EvtVector4R &p4)

void	setpart_num (EvtId particle_number)

Protected Attributes
bool	_validP4

Detailed Description

Definition at line 42 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

Constructor & Destructor Documentation

◆ EvtParticle() [1/2]

EvtParticle::EvtParticle ( )

Default constructor.

Definition at line 60 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                         {
   _ndaug=0;
   _parent=0;
   _channel=-10;
   _t=0.0;
   _genlifetime=1;
   _first=1;
   _isInit=false;
   _validP4=false;
   _isDecayed=false;
   _decayProb=0;
   //   _mix=false;
}

◆ ~EvtParticle() [1/2]

EvtParticle::~EvtParticle ( )

virtual

Destructor.

Definition at line 56 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                          {
  delete _decayProb;
}

◆ EvtParticle() [2/2]

EvtParticle::EvtParticle ( )

Default constructor.

◆ ~EvtParticle() [2/2]

virtual EvtParticle::~EvtParticle ( )

virtual

Destructor.

Member Function Documentation

◆ addDaug() [1/2]

void EvtParticle::addDaug ( EvtParticle * node )

Add another daughter to the particle

Definition at line 104 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                           {
  node->_daug[node->_ndaug++]=this;
  _ndaug=0;
  _parent=node; 
}

Referenced by EvtKstarstargamma::decay(), EvtPHOTOS::doRadCorr(), and makeDaughters().

◆ addDaug() [2/2]

void EvtParticle::addDaug ( EvtParticle * node )

Add another daughter to the particle

◆ compMassProb() [1/2]

double EvtParticle::compMassProb ( )

Definition at line 504 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                 {
 
  EvtParticle *p=this;
  //report(INFO,"EvtGen") << "compMassProb " << endl;
  //p->printTree();
  double mass=p->mass();
  double parMass=0.;
  if ( p->getParent()) { 
    parMass=p->getParent()->mass(); 
  }
 
  int nDaug=p->getNDaug();
  double *dMasses=0;
 
  int i;
  if ( nDaug>0 ) {
    dMasses=new double[nDaug];
    for (i=0; i<nDaug; i++) dMasses[i]=p->getDaug(i)->mass();
  }
 
  double temp=1.0;
  temp=EvtPDL::getMassProb(p->getId(), mass, parMass, nDaug, dMasses);
 
  //report(INFO,"EvtGen") << temp << " " << EvtPDL::name(p->getId()) << endl;
  //If the particle already has a mass, we dont need to include
  //it in the probability calculation
  if ( (!p->getParent() || _validP4 ) && temp>0.0 ) temp=1.; 
 
  delete [] dMasses;
  // if ( temp < 0.9999999 ) 
  for (i=0; i<nDaug; i++) {
    temp*=p->getDaug(i)->compMassProb();
  }
  return temp;
}

Referenced by compMassProb(), and generateMassTree().

◆ compMassProb() [2/2]

double EvtParticle::compMassProb ( )

◆ decay() [1/2]

void EvtParticle::decay ( )

Decay particle

Definition at line 404 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                       {
 
  //P is particle to decay, typically 'this' but sometime
  //modified by mixing 
  EvtParticle* p=this;
  //Did it mix?
  //if ( p->getMixed() ) {
    //should take C(p) - this should only
    //happen the first time we call decay for this
    //particle
  //p->takeCConj();
  // p->setUnMixed();
  //}
  EvtSpinDensity myRho;  //pingrg test code
  EvtDecayBase *decayer;
  decayer = EvtDecayTable::getDecayFunc(p);
  //  if ( decayer ) {
  //    report(INFO,"EvtGen") << "calling decay for " << EvtPDL::name(p->getId()) << " " << p->mass() << " " << p->getP4() << " " << p->getNDaug() << " " << p << endl;
  //    report(INFO,"EvtGen") << "NDaug= " << decayer->getNDaug() << endl;
  //    int ti;
  //    for ( ti=0; ti<decayer->getNDaug(); ti++) 
  //      report(INFO,"EvtGen") << "Daug " << ti << " " << EvtPDL::name(decayer->getDaug(ti)) << endl;
  //  }
  //if (p->_ndaug>0) {
  //      report(INFO,"EvtGen") <<"Is decaying particle with daughters!!!!!"<<endl;
  //     ::abort();
    //return;
    //call initdecay first - April 29,2002 - Lange
  //}
 
  //if there are already daughters, then this step is already done!
  // figure out the masses
  if ( _ndaug == 0 ) generateMassTree();
  static EvtId BS0=EvtPDL::getId("B_s0");
  static EvtId BSB=EvtPDL::getId("anti-B_s0");
  static EvtId BD0=EvtPDL::getId("B0");
  static EvtId BDB=EvtPDL::getId("anti-B0"); 
  if ( _ndaug==1 &&  (getId()==BS0||getId()==BSB||getId()==BD0||getId()==BDB) ) 
    {getDaug(0)->decay();
      std::cout<<"if ( _ndaug==1 &&  (getId()==BS0||getId()==BSB||getId()==BD0||getId()==BDB) )"<<endl;
    }
 
  else{
    //now we have accepted a set of masses - time
    if ( decayer ) {
      decayer->makeDecay(p);
      //p->printTree(); //for debugging
    }
    else{
      p->_rhoBackward.SetDiag(p->getSpinStates());
      
    }
  }  
  _isDecayed=true;
  return;  
}

Referenced by decay(), EvtJetSet::decay(), EvtLunda::decay(), EvtLundCharm::decay(), EvtPythia::decay(), EvtTauola::decay(), EvtGen::generateDecay(), EvtDecayAmp::makeDecay(), EvtDecayIncoherent::makeDecay(), and EvtDecayProb::makeDecay().

◆ decay() [2/2]

void EvtParticle::decay ( )

Decay particle

◆ decayProb() [1/2]

double * EvtParticle::decayProb ( )

inline

Definition at line 390 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

390{return _decayProb;}

◆ decayProb() [2/2]

double * EvtParticle::decayProb ( )

inline

Definition at line 390 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

390{return _decayProb;}

◆ deleteDaughters() [1/2]

void EvtParticle::deleteDaughters ( bool keepChannel = false )

Definition at line 540 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                                 {
  int i;
 
  for(i=0;i<_ndaug;i++){
    _daug[i]->deleteTree();
  }
  
  _ndaug=0;
  //if ( keepChannel ) report(INFO,"EvtGen") << "keeping \n";
  if ( !keepChannel) _channel=-10;
  //_t=0.0;
  //_genlifetime=1;
  _first=1;
  _isInit=false;
  //report(INFO,"EvtGen") << "calling deletedaughters " << EvtPDL::name(this->getId()) <<endl;
}

Referenced by EvtJetSet::decay(), EvtLunda::decay(), EvtLundCharm::decay(), EvtPythia::decay(), EvtTauola::decay(), EvtVSSBMixCPT::decay(), deleteTree(), EvtConExc::gamHXSection(), and initializePhaseSpace().

◆ deleteDaughters() [2/2]

void EvtParticle::deleteDaughters ( bool keepChannel = false )

◆ deleteTree() [1/2]

void EvtParticle::deleteTree ( )

Delete a decay chain

Definition at line 557 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                            {
 
  this->deleteDaughters();
  
  delete this;
  
}

Referenced by EvtbTosllAmp::CalcMaxProb(), EvtSemiLeptonicAmp::CalcMaxProb(), deleteDaughters(), EvtGen::generateDecay(), EvtGen::generateEvent(), EvtPsi3Sdecay::PHSPDecay(), and EvtConExc::~EvtConExc().

◆ deleteTree() [2/2]

void EvtParticle::deleteTree ( )

Delete a decay chain

◆ dumpTree() [1/2]

void EvtParticle::dumpTree ( ) const

Definition at line 978 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                 {  //pingrg, Mar. 25,2008
  
  report(INFO,"EvtGen") << "This is the current decay chain to dump"<<endl;
  report(INFO,"") << "This top particle is "<<
    EvtPDL::name(_id).c_str()<<endl;  
 
  this->dumpTreeRec(0,0);
  report(INFO,"EvtGen") << "End of decay chain."<<endl;
 
}

◆ dumpTree() [2/2]

void EvtParticle::dumpTree ( ) const

◆ dumpTreeRec() [1/2]

void EvtParticle::dumpTreeRec	(	int	level,
		int	dj
	)		const

Definition at line 948 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                                    {  //pingrg, Mar. 25,2008
 
  int newlevel,i;
  newlevel = level +1;
 
  
  if (_ndaug!=0) {
 
    int Ids= EvtPDL::getStdHep(_id); 
    std::string c1,cid;
    if(Ids>0) c1="p";
    if(Ids<0) {
      c1="m";Ids=-1*Ids;
     }   
    
    cid=c1+IntToStr(Ids);
 
    report(INFO,"") <<newlevel<<" "<< cid<<" "<<dj;
    report(INFO,"") << " =  ";
    
    int Nchannel=this->getChannel()+1;
    report(INFO,"") <<Nchannel<<endl;
 
    for(i=0;i<_ndaug;i++){
       _daug[i]->dumpTreeRec(newlevel,i);
    }
  }
}

Referenced by dumpTree(), and dumpTreeRec().

◆ dumpTreeRec() [2/2]

void EvtParticle::dumpTreeRec	(	int	level,
		int	dj
	)		const

◆ eps() [1/2]

EvtVector4C EvtParticle::eps ( int i ) const

virtual

Returns polarization vector in the particles own restframe.

Reimplemented in EvtVectorParticle, and EvtVectorParticle.

Definition at line 576 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                        {
  EvtVector4C temp;
  printParticle();
  report(ERROR,"EvtGen") << "and you have asked for the:"<<i
             <<"th polarization vector."
             <<" I.e. you thought it was a"
             <<" vector particle!" << endl;
  ::abort();
  return temp;
}

Referenced by EvtEDM::decay(), EvtJpipi::decay(), EvtOmegaDalitz::decay(), EvtSVPCP::decay(), EvtSVPHelAmp::decay(), EvtVll::decay(), EvtVPHOtoVISR::decay(), EvtVPHOtoVISRHi::decay(), EvtVSPPwave::decay(), EvtVSS::decay(), EvtVSSBMixCPT::decay(), EvtVSSMix::decay(), EvtVVP::decay(), EvtVVpipi::decay(), EvtVVPIPI_WEIGHTED::decay(), EvtVVSPwave::decay(), and EvtSVVHelAmp::SVVHel().

◆ eps() [2/2]

virtual EvtVector4C EvtParticle::eps ( int i ) const

virtual

Returns polarization vector in the particles own restframe.

Reimplemented in EvtVectorParticle, and EvtVectorParticle.

◆ epsParent() [1/2]

EvtVector4C EvtParticle::epsParent ( int i ) const

virtual

Returns polarization vector in the parents restframe.

Reimplemented in EvtVectorParticle, and EvtVectorParticle.

Definition at line 565 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                              {
  EvtVector4C temp;
  printParticle();
  report(ERROR,"EvtGen") << "and you have asked for the:"<<i
             <<"th polarization vector."
             <<" I.e. you thought it was a"
             <<" vector particle!" << endl;
  ::abort();
  return temp;
}

Referenced by EvtbTosllVectorAmp::CalcAmp(), EvtSemiLeptonicVectorAmp::CalcAmp(), EvtBHadronic::decay(), EvtKstarnunu::decay(), EvtSSDCP::decay(), EvtVectorIsr::decay(), EvtVPHOtoVISR::decay(), EvtVPHOtoVISRHi::decay(), and EvtVVPIPI_WEIGHTED::decay().

◆ epsParent() [2/2]

virtual EvtVector4C EvtParticle::epsParent ( int i ) const

virtual

Returns polarization vector in the parents restframe.

Reimplemented in EvtVectorParticle, and EvtVectorParticle.

◆ epsParentPhoton() [1/2]

EvtVector4C EvtParticle::epsParentPhoton ( int i )

virtual

Returns polarization vector in the parents restframe for a photon.

Reimplemented in EvtPhotonParticle, and EvtPhotonParticle.

Definition at line 587 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                             {
  EvtVector4C temp;
  printParticle();
  report(ERROR,"EvtGen") << "and you have asked for the:"<<i
             <<"th polarization vector of photon."
             <<" I.e. you thought it was a"
             <<" photon particle!" << endl;
  ::abort();
  return temp;
}

Referenced by EvtLNuGamma::decay(), EvtSVPCP::decay(), EvtSVPHelAmp::decay(), EvtVectorIsr::decay(), EvtVSPPwave::decay(), and EvtVVP::decay().

◆ epsParentPhoton() [2/2]

virtual EvtVector4C EvtParticle::epsParentPhoton ( int i )

virtual

Returns polarization vector in the parents restframe for a photon.

Reimplemented in EvtPhotonParticle, and EvtPhotonParticle.

◆ epsPhoton() [1/2]

EvtVector4C EvtParticle::epsPhoton ( int i )

virtual

Returns polarization vector in the particles own restframe for a photon.

Reimplemented in EvtPhotonParticle, and EvtPhotonParticle.

Definition at line 598 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                       {
  EvtVector4C temp;
  printParticle();
  report(ERROR,"EvtGen") << "and you have asked for the:"<<i
             <<"th polarization vector of a photon."
             <<" I.e. you thought it was a"
             <<" photon particle!" << endl;
  ::abort();
  return temp;
}

◆ epsPhoton() [2/2]

virtual EvtVector4C EvtParticle::epsPhoton ( int i )

virtual

Returns polarization vector in the particles own restframe for a photon.

Reimplemented in EvtPhotonParticle, and EvtPhotonParticle.

◆ epsTensor() [1/2]

EvtTensor4C EvtParticle::epsTensor ( int i ) const

virtual

Returns tensor in the particles own restframe for a spin 2 particle.

Reimplemented in EvtTensorParticle, and EvtTensorParticle.

Definition at line 670 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                              {
  int temp;
  temp = i;
  EvtTensor4C tempC; 
  printParticle();
  report(ERROR,"EvtGen") << "and you have asked for the:"<<i
             <<"th tensor."
             <<" I.e. you thought it was a"
             <<" Tensor particle!" << endl;
  ::abort();
  return tempC;
}

Referenced by EvtTSS::decay(), and EvtTVSPwave::decay().

◆ epsTensor() [2/2]

virtual EvtTensor4C EvtParticle::epsTensor ( int i ) const

virtual

Returns tensor in the particles own restframe for a spin 2 particle.

Reimplemented in EvtTensorParticle, and EvtTensorParticle.

◆ epsTensorParent() [1/2]

EvtTensor4C EvtParticle::epsTensorParent ( int i ) const

virtual

Returns tensor in the parents restframe for a spin 2 particle.

Reimplemented in EvtTensorParticle, and EvtTensorParticle.

Definition at line 657 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                                    {
  int temp;
  temp = i;
  EvtTensor4C tempC; 
  printParticle();
  report(ERROR,"EvtGen") << "and you have asked for the:"<<i
             <<"th tensor."
             <<" I.e. you thought it was a"
             <<" Tensor particle!" << endl;
  ::abort();
  return tempC;
}

Referenced by EvtSemiLeptonicTensorAmp::CalcAmp(), EvtBHadronic::decay(), EvtSSDCP::decay(), EvtSTS::decay(), and EvtSTSCP::decay().

◆ epsTensorParent() [2/2]

virtual EvtTensor4C EvtParticle::epsTensorParent ( int i ) const

virtual

Returns tensor in the parents restframe for a spin 2 particle.

Reimplemented in EvtTensorParticle, and EvtTensorParticle.

◆ firstornot() [1/2]

int EvtParticle::firstornot ( ) const

Used internally to decide if first time particle is decayed.

Definition at line 111 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

111{ return _first;}

Referenced by EvtDecayBase::findMasses().

◆ firstornot() [2/2]

int EvtParticle::firstornot ( ) const

Used internally to decide if first time particle is decayed.

◆ generateMassTree() [1/2]

void EvtParticle::generateMassTree ( )

Definition at line 461 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                   {
  double massProb=1.;
  double ranNum=2.;
  int counter=0;
  EvtParticle *p=this;
  while (massProb<ranNum) {
    //check it out the first time.
    p->initDecay();
    //report(INFO,"EvtGen") << "calling massProb \n";
    massProb=p->compMassProb();
    ranNum=EvtRandom::Flat();
   // report(INFO,"EvtGen") << "end of iter " << massProb <<endl;
    counter++; 
 
    if ( counter > 10000 ) {
      if ( counter == 10001 ) {
    report(INFO,"EvtGen") << "Too many iterations to determine the mass tree. Parent mass= "<< p->mass() << _p<<" " << massProb <<endl;
    p->printTree();
    report(INFO,"EvtGen") << "will take next combo with non-zero likelihood\n"; 
      }
      if ( massProb>0. ) massProb=2.0;
      if ( counter > 20000 ) {
    // one last try - take the minimum masses
    p->initDecay(true);
    p->printTree();
    massProb=p->compMassProb();
    if ( massProb>0. ) {
      massProb=2.0;
      report(INFO,"EvtGen") << "Taking the minimum mass of all particles in the chain\n";
    }
    else {
      report(INFO,"EvtGen") << "Sorry, no luck finding a valid set of masses.  This may be a pathological combo\n";
      std::cout<<EvtPDL::name(p->getId())<<": Parent mass "<<p->getP4().mass()<<" with momentum "<<p->getP4()<<std::endl;
      if(EvtGlobalSet::ConExcPythia){EvtGlobalSet::ConExcPythia=0;return;}else{abort();}
      //assert(0);
    }
      }
    }
  }
 // report(INFO,"EvtGen") << counter << endl;
 // p->printTree();
}

Referenced by decay(), EvtBtoKD3P::decay(), and initializePhaseSpace().

◆ generateMassTree() [2/2]

void EvtParticle::generateMassTree ( )

◆ get4Pos() [1/2]

EvtVector4R EvtParticle::get4Pos ( )

Returns the 4position of the particle in the lab frame.

Definition at line 706 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                 {
 
  EvtVector4R temp,mom;
  EvtParticle *ptemp;
  
  temp.set(0.0,0.0,0.0,0.0);
  ptemp=getParent();
 
  if (ptemp==0) return temp;
 
  temp=(ptemp->_t/ptemp->mass())*(ptemp->getP4());
 
  while (ptemp->getParent()!=0) {
    ptemp=ptemp->getParent();
    mom=ptemp->getP4();
    temp=boostTo(temp,mom);
    temp=temp+(ptemp->_t/ptemp->mass())*(ptemp->getP4());
  } 
  
  return temp;
}

Referenced by makeStdHep().

◆ get4Pos() [2/2]

EvtVector4R EvtParticle::get4Pos ( )

Returns the 4position of the particle in the lab frame.

◆ getChannel() [1/2]

int EvtParticle::getChannel ( ) const

Retunrs the decay channel.

Definition at line 123 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

123{ return _channel;}

Referenced by EvtJscont::decay(), EvtPycont::decay(), EvtPyGaGa::decay(), dumpTreeRec(), EvtParticleDecayList::getDecayModel(), EvtDecayAmp::makeDecay(), EvtDecayProb::makeDecay(), and EvtDecayTag::makeTag().

◆ getChannel() [2/2]

int EvtParticle::getChannel ( ) const

Retunrs the decay channel.

◆ getDaug() [1/2]

EvtParticle * EvtParticle::getDaug ( int i )

Get pointer the the i:th daugther.

Definition at line 85 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

85{ return _daug[i]; }

◆ getDaug() [2/2]

EvtParticle * EvtParticle::getDaug ( int i )

Get pointer the the i:th daugther.

◆ getGeneratorFlag() [1/2]

int EvtParticle::getGeneratorFlag ( )

inline

get generator information; pingrg-2011-1-6

Definition at line 146 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

146{ return _generatorFlag ;}

◆ getGeneratorFlag() [2/2]

int EvtParticle::getGeneratorFlag ( )

inline

get generator information; pingrg-2011-1-6

Definition at line 146 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

146{ return _generatorFlag ;}

◆ getId() [1/2]

EvtId EvtParticle::getId ( ) const

Returns Id of particle.

Definition at line 113 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

113{ return _id;}

◆ getId() [2/2]

EvtId EvtParticle::getId ( ) const

Returns Id of particle.

◆ getInclusiveMode() [1/2]

int EvtParticle::getInclusiveMode ( )

inline

Definition at line 394 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

394{return _inclusiveMode;}

◆ getInclusiveMode() [2/2]

int EvtParticle::getInclusiveMode ( )

inline

Definition at line 394 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

394{return _inclusiveMode;}

◆ getIntFlag() [1/2]

std::vector< int > EvtParticle::getIntFlag ( )

inline

get int flag for ConExc: pingrg-2015-2-7

Definition at line 158 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

                              {
    return _intFlag;
  }

Referenced by EvtDecayTag::getModeTag().

◆ getIntFlag() [2/2]

std::vector< int > EvtParticle::getIntFlag ( )

inline

get int flag for ConExc: pingrg-2015-2-7

Definition at line 158 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

                              {
    return _intFlag;
  }

◆ getLifetime() [1/2]

double EvtParticle::getLifetime ( )

Returns the lifetime.

Definition at line 99 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                               {
 
  return _t;
}

Referenced by EvtVSSBMixCPT::decay(), EvtVSSMix::decay(), and EvtCPUtil::OtherB().

◆ getLifetime() [2/2]

double EvtParticle::getLifetime ( )

Returns the lifetime.

◆ getNDaug() [1/2]

int EvtParticle::getNDaug ( ) const

Returns number of daugthers.

Definition at line 125 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

125{ return _ndaug;}

Referenced by EvtRexc::angularSampling(), EvtConExc::angularSampling(), EvtConExc::checkdecay(), EvtPsi3Sdecay::choseDecay(), compMassProb(), EvtSecondary::createSecondary(), EvtBtoKD3P::decay(), EvtDIY::decay(), EvtEtap2gpipi::decay(), EvtEtap2gpipiB::decay(), EvtIntervalDecayAmp< T >::decay(), EvtJetSet::decay(), EvtJscont::decay(), EvtLunda::decay(), EvtLundCharm::decay(), EvtmH2::decay(), EvtOpenCharm::decay(), EvtPycont::decay(), EvtPyGaGa::decay(), EvtPythia::decay(), EvtRhoPi::decay(), EvtSVSCPiso::decay(), EvtTauola::decay(), EvtVPHOtoVISRHi::decay(), EvtVSSBMixCPT::decay(), EvtPHOTOS::doRadCorr(), EvtPsi3Sdecay::EvtPsi3Sdecay(), EvtDecayBase::findMass(), EvtDecayBase::findMaxMass(), EvtFSPick::FSPick(), EvtGen::generateDecay(), EvtParticleDecayList::getDecayModel(), initDecay(), initializePhaseSpace(), EvtDecayAmp::makeDecay(), EvtDecayIncoherent::makeDecay(), EvtDecayProb::makeDecay(), EvtDecayTag::makeTag(), EvtConExc::photonSampling(), EvtPsi3Sdecay::PHSPDecay(), and treeStrRec().

◆ getNDaug() [2/2]

int EvtParticle::getNDaug ( ) const

Returns number of daugthers.

◆ getP4() [1/2]

const EvtVector4R & EvtParticle::getP4 ( ) const

Returns 4momentum in parents restframe.

Definition at line 121 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

121{ return _p;}

Referenced by EvtbTosllScalarAmp::CalcAmp(), EvtbTosllVectorAmp::CalcAmp(), EvtSemiLeptonicBaryonAmp::CalcAmp(), EvtSemiLeptonicTensorAmp::CalcAmp(), EvtSemiLeptonicVectorAmp::CalcAmp(), EvtPsi3Sdecay::choseDecay(), EvtAngSam::decay(), EvtAngSam3::decay(), EvtAngSamX::decay(), EvtAV2GV::decay(), EvtBHadronic::decay(), EvtBsquark::decay(), EvtBTo4piCP::decay(), EvtBtoKD3P::decay(), EvtCalHelAmp::decay(), EvtChi0BB1::decay(), EvtChi0BB2::decay(), EvtChi1BB1::decay(), EvtChi1BB2::decay(), EvtChi2BB1::decay(), EvtChi2BB2::decay(), EvtConExc::decay(), EvtD0mixDalitz::decay(), EvtD0ToKpipi0pi0::decay(), EvtD0ToKpipipi::decay(), EvtDDalitz::decay(), EvtDeBD::decay(), EvtDIY::decay(), EvtDToKpienu::decay(), EvtDToKSpipipi::decay(), EvtDTopipienu::decay(), EvtEDM::decay(), EvtEtaDalitz::decay(), EvtEtap2gpipi::decay(), EvtEtap2gpipiB::decay(), EvtFlatQ2::decay(), EvtHAngSam3::decay(), EvtHelPPJ::decay(), EvtHypWK::decay(), EvtJ2BB1::decay(), EvtJ2BB2::decay(), EvtJ2BB3::decay(), EvtJPE::decay(), EvtJpipi::decay(), EvtJTO3P::decay(), EvtKstarnunu::decay(), EvtKstarstargamma::decay(), EvtLNuGamma::decay(), EvtMBody3::decay(), EvtOmegaDalitz::decay(), EvtOpenCharm::decay(), EvtP2GC0::decay(), EvtP2GC1::decay(), EvtP2GC2::decay(), EvtPhiDalitz::decay(), EvtPhsp::decay(), EvtPi0Dalitz::decay(), EvtRexc::decay(), EvtRhoPi::decay(), EvtSPL::decay(), EvtSSDCP::decay(), EvtSTS::decay(), EvtSTSCP::decay(), EvtSVPCP::decay(), EvtSVPHelAmp::decay(), EvtSVS::decay(), EvtSVSCP::decay(), EvtSVSCPLH::decay(), EvtSVSNONCPEIGEN::decay(), EvtT2GV::decay(), EvtTauGamMu::decay(), EvtTauHadnu::decay(), EvtTauScalarnu::decay(), EvtTSS::decay(), EvtVSPPwave::decay(), EvtVSS::decay(), EvtVSSBMixCPT::decay(), EvtVSSMix::decay(), EvtVVpipi::decay(), EvtVVPIPI_WEIGHTED::decay(), EvtVVSPwave::decay(), EvtConExc::difgamXs(), EvtPHOTOS::doRadCorr(), EvtVectorParticle::epsParent(), EvtPhotonParticle::epsParentPhoton(), EvtTensorParticle::epsTensorParent(), EvtEvalHelAmp::evalAmp(), EvtConExc::findMaxXS(), EvtConExc::gamHXSection(), generateMassTree(), get4Pos(), getP4Lab(), EvtDecayAmp::makeDecay(), EvtPsi3Sdecay::PHSPDecay(), EvtConExc::Rad1difXs(), EvtConExc::Rad2difXs(), EvtDiracParticle::rotateToHelicityBasis(), EvtRaritaSchwingerParticle::rotateToHelicityBasis(), and EvtSVVHelAmp::SVVHel().

◆ getP4() [2/2]

const EvtVector4R & EvtParticle::getP4 ( ) const

Returns 4momentum in parents restframe.

◆ getP4Lab() [1/2]

EvtVector4R EvtParticle::getP4Lab ( )

Gets 4vector in the labframe, i.e., the frame in which the root particles momentum is measured.

Definition at line 685 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                  {
  EvtVector4R temp,mom;
  EvtParticle *ptemp;
  
  temp=this->getP4();
  ptemp=this;
  
  while (ptemp->getParent()!=0) {
    ptemp=ptemp->getParent();
    mom=ptemp->getP4();
    temp=boostTo(temp,mom);   
  } 
  return temp;
}

Referenced by EvtRexc::angularSampling(), EvtConExc::angularSampling(), EvtAngH2::decay(), EvtBody3::decay(), EvtDeBD::decay(), EvtDIY::decay(), EvtEtap2gpipi::decay(), EvtEtap2gpipiB::decay(), EvtLambdaP_BarGamma::decay(), EvtMassH1::decay(), EvtMassH2::decay(), EvtMBody3::decay(), EvtmH2::decay(), EvtmPhsp::decay(), EvtMultibody::decay(), EvtRhoPi::decay(), EvtSingleParticle::decay(), EvtSingleParticle2::decay(), EvtSPL::decay(), makeStdHep(), and EvtConExc::photonSampling().

◆ getP4Lab() [2/2]

EvtVector4R EvtParticle::getP4Lab ( )

Gets 4vector in the labframe, i.e., the frame in which the root particles momentum is measured.

◆ getP4Restframe() [1/2]

EvtVector4R EvtParticle::getP4Restframe ( )

Gets 4vector in the particles restframe, i.e. this functiont will return (m,0,0,0)

Definition at line 700 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                        {
 
  return EvtVector4R(mass(),0.0,0.0,0.0);
 
}

Referenced by EvtKstarstargamma::decay(), EvtSSDCP::decay(), EvtSVSCPLH::decay(), EvtVPHOtoVISR::decay(), and EvtVPHOtoVISRHi::decay().

◆ getP4Restframe() [2/2]

EvtVector4R EvtParticle::getP4Restframe ( )

Gets 4vector in the particles restframe, i.e. this functiont will return (m,0,0,0)

◆ getParent() [1/2]

EvtParticle * EvtParticle::getParent ( )

Returns pointer to parent particle.

Definition at line 87 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

87{ return _parent;}

Referenced by compMassProb(), EvtDDalitz::decay(), EvtDecayBase::findMass(), EvtDecayBase::findMaxMass(), get4Pos(), getP4Lab(), initDecay(), EvtDecayAmp::makeDecay(), and EvtCPUtil::OtherB().

◆ getParent() [2/2]

EvtParticle * EvtParticle::getParent ( )

Returns pointer to parent particle.

◆ getSpinDensityBackward() [1/2]

EvtSpinDensity EvtParticle::getSpinDensityBackward ( )

inline

Get backward spin density matrix.

Definition at line 357 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

357{return _rhoBackward;}

Referenced by EvtDecayAmp::makeDecay().

◆ getSpinDensityBackward() [2/2]

EvtSpinDensity EvtParticle::getSpinDensityBackward ( )

inline

Get backward spin density matrix.

Definition at line 357 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

357{return _rhoBackward;}

◆ getSpinDensityForward() [1/2]

EvtSpinDensity EvtParticle::getSpinDensityForward ( )

inline

Get forward spin density matrix.

Definition at line 347 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

347{return _rhoForward;}

Referenced by EvtDecayAmp::makeDecay(), and EvtDecayIncoherent::makeDecay().

◆ getSpinDensityForward() [2/2]

EvtSpinDensity EvtParticle::getSpinDensityForward ( )

inline

Get forward spin density matrix.

Definition at line 347 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

347{return _rhoForward;}

◆ getSpinStates() [1/2]

int EvtParticle::getSpinStates ( ) const

Returns number of spin states of the particle.

Definition at line 118 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

119 { return EvtSpinType::getSpinStates(EvtPDL::getSpinType(_id));}

EvtPDL::getSpinType

static EvtSpinType::spintype getSpinType(EvtId i)

Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtPDL.hh:61

EvtSpinType::getSpinStates

static int getSpinStates(spintype stype)

Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtSpinType.hh:64

Referenced by EvtbTosllAmp::CalcMaxProb(), EvtSemiLeptonicAmp::CalcMaxProb(), decay(), EvtDecayIncoherent::makeDecay(), EvtDecayProb::makeDecay(), setDiagonalSpinDensity(), setPolarizedSpinDensity(), and setVectorSpinDensity().

◆ getSpinStates() [2/2]

int EvtParticle::getSpinStates ( ) const

Returns number of spin states of the particle.

◆ getSpinType() [1/2]

EvtSpinType::spintype EvtParticle::getSpinType ( ) const

Returns particle type.

Definition at line 115 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

116 { return EvtPDL::getSpinType(_id);}

◆ getSpinType() [2/2]

EvtSpinType::spintype EvtParticle::getSpinType ( ) const

Returns particle type.

◆ hasValidP4() [1/2]

bool EvtParticle::hasValidP4 ( )

inline

Definition at line 383 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

383{return _validP4;}

Referenced by EvtDecayBase::findMaxMass(), EvtParticleDecayList::getDecayModel(), and initDecay().

◆ hasValidP4() [2/2]

bool EvtParticle::hasValidP4 ( )

inline

Definition at line 383 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

383{return _validP4;}

◆ init() [1/2]

virtual void EvtParticle::init	(	EvtId	part_n,
		const EvtVector4R &	p4
	)

pure virtual

Initialiaze particle with id and 4momentum.

Implemented in EvtHighSpinParticle, EvtHighSpinParticle, EvtRaritaSchwingerParticle, EvtStringParticle, EvtRaritaSchwingerParticle, EvtStringParticle, EvtScalarParticle, EvtVectorParticle, EvtScalarParticle, EvtVectorParticle, EvtDiracParticle, EvtNeutrinoParticle, EvtPhotonParticle, EvtTensorParticle, EvtDiracParticle, EvtNeutrinoParticle, EvtPhotonParticle, and EvtTensorParticle.

◆ init() [2/2]

virtual void EvtParticle::init	(	EvtId	part_n,
		const EvtVector4R &	p4
	)

pure virtual

Initialiaze particle with id and 4momentum.

Implemented in EvtHighSpinParticle, EvtHighSpinParticle, EvtRaritaSchwingerParticle, EvtStringParticle, EvtRaritaSchwingerParticle, EvtStringParticle, EvtScalarParticle, EvtVectorParticle, EvtScalarParticle, EvtVectorParticle, EvtDiracParticle, EvtNeutrinoParticle, EvtPhotonParticle, EvtTensorParticle, EvtDiracParticle, EvtNeutrinoParticle, EvtPhotonParticle, and EvtTensorParticle.

◆ initDecay() [1/2]

void EvtParticle::initDecay ( bool useMinMass = false )

Definition at line 237 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                           {
 
  EvtParticle* p=this;
  // carefull - the parent mass might be fixed in stone..
  EvtParticle* par=p->getParent();
  double parMass=-1.;
  if ( par != 0 ) {
    if ( par->hasValidP4() ) parMass=par->mass();
    int i;
    for ( i=0;i<par->getNDaug();i++) {
      EvtParticle *tDaug=par->getDaug(i);
      if ( p != tDaug )
    parMass-=EvtPDL::getMinMass(tDaug->getId());
    }
  }
 
  if ( _isInit ) {
    //we have already been here - just reroll the masses!
    if ( _ndaug>0) {
      int ii;
      for(ii=0;ii<_ndaug;ii++){
    if ( _ndaug==1 ||  EvtPDL::getWidth(p->getDaug(ii)->getId()) > 0.0000001)
      p->getDaug(ii)->initDecay(useMinMass);
    else p->getDaug(ii)->setMass(EvtPDL::getMeanMass(p->getDaug(ii)->getId()));
      }
    }
 
    int j;
    EvtId *dauId=0;
    double *dauMasses=0;
    if ( _ndaug > 0) {
      dauId=new EvtId[_ndaug];
      dauMasses=new double[_ndaug];
      for (j=0;j<_ndaug;j++) { 
    dauId[j]=p->getDaug(j)->getId();
    dauMasses[j]=p->getDaug(j)->mass();
      }
    }
    EvtId *parId=0;
    EvtId *othDauId=0;
    EvtParticle *tempPar=p->getParent();
    if (tempPar) {
      parId=new EvtId(tempPar->getId());
      if ( tempPar->getNDaug()==2 ) {
    if ( tempPar->getDaug(0) == this ) othDauId=new EvtId(tempPar->getDaug(1)->getId());
    else othDauId=new EvtId(tempPar->getDaug(0)->getId());
      }
    }
    if ( p->getParent() && _validP4==false ) {
      if ( !useMinMass ) p->setMass(EvtPDL::getRandMass(p->getId(),parId,_ndaug,dauId,othDauId,parMass,dauMasses));
      else p->setMass(EvtPDL::getMinMass(p->getId()));
    }
    if ( parId) delete parId;
    if ( othDauId) delete othDauId;
    if ( dauId) delete [] dauId;
    if ( dauMasses) delete [] dauMasses;
    return;
  }
 
  
  //Will include effects of mixing here
  //added by Lange Jan4,2000
  static EvtId BS0=EvtPDL::getId("B_s0");
  static EvtId BSB=EvtPDL::getId("anti-B_s0");
  static EvtId BD0=EvtPDL::getId("B0");
  static EvtId BDB=EvtPDL::getId("anti-B0");
 
  //only makes sense if there is no parent particle
  if ( (getNDaug()==0 && getParent()==0) && (getId()==BS0||getId()==BSB||getId()==BD0||getId()==BDB)){
    double t;
    int mix;
    EvtCPUtil::incoherentMix(getId(), t, mix);
    setLifetime(t);
    
    if (mix) {
 
      EvtScalarParticle* scalar_part;
    
      scalar_part=new EvtScalarParticle;
      if (getId()==BS0) {
    EvtVector4R p_init(EvtPDL::getMass(BSB),0.0,0.0,0.0);
    scalar_part->init(BSB,p_init);
      }
      else if (getId()==BSB) {
    EvtVector4R p_init(EvtPDL::getMass(BS0),0.0,0.0,0.0);
    scalar_part->init(BS0,p_init);
      }
      else if (getId()==BD0) {
    EvtVector4R p_init(EvtPDL::getMass(BDB),0.0,0.0,0.0);
    scalar_part->init(BDB,p_init);
      }
      else if (getId()==BDB) {
    EvtVector4R p_init(EvtPDL::getMass(BD0),0.0,0.0,0.0);
    scalar_part->init(BD0,p_init);
      }
 
      scalar_part->setLifetime(0);
 
      scalar_part->setDiagonalSpinDensity();      
    
      insertDaugPtr(0,scalar_part);
 
      _ndaug=1;
      _isInit=true;
      p=scalar_part;
      p->initDecay(useMinMass);
      return;
 
 
    }
  }
  if ( _ndaug==1 ) std::cout << "hi " << EvtPDL::name(this->getId()) << std::endl;
 
  EvtDecayBase *decayer;
  decayer = EvtDecayTable::getDecayFunc(p);
 
  if ( decayer ) {
    p->makeDaughters(decayer->nRealDaughters(),decayer->getDaugs());
    //    report(INFO,"EvtGen") << "has found decay " << decayer->nRealDaughters() << endl;
    //then loop over the daughters and init their decay
    int i;
    for(i=0;i<p->getNDaug();i++){
      if ( EvtPDL::getWidth(p->getDaug(i)->getId()) > 0.0000001)
    p->getDaug(i)->initDecay(useMinMass);
      else p->getDaug(i)->setMass(EvtPDL::getMeanMass(p->getDaug(i)->getId()));
      //  report(INFO,"EvtGen") << "has inited " << EvtPDL::name(p->getDaug(i)->getId()) <<
      //    " " <<  EvtPDL::name(p->getId()) << endl;
    }
  }
  //figure masses in separate step...
  //  if ( p->getParent() && _validP4==false ) EvtDecayBase::findMass(p);
  
  int j;
  EvtId *dauId=0;
  double *dauMasses=0;
  int nDaugT=p->getNDaug();
  if ( nDaugT > 0) {
    dauId=new EvtId[nDaugT];
    dauMasses=new double[nDaugT];
    for (j=0;j<nDaugT;j++) { 
      dauId[j]=p->getDaug(j)->getId();
      dauMasses[j]=p->getDaug(j)->mass();
    }
  }
 
  EvtId *parId=0;
  EvtId *othDauId=0;
  EvtParticle *tempPar=p->getParent();
  if (tempPar) {
    parId=new EvtId(tempPar->getId());
    if ( tempPar->getNDaug()==2 ) {
      if ( tempPar->getDaug(0) == this ) othDauId=new EvtId(tempPar->getDaug(1)->getId());
      else othDauId=new EvtId(tempPar->getDaug(0)->getId());
    }
  }
  if ( p->getParent() && p->hasValidP4()==false ) {
    if ( !useMinMass ) p->setMass(EvtPDL::getRandMass(p->getId(),parId,p->getNDaug(),dauId,othDauId,parMass,dauMasses));
    else p->setMass(EvtPDL::getMinMass(p->getId()));
  }
  if ( parId) delete parId;
  if ( othDauId) delete othDauId;
  if ( dauId) delete [] dauId;
  if ( dauMasses) delete [] dauMasses;
  _isInit=true;
}

Referenced by generateMassTree(), and initDecay().

◆ initDecay() [2/2]

void EvtParticle::initDecay ( bool useMinMass = false )

◆ initializePhaseSpace() [1/2]

double EvtParticle::initializePhaseSpace	(	int	numdaughter,
		EvtId *	daughters,
		double	poleSize = `-1.`,
		int	whichTwo1 = `0`,
		int	whichTwo2 = `1`
	)

Similar to the routine above except that here momentum is generated according to phase space daughters are filled with this momentum.

Definition at line 1071 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                         {
 
  double m_b;
  int i;
  //lange
  //  this->makeDaughters(numdaughter,daughters);
 
  static EvtVector4R p4[100];
  static double mass[100];
 
  m_b = this->mass();
 
  //lange - Jan2,2002 - Need to check to see if the daughters of the parent
  // have changed. If so, delete them and start over.
  //report(INFO,"EvtGen") << "the parent is\n";
  //if ( this->getParent() ) {
  //  if ( this->getParent()->getParent() ) this->getParent()->getParent()->printTree();
    //    this->getParent()->printTree();
  //}
  //report(INFO,"EvtGen") << "and this is\n";
  //if ( this) this->printTree();
  bool resetDaughters=false;
  if ( numdaughter != this->getNDaug() && this->getNDaug() > 0 ) resetDaughters=true;
  if ( numdaughter == this->getNDaug() ) 
    for (i=0; i<numdaughter;i++) {
      if ( this->getDaug(i)->getId() != daughters[i] ) resetDaughters=true;
      //report(INFO,"EvtGen") << this->getDaug(i)->getId() << " " << daughters[i] << endl;
    }
 
  if ( resetDaughters ) {
    //    report(INFO,"EvtGen") << "reseting daughters\n";
    //for (i=0; i<numdaughter;i++) {
    //  report(INFO,"EvtGen") << "reset " <<i<< " "<< EvtPDL::name(this->getDaug(i)->getId()) << " " << EvtPDL::name(daughters[i]) << endl;
    //}
    bool t1=true;
    //but keep the decay channel of the parent.
    this->deleteDaughters(t1);
    this->makeDaughters(numdaughter,daughters);
    this->generateMassTree();
  }
 
  double weight=0.;
  //  EvtDecayBase::findMasses( this, numdaughter, daughters, mass );
  //get the list of masses
  //report(INFO,"EvtGen") << "mpar= " << m_b << " " << this <<endl;
  for (i=0; i<numdaughter;i++) {
    mass[i]=this->getDaug(i)->mass();
    //    report(INFO,"EvtGen") << "mass " << i << " " << mass[i] << " " << this->getDaug(i) << endl;
  }
 
  if ( poleSize<-0.1) {
    EvtGenKine::PhaseSpace( numdaughter, mass, p4, m_b );
    for(i=0;i<numdaughter;i++){
      this->getDaug(i)->init(daughters[i],p4[i]);
    }
 
  }
  else  {
    if ( numdaughter != 3 ) {
      report(ERROR,"EvtGen") << "Only can generate pole phase space "
                 << "distributions for 3 body final states"
                 << endl<<"Will terminate."<<endl;
      ::abort();
    }
    bool ok=false;
    if ( (whichTwo1 == 1 && whichTwo2 == 0 ) ||
     (whichTwo1 == 0 && whichTwo2 == 1 ) ) {
      weight=EvtGenKine::PhaseSpacePole( m_b, mass[0], mass[1], mass[2], 
                      poleSize, p4);
      //report(INFO,"EvtGen") << "here " << weight << " " << poleSize << endl; 
      this->getDaug(0)->init(daughters[0],p4[0]);
      this->getDaug(1)->init(daughters[1],p4[1]);
      this->getDaug(2)->init(daughters[2],p4[2]);
      ok=true;
    }
    if ( (whichTwo1 == 1 && whichTwo2 == 2 ) ||
     (whichTwo1 == 2 && whichTwo2 == 1 ) ) {
      weight=EvtGenKine::PhaseSpacePole( m_b, mass[2], mass[1], mass[0], 
                      poleSize, p4);
      this->getDaug(0)->init(daughters[0],p4[2]);
      this->getDaug(1)->init(daughters[1],p4[1]);
      this->getDaug(2)->init(daughters[2],p4[0]);
      ok=true;
    }
    if ( (whichTwo1 == 0 && whichTwo2 == 2 ) ||
     (whichTwo1 == 2 && whichTwo2 == 0 ) ) {
      weight=EvtGenKine::PhaseSpacePole( m_b, mass[1], mass[0], mass[2], 
                      poleSize, p4);
      this->getDaug(0)->init(daughters[0],p4[1]);
      this->getDaug(1)->init(daughters[1],p4[0]);
      this->getDaug(2)->init(daughters[2],p4[2]);
      ok=true;
    }
    if ( !ok) {
      report(ERROR,"EvtGen") << "Invalid pair of particle to generate a pole dist"
                 << whichTwo1 << " " << whichTwo2
                 << endl<<"Will terminate."<<endl;
      ::abort();
    }
  }
 
  return weight;
}

Referenced by EvtPsi3Sdecay::choseDecay(), EvtAngH2::decay(), EvtAngSam::decay(), EvtAngSam3::decay(), EvtAngSamLab::decay(), EvtAngSamX::decay(), EvtAV2GV::decay(), EvtBHadronic::decay(), EvtBody3::decay(), EvtBsquark::decay(), EvtBto2piCPiso::decay(), EvtBTo4piCP::decay(), EvtBtoKD3P::decay(), EvtBtoKpiCPiso::decay(), EvtbTosllAli::decay(), EvtbTosllBall::decay(), EvtCalHelAmp::decay(), EvtChi0BB1::decay(), EvtChi0BB2::decay(), EvtChi1BB1::decay(), EvtChi1BB2::decay(), EvtChi2BB1::decay(), EvtChi2BB2::decay(), EvtConExc::decay(), EvtD0mixDalitz::decay(), EvtD0ToKpipi0pi0::decay(), EvtD0ToKpipipi::decay(), EvtDDalitz::decay(), EvtDeBD::decay(), EvtDIY::decay(), EvtDMix::decay(), EvtDToKpienu::decay(), EvtDToKSpipipi::decay(), EvtDTopipienu::decay(), EvtEDM::decay(), EvtEtaDalitz::decay(), EvtEtap2gpipi::decay(), EvtEtap2gpipiB::decay(), EvtFlatQ2::decay(), EvtHAngSam3::decay(), EvtHelAmp::decay(), EvtHelPPJ::decay(), EvtHQET::decay(), EvtHQET2::decay(), EvtHypNonLepton::decay(), EvtHypWK::decay(), EvtISGW::decay(), EvtISGW2::decay(), EvtJ2BB1::decay(), EvtJ2BB2::decay(), EvtJ2BB3::decay(), EvtJPE::decay(), EvtJpipi::decay(), EvtJTO3P::decay(), EvtKKLambdaC::decay(), EvtKstarnunu::decay(), EvtKstarstargamma::decay(), EvtLambdac2pKpi::decay(), EvtLambdaP_BarGamma::decay(), EvtLNuGamma::decay(), EvtMassH1::decay(), EvtMassH2::decay(), EvtMBody3::decay(), EvtMelikhov::decay(), EvtmH2::decay(), EvtmPhsp::decay(), EvtMultibody::decay(), EvtOmegaDalitz::decay(), EvtP2GC0::decay(), EvtP2GC1::decay(), EvtP2GC2::decay(), EvtPartWave::decay(), EvtPBB1::decay(), EvtPBB2::decay(), EvtPhiDalitz::decay(), EvtPhsp::decay(), EvtPi0Dalitz::decay(), EvtRexc::decay(), EvtRhoPi::decay(), EvtS2GV::decay(), EvtSLBKPole::decay(), EvtSll::decay(), EvtSLN::decay(), EvtSLPole::decay(), EvtSPL::decay(), EvtSSDCP::decay(), EvtSSSCP::decay(), EvtSSSCPpng::decay(), EvtSSSCPT::decay(), EvtSTS::decay(), EvtSTSCP::decay(), EvtSVPCP::decay(), EvtSVPHelAmp::decay(), EvtSVS::decay(), EvtSVSCP::decay(), EvtSVSCPiso::decay(), EvtSVSCPLH::decay(), EvtSVSNONCPEIGEN::decay(), EvtSVVNONCPEIGEN::decay(), EvtT2GV::decay(), EvtTauGamMu::decay(), EvtTauHadnu::decay(), EvtTaulnunu::decay(), EvtTauScalarnu::decay(), EvtTauVectornu::decay(), EvtTrackGen::decay(), EvtTSS::decay(), EvtTVSPwave::decay(), EvtVectorIsr::decay(), EvtVll::decay(), EvtVSPPwave::decay(), EvtVSS::decay(), EvtVSSBMixCPT::decay(), EvtVSSMix::decay(), EvtVub::decay(), EvtVubHybrid::decay(), EvtVubNLO::decay(), EvtVVP::decay(), EvtVVpipi::decay(), EvtVVPIPI_WEIGHTED::decay(), EvtVVSPwave::decay(), EvtLunda::ExclusiveDecay(), EvtConExc::findMaxXS(), EvtConExc::gamHXSection(), EvtPsi3Sdecay::PHSPDecay(), and EvtSVVHelAmp::SVVHel().

◆ initializePhaseSpace() [2/2]

double EvtParticle::initializePhaseSpace	(	int	numdaughter,
		EvtId *	daughters,
		double	poleSize = `-1.`,
		int	whichTwo1 = `0`,
		int	whichTwo2 = `1`
	)

Similar to the routine above except that here momentum is generated according to phase space daughters are filled with this momentum.

◆ insertDaugPtr() [1/2]

void EvtParticle::insertDaugPtr	(	int	idaug,
		EvtParticle *	partptr
	)

inline

Makes partptr the idaug:th daugther.

Definition at line 220 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

220 { _daug[idaug]=partptr;

221 partptr->_parent=this; }

Referenced by initDecay(), and EvtCPUtil::OtherB().

◆ insertDaugPtr() [2/2]

void EvtParticle::insertDaugPtr	(	int	idaug,
		EvtParticle *	partptr
	)

inline

Makes partptr the idaug:th daugther.

Definition at line 220 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

220 { _daug[idaug]=partptr;

221 partptr->_parent=this; }

◆ isDecayed() [1/2]

bool EvtParticle::isDecayed ( )

inline

Definition at line 384 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

384{return _isDecayed;}

◆ isDecayed() [2/2]

bool EvtParticle::isDecayed ( )

inline

Definition at line 384 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

384{return _isDecayed;}

◆ isInitialized() [1/2]

bool EvtParticle::isInitialized ( )

inline

Definition at line 382 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

382{return _isInit;}

Referenced by EvtDecayBase::findMaxMass().

◆ isInitialized() [2/2]

bool EvtParticle::isInitialized ( )

inline

Definition at line 382 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

382{return _isInit;}

◆ makeDaughters() [1/2]

void EvtParticle::makeDaughters	(	int	ndaug,
		EvtId *	id
	)

Creates the daughters in the list of ids and adds them to the parent. Note that momentum is left uninitialized, this is only creation.

Definition at line 1177 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                                         {
 
  int i;
  if ( _channel < 0 ) {
    //report(INFO,"EvtGen") << "setting channel " << EvtPDL::name(this->getId()) << endl;
    setChannel(0);
  }
  EvtParticle* pdaug;  
  if (_ndaug!=0 ){
    if (_ndaug!=ndaugstore){
      report(ERROR,"EvtGen") << "Asking to make a different number of "
                 << "daughters than what was previously created."
                 << endl<<"Will terminate."<<endl;
      ::abort();
    }
  } 
  else{
    for(i=0;i<ndaugstore;i++){
      pdaug=EvtParticleFactory::particleFactory(EvtPDL::getSpinType(id[i]));
      pdaug->setId(id[i]);
      pdaug->addDaug(this); 
    }
 
  } //else
} //makeDaughters

◆ makeDaughters() [2/2]

void EvtParticle::makeDaughters	(	int	ndaug,
		EvtId *	id
	)

Creates the daughters in the list of ids and adds them to the parent. Note that momentum is left uninitialized, this is only creation.

◆ makeStdHep() [1/4]

void EvtParticle::makeStdHep ( EvtStdHep & stdhep )

Definition at line 795 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                             {
 
  //first add particle to the stdhep list;
  stdhep.createParticle(getP4Lab(),get4Pos(),-1,-1,
            EvtPDL::getStdHep(getId()));
 
  int i;
  for(i=0;i<_ndaug;i++){
    stdhep.createParticle(_daug[i]->getP4Lab(),_daug[i]->get4Pos(),0,0,
              EvtPDL::getStdHep(_daug[i]->getId()));
  }
 
  for(i=0;i<_ndaug;i++){
    _daug[i]->makeStdHepRec(1+i,1+i,stdhep);
  }
  return;
 
}

◆ makeStdHep() [2/4]

void EvtParticle::makeStdHep ( EvtStdHep & stdhep )

◆ makeStdHep() [3/4]

void EvtParticle::makeStdHep	(	EvtStdHep &	stdhep,
		EvtSecondary &	secondary,
		EvtId *	stable_parent_ihep
	)

Makes stdhep list

Definition at line 759 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                       {
 
  //first add particle to the stdhep list;
  stdhep.createParticle(getP4Lab(),get4Pos(),-1,-1,
            EvtPDL::getStdHep(getId()));
 
  int ii=0;
 
  //lets see if this is a longlived particle and terminate the 
  //list building!
  
  while (list_of_stable[ii]!=EvtId(-1,-1)) {
    if (getId()==list_of_stable[ii]){
      secondary.createSecondary(0,this);
      return;
    }
    ii++;
  }
 
 
 
 
  int i;
  for(i=0;i<_ndaug;i++){
    stdhep.createParticle(_daug[i]->getP4Lab(),_daug[i]->get4Pos(),0,0,
              EvtPDL::getStdHep(_daug[i]->getId()));
  }
 
  for(i=0;i<_ndaug;i++){
    _daug[i]->makeStdHepRec(1+i,1+i,stdhep,secondary,list_of_stable);
  }
  return;
 
}

Referenced by EvtGen::generateDecay(), and EvtGen::generateEvent().

◆ makeStdHep() [4/4]

void EvtParticle::makeStdHep	(	EvtStdHep &	stdhep,
		EvtSecondary &	secondary,
		EvtId *	stable_parent_ihep
	)

Makes stdhep list

◆ mass() [1/2]

double EvtParticle::mass ( ) const

Returns mass of particle.

Definition at line 127 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                               {
 
     return _p.mass();
}

◆ mass() [2/2]

double EvtParticle::mass ( ) const

Returns mass of particle.

◆ nextIter() [1/2]

EvtParticle * EvtParticle::nextIter ( EvtParticle * rootOfTree = 0 )

Iterates over the particles in a decay chain.

Definition at line 729 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                                           {
 
  EvtParticle *bpart;
  EvtParticle *current;
 
  current=this;
  int i;
 
  if (_ndaug!=0) return _daug[0];
 
  do{
    bpart=current->_parent;
    if (bpart==0) return 0;
    i=0;
    while (bpart->_daug[i]!=current) {i++;}
 
    if ( bpart==rootOfTree ) {
      if ( i== bpart->_ndaug-1 ) return 0;
    }
 
    i++;
    current=bpart;
 
  }while(i>=bpart->_ndaug);
 
  return bpart->_daug[i];
 
}

◆ nextIter() [2/2]

EvtParticle * EvtParticle::nextIter ( EvtParticle * rootOfTree = 0 )

Iterates over the particles in a decay chain.

◆ noLifeTime() [1/2]

void EvtParticle::noLifeTime ( )

inline

Definition at line 362 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

362{ _genlifetime=0; }

Referenced by EvtbTosllAmp::CalcMaxProb(), and EvtSemiLeptonicAmp::CalcMaxProb().

◆ noLifeTime() [2/2]

void EvtParticle::noLifeTime ( )

inline

Definition at line 362 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

362{ _genlifetime=0; }

◆ printParticle() [1/2]

void EvtParticle::printParticle ( ) const

Prints information for the particle.

Definition at line 1000 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                      {
 
  switch (EvtPDL::getSpinType(_id)){ 
  case EvtSpinType::SCALAR:
    report(INFO,"EvtGen") << "This is a scalar particle:"<<EvtPDL::name(_id).c_str()<<"\n";
    break;     
  case EvtSpinType::VECTOR:
    report(INFO,"EvtGen") << "This is a vector particle:"<<EvtPDL::name(_id).c_str()<<"\n";
    break;     
  case EvtSpinType::TENSOR:
    report(INFO,"EvtGen") << "This is a tensor particle:"<<EvtPDL::name(_id).c_str()<<"\n";
    break;
  case EvtSpinType::DIRAC:
    report(INFO,"EvtGen") << "This is a dirac particle:"<<EvtPDL::name(_id).c_str()<<"\n";
    break;
  case EvtSpinType::PHOTON:
    report(INFO,"EvtGen") << "This is a photon:"<<EvtPDL::name(_id).c_str()<<"\n";
    break;
  case EvtSpinType::NEUTRINO:
    report(INFO,"EvtGen") << "This is a neutrino:"<<EvtPDL::name(_id).c_str()<<"\n";
    break;
  case EvtSpinType::RARITASCHWINGER:
    report(INFO,"EvtGen") << "This is a raritaschwinger:"<<EvtPDL::name(_id).c_str()<<"\n";
    break;
  case EvtSpinType::STRING:
    report(INFO,"EvtGen") << "This is a string:"<<EvtPDL::name(_id).c_str()<<"\n";
    break;
  default:
    report(INFO,"EvtGen") <<"Unknown particle type in EvtParticle::printParticle()"<<endl;
    break;
  }
  report(INFO,"EvtGen") << "Number of daughters:"<<_ndaug<<"\n";
 
 
}

Referenced by eps(), epsParent(), epsParentPhoton(), epsPhoton(), epsTensor(), epsTensorParent(), sp(), spNeutrino(), spParent(), and spParentNeutrino().

◆ printParticle() [2/2]

void EvtParticle::printParticle ( ) const

Prints information for the particle.

◆ printTree() [1/2]

void EvtParticle::printTree ( ) const

Prints out the particle "tree" of a given particle. The tree consists of all daughters (and their daughters, etc) and their properties.

Definition at line 897 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                  {
  
  report(INFO,"EvtGen") << "This is the current decay chain"<<endl;
  report(INFO,"") << "This top particle is "<<
    EvtPDL::name(_id).c_str()<<endl;  
 
  this->printTreeRec(0);
  report(INFO,"EvtGen") << "End of decay chain."<<endl;
 
}

Referenced by EvtConExc::decay(), EvtDecayBase::findMass(), and generateMassTree().

◆ printTree() [2/2]

void EvtParticle::printTree ( ) const

Prints out the particle "tree" of a given particle. The tree consists of all daughters (and their daughters, etc) and their properties.

◆ printTreeRec() [1/2]

void EvtParticle::printTreeRec ( int level ) const

Definition at line 870 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                              {
 
  int newlevel,i;
  newlevel = level +1;
 
  
  if (_ndaug!=0) {
    if ( level > 0 ) {
      for (i=0;i<(5*level);i++) {
    report(INFO,"") <<" ";
      }
    }
    report(INFO,"") << EvtPDL::name(_id).c_str();  
    report(INFO,"") << " -> ";
    for(i=0;i<_ndaug;i++){
      report(INFO,"") << EvtPDL::name(_daug[i]->getId()).c_str()<<" ";
    }
    for(i=0;i<_ndaug;i++){
      report(INFO,"") << _daug[i]->mass()<<" ";
    }
    report(INFO,"")<<endl;
    for(i=0;i<_ndaug;i++){
      _daug[i]->printTreeRec(newlevel);
    }
  }
}

Referenced by printTree(), and printTreeRec().

◆ printTreeRec() [2/2]

void EvtParticle::printTreeRec ( int level ) const

◆ resetFirstOrNot() [1/2]

void EvtParticle::resetFirstOrNot ( )

Definition at line 77 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                  {
  _first=1;
}

Referenced by EvtGen::generateDecay().

◆ resetFirstOrNot() [2/2]

void EvtParticle::resetFirstOrNot ( )

◆ resetNDaug() [1/2]

void EvtParticle::resetNDaug ( )

inline

Definition at line 269 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

269{_ndaug=0; return;}

Referenced by EvtOpenCharm::decay(), and EvtGen::generateDecay().

◆ resetNDaug() [2/2]

void EvtParticle::resetNDaug ( )

inline

Definition at line 269 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

269{_ndaug=0; return;}

◆ rotateToHelicityBasis() [1/4]

virtual EvtSpinDensity EvtParticle::rotateToHelicityBasis ( ) const

pure virtual

Returns a rotation matrix need to rotate the basis state to the helicity basis. The EvtSpinDensity matrix is just use as a matrix here. This function is to be implemented in each derived class.

Implemented in EvtDiracParticle, EvtHighSpinParticle, EvtNeutrinoParticle, EvtPhotonParticle, EvtRaritaSchwingerParticle, EvtScalarParticle, EvtStringParticle, EvtTensorParticle, EvtVectorParticle, EvtDiracParticle, EvtHighSpinParticle, EvtNeutrinoParticle, EvtPhotonParticle, EvtRaritaSchwingerParticle, EvtScalarParticle, EvtStringParticle, EvtTensorParticle, and EvtVectorParticle.

Referenced by EvtMultibody::decay(), and setSpinDensityForwardHelicityBasis().

◆ rotateToHelicityBasis() [2/4]

virtual EvtSpinDensity EvtParticle::rotateToHelicityBasis ( ) const

pure virtual

Returns a rotation matrix need to rotate the basis state to the helicity basis. The EvtSpinDensity matrix is just use as a matrix here. This function is to be implemented in each derived class.

Implemented in EvtDiracParticle, EvtHighSpinParticle, EvtNeutrinoParticle, EvtPhotonParticle, EvtRaritaSchwingerParticle, EvtScalarParticle, EvtStringParticle, EvtTensorParticle, EvtVectorParticle, EvtDiracParticle, EvtHighSpinParticle, EvtNeutrinoParticle, EvtPhotonParticle, EvtRaritaSchwingerParticle, EvtScalarParticle, EvtStringParticle, EvtTensorParticle, and EvtVectorParticle.

◆ rotateToHelicityBasis() [3/4]

virtual EvtSpinDensity EvtParticle::rotateToHelicityBasis	(	double	alpha,
		double	beta,
		double	gamma
	)		const

pure virtual

Implemented in EvtDiracParticle, EvtHighSpinParticle, EvtNeutrinoParticle, EvtPhotonParticle, EvtRaritaSchwingerParticle, EvtScalarParticle, EvtStringParticle, EvtTensorParticle, EvtVectorParticle, EvtDiracParticle, EvtHighSpinParticle, EvtNeutrinoParticle, EvtPhotonParticle, EvtRaritaSchwingerParticle, EvtScalarParticle, EvtStringParticle, EvtTensorParticle, and EvtVectorParticle.

◆ rotateToHelicityBasis() [4/4]

virtual EvtSpinDensity EvtParticle::rotateToHelicityBasis	(	double	alpha,
		double	beta,
		double	gamma
	)		const

pure virtual

Implemented in EvtDiracParticle, EvtHighSpinParticle, EvtNeutrinoParticle, EvtPhotonParticle, EvtRaritaSchwingerParticle, EvtScalarParticle, EvtStringParticle, EvtTensorParticle, EvtVectorParticle, EvtDiracParticle, EvtHighSpinParticle, EvtNeutrinoParticle, EvtPhotonParticle, EvtRaritaSchwingerParticle, EvtScalarParticle, EvtStringParticle, EvtTensorParticle, and EvtVectorParticle.

◆ setChannel() [1/2]

void EvtParticle::setChannel ( int i )

Should only be used internally.

Definition at line 81 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                    { 
  _channel=i;
}

Referenced by EvtParticleDecayList::getDecayModel(), and makeDaughters().

◆ setChannel() [2/2]

void EvtParticle::setChannel ( int i )

Should only be used internally.

◆ setDecayProb() [1/2]

void EvtParticle::setDecayProb ( double p )

Definition at line 1204 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                          {
 
  if ( _decayProb == 0 )  _decayProb=new double;
  *_decayProb=prob;
}

Referenced by EvtDecayAmp::makeDecay(), EvtDecayIncoherent::makeDecay(), and EvtDecayProb::makeDecay().

◆ setDecayProb() [2/2]

void EvtParticle::setDecayProb ( double p )

◆ setDiagonalSpinDensity() [1/2]

void EvtParticle::setDiagonalSpinDensity ( )

Set diagonal spindensity matrix.

Definition at line 133 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                        {
 
  _rhoForward.SetDiag(getSpinStates());
}

Referenced by EvtbTosllAmp::CalcMaxProb(), EvtSemiLeptonicAmp::CalcMaxProb(), EvtJetSet::decay(), EvtLunda::decay(), EvtLundCharm::decay(), EvtPythia::decay(), EvtTauola::decay(), and initDecay().

◆ setDiagonalSpinDensity() [2/2]

void EvtParticle::setDiagonalSpinDensity ( )

Set diagonal spindensity matrix.

◆ setFirstOrNot() [1/2]

void EvtParticle::setFirstOrNot ( )

Definition at line 74 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                {
  _first=0;
}

Referenced by EvtDecayBase::findMasses().

◆ setFirstOrNot() [2/2]

void EvtParticle::setFirstOrNot ( )

◆ setGeneratorFlag() [1/2]

void EvtParticle::setGeneratorFlag ( int flag )

inline

set generator information; pingrg-2011-1-6

Definition at line 141 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

141{_generatorFlag = flag;}

flag

long int flag

Definition: Eepipi/Eepipi-00-00-06/src/ee2eepp/basesv5.1/f2c.h:39

Referenced by EvtLundCharm::decay(), and EvtOpenCharm::decay().

◆ setGeneratorFlag() [2/2]

void EvtParticle::setGeneratorFlag ( int flag )

inline

set generator information; pingrg-2011-1-6

Definition at line 141 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

141{_generatorFlag = flag;}

◆ setId() [1/2]

void EvtParticle::setId ( EvtId id )

inline

Definition at line 365 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

365{ _id=id;}

Referenced by makeDaughters().

◆ setId() [2/2]

void EvtParticle::setId ( EvtId id )

inline

Definition at line 365 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

365{ _id=id;}

◆ setInclusiveMode() [1/2]

void EvtParticle::setInclusiveMode ( int im )

inline

Definition at line 393 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

393{_inclusiveMode=im ;}

◆ setInclusiveMode() [2/2]

void EvtParticle::setInclusiveMode ( int im )

inline

Definition at line 393 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

393{_inclusiveMode=im ;}

◆ setIntFlag() [1/2]

void EvtParticle::setIntFlag ( std::vector< int > vi )

inline

set int flag for ConExc: pingrg-2015-2-7

Definition at line 151 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

                                    {
    _intFlag=vi;
  }

Referenced by EvtConExc::decay().

◆ setIntFlag() [2/2]

void EvtParticle::setIntFlag ( std::vector< int > vi )

inline

set int flag for ConExc: pingrg-2015-2-7

Definition at line 151 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

                                    {
    _intFlag=vi;
  }

◆ setLifetime() [1/4]

void EvtParticle::setLifetime ( )

Generate lifetime according to pure exponential.

Definition at line 93 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                             {
  if (_genlifetime){
    _t=-log(EvtRandom::Flat())*EvtPDL::getctau(getId());
  }
}

Referenced by EvtHighSpinParticle::init(), EvtRaritaSchwingerParticle::init(), EvtDiracParticle::init(), EvtNeutrinoParticle::init(), EvtTensorParticle::init(), EvtPhotonParticle::init(), EvtScalarParticle::init(), EvtVectorParticle::init(), and initDecay().

◆ setLifetime() [2/4]

void EvtParticle::setLifetime ( )

Generate lifetime according to pure exponential.

◆ setLifetime() [3/4]

void EvtParticle::setLifetime ( double tau )

Set lifetime of the particle in parents restframe.

Definition at line 89 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                       {
  _t=tau;
}

Referenced by EvtD0mixDalitz::decay(), EvtDMix::decay(), EvtVSSBMixCPT::decay(), EvtVSSMix::decay(), EvtVub::decay(), EvtVubHybrid::decay(), initDecay(), and EvtCPUtil::OtherB().

◆ setLifetime() [4/4]

void EvtParticle::setLifetime ( double tau )

Set lifetime of the particle in parents restframe.

◆ setMass() [1/2]

void EvtParticle::setMass ( double m )

inline

Definition at line 372 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

372{ _p=EvtVector4R(m,0.0,0.0,0.0);}

Referenced by EvtPsi3Sdecay::choseDecay(), EvtConExc::decay(), EvtRexc::decay(), EvtLunda::ExclusiveDecay(), EvtDecayBase::findMass(), EvtConExc::findMaxXS(), EvtConExc::gamHXSection(), initDecay(), and EvtPsi3Sdecay::PHSPDecay().

◆ setMass() [2/2]

void EvtParticle::setMass ( double m )

inline

Definition at line 372 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

372{ _p=EvtVector4R(m,0.0,0.0,0.0);}

◆ setp() [1/4]

void EvtParticle::setp ( const EvtVector4R & p4 )

inlineprotected

Definition at line 399 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

399{ _p =p4; }

◆ setp() [2/4]

void EvtParticle::setp ( const EvtVector4R & p4 )

inlineprotected

Definition at line 399 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

399{ _p =p4; }

◆ setp() [3/4]

void EvtParticle::setp	(	double	e,
		double	px,
		double	py,
		double	pz
	)

inlineprotected

Definition at line 398 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

398{ _p.set(e,px,py,pz); }

Referenced by EvtHighSpinParticle::init(), EvtRaritaSchwingerParticle::init(), EvtStringParticle::init(), EvtDiracParticle::init(), EvtNeutrinoParticle::init(), EvtPhotonParticle::init(), EvtScalarParticle::init(), EvtTensorParticle::init(), and EvtVectorParticle::init().

◆ setp() [4/4]

void EvtParticle::setp	(	double	e,
		double	px,
		double	py,
		double	pz
	)

inlineprotected

Definition at line 398 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

398{ _p.set(e,px,py,pz); }

◆ setP4() [1/2]

void EvtParticle::setP4 ( const EvtVector4R & p4 )

inline

Sets the 4momentum in the parents restframe.

Definition at line 258 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

258{_p=p4;}

Referenced by EvtPHOTOS::doRadCorr().

◆ setP4() [2/2]

void EvtParticle::setP4 ( const EvtVector4R & p4 )

inline

Sets the 4momentum in the parents restframe.

Definition at line 258 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

258{_p=p4;}

◆ setpart_num() [1/2]

void EvtParticle::setpart_num ( EvtId particle_number )

inlineprotected

Definition at line 400 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

  { 
    assert(_channel==-10||
       _id.getId()==particle_number.getId()||
       _id.getId()==-1);
    _id = particle_number; 
  }

Referenced by EvtHighSpinParticle::init(), EvtRaritaSchwingerParticle::init(), EvtStringParticle::init(), EvtDiracParticle::init(), EvtNeutrinoParticle::init(), EvtPhotonParticle::init(), EvtScalarParticle::init(), EvtTensorParticle::init(), and EvtVectorParticle::init().

◆ setpart_num() [2/2]

void EvtParticle::setpart_num ( EvtId particle_number )

inlineprotected

Definition at line 400 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

  { 
    assert(_channel==-10||
       _id.getId()==particle_number.getId()||
       _id.getId()==-1);
    _id = particle_number; 
  }

◆ setPolarizedSpinDensity() [1/2]

void EvtParticle::setPolarizedSpinDensity	(	double	r00,
		double	r11,
		double	r22
	)

Definition at line 157 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                                                         {
 
  if (getSpinStates()!=3) {
    report(ERROR,"EvtGen")<<"Error in EvtParticle::setVectorSpinDensity"<<endl;
    report(ERROR,"EvtGen")<<"spin_states:"<<getSpinStates()<<endl;
    report(ERROR,"EvtGen")<<"particle:"<<EvtPDL::name(_id).c_str()<<endl;
    ::abort();
  }
 
  EvtSpinDensity rho;
 
  //Set helicity +1 and -1 to 1.
  rho.SetDiag(getSpinStates());
  rho.Set(0,0,EvtComplex(r00,0.));
  rho.Set(1,1,EvtComplex(r11,0.));
  rho.Set(2,2,EvtComplex(r22,0.));
  setSpinDensityForwardHelicityBasis(rho);
 
}

◆ setPolarizedSpinDensity() [2/2]

void EvtParticle::setPolarizedSpinDensity	(	double	r00,
		double	r11,
		double	r22
	)

◆ setSpinDensityBackward() [1/2]

void EvtParticle::setSpinDensityBackward ( const EvtSpinDensity & rho )

inline

Set backward spin density matrix.

Definition at line 352 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

352{_rhoBackward=rho;}

Referenced by EvtDecayAmp::makeDecay(), EvtDecayIncoherent::makeDecay(), and EvtDecayProb::makeDecay().

◆ setSpinDensityBackward() [2/2]

void EvtParticle::setSpinDensityBackward ( const EvtSpinDensity & rho )

inline

Set backward spin density matrix.

Definition at line 352 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

352{_rhoBackward=rho;}

◆ setSpinDensityForward() [1/2]

void EvtParticle::setSpinDensityForward ( const EvtSpinDensity & rho )

inline

Set forward spin density matrix.

Definition at line 321 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

321{_rhoForward=rho;}

Referenced by EvtVectorIsr::decay(), EvtDecayAmp::makeDecay(), EvtDecayIncoherent::makeDecay(), EvtDecayProb::makeDecay(), and EvtParticleFactory::particleFactory().

◆ setSpinDensityForward() [2/2]

void EvtParticle::setSpinDensityForward ( const EvtSpinDensity & rho )

inline

Set forward spin density matrix.

Definition at line 321 of file BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

321{_rhoForward=rho;}

◆ setSpinDensityForwardHelicityBasis() [1/4]

void EvtParticle::setSpinDensityForwardHelicityBasis ( const EvtSpinDensity & rho )

Set forward spin density matrix according to the density matrix rho in the helicity amplitude basis.

Definition at line 178 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                                                             {
 
  EvtSpinDensity R=rotateToHelicityBasis();
 
  assert(R.GetDim()==rho.GetDim());
 
  int n=rho.GetDim();
 
  _rhoForward.SetDim(n);
 
  int i,j,k,l;
 
  for(i=0;i<n;i++){
    for(j=0;j<n;j++){
      EvtComplex tmp=0.0;
      for(k=0;k<n;k++){
    for(l=0;l<n;l++){
      tmp+=R.Get(l,i)*rho.Get(l,k)*conj(R.Get(k,j));
    }
      }
      _rhoForward.Set(i,j,tmp);
    }
  }
 
  //report(INFO,"EvtGen") << "_rhoForward:"<<_rhoForward<<endl;
 
}

Referenced by setPolarizedSpinDensity(), and setVectorSpinDensity().

◆ setSpinDensityForwardHelicityBasis() [2/4]

void EvtParticle::setSpinDensityForwardHelicityBasis ( const EvtSpinDensity & rho )

Set forward spin density matrix according to the density matrix rho in the helicity amplitude basis.

◆ setSpinDensityForwardHelicityBasis() [3/4]

void EvtParticle::setSpinDensityForwardHelicityBasis	(	const EvtSpinDensity &	rho,
		double	alpha,
		double	beta,
		double	gamma
	)

Definition at line 206 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                          {
 
  EvtSpinDensity R=rotateToHelicityBasis(alpha,beta,gamma);
 
  assert(R.GetDim()==rho.GetDim());
 
  int n=rho.GetDim();
 
  _rhoForward.SetDim(n);
 
  int i,j,k,l;
 
  for(i=0;i<n;i++){
    for(j=0;j<n;j++){
      EvtComplex tmp=0.0;
      for(k=0;k<n;k++){
    for(l=0;l<n;l++){
      tmp+=R.Get(l,i)*rho.Get(l,k)*conj(R.Get(k,j));
    }
      }
      _rhoForward.Set(i,j,tmp);
    }
  }
 
  //report(INFO,"EvtGen") << "_rhoForward:"<<_rhoForward<<endl;
 
}

◆ setSpinDensityForwardHelicityBasis() [4/4]

void EvtParticle::setSpinDensityForwardHelicityBasis	(	const EvtSpinDensity &	rho,
		double	alpha,
		double	beta,
		double	gamma
	)

◆ setVectorSpinDensity() [1/2]

void EvtParticle::setVectorSpinDensity ( )

Set spindensity matrix for e+e- -> V

Definition at line 138 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                      {
 
  if (getSpinStates()!=3) {
    report(ERROR,"EvtGen")<<"Error in EvtParticle::setVectorSpinDensity"<<endl;
    report(ERROR,"EvtGen")<<"spin_states:"<<getSpinStates()<<endl;
    report(ERROR,"EvtGen")<<"particle:"<<EvtPDL::name(_id).c_str()<<endl;
    ::abort();
  }
 
  EvtSpinDensity rho;
 
  //Set helicity +1 and -1 to 1.
  rho.SetDiag(getSpinStates());
  rho.Set(1,1,EvtComplex(0.0,0.0));
  setSpinDensityForwardHelicityBasis(rho);
 
}

Referenced by EvtGen::generateEvent().

◆ setVectorSpinDensity() [2/2]

void EvtParticle::setVectorSpinDensity ( )

Set spindensity matrix for e+e- -> V

◆ sp() [1/2]

EvtDiracSpinor EvtParticle::sp ( int i ) const

virtual

Returns Dirac spinor in the particles own restframe for a Dirac particle.

Reimplemented in EvtDiracParticle, and EvtDiracParticle.

Definition at line 622 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                          {
  EvtDiracSpinor tempD;
  int temp;
  temp = i;
  printParticle();
  report(ERROR,"EvtGen") << "and you have asked for the:"<<i
             <<"th dirac spinor."
             <<" I.e. you thought it was a"
             <<" Dirac particle!" << endl;
  ::abort();
  return tempD;
}

Referenced by EvtHypNonLepton::calcAmp(), EvtSemiLeptonicBaryonAmp::CalcAmp(), EvtTauHadnu::decay(), EvtTaulnunu::decay(), EvtTauScalarnu::decay(), and EvtTauVectornu::decay().

◆ sp() [2/2]

virtual EvtDiracSpinor EvtParticle::sp ( int ) const

virtual

Returns Dirac spinor in the particles own restframe for a Dirac particle.

Reimplemented in EvtDiracParticle, and EvtDiracParticle.

◆ spNeutrino() [1/2]

EvtDiracSpinor EvtParticle::spNeutrino ( ) const

virtual

Returns Dirac spinor in the particles own restframe for a Neutrino particle.

Reimplemented in EvtNeutrinoParticle, and EvtNeutrinoParticle.

Definition at line 646 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                             {
  EvtDiracSpinor tempD;
  printParticle();
  report(ERROR,"EvtGen") << "and you have asked for the "
             <<"dirac spinor."
             <<" I.e. you thought it was a"
             <<" neutrino particle!" << endl;
  ::abort();
  return tempD;
}

◆ spNeutrino() [2/2]

virtual EvtDiracSpinor EvtParticle::spNeutrino ( ) const

virtual

Returns Dirac spinor in the particles own restframe for a Neutrino particle.

Reimplemented in EvtNeutrinoParticle, and EvtNeutrinoParticle.

◆ spParent() [1/2]

EvtDiracSpinor EvtParticle::spParent ( int i ) const

virtual

Returns Dirac spinor in the parents restframe for a Dirac particle.

Reimplemented in EvtDiracParticle, and EvtDiracParticle.

Definition at line 609 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                                {
  EvtDiracSpinor tempD;
  int temp;
  temp = i;
  printParticle();
  report(ERROR,"EvtGen") << "and you have asked for the:"<<i
             <<"th dirac spinor."
             <<" I.e. you thought it was a"
             <<" Dirac particle!" << endl;
  ::abort();
  return tempD;
}

Referenced by EvtHypNonLepton::calcAmp(), EvtbTosllScalarAmp::CalcAmp(), EvtbTosllVectorAmp::CalcAmp(), EvtSemiLeptonicBaryonAmp::CalcAmp(), EvtSemiLeptonicTensorAmp::CalcAmp(), EvtSemiLeptonicVectorAmp::CalcAmp(), EvtBsquark::decay(), EvtEDM::decay(), EvtLNuGamma::decay(), EvtSll::decay(), EvtSLN::decay(), EvtTaulnunu::decay(), and EvtVll::decay().

◆ spParent() [2/2]

virtual EvtDiracSpinor EvtParticle::spParent ( int ) const

virtual

Returns Dirac spinor in the parents restframe for a Dirac particle.

Reimplemented in EvtDiracParticle, and EvtDiracParticle.

◆ spParentNeutrino() [1/2]

EvtDiracSpinor EvtParticle::spParentNeutrino ( ) const

virtual

Returns Dirac spinor in the parents restframe for a Neutrino particle.

Reimplemented in EvtNeutrinoParticle, and EvtNeutrinoParticle.

Definition at line 635 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                                   {
  EvtDiracSpinor tempD;
  printParticle();
  report(ERROR,"EvtGen") << "and you have asked for the "
             <<"dirac spinor."
             <<" I.e. you thought it was a"
             <<" neutrino particle!" << endl;
  ::abort();
  return tempD;
}

Referenced by EvtSemiLeptonicBaryonAmp::CalcAmp(), EvtSemiLeptonicTensorAmp::CalcAmp(), EvtSemiLeptonicVectorAmp::CalcAmp(), EvtKstarnunu::decay(), EvtLNuGamma::decay(), EvtSLN::decay(), EvtTauHadnu::decay(), EvtTaulnunu::decay(), EvtTauScalarnu::decay(), and EvtTauVectornu::decay().

◆ spParentNeutrino() [2/2]

virtual EvtDiracSpinor EvtParticle::spParentNeutrino ( ) const

virtual

Returns Dirac spinor in the parents restframe for a Neutrino particle.

Reimplemented in EvtNeutrinoParticle, and EvtNeutrinoParticle.

◆ treeStr() [1/2]

std::string EvtParticle::treeStr ( ) const

Definition at line 990 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                     {
 
  std::string retval=EvtPDL::name(_id);
  retval+=" -> ";
 
  retval+=treeStrRec(0);
 
  return retval;
}

◆ treeStr() [2/2]

std::string EvtParticle::treeStr ( ) const

◆ treeStrRec() [1/2]

std::string EvtParticle::treeStrRec ( int level ) const

Definition at line 908 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                                 {
 
  int newlevel,i;
  newlevel = level +1;
 
  std::string retval="";
 
  for(i=0;i<_ndaug;i++){
    retval+=EvtPDL::name(_daug[i]->getId());
    if ( _daug[i]->getNDaug() > 0 ) {
      retval+= " (";
      retval+= _daug[i]->treeStrRec(newlevel);
      retval+= ") ";
    }
    else{
      if ( i!=_ndaug-1) retval+=" ";
    }
  }
 
  return retval;
}

Referenced by treeStr(), and treeStrRec().

◆ treeStrRec() [2/2]

std::string EvtParticle::treeStrRec ( int level ) const

◆ writeTreeRec() [1/2]

std::string EvtParticle::writeTreeRec ( std::string resonance ) const

Definition at line 930 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc.

                                                             {  //pingrg, Jun. 6, 2008
  std::string retval="";
  
  if (resonance ==  EvtPDL::name(_id).c_str() && _ndaug!=0) {
    retval=resonance+": "+IntToStr(_ndaug)+"  ";
    for(int i=0;i<_ndaug;i++){
      retval += EvtPDL::name(_daug[i]->getId()).c_str();
      retval += " ";
    }
  }
 
    for(int i=0;i<_ndaug;i++){
      _daug[i]->writeTreeRec(resonance);
    }
    std::cout<<retval;
  return retval;
}

Referenced by writeTreeRec().

◆ writeTreeRec() [2/2]

std::string EvtParticle::writeTreeRec ( std::string ) const

Member Data Documentation

◆ _validP4

bool EvtParticle::_validP4

protected

Definition at line 407 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.hh.

Referenced by compMassProb(), EvtParticle(), hasValidP4(), EvtHighSpinParticle::init(), EvtRaritaSchwingerParticle::init(), EvtStringParticle::init(), EvtDiracParticle::init(), EvtNeutrinoParticle::init(), EvtPhotonParticle::init(), EvtScalarParticle::init(), EvtTensorParticle::init(), EvtVectorParticle::init(), and initDecay().

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

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ EvtParticle() [1/2]

◆ ~EvtParticle() [1/2]

◆ EvtParticle() [2/2]

◆ ~EvtParticle() [2/2]

Member Function Documentation

◆ addDaug() [1/2]

◆ addDaug() [2/2]

◆ compMassProb() [1/2]

◆ compMassProb() [2/2]

◆ decay() [1/2]

◆ decay() [2/2]

◆ decayProb() [1/2]

◆ decayProb() [2/2]

◆ deleteDaughters() [1/2]

◆ deleteDaughters() [2/2]

◆ deleteTree() [1/2]

◆ deleteTree() [2/2]

◆ dumpTree() [1/2]

◆ dumpTree() [2/2]

◆ dumpTreeRec() [1/2]

◆ dumpTreeRec() [2/2]

◆ eps() [1/2]

◆ eps() [2/2]

◆ epsParent() [1/2]

◆ epsParent() [2/2]

◆ epsParentPhoton() [1/2]

◆ epsParentPhoton() [2/2]

◆ epsPhoton() [1/2]

◆ epsPhoton() [2/2]

◆ epsTensor() [1/2]

◆ epsTensor() [2/2]

◆ epsTensorParent() [1/2]

◆ epsTensorParent() [2/2]

◆ firstornot() [1/2]

◆ firstornot() [2/2]

◆ generateMassTree() [1/2]

◆ generateMassTree() [2/2]

◆ get4Pos() [1/2]

◆ get4Pos() [2/2]

◆ getChannel() [1/2]

◆ getChannel() [2/2]

◆ getDaug() [1/2]

◆ getDaug() [2/2]

◆ getGeneratorFlag() [1/2]

◆ getGeneratorFlag() [2/2]

◆ getId() [1/2]

◆ getId() [2/2]

◆ getInclusiveMode() [1/2]

◆ getInclusiveMode() [2/2]

◆ getIntFlag() [1/2]

◆ getIntFlag() [2/2]

◆ getLifetime() [1/2]

◆ getLifetime() [2/2]

◆ getNDaug() [1/2]

◆ getNDaug() [2/2]

◆ getP4() [1/2]

◆ getP4() [2/2]

◆ getP4Lab() [1/2]

◆ getP4Lab() [2/2]

◆ getP4Restframe() [1/2]

◆ getP4Restframe() [2/2]

◆ getParent() [1/2]

◆ getParent() [2/2]

◆ getSpinDensityBackward() [1/2]

◆ getSpinDensityBackward() [2/2]

◆ getSpinDensityForward() [1/2]

◆ getSpinDensityForward() [2/2]

◆ getSpinStates() [1/2]

◆ getSpinStates() [2/2]

◆ getSpinType() [1/2]

◆ getSpinType() [2/2]

◆ hasValidP4() [1/2]

◆ hasValidP4() [2/2]

◆ init() [1/2]

◆ init() [2/2]