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

None GenericIon <1GeV Problem 

Keywords: GenericIon, Missing Model
Forum: Hadronic Processes
Date: 08 Jun, 2006
From: Rick Foster <rickf@space.mit.edu>

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

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

1 None: Re: GenericIon <1GeV Problem   (Dennis H. Wright - 22 Jun, 2006)
1 None: Re: GenericIon <1GeV Problem   (Vladimir IVANTCHENKO - 23 Jun, 2006)
3 None: Re: GenericIon <1GeV Problem   (Vladimir IVANTCHENKO - 30 Jun, 2006)
 Add Message Add Message
to: "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 ]