Message: not getting neutron in volume Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None not getting neutron in volume 

Forum: Particles
Date: 20 Sep, 2011
From: Summit <Summit>

Hello Geant4Users,

In PhysicsList.cc, I have added

(name == "QGSP_BIC_HP") {
    hadronPhys.push_back(new HadronPhysicsQGSP_BIC_HP());
to get neutron.

1) During run, I am getting this error

physicsList->Construct() start. /u/ey/perl/geant4/geant4.9.4/data/G4NDL3.14

NeutronHP: /Capture file for Z = 8, A = 18 is not found and NeutronHP will use /u/ey/perl/geant4/geant4.9.4/data/G4NDL3.14/Capture//CrossSection/8_17_Oxygen
NeutronHP: /Elastic file for Z = 8, A = 18 is not found and NeutronHP will use /u/ey/perl/geant4/geant4.9.4/data/G4NDL3.14/Elastic//CrossSection/8_17_Oxygen
NeutronHP: /Inelastic file for Z = 8, A = 18 is not found and NeutronHP will use /u/ey/perl/geant4/geant4.9.4/data/G4NDL3.14/Inelastic//CrossSection/8_17_Oxygen
NeutronHP: /Capture/ file for Z = 8, A = 18 is not found and NeutronHP will use /u/ey/perl/geant4/geant4.9.4/data/G4NDL3.14/Capture///CrossSection/8_17_Oxygen
physicsList->CheckParticleList() start.

2) for neutron, I have added

else if (particleName == "neutron") {

	  // 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);

	  // inelastic scattering

	  G4NeutronInelasticProcess* theInelasticProcess = new G4NeutronInelasticProcess();

G4HENeutronInelastic* theHEInelasticModel = new G4HENeutronInelastic;

	  theInelasticProcess->RegisterMe(theHEInelasticModel);

	  G4NeutronHPInelastic * theLENeutronInelasticModel = new G4NeutronHPInelastic;

	  theInelasticProcess->RegisterMe(theLENeutronInelasticModel);

	  G4NeutronHPInelasticData * theNeutronData1 = new G4NeutronHPInelasticData;

	  theInelasticProcess->AddDataSet(theNeutronData1);

	  pmanager->AddDiscreteProcess(theInelasticProcess);

	  // capture

	  G4HadronCaptureProcess* theCaptureProcess = new G4HadronCaptureProcess;

	  G4LCapture* theCaptureModel = new G4LCapture;

	  theCaptureModel->SetMinEnergy(19*MeV);

	  theCaptureProcess->RegisterMe(theCaptureModel);

	  G4NeutronHPCapture * theLENeutronCaptureModel = new G4NeutronHPCapture;

	  theCaptureProcess->RegisterMe(theLENeutronCaptureModel);

	  G4NeutronHPCaptureData * theNeutronData3 = new G4NeutronHPCaptureData;

	  theCaptureProcess->AddDataSet(theNeutronData3);

	  pmanager->AddDiscreteProcess(theCaptureProcess);
}

 but neutrons were not observed during interaction.

Please suggest me what I do to observe neutrons.

Thanks in advance,

Summit

 Add Message Add Message
to: "not getting neutron in volume"

 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 ]