G4UnitsTable.cc File Reference

#include <iomanip>
#include <sstream>
#include "G4SystemOfUnits.hh"
#include "G4Threading.hh"
#include "G4UnitsTable.hh"

Go to the source code of this file.

Functions
std::ostream &	operator<< (std::ostream &flux, G4BestUnit a)

Function Documentation

operator<<()

std::ostream & operator<<	(	std::ostream &	flux,
		G4BestUnit	a
	)

Definition at line 601 of file G4UnitsTable.cc.

{
  G4UnitsTable& theUnitsTable = G4UnitDefinition::GetUnitsTable();
  G4UnitsContainer& List = theUnitsTable[a.IndexOfCategory]->GetUnitsList();
  G4int len              = theUnitsTable[a.IndexOfCategory]->GetSymbMxLen();
 
  G4int ksup(-1), kinf(-1);
  G4double umax(0.), umin(DBL_MAX);
  G4double rsup(DBL_MAX), rinf(0.);
 
  // for a ThreeVector, choose the best unit for the biggest value
  G4double value =
    std::max(std::max(std::fabs(a.Value[0]), std::fabs(a.Value[1])),
             std::fabs(a.Value[2]));
 
  //special treatement for Energy.
  if ((a.Category == "Energy") && (value == 0.)) {
    for (G4int j = 0; j < a.nbOfVals; ++j) {
       flux << a.Value[j] << " ";
    }
    std::ios::fmtflags oldform = flux.flags();
    flux.setf(std::ios::left, std::ios::adjustfield);
    flux << std::setw(len) << "eV";
    flux.flags(oldform);
    return flux;
  }      
         
  //here, value != 0.
  for(std::size_t k = 0; k < List.size(); ++k)
  {
    G4double unit = List[k]->GetValue();
    if(!(value != DBL_MAX))
    {
      if(unit > umax)
      {
        umax = unit;
        ksup = k;
      }
    }
    else if(value <= DBL_MIN)
    {
      if(unit < umin)
      {
        umin = unit;
        kinf = k;
      }
    }
    else
    {
      G4double ratio = value / unit;
      if((ratio >= 1.) && (ratio < rsup))
      {
        rsup = ratio;
        ksup = k;
      }
      if((ratio < 1.) && (ratio > rinf))
      {
        rinf = ratio;
        kinf = k;
      }
    }
  }
 
  G4int index = ksup;
  if(index == -1)
  {
    index = kinf;
  }
  if(index == -1)
  {
    index = 0;
  }
 
  for(G4int j = 0; j < a.nbOfVals; ++j)
  {
    flux << a.Value[j] / (List[index]->GetValue()) << " ";
  }
 
  std::ios::fmtflags oldform = flux.flags();
 
  flux.setf(std::ios::left, std::ios::adjustfield);
  flux << std::setw(len) << List[index]->GetSymbol();
  flux.flags(oldform);
 
  return flux;
}