EvtPythia Class Reference

#include <EvtPythia.hh>

Inheritance diagram for EvtPythia:

Public Member Functions
	EvtPythia ()
virtual	~EvtPythia ()
void	getName (std::string &name)
EvtDecayBase *	clone ()
void	decay (EvtParticle *p)
std::string	commandName ()
void	command (std::string cmd)
void	init ()
void	initProbMax ()
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
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 ()

Static Public Member Functions
static void	pythiaInit (int f)
static void	pythiacont (double *, int *, int *, double *, double *, double *, double *)
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)

Additional Inherited Members
Protected Member Functions inherited from EvtDecayBase
bool	daugsDecayedByParentModel ()
Protected Attributes inherited from EvtDecayBase
bool	_daugsDecayedByParentModel

Detailed Description

Definition at line 32 of file EvtPythia.hh.

Constructor & Destructor Documentation

EvtPythia()

EvtPythia::EvtPythia

(

)

Definition at line 91 of file EvtPythia.cc.

{}

Referenced by clone().

~EvtPythia()

EvtPythia::~EvtPythia ( )

virtual

Definition at line 93 of file EvtPythia.cc.

                      {
 
  int i;
 
  // the deletion of commands is really uggly!
 
  if ( njetsetdecays == 0 )
  {
    delete[] commands;
    commands = 0;
    return;
  }
 
  for ( i = 0; i < njetsetdecays; i++ )
  {
    if ( jetsetdecays[i] == this )
    {
      jetsetdecays[i] = jetsetdecays[njetsetdecays - 1];
      njetsetdecays--;
      if ( njetsetdecays == 0 )
      {
        delete[] commands;
        commands = 0;
      }
      return;
    }
  }
 
  report( ERROR, "EvtGen" ) << "Error in destroying Pythia model!" << endl;
}

Member Function Documentation

clone()

EvtDecayBase * EvtPythia::clone ( )

virtual

Implements EvtDecayBase.

Definition at line 126 of file EvtPythia.cc.

{ return new EvtPythia; }

command()

void EvtPythia::command ( std::string cmd )

virtual

Reimplemented from EvtDecayBase.

Definition at line 150 of file EvtPythia.cc.

                                       {
 
  if ( ncommand == lcommand )
  {
 
    lcommand = 10 + 2 * lcommand;
 
    std::string* newcommands = new std::string[lcommand];
 
    int i;
 
    for ( i = 0; i < ncommand; i++ ) { newcommands[i] = commands[i]; }
 
    delete[] commands;
 
    commands = newcommands;
  }
 
  commands[ncommand] = cmd;
 
  ncommand++;
}

commandName()

std::string EvtPythia::commandName ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 148 of file EvtPythia.cc.

{ return std::string( "JetSetPar" ); }

decay()

void EvtPythia::decay ( EvtParticle * p )

virtual

Implements EvtDecayBase.

Definition at line 178 of file EvtPythia.cc.

                                      {
 
  // added by Lange Jan4,2000
  static EvtId STRNG = EvtPDL::getId( "string" );
 
  int istdheppar = EvtPDL::getStdHep( p->getId() );
 
  if ( pycomp_( &istdheppar ) == 0 )
  {
    report( ERROR, "EvtGen" ) << "Pythia can not decay:" << EvtPDL::name( p->getId() ).c_str()
                              << endl;
    return;
  }
 
  double mp = p->mass();
 
  EvtVector4R p4[20];
 
  int i, more;
  int ip = EvtPDL::getStdHep( p->getId() );
 
  int   ndaugjs;
  int   kf[100];
  EvtId evtnumstable[100], evtnumparton[100];
  int   stableindex[100], partonindex[100];
  int   numstable;
  int   numparton;
  int   km[100];
  EvtId type[MAX_DAUG];
 
  cbbeam_.maximum = mp; // pingrg
  if ( mp == 0 )
  {
    std::cout << "Particle " << EvtPDL::name( p->getId() ) << " has zero mass" << std::endl;
    abort();
  }
  pythiaInit( 0 );
 
  double px[100], py[100], pz[100], e[100];
  if ( p->getNDaug() != 0 ) { p->deleteDaughters( true ); }
 
  int count = 0;
 
  do {
 
    pythiadec_( &ip, &mp, &ndaugjs, kf, km, px, py, pz, e );
 
    numstable = 0;
    numparton = 0;
 
    for ( i = 0; i < ndaugjs; i++ )
    {
      // std::cout<<"EvtPDL::evtIdFromStdHep(kf[i]),i= "<<i<<"
      // "<<EvtPDL::evtIdFromStdHep(kf[i])<<std::endl;
      if ( EvtPDL::evtIdFromStdHep( kf[i] ) == EvtId( -1, -1 ) )
      {
        report( ERROR, "EvtGen" ) << "Pythia returned particle:" << kf[i] << endl;
        report( ERROR, "EvtGen" ) << "This can not be translated to evt number" << endl;
        report( ERROR, "EvtGen" ) << "and the decay will be rejected!" << endl;
        report( ERROR, "EvtGen" ) << "The decay was of particle:" << ip << endl;
        int i = 1;
        pylist_( i );
      }
 
      // sort out the partons
      if ( abs( kf[i] ) <= 6 || kf[i] == 21 )
      {
        partonindex[numparton]  = i;
        evtnumparton[numparton] = EvtPDL::evtIdFromStdHep( kf[i] );
        numparton++;
      }
      else
      {
        stableindex[numstable]  = i;
        evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( kf[i] );
        numstable++;
      }
 
      // have to protect against negative mass^2 for massless particles
      // i.e. neutrinos and photons.
      // this is uggly but I need to fix it right now....
 
      if ( px[i] * px[i] + py[i] * py[i] + pz[i] * pz[i] >= e[i] * e[i] )
      { e[i] = sqrt( px[i] * px[i] + py[i] * py[i] + pz[i] * pz[i] ) + 0.0000000000001; }
 
      p4[i].set( e[i], px[i], py[i], pz[i] );
    }
 
    int channel = EvtDecayTable::inChannelList( p->getId(), numstable, evtnumstable );
 
    more = ( channel != -1 );
 
    count++;
 
  } while ( more && ( count < 10000 ) );
 
  if ( count > 9999 )
  {
    report( INFO, "EvtGen" ) << "Too many loops in EvtPythia!!!" << endl;
    report( INFO, "EvtGen" ) << "Parent:" << EvtPDL::name( getParentId() ).c_str() << endl;
    for ( i = 0; i < numstable; i++ )
    {
      report( INFO, "EvtGen" ) << "Daug(" << i << ")"
                               << EvtPDL::name( evtnumstable[i] ).c_str() << endl;
    }
  }
 
  if ( numparton == 0 )
  {
 
    p->makeDaughters( numstable, evtnumstable );
 
    for ( i = 0; i < numstable; i++ )
    { p->getDaug( i )->init( evtnumstable[i], p4[stableindex[i]] ); }
 
    fixPolarizations( p );
 
    return;
  }
  else
  {
 
    // have partons in JETSET
 
    EvtVector4R p4string( 0.0, 0.0, 0.0, 0.0 );
 
    for ( i = 0; i < numparton; i++ ) { p4string += p4[partonindex[i]]; }
 
    int nprimary = 1;
    type[0]      = STRNG;
    for ( i = 0; i < numstable; i++ )
    {
      if ( km[stableindex[i]] == 0 ) { type[nprimary++] = evtnumstable[i]; }
    }
 
    p->makeDaughters( nprimary, type );
 
    p->getDaug( 0 )->init( STRNG, p4string );
 
    EvtVector4R p4partons[10];
 
    for ( i = 0; i < numparton; i++ ) { p4partons[i] = p4[partonindex[i]]; }
 
    ( (EvtStringParticle*)p->getDaug( 0 ) )->initPartons( numparton, p4partons, evtnumparton );
 
    nprimary = 1;
 
    for ( i = 0; i < numstable; i++ )
    {
 
      if ( km[stableindex[i]] == 0 )
      { p->getDaug( nprimary++ )->init( evtnumstable[i], p4[stableindex[i]] ); }
    }
 
    int nsecond = 0;
    for ( i = 0; i < numstable; i++ )
    {
      if ( km[stableindex[i]] != 0 ) { type[nsecond++] = evtnumstable[i]; }
    }
 
    p->getDaug( 0 )->makeDaughters( nsecond, type );
 
    nsecond = 0;
    for ( i = 0; i < numstable; i++ )
    {
      if ( km[stableindex[i]] != 0 )
      {
        p4[stableindex[i]] = boostTo( p4[stableindex[i]], p4string );
        p->getDaug( 0 )->getDaug( nsecond )->init( evtnumstable[i], p4[stableindex[i]] );
        p->getDaug( 0 )->getDaug( nsecond )->setDiagonalSpinDensity();
        p->getDaug( 0 )->getDaug( nsecond )->decay();
        nsecond++;
      }
    }
 
    fixPolarizations( p );
 
    return;
  }
}

getName()

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

virtual

Implements EvtDecayBase.

Definition at line 124 of file EvtPythia.cc.

{ model_name = "PYTHIA"; }

init()

void EvtPythia::init ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 130 of file EvtPythia.cc.

                     {
 
  checkNArg( 1 );
 
  if ( getParentId().isAlias() )
  {
 
    report( ERROR, "EvtGen" ) << "EvtPythia finds that you are decaying the" << endl
                              << " aliased particle " << EvtPDL::name( getParentId() ).c_str()
                              << " with the Pythia model" << endl
                              << " this does not work, please modify decay table." << endl;
    report( ERROR, "EvtGen" ) << "Will terminate execution!" << endl;
    ::abort();
  }
 
  store( this );
}

initProbMax()

void EvtPythia::initProbMax ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 128 of file EvtPythia.cc.

{ noProbMax(); }

pythiacont()

void EvtPythia::pythiacont ( double * energy, int * ndaugjs, int * kf, double * px, double * py, double * pz, double * e )

static

Definition at line 173 of file EvtPythia.cc.

                                                    {
  pycontinuum_( energy, ndaugjs, kf, px, py, pz, e );
}

Referenced by EvtPycont::decay(), and EvtPyGaGa::decay().

pythiaInit()

void EvtPythia::pythiaInit ( int f )

static

Definition at line 909 of file EvtPythia.cc.

                                      {
 
  static int first = 1;
  if ( first )
  {
 
    first = 0;
 
    report( INFO, "EvtGen" ) << "Will initialize Pythia." << endl;
    for ( int i = 0; i < ncommand; i++ )
      pygive_( commands[i].c_str(), strlen( commands[i].c_str() ) );
 
    char fname[200];
 
    char hostBuffer[100];
 
    if ( gethostname( hostBuffer, 100 ) != 0 )
    {
      report( ERROR, "EvtGen" ) << " couldn't get hostname." << endl;
      strncpy( hostBuffer, "hostnameNotFound", 100 );
    }
 
    char pid[100];
 
    int thePid = getpid();
 
    if ( sprintf( pid, "%d", thePid ) == 0 )
    {
      report( ERROR, "EvtGen" ) << " couldn't get process ID." << endl;
      strncpy( pid, "666", 100 );
    }
 
    strcpy( fname, "jet.d-" );
    strcat( fname, hostBuffer );
    strcat( fname, "-" );
    strcat( fname, pid );
 
    MakePythiaFile( fname );
    evtpythiainit_( fname, strlen( fname ) );
    initpythia_( &dummy );
 
    if ( 0 == getenv( "EVTSAVEJETD" ) )
    {
      char delcmd[300];
      strcpy( delcmd, "rm -f " );
      strcat( delcmd, fname );
      system( delcmd );
    }
 
    report( INFO, "EvtGen" ) << "Done initializing Pythia." << endl;
  }
}

Referenced by EvtPycont::decay(), EvtPyGaGa::decay(), and decay().

---

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

• EvtPythia.hh
• EvtPythia.cc