#include <G4INCLDecayAvatar.hh>

Inheritance diagram for G4INCL::DecayAvatar:

Collaboration diagram for G4INCL::DecayAvatar:

Public Member Functions
	DecayAvatar (G4INCL::Particle aParticle, G4double time, G4INCL::Nucleus aNucleus, G4bool force=false)

virtual	~DecayAvatar ()

IChannel *	getChannel ()

void	fillFinalState (FinalState *fs)

virtual void	preInteraction ()

virtual void	postInteraction (FinalState *fs)

ParticleList	getParticles () const

std::string	dump () const

Public Member Functions inherited from G4INCL::InteractionAvatar
	InteractionAvatar (G4double, G4INCL::Nucleus , G4INCL::Particle )

	InteractionAvatar (G4double, G4INCL::Nucleus , G4INCL::Particle , G4INCL::Particle *)

virtual	~InteractionAvatar ()

Public Member Functions inherited from G4INCL::IAvatar
	IAvatar ()

	IAvatar (G4double time)

virtual	~IAvatar ()

FinalState *	getFinalState ()

void	fillFinalState (FinalState *fs)

G4double	getTime () const

AvatarType	getType () const

G4bool	isACollision () const

G4bool	isADecay () const

void	setType (AvatarType t)

long	getID () const

std::string	toString ()

Additional Inherited Members
Static Public Member Functions inherited from G4INCL::InteractionAvatar
static void	deleteBackupParticles ()
	Release the memory allocated for the backup particles. More...

Static Public Attributes inherited from G4INCL::InteractionAvatar
static const G4double	locEAccuracy = 1.E-4
	Target accuracy in the determination of the local-energy Q-value. More...

static const G4int	maxIterLocE = 50
	Max number of iterations for the determination of the local-energy Q-value. More...

Protected Member Functions inherited from G4INCL::InteractionAvatar
G4bool	bringParticleInside (Particle *const p)

void	preInteractionLocalEnergy (Particle *const p)
	Apply local-energy transformation, if appropriate. More...

void	preInteractionBlocking ()
	Store the state of the particles before the interaction. More...

void	preInteraction ()

void	postInteraction (FinalState *)

void	restoreParticles () const
	Restore the state of both particles. More...

G4bool	shouldUseLocalEnergy () const
	true if the given avatar should use local energy More...

G4bool	enforceEnergyConservation (FinalState *const fs)
	Enforce energy conservation. More...

Protected Attributes inherited from G4INCL::InteractionAvatar
Nucleus *	theNucleus

Particle *	particle1

Particle *	particle2

ThreeVector	boostVector

G4double	oldTotalEnergy

G4double	oldXSec

G4bool	isPiN

ParticleList	modified

ParticleList	created

ParticleList	modifiedAndCreated

Protected Attributes inherited from G4INCL::IAvatar
G4double	theTime

Static Protected Attributes inherited from G4INCL::InteractionAvatar
static G4ThreadLocal Particle *	backupParticle1 = NULL

static G4ThreadLocal Particle *	backupParticle2 = NULL

Detailed Description

Decay avatar

The reflection avatar is created when a particle reaches the boundary of the nucleus. At this point it can either be reflected from the boundary or exit the nucleus.

Definition at line 55 of file G4INCLDecayAvatar.hh.

Constructor & Destructor Documentation

G4INCL::DecayAvatar::DecayAvatar	(	G4INCL::Particle *	aParticle,
		G4double	time,
		G4INCL::Nucleus *	aNucleus,
		G4bool	force = `false`
	)

Definition at line 49 of file G4INCLDecayAvatar.cc.

     : InteractionAvatar(time, n, aParticle), forced(force),
       incidentDirection(aParticle->getMomentum())
   {
     setType(DecayAvatarType);
   }

Here is the call graph for this function:

G4INCL::DecayAvatar::~DecayAvatar ( )

virtual

Definition at line 56 of file G4INCLDecayAvatar.cc.

                             {
 
   }

Member Function Documentation

std::string G4INCL::DecayAvatar::dump ( ) const

virtual

Implements G4INCL::IAvatar.

Definition at line 174 of file G4INCLDecayAvatar.cc.

                                     {
     std::stringstream ss;
     ss << "(avatar " << theTime << " 'decay" << '\n'
       << "(list " << '\n'
       << particle1->dump()
       << "))" << '\n';
     return ss.str();
   }

Here is the call graph for this function:

void G4INCL::DecayAvatar::fillFinalState ( FinalState * fs )

G4INCL::IChannel * G4INCL::DecayAvatar::getChannel ( )

virtual

Implements G4INCL::InteractionAvatar.

Definition at line 60 of file G4INCLDecayAvatar.cc.

                                           {
     if(particle1->isDelta()) {
       INCL_DEBUG("DeltaDecayChannel chosen." << '\n');
       return new DeltaDecayChannel(particle1, incidentDirection);
     }
     else if(particle1->isEta() || particle1->isOmega()) {
       INCL_DEBUG("PionResonanceDecayChannel chosen." << '\n');
       return new PionResonanceDecayChannel(particle1, incidentDirection);
     }  
     else
       return NULL;
   }

Here is the call graph for this function:

ParticleList G4INCL::DecayAvatar::getParticles ( ) const

inlinevirtual

Implements G4INCL::IAvatar.

Definition at line 66 of file G4INCLDecayAvatar.hh.

                                       {
       ParticleList theParticleList;
       theParticleList.push_back(particle1);
       return theParticleList;
     }

void G4INCL::DecayAvatar::postInteraction ( FinalState * fs )

virtual

Implements G4INCL::IAvatar.

Definition at line 77 of file G4INCLDecayAvatar.cc.

                                                   {
     // Make sure we have at least two particles in the final state
 // assert(fs->getModifiedParticles().size() + fs->getCreatedParticles().size() - fs->getDestroyedParticles().size() >= 2);
 
     if(!forced) { // Normal decay
 
       // Call the postInteraction method of the parent class
       // (provides Pauli blocking and enforces energy conservation)
       InteractionAvatar::postInteraction(fs);
 
       if(fs->getValidity() == PauliBlockedFS)
         /* If the decay was Pauli-blocked, make sure the propagation model
          * generates a new decay avatar on the next call to propagate().
          *
          * \bug{Note that we don't generate new decay avatars for deltas that
          * could not satisfy energy conservation. This is in keeping with
          * INCL4.6, but doesn't seem to make much sense to me (DM), as energy
          * conservation can be impossible to satisfy due to weird local-energy
          * conditions, for example, that evolve with time.}
          */
         fs->addModifiedParticle(particle1);
 
     } else { // Forced decay
       created = fs->getCreatedParticles();
       modified = fs->getModifiedParticles();
 
       // Try to enforce energy conservation
       fs->setTotalEnergyBeforeInteraction(oldTotalEnergy);
       const G4bool success = enforceEnergyConservation(fs);
       if(!success) {
         INCL_DEBUG("Enforcing energy conservation: failed!" << '\n');
 
         if(theNucleus) {
           // Restore the state of the initial particles
           restoreParticles();
 
           // Delete newly created particles
           for(ParticleIter i=created.begin(), e=created.end(); i!=e; ++i )
             delete *i;
 
           fs->reset();
           fs->makeNoEnergyConservation();
           fs->setTotalEnergyBeforeInteraction(0.0);
 
           return; // Interaction is blocked. Return an empty final state.
         } else {
           // If there is no nucleus we have to continue anyway, even if energy
           // conservation failed. We cannot afford producing unphysical
           // remnants.
           INCL_DEBUG("No nucleus, continuing anyway." << '\n');
         }
       } else {
         INCL_DEBUG("Enforcing energy conservation: success!" << '\n');
       }
 
       if(theNucleus) {
         // Test CDPP blocking
         G4bool isCDPPBlocked = Pauli::isCDPPBlocked(created, theNucleus);
 
         if(isCDPPBlocked) {
           INCL_DEBUG("CDPP: Blocked!" << '\n');
 
           // Restore the state of both particles
           restoreParticles();
 
           // Delete newly created particles
           for(ParticleIter i=created.begin(), e=created.end(); i!=e; ++i )
             delete *i;
 
           fs->reset();
           fs->makePauliBlocked();
           fs->setTotalEnergyBeforeInteraction(0.0);
 
           return; // Interaction is blocked. Return an empty final state.
         }
         INCL_DEBUG("CDPP: Allowed!" << '\n');
 
       }
     }
 
     // If there is a nucleus, increment the counters
     if(theNucleus) {
       switch(fs->getValidity()) {
         case PauliBlockedFS:
           theNucleus->getStore()->getBook().incrementBlockedDecays();
           break;
         case NoEnergyConservationFS:
         case ParticleBelowFermiFS:
         case ParticleBelowZeroFS:
           break;
         case ValidFS:
           theNucleus->getStore()->getBook().incrementAcceptedDecays();
       }
     }
     return;
   }

Here is the call graph for this function:

void G4INCL::DecayAvatar::preInteraction ( )

virtual

Implements G4INCL::IAvatar.

Definition at line 73 of file G4INCLDecayAvatar.cc.

                                    {
     InteractionAvatar::preInteraction();
   }

Here is the call graph for this function:

The documentation for this class was generated from the following files:

source/geant4.10.03.p03/source/processes/hadronic/models/inclxx/incl_physics/include/G4INCLDecayAvatar.hh
source/geant4.10.03.p03/source/processes/hadronic/models/inclxx/incl_physics/src/G4INCLDecayAvatar.cc

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation