G4ConvergenceTester Class Reference

#include <G4ConvergenceTester.hh>

Public Member Functions
	G4ConvergenceTester (G4String theName="NONAME")
	~G4ConvergenceTester ()
	G4ConvergenceTester (G4double)
void	AddScore (G4double)
G4ConvergenceTester &	operator+= (G4double val)
void	ShowHistory (std::ostream &out=G4cout)
void	ShowResult (std::ostream &out=G4cout)
G4double	GetValueOfMinimizingFunction (std::vector< G4double > x)
void	ComputeStatistics ()
G4double	GetMean ()
G4double	GetStandardDeviation ()
G4double	GetVariance ()
G4double	GetR ()
G4double	GetEfficiency ()
G4double	GetR2eff ()
G4double	GetR2int ()
G4double	GetShift ()
G4double	GetVOV ()
G4double	GetFOM ()

Detailed Description

Definition at line 54 of file G4ConvergenceTester.hh.

Constructor & Destructor Documentation

G4ConvergenceTester() [1/2]

G4ConvergenceTester::G4ConvergenceTester

(

G4String

theName = "NONAME"

)

Definition at line 34 of file G4ConvergenceTester.cc.

  : name(theName)
{
  nonzero_histories.clear();
  largest_scores.clear();
  largest_scores.push_back(0.0);
 
  history_grid.resize(noBinOfHistory, 0);
  mean_history.resize(noBinOfHistory, 0.0);
  var_history.resize(noBinOfHistory, 0.0);
  sd_history.resize(noBinOfHistory, 0.0);
  r_history.resize(noBinOfHistory, 0.0);
  vov_history.resize(noBinOfHistory, 0.0);
  fom_history.resize(noBinOfHistory, 0.0);
  shift_history.resize(noBinOfHistory, 0.0);
  e_history.resize(noBinOfHistory, 0.0);
  r2eff_history.resize(noBinOfHistory, 0.0);
  r2int_history.resize(noBinOfHistory, 0.0);
 
  timer = new G4Timer();
  timer->Start();
  cpu_time.clear();
  cpu_time.push_back(0.0);
}

~G4ConvergenceTester()

G4ConvergenceTester::~G4ConvergenceTester

(

)

Definition at line 59 of file G4ConvergenceTester.cc.

{ delete timer; }

G4ConvergenceTester() [2/2]

G4ConvergenceTester::G4ConvergenceTester

(

G4double

)

Member Function Documentation

AddScore()

void G4ConvergenceTester::AddScore

(

G4double

x

)

Definition at line 61 of file G4ConvergenceTester.cc.

{
  // G4cout << x << G4endl;
 
  timer->Stop();
  cpu_time.push_back(timer->GetSystemElapsed() + timer->GetUserElapsed());
 
  if(x < 0.0)
  {
    G4cout << "Warning: G4convergenceTester expects zero or positive number as "
              "inputs, but received a negative number."
           << G4endl;
  }
 
  if(x == 0.0)
  {
  }
  else
  {
    nonzero_histories.insert(std::pair<G4int, G4double>(n, x));
    if(x > largest_scores.back())
    {
      //        Following serch should become faster if begin from bottom.
      std::vector<G4double>::iterator it;
      for(it = largest_scores.begin(); it != largest_scores.end(); it++)
      {
        if(x > *it)
        {
          largest_scores.insert(it, x);
          break;
        }
      }
 
      if(largest_scores.size() > 201)
      {
        largest_scores.pop_back();
      }
    }
    sum += x;
  }
 
  // Data has been added so statistics have not been updated to new values
  statsAreUpdated = false;
  n++;
  return;
}

Referenced by operator+=().

ComputeStatistics()

void G4ConvergenceTester::ComputeStatistics ( )

inline

Definition at line 79 of file G4ConvergenceTester.hh.

{ calStat(); }

GetEfficiency()

G4double G4ConvergenceTester::GetEfficiency ( )

inline

Definition at line 104 of file G4ConvergenceTester.hh.

  {
    CheckIsUpdated();
    return efficiency;
  }

GetFOM()

G4double G4ConvergenceTester::GetFOM ( )

inline

Definition at line 129 of file G4ConvergenceTester.hh.

  {
    CheckIsUpdated();
    return fom;
  }

GetMean()

G4double G4ConvergenceTester::GetMean ( )

inline

Definition at line 84 of file G4ConvergenceTester.hh.

  {
    CheckIsUpdated();
    return mean;
  }

GetR()

G4double G4ConvergenceTester::GetR ( )

inline

Definition at line 99 of file G4ConvergenceTester.hh.

  {
    CheckIsUpdated();
    return r;
  }

GetR2eff()

G4double G4ConvergenceTester::GetR2eff ( )

inline

Definition at line 109 of file G4ConvergenceTester.hh.

  {
    CheckIsUpdated();
    return r2eff;
  }

GetR2int()

G4double G4ConvergenceTester::GetR2int ( )

inline

Definition at line 114 of file G4ConvergenceTester.hh.

  {
    CheckIsUpdated();
    return r2int;
  }

GetShift()

G4double G4ConvergenceTester::GetShift ( )

inline

Definition at line 119 of file G4ConvergenceTester.hh.

  {
    CheckIsUpdated();
    return shift;
  }

GetStandardDeviation()

G4double G4ConvergenceTester::GetStandardDeviation ( )

inline

Definition at line 89 of file G4ConvergenceTester.hh.

  {
    CheckIsUpdated();
    return sd;
  }

GetValueOfMinimizingFunction()

G4double G4ConvergenceTester::GetValueOfMinimizingFunction ( std::vector< G4double >  x )

inline

Definition at line 73 of file G4ConvergenceTester.hh.

  {
    return slope_fitting_function(x);
  }

GetVariance()

G4double G4ConvergenceTester::GetVariance ( )

inline

Definition at line 94 of file G4ConvergenceTester.hh.

  {
    CheckIsUpdated();
    return var;
  }

GetVOV()

G4double G4ConvergenceTester::GetVOV ( )

inline

Definition at line 124 of file G4ConvergenceTester.hh.

  {
    CheckIsUpdated();
    return vov;
  }

operator+=()

G4ConvergenceTester & G4ConvergenceTester::operator+= ( G4double  val )

inline

Definition at line 63 of file G4ConvergenceTester.hh.

  {
    this->AddScore(val);
    return *this;
  }

ShowHistory()

void G4ConvergenceTester::ShowHistory

(

std::ostream &

out = G4cout

)

Definition at line 461 of file G4ConvergenceTester.cc.

{
  if(!showHistory)
  {
    out << "Number of events of this run is too small to show history."
        << G4endl;
    return;
  }
 
  out << std::setprecision(6);
 
  out << G4endl;
  out << "G4ConvergenceTester Output History of " << name << G4endl;
  out << "i/" << noBinOfHistory << " till_ith      mean" << std::setw(13)
      << "var" << std::setw(13) << "sd" << std::setw(13) << "r" << std::setw(13)
      << "vov" << std::setw(13) << "fom" << std::setw(13) << "shift"
      << std::setw(13) << "e" << std::setw(13) << "r2eff" << std::setw(13)
      << "r2int" << G4endl;
  for(G4int i = 1; i <= noBinOfHistory; i++)
  {
    out << std::setw(4) << i << " " << std::setw(5) << history_grid[i - 1]
        << std::setw(13) << mean_history[i - 1] << std::setw(13)
        << var_history[i - 1] << std::setw(13) << sd_history[i - 1]
        << std::setw(13) << r_history[i - 1] << std::setw(13)
        << vov_history[i - 1] << std::setw(13) << fom_history[i - 1]
        << std::setw(13) << shift_history[i - 1] << std::setw(13)
        << e_history[i - 1] << std::setw(13) << r2eff_history[i - 1]
        << std::setw(13) << r2int_history[i - 1] << G4endl;
  }
}

ShowResult()

void G4ConvergenceTester::ShowResult

(

std::ostream &

out = G4cout

)

Definition at line 353 of file G4ConvergenceTester.cc.

{
  // if data has been added since the last computation of the statistical values
  // (not statsAreUpdated) call calStat to recompute the statistical values
  if(!statsAreUpdated)
  {
    calStat();
  }
 
  out << std::setprecision(6);
 
  out << G4endl;
  out << "G4ConvergenceTester Output Result of " << name << G4endl;
  out << std::setw(20) << "EFFICIENCY = " << std::setw(13) << efficiency
      << G4endl;
  out << std::setw(20) << "MEAN = " << std::setw(13) << mean << G4endl;
  out << std::setw(20) << "VAR = " << std::setw(13) << var << G4endl;
  out << std::setw(20) << "SD = " << std::setw(13) << sd << G4endl;
  out << std::setw(20) << "R = " << std::setw(13) << r << G4endl;
  out << std::setw(20) << "SHIFT = " << std::setw(13) << shift << G4endl;
  out << std::setw(20) << "VOV = " << std::setw(13) << vov << G4endl;
  out << std::setw(20) << "FOM = " << std::setw(13) << fom << G4endl;
 
  out << std::setw(20) << "THE LARGEST SCORE = " << std::setw(13) << largest
      << " and it happened at " << largest_score_happened << "th event"
      << G4endl;
  if(mean != 0)
  {
    out << std::setw(20) << "Affected Mean = " << std::setw(13) << mean_1
        << " and its ratio to original is " << mean_1 / mean << G4endl;
  }
  else
  {
    out << std::setw(20) << "Affected Mean = " << std::setw(13) << mean_1
        << G4endl;
  }
  if(var != 0)
  {
    out << std::setw(20) << "Affected VAR = " << std::setw(13) << var_1
        << " and its ratio to original is " << var_1 / var << G4endl;
  }
  else
  {
    out << std::setw(20) << "Affected VAR = " << std::setw(13) << var_1
        << G4endl;
  }
  if(r != 0)
  {
    out << std::setw(20) << "Affected R = " << std::setw(13) << r_1
        << " and its ratio to original is " << r_1 / r << G4endl;
  }
  else
  {
    out << std::setw(20) << "Affected R = " << std::setw(13) << r_1 << G4endl;
  }
  if(shift != 0)
  {
    out << std::setw(20) << "Affected SHIFT = " << std::setw(13) << shift_1
        << " and its ratio to original is " << shift_1 / shift << G4endl;
  }
  else
  {
    out << std::setw(20) << "Affected SHIFT = " << std::setw(13) << shift_1
        << G4endl;
  }
  if(fom != 0)
  {
    out << std::setw(20) << "Affected FOM = " << std::setw(13) << fom_1
        << " and its ratio to original is " << fom_1 / fom << G4endl;
  }
  else
  {
    out << std::setw(20) << "Affected FOM = " << std::setw(13) << fom_1
        << G4endl;
  }
 
  if(!showHistory)
  {
    out << "Number of events of this run is too small to do convergence tests."
        << G4endl;
    return;
  }
 
  check_stat_history(out);
 
  // check SLOPE and output result
  if(calcSLOPE)
  {
    if(slope >= 3)
    {
      noPass++;
      out << "SLOPE is large enough" << G4endl;
    }
    else
    {
      out << "SLOPE is not large enough" << G4endl;
    }
  }
  else
  {
    out << "Number of non zero history too small to calculate SLOPE" << G4endl;
  }
 
  out << "This result passes " << noPass << " / " << noTotal
      << " Convergence Test." << G4endl;
  out << G4endl;
}

---

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

• G4ConvergenceTester.hh
• G4ConvergenceTester.cc