Message: avoid methastable daughter nuclid Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None avoid methastable daughter nuclid 

Forum: Physics List
Re: None using G4PhotonEvaporation (Karen Saucke)
Date: 03 May, 2005
From: Karen Saucke <Karen Saucke>

I solved the decay of the metastable Ba137 (as daughter of Cs137) in the following way:

my physics List is still the old one, still no G4PhotonEvaporation used.
while( (*theParticleIterator)() ){
		G4ParticleDefinition* particle = theParticleIterator->value();
		G4ProcessManager* pmanager = particle->GetProcessManager();
		if (theDecayProcess->IsApplicable(*particle)) { 
			pmanager ->AddProcess(theDecayProcess);
			// Setzen der Reihenfolge fuer PostStepDoIt und AtRestDoIt
			pmanager ->SetProcessOrdering(theDecayProcess, idxPostStep,1);
			pmanager ->SetProcessOrdering(theDecayProcess, idxAtRest,1);
		if (theRadioactiveDecay->IsApplicable(*particle)) { 
			pmanager ->AddProcess(theRadioactiveDecay);
			// Setzen der Reihenfolge fuer PostStepDoIt und AtRestDoIt
			pmanager ->SetProcessOrdering(theRadioactiveDecay, idxPostStep,2);
			pmanager ->SetProcessOrdering(theRadioactiveDecay, idxAtRest,2);


but I changed the life-lime of the nuclid by hand in 'G4UserStackingAction'

   TarDet01StackingAction::ClassifyNewTrack(const G4Track * aTrack)
	#ifdef TALK
		G4cout<<" G4ClassificationOfNewTrack TarDet01StackingAction::ClassifyNewTrack(const G4Track * aTrack) "<<"\r\n"; 

	G4TouchableHistory* theTouchable = (G4TouchableHistory*)(aTrack->GetTouchable());

	G4ParticleDefinition * particleType = aTrack->GetDefinition();   

	//  if(particleType->GetParticleName()=="Ba137[661.7]") {




first I set the lifetime to zero only for Ba137. With above version you may be carefull about long living radioactive nuclids.

That helped for me, may for someone else as well.


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

 Add Message Add Message
to: "avoid methastable daughter nuclid"

 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 ]