UrQMD generate zero kinematic energy secondries and induce wrong energy cut counting
In nucleus interaction, UrQMD will put all protons and neutrons of target nucleus into secondaries stack. Many of them do not actually participate in interaction and have zero kinematic energy.
Those zero kinematic energy particles will be removed in particle cut, but their total energy (mass) will be countted into total energy cut.
An example log:
[corsika:debug (Cascade.inl:316)] collide
[corsika:info (InteractionHistogram.inl:32)] pM=0.93957, tM=13.549178660449982, pid=n0, Elab=10.430479551095273, sqrtS=21.61277159601732, pdg=2112 a=6 z=0
[corsika:debug (UrQMD.inl:320)] 0 p+ (0 0 0) eV
[corsika:debug (UrQMD.inl:320)] 1 p+ (0 0 0) eV
[corsika:debug (UrQMD.inl:320)] 2 p+ (0 0 0) eV
[corsika:debug (UrQMD.inl:320)] 3 p+ (0 0 0) eV
[corsika:debug (UrQMD.inl:320)] 4 p+ (0 0 0) eV
[corsika:debug (UrQMD.inl:320)] 5 p+ (0 0 0) eV
[corsika:debug (UrQMD.inl:320)] 6 p+ (0 0 0) eV
[corsika:debug (UrQMD.inl:320)] 7 n0 (0 0 0) eV
[corsika:debug (UrQMD.inl:320)] 8 n0 (0 0 0) eV
[corsika:debug (UrQMD.inl:320)] 9 n0 (0 0 0) eV
[corsika:debug (UrQMD.inl:320)] 10 n0 (0 0 0) eV
[corsika:debug (UrQMD.inl:320)] 11 n0 (0 0 0) eV
[corsika:debug (UrQMD.inl:320)] 12 n0 (0 0 0) eV
[corsika:debug (UrQMD.inl:320)] 13 n0 (-2.54482e+08 1.3539e+08 -7.69883e+08) eV
[corsika:debug (UrQMD.inl:320)] 14 Sigma0 (-9.11716e+07 -1.8696e+08 -8.31646e+09) eV
[corsika:debug (UrQMD.inl:320)] 15 K+ (1.98633e+08 5.00449e+07 -4.75644e+08) eV
[corsika:debug (UrQMD.inl:320)] 16 pi- (-3.47611e+08 -4.11628e+07 -9.18623e+08) eV
[corsika:debug (UrQMD.inl:325)] UrQMD generated 17 secondaries!
[corsika:debug (Cascade.inl:259)] ===========CHECK SECONDARIES============
[corsika:debug (Cascade.inl:260)] Projectile check: n0 (2112), E_kin= 9.490909551095273 GeV, E=10.430479551095273 GeV:
[corsika:debug (Cascade.inl:269)] Secondaries check: p+ (2212), E_kin= 0 GeV, E=0.93827 GeV:
[corsika:debug (Cascade.inl:269)] Secondaries check: p+ (2212), E_kin= 0 GeV, E=0.93827 GeV:
[corsika:debug (Cascade.inl:269)] Secondaries check: p+ (2212), E_kin= 0 GeV, E=0.93827 GeV:
[corsika:debug (Cascade.inl:269)] Secondaries check: p+ (2212), E_kin= 0 GeV, E=0.93827 GeV:
[corsika:debug (Cascade.inl:269)] Secondaries check: p+ (2212), E_kin= 0 GeV, E=0.93827 GeV:
[corsika:debug (Cascade.inl:269)] Secondaries check: p+ (2212), E_kin= 0 GeV, E=0.93827 GeV:
[corsika:debug (Cascade.inl:269)] Secondaries check: p+ (2212), E_kin= 0 GeV, E=0.93827 GeV:
[corsika:debug (Cascade.inl:269)] Secondaries check: n0 (2112), E_kin= 0 GeV, E=0.93957 GeV:
[corsika:debug (Cascade.inl:269)] Secondaries check: n0 (2112), E_kin= 0 GeV, E=0.93957 GeV:
[corsika:debug (Cascade.inl:269)] Secondaries check: n0 (2112), E_kin= 0 GeV, E=0.93957 GeV:
[corsika:debug (Cascade.inl:269)] Secondaries check: n0 (2112), E_kin= 0 GeV, E=0.93957 GeV:
[corsika:debug (Cascade.inl:269)] Secondaries check: n0 (2112), E_kin= 0 GeV, E=0.93957 GeV:
[corsika:debug (Cascade.inl:269)] Secondaries check: n0 (2112), E_kin= 0 GeV, E=0.93957 GeV:
[corsika:debug (Cascade.inl:269)] Secondaries check: n0 (2112), E_kin= 0.3088703816162229 GeV, E=1.2484403816162228 GeV:
[corsika:debug (Cascade.inl:269)] Secondaries check: Sigma0 (3212), E_kin= 7.211473603523096 GeV, E=8.404113603523095 GeV:
[corsika:debug (Cascade.inl:269)] Secondaries check: K+ (321), E_kin= 0.2218038677297591 GeV, E=0.7154838677297591 GeV:
[corsika:debug (Cascade.inl:269)] Secondaries check: pi- (-211), E_kin= 0.8533429596920387 GeV, E=0.9929129596920386 GeV:
[corsika:debug (Cascade.inl:273)] ===========CHECK SECONDARIES============
[corsika:debug (ParticleCut.inl:127)] ParticleCut: checking p+ (2212), E_kin= 0 GeV, EcutTot=93.61082464847107 GeV, E=0.93827 GeV, m=0.93827 GeV
[corsika:debug (ParticleCut.inl:133)] p=[[[[particle: i=15, PID=p+, Ekin=0GeV][node=0x56490ad69a40]]][i_parent=0, [evt: hasParent=true, projIndex=14, Nsec=17]]]
[corsika:debug (ParticleCut.inl:147)] removing low en. particle...
[corsika:debug (ParticleCut.inl:127)] ParticleCut: checking p+ (2212), E_kin= 0 GeV, EcutTot=94.54909464847107 GeV, E=0.93827 GeV, m=0.93827 GeV
[corsika:debug (ParticleCut.inl:133)] p=[[[[particle: i=16, PID=p+, Ekin=0GeV][node=0x56490ad69a40]]][i_parent=1, [evt: hasParent=true, projIndex=14, Nsec=17]]]
[corsika:debug (ParticleCut.inl:147)] removing low en. particle...
[corsika:debug (ParticleCut.inl:127)] ParticleCut: checking p+ (2212), E_kin= 0 GeV, EcutTot=95.48736464847107 GeV, E=0.93827 GeV, m=0.93827 GeV
[corsika:debug (ParticleCut.inl:133)] p=[[[[particle: i=17, PID=p+, Ekin=0GeV][node=0x56490ad69a40]]][i_parent=2, [evt: hasParent=true, projIndex=14, Nsec=17]]]
[corsika:debug (ParticleCut.inl:147)] removing low en. particle...
[corsika:debug (ParticleCut.inl:127)] ParticleCut: checking p+ (2212), E_kin= 0 GeV, EcutTot=96.42563464847107 GeV, E=0.93827 GeV, m=0.93827 GeV
[corsika:debug (ParticleCut.inl:133)] p=[[[[particle: i=18, PID=p+, Ekin=0GeV][node=0x56490ad69a40]]][i_parent=3, [evt: hasParent=true, projIndex=14, Nsec=17]]]
[corsika:debug (ParticleCut.inl:147)] removing low en. particle...
[corsika:debug (ParticleCut.inl:127)] ParticleCut: checking p+ (2212), E_kin= 0 GeV, EcutTot=97.36390464847106 GeV, E=0.93827 GeV, m=0.93827 GeV
[corsika:debug (ParticleCut.inl:133)] p=[[[[particle: i=19, PID=p+, Ekin=0GeV][node=0x56490ad69a40]]][i_parent=4, [evt: hasParent=true, projIndex=14, Nsec=17]]]
[corsika:debug (ParticleCut.inl:147)] removing low en. particle...
[corsika:debug (ParticleCut.inl:127)] ParticleCut: checking p+ (2212), E_kin= 0 GeV, EcutTot=98.30217464847107 GeV, E=0.93827 GeV, m=0.93827 GeV
[corsika:debug (ParticleCut.inl:133)] p=[[[[particle: i=20, PID=p+, Ekin=0GeV][node=0x56490ad69a40]]][i_parent=5, [evt: hasParent=true, projIndex=14, Nsec=17]]]
[corsika:debug (ParticleCut.inl:147)] removing low en. particle...
[corsika:debug (ParticleCut.inl:127)] ParticleCut: checking p+ (2212), E_kin= 0 GeV, EcutTot=99.24044464847107 GeV, E=0.93827 GeV, m=0.93827 GeV
[corsika:debug (ParticleCut.inl:133)] p=[[[[particle: i=21, PID=p+, Ekin=0GeV][node=0x56490ad69a40]]][i_parent=6, [evt: hasParent=true, projIndex=14, Nsec=17]]]
[corsika:debug (ParticleCut.inl:147)] removing low en. particle...
[corsika:debug (ParticleCut.inl:127)] ParticleCut: checking n0 (2112), E_kin= 0 GeV, EcutTot=100.17871464847107 GeV, E=0.93957 GeV, m=0.93957 GeV
[corsika:debug (ParticleCut.inl:133)] p=[[[[particle: i=22, PID=n0, Ekin=0GeV][node=0x56490ad69a40]]][i_parent=7, [evt: hasParent=true, projIndex=14, Nsec=17]]]
[corsika:debug (ParticleCut.inl:147)] removing low en. particle...
[corsika:debug (ParticleCut.inl:127)] ParticleCut: checking n0 (2112), E_kin= 0 GeV, EcutTot=101.11828464847108 GeV, E=0.93957 GeV, m=0.93957 GeV
[corsika:debug (ParticleCut.inl:133)] p=[[[[particle: i=23, PID=n0, Ekin=0GeV][node=0x56490ad69a40]]][i_parent=8, [evt: hasParent=true, projIndex=14, Nsec=17]]]
[corsika:debug (ParticleCut.inl:147)] removing low en. particle...
[corsika:debug (ParticleCut.inl:127)] ParticleCut: checking n0 (2112), E_kin= 0 GeV, EcutTot=102.05785464847106 GeV, E=0.93957 GeV, m=0.93957 GeV
[corsika:debug (ParticleCut.inl:133)] p=[[[[particle: i=24, PID=n0, Ekin=0GeV][node=0x56490ad69a40]]][i_parent=9, [evt: hasParent=true, projIndex=14, Nsec=17]]]
[corsika:debug (ParticleCut.inl:147)] removing low en. particle...
[corsika:debug (ParticleCut.inl:127)] ParticleCut: checking n0 (2112), E_kin= 0 GeV, EcutTot=102.99742464847107 GeV, E=0.93957 GeV, m=0.93957 GeV
[corsika:debug (ParticleCut.inl:133)] p=[[[[particle: i=25, PID=n0, Ekin=0GeV][node=0x56490ad69a40]]][i_parent=10, [evt: hasParent=true, projIndex=14, Nsec=17]]]
[corsika:debug (ParticleCut.inl:147)] removing low en. particle...
[corsika:debug (ParticleCut.inl:127)] ParticleCut: checking n0 (2112), E_kin= 0 GeV, EcutTot=103.93699464847107 GeV, E=0.93957 GeV, m=0.93957 GeV
[corsika:debug (ParticleCut.inl:133)] p=[[[[particle: i=26, PID=n0, Ekin=0GeV][node=0x56490ad69a40]]][i_parent=11, [evt: hasParent=true, projIndex=14, Nsec=17]]]
[corsika:debug (ParticleCut.inl:147)] removing low en. particle...
[corsika:debug (ParticleCut.inl:127)] ParticleCut: checking n0 (2112), E_kin= 0 GeV, EcutTot=104.87656464847107 GeV, E=0.93957 GeV, m=0.93957 GeV
[corsika:debug (ParticleCut.inl:133)] p=[[[[particle: i=27, PID=n0, Ekin=0GeV][node=0x56490ad69a40]]][i_parent=12, [evt: hasParent=true, projIndex=14, Nsec=17]]]
[corsika:debug (ParticleCut.inl:147)] removing low en. particle...
[corsika:debug (ParticleCut.inl:127)] ParticleCut: checking n0 (2112), E_kin= 0.3088703816162229 GeV, EcutTot=105.81613464847106 GeV, E=1.2484403816162228 GeV, m=0.93957 GeV
[corsika:debug (ParticleCut.inl:133)] p=[[[[particle: i=28, PID=n0, Ekin=0.3088703816162229GeV][node=0x56490ad69a40]]][i_parent=13, [evt: hasParent=true, projIndex=14, Nsec=17]]]
[corsika:debug (ParticleCut.inl:147)] removing low en. particle...
[corsika:debug (ParticleCut.inl:127)] ParticleCut: checking Sigma0 (3212), E_kin= 7.211473603523096 GeV, EcutTot=107.0645750300873 GeV, E=8.404113603523095 GeV, m=1.19264 GeV
[corsika:debug (ParticleCut.inl:133)] p=[[[[particle: i=29, PID=Sigma0, Ekin=7.211473603523096GeV][node=0x56490ad69a40]]][i_parent=14, [evt: hasParent=true, projIndex=14, Nsec=17]]]
[corsika:debug (ParticleCut.inl:127)] ParticleCut: checking K+ (321), E_kin= 0.2218038677297591 GeV, EcutTot=107.0645750300873 GeV, E=0.7154838677297591 GeV, m=0.49368 GeV
[corsika:debug (ParticleCut.inl:133)] p=[[[[particle: i=30, PID=K+, Ekin=0.2218038677297591GeV][node=0x56490ad69a40]]][i_parent=15, [evt: hasParent=true, projIndex=14, Nsec=17]]]
[corsika:debug (ParticleCut.inl:147)] removing low en. particle...
[corsika:debug (ParticleCut.inl:127)] ParticleCut: checking pi- (-211), E_kin= 0.8533429596920387 GeV, EcutTot=107.78005889781704 GeV, E=0.9929129596920386 GeV, m=0.13957 GeV
[corsika:debug (ParticleCut.inl:133)] p=[[[[particle: i=31, PID=pi-, Ekin=0.8533429596920387GeV][node=0x56490ad69a40]]][i_parent=16, [evt: hasParent=true, projIndex=14, Nsec=17]]]
[corsika:debug (ParticleCut.inl:147)] removing low en. particle...
[corsika:debug (ParticleCut.inl:169)] Event cut: 15.162147209038022 GeV