Message: Re: GenericIon <1GeV Problem Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: GenericIon <1GeV Problem 

Forum: Hadronic Processes
Re: None GenericIon <1GeV Problem (Rick Foster)
Date: 30 Jun, 2006
From: Vladimir IVANTCHENKO <vnivanch@mail.cern.ch>

On Thu, 8 Jun 2006, Rick Foster wrote:

> *** Discussion title: Hadronic Processes
> Email replies to PublicHyperNews@slac.stanford.edu must include:
>   In-Reply-To: <"/hadronprocess/518"@geant4-hn.slac.stanford.edu>
>   Subject: ...change this to be about your reply.
> 
> When I run a simulation with Fe56 > 1 GeV , my model works o.k.
> 
> When I run with Fe56 (or a secondary ion fragment get produced) < ~1GeV,
> I get the following:
> 
> >>> Event 0
> Idle> /gps/energy .9 GeV
> Idle> /run/beamOn 1     
> ### Run 2 start.
> G4EnergyRangeManager:GetHadronicInteraction: counter=1, Ek=16.974333, Material = G4_Al, Element = Al
> *0* low=70, high=10100
> In src/G4EnergyRangeManager.cc, line 120: 
> ===> GetHadronicInteraction: No Model found
> Unrecoverable error for:
>  - Particle energy[GeV] = 0.89963967
>  - Material = G4_Al
>  - Particle type = Fe53[0.0]
> 
> *** G4Exception : 007
>       issued by : G4HadronicProcess
> ChooseHadronicInteraction failed.
> *** Fatal Exception *** core dump ***
> 
> For the GenericIon EM Process I use:
> 
>         else if (particleName == "deuteron" ||
>                  particleName == "triton" ||
>                  particleName == "He3" ||
>                  particleName == "alpha" ||
>                  particleName == "GenericIon" ) {
>         pmanager->AddProcess(new G4MultipleScattering,-1,1,1);
>         pmanager->AddProcess(new G4hIonisation(),-1, 2, 2);
> 
> And for Hadronic Process i use:
> 
>         else if (particleName == "GenericIon") {
>         //Elastic
>         pmanager->AddDiscreteProcess(theElasticProcess);
> 
>         //Inelastic
>         //Pass the ablation Model or use the default de-excitition models
>         G4WilsonAblationModel* AblModel = new G4WilsonAblationModel();
>         G4WilsonAbrasionModel* AbrModel = new G4WilsonAbrasionModel(AblModel);
>         G4IonInelasticProcess* ionInelProcess = new G4IonInelasticProcess();
>         ionInelProcess->RegisterMe(AbrModel);
>         pmanager->AddDiscreteProcess(ionInelProcess);
> 
>         // Load Cross Sections .. Last loaded is used first for energy range
>         // Shen 0-10 GeV/nuc .. Tripathi 0-1 GeV/nuc
>         G4IonsShenCrossSection* shenCS = new G4IonsShenCrossSection();
>         ionInelProcess->AddDataSet(shenCS);
>         G4TripathiCrossSection* tripathiData = new G4TripathiCrossSection();
>         ionInelProcess->AddDataSet(tripathiData);
>         }
> 
> But for some reason, it does not seem to find a model for these <1GeV
> ions.
> 
> Thanks for any insight.
> 
> Rick
> 

Hello,

The problem is that you are trying to use Wilson model, which by default 
is working in limited interval of energies. It seems that you need 
BinaryIonCascade at least for low energies.

Look into example/extended/electromagnetic/TestEm7, 9 
how it can be activated. I would also mension that until now we have no 
clear evidences which of to models is more adequate to your tasks.

VI

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

 Add Message Add Message
to: "Re: GenericIon <1GeV Problem"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]