#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 46 of file G4StatMFChannel.cc.

                                 : 
  _NumOfNeutralFragments(0), 
  _NumOfChargedFragments(0)
{
  Pos.resize(8);
  Vel.resize(8);
  Accel.resize(8);
}

◆ ~G4StatMFChannel()

G4StatMFChannel::~G4StatMFChannel ( )

Definition at line 55 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 63 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 76 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 122 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 92 of file G4StatMFChannel.cc.

{
  G4double Coulomb =
    std::accumulate(_theFragments.begin(),_theFragments.end(),
                    0.0,
                    [](const G4double& running_total,
                       G4StatMFFragment*& fragment)
                    {
                      return running_total + fragment->GetCoulombEnergy();
                    } );
  //      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 108 of file G4StatMFChannel.cc.

{
  G4double Energy = 0.0;
    
  G4double TranslationalEnergy = 1.5*T*_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();}