Decayer to simulate J/Psi polarization. (E. Scomparini)
[u/mrichter/AliRoot.git] / EVGEN / AliDecayerPolarized.h
CommitLineData
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
12#include "AliDecayer.h"
13
14
15class 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