Decayer to simulate J/Psi polarization. (E. Scomparini)
[u/mrichter/AliRoot.git] / EVGEN / AliDecayerPolarized.h
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
12 #include "AliDecayer.h"
13
14
15 class AliDecayerPolarized : public AliDecayer
16 {
17  public:
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   virtual ~AliDecayerPolarized();
23   void SetPolDec(Double_t alpha=0) {fAlpha=alpha;}
24   void SetPolRefSys(Polar_t systref=kColSop) {fSystRef=systref;}
25   void SetDecProd(FinState_t decprod=kMuon) {fDecProd=decprod;}
26   virtual void Init(){;}
27   virtual void Decay(Int_t ipart, TLorentzVector *p);
28   virtual Int_t ImportParticles(TClonesArray *part);
29   AliDecayerPolarized(const AliDecayerPolarized &decayer):AliDecayer(decayer)
30       {decayer.Copy(*this);}
31   virtual AliDecayerPolarized &operator=(const AliDecayerPolarized &decayer) 
32       {decayer.Copy(*this);return(*this);}
33     
34  protected:
35   Double_t fAlpha;       // Polarization parameter
36   Polar_t fSystRef;      // Reference system for polarization
37   FinState_t fDecProd;   // Choice of decay products
38   TF1 *fPol;             // ! Angular distribution for decay products
39   TParticle *fMother;    // ! Particle that has to be decayed
40   TParticle *fDaughter1; // ! Decay product no. 1
41   TParticle *fDaughter2; // ! Decay product no. 2
42   
43   ClassDef(AliDecayerPolarized,1) // Polarized 2-body quarkonium decay
44 };
45 #endif
46
47