#include <G4FissionStore.hh>

Public Member Functions
void	addConfig (G4double a, G4double z, G4double ez, G4double ek, G4double ev)

void	clear ()

	G4FissionStore ()

G4FissionConfiguration	generateConfiguration (G4double amax, G4double rand) const

void	setVerboseLevel (G4int verbose=1)

size_t	size () const

Private Attributes
std::vector< G4double >	configProbs

std::vector< G4FissionConfiguration >	configurations

G4int	verboseLevel

Detailed Description

Definition at line 36 of file G4FissionStore.hh.

Constructor & Destructor Documentation

◆ G4FissionStore()

G4FissionStore::G4FissionStore ( )

Definition at line 38 of file G4FissionStore.cc.

                               : verboseLevel(0) {
  if (verboseLevel > 1) 
    G4cout << " >>> G4FissionStore::G4FissionStore" << G4endl;
}

References G4cout, G4endl, and verboseLevel.

Member Function Documentation

◆ addConfig()

void G4FissionStore::addConfig	(	G4double	a,
		G4double	z,
		G4double	ez,
		G4double	ek,
		G4double	ev
	)

Definition at line 43 of file G4FissionStore.cc.

                                             {
  G4FissionConfiguration config(a, z, ez, ek, ev);
  configurations.push_back(config);
  if (verboseLevel > 2) G4cout << config << G4endl;
}

References configurations, G4cout, G4endl, and verboseLevel.

Referenced by G4Fissioner::deExcite().

◆ clear()

void G4FissionStore::clear ( )

inline

Definition at line 44 of file G4FissionStore.hh.

44{ configurations.clear(); }

References configurations.

Referenced by G4Fissioner::deExcite().

◆ generateConfiguration()

G4FissionConfiguration G4FissionStore::generateConfiguration	(	G4double	amax,
		G4double	rand
	)		const

Definition at line 50 of file G4FissionStore.cc.

                                                      {
  if (verboseLevel > 1)
    G4cout << " >>> G4FissionStore::generateConfiguration" << G4endl;
  
  const G4double small = -30.0;
 
  G4double totProb = 0.0;
  configProbs.resize(size(),0.);
 
  if (verboseLevel > 3)
    G4cout << " amax " << amax << " ic " << size() << G4endl;
 
  for (size_t i = 0; i < size(); i++) {
    G4double ez = configurations[i].ezet;
    G4double pr = ez - amax;
 
    if (pr < small) pr = small;
    pr = G4Exp(pr);
    if (verboseLevel > 2) {
      G4cout << configurations[i] << "\n probability " << pr << G4endl; 
    }
    totProb += pr;
    configProbs[i] = totProb;  
  };
 
  G4double st = totProb * rand;
 
  size_t igen = 0;
  /* Loop checking 08.06.2015 MHK */
  while (configProbs[igen] <= st && igen < size()) igen++;
 
  if (verboseLevel > 3) G4cout << " igen " << igen << G4endl;
 
  return configurations[igen];
}           