IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 34f29d46 authored by Maximilian Reininghaus's avatar Maximilian Reininghaus :vulcan:
Browse files

copy particles from eventMain to C8

parent 444d19b5
No related branches found
No related tags found
1 merge request!427Resolve "upgrade pythia to version 8.3xx"
...@@ -43,7 +43,7 @@ namespace corsika::pythia8 { ...@@ -43,7 +43,7 @@ namespace corsika::pythia8 {
pythiaMain_.readString("321:mayDecay = on"); pythiaMain_.readString("321:mayDecay = on");
pythiaMain_.readString("130:mayDecay = on"); pythiaMain_.readString("130:mayDecay = on");
// Redure statistics printout to relevant ones. // Reduce statistics printout to relevant ones.
pythiaMain_.readString("Stat:showProcessLevel = off"); pythiaMain_.readString("Stat:showProcessLevel = off");
pythiaMain_.readString("Stat:showPartonLevel = off"); pythiaMain_.readString("Stat:showPartonLevel = off");
...@@ -349,7 +349,7 @@ namespace corsika::pythia8 { ...@@ -349,7 +349,7 @@ namespace corsika::pythia8 {
sizeOld = eventMain.size(); sizeOld = eventMain.size();
for (int iSub = 3; iSub < eventColl.size(); ++iSub) { for (int iSub = 3; iSub < eventColl.size(); ++iSub) {
if (!eventColl[iSub].isFinal()) continue; if (!eventColl[iSub].isFinal()) continue;
int iNew = eventMain.append(eventColl[iSub]); int const iNew = eventMain.append(eventColl[iSub]);
eventMain[iNew].mothers(iNuc, iProj); eventMain[iNew].mothers(iNuc, iProj);
eventMain[iNew].vProdAdd(vNow); eventMain[iNew].vProdAdd(vNow);
} }
...@@ -377,9 +377,9 @@ namespace corsika::pythia8 { ...@@ -377,9 +377,9 @@ namespace corsika::pythia8 {
MomentumVector Plab_final{labFrameBoost.getOriginalCS()}; MomentumVector Plab_final{labFrameBoost.getOriginalCS()};
auto Elab_final = HEPEnergyType::zero(); auto Elab_final = HEPEnergyType::zero();
for (Pythia8::Particle const& p8p: eventColl) { for (Pythia8::Particle const& p8p : eventMain) {
// skip particles that have decayed / are initial particles in pythia's event record // skip particles that have decayed / are initial particles in pythia's event record
if (!p8p.isFinal()) continue; if (!p8p.isFinal()) continue;
try { try {
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment