Skip to content

Commit

Permalink
Merge pull request #44264 from hjbossi/PhotonFluxFix
Browse files Browse the repository at this point in the history
Fix to PYTHIA 8 photon flux.
  • Loading branch information
cmsbuild authored Mar 6, 2024
2 parents 4c67638 + bca208b commit 5fe7484
Showing 1 changed file with 4 additions and 4 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -84,16 +84,16 @@ class Nucleus2gamma2 : public Pythia8::PDF {
Nucleus2gamma2(int idBeamIn) : Pythia8::PDF(idBeamIn) {}

// Update the photon flux.
void xfUpdate(int idBeamIn, double x, double) override {
void xfUpdate(int, double x, double) override {
// lead
double radius = 0; // radius in [fm]
double z = 0;
if (idBeamIn == 1000822080) {
if (idBeam == 1000822080) {
radius = 6.636;
z = 82;
}
// oxygen
else if (idBeamIn == 80160) {
else if (idBeam == 80160) {
radius = 3.02;
z = 8;
}
Expand Down Expand Up @@ -413,7 +413,7 @@ bool Pythia8Hadronizer::initializeForInternalPartons() {
}

if (doProtonPhotonFlux) {
photonFlux = make_shared<Nucleus2gamma2>(2212);
photonFlux = make_shared<Nucleus2gamma2>(1000822080);
fMasterGen->setPhotonFluxPtr(photonFlux, nullptr);
}

Expand Down

0 comments on commit 5fe7484

Please sign in to comment.