G4ConvergenceTester Class Reference

#include <G4ConvergenceTester.hh>

Public Member Functions
	G4ConvergenceTester ()
	~G4ConvergenceTester ()
	G4ConvergenceTester (G4double)
void	AddScore (G4double)
void	ShowHistory ()
void	ShowResult ()
G4double	GetValueOfMinimizingFunction (std::vector< G4double > x)

Detailed Description

Definition at line 58 of file G4ConvergenceTester.hh.

Constructor & Destructor Documentation

G4ConvergenceTester() [1/2]

G4ConvergenceTester::G4ConvergenceTester

(

)

Definition at line 45 of file G4ConvergenceTester.cc.

 : n(0), sum(0.), mean(0.), var(0.), sd(0.), r(0.), efficiency(0.),
   r2eff(0.), r2int(0.), shift(0.), vov(0.), fom(0.), largest(0.),
   largest_score_happened(0), mean_1(0.), var_1(0.), sd_1(0.), r_1(0.),
   shift_1(0.), vov_1(0.), fom_1(0.), noBinOfHistory(16), slope(0.),
   noBinOfPDF(10), minimizer(0), noPass(0), noTotal(8)
{
   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 76 of file G4ConvergenceTester.cc.

{
   delete timer;
}

G4ConvergenceTester() [2/2]

G4ConvergenceTester::G4ConvergenceTester

(

G4double

)

Member Function Documentation

AddScore()

void G4ConvergenceTester::AddScore

(

G4double

x

)

Definition at line 83 of file G4ConvergenceTester.cc.

{ 
 
   //G4cout << x << G4endl;
 
   timer->Stop();
   cpu_time.push_back( timer->GetSystemElapsed() + timer->GetUserElapsed() );
 
   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();
          } 
          //G4cout << largest_scores.size() << " " << largest_scores.front() << " " << largest_scores.back() << G4endl;
       }
       sum += x; 
   }
 
   n++;
   return; 
}

GetValueOfMinimizingFunction()

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

inline

Definition at line 73 of file G4ConvergenceTester.hh.

             { return slope_fitting_function( x ); }

ShowHistory()

void G4ConvergenceTester::ShowHistory

(

)

Definition at line 377 of file G4ConvergenceTester.cc.

{
   G4cout << "i/" << noBinOfHistory << " till_ith  mean  var  sd  r  vov  fom  shift  e  r2eff  r2int" << G4endl;
   for ( G4int i = 1 ; i <=  noBinOfHistory  ; i++ )
   {
      G4cout << i << " " 
             << history_grid [ i-1 ] << " "
             << mean_history [ i-1 ] << " "  
             << var_history [ i-1 ] << " "  
             << sd_history [ i-1 ] << " "  
             << r_history [ i-1 ] << " "  
             << vov_history [ i-1 ] << " "  
             << fom_history [ i-1 ] << " " 
             << shift_history [ i-1 ] << " "  
             << e_history [ i-1 ] << " " 
             << r2eff_history [ i-1 ] << " " 
             << r2int_history [ i-1 ] << " " 
             << G4endl;
   }
}

ShowResult()

void G4ConvergenceTester::ShowResult

(

)

Definition at line 339 of file G4ConvergenceTester.cc.

{
   calStat();
 
   G4cout << "EFFICIENCY = " << efficiency << G4endl;
   G4cout << "MEAN = " << mean << G4endl;
   G4cout << "VAR = " << var << G4endl;
   G4cout << "SD = " << sd << G4endl;
   G4cout << "R = "<< r << G4endl;
   G4cout << "SHIFT = "<< shift << G4endl;
   G4cout << "VOV = "<< vov << G4endl;
   G4cout << "FOM = "<< fom << G4endl;
 
   G4cout << "THE LARGEST SCORE = " << largest << " and it happend at " << largest_score_happened << "th event" << G4endl;
   G4cout << "Affected Mean = " << mean_1 << " and its ratio to orignal is " << mean_1/mean << G4endl;
   G4cout << "Affected VAR = " << var_1 << " and its ratio to orignal is " << var_1/var << G4endl;
   G4cout << "Affected R = " << r_1 << " and its ratio to orignal is " << r_1/r << G4endl;
   G4cout << "Affected SHIFT = " << shift_1 << " and its ratio to orignal is " << shift_1/shift << G4endl;
   G4cout << "Affected FOM = " << fom_1 << " and its ratio to orignal is " << fom_1/fom << G4endl;
 
   check_stat_history();
 
// check SLOPE and output result
   if ( slope >= 3 )
   {    
      noPass++;
      G4cout << "SLOPE is large enough" << G4endl; 
   }
   else
   {
      G4cout << "SLOPE is not large enough" << G4endl; 
   }
 
   G4cout << "This result passes " << noPass << " / "<< noTotal << " Convergence Test." << G4endl; 
   G4cout << G4endl;
 
}

---

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

• G4ConvergenceTester.hh
• G4ConvergenceTester.cc