#include <G4StatMFChannel.hh>

Public Member Functions
	G4StatMFChannel ()

	~G4StatMFChannel ()

void	CreateFragment (G4int A, G4int Z)

size_t	GetMultiplicity (void)

G4bool	CheckFragments (void)

G4double	GetFragmentsCoulombEnergy (void)

G4double	GetFragmentsEnergy (G4double T) const

G4FragmentVector *	GetFragments (G4int anA, G4int anZ, G4double T)

Detailed Description

Definition at line 41 of file G4StatMFChannel.hh.

Constructor & Destructor Documentation

◆ G4StatMFChannel()

G4StatMFChannel::G4StatMFChannel ( )

Definition at line 59 of file G4StatMFChannel.cc.

                                 : 
  _NumOfNeutralFragments(0), 
  _NumOfChargedFragments(0)
{}

◆ ~G4StatMFChannel()

G4StatMFChannel::~G4StatMFChannel ( )

Definition at line 64 of file G4StatMFChannel.cc.

{ 
  if (!_theFragments.empty()) {
    std::for_each(_theFragments.begin(),_theFragments.end(),
          DeleteFragment());
  }
}

Member Function Documentation

◆ CheckFragments()

G4bool G4StatMFChannel::CheckFragments ( void )

Definition at line 72 of file G4StatMFChannel.cc.

{
    std::deque<G4StatMFFragment*>::iterator i;
    for (i = _theFragments.begin(); 
     i != _theFragments.end(); ++i) 
      {
    G4int A = (*i)->GetA();
    G4int Z = (*i)->GetZ();
        if ( (A > 1 && (Z > A || Z <= 0)) || (A==1 && Z > A) || A <= 0 ) return false;
    }
    return true;
}

Referenced by G4StatMF::BreakItUp().

◆ CreateFragment()

void G4StatMFChannel::CreateFragment	(	G4int	A,
		G4int	Z
	)

Definition at line 85 of file G4StatMFChannel.cc.

{
    if (Z <= 0.5) {
    _theFragments.push_back(new G4StatMFFragment(A,Z));
    _NumOfNeutralFragments++;
    } else {
    _theFragments.push_front(new G4StatMFFragment(A,Z));
    _NumOfChargedFragments++;
    }
    
    return;
}

Referenced by G4StatMFMicroCanonical::ChooseAandZ(), and G4StatMFMicroPartition::ChooseZ().

◆ GetFragments()

G4FragmentVector * G4StatMFChannel::GetFragments	(	G4int	anA,
		G4int	anZ,
		G4double	T
	)

Definition at line 127 of file G4StatMFChannel.cc.

{
    // calculate momenta of charged fragments  
    CoulombImpulse(anA,anZ,T);
    
    // calculate momenta of neutral fragments
    FragmentsMomenta(_NumOfNeutralFragments, _NumOfChargedFragments, T);
 
    G4FragmentVector * theResult = new G4FragmentVector;
    std::deque<G4StatMFFragment*>::iterator i;
    for (i = _theFragments.begin(); i != _theFragments.end(); ++i)
    theResult->push_back((*i)->GetFragment(T));
 
    return theResult;
}

◆ GetFragmentsCoulombEnergy()

G4double G4StatMFChannel::GetFragmentsCoulombEnergy ( void )

Definition at line 101 of file G4StatMFChannel.cc.

{
    G4double Coulomb = std::accumulate(_theFragments.begin(),_theFragments.end(),
                     0.0,SumCoulombEnergy());
//      G4double Coulomb = 0.0;
//      for (unsigned int i = 0;i < _theFragments.size(); i++)
//      Coulomb += _theFragments[i]->GetCoulombEnergy();
    return Coulomb;
}

◆ GetFragmentsEnergy()

G4double G4StatMFChannel::GetFragmentsEnergy ( G4double T ) const

Definition at line 113 of file G4StatMFChannel.cc.

{
    G4double Energy = 0.0;
    
    G4double TranslationalEnergy = (3./2.)*T*static_cast<G4double>(_theFragments.size());
 
    std::deque<G4StatMFFragment*>::const_iterator i;
    for (i = _theFragments.begin(); i != _theFragments.end(); ++i)
      {
    Energy += (*i)->GetEnergy(T);
      }
    return Energy + TranslationalEnergy;    
}

◆ GetMultiplicity()

size_t G4StatMFChannel::GetMultiplicity ( void )

inline

Definition at line 65 of file G4StatMFChannel.hh.

65{ return _theFragments.size();}