Pythia6 Class Reference

#include <Pythia6.hh>

Public Member Functions
	Pythia6 ()
virtual	~Pythia6 ()
int	Pycomp (int kf)
void	Py1ent (int line, int kf, double pe, double theta, double phi)
ParticleVector *	ImportParticles ()
int	ImportParticles (ParticleVector *particles, const char *option="")
void	SetMSTJ (int i, int m)
int	GetMDCY (int i, int j)
int	GetKFDP (int i, int j)
void	SetMDCY (int i, int j, int m)
void	SetMDME (int i, int j, int m)

Static Public Member Functions
static Pythia6 *	Instance ()

Detailed Description

The C++ interface class to Pythia6

According to TPythia6 class from Root: (The TPythia6 class is an interface class to F77 routines in Pythia6 // CERNLIB event generators, written by T.Sjostrand.) http://root.cern.ch/ see http://root.cern.ch/root/License.html

The complete Pythia6 documentation can be found at: http://home.thep.lu.se/~torbjorn/pythiaaux/recent.html

Definition at line 163 of file Pythia6.hh.

Constructor & Destructor Documentation

Pythia6::Pythia6

(

)

Pythia6 constructor: creates a vector of Pythia6Particle in which it will store all particles. Note that there may be only one functional Pythia6 object at a time, so it's not use to create more than one instance of it.

Definition at line 128 of file Pythia6.cc.

References pythia6_common_address().

Referenced by Instance().

  : fParticles(0),
    fPyjets(0),
    fPydat1(0),
    fPydat3(0)
{
/// Pythia6 constructor: creates a vector of Pythia6Particle in which it will 
/// store all particles. Note that there may be only one functional Pythia6 
/// object at a time, so it's not use to create more than one instance of it.
  
   // Protect against multiple objects.   All access should be via the
   // Instance member function. 
   if ( fgInstance ) {
      std::cerr << "There's already an instance of Pythia6" << std::endl;
      exit (1);
   }   
  
   fParticles = new ParticleVector();

   // initialize common-blocks
   fPyjets = (Pyjets_t*) pythia6_common_address("PYJETS");
   fPydat1 = (Pydat1_t*) pythia6_common_address("PYDAT1");
   fPydat3 = (Pydat3_t*) pythia6_common_address("PYDAT3");
}

Pythia6::~Pythia6 ( )

virtual

Destroy the object, delete and dispose all Pythia6Particles currently on list.

Definition at line 155 of file Pythia6.cc.

{
/// Destroy the object, delete and dispose all Pythia6Particles currently on 
/// list.

   if ( fParticles ) {
      ParticleVector::const_iterator it;
      for ( it = fParticles->begin(); it != fParticles->end(); it++ )
        delete  *it;
      delete fParticles;
   }
}

Member Function Documentation

int Pythia6::GetKFDP ( int  i, int  j  )

inline

Definition at line 187 of file Pythia6.hh.

References Pydat3_t::KFDP.

{ return fPydat3->KFDP[j-1][i-1]; }

int Pythia6::GetMDCY ( int  i, int  j  )

inline

Definition at line 186 of file Pythia6.hh.

References Pydat3_t::MDCY.

{ return fPydat3->MDCY[j-1][i-1]; }

ParticleVector* Pythia6::ImportParticles

(

)

int Pythia6::ImportParticles	(	ParticleVector *	particles,
		const char *	option = ""
	)

Default primary creation method. It reads the /HEPEVT/ common block which has been filled by the GenerateEvent method. If the event generator does not use the HEPEVT common block, This routine has to be overloaded by the subclasses. The function loops on the generated particles and store them in the TClonesArray pointed by the argument particles. The default action is to store only the stable particles (ISTHEP = 1) This can be demanded explicitly by setting the option = "Final" If the option = "All", all the particles are stored.

Definition at line 204 of file Pythia6.cc.

References Pyjets_t::K, Pyjets_t::N, Pyjets_t::P, and Pyjets_t::V.

{
///  Default primary creation method. It reads the /HEPEVT/ common block which
///  has been filled by the GenerateEvent method. If the event generator does
///  not use the HEPEVT common block, This routine has to be overloaded by
///  the subclasses.
///  The function loops on the generated particles and store them in
///  the TClonesArray pointed by the argument particles.
///  The default action is to store only the stable particles (ISTHEP = 1)
///  This can be demanded explicitly by setting the option = "Final"
///  If the option = "All", all the particles are stored.

   if ( particles == 0 ) return 0;
   
   ParticleVector::const_iterator it;
   for ( it = particles->begin(); it != particles->end(); it++ )
     delete  *it;
   particles->clear();
   
   int numpart = fPyjets->N;
   int nparts=0;
   if (!strcmp(option,"") || !strcmp(option,"Final")) {
      for (int i = 0; i<numpart; i++) {

        if (fPyjets->K[0][i] == 1) {
          //
          //  Use the common block values for the TParticle constructor
          //
          particles->push_back(
            new Pythia6Particle(
                            fPyjets->K[0][i] ,
                            fPyjets->K[1][i] ,
                            fPyjets->K[2][i] ,
                            fPyjets->K[3][i] ,
                            fPyjets->K[4][i] ,
                            fPyjets->P[0][i] ,
                            fPyjets->P[1][i] ,
                            fPyjets->P[2][i] ,
                            fPyjets->P[3][i] ,
                            fPyjets->P[4][i] ,
                            fPyjets->V[0][i] ,
                            fPyjets->V[1][i] ,
                            fPyjets->V[2][i] ,
                            fPyjets->V[3][i] ,
                            fPyjets->V[4][i]));

          //     if(gDebug) printf("%d %d %d! ",i,fPyjets->K[1][i],numpart);
          nparts++;
       }
     }
   } 
   else if (!strcmp(option,"All")) {
      for (int i = 0; i<numpart; i++) {
          particles->push_back(
            new Pythia6Particle(
                            fPyjets->K[0][i] ,
                            fPyjets->K[1][i] ,
                            fPyjets->K[2][i] ,
                            fPyjets->K[3][i] ,
                            fPyjets->K[4][i] ,
                            fPyjets->P[0][i] ,
                            fPyjets->P[1][i] ,
                            fPyjets->P[2][i] ,
                            fPyjets->P[3][i] ,
                            fPyjets->P[4][i] ,
                            fPyjets->V[0][i] ,
                            fPyjets->V[1][i] ,
                            fPyjets->V[2][i] ,
                            fPyjets->V[3][i] ,
                            fPyjets->V[4][i]));
      }
      nparts=numpart;
   }

   return nparts;
}

Pythia6 * Pythia6::Instance ( void  )

static

Static access method

Definition at line 117 of file Pythia6.cc.

References Pythia6().

{
/// Static access method

   if ( ! fgInstance ) fgInstance = new Pythia6();

   return fgInstance;
}

void Pythia6::Py1ent	(	int	line,
		int	kf,
		double	pe,
		double	theta,
		double	phi
	)

Add one entry to the event record, i.e. either a parton or a particle.

IP: normally line number for the parton/particle. There are two exceptions:

If IP = 0: line number 1 is used and PYEXEC is called. If IP < 0: line -IP is used, with status code K(-IP,2)=2 rather than 1; thus a parton system may be built up by filling all but the last parton of the system with IP < 0. KF: parton/particle flavour code (PDG code) PE: parton/particle energy. If PE is smaller than the mass, the parton/particle is taken to be at rest. THETA: PHI: polar and azimuthal angle for the momentum vector of the parton/particle.

Definition at line 179 of file Pythia6.cc.

References py1ent.

{
/// Add one entry to the event record, i.e. either a parton or a
/// particle. 
///
///  IP:   normally line number for the parton/particle. There are two
///        exceptions:
/// 
///        If IP = 0: line number 1 is used and PYEXEC is called. 
///        If IP < 0: line -IP is used, with status code K(-IP,2)=2
///                   rather than 1; thus a parton system may be built
///                   up by filling all but the last parton of the
///                   system with IP < 0.   
///  KF:   parton/particle flavour code (PDG code)
///  PE:   parton/particle energy. If PE is smaller than the mass,
///        the parton/particle is taken to be at rest.  
///  THETA:
///  PHI:  polar and azimuthal angle for the momentum vector of the
///        parton/particle. 

   py1ent(ip, kf, pe, theta, phi);
}

int Pythia6::Pycomp

(

int

kf

)

Interface with fortran routine pycomp

Definition at line 170 of file Pythia6.cc.

References pycomp.

{
/// Interface with fortran routine pycomp

   return pycomp(&kf);
}

void Pythia6::SetMDCY ( int  i, int  j, int  m  )

inline

Definition at line 188 of file Pythia6.hh.

References python.hepunit::m, and Pydat3_t::MDCY.

{ fPydat3->MDCY[j-1][i-1] = m; }

void Pythia6::SetMDME ( int  i, int  j, int  m  )

inline

Definition at line 189 of file Pythia6.hh.

References python.hepunit::m, and Pydat3_t::MDME.

{ fPydat3->MDME[j-1][i-1] = m; }

void Pythia6::SetMSTJ ( int  i, int  m  )

inline

Definition at line 182 of file Pythia6.hh.

References python.hepunit::m, and Pydat1_t::MSTJ.

{ fPydat1->MSTJ[i-1] = m; }

---

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

• Pythia6.hh
• Pythia6.cc