55#include "FairPrimaryGenerator.h"
66#include "Generators/GeneratorPythia8.h"
77#include "TRandom3.h"
8+ #include "TF1.h"
89#include "TParticlePDG.h"
910#include "TDatabasePDG.h"
1011
@@ -48,11 +49,11 @@ public:
4849
4950 fSpectrumOm = std ::make_unique < TF1 > ("fSpectrumOm" , this , & GeneratorPythia8ExtraStrangeness ::boltzPlusPower , 0. , genMaxPt , 5 , "GeneratorPythia8ExtraStrangeness" , "boltzPlusPower" );
5051
51- fSpectrumXi -> FixParameter (0 , 1.67245e+00 );
52- fSpectrumXi -> FixParameter (1 , 5.18174e-01 );
53- fSpectrumXi -> FixParameter (2 , 1.73747e+01 );
54- fSpectrumXi -> FixParameter (3 , -2.56681e+00 );
55- fSpectrumXi -> FixParameter (4 , 1.87513e-04 );
52+ fSpectrumOm -> FixParameter (0 , 1.67245e+00 );
53+ fSpectrumOm -> FixParameter (1 , 5.18174e-01 );
54+ fSpectrumOm -> FixParameter (2 , 1.73747e+01 );
55+ fSpectrumOm -> FixParameter (3 , -2.56681e+00 );
56+ fSpectrumOm -> FixParameter (4 , 1.87513e-04 );
5657 }
5758
5859 Double_t y2eta (Double_t pt , Double_t mass , Double_t y ){
@@ -101,7 +102,7 @@ public:
101102 ranGenerator -> SetSeed (0 );
102103
103104 // generate transverse momentum
104- const double gen_pT = fSpectraXi -> GetRandom (genMinPt ,genMaxPt );
105+ const double gen_pT = fSpectrumXi -> GetRandom (genMinPt ,genMaxPt );
105106
106107 //Actually could be something else without loss of generality but okay
107108 const double gen_phi = ranGenerator -> Uniform (0 ,2 * TMath ::Pi ());
@@ -127,7 +128,7 @@ public:
127128 ranGenerator -> SetSeed (0 );
128129
129130 // generate transverse momentum
130- const double gen_pT = fSpectraOm -> GetRandom (genMinPt ,genMaxPt );
131+ const double gen_pT = fSpectrumOm -> GetRandom (genMinPt ,genMaxPt );
131132
132133 //Actually could be something else without loss of generality but okay
133134 const double gen_phi = ranGenerator -> Uniform (0 ,2 * TMath ::Pi ());
0 commit comments