IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 8fc60ea0 authored by Maximilian Reininghaus's avatar Maximilian Reininghaus :vulcan: Committed by Maximilian Reininghaus
Browse files

fixed error in PDG/UrQMD conversion

parent 954acac1
No related branches found
No related tags found
No related merge requests found
...@@ -176,9 +176,6 @@ corsika::process::EProcessReturn UrQMD::DoInteraction(SetupProjectile& vProjecti ...@@ -176,9 +176,6 @@ corsika::process::EProcessReturn UrQMD::DoInteraction(SetupProjectile& vProjecti
std::cout << "UrQMD generated " << sys_.npart << " secondaries!" << std::endl; std::cout << "UrQMD generated " << sys_.npart << " secondaries!" << std::endl;
if (sys_.npart > 0) // delete only in case of inelastic collision, otherwise keep
vProjectile.Delete();
return process::EProcessReturn::eOk; return process::EProcessReturn::eOk;
} }
...@@ -206,13 +203,15 @@ corsika::particles::Code corsika::process::UrQMD::ConvertFromUrQMD(int vItyp, in ...@@ -206,13 +203,15 @@ corsika::particles::Code corsika::process::UrQMD::ConvertFromUrQMD(int vItyp, in
std::pair<int, int> corsika::process::UrQMD::ConvertToUrQMD( std::pair<int, int> corsika::process::UrQMD::ConvertToUrQMD(
corsika::particles::Code code) { corsika::particles::Code code) {
static const std::map<int, std::pair<int, int>> mapPDGToUrQMD{ static const std::map<int, std::pair<int, int>> mapPDGToUrQMD{
// data from github.com/afedynitch/ParticleDataTool // data mostly from github.com/afedynitch/ParticleDataTool
{22, {100, 0}}, // gamma {22, {100, 0}}, // gamma
{111, {101, 0}}, // pi0 {111, {101, 0}}, // pi0
{211, {101, 2}}, // pi+ {211, {101, 2}}, // pi+
{-211, {101, -2}}, // pi- {-211, {101, -2}}, // pi-
{321, {106, 1}}, // K+ {321, {106, 1}}, // K+
{-321, {106, -1}}, // K- {-321, {-106, -1}}, // K-
{311, {106, -1}}, // K0
{-311, {-106, 1}}, // K0bar
{2212, {1, 1}}, // p {2212, {1, 1}}, // p
{2112, {1, -1}}, // n {2112, {1, -1}}, // n
{221, {102, 0}}, // eta {221, {102, 0}}, // eta
......
...@@ -216,7 +216,8 @@ cl check if we found the correct values in IDTAB ...@@ -216,7 +216,8 @@ cl check if we found the correct values in IDTAB
pdgid = idtab(3,next) pdgid = idtab(3,next)
endif endif
else else
call error ('pdgid','Error in tablelookup',dble(next),3) ! M.R 2019-04-27, error to be handled in CORSIKA instead of here
! call error ('pdgid','Error in tablelookup',dble(next),3)
pdgid = 0 pdgid = 0
endif endif
......
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