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