TuneA settings removed from kPyJets
[u/mrichter/AliRoot.git] / EVGEN / AliGenLcLib.cxx
CommitLineData
feb24caa 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
16
17// Library class for particle pt and y distributions used for
18// LambdaC simulations.
19// To be used with AliGenParam.
20//
21// Author: Annalisa Mastroserio <Annalisa.Mastroserio@cern.ch>
22//
23
24#include <TPDGCode.h>
25#include <TMath.h>
26#include <TRandom.h>
27#include <TString.h>
28
29#include "AliGenLcLib.h"
30#include "AliLog.h"
31
32ClassImp(AliGenLcLib)
33
34
35//---------------------------------------------
36// LambdaC
37//---------------------------------------------
38 Int_t AliGenLcLib::IpLcPlus(TRandom *)
39{
40 //PDG code
41 return 4122;
42}
43
44Int_t AliGenLcLib::IpLcMinus(TRandom *)
45{
46 //PDG code
47 return -4122;
48}
49
50Double_t AliGenLcLib::PtLcFlat( const Double_t *, const Double_t *)
51{
52 // FLAT pt-distribution
53 return 1;
54}
55
56Double_t AliGenLcLib::PtLcExp( const Double_t *x, const Double_t *)
57{
0f8eca74 58 // pt-distribution
59 //return x[0]*TMath::Exp(-x[0]/0.16); // distribution used in LHC11f1 for the anchor runs : 139441, 139510, 139511, 130513, 130514, 130517.
60 return TMath::GammaDist(x[0],2,0,1.7); //distribution as in LHC11a10a of the prompt Lc whose daughters are in |eta|<0.9. Used for Lb as well.
feb24caa 61}
62
0f8eca74 63Double_t AliGenLcLib::YLcFlat(const Double_t *,const Double_t *)
feb24caa 64{
65 //LambdaC y-distribution
0f8eca74 66 return 1;
feb24caa 67}
68
69
70
71typedef Double_t (*GenFunc) (const Double_t*, const Double_t*);
72typedef Int_t (*GenFuncIp) (TRandom *);
73
74GenFunc AliGenLcLib::GetPt(Int_t iPID, const char * sForm) const
75{
76 // Return pointer to Pt parameterisation
77 printf("PID: %i, form: %s \n",iPID,sForm);
78 TString type(sForm);
79 GenFunc func;
80
81 switch(iPID) {
82
83 case kLcPlus:
84 if (type=="FLAT") {func=PtLcFlat; break;}
85 else if(type=="EXP") {func=PtLcExp; break;}
86 else {
87 AliFatal(Form("Unknown Pt distribution form: %s",sForm)); func=0;
88 }
89
90 case kLcMinus:
91 if (type=="FLAT") {func=PtLcFlat; break;}
92 else if(type=="EXP") {func=PtLcExp; break;}
93 else {
94 AliFatal(Form("Unknown Pt distribution form: %s",sForm)); func=0;
95 }
96
97 default : AliFatal(Form("Unknown particle type: %i",iPID)); func=0;
98 }//switch
99
100 return func;
101}
102
103GenFunc AliGenLcLib::GetY(Int_t iPID, const char *sForm) const
104{
105 AliDebug(1,Form("PID: %i, form: %s",iPID,sForm));
106 GenFunc func;
107 switch (iPID) {
108
109 case kLcPlus: func=YLcFlat; break;
110 case kLcMinus: func=YLcFlat; break;
111
112 default : AliFatal(Form("Unknown particle type: %i",iPID)); func=0; break;
113
114 }//switch
115 return func;
116}
117
118GenFuncIp AliGenLcLib::GetIp(Int_t iPID, const char *sForm) const
119{
120 // Return pointer to particle type parameterisation
121 AliDebug(1,Form("PID: %i, form: %s",iPID,sForm)); //////////
122
123 switch (iPID){
124
125 case kLcPlus: return IpLcPlus;
126 case kLcMinus: return IpLcMinus;
127
128 default : AliFatal(Form("Unknown particle type: %i",iPID)); return 0;
129 }
130}