Option to switch pi0 decay off
[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    HeavyFlavourOff() {fHeavyFlavour = kFALSE;}
33     virtual void    DecayLongLivedParticles()  {fLongLived    = kTRUE;}
34     virtual Float_t GetPartialBranchingRatio(Int_t ipart);
35     virtual Float_t GetLifetime(Int_t kf);
36     virtual void    SwitchOffBDecay();
37     virtual void    SwitchOffPi0() {fPi0 = 0;}
38     virtual void    WriteDecayTable();
39     virtual void    ReadDecayTable();
40     
41  private:
42     void     DefineParticles();
43     void     Lu1Ent(Int_t flag, Int_t idpart, 
44                     Double_t mom, Double_t theta, Double_t phi);
45     Int_t    CountProducts(Int_t channel, Int_t particle);
46     void     ForceParticleDecay(Int_t particle, Int_t product, Int_t mult);
47     void     ForceParticleDecay(Int_t particle, const Int_t* products, Int_t* mult, Int_t npart, Bool_t flag = 0);
48     void     ForceHadronicD(Int_t optUser4Bodies=1);
49     void     ForceOmega();
50     void     ForceLambda();
51     void     SwitchOffHeavyFlavour();
52     void     ForceBeautyUpgrade();
53     Float_t  GetBraPart(Int_t kf);
54     void     Copy(TObject &decayer) const;
55
56     AliDecayerPythia &operator=(const AliDecayerPythia &decayer) 
57     {decayer.Copy(*this);return(*this);}
58     
59     
60  private:
61     AliPythia*  fPythia;          //! Pointer to AliPythia
62     Decay_t     fDecay;           //  Forced decay mode
63     Float_t     fBraPart[501];    //! Branching ratios
64     Bool_t      fHeavyFlavour;    //! Flag for heavy flavors
65     Bool_t      fLongLived;       //! Flag for long lived particle decay
66     Bool_t      fPatchOmegaDalitz;//! Flag to patch the omega Dalitz decays 
67     Bool_t      fPi0;             //! Flag for pi0 decay 
68     static Bool_t fgInit;         //! initialization flag 
69     
70     ClassDef(AliDecayerPythia, 4) // AliDecayer implementation using Pythia  
71 };
72 #endif
73
74
75
76
77
78
79