[u/mrichter/AliRoot.git] / EVGEN / AliDecayerPolarized.h

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

// Class to generate decay products for polarized heavy quarkonia

#include <TLorentzVector.h>
#include <TClonesArray.h>
#include <TF1.h>

#include "AliDecayer.h"

class TParticle;

class AliDecayerPolarized : public AliDecayer
{
 public:
    typedef enum { kNoPol = 0, kColSop = 1, kHelicity = 2} Polar_t;
    typedef enum { kElectron = 1, kMuon = 2} FinState_t;
  AliDecayerPolarized();
  AliDecayerPolarized(Double_t alpha, Polar_t systref, FinState_t decprod);
  AliDecayerPolarized(const AliDecayerPolarized &decayer);
  virtual ~AliDecayerPolarized();
  void SetPolDec(Double_t alpha=0) {fAlpha=alpha;}
  void SetPolRefSys(Polar_t systref=kColSop) {fSystRef=systref;}
  void SetDecProd(FinState_t decprod=kMuon) {fDecProd=decprod;}
  virtual void  Init(){;}
  virtual void  Decay(Int_t ipart, TLorentzVector *p);
  virtual Int_t ImportParticles(TClonesArray *part);
  // The following methods are dummy
  virtual void    SetForceDecay(Int_t type);
  virtual void    ForceDecay();
  virtual Float_t GetPartialBranchingRatio(Int_t ipart);
  virtual Float_t GetLifetime(Int_t kf);
  virtual void    ReadDecayTable();
    
 protected:
  Double_t fAlpha;       // Polarization parameter
  Polar_t fSystRef;      // Reference system for polarization
  FinState_t fDecProd;   // Choice of decay products
  TF1 *fPol;             // ! Angular distribution for decay products
  TParticle *fMother;    // ! Particle that has to be decayed
  TParticle *fDaughter1; // ! Decay product no. 1
  TParticle *fDaughter2; // ! Decay product no. 2

 private:
  void  Copy(TObject &decayer) const;
  AliDecayerPolarized &operator=(const AliDecayerPolarized &decayer) 
      {decayer.Copy(*this);return(*this);}
  
  ClassDef(AliDecayerPolarized,1) // Polarized 2-body quarkonium decay
};
#endif
Commit	Line	Data
c1e10a2a	1	#ifndef ALIDECAYERPOLARIZED_H
	2	#define ALIDECAYERPOLARIZED_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6	// Class to generate decay products for polarized heavy quarkonia
	7
	8	#include <TLorentzVector.h>
	9	#include <TClonesArray.h>
	10	#include <TF1.h>
	11
3cad937a	12	#include "AliDecayer.h"
c1e10a2a	13
e9912df8	14	class TParticle;
c1e10a2a	15
3cad937a	16	class AliDecayerPolarized : public AliDecayer
c1e10a2a	17	{
	18	public:
	19	typedef enum { kNoPol = 0, kColSop = 1, kHelicity = 2} Polar_t;
	20	typedef enum { kElectron = 1, kMuon = 2} FinState_t;
	21	AliDecayerPolarized();
	22	AliDecayerPolarized(Double_t alpha, Polar_t systref, FinState_t decprod);
1c56e311	23	AliDecayerPolarized(const AliDecayerPolarized &decayer);
c1e10a2a	24	virtual ~AliDecayerPolarized();
	25	void SetPolDec(Double_t alpha=0) {fAlpha=alpha;}
	26	void SetPolRefSys(Polar_t systref=kColSop) {fSystRef=systref;}
	27	void SetDecProd(FinState_t decprod=kMuon) {fDecProd=decprod;}
3cad937a	28	virtual void Init(){;}
3cad937a	29	virtual void Decay(Int_t ipart, TLorentzVector *p);
c1e10a2a	30	virtual Int_t ImportParticles(TClonesArray *part);
3cad937a	31	// The following methods are dummy
	32	virtual void SetForceDecay(Int_t type);
	33	virtual void ForceDecay();
	34	virtual Float_t GetPartialBranchingRatio(Int_t ipart);
	35	virtual Float_t GetLifetime(Int_t kf);
	36	virtual void ReadDecayTable();
c1e10a2a	37
	38	protected:
	39	Double_t fAlpha; // Polarization parameter
	40	Polar_t fSystRef; // Reference system for polarization
	41	FinState_t fDecProd; // Choice of decay products
	42	TF1 *fPol; // ! Angular distribution for decay products
	43	TParticle *fMother; // ! Particle that has to be decayed
	44	TParticle *fDaughter1; // ! Decay product no. 1
	45	TParticle *fDaughter2; // ! Decay product no. 2
78fd4e85	46
	47	private:
	48	void Copy(TObject &decayer) const;
	49	AliDecayerPolarized &operator=(const AliDecayerPolarized &decayer)
	50	{decayer.Copy(this);return(this);}
c1e10a2a	51
	52	ClassDef(AliDecayerPolarized,1) // Polarized 2-body quarkonium decay
	53	};
	54	#endif
	55
	56