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 */
6 // Class to generate decay products for polarized heavy quarkonia
8 #include <TLorentzVector.h>
9 #include <TClonesArray.h>
12 #include "AliDecayer.h"
15 class AliDecayerPolarized : public AliDecayer
18 typedef enum { kNoPol = 0, kColSop = 1, kHelicity = 2} Polar_t;
19 typedef enum { kElectron = 1, kMuon = 2} FinState_t;
20 AliDecayerPolarized();
21 AliDecayerPolarized(Double_t alpha, Polar_t systref, FinState_t decprod);
22 AliDecayerPolarized(const AliDecayerPolarized &decayer);
23 virtual ~AliDecayerPolarized();
24 void SetPolDec(Double_t alpha=0) {fAlpha=alpha;}
25 void SetPolRefSys(Polar_t systref=kColSop) {fSystRef=systref;}
26 void SetDecProd(FinState_t decprod=kMuon) {fDecProd=decprod;}
27 virtual void Init(){;}
28 virtual void Decay(Int_t ipart, TLorentzVector *p);
29 virtual Int_t ImportParticles(TClonesArray *part);
30 // The following methods are dummy
31 virtual void SetForceDecay(Int_t type);
32 virtual void ForceDecay();
33 virtual Float_t GetPartialBranchingRatio(Int_t ipart);
34 virtual Float_t GetLifetime(Int_t kf);
35 virtual void ReadDecayTable();
37 void Copy(TObject &decayer) const;
38 AliDecayerPolarized &operator=(const AliDecayerPolarized &decayer)
39 {decayer.Copy(*this);return(*this);}
42 Double_t fAlpha; // Polarization parameter
43 Polar_t fSystRef; // Reference system for polarization
44 FinState_t fDecProd; // Choice of decay products
45 TF1 *fPol; // ! Angular distribution for decay products
46 TParticle *fMother; // ! Particle that has to be decayed
47 TParticle *fDaughter1; // ! Decay product no. 1
48 TParticle *fDaughter2; // ! Decay product no. 2
50 ClassDef(AliDecayerPolarized,1) // Polarized 2-body quarkonium decay