Extended particle list in TGeant3

From: Yuri Kharlov <yuri.kharlov_at_ihep.ru>
Date: Wed, 05 Jan 2005 20:25:49 +0300

Dear VMC developpers,

The list of particles defined in TGeant3::DefineParticles() seems to require some corrections. J/psi is defined as a particle 66:

  Gspart(fNG3Particles++, "J/Psi",       3, 3.09688, 0., 0.);
  fPDGCode[fNPDGCodes++]=443;          // 66 = J/Psi

First, its lifetime (last paratemer) is set to 0 which causes fpe in gtrak/gtneut.F:234.
Second, decay modes which were intended for J/psi, are defined for a particle 113, not 66:
// --- jpsi ---

    for (kz = 0; kz < 6; ++kz) {

        bratio[kz] = 0.;
        mode[kz] = 0;

    }
    ipa = 113;
    bratio[0] = 50.;
    bratio[1] = 50.;
    mode[0] = 506;
    mode[1] = 605;
    Gsdk(ipa, bratio, mode);
I expect to have ipa=66 instead of ipa=66. The same concerns some other particles (upsilon, phi) which are defined in TGeant3 beyond the standard geant3.

I use root 4.01/02, geant3 checked out CVS of late September 2004 (TGeant3 has cvs version 1.34), gcc 3.2.3.

Best regards,
Yuri Kharlov. Received on Wed Jan 05 2005 - 18:25:55 MET