[u/mrichter/AliRoot.git] / PYTHIA6 / 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>
#include <AliStructFuncType.h>
typedef enum
{kPyCharm, kPyBeauty, kPyCharmUnforced, kPyBeautyUnforced,
 kPyJpsi, kPyJpsiChi, kPyMb, kPyMbNonDiffr, kPyJets, kPyDirectGamma,
 kPyCharmPbPbMNR, kPyD0PbPbMNR, kPyDPlusPbPbMNR, kPyBeautyPbPbMNR,
 kPyCharmpPbMNR, kPyD0pPbMNR, kPyDPluspPbMNR, kPyBeautypPbMNR,
 kPyCharmppMNR, kPyD0ppMNR, kPyDPlusppMNR, kPyBeautyppMNR, kPyW,
 kPyOldUEQ2ordered, kPyOldUEQ2ordered2, kPyOldPopcorn}
Process_t;
/*
typedef enum
{
    kDOSet1     = 1006,
    kGRVLO      = 5005,
    kGRVHO      = 5006,
    kMRSDminus  = 3031,
    kMRSD0      = 3030,
    kMRSG       = 3041,
    kCTEQ2pM    = 4024,
    kCTEQ4L     = 4032,
    kCTEQ4M     = 4034,
    kMRSTcgLO   = 3072,
    kCTEQ5L     = 4046,
    kGRVLO98    = 5012
}
StrucFunc_t;
*/
class AliFastGlauber;
class AliQuenchingWeights;

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();
    virtual void Pycell(Int_t& nclus);
    virtual void Pyclus(Int_t& nclus);
    virtual void Pyshow(Int_t ip1, Int_t ip2, Double_t qmax);
    virtual void Pyrobo(Int_t imi, Int_t ima, Double_t the, Double_t phi, Double_t bex, Double_t bey, Double_t bez);
    virtual void InitQuenching(Float_t bmin, Float_t bmax, Float_t k, Int_t iECMethod);
    virtual void Pyquen(Double_t a, Int_t ibf, Double_t b);
    virtual void GetQuenchingParameters(Double_t& xp, Double_t& yp, Double_t z[4]);
    // return instance of the singleton
    static  AliPythia* Instance();
    virtual void Quench();
 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[2001];     //  ! Default decay switches per mode
    Double_t              fZQuench[4];        //  ! Quenching fractions for this even
    Double_t              fXJet;              //  ! Jet production point X
    Double_t              fYJet;              //  ! Jet production point Y
    AliFastGlauber*       fGlauber;           //  ! The Glauber model
    AliQuenchingWeights*  fQuenchingWeights;  //  ! The Quenching Weights model
    static AliPythia*     fgAliPythia;        // Pointer to single instance
 private: 
    AliPythia();

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

#endif
Commit	Line	Data
8d2cd130	1	#ifndef ALIPYTHIA_H
	2	#define ALIPYTHIA_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6	/* $Id$ */
	7
	8	#include <TPythia6.h>
	9	#include <AliRndm.h>
	10	#include <AliStructFuncType.h>
	11	typedef enum
	12	{kPyCharm, kPyBeauty, kPyCharmUnforced, kPyBeautyUnforced,
	13	kPyJpsi, kPyJpsiChi, kPyMb, kPyMbNonDiffr, kPyJets, kPyDirectGamma,
90d7b703	14	kPyCharmPbPbMNR, kPyD0PbPbMNR, kPyDPlusPbPbMNR, kPyBeautyPbPbMNR,
90d7b703	15	kPyCharmpPbMNR, kPyD0pPbMNR, kPyDPluspPbMNR, kPyBeautypPbMNR,
65f2626c	16	kPyCharmppMNR, kPyD0ppMNR, kPyDPlusppMNR, kPyBeautyppMNR, kPyW,
65f2626c	17	kPyOldUEQ2ordered, kPyOldUEQ2ordered2, kPyOldPopcorn}
8d2cd130	18	Process_t;
	19	/*
	20	typedef enum
	21	{
	22	kDOSet1 = 1006,
	23	kGRVLO = 5005,
	24	kGRVHO = 5006,
	25	kMRSDminus = 3031,
	26	kMRSD0 = 3030,
	27	kMRSG = 3041,
	28	kCTEQ2pM = 4024,
	29	kCTEQ4L = 4032,
	30	kCTEQ4M = 4034,
	31	kMRSTcgLO = 3072,
	32	kCTEQ5L = 4046,
	33	kGRVLO98 = 5012
	34	}
	35	StrucFunc_t;
	36	*/
0f482ae4	37	class AliFastGlauber;
	38	class AliQuenchingWeights;
	39
8d2cd130	40	class AliPythia : public TPythia6, public AliRndm
	41	{
	42
	43	public:
	44	virtual ~AliPythia(){;}
	45	// convert to compressed code and print result (for debugging only)
	46	virtual Int_t CheckedLuComp(Int_t kf);
	47	// Pythia initialisation for selected processes
	48	virtual void ProcInit
	49	(Process_t process, Float_t energy, StrucFunc_t strucfunc);
	50	// treat protons as inside nuclei
	51	virtual void SetNuclei(Int_t a1, Int_t a2);
	52	// Print particle properties
	53	virtual void PrintParticles();
	54	virtual void ResetDecayTable();
	55	virtual void SetDecayTable();
	56	virtual void Pycell(Int_t& nclus);
	57	virtual void Pyclus(Int_t& nclus);
452af8c7	58	virtual void Pyshow(Int_t ip1, Int_t ip2, Double_t qmax);
452af8c7	59	virtual void Pyrobo(Int_t imi, Int_t ima, Double_t the, Double_t phi, Double_t bex, Double_t bey, Double_t bez);
86b6ad68	60	virtual void InitQuenching(Float_t bmin, Float_t bmax, Float_t k, Int_t iECMethod);
992f2843	61	virtual void Pyquen(Double_t a, Int_t ibf, Double_t b);
b280c4cc	62	virtual void GetQuenchingParameters(Double_t& xp, Double_t& yp, Double_t z[4]);
8d2cd130	63	// return instance of the singleton
8d2cd130	64	static AliPythia* Instance();
452af8c7	65	virtual void Quench();
8d2cd130	66	protected:
0f482ae4	67	Process_t fProcess; // Process type
	68	Float_t fEcms; // Centre of mass energy
	69	StrucFunc_t fStrucFunc; // Structure function
	70	Int_t fDefMDCY[501]; // ! Default decay switches per particle
	71	Int_t fDefMDME[2001]; // ! Default decay switches per mode
b280c4cc	72	Double_t fZQuench[4]; // ! Quenching fractions for this even
	73	Double_t fXJet; // ! Jet production point X
	74	Double_t fYJet; // ! Jet production point Y
0f482ae4	75	AliFastGlauber* fGlauber; // ! The Glauber model
	76	AliQuenchingWeights* fQuenchingWeights; // ! The Quenching Weights model
	77	static AliPythia* fgAliPythia; // Pointer to single instance
8d2cd130	78	private:
	79	AliPythia();
	80
	81	ClassDef(AliPythia,1) //ALICE UI to PYTHIA
	82	};
	83
	84	#endif
	85
	86
	87
90d7b703	88
589380c6	89