[u/mrichter/AliRoot.git] / PYTHIA6 / AliDecayerPythia.h

#ifndef ALIDECAYERPYTHIA_H
#define ALIDECAYERPYTHIA_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

/* $Id$ */

// Implementation of AliDecayer using Pythia
// Method forwarding to the AliPythia instance.
// Author: andreas.morsch@cern.ch

#include "AliDecayer.h"
class AliPythia;
class TClonesArrray;

class AliDecayerPythia :
public AliDecayer
{
 public:
    AliDecayerPythia();
    AliDecayerPythia(const AliDecayerPythia &decayer);
    //    
    virtual         ~AliDecayerPythia(){;}
    virtual void    Init();
    virtual void    Decay(Int_t idpart, TLorentzVector *p);
    virtual Int_t   ImportParticles(TClonesArray *particles);
    virtual void    SetForceDecay(Decay_t decay) {fDecay=decay;}
    virtual void    SetForceDecay(Int_t decay)
      {SetForceDecay((Decay_t) decay);}
    virtual void    ForceDecay();
    virtual void    SetPatchOmegaDalitz() {fPatchOmegaDalitz = 1;}
    virtual void    HeavyFlavourOff() {fHeavyFlavour = kFALSE;}
    virtual void    DecayLongLivedParticles()  {fLongLived    = kTRUE;}
    virtual Float_t GetPartialBranchingRatio(Int_t ipart);
    virtual Float_t GetLifetime(Int_t kf);
    virtual void    SwitchOffBDecay();
    virtual void    SwitchOffPi0() {fPi0 = 0;}
    virtual void    WriteDecayTable();
    virtual void    ReadDecayTable();
    
 private:
    void     DefineParticles();
    void     Lu1Ent(Int_t flag, Int_t idpart, 
		    Double_t mom, Double_t theta, Double_t phi);
    Int_t    CountProducts(Int_t channel, Int_t particle);
    void     ForceParticleDecay(Int_t particle, Int_t product, Int_t mult);
    void     ForceParticleDecay(Int_t particle, const Int_t* products, Int_t* mult, Int_t npart, Bool_t flag = 0);
    void     ForceHadronicD(Int_t optUser4Bodies=1);
    void     ForceOmega();
    void     ForceLambda();
    void     SwitchOffHeavyFlavour();
    void     ForceBeautyUpgrade();
    Float_t  GetBraPart(Int_t kf);
    void     Copy(TObject &decayer) const;

    AliDecayerPythia &operator=(const AliDecayerPythia &decayer) 
    {decayer.Copy(*this);return(*this);}
    
    
 private:
    AliPythia*  fPythia;          //! Pointer to AliPythia
    Decay_t     fDecay;           //  Forced decay mode
    Float_t     fBraPart[501];    //! Branching ratios
    Bool_t      fHeavyFlavour;    //! Flag for heavy flavors
    Bool_t      fLongLived;       //! Flag for long lived particle decay
    Bool_t      fPatchOmegaDalitz;//! Flag to patch the omega Dalitz decays 
    Bool_t      fPi0;             //! Flag for pi0 decay 
    static Bool_t fgInit;         //! initialization flag 
    
    ClassDef(AliDecayerPythia, 4) // AliDecayer implementation using Pythia  
};
#endif
Commit	Line	Data
	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