Message: gold activation with neutrons Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question gold activation with neutrons 

Keywords: neutron gold activation isotope production
Forum: Hadronic Processes
Date: 25 Jan, 2008
From: Sylvia Studeny <sylvia.studeny@helmholtz-muenchen.de>

Hi,

I have a problem with counting Au-197(n,gamma)Au-198 reaction: I do this in two ways:

1) scoring the neutron flux and folding this with the capture crosssection

2) counting the number of Au-198 isotopes directly in Stepping action.

The fist way is ok.
The second way give to few Au-198 (only 30% of the expected). I count them in the following way:
In UserSteppingAction:
 //Energiedeposition im Detektor:
 physiVol = aStep->GetPreStepPoint()->GetPhysicalVolume();
//  physiVol = aStep->GetPostStepPoint()->GetPhysicalVolume();
 if (physiVol->GetName()!= "Voxel"){return;}

  aTrack = aStep->GetTrack();
//- check for newly produced particle. e.g. first step.
  if ( aTrack->GetCurrentStepNumber() != 1) return;
  //- check for this is not a primary particle. e.g. ParentID > 0 .
  if ( aTrack->GetParentID() == 0 ) return;

  Particle = aTrack->GetDefinition();
  particleName = Particle -> GetParticleName();
  //G4cout << "Stepping: TeilchenName: " << particleName << G4endl;
  if(particleName=="Au198[0.0]"){au198No++;}
  if(particleName=="Au197[0.0]"){au197No++;}
  return;

My neutron physics looks like that:

   pManager = G4Neutron::Neutron()->GetProcessManager();

   G4HadronElasticProcess* thenElasticProcess 
                         = new G4HadronElasticProcess();
   thenElasticProcess->AddDataSet(new G4NeutronHPElasticData);

   G4NeutronHPThermalScatteringData * theNeutronThermalData =
         new G4NeutronHPThermalScatteringData;
   thenElasticProcess->AddDataSet(theNeutronThermalData);

   G4LElastic* thenElasticModel = new G4LElastic();
   thenElasticModel->SetMinEnergy(19.*MeV);
   G4NeutronHPElastic* theNeutronHPElasticModel= new G4NeutronHPElastic();
   theNeutronHPElasticModel ->SetMaxEnergy(20.*MeV);
   theNeutronHPElasticModel ->SetMinEnergy(4.*eV);
   thenElasticProcess->RegisterMe(thenElasticModel);
   thenElasticProcess->RegisterMe(theNeutronHPElasticModel);
   G4NeutronHPThermalScattering* theNeutronHPThermalModel = new G4NeutronHPThermalScattering();
   theNeutronHPThermalModel->SetMaxEnergy(4.*eV);
   thenElasticProcess->RegisterMe(theNeutronHPThermalModel);
   pManager->AddDiscreteProcess(thenElasticProcess);

   G4NeutronInelasticProcess* theNeutronInelasticProcess 
                         = new G4NeutronInelasticProcess(); 
   theNeutronInelasticProcess->AddDataSet(new G4NeutronHPInelasticData);

   G4NeutronHPInelastic* theNeutronHPInelasticModel= new G4NeutronHPInelastic();
   theNeutronHPInelasticModel ->SetMaxEnergy(20.*MeV);
   G4BinaryCascade* theNeutronBCModel = new G4BinaryCascade();
   theNeutronBCModel->SetMinEnergy(19.*MeV);
   theNeutronInelasticProcess->RegisterMe(theNeutronHPInelasticModel);
   theNeutronInelasticProcess->RegisterMe(theNeutronBCModel);

   G4NeutronIsotopeProduction *aProductionModel = new G4NeutronIsotopeProduction;
  theNeutronInelasticProcess->RegisterIsotopeProductionModel(aProductionModel);
  G4HadronicProcess:: EnableIsotopeProductionGlobally();

   pManager->AddDiscreteProcess(theNeutronInelasticProcess);

   G4HadronFissionProcess* thenFission
                         = new G4HadronFissionProcess();
   thenFission->AddDataSet(new G4NeutronHPFissionData);

   G4LFission* thenFissionModel = new G4LFission();
   thenFissionModel->SetMinEnergy(19.*MeV);
   G4NeutronHPFission* theNeutronHPFissionModel= new G4NeutronHPFission();
   theNeutronHPFissionModel ->SetMaxEnergy(20.*MeV);
   thenFission->RegisterMe(thenFissionModel);
   thenFission->RegisterMe(theNeutronHPFissionModel);
   pManager->AddDiscreteProcess(thenFission);

   G4HadronCaptureProcess* thenCapture
                         = new G4HadronCaptureProcess();
   thenCapture->AddDataSet(new G4NeutronHPCaptureData);

   G4LCapture* thenCaptureModel = new G4LCapture();
   thenCaptureModel->SetMinEnergy(19.*MeV);
   G4NeutronHPCapture* theNeutronHPCaptureModel= new G4NeutronHPCapture();
   theNeutronHPCaptureModel ->SetMaxEnergy(20.*MeV);
   thenCapture->RegisterMe(thenCaptureModel);
   thenCapture->RegisterMe(theNeutronHPCaptureModel);
   pManager->AddDiscreteProcess(thenCapture);

I read the older forum-discussion about gold activation, but I could not detect my fault there. Do I do something wrong? I assume that all particles created by elastic/inelastic interaction of the neutron are at least once in Stepping action, where I can count them. Is that right? If not, where should I count the activated gold atoms?

Thanks for any suggestions, Sylvia

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

1 None: Re: gold activation with neutrons   (Luciano - 29 Jan, 2008)
1 None: Re: gold activation with neutrons   (Sylvia Studeny - 29 Jan, 2008)
3 None: Re: gold activation with neutrons   (Luciano - 29 Jan, 2008)
1 None: Re: gold activation with neutrons   (Sylvia Studeny - 29 Jan, 2008)
(_ Question: Re: gold activation with neutrons   (Franciole Marinho - 11 Mar, 2009)
(_ Question: Re: gold activation with neutrons   (Florestan - 24 Nov, 2010)
(_ None: Re: gold activation with neutrons   (Colin Paul Gloster - 24 Nov, 2010)
(_ None: Re: gold activation with neutrons   (Eleni - 10 Mar, 2012)
(_ Feedback: Re: gold activation with neutrons   (Sean Turnbull - 11 Mar, 2012)
5 None: Re: gold activation with neutrons   (Eleni - 12 Mar, 2012)
 Add Message Add Message
to: "gold activation with 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 ]