search

CRPropa / CRPropa3

CRPropa is a public astrophysical simulation framework for propagating extraterrestrial ultra-high energy particles. https://crpropa.github.io/CRPropa3/
https://crpropa.desy.de
GNU General Public License v3.0
71 stars 69 forks source link

some photons in secondary neutrino output #159

Closed saikatxdas closed 6 years ago

saikatxdas commented 6 years ago

In the example for secondary neutrinos, I am keeping all other modules same and change the particle type for propagation to iron nuclei. source.add(SourceParticleType(nucleusId(56,26))) Thus the output file for neutrino should store only neutrinos & antineutrinos. But in the output, I obtain few but a certain number of particles with ID 22, i.e. photons. What can be the possible reason? Is ElectronPairProduction( ) causing the trouble? Should I use ObserverPhotonVeto( ) too in case of propagating heavier nuclei? There is no issue with protons.

Can you please look at it?

avvliet commented 6 years ago

This was also due to the missing electrons in m.add(NuclearDecay(electrons, photons, neutrinos)). I've fixed the example. It should be correct now.

saikatxdas commented 6 years ago

Okay. Thank you very much for the response. Yes, now it's working completely fine. Thank you once again.

Saikat