Message: Problems including physics process list for neutrons Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Problems including physics process list for neutrons  

Forum: Physics List
Date: 18 Nov, 2010
From: Manuel Castro Avila <Manuel Castro Avila>

Hi everyone,

I am creating a physics process list for neutrons. So, I am using the following code:

 else if (particleName == "neutron") {
        // Inelastic
        G4NeutronInelasticProcess* theInelasticProcess = new G4NeutronInelasticProcess("inelastic");

        //PreEquilibrium
        G4PreCompoundModel * thePreEquilib = new G4PreCompoundModel(new G4ExcitationHandler);
        theInelasticProcess->RegisterMe(thePreEquilib);
        G4NeutronHPInelasticData* theNeutronData1 = new G4NeutronHPInelasticData;
        theInelasticProcess->AddDataSet(theNeutronData1);
	// Isotope production
        G4NeutronIsotopeProduction* theNeutronIsotopeProduction = new G4NeutronIsotopeProduction;
        theInelasticProcess->RegisterIsotopeProductionModel(theNeutronIsotopeProduction);
        pmanager->AddDiscreteProcess(theInelasticProcess);
	// Capture
        G4HadronCaptureProcess* theCaptureProcess = new G4HadronCaptureProcess;
        G4NeutronHPCapture* theLENeutronCaptureModel = new G4NeutronHPCapture;
        theCaptureProcess->RegisterMe(theLENeutronCaptureModel);
        G4NeutronHPCaptureData* theNeutronData3 = new G4NeutronHPCaptureData;
        theCaptureProcess->AddDataSet(theNeutronData3);
        pmanager->AddDiscreteProcess(theCaptureProcess);

      // elastic scattering
         G4HadronElasticProcess* theNeutronElasticProcess = new G4HadronElasticProcess;
         G4LElastic* theElasticModel1 = new G4LElastic;
         G4NeutronHPElastic * theElasticNeutron = new G4NeutronHPElastic;
         theNeutronElasticProcess->RegisterMe(theElasticModel1);
         theElasticModel1->SetMinEnergy(19*MeV);
         theNeutronElasticProcess->RegisterMe(theElasticNeutron);
         G4NeutronHPElasticData * theNeutronData = new G4NeutronHPElasticData;
         theNeutronElasticProcess->AddDataSet(theNeutronData);
         pmanager->AddDiscreteProcess(theNeutronElasticProcess);
      }

But when I execute the code appear a error:

*** G4Exception : No GenericIon
      issued by : G4VUserPhysicsList::AddProcessManager
GenericIon does not exist
*** Run Must Be Aborted 

I don't understand what is happening or if I need to implement the physics process list for ions, too?

Do anyone know where is the problem??

I thank a lot your help

Manuel

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

1 Idea: Re: Problems including physics process list for neutrons   (Aatif Imtiaz Butt - 18 Nov, 2010)
(_ None: Re: Problems including physics process list for neutrons   (Manuel Castro Avila - 19 Nov, 2010)
 Add Message Add Message
to: "Problems including physics process list for neutrons "

 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 ]