Hi @xzy66
I assume that you are referring to the following thread:
where indeed the cross section for photonuclear reaction of 68Zn to 67Cu was reported to be lower. This was the case for FLUKA v4-1. However, it seems that this particular discrepancy could still be there in the current version as I did not find any direct mention in the release notes of any of the newer versions up to the current v4-5.
About the plots that you are showing - I think that it’s hard to compare the spectra if you are using different scales, try putting a logarithmic scale on the plot that you are producing and then compare with your reference.
The other thing is that for what you are trying to obtain, you don’t really need the RADDECAY card. It is only necessary when you want to look at residual dose/isotopes after a certain cooling time. Following this remark, I think that you should not set DCYSCORE for your RESNUCLEI cards as, if I understand correctly, what you are really trying to see is the prompt production of 67Cu.
I tried to run FLUKA with the input file that you provided (I just increased the number of bins in USRTRACK cards and added one more) and I can clearly see a cut-off at 9 MeV in the photon spectra as the minimum energy registered. Then, if the reaction threshold is also 9 MeV, then you would not really see much difference in the photoproduction of 67Cu.
Let me know if you have any other questions.
Cheers,
Jerzy