G4StatMFMicroManager Class Reference

#include <G4StatMFMicroManager.hh>

Data Structures
struct	DeleteFragment

Public Member Functions
G4StatMFChannel *	ChooseChannel (G4int A0, G4int Z0, G4double MeanT)
	G4StatMFMicroManager (const G4Fragment &theFragment, G4int multiplicity, G4double FreeIntE, G4double SCompNuc)
G4double	GetMeanEntropy (void) const
G4double	GetMeanMultiplicity (void) const
G4double	GetMeanTemperature (void) const
G4double	GetProbability (void) const
void	Normalize (G4double Norm)
G4bool	operator!= (const G4StatMFMicroManager &right) const
G4bool	operator== (const G4StatMFMicroManager &right) const
	~G4StatMFMicroManager ()

Private Member Functions
	G4StatMFMicroManager ()
	G4StatMFMicroManager (const G4StatMFMicroManager &right)
void	Initialize (const G4Fragment &theFragment, G4int m, G4double FreeIntE, G4double SCompNuc)
G4bool	MakePartition (G4int k, G4int *ANumbers)
G4StatMFMicroManager &	operator= (const G4StatMFMicroManager &right)

Private Attributes
G4double	_MeanEntropy
G4double	_MeanMultiplicity
G4double	_MeanTemperature
G4double	_Normalization
std::vector< G4StatMFMicroPartition * >	_Partition
G4double	_WW

Detailed Description

Definition at line 43 of file G4StatMFMicroManager.hh.

Constructor & Destructor Documentation

G4StatMFMicroManager() [1/3]

G4StatMFMicroManager::G4StatMFMicroManager	(	const G4Fragment &	theFragment,
		G4int	multiplicity,
		G4double	FreeIntE,
		G4double	SCompNuc
	)

Definition at line 62 of file G4StatMFMicroManager.cc.

                                                             : 
  _Normalization(0.0)
{
  // Perform class initialization
  Initialize(theFragment,multiplicity,FreeIntE,SCompNuc);
}

References Initialize().

~G4StatMFMicroManager()

G4StatMFMicroManager::~G4StatMFMicroManager

(

)

Definition at line 72 of file G4StatMFMicroManager.cc.

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

References _Partition.

G4StatMFMicroManager() [2/3]

G4StatMFMicroManager::G4StatMFMicroManager ( )

inlineprivate

Definition at line 57 of file G4StatMFMicroManager.hh.

{};

G4StatMFMicroManager() [3/3]

G4StatMFMicroManager::G4StatMFMicroManager ( const G4StatMFMicroManager &  right )

private

Definition at line 35 of file G4StatMFMicroManager.cc.

{
    throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMicroManager::copy_constructor meant to not be accessible");
}

Member Function Documentation

ChooseChannel()

G4StatMFChannel * G4StatMFMicroManager::ChooseChannel	(	G4int	A0,
		G4int	Z0,
		G4double	MeanT
	)

Definition at line 170 of file G4StatMFMicroManager.cc.

{
  G4double RandNumber = _Normalization * _WW * G4UniformRand();
  G4double AccumWeight = 0.0;
    
  for (std::vector<G4StatMFMicroPartition*>::iterator i = _Partition.begin();
       i != _Partition.end(); ++i)
    {
    AccumWeight += (*i)->GetProbability();
    if (RandNumber < AccumWeight)
      return (*i)->ChooseZ(A0,Z0,MeanT);
    }
 
  throw G4HadronicException(__FILE__, __LINE__, 
                "G4StatMFMicroCanonical::ChooseChannel: Couldn't find a channel.");
  return 0;
}

References _Normalization, _Partition, _WW, and G4UniformRand.

GetMeanEntropy()

G4double G4StatMFMicroManager::GetMeanEntropy ( void  ) const

inline

Definition at line 85 of file G4StatMFMicroManager.hh.

{return _MeanEntropy; }

References _MeanEntropy.

GetMeanMultiplicity()

G4double G4StatMFMicroManager::GetMeanMultiplicity ( void  ) const

inline

Definition at line 81 of file G4StatMFMicroManager.hh.

{return _MeanMultiplicity; }

References _MeanMultiplicity.

GetMeanTemperature()

G4double G4StatMFMicroManager::GetMeanTemperature ( void  ) const

inline

Definition at line 83 of file G4StatMFMicroManager.hh.

{return _MeanTemperature; }

References _MeanTemperature.

GetProbability()

G4double G4StatMFMicroManager::GetProbability ( void  ) const

inline

Definition at line 77 of file G4StatMFMicroManager.hh.

{return _WW;}

References _WW.

Initialize()

void G4StatMFMicroManager::Initialize ( const G4Fragment &  theFragment, G4int  m, G4double  FreeIntE, G4double  SCompNuc  )

private

Definition at line 81 of file G4StatMFMicroManager.cc.

{
  G4int i;
 
  G4double U = theFragment.GetExcitationEnergy();
 
  G4int A = theFragment.GetA_asInt();
  G4int Z = theFragment.GetZ_asInt();
    
  // Statistical weights
  _WW = 0.0;
 
  // Mean breakup multiplicity
  _MeanMultiplicity = 0.0;
 
  // Mean channel temperature
  _MeanTemperature = 0.0;
 
  // Mean channel entropy
  _MeanEntropy = 0.0;   
    
  // Keep fragment atomic numbers
  //    G4int * FragmentAtomicNumbers = new G4int(static_cast<G4int>(A+0.5));
  //    G4int * FragmentAtomicNumbers = new G4int(m);
  G4int FragmentAtomicNumbers[4];
    
  // We distribute A nucleons between m fragments mantaining the order
  // FragmentAtomicNumbers[m-1]>FragmentAtomicNumbers[m-2]>...>FragmentAtomicNumbers[0]
  // Our initial distribution is 
  // FragmentAtomicNumbers[m-1]=A, FragmentAtomicNumbers[m-2]=0, ..., FragmentAtomicNumbers[0]=0
  FragmentAtomicNumbers[im-1] = A;
  for (i = 0; i <  (im - 1); i++) FragmentAtomicNumbers[i] = 0;
 
  // We try to distribute A nucleons in partitions of m fragments
  // MakePartition return true if it is possible 
  // and false if it is not 
 
  // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
  while (MakePartition(im,FragmentAtomicNumbers)) {
    // Allowed partitions are stored and its probability calculated
            
    G4StatMFMicroPartition * aPartition = new G4StatMFMicroPartition(A,Z);
    G4double PartitionProbability = 0.0;
            
    for (i = im-1; i >= 0; i--) aPartition->SetPartitionFragment(FragmentAtomicNumbers[i]);
    PartitionProbability = aPartition->CalcPartitionProbability(U,FreeIntE,SCompNuc);
    _Partition.push_back(aPartition);
            
    _WW += PartitionProbability;
    _MeanMultiplicity += im*PartitionProbability;
    _MeanTemperature += aPartition->GetTemperature() * PartitionProbability;
    if (PartitionProbability > 0.0) 
      _MeanEntropy += PartitionProbability * aPartition->GetEntropy();
  }
}

References _MeanEntropy, _MeanMultiplicity, _MeanTemperature, _Partition, _WW, A, G4StatMFMicroPartition::CalcPartitionProbability(), G4Fragment::GetA_asInt(), G4StatMFMicroPartition::GetEntropy(), G4Fragment::GetExcitationEnergy(), G4StatMFMicroPartition::GetTemperature(), G4Fragment::GetZ_asInt(), MakePartition(), G4StatMFMicroPartition::SetPartitionFragment(), and Z.

Referenced by G4StatMFMicroManager().

MakePartition()

G4bool G4StatMFMicroManager::MakePartition ( G4int  k, G4int *  ANumbers  )

private

Definition at line 138 of file G4StatMFMicroManager.cc.

{
  G4int l = 1;
  // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
  while (l < k) {
    G4int tmp = ANumbers[l-1] + ANumbers[k-1];
    ANumbers[l-1] += 1;
    ANumbers[k-1] -= 1;
    if (ANumbers[l-1] > ANumbers[l] || ANumbers[k-2] > ANumbers[k-1]) {
      ANumbers[l-1] = 1;
      ANumbers[k-1] = tmp - 1;
      l++;
    } else return true;
  }
  return false;
}

Referenced by Initialize().

Normalize()

void G4StatMFMicroManager::Normalize

(

G4double

Norm

)

Definition at line 158 of file G4StatMFMicroManager.cc.

{
  _Normalization = Norm;
  _WW /= Norm;
  _MeanMultiplicity /= Norm;
  _MeanTemperature /= Norm;
  _MeanEntropy /= Norm; 
    
  return;
}

References _MeanEntropy, _MeanMultiplicity, _MeanTemperature, _Normalization, and _WW.

operator!=()

G4bool G4StatMFMicroManager::operator!=

(

const G4StatMFMicroManager &

right

)

const

Definition at line 56 of file G4StatMFMicroManager.cc.

{
    return true;
}

operator=()

G4StatMFMicroManager & G4StatMFMicroManager::operator= ( const G4StatMFMicroManager &  right )

private

Definition at line 42 of file G4StatMFMicroManager.cc.

{
    throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMicroManager::operator= meant to not be accessible");
    return *this;
}

operator==()

G4bool G4StatMFMicroManager::operator==

(

const G4StatMFMicroManager &

right

)

const

Definition at line 50 of file G4StatMFMicroManager.cc.

{
    return false;
}

Field Documentation

_MeanEntropy

G4double G4StatMFMicroManager::_MeanEntropy

private

Definition at line 115 of file G4StatMFMicroManager.hh.

Referenced by GetMeanEntropy(), Initialize(), and Normalize().

_MeanMultiplicity

G4double G4StatMFMicroManager::_MeanMultiplicity

private

Definition at line 111 of file G4StatMFMicroManager.hh.

Referenced by GetMeanMultiplicity(), Initialize(), and Normalize().

_MeanTemperature

G4double G4StatMFMicroManager::_MeanTemperature

private

Definition at line 113 of file G4StatMFMicroManager.hh.

Referenced by GetMeanTemperature(), Initialize(), and Normalize().

_Normalization

G4double G4StatMFMicroManager::_Normalization

private

Definition at line 109 of file G4StatMFMicroManager.hh.

Referenced by ChooseChannel(), and Normalize().

_Partition

std::vector<G4StatMFMicroPartition*> G4StatMFMicroManager::_Partition

private

Definition at line 103 of file G4StatMFMicroManager.hh.

Referenced by ChooseChannel(), Initialize(), and ~G4StatMFMicroManager().

_WW

G4double G4StatMFMicroManager::_WW

private

Definition at line 107 of file G4StatMFMicroManager.hh.

Referenced by ChooseChannel(), GetProbability(), Initialize(), and Normalize().

---

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

• source/processes/hadronic/models/de_excitation/multifragmentation/include/G4StatMFMicroManager.hh
• source/processes/hadronic/models/de_excitation/multifragmentation/src/G4StatMFMicroManager.cc