From 9fd8732dcc6d7928407fbf113d9c19eb46053f5b Mon Sep 17 00:00:00 2001
From: Felix Riehn <friehn@lip.pt>
Date: Fri, 16 Feb 2024 10:47:08 +0100
Subject: [PATCH] fix frame
---
corsika/detail/modules/pythia8/NeutrinoInteraction.inl | 9 +++++----
1 file changed, 5 insertions(+), 4 deletions(-)
diff --git a/corsika/detail/modules/pythia8/NeutrinoInteraction.inl b/corsika/detail/modules/pythia8/NeutrinoInteraction.inl
index 183d360c2..9bb135ad1 100644
--- a/corsika/detail/modules/pythia8/NeutrinoInteraction.inl
+++ b/corsika/detail/modules/pythia8/NeutrinoInteraction.inl
@@ -116,16 +116,17 @@ namespace corsika::pythia8 {
double Q2min = minQ2_ / 1_GeV;
// Set up incoming beams, for frame with unequal beam energies.
+ // projectile is along -z
pythiaMain_.readString("Beams:frameType = 2");
// BeamA = proton.
- pythiaMain_.readString("Beams:idB = 2212");
- pythiaMain_.settings.parm("Beams:eB", eProton);
+ pythiaMain_.readString("Beams:idA = 2212");
+ pythiaMain_.settings.parm("Beams:eA", eProton);
// BeamB = electron.
int const idNow = static_cast<int>(get_PDG(projectileId));
- pythiaMain_.readString("Beams:idA = 12");
+ pythiaMain_.readString("Beams:idB = 12");
// pythiaMain_.settings.parm("Beams:idA", idNow);
- pythiaMain_.settings.parm("Beams:eA", eElectron);
+ pythiaMain_.settings.parm("Beams:eB", eElectron);
// Set up DIS process within some phase space.
// Neutral current (with gamma/Z interference).
--
GitLab