Skip to content

PWGLF: improve Xi/Om injected spectra #1784

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
sawenzel merged 2 commits into from
Nov 5, 2024
Merged

PWGLF: improve Xi/Om injected spectra #1784

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
a178e8f
PWGLF: Change Xi/Om spectra for injection
ddobrigk Nov 4, 2024
8bd2fdd
Fix bugs
ddobrigk Nov 4, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
43 changes: 41 additions & 2 deletions MC/config/PWGLF/pythia8/generator_pythia8_extraStrangeness.C
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -5,6 +5,7 @@
#include "FairPrimaryGenerator.h"
#include "Generators/GeneratorPythia8.h"
#include "TRandom3.h"
#include "TF1.h"
#include "TParticlePDG.h"
#include "TDatabasePDG.h"

Expand Down Expand Up @@ -37,6 +38,22 @@ public:
xProd=0.; yProd=0.; zProd=0.;

fLVHelper = std::make_unique<TLorentzVector>();

fSpectrumXi = std::make_unique<TF1>("fSpectrumXi", this, &GeneratorPythia8ExtraStrangeness::boltzPlusPower, 0., genMaxPt, 5, "GeneratorPythia8ExtraStrangeness", "boltzPlusPower");

fSpectrumXi->FixParameter(0, 1.32171);
fSpectrumXi->FixParameter(1, 4.84e-1);
fSpectrumXi->FixParameter(2, 111.9);
fSpectrumXi->FixParameter(3, -2.56511e+00);
fSpectrumXi->FixParameter(4, 1.14011e-04);

fSpectrumOm = std::make_unique<TF1>("fSpectrumOm", this, &GeneratorPythia8ExtraStrangeness::boltzPlusPower, 0., genMaxPt, 5, "GeneratorPythia8ExtraStrangeness", "boltzPlusPower");

fSpectrumOm->FixParameter(0, 1.67245e+00);
fSpectrumOm->FixParameter(1, 5.18174e-01);
fSpectrumOm->FixParameter(2, 1.73747e+01);
fSpectrumOm->FixParameter(3, -2.56681e+00);
fSpectrumOm->FixParameter(4, 1.87513e-04);
}

Double_t y2eta(Double_t pt, Double_t mass, Double_t y){
Expand Down Expand Up @@ -85,7 +102,7 @@ public:
ranGenerator->SetSeed(0);

// generate transverse momentum
const double gen_pT = ranGenerator->Uniform(0,5);
const double gen_pT = fSpectrumXi->GetRandom(genMinPt,genMaxPt);

//Actually could be something else without loss of generality but okay
const double gen_phi = ranGenerator->Uniform(0,2*TMath::Pi());
Expand All @@ -111,7 +128,7 @@ public:
ranGenerator->SetSeed(0);

// generate transverse momentum
const double gen_pT = ranGenerator->Uniform(0,5);
const double gen_pT = fSpectrumOm->GetRandom(genMinPt,genMaxPt);

//Actually could be something else without loss of generality but okay
const double gen_phi = ranGenerator->Uniform(0,2*TMath::Pi());
Expand All @@ -129,6 +146,26 @@ public:
set4momentum(fLVHelper->Px(),fLVHelper->Py(),fLVHelper->Pz());
}

Double_t boltzPlusPower(const Double_t *x, const Double_t *p)
{
// a plain parametrization. not meant to be physics worthy.
// adjusted to match preliminary 5 TeV shape.

Double_t pt = x[0];
Double_t mass = p[0];
Double_t mt = TMath::Sqrt(pt * pt + mass * mass);
Double_t T = p[1];
Double_t norm = p[2];

Double_t lowptpart = mt * TMath::Exp(-mt / T);
Double_t highptpart = p[4]*TMath::Power(x[0], p[3]);

Double_t mixup = 1./(1.+TMath::Exp((x[0]-4.5)/.1));

//return pt * norm * (mixup * mt * TMath::Exp(-mt / T) + (1.-mixup)*highptpart) ;
return pt * norm * (mt * TMath::Exp(-mt / T) + (1.-mixup)*highptpart) ;
}

//__________________________________________________________________
Bool_t generateEvent() override {

Expand Down Expand Up @@ -235,6 +272,8 @@ private:
double zProd; /// z-coordinate position production vertex [cm]

std::unique_ptr<TLorentzVector> fLVHelper;
std::unique_ptr<TF1> fSpectrumXi;
std::unique_ptr<TF1> fSpectrumOm;
};

FairGenerator *generator_extraStrangeness()
Expand Down