Modifications needed to use EXODUS Decayer
[u/mrichter/AliRoot.git] / PYTHIA6 / AliDecayerPythia.h
1 #ifndef ALIDECAYERPYTHIA_H
2 #define ALIDECAYERPYTHIA_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 // Implementation of AliDecayer using Pythia
9 // Method forwarding to the AliPythia instance.
10 // Author: andreas.morsch@cern.ch
11
12 #include "AliDecayer.h"
13 class AliPythia;
14 class TClonesArrray;
15
16 class AliDecayerPythia :
17 public AliDecayer
18 {
19  public:
20     AliDecayerPythia();
21     AliDecayerPythia(const AliDecayerPythia &decayer);
22     //    
23     virtual         ~AliDecayerPythia(){;}
24     virtual void    Init();
25     virtual void    Decay(Int_t idpart, TLorentzVector *p);
26     virtual Int_t   ImportParticles(TClonesArray *particles);
27     virtual void    SetForceDecay(Decay_t decay) {fDecay=decay;}
28     virtual void    SetForceDecay(Int_t decay)
29       {SetForceDecay((Decay_t) decay);}
30     virtual void    ForceDecay();
31     virtual void    SetPatchOmegaDalitz() {fPatchOmegaDalitz = 1;}
32     virtual void    SetDecayerExodus()    {fDecayerExodus = 1;}
33     virtual void    HeavyFlavourOff() {fHeavyFlavour = kFALSE;}
34     virtual void    DecayLongLivedParticles()  {fLongLived    = kTRUE;}
35     virtual Float_t GetPartialBranchingRatio(Int_t ipart);
36     virtual Float_t GetLifetime(Int_t kf);
37     virtual void    SwitchOffBDecay();
38     virtual void    SwitchOffPi0() {fPi0 = 0;}
39     virtual void    WriteDecayTable();
40     virtual void    ReadDecayTable();
41     
42  private:
43     void     DefineParticles();
44     void     Lu1Ent(Int_t flag, Int_t idpart, 
45                     Double_t mom, Double_t theta, Double_t phi);
46     Int_t    CountProducts(Int_t channel, Int_t particle);
47     void     ForceParticleDecay(Int_t particle, Int_t product, Int_t mult);
48     void     ForceParticleDecay(Int_t particle, const Int_t* products, Int_t* mult, Int_t npart, Bool_t flag = 0);
49     void     ForceHadronicD(Int_t optUser4Bodies=1);
50     void     ForceOmega();
51     void     ForceLambda();
52     void     SwitchOffHeavyFlavour();
53     void     ForceBeautyUpgrade();
54     Float_t  GetBraPart(Int_t kf);
55     void     Copy(TObject &decayer) const;
56
57     AliDecayerPythia &operator=(const AliDecayerPythia &decayer) 
58     {decayer.Copy(*this);return(*this);}
59     
60     
61  private:
62     AliPythia*  fPythia;          //! Pointer to AliPythia
63     Decay_t     fDecay;           //  Forced decay mode
64     Float_t     fBraPart[501];    //! Branching ratios
65     Bool_t      fHeavyFlavour;    //! Flag for heavy flavors
66     Bool_t      fLongLived;       //! Flag for long lived particle decay
67     Bool_t      fPatchOmegaDalitz;//! Flag to patch the omega Dalitz decays 
68     Bool_t      fDecayerExodus;  //! Flag for EXODUS decayer
69     Bool_t      fPi0;             //! Flag for pi0 decay 
70     static Bool_t fgInit;         //! initialization flag 
71     
72     ClassDef(AliDecayerPythia, 5) // AliDecayer implementation using Pythia  
73 };
74 #endif
75
76
77
78
79
80
81