EvtFromParent Class Reference

#include <EvtFromParent.hh>

Inheritance diagram for EvtFromParent:

Public Member Functions
	EvtFromParent ()
virtual	~EvtFromParent ()
void	getName (std::string &name)
EvtDecayBase *	clone ()
void	initProbMax ()
void	init ()
void	decay (EvtParticle *p)
Public Member Functions inherited from EvtDecayIncoherent
void	makeDecay (EvtParticle *p)
virtual	~EvtDecayIncoherent ()
void	setDaughterSpinDensity (int daughter)
int	isDaughterSpinDensitySet (int daughter)
Public Member Functions inherited from EvtDecayBase
virtual std::string	commandName ()
virtual void	command (std::string cmd)
double	getProbMax (double prob)
double	resetProbMax (double prob)
	EvtDecayBase ()
virtual	~EvtDecayBase ()
virtual bool	matchingDecay (const EvtDecayBase &other) const
EvtId	getParentId ()
double	getBranchingFraction ()
void	disableCheckQ ()
void	checkQ ()
int	getNDaug ()
EvtId *	getDaugs ()
EvtId	getDaug (int i)
int	getNArg ()
int	getPHOTOS ()
void	setPHOTOS ()
void	setVerbose ()
void	setSummary ()
double *	getArgs ()
std::string *	getArgsStr ()
double	getArg (int j)
std::string	getArgStr (int j)
std::string	getModelName ()
int	getDSum ()
int	summary ()
int	verbose ()
void	saveDecayInfo (EvtId ipar, int ndaug, EvtId *daug, int narg, std::vector< std::string > &args, std::string name, double brfr)
void	printSummary ()
void	setProbMax (double prbmx)
void	noProbMax ()
void	checkNArg (int a1, int a2=-1, int a3=-1, int a4=-1)
void	checkNDaug (int d1, int d2=-1)
void	checkSpinParent (EvtSpinType::spintype sp)
void	checkSpinDaughter (int d1, EvtSpinType::spintype sp)
virtual int	nRealDaughters ()

Additional Inherited Members
Static Public Member Functions inherited from EvtDecayBase
static void	findMasses (EvtParticle *p, int ndaugs, EvtId daugs[10], double masses[10])
static void	findMass (EvtParticle *p)
static double	findMaxMass (EvtParticle *p)
Protected Member Functions inherited from EvtDecayBase
bool	daugsDecayedByParentModel ()
Protected Attributes inherited from EvtDecayBase
bool	_daugsDecayedByParentModel

Detailed Description

Definition at line 30 of file EvtFromParent.hh.

Constructor & Destructor Documentation

EvtFromParent()

EvtFromParent::EvtFromParent ( )

inline

Definition at line 33 of file EvtFromParent.hh.

{}

Referenced by clone().

~EvtFromParent()

EvtFromParent::~EvtFromParent ( )

virtual

Definition at line 39 of file EvtFromParent.cc.

{}

Member Function Documentation

clone()

EvtDecayBase * EvtFromParent::clone ( )

virtual

Implements EvtDecayBase.

Definition at line 43 of file EvtFromParent.cc.

{ return new EvtFromParent; }

decay()

void EvtFromParent::decay ( EvtParticle * p )

virtual

Implements EvtDecayBase.

Definition at line 53 of file EvtFromParent.cc.

                                          {
  bool debug = 0;
  // p->printTree();
  // std::cout<<"FromParent:p->getId() = "<<p->getId()<<std::endl;
  EvtId LambdaId = EvtPDL::getId( "Lambda0" );
  // std::cout<<"LambdaId = "<<LambdaId<<std::endl;
  EvtId antiLambdaId = EvtPDL::getId( "anti-Lambda0" );
  // std::cout<<"antiLambdaId = "<<antiLambdaId<<std::endl;
  if ( p->getId() != LambdaId && p->getId() != antiLambdaId )
  {
    std::cout << "Parent particle is required to be Lambda0 or antiLambda0" << std::endl;
    abort();
  }
  int         more      = 0;
  int         numstable = 0;
  EvtId       evtnumstable[100];
  EvtVector4R p4[20];
  /*if(EvtParticle::_NextLevelDauNum == 0) {
    report(ERROR,"EvtGen") << "EvtParticle failed to store the info of NextLevel";
  ::abort();
  }*/
  if ( EvtParticle::_NextLevelDauNum != 0 )
  {
    if ( p->getId() == LambdaId )
    {
      // pi-
      evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( -211 );
      p4[numstable]           = EvtParticle::_NextLevelP4[2];
      if ( debug )
        std::cout << "FromParent: pi- p4[numstable] = " << p4[numstable] << std::endl;
      numstable++;
      // p
      evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 2212 );
      p4[numstable]           = EvtParticle::_NextLevelP4[3];
      if ( debug )
        std::cout << "FromParent: p+ p4[numstable] = " << p4[numstable] << std::endl;
      numstable++;
    }
    if ( p->getId() == antiLambdaId )
    {
      if ( debug )
        std::cout << "anti-Lambda0's parent id:" << p->getParent()->getId()
                  << " p4 =" << p->getParent()->getP4() << std::endl;
      // pi+
      evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 211 );
      p4[numstable]           = EvtParticle::_NextLevelP4[0];
      if ( debug )
        std::cout << "FromParent: pi+ p4[numstable] = " << p4[numstable] << std::endl;
      numstable++;
      // pbar
      evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( -2212 );
      p4[numstable]           = EvtParticle::_NextLevelP4[1];
      if ( debug )
        std::cout << "FromParent: anti-p- p4[numstable] = " << p4[numstable] << std::endl;
      numstable++;
    }
    p->makeDaughters( numstable, evtnumstable );
    int ndaugFound = 0;
    for ( int i = 0; i < numstable; i++ )
    {
      p->getDaug( i )->init( evtnumstable[i], p4[i] );
      ndaugFound++;
    }
    if ( ndaugFound == 0 )
    {
      report( ERROR, "EvtGen" ) << "Phokhara has failed to do a decay ";
      report( ERROR, "EvtGen" ) << EvtPDL::name( p->getId() ).c_str() << " " << p->mass()
                                << endl;
      assert( 0 );
    }
    if ( p->getId() == LambdaId )
    {
      // int channel=EvtDecayTable::inChannelList(p->getId(),numstable,evtnumstable);
      // more=(channel!=-1)&&(channel!=0);
      // if(more) { std::cout<<"EvtFromParent:Existence of mode "<<channel<<" in exclusive
      // decay list has the same final state as this one"<<std::endl;abort(); }
 
      if ( debug )
        std::cout << "EvtFromParent Lambda0's parent id:" << p->getParent()->getId()
                  << " p4 =" << p->getParent()->getP4Lab() << std::endl;
      if ( debug ) std::cout << "EvtFromParent SUMMARY: part p4" << p->getP4Lab() << std::endl;
      if ( debug )
        std::cout << "EvtFromParent SUMMARY: Daug0 p4" << p->getDaug( 0 )->getP4Lab()
                  << std::endl;
      if ( debug )
        std::cout << "EvtFromParent SUMMARY: Daug1 p4" << p->getDaug( 1 )->getP4Lab()
                  << std::endl;
    }
    if ( p->getId() == antiLambdaId )
    {
      // int channel=EvtDecayTable::inChannelList(p->getId(),numstable,evtnumstable);
      // more=(channel!=0);
      // if(more) { std::cout<<"EvtFromParent:Existence of mode "<<channel<<" in exclusive
      // decay list has the same final state as this one"<<std::endl;abort(); }
 
      if ( debug )
        std::cout << "EvtFromParent antiLambda0's parent id:" << p->getParent()->getId()
                  << " p4 =" << p->getParent()->getP4Lab() << std::endl;
      if ( debug ) std::cout << "EvtFromParent SUMMARY: part p4" << p->getP4Lab() << std::endl;
      if ( debug )
        std::cout << "EvtFromParent SUMMARY: Daug0 p4" << p->getDaug( 0 )->getP4Lab()
                  << std::endl;
      if ( debug )
        std::cout << "EvtFromParent SUMMARY: Daug1 p4" << p->getDaug( 1 )->getP4Lab()
                  << std::endl;
    }
  }
  return;
}

getName()

void EvtFromParent::getName ( std::string & name )

virtual

Implements EvtDecayBase.

Definition at line 41 of file EvtFromParent.cc.

{ model_name = "FromParent"; }

init()

void EvtFromParent::init ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 45 of file EvtFromParent.cc.

                         {
 
  // check that there are 0 arguments
  checkNArg( 0 );
}

initProbMax()

void EvtFromParent::initProbMax ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 51 of file EvtFromParent.cc.

{ noProbMax(); }

---

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

• EvtFromParent.hh
• EvtFromParent.cc