Skip to content

[PWGDQ] Add pythia8 prompt J/psi generator at forward rapidity #1902

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
lucamicheletti93 merged 2 commits into from
Feb 14, 2025
Merged

[PWGDQ] Add pythia8 prompt J/psi generator at forward rapidity #1902

Show file tree
Hide file tree
Changes from 1 commit
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
34f5956
[PWGDQ] Add pythia8 prompt J/psi generator at forward rapidity
mcoquet642 Feb 11, 2025
0f964a1
Adding test macro
mcoquet642 Feb 12, 2025
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
35 changes: 34 additions & 1 deletion MC/config/PWGDQ/external/generator/generator_pythia8Onia_PromptSignals_gaptriggered.C
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -173,4 +173,37 @@ FairGenerator*
// gen->PrintDebug();

return gen;
}
}

FairGenerator*
GeneratorPromptJpsi_EvtGenFwdy(int triggerGap, double rapidityMin = -4.3, double rapidityMax = -2.3, bool verbose = false)
{
auto gen = new o2::eventgen::GeneratorEvtGen<o2::eventgen::GeneratorPythia8OniaPromptSignalsGapTriggered>();
gen->setTriggerGap(triggerGap);
gen->setRapidityRange(rapidityMin, rapidityMax);
gen->addHadronPDGs(443);
gen->setVerbose(verbose);

TString pathO2table = gSystem->ExpandPathName("${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGDQ/pythia8/decayer/switchOffJpsi.cfg");
gen->readFile(pathO2table.Data());
gen->setConfigMBdecays(pathO2table);
gen->PrintDebug(true);

gen->SetSizePdg(1);
gen->AddPdg(443, 0);

gen->SetForceDecay(kEvtDiMuon);

// set random seed
gen->readString("Random:setSeed on");
uint random_seed;
unsigned long long int random_value = 0;
ifstream urandom("/dev/urandom", ios::in|ios::binary);
urandom.read(reinterpret_cast<char*>(&random_value), sizeof(random_seed));
gen->readString(Form("Random:seed = %llu", random_value % 900000001));

// print debug
// gen->PrintDebug();

return gen;
}
7 changes: 7 additions & 0 deletions MC/config/PWGDQ/ini/Generator_InjectedPromptJpsiFwdy_PythiaOnia_TriggerGap.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,7 @@
### The external generator derives from GeneratorPythia8.
[GeneratorExternal]
fileName=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGDQ/external/generator/generator_pythia8Onia_PromptSignals_gaptriggered.C
funcName=GeneratorPromptJpsi_EvtGenFwdy(5,-4.3,-2.3)

[GeneratorPythia8]
config=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGDQ/pythia8/generator/pythia8_onia_triggerGap.cfg
2 changes: 1 addition & 1 deletion MC/config/common/pythia8/generator/pythia8_hf.cfg
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
### beams
Beams:idA 2212 # proton
Beams:idB 2212 # proton
Beams:eCM 14000. # GeV
Beams:eCM 13600. # GeV

### processes
HardQCD:hardccbar on # scatterings g-g / q-qbar -> c-cbar
Expand Down