No code takes responsibility for decaying charm (and B) mesons
GENIE currently leaves charmed mesons/baryons (e.g. D, etc) as part of the final state of the interaction (status code = 1) to allow then to decay during propagation. Geant4, for reasons to do with correlated production and spin kinematics important for other experiments, does not come configured out-of-the-box to decay these heavy mesons/baryons. It was suggested by the G4 experts that the code should interface to an existing package such as pythia6 (or pythia8).
G4Decay::DoIt : decay table not defined for lambda_c+ *** G4Exception : No Decay Table issued by : G4Decay::DecayIt Decay table is not defined *** This is just a warning message.
Sometimes the particle is a
#1 Updated by Robert Hatcher over 7 years ago
- Status changed from New to Assigned
- Assignee set to Robert Hatcher
- % Done changed from 0 to 50
A GENIE based (vs. Geant4) solution has been implemented through the use of a modified
UserPhysicsOptions.xml file. The file is made available to GENIE through the addition of the file location to the
GXMLPATH environment variable. To enable this, use the corresponding UPS product
setup genie v3665 -q e2:debug setup genie_phyopt v3665 -q dkcharmtau