Ion defnition Fluka / G3

From: Bertini, Denis Dr. <D.Bertini_at_gsi.de>
Date: Thu, 17 Sep 2009 11:58:37 +0200

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 Received on Thu Sep 17 2009 - 11:58:44 CEST