How to simulate D mesons in GSTAR
By Pavel Nevski
' ====================================================== '
D meason has many decays (4 pages in PDG handbook for D+, 5 pages for D0)
and GEANT may have not more than 6 decay channels. Also, the number
of products in each decay in geant is limited by 3.
If you want to have a general (more or less complete) simulation of a
D-meson decay have to generate and decay D meson in an external program
(i.e. PYTHIA). Then you have to write the decay product (including D-meason
itself with it's PDG code) in a txt file as shown at GSTAR web or in an
n-tuple using interface in $STAR/pams/gen/heptup.f
If you are interested in few selected decay channels, you can add D meson
and selected channels to $STAR/pams/sim/gstar/gstar_part.g
There is already an example put there as a comment:
* For D+, include only K0_S + pi+ decay mode pmj 16/2/95
* Particle D_plus code=35 TrkTyp=4 mass=1.869 charge=1 tlife=1.057e-12,
* pdg=411 bratio= { 1, } mode= { 1608, }
* For D0, include only K+ + pi- decay mode pmj 16/2/95
* Particle D0 code=37 TrkTyp=3 mass=1.865 charge=0 tlife=0.415e-12,
* pdg=421 bratio= { 1, } mode= { 1208, }
It is up to you to uncomment it and to tune the parameters as you like