Message: Photonuclear reactions Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Photonuclear reactions 

Forum: Hadronic Processes
Date: 25 Jan, 2010
From: Ralph Massarczyk <Ralph Massarczyk>

Is there a possibility to describe the following reactions? - Ag + gamma -> Ag* , and after that: - Ag* -> Ag + gamma I implemented in my Physicslist

  G4GammaNuclearReaction* gammaNuclear = new G4GammaNuclearReaction();
  G4PhotoNuclearProcess* photoNuclearProcess = new G4PhotoNuclearProcess();
  photoNuclearProcess->RegisterMe(gammaNuclear);

and for Decay:

#include "G4Decay.hh"
#include "G4RadioactiveDecay.hh"
#include "G4IonTable.hh"
#include "G4Ions.hh"

void PhysicsList::ConstructGeneral() {

  // Add Decay Process
  G4Decay* theDecayProcess = new G4Decay();
  theParticleIterator->reset();
  while( (*theParticleIterator)() )
    {
      G4ParticleDefinition* particle = theParticleIterator->value();
      G4ProcessManager* pmanager = particle->GetProcessManager();

      if (theDecayProcess->IsApplicable(*particle) && !particle->IsShortLived()) 
	{ 
	  pmanager ->AddProcess(theDecayProcess);
	  // set ordering for PostStepDoIt and AtRestDoIt
	  pmanager ->SetProcessOrdering(theDecayProcess, idxPostStep);
	  pmanager ->SetProcessOrdering(theDecayProcess, idxAtRest);
	}
    }

  // Declare radioactive decay to the GenericIon in the IonTable.
  const G4IonTable *theIonTable = 
  G4ParticleTable::GetParticleTable()->GetIonTable();
  G4RadioactiveDecay *theRadioactiveDecay = new G4RadioactiveDecay();

  for (G4int i=0; i<theIonTable->Entries(); i++) 
    {
      G4String particleName = theIonTable->GetParticle(i)->GetParticleName();
      G4String particleType = theIonTable->GetParticle(i)->GetParticleType();

      if (particleName == "GenericIon") 
	{
	  G4ProcessManager* pmanager = 
	    theIonTable->GetParticle(i)->GetProcessManager();
	  pmanager->SetVerboseLevel(VerboseLevel);
	  pmanager ->AddProcess(theRadioactiveDecay);

	  pmanager ->SetProcessOrdering(theRadioactiveDecay, idxPostStep);
	  pmanager ->SetProcessOrdering(theRadioactiveDecay, idxAtRest);

	} 
    }
}

But unfortunally I see no gammas from excited states, only if I use excited nuclei oder radioactive nuclei as PrimaryParticles. But I need the excitation also. Is this case described somewhere in one of the examples?

I read something about a G4NRF package, is this still available?

Thanks for your help

Ralph Massarczyk

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

1 News: Re: Photonuclear reactions   (Ralph Massarczyk - 26 Jan, 2010)
(_ Idea: Re: Photonuclear reactions   (Zach Hartwig - 26 Jan, 2010)
(_ Question: Re: Photonuclear reactions   (Ralph Massarczyk - 29 Jan, 2010)
 Add Message Add Message
to: "Photonuclear reactions"

 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 ]