Pythia6 Class Reference

#include <Pythia6.hh>

Collaboration diagram for Pythia6:

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

Private Attributes
ParticleVector *	fParticles
Pyjets_t *	fPyjets
Pydat1_t *	fPydat1
Pydat3_t *	fPydat3

Static Private Attributes
static Pythia6 *	fgInstance = 0

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

(

)

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.

  : fParticles(0),
    fPyjets(0),
    fPydat1(0),
    fPydat3(0)
{
  
   // 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");
}

Here is the call graph for this function:

Here is the caller graph for this function:

~Pythia6()

Pythia6::~Pythia6 ( )

virtual

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

Definition at line 155 of file Pythia6.cc.

{

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

Member Function Documentation

GetKFDP()

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

inline

Definition at line 187 of file Pythia6.hh.

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

Here is the caller graph for this function:

GetMDCY()

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

inline

Definition at line 186 of file Pythia6.hh.

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

Here is the caller graph for this function:

ImportParticles() [1/2]

ParticleVector* Pythia6::ImportParticles

(

)

Here is the caller graph for this function:

ImportParticles() [2/2]

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.

{

   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;
}

Instance()

Pythia6 * Pythia6::Instance ( void  )

static

Static access method

Definition at line 117 of file Pythia6.cc.

{

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

   return fgInstance;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Py1ent()

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.

{

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

Here is the caller graph for this function:

Pycomp()

int Pythia6::Pycomp

(

int

kf

)

Interface with fortran routine pycomp

Definition at line 170 of file Pythia6.cc.

{

   return pycomp(&kf);
}

Here is the caller graph for this function:

SetMDCY()

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

inline

Definition at line 188 of file Pythia6.hh.

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

Here is the caller graph for this function:

SetMDME()

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

inline

Definition at line 189 of file Pythia6.hh.

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

Here is the caller graph for this function:

SetMSTJ()

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

inline

Definition at line 182 of file Pythia6.hh.

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

Here is the caller graph for this function:

Member Data Documentation

fgInstance

Pythia6 * Pythia6::fgInstance = 0

staticprivate

Definition at line 192 of file Pythia6.hh.

fParticles

ParticleVector* Pythia6::fParticles

private

Definition at line 194 of file Pythia6.hh.

fPydat1

Pydat1_t* Pythia6::fPydat1

private

Definition at line 196 of file Pythia6.hh.

fPydat3

Pydat3_t* Pythia6::fPydat3

private

Definition at line 197 of file Pythia6.hh.

fPyjets

Pyjets_t* Pythia6::fPyjets

private

Definition at line 195 of file Pythia6.hh.

---

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

• Pythia6.hh
• Pythia6.cc