#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	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)

Protected Attributes
bool	_validP4

Detailed Description

Definition at line 42 of file EvtParticle.hh.

Constructor & Destructor Documentation

◆ EvtParticle()

EvtParticle::EvtParticle ( )

Default constructor.

Definition at line 60 of file 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()

EvtParticle::~EvtParticle ( )

virtual

Destructor.

Definition at line 56 of file EvtParticle.cc.

                          {
  delete _decayProb;
}

Member Function Documentation

◆ addDaug()

void EvtParticle::addDaug ( EvtParticle * node )

Add another daughter to the particle

Definition at line 104 of file EvtParticle.cc.

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

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

◆ compMassProb()

double EvtParticle::compMassProb ( )

Definition at line 502 of file 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().

◆ decay()

void EvtParticle::decay ( )

Decay particle

Definition at line 404 of file 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);
 
   }
    else{
      p->_rhoBackward.SetDiag(p->getSpinStates());
 
    }
  }  
  _isDecayed=true;
  return;  
}

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

◆ decayProb()

double * EvtParticle::decayProb ( )

inline

Definition at line 376 of file EvtParticle.hh.

376{return _decayProb;}

◆ deleteDaughters()

void EvtParticle::deleteDaughters ( bool keepChannel = false )

Definition at line 538 of file 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().

◆ deleteTree()

void EvtParticle::deleteTree ( )

Delete a decay chain

Definition at line 555 of file EvtParticle.cc.

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

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

◆ dumpTree()

void EvtParticle::dumpTree ( ) const

Definition at line 976 of file 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;
 
}

◆ dumpTreeRec()

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

Definition at line 946 of file 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().

◆ eps()

EvtVector4C EvtParticle::eps ( int i ) const

virtual

Returns polarization vector in the particles own restframe.

Reimplemented in EvtVectorParticle.

Definition at line 574 of file 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 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().

◆ epsParent()

EvtVector4C EvtParticle::epsParent ( int i ) const

virtual

Returns polarization vector in the parents restframe.

Reimplemented in EvtVectorParticle.

Definition at line 563 of file 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().

◆ epsParentPhoton()

EvtVector4C EvtParticle::epsParentPhoton ( int i )

virtual

Returns polarization vector in the parents restframe for a photon.

Reimplemented in EvtPhotonParticle.

Definition at line 585 of file 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().

◆ epsPhoton()

EvtVector4C EvtParticle::epsPhoton ( int i )

virtual

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

Reimplemented in EvtPhotonParticle.

Definition at line 596 of file 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;
}

◆ epsTensor()

EvtTensor4C EvtParticle::epsTensor ( int i ) const

virtual

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

Reimplemented in EvtTensorParticle.

Definition at line 668 of file 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().

◆ epsTensorParent()

EvtTensor4C EvtParticle::epsTensorParent ( int i ) const

virtual

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

Reimplemented in EvtTensorParticle.

Definition at line 655 of file 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 EvtSSDCP::decay(), EvtSTS::decay(), and EvtSTSCP::decay().

◆ firstornot()

int EvtParticle::firstornot ( ) const

Used internally to decide if first time particle is decayed.

Definition at line 111 of file EvtParticle.cc.

111{ return _first;}

Referenced by EvtDecayBase::findMasses().

◆ generateMassTree()

void EvtParticle::generateMassTree ( )

Definition at line 461 of file 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() << " " << 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";
      assert(0);
    }
      }
    }
  }
 // report(INFO,"EvtGen") << counter << endl;
 // p->printTree();
}

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

◆ get4Pos()

EvtVector4R EvtParticle::get4Pos ( )

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

Definition at line 704 of file 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(), and makeStdHep().

◆ getChannel()

int EvtParticle::getChannel ( ) const

Retunrs the decay channel.

Definition at line 123 of file EvtParticle.cc.

123{ return _channel;}

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

◆ getDaug()

EvtParticle * EvtParticle::getDaug ( int i )

Get pointer the the i:th daugther.

Definition at line 85 of file EvtParticle.cc.

85{ return _daug[i]; }

Referenced by EvtbTosllScalarAmp::CalcAmp(), EvtbTosllVectorAmp::CalcAmp(), EvtSemiLeptonicBaryonAmp::CalcAmp(), EvtSemiLeptonicTensorAmp::CalcAmp(), EvtSemiLeptonicVectorAmp::CalcAmp(), EvtHypNonLepton::calcAmp(), EvtbTosllAmp::CalcMaxProb(), EvtSemiLeptonicAmp::CalcMaxProb(), EvtPsi3Sdecay::choseDecay(), compMassProb(), EvtSecondary::createSecondary(), EvtAngH2::decay(), EvtAngSam3::decay(), EvtAngSam::decay(), EvtAngSamLab::decay(), EvtAngSamX::decay(), EvtAV2GV::decay(), EvtBHadronic::decay(), EvtBody3::decay(), EvtBsquark::decay(), EvtBTo3piCP::decay(), EvtBTo4piCP::decay(), EvtBtoKD3P::decay(), EvtBToKpipiCP::decay(), EvtBtoXsEtap::decay(), EvtBtoXsgamma::decay(), EvtBtoXsll::decay(), EvtCBTo3piMPP::decay(), EvtCBTo3piP00::decay(), EvtChi0BB1::decay(), EvtChi0BB2::decay(), EvtChi1BB1::decay(), EvtChi1BB2::decay(), EvtChi2BB1::decay(), EvtChi2BB2::decay(), EvtConExc::decay(), EvtD0mixDalitz::decay(), EvtDDalitz::decay(), EvtDeBD::decay(), EvtDIY::decay(), EvtEtaDalitz::decay(), EvtFlatQ2::decay(), EvtHAngSam3::decay(), EvtHelPPJ::decay(), EvtHypWK::decay(), EvtIntervalDecayAmp< T >::decay(), EvtJ2BB1::decay(), EvtJ2BB2::decay(), EvtJ2BB3::decay(), EvtJetSet::decay(), EvtJPE::decay(), EvtJpipi::decay(), EvtJscont::decay(), EvtJTO3P::decay(), EvtKstarnunu::decay(), EvtKstarstargamma::decay(), EvtLambdaP_BarGamma::decay(), EvtLNuGamma::decay(), EvtLunda::decay(), EvtLundCharm::decay(), EvtMassH1::decay(), EvtMassH2::decay(), EvtMBody3::decay(), EvtmH2::decay(), EvtmPhsp::decay(), EvtMultibody::decay(), EvtOmegaDalitz::decay(), EvtOpenCharm::decay(), EvtP2GC0::decay(), EvtP2GC1::decay(), EvtP2GC2::decay(), decay(), EvtPhiDalitz::decay(), EvtPi0Dalitz::decay(), EvtPycont::decay(), EvtPyGaGa::decay(), EvtPythia::decay(), EvtRhoPi::decay(), EvtSingleParticle2::decay(), EvtSingleParticle::decay(), EvtSinglePoint::decay(), EvtSll::decay(), EvtSLN::decay(), EvtSPL::decay(), EvtSSDCP::decay(), EvtSTS::decay(), EvtSTSCP::decay(), EvtSVPCP::decay(), EvtSVPHelAmp::decay(), EvtSVS::decay(), EvtSVSCP::decay(), EvtSVSCPiso::decay(), EvtSVSCPLH::decay(), EvtSVSNONCPEIGEN::decay(), EvtT2GV::decay(), EvtTauHadnu::decay(), EvtTaulnunu::decay(), EvtTauola::decay(), EvtTauScalarnu::decay(), EvtTauVectornu::decay(), EvtTSS::decay(), EvtTVSPwave::decay(), EvtVectorIsr::decay(), EvtVll::decay(), EvtVPHOtoVISR::decay(), EvtVPHOtoVISRHi::decay(), EvtVSPPwave::decay(), EvtVSS::decay(), EvtVSSBMixCPT::decay(), EvtVSSMix::decay(), EvtVub::decay(), EvtVubHybrid::decay(), EvtVubNLO::decay(), EvtVVP::decay(), EvtVVpipi::decay(), EvtVVPIPI_WEIGHTED::decay(), EvtVVSPwave::decay(), EvtConExc::difgamXs(), EvtPHOTOS::doRadCorr(), EvtEvalHelAmp::evalAmp(), EvtDecayBase::findMass(), EvtDecayBase::findMasses(), EvtDecayBase::findMaxMass(), EvtConExc::findMaxXS(), EvtPsi3Sdecay::findMode(), EvtFSPick::FSPick(), EvtConExc::gamHXSection(), EvtGen::generateDecay(), initDecay(), initializePhaseSpace(), EvtDecayAmp::makeDecay(), EvtDecayIncoherent::makeDecay(), EvtDecayProb::makeDecay(), EvtDecayTag::makeTag(), EvtCPUtil::OtherB(), EvtCPUtil::OtherB(), EvtPsi3Sdecay::PHSPDecay(), EvtConExc::Rad1difXs(), EvtConExc::Rad2difXs(), and EvtSVVHelAmp::SVVHel().

◆ getGeneratorFlag()

int EvtParticle::getGeneratorFlag ( )

inline

get generator information; pingrg-2011-1-6

Definition at line 146 of file EvtParticle.hh.

146{ return _generatorFlag ;}

◆ getId()

EvtId EvtParticle::getId ( ) const

Returns Id of particle.

Definition at line 113 of file EvtParticle.cc.

113{ return _id;}

Referenced by EvtbTosllScalarAmp::CalcAmp(), EvtbTosllVectorAmp::CalcAmp(), EvtSemiLeptonicBaryonAmp::CalcAmp(), EvtSemiLeptonicTensorAmp::CalcAmp(), EvtSemiLeptonicVectorAmp::CalcAmp(), EvtPsi3Sdecay::choseDecay(), compMassProb(), EvtSecondary::createSecondary(), EvtAngH2::decay(), EvtBsquark::decay(), EvtBTo3piCP::decay(), EvtBtoXsll::decay(), EvtCBTo3piMPP::decay(), EvtCBTo3piP00::decay(), EvtConExc::decay(), EvtDDalitz::decay(), EvtDeBD::decay(), EvtDMix::decay(), EvtHypWK::decay(), EvtJetSet::decay(), EvtJscont::decay(), EvtKstarstargamma::decay(), EvtLNuGamma::decay(), EvtLunda::decay(), EvtLundCharm::decay(), EvtOpenCharm::decay(), decay(), EvtPycont::decay(), EvtPyGaGa::decay(), EvtPythia::decay(), EvtSVSCPiso::decay(), EvtSVSNONCPEIGEN::decay(), EvtSVVNONCPEIGEN::decay(), EvtTauHadnu::decay(), EvtTaulnunu::decay(), EvtTauola::decay(), EvtTauScalarnu::decay(), EvtTauVectornu::decay(), EvtVPHOtoVISR::decay(), EvtVPHOtoVISRHi::decay(), EvtVSSBMixCPT::decay(), EvtVSSMix::decay(), EvtPHOTOS::doRadCorr(), EvtDecayBase::findMass(), EvtDecayBase::findMasses(), EvtDecayBase::findMaxMass(), EvtConExc::findMaxXS(), EvtPsi3Sdecay::findMode(), EvtFSPick::FSPick(), EvtConExc::gamHXSection(), EvtDecayTable::getDecayFunc(), EvtParticleDecayList::getDecayModel(), initDecay(), initializePhaseSpace(), EvtDecayAmp::makeDecay(), EvtDecayProb::makeDecay(), makeStdHep(), makeStdHep(), EvtDecayTag::makeTag(), EvtCPUtil::OtherB(), EvtCPUtil::OtherB(), EvtPsi3Sdecay::PHSPDecay(), printTreeRec(), EvtDiracParticle::rotateToHelicityBasis(), EvtDiracParticle::rotateToHelicityBasis(), EvtHighSpinParticle::rotateToHelicityBasis(), EvtHighSpinParticle::rotateToHelicityBasis(), EvtRaritaSchwingerParticle::rotateToHelicityBasis(), EvtRaritaSchwingerParticle::rotateToHelicityBasis(), setLifetime(), treeStrRec(), and writeTreeRec().

◆ getInclusiveMode()

int EvtParticle::getInclusiveMode ( )

inline

Definition at line 380 of file EvtParticle.hh.

380{return _inclusiveMode;}

◆ getLifetime()

double EvtParticle::getLifetime ( )

Returns the lifetime.

Definition at line 99 of file EvtParticle.cc.

                               {
 
  return _t;
}

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

◆ getNDaug()

int EvtParticle::getNDaug ( ) const

Returns number of daugthers.

Definition at line 125 of file EvtParticle.cc.

125{ return _ndaug;}

Referenced by EvtPsi3Sdecay::choseDecay(), compMassProb(), EvtSecondary::createSecondary(), EvtBtoKD3P::decay(), EvtDIY::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(), EvtPsi3Sdecay::PHSPDecay(), and treeStrRec().

◆ getP4()

const EvtVector4R & EvtParticle::getP4 ( ) const

Returns 4momentum in parents restframe.

Definition at line 121 of file EvtParticle.cc.

121{ return _p;}

Referenced by EvtbTosllScalarAmp::CalcAmp(), EvtbTosllVectorAmp::CalcAmp(), EvtSemiLeptonicBaryonAmp::CalcAmp(), EvtSemiLeptonicTensorAmp::CalcAmp(), EvtSemiLeptonicVectorAmp::CalcAmp(), EvtPsi3Sdecay::choseDecay(), EvtAngSam3::decay(), EvtAngSam::decay(), EvtAngSamX::decay(), EvtAV2GV::decay(), EvtBHadronic::decay(), EvtBsquark::decay(), EvtBTo4piCP::decay(), EvtBtoKD3P::decay(), EvtChi0BB1::decay(), EvtChi0BB2::decay(), EvtChi1BB1::decay(), EvtChi1BB2::decay(), EvtChi2BB1::decay(), EvtChi2BB2::decay(), EvtConExc::decay(), EvtD0mixDalitz::decay(), EvtDDalitz::decay(), EvtDeBD::decay(), EvtDIY::decay(), EvtEtaDalitz::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(), EvtPi0Dalitz::decay(), EvtRhoPi::decay(), EvtSPL::decay(), EvtSSDCP::decay(), EvtSTS::decay(), EvtSTSCP::decay(), EvtSVPCP::decay(), EvtSVPHelAmp::decay(), EvtSVS::decay(), EvtSVSCP::decay(), EvtSVSNONCPEIGEN::decay(), EvtT2GV::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(), get4Pos(), getP4Lab(), EvtConExc::init(), EvtDecayAmp::makeDecay(), EvtPsi3Sdecay::PHSPDecay(), EvtConExc::Rad1difXs(), EvtConExc::Rad2difXs(), EvtDiracParticle::rotateToHelicityBasis(), EvtDiracParticle::rotateToHelicityBasis(), EvtRaritaSchwingerParticle::rotateToHelicityBasis(), EvtRaritaSchwingerParticle::rotateToHelicityBasis(), and EvtSVVHelAmp::SVVHel().

◆ getP4Lab()

EvtVector4R EvtParticle::getP4Lab ( )

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

Definition at line 683 of file 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 EvtAngH2::decay(), EvtBody3::decay(), EvtDeBD::decay(), EvtDIY::decay(), EvtLambdaP_BarGamma::decay(), EvtMassH1::decay(), EvtMassH2::decay(), EvtMBody3::decay(), EvtmH2::decay(), EvtmPhsp::decay(), EvtMultibody::decay(), EvtRhoPi::decay(), EvtSingleParticle2::decay(), EvtSingleParticle::decay(), EvtSPL::decay(), makeStdHep(), and makeStdHep().

◆ getP4Restframe()

EvtVector4R EvtParticle::getP4Restframe ( )

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

Definition at line 698 of file EvtParticle.cc.

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

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

◆ getParent()

EvtParticle * EvtParticle::getParent ( )

Returns pointer to parent particle.

Definition at line 87 of file EvtParticle.cc.

87{ return _parent;}

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

◆ getSpinDensityBackward()

EvtSpinDensity EvtParticle::getSpinDensityBackward ( )

inline

Get backward spin density matrix.

Definition at line 343 of file EvtParticle.hh.

343{return _rhoBackward;}

Referenced by EvtDecayAmp::makeDecay().

◆ getSpinDensityForward()

EvtSpinDensity EvtParticle::getSpinDensityForward ( )

inline

Get forward spin density matrix.

Definition at line 333 of file EvtParticle.hh.

333{return _rhoForward;}

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

◆ getSpinStates()

int EvtParticle::getSpinStates ( ) const

Returns number of spin states of the particle.

Definition at line 118 of file EvtParticle.cc.

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

EvtPDL::getSpinType

static EvtSpinType::spintype getSpinType(EvtId i)

Definition EvtPDL.hh:61

EvtSpinType::getSpinStates

static int getSpinStates(spintype stype)

Definition EvtSpinType.hh:64

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

◆ getSpinType()

EvtSpinType::spintype EvtParticle::getSpinType ( ) const

Returns particle type.

Definition at line 115 of file EvtParticle.cc.

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

◆ hasValidP4()

bool EvtParticle::hasValidP4 ( )

inline

Definition at line 369 of file EvtParticle.hh.

369{return _validP4;}

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

◆ init()

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

pure virtual

Initialiaze particle with id and 4momentum.

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

Referenced by EvtbTosllAmp::CalcMaxProb(), EvtSemiLeptonicAmp::CalcMaxProb(), EvtBTo3piCP::decay(), EvtBtoKD3P::decay(), EvtBToKpipiCP::decay(), EvtBtoXsEtap::decay(), EvtBtoXsgamma::decay(), EvtBtoXsll::decay(), EvtCBTo3piMPP::decay(), EvtCBTo3piP00::decay(), EvtIntervalDecayAmp< T >::decay(), EvtJetSet::decay(), EvtJscont::decay(), EvtLunda::decay(), EvtLundCharm::decay(), EvtOpenCharm::decay(), EvtPycont::decay(), EvtPyGaGa::decay(), EvtPythia::decay(), EvtSingleParticle2::decay(), EvtSingleParticle::decay(), EvtSinglePoint::decay(), EvtTauola::decay(), EvtVectorIsr::decay(), EvtVPHOtoVISR::decay(), EvtVPHOtoVISRHi::decay(), EvtVSSBMixCPT::decay(), EvtVub::decay(), EvtVubHybrid::decay(), EvtVubNLO::decay(), and initializePhaseSpace().

◆ initDecay()

void EvtParticle::initDecay ( bool useMinMass = false )

Definition at line 237 of file 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().

◆ initializePhaseSpace()

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 1069 of file 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(), EvtAngSam3::decay(), EvtAngSam::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(), EvtChi0BB1::decay(), EvtChi0BB2::decay(), EvtChi1BB1::decay(), EvtChi1BB2::decay(), EvtChi2BB1::decay(), EvtChi2BB2::decay(), EvtConExc::decay(), EvtD0mixDalitz::decay(), EvtDDalitz::decay(), EvtDeBD::decay(), EvtDIY::decay(), EvtDMix::decay(), EvtEtaDalitz::decay(), EvtFlatQ2::decay(), EvtHAngSam3::decay(), EvtHelAmp::decay(), EvtHelPPJ::decay(), EvtHQET2::decay(), EvtHQET::decay(), EvtHypNonLepton::decay(), EvtHypWK::decay(), EvtISGW2::decay(), EvtISGW::decay(), EvtJ2BB1::decay(), EvtJ2BB2::decay(), EvtJ2BB3::decay(), EvtJPE::decay(), EvtJpipi::decay(), EvtJTO3P::decay(), EvtKKLambdaC::decay(), EvtKstarnunu::decay(), EvtKstarstargamma::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(), 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(), EvtTauHadnu::decay(), EvtTaulnunu::decay(), EvtTauScalarnu::decay(), EvtTauVectornu::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(), EvtConExc::findMaxXS(), EvtConExc::gamHXSection(), EvtPsi3Sdecay::PHSPDecay(), and EvtSVVHelAmp::SVVHel().

◆ insertDaugPtr()

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

inline

Makes partptr the idaug:th daugther.

Definition at line 206 of file EvtParticle.hh.

206 { _daug[idaug]=partptr;

207 partptr->_parent=this; }

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

◆ isDecayed()

bool EvtParticle::isDecayed ( )

inline

Definition at line 370 of file EvtParticle.hh.

370{return _isDecayed;}

◆ isInitialized()

bool EvtParticle::isInitialized ( )

inline

Definition at line 368 of file EvtParticle.hh.

368{return _isInit;}

Referenced by EvtDecayBase::findMaxMass().

◆ makeDaughters()

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 1175 of file 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

Referenced by EvtbTosllAmp::CalcMaxProb(), EvtSemiLeptonicAmp::CalcMaxProb(), EvtPsi3Sdecay::choseDecay(), EvtBTo3piCP::decay(), EvtBToKpipiCP::decay(), EvtBtoXsEtap::decay(), EvtBtoXsgamma::decay(), EvtBtoXsll::decay(), EvtCBTo3piMPP::decay(), EvtCBTo3piP00::decay(), EvtConExc::decay(), EvtJetSet::decay(), EvtJscont::decay(), EvtLunda::decay(), EvtLundCharm::decay(), EvtOpenCharm::decay(), EvtPycont::decay(), EvtPyGaGa::decay(), EvtPythia::decay(), EvtSingleParticle2::decay(), EvtSingleParticle::decay(), EvtSinglePoint::decay(), EvtTauola::decay(), EvtConExc::findMaxXS(), EvtConExc::gamHXSection(), initDecay(), initializePhaseSpace(), and EvtPsi3Sdecay::PHSPDecay().

◆ makeStdHep() [1/2]

void EvtParticle::makeStdHep ( EvtStdHep & stdhep )

Definition at line 793 of file 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/2]

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

Makes stdhep list

Definition at line 757 of file 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().

◆ mass()

double EvtParticle::mass ( ) const

Returns mass of particle.

Definition at line 127 of file EvtParticle.cc.

                               {
 
     return _p.mass();
}

Referenced by EvtbTosllScalarAmp::CalcAmp(), EvtbTosllVectorAmp::CalcAmp(), EvtSemiLeptonicBaryonAmp::CalcAmp(), EvtSemiLeptonicTensorAmp::CalcAmp(), EvtSemiLeptonicVectorAmp::CalcAmp(), EvtbTosllAmp::CalcMaxProb(), EvtSemiLeptonicAmp::CalcMaxProb(), EvtPsi3Sdecay::choseDecay(), compMassProb(), EvtBHadronic::decay(), EvtBsquark::decay(), EvtBtoXsEtap::decay(), EvtBtoXsgamma::decay(), EvtBtoXsll::decay(), EvtDDalitz::decay(), EvtEtaDalitz::decay(), EvtHypWK::decay(), EvtJetSet::decay(), EvtJpipi::decay(), EvtJscont::decay(), EvtKstarnunu::decay(), EvtLunda::decay(), EvtLundCharm::decay(), EvtOmegaDalitz::decay(), EvtOpenCharm::decay(), EvtPycont::decay(), EvtPyGaGa::decay(), EvtPythia::decay(), EvtSll::decay(), EvtSLN::decay(), EvtSSDCP::decay(), EvtSTS::decay(), EvtSTSCP::decay(), EvtSVS::decay(), EvtSVSCPLH::decay(), EvtSVSNONCPEIGEN::decay(), EvtTauola::decay(), EvtTauVectornu::decay(), EvtTVSPwave::decay(), EvtVectorIsr::decay(), EvtVll::decay(), EvtVPHOtoVISR::decay(), EvtVPHOtoVISRHi::decay(), EvtVSPPwave::decay(), EvtVub::decay(), EvtVubHybrid::decay(), EvtVubNLO::decay(), EvtVVpipi::decay(), EvtVVPIPI_WEIGHTED::decay(), EvtPHOTOS::doRadCorr(), EvtDecayBase::findMass(), EvtDecayBase::findMasses(), EvtDecayBase::findMaxMass(), generateMassTree(), get4Pos(), EvtParticleDecayList::getDecayModel(), getP4Restframe(), EvtDiracParticle::init(), EvtRaritaSchwingerParticle::init(), initDecay(), initializePhaseSpace(), EvtDecayAmp::makeDecay(), EvtDecayProb::makeDecay(), EvtPsi3Sdecay::PHSPDecay(), printTreeRec(), EvtDiracParticle::rotateToHelicityBasis(), EvtDiracParticle::rotateToHelicityBasis(), EvtRaritaSchwingerParticle::rotateToHelicityBasis(), and EvtRaritaSchwingerParticle::rotateToHelicityBasis().

◆ nextIter()

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

Iterates over the particles in a decay chain.

Definition at line 727 of file 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];
 
}

◆ noLifeTime()

void EvtParticle::noLifeTime ( )

inline

Definition at line 348 of file EvtParticle.hh.

348{ _genlifetime=0; }

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

◆ printParticle()

void EvtParticle::printParticle ( ) const

Prints information for the particle.

Definition at line 998 of file 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().

◆ printTree()

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 895 of file 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().

◆ printTreeRec()

void EvtParticle::printTreeRec ( int level ) const

Definition at line 868 of file 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().

◆ resetFirstOrNot()

void EvtParticle::resetFirstOrNot ( )

Definition at line 77 of file EvtParticle.cc.

                                  {
  _first=1;
}

Referenced by EvtGen::generateDecay().

◆ resetNDaug()

void EvtParticle::resetNDaug ( )

inline

Definition at line 255 of file EvtParticle.hh.

255{_ndaug=0; return;}

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

◆ rotateToHelicityBasis() [1/2]

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, and EvtVectorParticle.

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

◆ rotateToHelicityBasis() [2/2]

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

pure virtual

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

◆ setChannel()

void EvtParticle::setChannel ( int i )

Should only be used internally.

Definition at line 81 of file EvtParticle.cc.

                                    { 
  _channel=i;
}

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

◆ setDecayProb()

void EvtParticle::setDecayProb ( double p )

Definition at line 1202 of file EvtParticle.cc.

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

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

◆ setDiagonalSpinDensity()

void EvtParticle::setDiagonalSpinDensity ( )

Set diagonal spindensity matrix.

Definition at line 133 of file EvtParticle.cc.

                                        {
 
  _rhoForward.SetDiag(getSpinStates());
}

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

◆ setFirstOrNot()

void EvtParticle::setFirstOrNot ( )

Definition at line 74 of file EvtParticle.cc.

                                {
  _first=0;
}

Referenced by EvtDecayBase::findMasses().

◆ setGeneratorFlag()

void EvtParticle::setGeneratorFlag ( int flag )

inline

set generator information; pingrg-2011-1-6

Definition at line 141 of file EvtParticle.hh.

141{_generatorFlag = flag;}

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

◆ setId()

void EvtParticle::setId ( EvtId id )

inline

Definition at line 351 of file EvtParticle.hh.

351{ _id=id;}

Referenced by makeDaughters().

◆ setInclusiveMode()

void EvtParticle::setInclusiveMode ( int im )

inline

Definition at line 379 of file EvtParticle.hh.

379{_inclusiveMode=im ;}

◆ setLifetime() [1/2]

void EvtParticle::setLifetime ( )

Generate lifetime according to pure exponential.

Definition at line 93 of file EvtParticle.cc.

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

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

◆ setLifetime() [2/2]

void EvtParticle::setLifetime ( double tau )

Set lifetime of the particle in parents restframe.

Definition at line 89 of file EvtParticle.cc.

                                       {
  _t=tau;
}

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

◆ setMass()

void EvtParticle::setMass ( double m )

inline

Definition at line 358 of file EvtParticle.hh.

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

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

◆ setp() [1/2]

void EvtParticle::setp ( const EvtVector4R & p4 )

inlineprotected

Definition at line 385 of file EvtParticle.hh.

385{ _p =p4; }

◆ setp() [2/2]

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

inlineprotected

Definition at line 384 of file EvtParticle.hh.

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

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

◆ setP4()

void EvtParticle::setP4 ( const EvtVector4R & p4 )

inline

Sets the 4momentum in the parents restframe.

Definition at line 244 of file EvtParticle.hh.

244{_p=p4;}

Referenced by EvtPHOTOS::doRadCorr().

◆ setpart_num()

void EvtParticle::setpart_num ( EvtId particle_number )

inlineprotected

Definition at line 386 of file EvtParticle.hh.

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

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

◆ setPolarizedSpinDensity()

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

Definition at line 157 of file 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);
 
}

◆ setSpinDensityBackward()

void EvtParticle::setSpinDensityBackward ( const EvtSpinDensity & rho )

inline

Set backward spin density matrix.

Definition at line 338 of file EvtParticle.hh.

338{_rhoBackward=rho;}

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

◆ setSpinDensityForward()

void EvtParticle::setSpinDensityForward ( const EvtSpinDensity & rho )

inline

Set forward spin density matrix.

Definition at line 307 of file EvtParticle.hh.

307{_rhoForward=rho;}

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

◆ setSpinDensityForwardHelicityBasis() [1/2]

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 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/2]

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

Definition at line 206 of file 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;
 
}

◆ setVectorSpinDensity()

void EvtParticle::setVectorSpinDensity ( )

Set spindensity matrix for e+e- -> V

Definition at line 138 of file 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().

◆ sp()

EvtDiracSpinor EvtParticle::sp ( int i ) const

virtual

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

Reimplemented in EvtDiracParticle.

Definition at line 620 of file 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 EvtSemiLeptonicBaryonAmp::CalcAmp(), EvtHypNonLepton::calcAmp(), EvtTauHadnu::decay(), EvtTaulnunu::decay(), EvtTauScalarnu::decay(), and EvtTauVectornu::decay().

◆ spNeutrino()

EvtDiracSpinor EvtParticle::spNeutrino ( ) const

virtual

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

Reimplemented in EvtNeutrinoParticle.

Definition at line 644 of file 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;
}

◆ spParent()

EvtDiracSpinor EvtParticle::spParent ( int i ) const

virtual

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

Reimplemented in EvtDiracParticle.

Definition at line 607 of file 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 EvtbTosllScalarAmp::CalcAmp(), EvtbTosllVectorAmp::CalcAmp(), EvtSemiLeptonicBaryonAmp::CalcAmp(), EvtSemiLeptonicTensorAmp::CalcAmp(), EvtSemiLeptonicVectorAmp::CalcAmp(), EvtHypNonLepton::calcAmp(), EvtBsquark::decay(), EvtLNuGamma::decay(), EvtSll::decay(), EvtSLN::decay(), EvtTaulnunu::decay(), and EvtVll::decay().

◆ spParentNeutrino()

EvtDiracSpinor EvtParticle::spParentNeutrino ( ) const

virtual

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

Reimplemented in EvtNeutrinoParticle.

Definition at line 633 of file 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().

◆ treeStr()

std::string EvtParticle::treeStr ( ) const

Definition at line 988 of file EvtParticle.cc.

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

◆ treeStrRec()

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

Definition at line 906 of file 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().

◆ writeTreeRec()

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

Definition at line 928 of file 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().

Member Data Documentation

◆ _validP4

bool EvtParticle::_validP4

protected

Definition at line 393 of file EvtParticle.hh.

Referenced by compMassProb(), EvtParticle(), hasValidP4(), EvtDiracParticle::init(), EvtHighSpinParticle::init(), EvtNeutrinoParticle::init(), EvtPhotonParticle::init(), EvtRaritaSchwingerParticle::init(), EvtScalarParticle::init(), EvtScalarParticle::init(), EvtStringParticle::init(), EvtTensorParticle::init(), EvtVectorParticle::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()

◆ ~EvtParticle()

Member Function Documentation

◆ addDaug()

◆ compMassProb()

◆ decay()

◆ decayProb()

◆ deleteDaughters()

◆ deleteTree()

◆ dumpTree()

◆ dumpTreeRec()

◆ eps()

◆ epsParent()

◆ epsParentPhoton()

◆ epsPhoton()

◆ epsTensor()

◆ epsTensorParent()

◆ firstornot()

◆ generateMassTree()

◆ get4Pos()

◆ getChannel()

◆ getDaug()

◆ getGeneratorFlag()

◆ getId()

◆ getInclusiveMode()

◆ getLifetime()

◆ getNDaug()

◆ getP4()

◆ getP4Lab()

◆ getP4Restframe()

◆ getParent()

◆ getSpinDensityBackward()

◆ getSpinDensityForward()

◆ getSpinStates()

◆ getSpinType()

◆ hasValidP4()

◆ init()

◆ initDecay()

◆ initializePhaseSpace()

◆ insertDaugPtr()

◆ isDecayed()

◆ isInitialized()

◆ makeDaughters()

◆ makeStdHep() [1/2]

◆ makeStdHep() [2/2]

◆ mass()

◆ nextIter()

◆ noLifeTime()

◆ printParticle()

◆ printTree()

◆ printTreeRec()

◆ resetFirstOrNot()

◆ resetNDaug()

◆ rotateToHelicityBasis() [1/2]

◆ rotateToHelicityBasis() [2/2]

◆ setChannel()

◆ setDecayProb()

◆ setDiagonalSpinDensity()

◆ setFirstOrNot()

◆ setGeneratorFlag()

◆ setId()

◆ setInclusiveMode()

◆ setLifetime() [1/2]

◆ setLifetime() [2/2]

◆ setMass()

◆ setp() [1/2]

◆ setp() [2/2]

◆ setP4()

◆ setpart_num()

◆ setPolarizedSpinDensity()

◆ setSpinDensityBackward()

◆ setSpinDensityForward()

◆ setSpinDensityForwardHelicityBasis() [1/2]

◆ setSpinDensityForwardHelicityBasis() [2/2]

◆ setVectorSpinDensity()