G4StatMFMicroManager Class Reference

#include <G4StatMFMicroManager.hh>

Collaboration diagram for G4StatMFMicroManager:

Classes
struct	DeleteFragment

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

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

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

Detailed Description

Definition at line 44 of file G4StatMFMicroManager.hh.

Constructor & Destructor Documentation

G4StatMFMicroManager() [1/3]

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

Definition at line 63 of file G4StatMFMicroManager.cc.

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

Here is the call graph for this function:

~G4StatMFMicroManager()

G4StatMFMicroManager::~G4StatMFMicroManager

(

)

Definition at line 73 of file G4StatMFMicroManager.cc.

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

G4StatMFMicroManager() [2/3]

G4StatMFMicroManager::G4StatMFMicroManager ( )

inlineprivate

Definition at line 58 of file G4StatMFMicroManager.hh.

{};

Here is the call graph for this function:

G4StatMFMicroManager() [3/3]

G4StatMFMicroManager::G4StatMFMicroManager ( const G4StatMFMicroManager &  right )

private

Definition at line 36 of file G4StatMFMicroManager.cc.

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

Here is the call graph for this function:

Member Function Documentation

ChooseChannel()

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

Definition at line 171 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;
}

Here is the caller graph for this function:

GetMeanEntropy()

G4double G4StatMFMicroManager::GetMeanEntropy ( void  ) const

inline

Definition at line 86 of file G4StatMFMicroManager.hh.

{return _MeanEntropy; }

Here is the call graph for this function:

GetMeanMultiplicity()

G4double G4StatMFMicroManager::GetMeanMultiplicity ( void  ) const

inline

Definition at line 82 of file G4StatMFMicroManager.hh.

{return _MeanMultiplicity; }

GetMeanTemperature()

G4double G4StatMFMicroManager::GetMeanTemperature ( void  ) const

inline

Definition at line 84 of file G4StatMFMicroManager.hh.

{return _MeanTemperature; }

GetProbability()

G4double G4StatMFMicroManager::GetProbability ( void  ) const

inline

Definition at line 78 of file G4StatMFMicroManager.hh.

{return _WW;}

Here is the call graph for this function:

Here is the caller graph for this function:

Initialize()

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

private

Definition at line 82 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();
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

MakePartition()

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

private

Definition at line 139 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;
}

Here is the caller graph for this function:

Normalize()

void G4StatMFMicroManager::Normalize

(

G4double

Norm

)

Definition at line 159 of file G4StatMFMicroManager.cc.

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

Here is the caller graph for this function:

operator!=()

G4bool G4StatMFMicroManager::operator!=

(

const G4StatMFMicroManager &

right

)

const

Definition at line 57 of file G4StatMFMicroManager.cc.

{
    return true;
}

Here is the caller graph for this function:

operator=()

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

private

Definition at line 44 of file G4StatMFMicroManager.cc.

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

Here is the caller graph for this function:

operator==()

G4bool G4StatMFMicroManager::operator==

(

const G4StatMFMicroManager &

right

)

const

Definition at line 51 of file G4StatMFMicroManager.cc.

{
    return false;
}

Here is the caller graph for this function:

Member Data Documentation

_MeanEntropy

G4double G4StatMFMicroManager::_MeanEntropy

private

Definition at line 116 of file G4StatMFMicroManager.hh.

_MeanMultiplicity

G4double G4StatMFMicroManager::_MeanMultiplicity

private

Definition at line 112 of file G4StatMFMicroManager.hh.

_MeanTemperature

G4double G4StatMFMicroManager::_MeanTemperature

private

Definition at line 114 of file G4StatMFMicroManager.hh.

_Normalization

G4double G4StatMFMicroManager::_Normalization

private

Definition at line 110 of file G4StatMFMicroManager.hh.

_Partition

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

private

Definition at line 104 of file G4StatMFMicroManager.hh.

_WW

G4double G4StatMFMicroManager::_WW

private

Definition at line 108 of file G4StatMFMicroManager.hh.

---

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

• G4StatMFMicroManager.hh
• G4StatMFMicroManager.cc