RE: Ion defnition Fluka / G3

Dear Andreas,

This is an G3 implementation problem and not VMC ! Look at line 71 in gthion.F in G3:
  FACFLU = DME*(Z*DENS/A) when A<=1-6 (vacuum ) then fortran crashes: division per zero !

                                      Cheers, 
                                       Denis    

-----Original Message-----
From: Andreas Morsch [mailto:morsch_at_mail.cern.ch] Sent: Thu 9/17/2009 12:15 PM
To: Bertini, Denis Dr.
Cc: vmc_at_root.cern.ch
Subject: Re: [VMC] Ion defnition Fluka / G3  

Dear Denis,
independently of this crash which we should investigate: Geant3 does not handle any heavy ion interactions.

Best regards - Andreas.

>
>
> Dear All,
>
> I am using the newest version of fluka_vmc together with the examples
> that i foung in geant4_vmc package.
> I am suprised that the way i used for defining an Heavy Ion works in
> the case of Fluka interface but not in the case of Geant3 interface !
>
>
> So what i have done is first to define my User_Ion in the
> Ex03MCApplication.cxx class using the appropriate AddIons() function:
>
> //////////////////////////
> void Ex03MCApplication::AddIons() {
> // Define a Gold Ion
> gMC->DefineIon("MyGold", 79, 197, 79, 0.);
>
> }
> /////////////////////////////////////////////////.
>
> then in the primary generator class Ex03PrimaryGenerator.cxx i add
> the ions as primary to the main VMC stack
>
>
> particle = TDatabasePDG::Instance()->GetParticle("MyGold");
>
> pdg = particle->PdgCode ();
> kinEnergy = 1.0*197; // 1GeV/n
> mass = particle->Mass();
> e = mass + kinEnergy;
>
> pz = sqrt(e*e - mass*mass);
> py = 0.;
> px = 0.;
>
> vx = 0.;
> vx = 0.;
> vx = 0.;
>
>
> cout << "-I- Generate Gold Ion PDG: " << pdg << " Total E: (Gev) " << e << endl;
> fStack->PushTrack(toBeDone, -1, pdg, px, py, pz, e, vx, vy, vz, tof, polx, poly, polz,
> kPPrimary, ntr, 1., 0);
>
> ///////////////////////////////////
>
> the run_fl.C macros accept this and run fine, the run_g3.C crashes in
> the initialisation phase.
> Do you have any idea what is here going on ?
>
> thanks in advance
>
> Denis Bertini
>

-- 
---------------------------------------------------------------------
  Andreas Morsch 
  CERN PH/AIP    
  1211 Geneva 23 
  Switzerland    
tel. +41.22.767.8617
fax  +41.22.767.9480
---------------------------------------------------------------------