neutrino interactions with pythia

switch on neutral and charged current interactions in pythia. either by extending the existing interaction process or adding a new corsika::pythia::NeutrinoInteraction with different initialization of pythia.

changed milestone to %neutrino showers

added Feature request label

marked this issue as related to #369 (closed)

@Fan could you create a branch with the solution you presented during the technical call. Or a link to your repository?
I think if we adjust the initialization of the existing pythia in c8 it should be straight forward. We should be able to treat neutrino's via the process sequence as all other particles. The workaround you presented with creating the interaction externally and feeding them into c8 is super cool though. we should keep it as our stack-in solution.

created branch 569-neutrino-interactions-with-pythia to address this issue

mentioned in merge request !548 (closed)

mentioned in merge request !427 (merged)

explicitly allow weak exchange:
Pythia8::Pythia::readString("WeakBosonExchange:ff2ff(t:W) = on");

see presentation by Fan Hu on neutrino interactions in C8

@mreininghaus I think we may not get away without a cross section for neutrinos after all. Looking at cascade.inl even when an interaction is forced, the interaction function re-evaluates the cross section (because energy loss?). Also the cross section enters again when the target is selected in selectInteraction. Or do I read this wrong?

Yes, you're right. But the actual value can be arbitrarily small to make it practically irrelevant. We could put a quectobarn .

Ah yes. I knew it was going to pay off to add those ultrasmall and ultralarge SI prefixes(?)

A primitive test confirms this. Adding if (is_neutrino(projectileId)) return 1_mb; to getCrossSection in pythia interface makes corsika.cpp produce a shower when called with a primary neutrino and forcing the interaction

changed milestone to %Beta Release

assigned to @riehn

see branch neutrino-interactions-with-pythia-8245. It implements interactions of primary neutrinos with pythia 8.245 based on the example program for DIS main36.cc. Curiously enough that example does not work with neutrinos since pythia 8.3

mentioned in merge request !593 (merged)

mentioned in merge request !594 (merged)

Closed by !594 (merged)

closed