Modifications needed to use EXODUS Decayer
[u/mrichter/AliRoot.git] / PYTHIA6 / AliPythia.h
CommitLineData
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>
e2de0ce1 11#include "PythiaProcesses.h"
03358a32 12#include "AliOmegaDalitz.h"
d7d0c637 13#include "AliDecayerExodus.h"
0f482ae4 14class AliFastGlauber;
15class AliQuenchingWeights;
16
8d2cd130 17class AliPythia : public TPythia6, public AliRndm
18{
19
20 public:
21 virtual ~AliPythia(){;}
22 // convert to compressed code and print result (for debugging only)
23 virtual Int_t CheckedLuComp(Int_t kf);
24 // Pythia initialisation for selected processes
25 virtual void ProcInit
efe3b1cd 26 (Process_t process, Float_t energy, StrucFunc_t strucfunc, Int_t itune = -1);
8d2cd130 27 // treat protons as inside nuclei
20e47f08 28 virtual void SetNuclei(Int_t a1, Int_t a2, Int_t pdf);
e28ccdaf 29 // Set colliding nuclei ("p","n",...)
30 virtual void SetCollisionSystem(TString projectile, TString target) { fProjectile = projectile; fTarget = target; }
8d2cd130 31 // Print particle properties
32 virtual void PrintParticles();
33 virtual void ResetDecayTable();
34 virtual void SetDecayTable();
16a82508 35 virtual void Pyevnw();
8d2cd130 36 virtual void Pycell(Int_t& nclus);
37 virtual void Pyclus(Int_t& nclus);
452af8c7 38 virtual void Pyshow(Int_t ip1, Int_t ip2, Double_t qmax);
cd07c39b 39 virtual void Pyshowq(Int_t ip1, Int_t ip2, Double_t qmax);
452af8c7 40 virtual void Pyrobo(Int_t imi, Int_t ima, Double_t the, Double_t phi, Double_t bex, Double_t bey, Double_t bez);
694b39f9 41 virtual void Pytune(Int_t itune);
9b61ba2a 42 virtual void Py2ent(Int_t idx, Int_t pdg1, Int_t pdg2, Double_t p);
32c8e463 43 virtual void InitQuenching(Float_t bmin, Float_t bmax, Float_t k, Int_t iECMethod, Float_t zmax = 0.97, Int_t ngmax = 30);
12cb0bc0 44 virtual void SetPyquenParameters(Double_t t0, Double_t tau0, Int_t nf, Int_t iengl, Int_t iangl);
992f2843 45 virtual void Pyquen(Double_t a, Int_t ibf, Double_t b);
6c43eccb 46 virtual void Qpygin0();
b280c4cc 47 virtual void GetQuenchingParameters(Double_t& xp, Double_t& yp, Double_t z[4]);
8d2cd130 48 // return instance of the singleton
49 static AliPythia* Instance();
452af8c7 50 virtual void Quench();
03358a32 51 void DalitzDecays();
d7d0c637 52 // Dalitz and resonance decays from EXODUS
53 void PizeroDalitz();
54 void EtaDalitz();
55 void RhoDirect();
56 void OmegaDalitz();
57 void OmegaDirect();
58 void EtaprimeDalitz();
59 void PhiDalitz();
60 void PhiDirect();
61 void JPsiDirect();
e8a8adcd 62 // Assignment Operator
63 AliPythia & operator=(const AliPythia & rhs);
64 void Copy(TObject&) const;
8d2cd130 65 protected:
0f482ae4 66 Process_t fProcess; // Process type
67 Float_t fEcms; // Centre of mass energy
68 StrucFunc_t fStrucFunc; // Structure function
e28ccdaf 69 TString fProjectile; // Projectile
70 TString fTarget; // Target
0f482ae4 71 Int_t fDefMDCY[501]; // ! Default decay switches per particle
72 Int_t fDefMDME[2001]; // ! Default decay switches per mode
b280c4cc 73 Double_t fZQuench[4]; // ! Quenching fractions for this even
74 Double_t fXJet; // ! Jet production point X
75 Double_t fYJet; // ! Jet production point Y
32c8e463 76 Int_t fNGmax; // Maximum number of radiated gluons in quenching
77 Float_t fZmax; // Maximum energy loss in quenching
0f482ae4 78 AliFastGlauber* fGlauber; // ! The Glauber model
79 AliQuenchingWeights* fQuenchingWeights; // ! The Quenching Weights model
0bd3d7c5 80 Int_t fItune; // ! Pythia tune
03358a32 81 AliOmegaDalitz fOmegaDalitz; // ! omega dalitz decayer
d7d0c637 82 AliDecayerExodus fExodus; // ! EXODUS decayer
0f482ae4 83 static AliPythia* fgAliPythia; // Pointer to single instance
8d2cd130 84 private:
85 AliPythia();
e8a8adcd 86 AliPythia(const AliPythia& pythia);
3dc3ec94 87 void ConfigHeavyFlavor();
e0e89f40 88 void AtlasTuning();
b4d68b2a 89 void AtlasTuningMC09();
8d2cd130 90 ClassDef(AliPythia,1) //ALICE UI to PYTHIA
91};
92
93#endif
94
95
96
90d7b703 97
589380c6 98