From 0a7077f11dee923548a2ec24c829e2937e7f5352 Mon Sep 17 00:00:00 2001 From: Felix Riehn <friehn@lip.pt> Date: Wed, 31 Jul 2024 20:41:32 +0200 Subject: [PATCH] shut up pythia --- corsika/detail/modules/pythia8/Interaction.inl | 17 +++++++++-------- 1 file changed, 9 insertions(+), 8 deletions(-) diff --git a/corsika/detail/modules/pythia8/Interaction.inl b/corsika/detail/modules/pythia8/Interaction.inl index 470cc66b2..a6ad2771e 100644 --- a/corsika/detail/modules/pythia8/Interaction.inl +++ b/corsika/detail/modules/pythia8/Interaction.inl @@ -50,7 +50,7 @@ namespace corsika::pythia8 { pythia_.readString("Beams:allowVariableEnergy = on"); pythia_.readString("Beams:frameType = 1"); // CoM - pythia_.settings.parm("Beams:eCM", 400_TeV/1_GeV); + pythia_.settings.parm("Beams:eCM", 400_TeV / 1_GeV); // needed when regular Pythia (not Angantyr) takes over for pp pythia_.readString("SoftQCD:all = on"); @@ -74,7 +74,7 @@ namespace corsika::pythia8 { pythia_.readString("HadronLevel:Decay = off"); // Reduce printout and relax energy-momentum conservation. - //pythia_.readString("Print:quiet = on"); + pythia_.readString("Print:quiet = on"); pythia_.readString("Check:epTolErr = 0.1"); pythia_.readString("Check:epTolWarn = 0.0001"); pythia_.readString("Check:mTolErr = 0.01"); @@ -218,10 +218,11 @@ namespace corsika::pythia8 { projectileId, corsika::pythia8::Interaction::canInteract(projectileId)); // define system (this is Nucleus-Nucleus!!) - auto const sqrtS2NN = (projectileP4 + targetP4/get_nucleus_A(targetId)).getNormSqr(); + auto const sqrtS2NN = + (projectileP4 + targetP4 / get_nucleus_A(targetId)).getNormSqr(); HEPEnergyType const sqrtSNN = sqrt(sqrtS2NN); - HEPEnergyType const eProjectileLab = - calculate_lab_energy(sqrtS2NN, get_mass(projectileId), get_mass(targetId)/get_nucleus_A(targetId)); + HEPEnergyType const eProjectileLab = calculate_lab_energy( + sqrtS2NN, get_mass(projectileId), get_mass(targetId) / get_nucleus_A(targetId)); if (!isValid(projectileId, targetId, sqrtSNN)) { throw std::runtime_error("invalid target,projectile,energy combination."); @@ -240,11 +241,11 @@ namespace corsika::pythia8 { // auto const pProjLab = proj4MomLab.getSpaceLikeComponents().getComponents(rotCS); // variables to talk to Pythia - double const eCM_pythia = - sqrtSNN / 1_GeV; // CoM energy hadron-Nucleon + double const eCM_pythia = sqrtSNN / 1_GeV; // CoM energy hadron-Nucleon double const idA_pythia = static_cast<double>(get_PDG(projectileId)); double const idB_pythia = static_cast<double>(get_PDG(targetId)); - CORSIKA_LOG_INFO("re-configuring pythia idA={}, idB={}, ecm={} GeV", idA_pythia, idB_pythia, eCM_pythia); + CORSIKA_LOG_INFO("re-configuring pythia idA={}, idB={}, ecm={} GeV", idA_pythia, + idB_pythia, eCM_pythia); // configure event pythia_.setBeamIDs(idA_pythia, idB_pythia); pythia_.setKinematics(eCM_pythia); -- GitLab