EvtParticle Class Reference

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

Static Public Attributes
static int	_NextLevelDauNum =0
static EvtId	_NextLevelId [20]
static EvtVector4R	_NextLevelP4 [20]

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 59 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 55 of file EvtParticle.cc.

                          {
  delete _decayProb;
}

Member Function Documentation

addDaug()

void EvtParticle::addDaug

(

EvtParticle *

node

)

Add another daughter to the particle

Definition at line 103 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 503 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 403 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);
      //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().

decayProb()

double * EvtParticle::decayProb ( )

inline

Definition at line 395 of file EvtParticle.hh.

{return _decayProb;}

deleteDaughters()

void EvtParticle::deleteDaughters

(

bool

keepChannel = false

)

Definition at line 539 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 556 of file EvtParticle.cc.

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

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

dumpTree()

void EvtParticle::dumpTree

(

)

const

Definition at line 977 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 947 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 575 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 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().

epsParent()

EvtVector4C EvtParticle::epsParent ( int  i ) const

virtual

Returns polarization vector in the parents restframe.

Reimplemented in EvtVectorParticle.

Definition at line 564 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 586 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 597 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 669 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 656 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 EvtSemiLeptonicTensorAmp::CalcAmp(), EvtBHadronic::decay(), 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 110 of file EvtParticle.cc.

{ return _first;}

Referenced by EvtDecayBase::findMasses().

generateMassTree()

void EvtParticle::generateMassTree

(

)

Definition at line 460 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() << _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().

get4Pos()

EvtVector4R EvtParticle::get4Pos

(

)

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

Definition at line 705 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().

getChannel()

int EvtParticle::getChannel

(

)

const

Retunrs the decay channel.

Definition at line 122 of file EvtParticle.cc.

{ 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 84 of file EvtParticle.cc.

{ return _daug[i]; }

Referenced by EvtRexc::angularSampling(), EvtConExc::angularSampling(), EvtHypNonLepton::calcAmp(), EvtbTosllScalarAmp::CalcAmp(), EvtbTosllVectorAmp::CalcAmp(), EvtSemiLeptonicBaryonAmp::CalcAmp(), EvtSemiLeptonicTensorAmp::CalcAmp(), EvtSemiLeptonicVectorAmp::CalcAmp(), EvtbTosllAmp::CalcMaxProb(), EvtSemiLeptonicAmp::CalcMaxProb(), EvtConExc::checkdecay(), EvtPsi3Sdecay::choseDecay(), compMassProb(), EvtSecondary::createSecondary(), decay(), EvtAngH2::decay(), EvtAngSam::decay(), EvtAngSam3::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(), EvtCalHelAmp::decay(), EvtCBTo3piMPP::decay(), EvtCBTo3piP00::decay(), EvtChi0BB1::decay(), EvtChi0BB2::decay(), EvtChi1BB1::decay(), EvtChi1BB2::decay(), EvtChi2BB1::decay(), EvtChi2BB2::decay(), EvtConExc::decay(), EvtD0mixDalitz::decay(), EvtD0To2pip2pim::decay(), EvtD0Toa0enu::decay(), EvtD0ToKKpipi::decay(), EvtD0toKLpipiPlot::decay(), EvtD0ToKpiEtap::decay(), EvtD0toKpietaPlot::decay(), EvtD0toKpiomegaPlot::decay(), EvtD0ToKpipi0::decay(), EvtD0ToKpipi0pi0::decay(), EvtD0toKpipi0Plot::decay(), EvtD0ToKpipipi::decay(), EvtD0ToKSKK::decay(), EvtD0toKSpi0etaPlot::decay(), EvtD0toKSpi0omegaPlot::decay(), EvtD0ToKSpi0pi0::decay(), EvtD0ToKSpipi::decay(), EvtD0ToKSpipipi0::decay(), EvtD0ToKSpipipi0pi0::decay(), EvtD0toKSpipiPlot::decay(), EvtD0Topipienu::decay(), EvtD0TopipiEta::decay(), EvtD0topipietaPlot::decay(), EvtD0Topippim2pi0::decay(), EvtDDalitz::decay(), EvtDeBD::decay(), EvtDIY::decay(), EvtDToa0enu::decay(), EvtDtoKKpiPlot::decay(), EvtDToKmPipPipPi0::decay(), EvtDToKpienu::decay(), EvtDtoKpipiPlot::decay(), EvtDToKppipipi0::decay(), EvtDToKSKpi0::decay(), EvtDToKSKSpi::decay(), EvtDToKSpieta::decay(), EvtDtoKSpietaPlot::decay(), EvtDtoKSpiomegaPlot::decay(), EvtDToKSpipi0pi0::decay(), EvtDToKSpipipi::decay(), EvtDTopi0pi0enu::decay(), EvtDTopiEtaEta::decay(), EvtDtopietaetaPlot::decay(), EvtDTopipi0Eta::decay(), EvtDtopipi0etaPlot::decay(), EvtDTopipi0pi0::decay(), EvtDtopipi0pi0Plot::decay(), EvtDTopipienu::decay(), EvtDtopipipiPlot::decay(), EvtEDM::decay(), EvtEtaDalitz::decay(), EvtEtap2gpipi::decay(), EvtFDC::decay(), EvtFlatQ2::decay(), EvtFromParent::decay(), EvtFromParent_antiLambda::decay(), EvtFromParent_Lambda::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(), EvtLambdac2pKpi::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(), EvtPhiDalitz::decay(), EvtPhokhara::decay(), EvtPhokhara_4pi::decay(), EvtPhokhara_etagam::decay(), EvtPhokhara_etapgam::decay(), EvtPhokhara_K0K0::decay(), EvtPhokhara_KK::decay(), EvtPhokhara_Lambda::decay(), EvtPhokhara_LLB::decay(), EvtPhokhara_nnbar::decay(), EvtPhokhara_pi0gam::decay(), EvtPhokhara_pi0pi0pipi::decay(), EvtPhokhara_pipi::decay(), EvtPhokhara_pipieta::decay(), EvtPhokhara_pipipi0::decay(), EvtPhokhara_ppbar::decay(), EvtPhsp::decay(), EvtPhspStepf::decay(), EvtPi0Dalitz::decay(), EvtPycont::decay(), EvtPyGaGa::decay(), EvtPythia::decay(), EvtRexc::decay(), EvtRhoPi::decay(), EvtSingleParticle::decay(), EvtSingleParticle2::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(), EvtTauGamMu::decay(), EvtTauHadnu::decay(), EvtTaulnunu::decay(), EvtTauola::decay(), EvtTauScalarnu::decay(), EvtTauVectornu::decay(), EvtTrackGen::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(), EvtLunda::ExclusiveDecay(), 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(), EvtConExc::photonSampling(), EvtPsi3Sdecay::PHSPDecay(), EvtConExc::Rad1difXs(), EvtConExc::Rad2difXs(), EvtConExc::SetP4(), EvtConExc::SetP4Rvalue(), and EvtSVVHelAmp::SVVHel().

getGeneratorFlag()

int EvtParticle::getGeneratorFlag ( )

inline

get generator information; pingrg-2011-1-6

Definition at line 146 of file EvtParticle.hh.

{ return _generatorFlag ;}

getId()

EvtId EvtParticle::getId

(

)

const

Returns Id of particle.

Definition at line 112 of file EvtParticle.cc.

{ return _id;}

Referenced by EvtRexc::angularSampling(), EvtConExc::angularSampling(), EvtbTosllScalarAmp::CalcAmp(), EvtbTosllVectorAmp::CalcAmp(), EvtSemiLeptonicBaryonAmp::CalcAmp(), EvtSemiLeptonicTensorAmp::CalcAmp(), EvtSemiLeptonicVectorAmp::CalcAmp(), EvtConExc::checkdecay(), EvtPsi3Sdecay::choseDecay(), compMassProb(), EvtSecondary::createSecondary(), decay(), EvtAngH2::decay(), EvtBsquark::decay(), EvtBTo3piCP::decay(), EvtBtoXsll::decay(), EvtCalHelAmp::decay(), EvtCBTo3piMPP::decay(), EvtCBTo3piP00::decay(), EvtConExc::decay(), EvtD0Toa0enu::decay(), EvtD0ToKKpipi::decay(), EvtD0Topipienu::decay(), EvtDDalitz::decay(), EvtDeBD::decay(), EvtDMix::decay(), EvtDToa0enu::decay(), EvtDToKmPipPipPi0::decay(), EvtDToKpienu::decay(), EvtDToKSpipi0pi0::decay(), EvtDTopi0pi0enu::decay(), EvtDTopipienu::decay(), EvtFromParent::decay(), EvtFromParent_antiLambda::decay(), EvtFromParent_Lambda::decay(), EvtHypWK::decay(), EvtJetSet::decay(), EvtJscont::decay(), EvtKstarstargamma::decay(), EvtLambdac2pKpi::decay(), EvtLNuGamma::decay(), EvtLunda::decay(), EvtLundCharm::decay(), EvtOpenCharm::decay(), EvtPhokhara::decay(), EvtPhokhara_4pi::decay(), EvtPhokhara_etagam::decay(), EvtPhokhara_etapgam::decay(), EvtPhokhara_K0K0::decay(), EvtPhokhara_KK::decay(), EvtPhokhara_Lambda::decay(), EvtPhokhara_LLB::decay(), EvtPhokhara_nnbar::decay(), EvtPhokhara_pi0gam::decay(), EvtPhokhara_pi0pi0pipi::decay(), EvtPhokhara_pipi::decay(), EvtPhokhara_pipieta::decay(), EvtPhokhara_pipipi0::decay(), EvtPhokhara_ppbar::decay(), EvtPycont::decay(), EvtPyGaGa::decay(), EvtPythia::decay(), EvtRexc::decay(), EvtSVSCPiso::decay(), EvtSVSNONCPEIGEN::decay(), EvtSVVNONCPEIGEN::decay(), EvtTauHadnu::decay(), EvtTaulnunu::decay(), EvtTauola::decay(), EvtTauScalarnu::decay(), EvtTauVectornu::decay(), EvtTrackGen::decay(), EvtVPHOtoVISR::decay(), EvtVPHOtoVISRHi::decay(), EvtVSSBMixCPT::decay(), EvtVSSMix::decay(), EvtPHOTOS::doRadCorr(), EvtLunda::ExclusiveDecay(), EvtDecayBase::findMass(), EvtDecayBase::findMasses(), EvtDecayBase::findMaxMass(), EvtConExc::findMaxXS(), EvtPsi3Sdecay::findMode(), EvtFSPick::FSPick(), EvtConExc::gamHXSection(), generateMassTree(), EvtDecayTable::getDecayFunc(), EvtParticleDecayList::getDecayModel(), initDecay(), initializePhaseSpace(), EvtDecayAmp::makeDecay(), EvtDecayProb::makeDecay(), makeStdHep(), EvtDecayTag::makeTag(), EvtCPUtil::OtherB(), EvtConExc::photonSampling(), EvtPsi3Sdecay::PHSPDecay(), printTreeRec(), EvtDiracParticle::rotateToHelicityBasis(), EvtHighSpinParticle::rotateToHelicityBasis(), EvtRaritaSchwingerParticle::rotateToHelicityBasis(), EvtConExc::selectMode(), setLifetime(), EvtConExc::SetP4(), EvtConExc::SetP4Rvalue(), treeStrRec(), and writeTreeRec().

getInclusiveMode()

int EvtParticle::getInclusiveMode ( )

inline

Definition at line 399 of file EvtParticle.hh.

{return _inclusiveMode;}

getIntFlag()

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

inline

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

Definition at line 163 of file EvtParticle.hh.

                              {
    return _intFlag;
  }

Referenced by EvtDecayTag::getModeTag().

getLifetime()

double EvtParticle::getLifetime

(

)

Returns the lifetime.

Definition at line 98 of file EvtParticle.cc.

                               {
 
  return _t;
}

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

getNDaug()

int EvtParticle::getNDaug

(

)

const

Returns number of daugthers.

Definition at line 124 of file EvtParticle.cc.

{ return _ndaug;}

Referenced by EvtRexc::angularSampling(), EvtConExc::angularSampling(), EvtConExc::checkdecay(), EvtPsi3Sdecay::choseDecay(), compMassProb(), EvtSecondary::createSecondary(), EvtBtoKD3P::decay(), EvtDIY::decay(), EvtFDC::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().

getP4()

const EvtVector4R & EvtParticle::getP4

(

)

const

Returns 4momentum in parents restframe.

Definition at line 120 of file EvtParticle.cc.

{ 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(), EvtD0To2pip2pim::decay(), EvtD0Toa0enu::decay(), EvtD0ToKKpipi::decay(), EvtD0ToKpiEtap::decay(), EvtD0ToKpipi0::decay(), EvtD0ToKpipi0pi0::decay(), EvtD0ToKpipipi::decay(), EvtD0ToKSKK::decay(), EvtD0ToKSpi0pi0::decay(), EvtD0ToKSpipi::decay(), EvtD0ToKSpipipi0::decay(), EvtD0ToKSpipipi0pi0::decay(), EvtD0Topipienu::decay(), EvtD0TopipiEta::decay(), EvtD0Topippim2pi0::decay(), EvtDDalitz::decay(), EvtDeBD::decay(), EvtDIY::decay(), EvtDToa0enu::decay(), EvtDToKmPipPipPi0::decay(), EvtDToKpienu::decay(), EvtDToKppipipi0::decay(), EvtDToKSKpi0::decay(), EvtDToKSKSpi::decay(), EvtDToKSpieta::decay(), EvtDToKSpipi0pi0::decay(), EvtDToKSpipipi::decay(), EvtDTopi0pi0enu::decay(), EvtDTopiEtaEta::decay(), EvtDTopipi0Eta::decay(), EvtDTopipi0pi0::decay(), EvtDTopipienu::decay(), EvtEDM::decay(), EvtEtaDalitz::decay(), EvtEtap2gpipi::decay(), EvtFlatQ2::decay(), EvtFromParent::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(), EvtPhspStepf::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().

getP4Lab()

EvtVector4R EvtParticle::getP4Lab

(

)

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

Definition at line 684 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 EvtRexc::angularSampling(), EvtConExc::angularSampling(), EvtAngH2::decay(), EvtBody3::decay(), EvtD0To2pip2pim::decay(), EvtD0toKLpipiPlot::decay(), EvtD0toKpietaPlot::decay(), EvtD0toKpiomegaPlot::decay(), EvtD0toKpipi0Plot::decay(), EvtD0toKSpi0etaPlot::decay(), EvtD0toKSpi0omegaPlot::decay(), EvtD0ToKSpipi::decay(), EvtD0ToKSpipipi0::decay(), EvtD0toKSpipiPlot::decay(), EvtD0topipietaPlot::decay(), EvtD0Topippim2pi0::decay(), EvtDeBD::decay(), EvtDIY::decay(), EvtDtoKKpiPlot::decay(), EvtDtoKpipiPlot::decay(), EvtDToKSKpi0::decay(), EvtDtoKSpietaPlot::decay(), EvtDtoKSpiomegaPlot::decay(), EvtDtopietaetaPlot::decay(), EvtDtopipi0etaPlot::decay(), EvtDtopipi0pi0Plot::decay(), EvtDtopipipiPlot::decay(), EvtFDC::decay(), EvtFromParent::decay(), EvtFromParent_antiLambda::decay(), EvtFromParent_Lambda::decay(), EvtLambdaP_BarGamma::decay(), EvtMassH1::decay(), EvtMassH2::decay(), EvtMBody3::decay(), EvtmH2::decay(), EvtmPhsp::decay(), EvtMultibody::decay(), EvtPhokhara_Lambda::decay(), EvtRhoPi::decay(), EvtSingleParticle::decay(), EvtSingleParticle2::decay(), EvtSPL::decay(), makeStdHep(), and EvtConExc::photonSampling().

getP4Restframe()

EvtVector4R EvtParticle::getP4Restframe

(

)

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

Definition at line 699 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 86 of file EvtParticle.cc.

{ return _parent;}

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

getSpinDensityBackward()

EvtSpinDensity EvtParticle::getSpinDensityBackward ( )

inline

Get backward spin density matrix.

Definition at line 362 of file EvtParticle.hh.

{return _rhoBackward;}

Referenced by EvtDecayAmp::makeDecay().

getSpinDensityForward()

EvtSpinDensity EvtParticle::getSpinDensityForward ( )

inline

Get forward spin density matrix.

Definition at line 352 of file EvtParticle.hh.

{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 117 of file EvtParticle.cc.

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

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 114 of file EvtParticle.cc.

      { return EvtPDL::getSpinType(_id);}

hasValidP4()

bool EvtParticle::hasValidP4 ( )

inline

Definition at line 388 of file EvtParticle.hh.

{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 EvtHighSpinParticle, EvtRaritaSchwingerParticle, EvtStringParticle, EvtScalarParticle, EvtVectorParticle, EvtDiracParticle, EvtNeutrinoParticle, EvtPhotonParticle, and EvtTensorParticle.

Referenced by EvtbTosllAmp::CalcMaxProb(), EvtSemiLeptonicAmp::CalcMaxProb(), EvtBTo3piCP::decay(), EvtBtoKD3P::decay(), EvtBToKpipiCP::decay(), EvtBtoXsEtap::decay(), EvtBtoXsgamma::decay(), EvtBtoXsll::decay(), EvtCBTo3piMPP::decay(), EvtCBTo3piP00::decay(), EvtFromParent::decay(), EvtFromParent_antiLambda::decay(), EvtFromParent_Lambda::decay(), EvtIntervalDecayAmp< T >::decay(), EvtJetSet::decay(), EvtJscont::decay(), EvtLambdac2pKpi::decay(), EvtLunda::decay(), EvtLundCharm::decay(), EvtOpenCharm::decay(), EvtPhokhara::decay(), EvtPhokhara_4pi::decay(), EvtPhokhara_etagam::decay(), EvtPhokhara_etapgam::decay(), EvtPhokhara_K0K0::decay(), EvtPhokhara_KK::decay(), EvtPhokhara_Lambda::decay(), EvtPhokhara_LLB::decay(), EvtPhokhara_nnbar::decay(), EvtPhokhara_pi0gam::decay(), EvtPhokhara_pi0pi0pipi::decay(), EvtPhokhara_pipi::decay(), EvtPhokhara_pipieta::decay(), EvtPhokhara_pipipi0::decay(), EvtPhokhara_ppbar::decay(), EvtPycont::decay(), EvtPyGaGa::decay(), EvtPythia::decay(), EvtSingleParticle::decay(), EvtSingleParticle2::decay(), EvtSinglePoint::decay(), EvtTauola::decay(), EvtTrackGen::decay(), EvtVectorIsr::decay(), EvtVPHOtoVISR::decay(), EvtVPHOtoVISRHi::decay(), EvtVSSBMixCPT::decay(), EvtVub::decay(), EvtVubHybrid::decay(), EvtVubNLO::decay(), initializePhaseSpace(), EvtConExc::SetP4(), and EvtConExc::SetP4Rvalue().

initDecay()

void EvtParticle::initDecay

(

bool

useMinMass = false

)

Definition at line 236 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 1070 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(), 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(), EvtD0To2pip2pim::decay(), EvtD0Toa0enu::decay(), EvtD0ToKKpipi::decay(), EvtD0toKLpipiPlot::decay(), EvtD0ToKpiEtap::decay(), EvtD0toKpietaPlot::decay(), EvtD0toKpiomegaPlot::decay(), EvtD0ToKpipi0::decay(), EvtD0ToKpipi0pi0::decay(), EvtD0toKpipi0Plot::decay(), EvtD0ToKpipipi::decay(), EvtD0ToKSKK::decay(), EvtD0toKSpi0etaPlot::decay(), EvtD0toKSpi0omegaPlot::decay(), EvtD0ToKSpi0pi0::decay(), EvtD0ToKSpipi::decay(), EvtD0ToKSpipipi0::decay(), EvtD0ToKSpipipi0pi0::decay(), EvtD0toKSpipiPlot::decay(), EvtD0Topipienu::decay(), EvtD0TopipiEta::decay(), EvtD0topipietaPlot::decay(), EvtD0Topippim2pi0::decay(), EvtDDalitz::decay(), EvtDeBD::decay(), EvtDIY::decay(), EvtDMix::decay(), EvtDToa0enu::decay(), EvtDtoKKpiPlot::decay(), EvtDToKmPipPipPi0::decay(), EvtDToKpienu::decay(), EvtDtoKpipiPlot::decay(), EvtDToKppipipi0::decay(), EvtDToKSKpi0::decay(), EvtDToKSKSpi::decay(), EvtDToKSpieta::decay(), EvtDtoKSpietaPlot::decay(), EvtDtoKSpiomegaPlot::decay(), EvtDToKSpipi0pi0::decay(), EvtDToKSpipipi::decay(), EvtDTopi0pi0enu::decay(), EvtDTopiEtaEta::decay(), EvtDtopietaetaPlot::decay(), EvtDTopipi0Eta::decay(), EvtDtopipi0etaPlot::decay(), EvtDTopipi0pi0::decay(), EvtDtopipi0pi0Plot::decay(), EvtDTopipienu::decay(), EvtDtopipipiPlot::decay(), EvtEDM::decay(), EvtEtaDalitz::decay(), EvtEtap2gpipi::decay(), EvtFDC::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(), EvtPhspStepf::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().

insertDaugPtr()

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

inline

Makes partptr the idaug:th daugther.

Definition at line 225 of file EvtParticle.hh.

                                                    { _daug[idaug]=partptr;
                                                 partptr->_parent=this; }

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

isDecayed()

bool EvtParticle::isDecayed ( )

inline

Definition at line 389 of file EvtParticle.hh.

{return _isDecayed;}

isInitialized()

bool EvtParticle::isInitialized ( )

inline

Definition at line 387 of file EvtParticle.hh.

{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 1176 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(), EvtFromParent::decay(), EvtFromParent_antiLambda::decay(), EvtFromParent_Lambda::decay(), EvtJetSet::decay(), EvtJscont::decay(), EvtLunda::decay(), EvtLundCharm::decay(), EvtOpenCharm::decay(), EvtPhokhara::decay(), EvtPhokhara_4pi::decay(), EvtPhokhara_etagam::decay(), EvtPhokhara_etapgam::decay(), EvtPhokhara_K0K0::decay(), EvtPhokhara_KK::decay(), EvtPhokhara_Lambda::decay(), EvtPhokhara_LLB::decay(), EvtPhokhara_nnbar::decay(), EvtPhokhara_pi0gam::decay(), EvtPhokhara_pi0pi0pipi::decay(), EvtPhokhara_pipi::decay(), EvtPhokhara_pipieta::decay(), EvtPhokhara_pipipi0::decay(), EvtPhokhara_ppbar::decay(), EvtPycont::decay(), EvtPyGaGa::decay(), EvtPythia::decay(), EvtRexc::decay(), EvtSingleParticle::decay(), EvtSingleParticle2::decay(), EvtSinglePoint::decay(), EvtTauola::decay(), EvtLunda::ExclusiveDecay(), EvtConExc::findMaxXS(), EvtConExc::gamHXSection(), initDecay(), initializePhaseSpace(), and EvtPsi3Sdecay::PHSPDecay().

makeStdHep() [1/2]

void EvtParticle::makeStdHep

(

EvtStdHep &

stdhep

)

Definition at line 794 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 758 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 126 of file EvtParticle.cc.

                               {
 
     return _p.mass();
}

Referenced by EvtbTosllScalarAmp::CalcAmp(), EvtbTosllVectorAmp::CalcAmp(), EvtSemiLeptonicBaryonAmp::CalcAmp(), EvtSemiLeptonicTensorAmp::CalcAmp(), EvtSemiLeptonicVectorAmp::CalcAmp(), EvtbTosllAmp::CalcMaxProb(), EvtSemiLeptonicAmp::CalcMaxProb(), EvtConExc::checkdecay(), EvtPsi3Sdecay::choseDecay(), compMassProb(), EvtBHadronic::decay(), EvtBsquark::decay(), EvtBtoXsEtap::decay(), EvtBtoXsgamma::decay(), EvtBtoXsll::decay(), EvtConExc::decay(), EvtDDalitz::decay(), EvtEDM::decay(), EvtEtaDalitz::decay(), EvtFromParent::decay(), EvtFromParent_antiLambda::decay(), EvtFromParent_Lambda::decay(), EvtHypWK::decay(), EvtJetSet::decay(), EvtJpipi::decay(), EvtJscont::decay(), EvtKstarnunu::decay(), EvtLunda::decay(), EvtLundCharm::decay(), EvtOmegaDalitz::decay(), EvtOpenCharm::decay(), EvtPhokhara::decay(), EvtPhokhara_4pi::decay(), EvtPhokhara_etagam::decay(), EvtPhokhara_etapgam::decay(), EvtPhokhara_K0K0::decay(), EvtPhokhara_KK::decay(), EvtPhokhara_Lambda::decay(), EvtPhokhara_LLB::decay(), EvtPhokhara_nnbar::decay(), EvtPhokhara_pi0gam::decay(), EvtPhokhara_pi0pi0pipi::decay(), EvtPhokhara_pipi::decay(), EvtPhokhara_pipieta::decay(), EvtPhokhara_pipipi0::decay(), EvtPhokhara_ppbar::decay(), EvtPycont::decay(), EvtPyGaGa::decay(), EvtPythia::decay(), EvtRexc::decay(), EvtSll::decay(), EvtSLN::decay(), EvtSSDCP::decay(), EvtSTS::decay(), EvtSTSCP::decay(), EvtSVS::decay(), EvtSVSCPLH::decay(), EvtSVSNONCPEIGEN::decay(), EvtTauola::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(), EvtLunda::ExclusiveDecay(), EvtDecayBase::findMass(), EvtDecayBase::findMasses(), EvtDecayBase::findMaxMass(), generateMassTree(), get4Pos(), EvtParticleDecayList::getDecayModel(), getP4Restframe(), EvtRaritaSchwingerParticle::init(), EvtDiracParticle::init(), EvtPhokhara::init_evt(), EvtPhokhara::init_mode(), initDecay(), initializePhaseSpace(), EvtDecayAmp::makeDecay(), EvtDecayProb::makeDecay(), EvtPsi3Sdecay::PHSPDecay(), printTreeRec(), EvtDiracParticle::rotateToHelicityBasis(), and EvtRaritaSchwingerParticle::rotateToHelicityBasis().

nextIter()

EvtParticle * EvtParticle::nextIter

(

EvtParticle *

rootOfTree = 0

)

Iterates over the particles in a decay chain.

Definition at line 728 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 367 of file EvtParticle.hh.

{ _genlifetime=0; }

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

printParticle()

void EvtParticle::printParticle

(

)

const

Prints information for the particle.

Definition at line 999 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 896 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 869 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 76 of file EvtParticle.cc.

                                  {
  _first=1;
}

Referenced by EvtGen::generateDecay().

resetNDaug()

void EvtParticle::resetNDaug ( )

inline

Definition at line 274 of file EvtParticle.hh.

{_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(), 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 80 of file EvtParticle.cc.

                                    { 
  _channel=i;
}

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

setDecayProb()

void EvtParticle::setDecayProb

(

double

p

)

Definition at line 1203 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 132 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 73 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.

{_generatorFlag = flag;}

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

setId()

void EvtParticle::setId ( EvtId  id )

inline

Definition at line 370 of file EvtParticle.hh.

{ _id=id;}

Referenced by makeDaughters().

setInclusiveMode()

void EvtParticle::setInclusiveMode ( int  im )

inline

Definition at line 398 of file EvtParticle.hh.

{_inclusiveMode=im ;}

setIntFlag()

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

inline

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

Definition at line 156 of file EvtParticle.hh.

                                    {
    _intFlag=vi;
  }

Referenced by EvtConExc::decay().

setLifetime() [1/2]

void EvtParticle::setLifetime

(

)

Generate lifetime according to pure exponential.

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

void EvtParticle::setLifetime

(

double

tau

)

Set lifetime of the particle in parents restframe.

Definition at line 88 of file EvtParticle.cc.

                                       {
  _t=tau;
}

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

setMass()

void EvtParticle::setMass ( double  m )

inline

Definition at line 377 of file EvtParticle.hh.

{ _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().

setp() [1/2]

void EvtParticle::setp ( const EvtVector4R &  p4 )

inlineprotected

Definition at line 404 of file EvtParticle.hh.

{ _p =p4; }

setp() [2/2]

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

inlineprotected

Definition at line 403 of file EvtParticle.hh.

{ _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().

setP4()

void EvtParticle::setP4 ( const EvtVector4R &  p4 )

inline

Sets the 4momentum in the parents restframe.

Definition at line 263 of file EvtParticle.hh.

{_p=p4;}

Referenced by EvtPHOTOS::doRadCorr().

setpart_num()

void EvtParticle::setpart_num ( EvtId  particle_number )

inlineprotected

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

setPolarizedSpinDensity()

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

Definition at line 156 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 357 of file EvtParticle.hh.

{_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 326 of file EvtParticle.hh.

{_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 177 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 205 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 137 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 621 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 EvtHypNonLepton::calcAmp(), EvtSemiLeptonicBaryonAmp::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 645 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 608 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 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().

spParentNeutrino()

EvtDiracSpinor EvtParticle::spParentNeutrino ( ) const

virtual

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

Reimplemented in EvtNeutrinoParticle.

Definition at line 634 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 989 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 907 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 929 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

_NextLevelDauNum

int EvtParticle::_NextLevelDauNum =0

static

Definition at line 149 of file EvtParticle.hh.

Referenced by EvtFromParent::decay(), EvtFromParent_antiLambda::decay(), EvtFromParent_Lambda::decay(), EvtPhokhara_Lambda::decay(), and EvtDecay::execute().

_NextLevelId

EvtId EvtParticle::_NextLevelId

static

Definition at line 150 of file EvtParticle.hh.

Referenced by EvtPhokhara_Lambda::decay().

_NextLevelP4

EvtVector4R EvtParticle::_NextLevelP4

static

Definition at line 151 of file EvtParticle.hh.

Referenced by EvtFromParent::decay(), EvtFromParent_antiLambda::decay(), EvtFromParent_Lambda::decay(), and EvtPhokhara_Lambda::decay().

_validP4

bool EvtParticle::_validP4

protected

Definition at line 412 of file 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:

• EvtParticle.hh
• EvtDecay.cxx
• EvtParticle.cc