#include <G4NuclWatcher.hh>

Public Member Functions
	G4NuclWatcher (G4int z, const std::vector< G4double > &expa, const std::vector< G4double > &expcs, const std::vector< G4double > &experr, G4bool check, G4bool nucl)

std::pair< G4double, G4double >	getAverageRatio () const

G4double	getChsq () const

std::pair< G4double, G4double >	getExpCs () const

std::pair< G4double, G4double >	getInuclCs () const

G4double	getLhood () const

G4double	getNmatched () const

G4bool	look_forNuclei () const

void	print ()

void	setInuclCs (G4double csec, G4int nev)

G4bool	to_check () const

void	watch (G4int a, G4int z)

	~G4NuclWatcher ()

Private Attributes
G4double	aver_lhood

G4double	aver_matched

G4double	aver_rat_err

G4double	average_ratio

G4bool	checkable

std::vector< G4double >	exper_as

std::vector< G4double >	exper_cs

std::vector< G4double >	exper_err

G4double	izotop_chsq

G4bool	nucleable

G4int	nuclz

std::vector< G4double >	simulated_as

std::vector< G4double >	simulated_cs

std::vector< G4double >	simulated_errors

std::vector< G4double >	simulated_prob

Detailed Description

Definition at line 40 of file G4NuclWatcher.hh.

Constructor & Destructor Documentation

◆ G4NuclWatcher()

G4NuclWatcher::G4NuclWatcher	(	G4int	z,
		const std::vector< G4double > &	expa,
		const std::vector< G4double > &	expcs,
		const std::vector< G4double > &	experr,
		G4bool	check,
		G4bool	nucl
	)

Definition at line 39 of file G4NuclWatcher.cc.

  : nuclz(z), izotop_chsq(0.), average_ratio(0.), aver_rat_err(0.),
    aver_lhood(0.), aver_matched(0.), exper_as(expa), exper_cs(expcs),
    exper_err(experr), checkable(check), nucleable(nucl) {}

◆ ~G4NuclWatcher()

G4NuclWatcher::~G4NuclWatcher ( )

inline

Definition at line 49 of file G4NuclWatcher.hh.

49{}

Member Function Documentation

◆ getAverageRatio()

std::pair< G4double, G4double > G4NuclWatcher::getAverageRatio ( ) const

inline

Definition at line 63 of file G4NuclWatcher.hh.

                                                      { 
    return std::pair<G4double, G4double>(average_ratio, aver_rat_err); 
  }

References aver_rat_err, and average_ratio.

◆ getChsq()

G4double G4NuclWatcher::getChsq ( ) const

inline

Definition at line 54 of file G4NuclWatcher.hh.

54{ return izotop_chsq; }

References izotop_chsq.

◆ getExpCs()

std::pair< G4double, G4double > G4NuclWatcher::getExpCs ( ) const

Definition at line 80 of file G4NuclWatcher.cc.

                                                          {
  G4double cs = 0.0;
  G4double err = 0.0;
  
  G4int experAsSize = exper_as.size();
  for(G4int iz = 0; iz < experAsSize; iz++) {
    cs += exper_cs[iz];
    err += exper_err[iz];
  }
  
  return std::pair<G4double, G4double>(cs, err);
}

References exper_as, exper_cs, and exper_err.

◆ getInuclCs()

std::pair< G4double, G4double > G4NuclWatcher::getInuclCs ( ) const

Definition at line 93 of file G4NuclWatcher.cc.

                                                            {
  G4double cs = 0.0;
  G4double err = 0.0;
  G4int  simulatedAsSize = simulated_as.size();
  for(G4int iz = 0; iz < simulatedAsSize; iz++) {
    cs += simulated_cs[iz];
    err += simulated_errors[iz];
  }
  
  return std::pair<G4double, G4double>(cs, err);
}

References simulated_as, simulated_cs, and simulated_errors.

◆ getLhood()

G4double G4NuclWatcher::getLhood ( ) const

inline

Definition at line 57 of file G4NuclWatcher.hh.

57{ return aver_lhood; }

References aver_lhood.

◆ getNmatched()

G4double G4NuclWatcher::getNmatched ( ) const

inline

Definition at line 58 of file G4NuclWatcher.hh.

58{ return aver_matched; }

References aver_matched.

◆ look_forNuclei()

G4bool G4NuclWatcher::look_forNuclei ( ) const

inline

Definition at line 56 of file G4NuclWatcher.hh.

56{ return nucleable; }

References nucleable.

◆ print()

void G4NuclWatcher::print ( )

Definition at line 105 of file G4NuclWatcher.cc.

                          {
  const G4double small = 0.001;
  
  G4cout << "\n ++++++++++++++++++++++++++++++++++++++++++++++++++++++++ "
     << "\n **** izotop Z **** " << nuclz << G4endl;
  
  izotop_chsq = 0.0;
  average_ratio = 0.0;
  aver_rat_err = 0.0;
  G4double exp_cs = 0.0;
  G4double exp_cs_err = 0.0;
  G4double inucl_cs = 0.0;
  G4double inucl_cs_err = 0.0;
  std::vector<G4bool> not_used(simulated_cs.size(), true);
  G4int nmatched = exper_as.size();
  G4int nused = simulated_cs.size();
  G4double lhood = 0.0;
  G4int experAsSize = exper_as.size();
 
  for (G4int iz = 0; iz < experAsSize; iz++) {
    G4double a = exper_as[iz];
    
    exp_cs += exper_cs[iz];
    exp_cs_err += exper_err[iz];
    
    G4bool found = false;
    G4int  simulatedAsSize = simulated_as.size();
    for (G4int i = 0; i<simulatedAsSize && !found; i++) {
      if (std::fabs(simulated_as[i] - a) < small) {
    G4double rat = simulated_cs[i] / exper_cs[iz];
    
    lhood += std::log10(rat) * std::log10(rat);
    
    G4double rat_err
      = std::sqrt(simulated_errors[i]*simulated_errors[i] +
              exper_err[iz]*exper_err[iz] * rat*rat) / exper_cs[iz];
    average_ratio += rat;
    aver_rat_err += rat_err; 
    
    G4cout << " A " << a << " exp.cs " << exper_cs[iz] << " err " 
           << exper_err[iz] << G4endl << " sim. cs " << simulated_cs[i]
           << " err " << simulated_errors[i] << G4endl
           << " ratio " << rat << " err " << rat_err << G4endl
           << " simulated production rate " << simulated_prob[i] << G4endl;
    
    not_used[i] = false;
    izotop_chsq += (rat - 1.0) * (rat - 1.0) / rat_err / rat_err; 
    found = true;
    nused--;
      }
    }
 
    if (found) nmatched--;
    else
      G4cout << " not found exper.: A " << a << " exp.cs " << exper_cs[iz] 
         << " err " << exper_err[iz] << G4endl;
  }
  
  G4cout << " not found in simulations " << nmatched << G4endl
     << " not found in exper: " << nused << G4endl;
 
  G4int  simulatedAsSize = simulated_as.size();
  for(G4int i = 0; i < simulatedAsSize; i++) {
    inucl_cs += simulated_cs[i];
    inucl_cs_err += simulated_errors[i];
    
    if(not_used[i]) 
      G4cout << " extra simul.: A " << simulated_as[i] << " sim. cs "
         << simulated_cs[i] << " err " << simulated_errors[i] << G4endl;
 
    G4cout << " simulated production rate " << simulated_prob[i] << G4endl;
  }
 
  G4int matched = exper_as.size() - nmatched;
  
  if (matched > 0) {
    aver_lhood = lhood;
    aver_matched = matched;    
    lhood = std::pow(10.0, std::sqrt(lhood/matched));
    
    G4cout << " matched " << matched << " CHSQ " << std::sqrt(izotop_chsq) / matched
       << G4endl
       << " raw chsq " << izotop_chsq << G4endl
       << " average ratio " << average_ratio / matched 
       << " err " << aver_rat_err / matched << G4endl 
       << " lhood " << lhood << G4endl;
  }
  else {
    izotop_chsq = 0.0;
    aver_lhood = 0.0;    
  } 
  
  G4cout << " exper. cs " << exp_cs << " err " << exp_cs_err << G4endl
     << " inucl. cs " << inucl_cs << " err " << inucl_cs_err << G4endl
     <<  " ++++++++++++++++++++++++++++++++++++++++++++++++++++++++ "
     << G4endl;
}

References aver_lhood, aver_matched, aver_rat_err, average_ratio, exper_as, exper_cs, exper_err, G4cout, G4endl, izotop_chsq, nuclz, simulated_as, simulated_cs, simulated_errors, and simulated_prob.

◆ setInuclCs()

void G4NuclWatcher::setInuclCs	(	G4double	csec,
		G4int	nev
	)

Definition at line 69 of file G4NuclWatcher.cc.

                                                       { 
  G4int  simulatedAsSize = simulated_as.size();
  for(G4int i = 0; i < simulatedAsSize ; i++) {
    double err = std::sqrt(simulated_cs[i]) / simulated_cs[i];
    
    simulated_prob.push_back(simulated_cs[i] / nev);
    simulated_cs[i] *= csec / nev;
    simulated_errors.push_back(simulated_cs[i] * err);    
  }
}

References simulated_as, simulated_cs, simulated_errors, and simulated_prob.

◆ to_check()

G4bool G4NuclWatcher::to_check ( ) const

inline

Definition at line 55 of file G4NuclWatcher.hh.

55{ return checkable; }

References checkable.

◆ watch()

void G4NuclWatcher::watch	(	G4int	a,
		G4int	z
	)

Definition at line 49 of file G4NuclWatcher.cc.

                                          {
  const G4double small = 0.001;
  
  if (std::abs(z-nuclz) >= small) return;
 
  G4bool here = false;      // Increment specified nucleus count
  G4int  simulatedAsSize = simulated_as.size();
  for (G4int i = 0; i<simulatedAsSize && !here; i++) {
    if (std::abs(simulated_as[i] - a) < small) {
      simulated_cs[i] += 1.0;
      here = true;      // Terminates loop
    }
  }
 
  if (!here) {          // Nothing found, so add new entry
    simulated_as.push_back(a); 
    simulated_cs.push_back(1.0); 
  }
}

References nuclz, simulated_as, and simulated_cs.

Field Documentation

◆ aver_lhood

G4double G4NuclWatcher::aver_lhood

private

Definition at line 74 of file G4NuclWatcher.hh.

Referenced by getLhood(), and print().

◆ aver_matched

G4double G4NuclWatcher::aver_matched

private

Definition at line 75 of file G4NuclWatcher.hh.

Referenced by getNmatched(), and print().

◆ aver_rat_err

G4double G4NuclWatcher::aver_rat_err

private

Definition at line 73 of file G4NuclWatcher.hh.

Referenced by getAverageRatio(), and print().

◆ average_ratio

G4double G4NuclWatcher::average_ratio

private

Definition at line 72 of file G4NuclWatcher.hh.

Referenced by getAverageRatio(), and print().

◆ checkable

G4bool G4NuclWatcher::checkable

private

Definition at line 83 of file G4NuclWatcher.hh.

Referenced by to_check().

◆ exper_as

std::vector<G4double> G4NuclWatcher::exper_as

private

Definition at line 76 of file G4NuclWatcher.hh.

Referenced by getExpCs(), and print().

◆ exper_cs

std::vector<G4double> G4NuclWatcher::exper_cs

private

Definition at line 77 of file G4NuclWatcher.hh.

Referenced by getExpCs(), and print().

◆ exper_err

std::vector<G4double> G4NuclWatcher::exper_err

private

Definition at line 78 of file G4NuclWatcher.hh.

Referenced by getExpCs(), and print().

◆ izotop_chsq

G4double G4NuclWatcher::izotop_chsq

private

Definition at line 71 of file G4NuclWatcher.hh.

Referenced by getChsq(), and print().

◆ nucleable

G4bool G4NuclWatcher::nucleable

private

Definition at line 84 of file G4NuclWatcher.hh.

Referenced by look_forNuclei().

◆ nuclz

G4int G4NuclWatcher::nuclz

private

Definition at line 70 of file G4NuclWatcher.hh.

Referenced by print(), and watch().

◆ simulated_as

std::vector<G4double> G4NuclWatcher::simulated_as

private

Definition at line 79 of file G4NuclWatcher.hh.

Referenced by getInuclCs(), print(), setInuclCs(), and watch().

◆ simulated_cs

std::vector<G4double> G4NuclWatcher::simulated_cs

private

Definition at line 80 of file G4NuclWatcher.hh.

Referenced by getInuclCs(), print(), setInuclCs(), and watch().

◆ simulated_errors

std::vector<G4double> G4NuclWatcher::simulated_errors

private

Definition at line 81 of file G4NuclWatcher.hh.

Referenced by getInuclCs(), print(), and setInuclCs().

◆ simulated_prob

std::vector<G4double> G4NuclWatcher::simulated_prob

private

Definition at line 82 of file G4NuclWatcher.hh.

Referenced by print(), and setInuclCs().

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

source/processes/hadronic/models/cascade/cascade/include/G4NuclWatcher.hh
source/processes/hadronic/models/cascade/cascade/src/G4NuclWatcher.cc

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4NuclWatcher()

◆ ~G4NuclWatcher()

Member Function Documentation

◆ getAverageRatio()

◆ getChsq()

◆ getExpCs()

◆ getInuclCs()

◆ getLhood()

◆ getNmatched()

◆ look_forNuclei()

◆ print()

◆ setInuclCs()

◆ to_check()

◆ watch()

Field Documentation

◆ aver_lhood

◆ aver_matched

◆ aver_rat_err

◆ average_ratio

◆ checkable

◆ exper_as

◆ exper_cs

◆ exper_err

◆ izotop_chsq

◆ nucleable

◆ nuclz

◆ simulated_as

◆ simulated_cs

◆ simulated_errors

◆ simulated_prob