| | 1 | #ifndef ALIGENMUONLIB_H |
| | 2 | #define ALIGENMUONLIB_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 | // Library class for particle pt and y distributions used for |
| | 9 | // muon spectrometer simulations. |
| | 10 | // To be used with AliGenParam. |
| | 11 | // |
| | 12 | // andreas.morsch@cern.ch |
| | 13 | |
| | 14 | #include "AliGenLib.h" |
| | 15 | |
| | 16 | class AliGenMUONlib : |
| | 17 | public AliGenLib |
| | 18 | { |
| | 19 | public: |
| | 20 | enum constants{kPhi, kOmega, kEta, kJpsi, kJpsiFamily, kPsiP, kJpsiFromB, kUpsilon, kUpsilonFamily, |
| | 21 | kUpsilonP, kUpsilonPP, kCharm, kBeauty, kPion, kKaon}; |
| | 22 | |
| | 23 | GenFunc GetPt(Int_t param, const char* tname=0) const; |
| | 24 | GenFunc GetY (Int_t param, const char* tname=0) const; |
| | 25 | GenFuncIp GetIp(Int_t param, const char* tname=0) const; |
| | 26 | private: |
| | 27 | |
| | 28 | // pions |
| | 29 | static Double_t PtPion(Double_t *px, Double_t *dummy); |
| | 30 | static Double_t PtScal(Double_t pt, Int_t np); |
| | 31 | static Double_t YPion( Double_t *py, Double_t *dummy); |
| | 32 | static Int_t IpPion(TRandom *ran); |
| | 33 | // kaons |
| | 34 | static Double_t PtKaon(Double_t *px, Double_t *dummy); |
| | 35 | static Double_t YKaon( Double_t *py, Double_t *dummy); |
| | 36 | static Int_t IpKaon(TRandom *ran); |
| | 37 | // Phi |
| | 38 | static Double_t PtPhi( Double_t *px, Double_t *dummy); |
| | 39 | static Double_t YPhi( Double_t *px, Double_t *dummy); |
| | 40 | static Int_t IpPhi(TRandom *ran); |
| | 41 | // Omega |
| | 42 | static Double_t PtOmega( Double_t *px, Double_t *dummy); |
| | 43 | static Double_t YOmega( Double_t *px, Double_t *dummy); |
| | 44 | static Int_t IpOmega(TRandom *ran); |
| | 45 | // Eta |
| | 46 | static Double_t PtEta( Double_t *px, Double_t *dummy); |
| | 47 | static Double_t YEta( Double_t *px, Double_t *dummy); |
| | 48 | static Int_t IpEta(TRandom *ran); |
| | 49 | // J/Psi |
| | 50 | static Double_t PtJpsi( Double_t *px, Double_t *dummy); |
| | 51 | static Double_t PtJpsiCDFscaled( Double_t *px, Double_t *dummy); |
| | 52 | static Double_t PtJpsiCDFscaledPP( Double_t *px, Double_t *dummy); |
| | 53 | static Double_t YJpsi(Double_t *py, Double_t *dummy); |
| | 54 | static Double_t PtJpsiPbPb( Double_t *px, Double_t *dummy); |
| | 55 | static Double_t PtJpsiBPbPb( Double_t *px, Double_t *dummy); |
| | 56 | static Double_t YJpsiPbPb(Double_t *py, Double_t *dummy); |
| | 57 | static Double_t YJpsiCDFscaled(Double_t *py, Double_t *dummy); |
| | 58 | static Double_t YJpsiCDFscaledPP( Double_t *px, Double_t *dummy); |
| | 59 | static Double_t PtJpsiPP( Double_t *px, Double_t *dummy); |
| | 60 | static Double_t YJpsiPP(Double_t *py, Double_t *dummy); |
| | 61 | static Double_t YJpsiBPbPb(Double_t *py, Double_t *dummy); |
| | 62 | static Int_t IpJpsi(TRandom *ran); |
| | 63 | static Int_t IpJpsiFamily(TRandom *ran); |
| | 64 | static Int_t IpPsiP(TRandom *ran); |
| | 65 | static Double_t PtJpsiFlat( Double_t *px, Double_t *dummy ); |
| | 66 | static Double_t YJpsiFlat(Double_t *py, Double_t *dummy); |
| | 67 | |
| | 68 | // Upsilon |
| | 69 | static Double_t PtUpsilon( Double_t *px, Double_t *dummy ); |
| | 70 | static Double_t PtUpsilonCDFscaled( Double_t *px, Double_t *dummy ); |
| | 71 | static Double_t PtUpsilonCDFscaledPP( Double_t *px, Double_t *dummy ); |
| | 72 | static Double_t YUpsilon(Double_t *py, Double_t *dummy); |
| | 73 | static Double_t YUpsilonCDFscaled(Double_t *py, Double_t *dummy); |
| | 74 | static Double_t YUpsilonCDFscaledPP( Double_t *px, Double_t *dummy ); |
| | 75 | static Double_t PtUpsilonPbPb( Double_t *px, Double_t *dummy ); |
| | 76 | static Double_t YUpsilonPbPb(Double_t *py, Double_t *dummy); |
| | 77 | static Double_t PtUpsilonPP( Double_t *px, Double_t *dummy ); |
| | 78 | static Double_t YUpsilonPP(Double_t *py, Double_t *dummy); |
| | 79 | static Int_t IpUpsilon(TRandom *ran); |
| | 80 | static Int_t IpUpsilonFamily(TRandom *ran); |
| | 81 | static Int_t IpUpsilonP(TRandom *ran); |
| | 82 | static Int_t IpUpsilonPP(TRandom *ran); |
| | 83 | static Double_t PtUpsilonFlat( Double_t *px, Double_t *dummy ); |
| | 84 | static Double_t YUpsilonFlat(Double_t *py, Double_t *dummy); |
| | 85 | // |
| | 86 | // Charm |
| | 87 | static Double_t PtCharm( Double_t *px, Double_t *dummy ); |
| | 88 | static Double_t PtCharmCentral( Double_t *px, Double_t *dummy ); |
| | 89 | static Double_t YCharm(Double_t *py, Double_t *dummy); |
| | 90 | static Int_t IpCharm(TRandom *ran); |
| | 91 | // |
| | 92 | // Beauty |
| | 93 | static Double_t PtBeauty( Double_t *px, Double_t *dummy ); |
| | 94 | static Double_t PtBeautyCentral( Double_t *px, Double_t *dummy ); |
| | 95 | static Double_t YBeauty(Double_t *py, Double_t *dummy); |
| | 96 | static Int_t IpBeauty(TRandom *ran); |
| | 97 | // |
| | 98 | |
| | 99 | static Float_t Interpolate(Float_t x, Float_t* y, Float_t x0, |
| | 100 | Float_t dx, |
| | 101 | Int_t n, Int_t no); |
| | 102 | |
| | 103 | ClassDef(AliGenMUONlib,0) // Library providing y and pT parameterisations |
| | 104 | }; |
| | 105 | #endif |
| | 106 | |
| | 107 | |
| | 108 | |
| | 109 | |
| | 110 | |
| | 111 | |
| | 112 | |
git://git.uio.no / u / mrichter / AliRoot.git / blame_incremental