[u/mrichter/AliRoot.git] / EVGEN / AliPythia.h

#ifndef ALIPYTHIA_H
#define ALIPYTHIA_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

/* $Id$ */

#include <TPythia6.h>
#include <AliRndm.h>

typedef enum
{kPyCharm, kPyBeauty, kPyCharmUnforced, kPyBeautyUnforced, kPyJpsi, kPyJpsiChi, kPyMb, kPyJets, kPyDirectGamma, kPyCharmPbMNR}
Process_t;

typedef enum
{
    kDO_Set_1=1006,
    kGRV_LO=5005,
    kGRV_HO=5006,
    kMRS_D_minus=3031,
    kMRS_D0=3030,
    kMRS_G=3041,
    kCTEQ_2pM=4024,
    kCTEQ_4L=4032,
    kCTEQ_4M=4034,
    kMRST_cg_LO=3072,
    kCTEQ_5L=4046,
    kGRV_LO_98=5012
}
StrucFunc_t;

class AliPythia : public TPythia6, public AliRndm
{

 public:
    virtual ~AliPythia(){;}
    // convert to compressed code and print result (for debugging only)
    virtual Int_t CheckedLuComp(Int_t kf);
    // Pythia initialisation for selected processes
    virtual void ProcInit
	(Process_t process, Float_t energy, StrucFunc_t strucfunc);
    // treat protons as inside nuclei
    virtual void    SetNuclei(Int_t a1, Int_t a2);
    // Print particle properties
    virtual void PrintParticles();
    virtual void ResetDecayTable();
    virtual void SetDecayTable();
    // return instance of the singleton
    static  AliPythia* Instance();

 protected:
    Process_t     fProcess;           // Process type
    Float_t       fEcms;              // Centre of mass energy
    StrucFunc_t   fStrucFunc;         // Structure function
    Int_t         fDefMDCY[501];      //  ! Default decay switches per particle
    Int_t         fDefMDME[2000];     //  ! Default decay switches per mode
    static AliPythia*    fgAliPythia; // Pointer to single instance
 private: 
    AliPythia();

    ClassDef(AliPythia,1) //ALICE UI to PYTHIA
};

#endif
Commit	Line	Data
f87cfe57	1	#ifndef ALIPYTHIA_H
f87cfe57	2	#define ALIPYTHIA_H
3da30618	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6	/* $Id$ */
fe4da5cc	7
95b811fe	8	#include <TPythia6.h>
3356c022	9	#include <AliRndm.h>
f1a48a38	10
f1a48a38	11	typedef enum
fbd1348b	12	{kPyCharm, kPyBeauty, kPyCharmUnforced, kPyBeautyUnforced, kPyJpsi, kPyJpsiChi, kPyMb, kPyJets, kPyDirectGamma, kPyCharmPbMNR}
f1a48a38	13	Process_t;
	14
	15	typedef enum
	16	{
	17	kDO_Set_1=1006,
	18	kGRV_LO=5005,
	19	kGRV_HO=5006,
	20	kMRS_D_minus=3031,
	21	kMRS_D0=3030,
	22	kMRS_G=3041,
	23	kCTEQ_2pM=4024,
2afdd95f	24	kCTEQ_4L=4032,
c1562ef9	25	kCTEQ_4M=4034,
	26	kMRST_cg_LO=3072,
	27	kCTEQ_5L=4046,
	28	kGRV_LO_98=5012
f1a48a38	29	}
f1a48a38	30	StrucFunc_t;
fe4da5cc	31
3356c022	32	class AliPythia : public TPythia6, public AliRndm
fe4da5cc	33	{
f1a48a38	34
fe4da5cc	35	public:
75c6d54e	36	virtual ~AliPythia(){;}
fe4da5cc	37	// convert to compressed code and print result (for debugging only)
f87cfe57	38	virtual Int_t CheckedLuComp(Int_t kf);
fe4da5cc	39	// Pythia initialisation for selected processes
	40	virtual void ProcInit
	41	(Process_t process, Float_t energy, StrucFunc_t strucfunc);
811826d8	42	// treat protons as inside nuclei
811826d8	43	virtual void SetNuclei(Int_t a1, Int_t a2);
14ee1cd0	44	// Print particle properties
	45	virtual void PrintParticles();
	46	virtual void ResetDecayTable();
	47	virtual void SetDecayTable();
811826d8	48	// return instance of the singleton
95b811fe	49	static AliPythia* Instance();
95b811fe	50
f87cfe57	51	protected:
95b811fe	52	Process_t fProcess; // Process type
	53	Float_t fEcms; // Centre of mass energy
	54	StrucFunc_t fStrucFunc; // Structure function
14ee1cd0	55	Int_t fDefMDCY[501]; // ! Default decay switches per particle
14ee1cd0	56	Int_t fDefMDME[2000]; // ! Default decay switches per mode
95b811fe	57	static AliPythia* fgAliPythia; // Pointer to single instance
	58	private:
	59	AliPythia();
f87cfe57	60
3356c022	61	ClassDef(AliPythia,1) //ALICE UI to PYTHIA
fe4da5cc	62	};
	63
	64	#endif
	65
	66
	67