ALIUnitsTable.cc File Reference

#include "Alignment/CocoaModel/interface/ALIUnitsTable.h"
#include "CLHEP/Units/GlobalSystemOfUnits.h"
#include <iomanip>
#include <cstdlib>
#include <cmath>

Go to the source code of this file.

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

Function Documentation

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

Definition at line 329 of file ALIUnitsTable.cc.

References funct::abs(), gather_cfg::cout, ALIUnitDefinition::GetUnitsTable(), ALIBestUnit::IndexOfCategory, dqmiolumiharvest::j, dqmdumpme::k, SiStripPI::max, ALIBestUnit::nbOfVals, particleFlowDisplacedVertex_cfi::ratio, unit(), and ALIBestUnit::Value.

                                                        {
  ALIUnitsTable& theUnitsTable = ALIUnitDefinition::GetUnitsTable();
  ALIUnitsContainer& List = theUnitsTable[a.IndexOfCategory]->GetUnitsList();
  ALIint len = theUnitsTable[a.IndexOfCategory]->GetSymbMxLen();

  ALIint ksup(-1), kinf(-1);
  ALIdouble umax(0.), umin(1.E12);
  ALIdouble rsup(1.E12), rinf(0.);

  //for a ThreeVector, choose the best unit for the biggest value
  ALIdouble value = std::max(std::max(std::abs(a.Value[0]), std::abs(a.Value[1])), std::abs(a.Value[2]));

  for (size_t k = 0; k < List.size(); k++) {
    ALIdouble unit = List[k]->GetValue();
    if (value == 1.E12) {
      if (unit > umax) {
        umax = unit;
        ksup = k;
      }
    } else if (value <= -1.E12) {
      if (unit < umin) {
        umin = unit;
        kinf = k;
      }
    }

    else {
      ALIdouble ratio = value / unit;
      if ((ratio >= 1.) && (ratio < rsup)) {
        rsup = ratio;
        ksup = k;
      }
      if ((ratio < 1.) && (ratio > rinf)) {
        rinf = ratio;
        kinf = k;
      }
    }
  }

  ALIint index = ksup;
  if (index == -1)
    index = kinf;
  if (index == -1)
    index = 0;

  for (ALIint j = 0; j < a.nbOfVals; j++) {
    flux << a.Value[j] / (List[index]->GetValue()) << " ";
  }

#ifdef ALIUSE_STD_NAMESPACE
  std::ios::fmtflags oldform = std::cout.flags();
#else
  //    ALIint oldform = std::cout.flags();
#endif

  flux.setf(std::ios::left, std::ios::adjustfield);
  flux << std::setw(len) << List[index]->GetSymbol();
  //??  flux.flags(oldform);

  return flux;
}