| | 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, kChic, kChic0, kChic1, kChic2 }; |
| | 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(const Double_t *px, const Double_t *dummy); |
| | 30 | static Double_t PtScal(Double_t pt, Int_t np); |
| | 31 | static Double_t YPion( const Double_t *py, const Double_t *dummy); |
| | 32 | static Int_t IpPion(TRandom *ran); |
| | 33 | // kaons |
| | 34 | static Double_t PtKaon(const Double_t *px, const Double_t *dummy); |
| | 35 | static Double_t YKaon( const Double_t *py, const Double_t *dummy); |
| | 36 | static Int_t IpKaon(TRandom *ran); |
| | 37 | // Phi |
| | 38 | static Double_t PtPhi( const Double_t *px, const Double_t *dummy); |
| | 39 | static Double_t YPhi( const Double_t *px, const Double_t *dummy); |
| | 40 | static Int_t IpPhi(TRandom *ran); |
| | 41 | // Omega |
| | 42 | static Double_t PtOmega( const Double_t *px, const Double_t *dummy); |
| | 43 | static Double_t YOmega( const Double_t *px, const Double_t *dummy); |
| | 44 | static Int_t IpOmega(TRandom *ran); |
| | 45 | // Eta |
| | 46 | static Double_t PtEta( const Double_t *px, const Double_t *dummy); |
| | 47 | static Double_t YEta( const Double_t *px, const Double_t *dummy); |
| | 48 | static Int_t IpEta(TRandom *ran); |
| | 49 | // J/Psi |
| | 50 | static Double_t PtJpsi( const Double_t *px, const Double_t *dummy); |
| | 51 | static Double_t PtJpsiCDFscaled( const Double_t *px, const Double_t *dummy); |
| | 52 | static Double_t PtJpsiCDFscaledPP( const Double_t *px, const Double_t *dummy); |
| | 53 | static Double_t PtJpsiCDFscaledPP10( const Double_t *px, const Double_t *dummy); |
| | 54 | static Double_t PtJpsiCDFscaledold( const Double_t *px, const Double_t *dummy); |
| | 55 | static Double_t PtJpsiCDFscaledPPold( const Double_t *px, const Double_t *dummy); |
| | 56 | static Double_t YJpsi(const Double_t *py, const Double_t *dummy); |
| | 57 | static Double_t PtJpsiPbPb( const Double_t *px, const Double_t *dummy); |
| | 58 | static Double_t PtJpsiBPbPb( const Double_t *px, const Double_t *dummy); |
| | 59 | static Double_t YJpsiPbPb(const Double_t *py, const Double_t *dummy); |
| | 60 | static Double_t YJpsiCDFscaled(const Double_t *py, const Double_t *dummy); |
| | 61 | static Double_t YJpsiCDFscaledPP( const Double_t *px, const Double_t *dummy); |
| | 62 | static Double_t YJpsiCDFscaledPP10( const Double_t *px, const Double_t *dummy); |
| | 63 | static Double_t YJpsiCDFscaledold(const Double_t *py, const Double_t *dummy); |
| | 64 | static Double_t YJpsiCDFscaledPPold( const Double_t *px, const Double_t *dummy); |
| | 65 | static Double_t PtJpsiPP( const Double_t *px, const Double_t *dummy); |
| | 66 | static Double_t YJpsiPP(const Double_t *py, const Double_t *dummy); |
| | 67 | static Double_t YJpsiBPbPb(const Double_t *py, const Double_t *dummy); |
| | 68 | static Int_t IpJpsi(TRandom *ran); |
| | 69 | static Int_t IpJpsiFamily(TRandom *ran); |
| | 70 | static Int_t IpPsiP(TRandom *ran); |
| | 71 | static Double_t PtJpsiFlat( const Double_t *px, const Double_t *dummy ); |
| | 72 | static Double_t YJpsiFlat(const Double_t *py, const Double_t *dummy); |
| | 73 | |
| | 74 | // Upsilon |
| | 75 | static Double_t PtUpsilon( const Double_t *px, const Double_t *dummy ); |
| | 76 | static Double_t PtUpsilonCDFscaled( const Double_t *px, const Double_t *dummy ); |
| | 77 | static Double_t PtUpsilonCDFscaledPP( const Double_t *px, const Double_t *dummy ); |
| | 78 | static Double_t YUpsilon(const Double_t *py, const Double_t *dummy); |
| | 79 | static Double_t YUpsilonCDFscaled(const Double_t *py, const Double_t *dummy); |
| | 80 | static Double_t YUpsilonCDFscaledPP( const Double_t *px, const Double_t *dummy ); |
| | 81 | static Double_t PtUpsilonPbPb( const Double_t *px, const Double_t *dummy ); |
| | 82 | static Double_t YUpsilonPbPb(const Double_t *py, const Double_t *dummy); |
| | 83 | static Double_t PtUpsilonPP( const Double_t *px, const Double_t *dummy ); |
| | 84 | static Double_t YUpsilonPP(const Double_t *py, const Double_t *dummy); |
| | 85 | static Int_t IpUpsilon(TRandom *ran); |
| | 86 | static Int_t IpUpsilonFamily(TRandom *ran); |
| | 87 | static Int_t IpUpsilonP(TRandom *ran); |
| | 88 | static Int_t IpUpsilonPP(TRandom *ran); |
| | 89 | static Double_t PtUpsilonFlat( const Double_t *px, const Double_t *dummy ); |
| | 90 | static Double_t YUpsilonFlat(const Double_t *py, const Double_t *dummy); |
| | 91 | // |
| | 92 | // Charm |
| | 93 | static Double_t PtCharm( const Double_t *px, const Double_t *dummy ); |
| | 94 | static Double_t PtCharmCentral( const Double_t *px, const Double_t *dummy ); |
| | 95 | static Double_t PtCharmF0M0S0PP( const Double_t *px, const Double_t *dummy ); |
| | 96 | static Double_t PtCharmF1M0S0PP( const Double_t *px, const Double_t *dummy ); |
| | 97 | static Double_t PtCharmF2M0S0PP( const Double_t *px, const Double_t *dummy ); |
| | 98 | static Double_t PtCharmF0M1S0PP( const Double_t *px, const Double_t *dummy ); |
| | 99 | static Double_t PtCharmF0M2S0PP( const Double_t *px, const Double_t *dummy ); |
| | 100 | static Double_t PtCharmF0M0S1PP( const Double_t *px, const Double_t *dummy ); |
| | 101 | static Double_t PtCharmF0M0S2PP( const Double_t *px, const Double_t *dummy ); |
| | 102 | static Double_t PtCharmF0M0S3PP( const Double_t *px, const Double_t *dummy ); |
| | 103 | static Double_t PtCharmF0M0S4PP( const Double_t *px, const Double_t *dummy ); |
| | 104 | static Double_t PtCharmF0M0S5PP( const Double_t *px, const Double_t *dummy ); |
| | 105 | static Double_t PtCharmF0M0S6PP( const Double_t *px, const Double_t *dummy ); |
| | 106 | static Double_t YCharm(const Double_t *py, const Double_t *dummy); |
| | 107 | static Double_t YCharmF0M0S0PP(const Double_t *py, const Double_t *dummy); |
| | 108 | static Double_t YCharmF1M0S0PP(const Double_t *py, const Double_t *dummy); |
| | 109 | static Double_t YCharmF2M0S0PP(const Double_t *py, const Double_t *dummy); |
| | 110 | static Double_t YCharmF0M1S0PP(const Double_t *py, const Double_t *dummy); |
| | 111 | static Double_t YCharmF0M2S0PP(const Double_t *py, const Double_t *dummy); |
| | 112 | static Double_t YCharmF0M0S1PP(const Double_t *py, const Double_t *dummy); |
| | 113 | static Double_t YCharmF0M0S2PP(const Double_t *py, const Double_t *dummy); |
| | 114 | static Double_t YCharmF0M0S3PP(const Double_t *py, const Double_t *dummy); |
| | 115 | static Double_t YCharmF0M0S4PP(const Double_t *py, const Double_t *dummy); |
| | 116 | static Double_t YCharmF0M0S5PP(const Double_t *py, const Double_t *dummy); |
| | 117 | static Double_t YCharmF0M0S6PP(const Double_t *py, const Double_t *dummy); |
| | 118 | static Int_t IpCharm(TRandom *ran); |
| | 119 | // |
| | 120 | // Beauty |
| | 121 | static Double_t PtBeauty( const Double_t *px, const Double_t *dummy ); |
| | 122 | static Double_t PtBeautyF0M0S0PP( const Double_t *px, const Double_t *dummy ); |
| | 123 | static Double_t PtBeautyF1M0S0PP( const Double_t *px, const Double_t *dummy ); |
| | 124 | static Double_t PtBeautyF2M0S0PP( const Double_t *px, const Double_t *dummy ); |
| | 125 | static Double_t PtBeautyF0M1S0PP( const Double_t *px, const Double_t *dummy ); |
| | 126 | static Double_t PtBeautyF0M2S0PP( const Double_t *px, const Double_t *dummy ); |
| | 127 | static Double_t PtBeautyF0M0S1PP( const Double_t *px, const Double_t *dummy ); |
| | 128 | static Double_t PtBeautyF0M0S2PP( const Double_t *px, const Double_t *dummy ); |
| | 129 | static Double_t PtBeautyF0M0S3PP( const Double_t *px, const Double_t *dummy ); |
| | 130 | static Double_t PtBeautyF0M0S4PP( const Double_t *px, const Double_t *dummy ); |
| | 131 | static Double_t PtBeautyF0M0S5PP( const Double_t *px, const Double_t *dummy ); |
| | 132 | static Double_t PtBeautyF0M0S6PP( const Double_t *px, const Double_t *dummy ); |
| | 133 | static Double_t YBeauty(const Double_t *py, const Double_t *dummy); |
| | 134 | static Double_t YBeautyF0M0S0PP(const Double_t *py, const Double_t *dummy); |
| | 135 | static Double_t YBeautyF1M0S0PP(const Double_t *py, const Double_t *dummy); |
| | 136 | static Double_t YBeautyF2M0S0PP(const Double_t *py, const Double_t *dummy); |
| | 137 | static Double_t YBeautyF0M1S0PP(const Double_t *py, const Double_t *dummy); |
| | 138 | static Double_t YBeautyF0M2S0PP(const Double_t *py, const Double_t *dummy); |
| | 139 | static Double_t YBeautyF0M0S1PP(const Double_t *py, const Double_t *dummy); |
| | 140 | static Double_t YBeautyF0M0S2PP(const Double_t *py, const Double_t *dummy); |
| | 141 | static Double_t YBeautyF0M0S3PP(const Double_t *py, const Double_t *dummy); |
| | 142 | static Double_t YBeautyF0M0S4PP(const Double_t *py, const Double_t *dummy); |
| | 143 | static Double_t YBeautyF0M0S5PP(const Double_t *py, const Double_t *dummy); |
| | 144 | static Double_t YBeautyF0M0S6PP(const Double_t *py, const Double_t *dummy); |
| | 145 | static Double_t PtBeautyCentral( const Double_t *px, const Double_t *dummy ); |
| | 146 | static Int_t IpBeauty(TRandom *ran); |
| | 147 | // |
| | 148 | |
| | 149 | // Chi 1c 2c |
| | 150 | static Double_t PtChic0( const Double_t *px, const Double_t *dummy); |
| | 151 | static Double_t YChic0(const Double_t *py, const Double_t *dummy); |
| | 152 | static Int_t IpChic0(TRandom *ran); |
| | 153 | |
| | 154 | static Double_t PtChic1( const Double_t *px, const Double_t *dummy); |
| | 155 | static Double_t YChic1(const Double_t *py, const Double_t *dummy); |
| | 156 | static Int_t IpChic1(TRandom *ran); |
| | 157 | |
| | 158 | static Double_t PtChic2( const Double_t *px, const Double_t *dummy); |
| | 159 | static Double_t YChic2(const Double_t *py, const Double_t *dummy); |
| | 160 | static Int_t IpChic2(TRandom *ran); |
| | 161 | |
| | 162 | static Double_t PtChic( const Double_t *px, const Double_t *dummy); |
| | 163 | static Double_t YChic(const Double_t *py, const Double_t *dummy); |
| | 164 | static Int_t IpChic(TRandom *ran); |
| | 165 | |
| | 166 | // |
| | 167 | |
| | 168 | |
| | 169 | static Float_t Interpolate(Float_t x, Float_t* y, Float_t x0, |
| | 170 | Float_t dx, |
| | 171 | Int_t n, Int_t no); |
| | 172 | |
| | 173 | ClassDef(AliGenMUONlib,0) // Library providing y and pT parameterisations |
| | 174 | }; |
| | 175 | #endif |
| | 176 | |
| | 177 | |
| | 178 | |
| | 179 | |
| | 180 | |
| | 181 | |
| | 182 | |
git://git.uio.no / u / mrichter / AliRoot.git / blame_incremental