ALIUtils Class Reference

#include <ALIUtils.h>

Public Member Functions
	ALIUtils ()
	~ALIUtils ()

Static Public Member Functions
static double	addPii (double val)
static ALIdouble	AngleSigmaDimensionFactor ()
static ALIdouble	AngleValueDimensionFactor ()
static double	approxTo0 (double val)
static ALIdouble	CalculateAngleDimensionFactorFromInt (ALIint ad)
static ALIdouble	CalculateAngleDimensionFactorFromString (ALIstring dimstr)
static ALIdouble	CalculateLengthDimensionFactorFromInt (ALIint ad)
static ALIdouble	CalculateLengthDimensionFactorFromString (ALIstring dimstr)
static std::string	changeName (const std::string &oldName, const std::string &subsstr1, const std::string &subsstr2)
static int	checkMatrixEquations (double angleX, double angleY, double angleZ, CLHEP::HepRotation *rot)
static double	diff2pi (double ang1, double ang2)
static void	dump3v (const CLHEP::Hep3Vector &vec, const std::string &msg)
static void	dumpDimensions (std::ofstream &fout)
static void	dumprm (const CLHEP::HepRotation &rm, const std::string &msg, std::ostream &out=std::cout)
static void	dumpVS (const std::vector< ALIstring > &wl, const std::string &msg, std::ostream &outs=std::cout)
	dumps a vector of strings with a message to outs More...
static bool	eq2ang (double ang1, double ang2)
static bool	getBool (const ALIstring &str)
	Convert a bool to an integer, checking that it is really a bool. More...
static ALIdouble	getDimensionValue (const ALIstring &dim, const ALIstring &dimType)
static ALIbool	getFirstTime ()
static double	getFloat (const ALIstring &str)
	Convert a string to an float, checking that it is really a number. More...
static int	getInt (const ALIstring &str)
	Convert a string to an integer, checking that it is really an integer. More...
static ALIdouble	getMaximumDeviationDerivative ()
static std::vector< double >	getRotationAnglesFromMatrix (CLHEP::HepRotation &rmLocal, double origAngleX, double origAngleY, double origAngleZ)
static int	IsNumber (const ALIstring &str)
static ALIdouble	LengthSigmaDimensionFactor ()
static ALIdouble	LengthValueDimensionFactor ()
static ALIdouble	OutputAngleSigmaDimensionFactor ()
static ALIdouble	OutputAngleValueDimensionFactor ()
static ALIdouble	OutputLengthSigmaDimensionFactor ()
static ALIdouble	OutputLengthValueDimensionFactor ()
static void	set_time_now (time_t now)
static void	SetAngleDimensionFactors ()
static void	setDebugVerbosity (ALIint val)
static void	setFirstTime (ALIbool val)
static void	SetLengthDimensionFactors ()
static void	setMaximumDeviationDerivative (ALIdouble val)
static void	SetOutputAngleDimensionFactors ()
static void	SetOutputLengthDimensionFactors ()
static void	setReportVerbosity (ALIint val)
static ALIstring	subQuotes (const ALIstring &str)
static time_t	time_now ()
static ALIdouble	val0 (ALIdouble val)

Static Public Attributes
static ALIint	debug = -1
static ALIdouble	deg = 0.017453293
static ALIint	report = 1

Static Private Attributes
static ALIdouble	_AngleSigmaDimensionFactor = 1.
static ALIdouble	_AngleValueDimensionFactor = 1.
static ALIdouble	_LengthSigmaDimensionFactor = 1.
static ALIdouble	_LengthValueDimensionFactor = 1.
static ALIdouble	_OutputAngleSigmaDimensionFactor = 1.
static ALIdouble	_OutputAngleValueDimensionFactor = 1.
static ALIdouble	_OutputLengthSigmaDimensionFactor = 1.
static ALIdouble	_OutputLengthValueDimensionFactor = 1.
static time_t	_time_now
static ALIbool	firstTime = false
static ALIdouble	maximum_deviation_derivative = 1.E-6

Detailed Description

Definition at line 23 of file ALIUtils.h.

Constructor & Destructor Documentation

ALIUtils::ALIUtils ( )

inline

Definition at line 25 of file ALIUtils.h.

{};

ALIUtils::~ALIUtils ( )

inline

Definition at line 26 of file ALIUtils.h.

References gather_cfg::cout, dump3v(), dumprm(), IsNumber(), mps_check::msg, MillePedeFileConverter_cfg::out, submit::rm, str, and AlCaHLTBitMon_QueryRunRegistry::string.

{};

Member Function Documentation

double ALIUtils::addPii ( double  val )

static

Definition at line 751 of file ALIUtils.cc.

References M_PI, and heppy_batch::val.

Referenced by getRotationAnglesFromMatrix(), and setMaximumDeviationDerivative().

                                  {
  if (val < M_PI) {
    val += M_PI;
  } else {
    val -= M_PI;
  }

  return val;
}

static ALIdouble ALIUtils::AngleSigmaDimensionFactor ( )

inlinestatic

Definition at line 65 of file ALIUtils.h.

References _AngleSigmaDimensionFactor.

Referenced by EntryAngle::SigmaDimensionFactor().

{ return _AngleSigmaDimensionFactor; }

static ALIdouble ALIUtils::AngleValueDimensionFactor ( )

inlinestatic

Definition at line 64 of file ALIUtils.h.

References _AngleValueDimensionFactor.

Referenced by EntryMgr::readEntryFromReportOut(), OpticalObject::transformCylindrical2Cartesian(), and EntryAngle::ValueDimensionFactor().

{ return _AngleValueDimensionFactor; }

double ALIUtils::approxTo0 ( double  val )

static

Definition at line 743 of file ALIUtils.cc.

References common_cff::precision, and heppy_batch::val.

Referenced by getRotationAnglesFromMatrix(), and setMaximumDeviationDerivative().

                                     {
  double precision = 1.e-9;
  if (fabs(val) < precision)
    val = 0;
  return val;
}

ALIdouble ALIUtils::CalculateAngleDimensionFactorFromInt ( ALIint  ad )

static

Definition at line 271 of file ALIUtils.cc.

References CalculateAngleDimensionFactorFromString(), beam_dqm_sourceclient-live_cfg::cerr, and beamvalidation::exit().

Referenced by set_time_now(), SetAngleDimensionFactors(), and SetOutputAngleDimensionFactors().

                                                                  {
  ALIdouble valsig;
  switch (ad) {
    case 0:  //----- radians
      valsig = CalculateAngleDimensionFactorFromString("rad");
      break;
    case 1:  //----- miliradians
      valsig = CalculateAngleDimensionFactorFromString("mrad");
      break;
    case 2:  //----- microradians
      valsig = CalculateAngleDimensionFactorFromString("murad");
      break;
    case 3:  //----- degrees
      valsig = CalculateAngleDimensionFactorFromString("deg");
      break;
    case 4:  //----- grads
      valsig = CalculateAngleDimensionFactorFromString("grad");
      break;
    default:
      std::cerr << "!!! UNKNOWN DIMENSION SCALING " << ad << std::endl << "VALUE MUST BE BETWEEN 0 AND 3 " << std::endl;
      exit(1);
  }

  // use microradinas instead of radians
  //-  valsig *= 1000000.;

  return valsig;
}

ALIdouble ALIUtils::CalculateAngleDimensionFactorFromString ( ALIstring  dimstr )

static

Definition at line 218 of file ALIUtils.cc.

References beam_dqm_sourceclient-live_cfg::cerr, beamvalidation::exit(), and M_PI.

Referenced by CalculateAngleDimensionFactorFromInt(), FittedEntriesReader::FittedEntriesReader(), and set_time_now().

                                                                            {
  ALIdouble valsig;
  if (dimstr == "rad") {
    valsig = 1.;
  } else if (dimstr == "mrad") {
    valsig = 1.E-3;
  } else if (dimstr == "murad") {
    valsig = 1.E-6;
  } else if (dimstr == "deg") {
    valsig = M_PI / 180.;
  } else if (dimstr == "grad") {
    valsig = M_PI / 200.;
  } else {
    std::cerr << "!!! UNKNOWN DIMENSION SCALING: " << dimstr << std::endl
              << "VALUE MUST BE AN ANGLE DIMENSION " << std::endl;
    exit(1);
  }

  return valsig;
}

ALIdouble ALIUtils::CalculateLengthDimensionFactorFromInt ( ALIint  ad )

static

Definition at line 242 of file ALIUtils.cc.

References CalculateLengthDimensionFactorFromString(), beam_dqm_sourceclient-live_cfg::cerr, and beamvalidation::exit().

Referenced by set_time_now(), SetLengthDimensionFactors(), and SetOutputLengthDimensionFactors().

                                                                   {
  ALIdouble valsig;
  switch (ad) {
    case 0:  //----- metres
      valsig = CalculateLengthDimensionFactorFromString("m");
      break;
    case 1:  //----- milimetres
      valsig = CalculateLengthDimensionFactorFromString("mm");
      break;
    case 2:  //----- micrometres
      valsig = CalculateLengthDimensionFactorFromString("mum");
      break;
    case 3:  //----- centimetres
      valsig = CalculateLengthDimensionFactorFromString("cm");
      break;
    default:
      std::cerr << "!!! UNKNOWN DIMENSION SCALING " << ad << std::endl << "VALUE MUST BE BETWEEN 0 AND 3 " << std::endl;
      exit(1);
  }

  // use microradinas instead of radians
  //-  valsig *= 1000000.;

  return valsig;
}

ALIdouble ALIUtils::CalculateLengthDimensionFactorFromString ( ALIstring  dimstr )

static

Definition at line 179 of file ALIUtils.cc.

References beam_dqm_sourceclient-live_cfg::cerr, and beamvalidation::exit().

Referenced by CalculateLengthDimensionFactorFromInt(), FittedEntriesReader::FittedEntriesReader(), and set_time_now().

                                                                             {
  ALIdouble valsig = 1.;
  ALIstring internalDim = "m";
  if (internalDim == "m") {
    if (dimstr == "m") {
      valsig = 1.;
    } else if (dimstr == "mm") {
      valsig = 1.E-3;
    } else if (dimstr == "mum") {
      valsig = 1.E-6;
    } else if (dimstr == "cm") {
      valsig = 1.E-2;
    } else {
      std::cerr << "!!! UNKNOWN DIMENSION SCALING " << dimstr << std::endl
                << "VALUE MUST BE BETWEEN 0 AND 3 " << std::endl;
      exit(1);
    }
  } else if (internalDim == "mm") {
    if (dimstr == "m") {
      valsig = 1.E3;
    } else if (dimstr == "mm") {
      valsig = 1.;
    } else if (dimstr == "mum") {
      valsig = 1.E-3;
    } else if (dimstr == "cm") {
      valsig = 1.E+1;
    } else {
      std::cerr << "!!! UNKNOWN DIMENSION SCALING: " << dimstr << std::endl
                << "VALUE MUST BE A LENGTH DIMENSION " << std::endl;
      exit(1);
    }
  }

  return valsig;
}

std::string ALIUtils::changeName ( const std::string &  oldName, const std::string &  subsstr1, const std::string &  subsstr2  )

static

Definition at line 535 of file ALIUtils.cc.

References mps_merge::newName, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by Measurement::copyMeas(), Model::copyMeasurements(), and val0().

                                                                                                             {
  std::string newName = oldName;
  int il = oldName.find(subsstr1);
  //  std::cout << " il " << il << " oldname " << oldName << " " << subsstr1 << std::endl;
  while (il >= 0) {
    newName = newName.substr(0, il) + subsstr2 + newName.substr(il + subsstr1.length(), newName.length());
    //    std::cout << " dnewName " << newName << " " << newName.substr( 0, il ) << " " << subsstr2 << " " << newName.substr( il+subsstr1.length(), newName.length() ) << std::endl;
    il = oldName.find(subsstr1, il + 1);
  }

  return newName;
}

int ALIUtils::checkMatrixEquations ( double  angleX, double  angleY, double  angleZ, CLHEP::HepRotation *  rot  )

static

Definition at line 762 of file ALIUtils.cc.

References beam_dqm_sourceclient-live_cfg::cerr, funct::cos(), eq2ang(), funct::sin(), fftjetcommon_cfi::sx, and fftjetcommon_cfi::sy.

Referenced by getRotationAnglesFromMatrix(), and setMaximumDeviationDerivative().

                                                                                                     {
  //-  std::cout << " cme " << angleX << " " << angleY << " " << angleZ << std::endl;
  if (rot == nullptr) {
    rot = new CLHEP::HepRotation();
    rot->rotateX(angleX);
    rot->rotateY(angleY);
    rot->rotateZ(angleZ);
  }
  double sx = sin(angleX);
  double cx = cos(angleX);
  double sy = sin(angleY);
  double cy = cos(angleY);
  double sz = sin(angleZ);
  double cz = cos(angleZ);

  double rotxx = cy * cz;
  double rotxy = sx * sy * cz - cx * sz;
  double rotxz = cx * sy * cz + sx * sz;
  double rotyx = cy * sz;
  double rotyy = sx * sy * sz + cx * cz;
  double rotyz = cx * sy * sz - sx * cz;
  double rotzx = -sy;
  double rotzy = sx * cy;
  double rotzz = cx * cy;

  int matrixElemBad = 0;
  if (!eq2ang(rot->xx(), rotxx)) {
    std::cerr << " EQUATION for xx() IS BAD " << rot->xx() << " <> " << rotxx << std::endl;
    matrixElemBad++;
  }
  if (!eq2ang(rot->xy(), rotxy)) {
    std::cerr << " EQUATION for xy() IS BAD " << rot->xy() << " <> " << rotxy << std::endl;
    matrixElemBad++;
  }
  if (!eq2ang(rot->xz(), rotxz)) {
    std::cerr << " EQUATION for xz() IS BAD " << rot->xz() << " <> " << rotxz << std::endl;
    matrixElemBad++;
  }
  if (!eq2ang(rot->yx(), rotyx)) {
    std::cerr << " EQUATION for yx() IS BAD " << rot->yx() << " <> " << rotyx << std::endl;
    matrixElemBad++;
  }
  if (!eq2ang(rot->yy(), rotyy)) {
    std::cerr << " EQUATION for yy() IS BAD " << rot->yy() << " <> " << rotyy << std::endl;
    matrixElemBad++;
  }
  if (!eq2ang(rot->yz(), rotyz)) {
    std::cerr << " EQUATION for yz() IS BAD " << rot->yz() << " <> " << rotyz << std::endl;
    matrixElemBad++;
  }
  if (!eq2ang(rot->zx(), rotzx)) {
    std::cerr << " EQUATION for zx() IS BAD " << rot->zx() << " <> " << rotzx << std::endl;
    matrixElemBad++;
  }
  if (!eq2ang(rot->zy(), rotzy)) {
    std::cerr << " EQUATION for zy() IS BAD " << rot->zy() << " <> " << rotzy << std::endl;
    matrixElemBad++;
  }
  if (!eq2ang(rot->zz(), rotzz)) {
    std::cerr << " EQUATION for zz() IS BAD " << rot->zz() << " <> " << rotzz << std::endl;
    matrixElemBad++;
  }

  //-  std::cout << " cme: matrixElemBad " << matrixElemBad << std::endl;
  return matrixElemBad;
}

double ALIUtils::diff2pi ( double  ang1, double  ang2  )

static

Definition at line 717 of file ALIUtils.cc.

References change_name::diff, and createfilelist::int.

Referenced by eq2ang(), getRotationAnglesFromMatrix(), and setMaximumDeviationDerivative().

                                                 {
  double pii = acos(0.) * 2;
  double diff = fabs(ang1 - ang2);
  diff = diff - int(diff / 2. / pii) * 2 * pii;
  return diff;
}

void ALIUtils::dump3v ( const CLHEP::Hep3Vector &  vec, const std::string &  msg  )

static

Definition at line 58 of file ALIUtils.cc.

References gather_cfg::cout.

Referenced by ALILine::ALILine(), ALIPlane::ALIPlane(), OptOOpticalSquare::calculateFaces(), OptOCOPS::convertPointToLocalCoordinates(), OptOSensor2D::convertPointToLocalCoordinates(), OptOPinhole::defaultBehaviour(), OptOPlateSplitter::detailedDeviatesLightRay(), OptOMirror::detailedDeviatesLightRay(), OptOCubeSplitter::detailedDeviatesLightRay(), OptOModifiedRhomboidPrism::detailedDeviatesLightRay(), OpticalObject::displaceCentreGlobOriginal(), OpticalObject::displaceCentreGlobOriginalOriginal(), OpticalObject::displaceRmGlobAroundGlobal(), OpticalObject::displaceRmGlobAroundLocal(), OpticalObject::displaceRmGlobOriginal(), OpticalObject::displaceRmGlobOriginalOriginal(), LightRay::dumpData(), OpticalObjectMgr::dumpOptOs(), OpticalObject::getDispVec(), OptOCubeSplitter::getMiddlePlate(), OpticalObject::getPlate(), OptOModifiedRhomboidPrism::getRotatedPlate(), OptOCubeSplitter::getUpperPlate(), ALILine::intersect(), LightRay::intersect(), OptOTiltmeter::makeMeasurement(), OptODistancemeter::makeMeasurement(), OptODistancemeter3dim::makeMeasurement(), OptOSensor2D::makeMeasurement(), OptOCOPS::makeMeasurement(), LightRay::reflect(), LightRay::refract(), OpticalObject::SetCentreGlobFromCentreLocal(), OpticalObject::SetCentreLocalFromEntryValues(), LightRay::shiftAndDeviateWhileTraversing(), OptOUserDefined::userDefinedBehaviour(), and ~ALIUtils().

                                                                      {
  //  double phicyl = atan( vec.y()/vec.x() );
  std::cout << msg << std::setprecision(8) << vec;
  std::cout << std::endl;
  //  std::cout << " " << vec.theta()/deg << " " << vec.phi()/deg << " " << vec.perp() << " " << phicyl/deg << std::endl;
  //  setw(10);
  //  std::cout << msg << " x=" << std::setprecision(8) << vec.x() << " y=" << setprecision(8) <<vec.y() << " z=" << std::setprecision(8) << vec.z() << std::endl;
  // std::setprecision(8);
}

void ALIUtils::dumpDimensions ( std::ofstream &  fout )

static

Definition at line 308 of file ALIUtils.cc.

References _OutputAngleSigmaDimensionFactor, _OutputAngleValueDimensionFactor, _OutputLengthSigmaDimensionFactor, _OutputLengthValueDimensionFactor, beam_dqm_sourceclient-live_cfg::cerr, beamvalidation::exit(), and M_PI.

Referenced by Fit::fitNextEvent(), FittedEntriesManager::MakeHistos(), and set_time_now().

                                               {
  fout << "DIMENSIONS: lengths = ";
  ALIstring internalDim = "m";
  if (_OutputLengthValueDimensionFactor == 1.) {
    fout << "m";
  } else if (_OutputLengthValueDimensionFactor == 1.E-3) {
    fout << "mm";
  } else if (_OutputLengthValueDimensionFactor == 1.E-6) {
    fout << "mum";
  } else if (_OutputLengthValueDimensionFactor == 1.E-2) {
    fout << "cm";
  } else {
    std::cerr << " !! unknown OutputLengthValueDimensionFactor " << _OutputLengthValueDimensionFactor << std::endl;
    exit(1);
  }

  fout << " +- ";
  if (_OutputLengthSigmaDimensionFactor == 1.) {
    fout << "m";
  } else if (_OutputLengthSigmaDimensionFactor == 1.E-3) {
    fout << "mm";
  } else if (_OutputLengthSigmaDimensionFactor == 1.E-6) {
    fout << "mum";
  } else if (_OutputLengthSigmaDimensionFactor == 1.E-2) {
    fout << "cm";
  } else {
    std::cerr << " !! unknown OutputLengthSigmaDimensionFactor " << _OutputLengthSigmaDimensionFactor << std::endl;
    exit(1);
  }

  fout << "  angles = ";
  if (_OutputAngleValueDimensionFactor == 1.) {
    fout << "rad";
  } else if (_OutputAngleValueDimensionFactor == 1.E-3) {
    fout << "mrad";
  } else if (_OutputAngleValueDimensionFactor == 1.E-6) {
    fout << "murad";
  } else if (_OutputAngleValueDimensionFactor == M_PI / 180.) {
    fout << "deg";
  } else if (_OutputAngleValueDimensionFactor == M_PI / 200.) {
    fout << "grad";
  } else {
    std::cerr << " !! unknown OutputAngleValueDimensionFactor " << _OutputAngleValueDimensionFactor << std::endl;
    exit(1);
  }

  fout << " +- ";
  if (_OutputAngleSigmaDimensionFactor == 1.) {
    fout << "rad";
  } else if (_OutputAngleSigmaDimensionFactor == 1.E-3) {
    fout << "mrad";
  } else if (_OutputAngleSigmaDimensionFactor == 1.E-6) {
    fout << "murad";
  } else if (_OutputAngleSigmaDimensionFactor == M_PI / 180.) {
    fout << "deg";
  } else if (_OutputAngleSigmaDimensionFactor == M_PI / 200.) {
    fout << "grad";
  } else {
    std::cerr << " !! unknown OutputAngleSigmaDimensionFactor " << _OutputAngleSigmaDimensionFactor << std::endl;
    exit(1);
  }
  fout << std::endl;
}

void ALIUtils::dumprm ( const CLHEP::HepRotation &  rm, const std::string &  msg, std::ostream &  out = std::cout  )

static

Definition at line 71 of file ALIUtils.cc.

Referenced by OpticalObject::buildRmFromEntryValuesOriginalOriginal(), OptOCOPS::convertPointToLocalCoordinates(), OptOSensor2D::convertPointToLocalCoordinates(), OpticalObject::displaceRmGlobAroundGlobal(), OpticalObject::displaceRmGlobAroundLocal(), OpticalObject::displaceRmGlobOriginal(), OpticalObject::displaceRmGlobOriginalOriginal(), OpticalObject::getDispVec(), OptOCubeSplitter::getMiddlePlate(), OpticalObject::getPlate(), OptOModifiedRhomboidPrism::getRotatedPlate(), getRotationAnglesFromMatrix(), OpticalObject::getRotationAnglesInOptOFrame(), OptOCubeSplitter::getUpperPlate(), OptOTiltmeter::makeMeasurement(), Fit::printCentreInOptOFrame(), OpticalObject::propagateGlobalRMOriginalOriginalChangeToChildren(), OpticalObject::SetCentreGlobFromCentreLocal(), OpticalObject::setGlobalRMOriginalOriginal(), OpticalObject::SetRMGlobFromRMLocal(), OpticalObject::SetRMLocalFromEntryValues(), LightRay::startLightRay(), EntryAngleAffAngles::valueDisplaced(), and ~ALIUtils().

                                                                                       {
  out << msg << " xx=" << rm.xx() << " xy=" << rm.xy() << " xz=" << rm.xz() << std::endl;
  out << msg << " yx=" << rm.yx() << " yy=" << rm.yy() << " yz=" << rm.yz() << std::endl;
  out << msg << " zx=" << rm.zx() << " zy=" << rm.zy() << " zz=" << rm.zz() << std::endl;
}

void ALIUtils::dumpVS ( const std::vector< ALIstring > &  wl, const std::string &  msg, std::ostream &  outs = std::cout  )

static

dumps a vector of strings with a message to outs

Definition at line 465 of file ALIUtils.cc.

References cuy::ii.

Referenced by OpticalObject::buildWordList(), Measurement::constructFromOA(), OpticalObject::copyData(), EntryData::fill(), Entry::fill(), DeviationSensor2D::fillData(), OpticalObject::fillExtraEntry(), FittedEntriesReader::FittedEntriesReader(), OpticalObject::getCoordinateFromOptAlignParam(), DeviationsFromFileSensor2D::readFile(), ErrorCorrelationMgr::readFromReportFile(), Model::readMeasurementsFromFile(), CocoaDaqReaderText::ReadNextEvent(), Model::readSystemDescription(), set_time_now(), MeasurementDistancemeter::setConversionFactor(), MeasurementDistancemeter3dim::setConversionFactor(), MeasurementSensor2D::setConversionFactor(), MeasurementCOPS::setConversionFactor(), MeasurementTiltmeter::setConversionFactor(), and Measurement::setCurrentDate().

                                                                                            {
  outs << msg << std::endl;
  ALIuint siz = wl.size();
  for (ALIuint ii = 0; ii < siz; ii++) {
    outs << wl[ii] << " ";
    /*  ostream_iterator<ALIstring> os(outs," ");
    copy(wl.begin(), wl.end(), os);*/
  }
  outs << std::endl;
}

bool ALIUtils::eq2ang ( double  ang1, double  ang2  )

static

Definition at line 725 of file ALIUtils.cc.

References change_name::diff, and diff2pi().

Referenced by checkMatrixEquations(), getRotationAnglesFromMatrix(), and setMaximumDeviationDerivative().

                                              {
  bool beq = true;

  double pii = acos(0.) * 2;
  double diff = diff2pi(ang1, ang2);
  if (diff > 0.00001) {
    if (fabs(diff - 2 * pii) > 0.00001) {
      //-      std::cout << " diff " << diff << " " << ang1 << " " << ang2 << std::endl;
      beq = false;
    }
  } else {
    beq = true;
  }

  return beq;
}

bool ALIUtils::getBool ( const ALIstring &  str )

static

Convert a bool to an integer, checking that it is really a bool.

Definition at line 420 of file ALIUtils.cc.

References beam_dqm_sourceclient-live_cfg::cerr, beamvalidation::exit(), and heppy_batch::val.

Referenced by set_time_now().

                                           {
  bool val;

  //t str = upper( str );
  //----------- first check that it is a not number
  if (str == "ON" || str == "TRUE") {
    val = true;
  } else if (str == "OFF" || str == "FALSE") {
    val = false;
  } else {
    std::cerr << "!!!! EXITING: trying to get the float from a string that is not 'ON'/'OFF'/'TRUE'/'FALSE' " << str
              << std::endl;
    exit(1);
  }

  return val;
}

ALIdouble ALIUtils::getDimensionValue ( const ALIstring &  dim, const ALIstring &  dimType  )

static

Definition at line 477 of file ALIUtils.cc.

References beam_dqm_sourceclient-live_cfg::cerr, M_PI, and relativeConstraints::value.

Referenced by EntryMgr::readEntryFromReportOut(), and val0().

                                                                                    {
  ALIdouble value;
  if (dimType == "Length") {
    if (dim == "mm") {
      value = 1.E-3;
    } else if (dim == "cm") {
      value = 1.E-2;
    } else if (dim == "m") {
      value = 1.;
    } else if (dim == "mum") {
      value = 1.E-6;
    } else if (dim == "dm") {
      value = 1.E-1;
    } else if (dim == "nm") {
      value = 1.E-9;
    } else {
      std::cerr << "!!!!FATAL ERROR:  ALIUtils::GetDimensionValue. " << dim
                << " is a dimension not supported for dimension type " << dimType << std::endl;
      abort();
    }
  } else if (dimType == "Angle") {
    if (dim == "rad") {
      value = 1.;
    } else if (dim == "mrad") {
      value = 1.E-3;
    } else if (dim == "murad") {
      value = 1.E-6;
    } else if (dim == "deg") {
      value = M_PI / 180.;
    } else if (dim == "grad") {
      value = M_PI / 200.;
    } else {
      std::cerr << "!!!!FATAL ERROR:  ALIUtils::GetDimensionValue. " << dim
                << " is a dimension not supported for dimension type " << dimType << std::endl;
      abort();
    }
  } else {
    std::cerr << "!!!!FATAL ERROR:  ALIUtils::GetDimensionValue. " << dimType << " is a dimension type not supported "
              << std::endl;
    abort();
  }

  return value;
}

static ALIbool ALIUtils::getFirstTime ( )

inlinestatic

Definition at line 87 of file ALIUtils.h.

References firstTime.

Referenced by Measurement::calculateOriginalSimulatedValue(), MeasurementSensor2D::calculateSimulatedValue(), MeasurementDistancemeter::calculateSimulatedValue(), MeasurementDistancemeter3dim::calculateSimulatedValue(), MeasurementCOPS::calculateSimulatedValue(), MeasurementDiffEntry::calculateSimulatedValue(), Fit::PropagateErrors(), OptOUserDefined::userDefinedBehaviour(), and Fit::WriteVisualisationFiles().

{ return firstTime; }

double ALIUtils::getFloat ( const ALIstring &  str )

static

Convert a string to an float, checking that it is really a number.

Definition at line 373 of file ALIUtils.cc.

References beam_dqm_sourceclient-live_cfg::cerr, beamvalidation::exit(), and IsNumber().

Referenced by EntryData::fill(), DeviationSensor2D::fillData(), FittedEntriesSet::FillEntriesFromFile(), OptOSensor2D::fillExtraEntry(), FittedEntriesReader::readFittedEntriesFromFile(), ErrorCorrelationMgr::readFromReportFile(), Model::readSystemDescription(), and set_time_now().

                                              {
  //----------- first check that it is a number
  if (!IsNumber(str)) {
    std::cerr << "!!!! EXITING: trying to get the float from a string that is not a number " << str << std::endl;
    exit(1);
  }

  return atof(str.c_str());
}

int ALIUtils::getInt ( const ALIstring &  str )

static

Convert a string to an integer, checking that it is really an integer.

Definition at line 384 of file ALIUtils.cc.

References beam_dqm_sourceclient-live_cfg::cerr, beamvalidation::exit(), cuy::ii, createfilelist::int, and IsNumber().

Referenced by OpticalObject::createComponentOptOs(), ErrorCorrelationMgr::readFromReportFile(), Model::readSystemDescription(), and set_time_now().

                                         {
  //----------- first check that it is an integer
  if (!IsNumber(str)) {
    //----- Check that it is a number
    std::cerr << "!!!! EXITING: trying to get the integer from a string that is not a number " << str << std::endl;
    exit(1);
  } else {
    //----- Check that it is not a float, no decimal or E-n
    bool isFloat = false;
    int ch = str.find('.');
    ALIuint ii = 0;
    if (ch != -1) {
      for (ii = ch + 1; ii < str.size(); ii++) {
        if (str[ii] != '0')
          isFloat = true;
      }
    }

    ch = str.find('E');
    if (ch != -1)
      ch = str.find('e');
    if (ch != -1) {
      if (str[ch + 1] == '-')
        isFloat = true;
    }

    if (isFloat) {
      std::cerr << "!!!! EXITING: trying to get the integer from a string that is a float: " << str << std::endl;
      std::cerr << ii << " ii " << ch << std::endl;
      exit(1);
    }
  }
  return int(atof(str.c_str()));
}

static ALIdouble ALIUtils::getMaximumDeviationDerivative ( )

inlinestatic

Definition at line 89 of file ALIUtils.h.

References maximum_deviation_derivative.

Referenced by Measurement::DerivativeRespectEntry(), and Fit::getInstance().

{ return maximum_deviation_derivative; }

std::vector< double > ALIUtils::getRotationAnglesFromMatrix ( CLHEP::HepRotation &  rmLocal, double  origAngleX, double  origAngleY, double  origAngleZ  )

static

Definition at line 549 of file ALIUtils.cc.

References addPii(), approxTo0(), beam_dqm_sourceclient-live_cfg::cerr, checkMatrixEquations(), funct::cos(), gather_cfg::cout, debug, diff2pi(), dumprm(), eq2ang(), cuy::ii, makeMuonMisalignmentScenario::rot, and funct::sin().

Referenced by OpticalObject::getRotationAnglesFromMatrix(), and setMaximumDeviationDerivative().

                                                                             {
  double pii = acos(0.) * 2;
  std::vector<double> newang(3);
  double angleX = origAngleX;
  double angleY = origAngleY;
  double angleZ = origAngleZ;

  if (ALIUtils::debug >= 4) {
    std::cout << " angles as value entries: X= " << angleX << " Y= " << angleY << " Z " << angleZ << std::endl;
  }

  //-  std::cout << name () << " vdbf " << angleX << " " << angleY << " " << angleZ << std::endl;
  double rotzx = approxTo0(rmLocal.zx());
  double rotzy = approxTo0(rmLocal.zy());
  double rotzz = approxTo0(rmLocal.zz());
  double rotyx = approxTo0(rmLocal.yx());
  double rotxx = approxTo0(rmLocal.xx());
  if (rotzy == 0. && rotzz == 0.) {
    //check that entry is z angle
    newang[0] = angleX;
    //beware of aa <==> pii - aa
    if (eq2ang(rmLocal.zx(), -1.)) {
      double aa = asin(rmLocal.xy());
      if (diff2pi(angleZ, -aa + newang[0]) < diff2pi(angleZ, -pii + aa + newang[0])) {
        newang[2] = -aa + newang[0];
        if (ALIUtils::debug >= 5)
          std::cout << " newang[0] = -aa + newang[0] " << std::endl;
      } else {
        newang[2] = -pii + aa + newang[0];
        if (ALIUtils::debug >= 5)
          std::cout << " newang[0] = -pii + aa + newang[0] " << newang[0] << " " << aa << " " << newang[2] << std::endl;
      }
    } else {
      double aa = asin(-rmLocal.xy());
      if (diff2pi(angleZ, aa - newang[0]) < diff2pi(angleZ, pii - aa - newang[0])) {
        newang[2] = aa - newang[0];
        if (ALIUtils::debug >= 5)
          std::cout << " newang[0] = aa - newang[2] " << std::endl;
      } else {
        newang[2] = pii - aa - newang[0];
        if (ALIUtils::debug >= 5)
          std::cout << " newang[0] = pii - aa - newang[2] " << newang[0] << " " << aa << " " << newang[2] << std::endl;
      }
    }
  } else {
    newang[0] = atan(rotzy / rotzz);
    newang[2] = atan(rotyx / rotxx);
  }
  if (rotzx < -1.) {
    //-    std::cerr << " rotzx too small " << rotzx << " = " << rmLocal.zx() << " " << rotzx-rmLocal.zx() << std::endl;
    rotzx = -1.;
  } else if (rotzx > 1.) {
    //-    std::cerr << " rotzx too big " << rotzx << " = " << rmLocal.zx() << " " << rotzx-rmLocal.zx() << std::endl;
    rotzx = 1.;
  }
  newang[1] = -asin(rotzx);
  if (ALIUtils::debug >= 5)
    std::cout << "First calculation of angles: " << std::endl
              << " newang[0] " << newang[0] << " " << rotzy << " " << rotzz << std::endl
              << " newang[1] " << newang[1] << " " << rotzx << std::endl
              << " newang[2] " << newang[2] << " " << rotyx << " " << rotxx << std::endl;

  //    newang[2] = acos( rmLocal.xx() / cos( newang[1] ) );
  //----- CHECK if the angles are OK (there are several symmetries)
  //--- Check if the predictions with the angles obtained match the values of the rotation matrix (they may differ for exampole by a sign or more in complicated formulas)
  double rotnewxx = cos(newang[1]) * cos(newang[2]);
  double rotnewzz = cos(newang[0]) * cos(newang[1]);
  double rotnewxy = sin(newang[0]) * sin(newang[1]) * cos(newang[2]) - cos(newang[0]) * sin(newang[2]);
  double rotnewxz = cos(newang[0]) * sin(newang[1]) * cos(newang[2]) + sin(newang[0]) * sin(newang[2]);
  double rotnewyy = sin(newang[0]) * sin(newang[1]) * sin(newang[2]) + cos(newang[0]) * cos(newang[2]);
  double rotnewyz = cos(newang[0]) * sin(newang[1]) * sin(newang[2]) - sin(newang[0]) * cos(newang[2]);

  bool eqxx = eq2ang(rotnewxx, rmLocal.xx());
  bool eqzz = eq2ang(rotnewzz, rmLocal.zz());
  bool eqxy = eq2ang(rotnewxy, rmLocal.xy());
  bool eqxz = eq2ang(rotnewxz, rmLocal.xz());
  bool eqyy = eq2ang(rotnewyy, rmLocal.yy());
  bool eqyz = eq2ang(rotnewyz, rmLocal.yz());

  //--- Check if one of the tree angles should be changed
  if (ALIUtils::debug >= 5) {
    std::cout << " pred rm.xx " << rotnewxx << " =? " << rmLocal.xx() << " pred rm.zz " << rotnewzz << " =? "
              << rmLocal.zz() << std::endl;
    std::cout << " eqxx " << eqxx << " eqzz " << eqzz << std::endl;
    //-    std::cout << " rotnewxx " << rotnewxx << " = " << rmLocal.xx() << " " << fabs( rotnewxx - rmLocal.xx() ) << " " <<(fabs( rotnewxx - rmLocal.xx() ) < 0.0001) << std::endl;
  }

  if (eqxx & !eqzz) {
    newang[0] = pii + newang[0];
    if (ALIUtils::debug >= 5)
      std::cout << " change newang[0] " << newang[0] << std::endl;
  } else if (!eqxx & !eqzz) {
    newang[1] = pii - newang[1];
    if (ALIUtils::debug >= 5)
      std::cout << " change newang[1] " << newang[1] << std::endl;
  } else if (!eqxx & eqzz) {
    newang[2] = pii + newang[2];
    if (ALIUtils::debug >= 5)
      std::cout << " change newang[2] " << newang[2] << std::endl;
  }

  //--- Check if the 3 angles should be changed (previous check is invariant to the 3 changing)
  if (ALIUtils::debug >= 5) {
    std::cout << " pred rm.xy " << rotnewxy << " =? " << rmLocal.xy() << " pred rm.xz " << rotnewxz << " =? "
              << rmLocal.xz() << " pred rm.yy " << rotnewyy << " =? " << rmLocal.yy() << " pred rm.yz " << rotnewyz
              << " =? " << rmLocal.yz() << std::endl;
    std::cout << " eqxy " << eqxy << " eqxz " << eqxz << " eqyy " << eqyy << " eqyz " << eqyz << std::endl;
  }

  if (!eqxy || !eqxz || !eqyy || !eqyz) {
    // check also cases where one of the above 'eq' is OK because it is = 0
    if (ALIUtils::debug >= 5)
      std::cout << " change the 3 newang " << std::endl;
    newang[0] = addPii(newang[0]);
    newang[1] = pii - newang[1];
    newang[2] = addPii(newang[2]);
    double rotnewxy = -sin(newang[0]) * sin(newang[1]) * cos(newang[2]) - cos(newang[0]) * sin(newang[2]);
    double rotnewxz = -cos(newang[0]) * sin(newang[1]) * cos(newang[2]) - sin(newang[0]) * sin(newang[2]);
    if (ALIUtils::debug >= 5)
      std::cout << " rotnewxy " << rotnewxy << " = " << rmLocal.xy() << " rotnewxz " << rotnewxz << " = "
                << rmLocal.xz() << std::endl;
  }
  if (diff2pi(angleX, newang[0]) + diff2pi(angleY, newang[1]) + diff2pi(angleZ, newang[2]) >
      diff2pi(angleX, pii + newang[0]) + diff2pi(angleY, pii - newang[1]) + diff2pi(angleZ, pii + newang[2])) {
    // check also cases where one of the above 'eq' is OK because it is = 0
    if (ALIUtils::debug >= 5)
      std::cout << " change the 3 newang " << std::endl;
    newang[0] = addPii(newang[0]);
    newang[1] = pii - newang[1];
    newang[2] = addPii(newang[2]);
    double rotnewxy = -sin(newang[0]) * sin(newang[1]) * cos(newang[2]) - cos(newang[0]) * sin(newang[2]);
    double rotnewxz = -cos(newang[0]) * sin(newang[1]) * cos(newang[2]) - sin(newang[0]) * sin(newang[2]);
    if (ALIUtils::debug >= 5)
      std::cout << " rotnewxy " << rotnewxy << " = " << rmLocal.xy() << " rotnewxz " << rotnewxz << " = "
                << rmLocal.xz() << std::endl;
  }

  for (int ii = 0; ii < 3; ii++) {
    newang[ii] = approxTo0(newang[ii]);
  }
  //  double rotnewyx = cos( newang[1] ) * sin( newang[2] );

  if (checkMatrixEquations(newang[0], newang[1], newang[2], &rmLocal) != 0) {
    std::cerr << " wrong rotation matrix " << newang[0] << " " << newang[1] << " " << newang[2] << std::endl;
    ALIUtils::dumprm(rmLocal, " matrix is ");
  }
  if (ALIUtils::debug >= 5) {
    std::cout << "Final angles:  newang[0] " << newang[0] << " newang[1] " << newang[1] << " newang[2] " << newang[2]
              << std::endl;
    CLHEP::HepRotation rot;
    rot.rotateX(newang[0]);
    ALIUtils::dumprm(rot, " new rot after X ");
    rot.rotateY(newang[1]);
    ALIUtils::dumprm(rot, " new rot after Y ");
    rot.rotateZ(newang[2]);
    ALIUtils::dumprm(rot, " new rot ");
    ALIUtils::dumprm(rmLocal, " rmLocal ");
    //-    ALIUtils::dumprm( theRmGlobOriginal, " theRmGlobOriginal " );
  }

  //-  std::cout << " before return newang[0] " << newang[0] << " newang[1] " << newang[1] << " newang[2] " << newang[2] << std::endl;
  return newang;
}

int ALIUtils::IsNumber ( const ALIstring &  str )

static

Definition at line 33 of file ALIUtils.cc.

References cuy::ii.

Referenced by Measurement::fillData(), Entry::fillFromInputFileSigma(), Entry::fillFromInputFileValue(), getFloat(), getInt(), ParameterMgr::getVal(), Model::readSystemDescription(), MeasurementDistancemeter3dim::setConversionFactor(), MeasurementDistancemeter::setConversionFactor(), MeasurementSensor2D::setConversionFactor(), MeasurementCOPS::setConversionFactor(), MeasurementTiltmeter::setConversionFactor(), and ~ALIUtils().

                                           {
  int isnum = 1;
  int numE = 0;
  for (ALIuint ii = 0; ii < str.length(); ii++) {
    if (!isdigit(str[ii]) && str[ii] != '.' && str[ii] != '-' && str[ii] != '+') {
      //--- check for E(xponential)
      if (str[ii] == 'E' || str[ii] == 'e') {
        if (numE != 0 || ii == str.length() - 1) {
          isnum = 0;
          break;
        }
        numE++;
      } else {
        isnum = 0;
        break;
      }
    }
  }

  return isnum;
}

static ALIdouble ALIUtils::LengthSigmaDimensionFactor ( )

inlinestatic

Definition at line 63 of file ALIUtils.h.

References _LengthSigmaDimensionFactor.

Referenced by EntryLength::SigmaDimensionFactor(), and Measurement::sigmaDimensionFactor().

{ return _LengthSigmaDimensionFactor; }

static ALIdouble ALIUtils::LengthValueDimensionFactor ( )

inlinestatic

Definition at line 62 of file ALIUtils.h.

References _LengthValueDimensionFactor.

Referenced by OptOSensor2D::fastTraversesLightRay(), EntryMgr::readEntryFromReportOut(), OpticalObject::transformCylindrical2Cartesian(), EntryLength::ValueDimensionFactor(), and Measurement::valueDimensionFactor().

{ return _LengthValueDimensionFactor; }

static ALIdouble ALIUtils::OutputAngleSigmaDimensionFactor ( )

inlinestatic

Definition at line 69 of file ALIUtils.h.

References _OutputAngleSigmaDimensionFactor.

Referenced by EntryAngle::OutputSigmaDimensionFactor().

{ return _OutputAngleSigmaDimensionFactor; }

static ALIdouble ALIUtils::OutputAngleValueDimensionFactor ( )

inlinestatic

Definition at line 68 of file ALIUtils.h.

References _OutputAngleValueDimensionFactor.

Referenced by EntryAngle::OutputValueDimensionFactor().

{ return _OutputAngleValueDimensionFactor; }

static ALIdouble ALIUtils::OutputLengthSigmaDimensionFactor ( )

inlinestatic

Definition at line 67 of file ALIUtils.h.

References _OutputLengthSigmaDimensionFactor.

Referenced by EntryLength::OutputSigmaDimensionFactor().

{ return _OutputLengthSigmaDimensionFactor; }

static ALIdouble ALIUtils::OutputLengthValueDimensionFactor ( )

inlinestatic

Definition at line 66 of file ALIUtils.h.

References _OutputLengthValueDimensionFactor.

Referenced by EntryLength::OutputValueDimensionFactor().

{ return _OutputLengthValueDimensionFactor; }

static void ALIUtils::set_time_now ( time_t  now )

inlinestatic

Definition at line 40 of file ALIUtils.h.

References _time_now, CalculateAngleDimensionFactorFromInt(), CalculateAngleDimensionFactorFromString(), CalculateLengthDimensionFactorFromInt(), CalculateLengthDimensionFactorFromString(), gather_cfg::cout, dumpDimensions(), dumpVS(), groupFilesInBlocks::fout, getBool(), getFloat(), getInt(), fileCollector::now, SetAngleDimensionFactors(), SetLengthDimensionFactors(), SetOutputAngleDimensionFactors(), SetOutputLengthDimensionFactors(), AlCaHLTBitMon_QueryRunRegistry::string, and LaserClient_cfi::wl.

Referenced by Fit::getFitQuality(), Fit::multiplyMatrices(), and Fit::PropagateErrors().

{ _time_now = now; }

void ALIUtils::SetAngleDimensionFactors ( )

static

Definition at line 104 of file ALIUtils.cc.

References _AngleSigmaDimensionFactor, _AngleValueDimensionFactor, CalculateAngleDimensionFactorFromInt(), gather_cfg::cout, debug, GlobalOptionMgr::getInstance(), and GlobalOptionMgr::GlobalOptions().

Referenced by Model::readSystemDescription(), and set_time_now().

                                        {
  //--------------------- if it doesn exist, GlobalOptions is 0
  //---------- Calculate factors to convert to radians
  GlobalOptionMgr* gomgr = GlobalOptionMgr::getInstance();
  ALIint ad = ALIint(gomgr->GlobalOptions()[ALIstring("angle_value_dimension")]);
  _AngleValueDimensionFactor = CalculateAngleDimensionFactorFromInt(ad);

  ad = ALIint(gomgr->GlobalOptions()[ALIstring("angle_error_dimension")]);
  _AngleSigmaDimensionFactor = CalculateAngleDimensionFactorFromInt(ad);

  //---------- Change factor to convert to error dimensions
  //  _AngleValueDimensionFactor /= _AngleSigmaDimensionFactor;
  //_AngleSigmaDimensionFactor = 1;

  if (ALIUtils::debug >= 6)
    std::cout << _AngleValueDimensionFactor << "Set Angle DimensionFactors" << _AngleSigmaDimensionFactor << std::endl;
}

static void ALIUtils::setDebugVerbosity ( ALIint  val )

inlinestatic

Definition at line 38 of file ALIUtils.h.

References heppy_batch::val.

Referenced by GlobalOptionMgr::setGlobalOption().

{ debug = val; }

static void ALIUtils::setFirstTime ( ALIbool  val )

inlinestatic

Definition at line 88 of file ALIUtils.h.

References firstTime, and heppy_batch::val.

Referenced by Fit::PropagateErrors(), and Fit::startFit().

{ firstTime = val; }

void ALIUtils::SetLengthDimensionFactors ( )

static

Definition at line 81 of file ALIUtils.cc.

References _LengthSigmaDimensionFactor, _LengthValueDimensionFactor, CalculateLengthDimensionFactorFromInt(), gather_cfg::cout, debug, GlobalOptionMgr::getGlobalOption(), GlobalOptionMgr::getInstance(), and GlobalOptionMgr::GlobalOptions().

Referenced by Model::readSystemDescription(), and set_time_now().

                                         {
  //---------------------------------------- if it doesn exist, GlobalOptions is 0
  //---------- Calculate factors to convert to meters
  GlobalOptionMgr* gomgr = GlobalOptionMgr::getInstance();
  ALIint ad = ALIint(gomgr->getGlobalOption("length_value_dimension"));

  _LengthValueDimensionFactor = CalculateLengthDimensionFactorFromInt(ad);
  ad = ALIint(gomgr->GlobalOptions()[ALIstring("length_error_dimension")]);
  _LengthSigmaDimensionFactor = CalculateLengthDimensionFactorFromInt(ad);

  //---------- Change factor to convert to error dimensions
  //  _LengthValueDimensionFactor /= _LengthSigmaDimensionFactor;
  //_LengthSigmaDimensionFactor = 1;

  if (ALIUtils::debug >= 6)
    std::cout << _LengthValueDimensionFactor << " Set Length DimensionFactors " << _LengthSigmaDimensionFactor
              << std::endl;
}

static void ALIUtils::setMaximumDeviationDerivative ( ALIdouble  val )

inlinestatic

Definition at line 90 of file ALIUtils.h.

References addPii(), approxTo0(), checkMatrixEquations(), diff2pi(), eq2ang(), getRotationAnglesFromMatrix(), maximum_deviation_derivative, makeMuonMisalignmentScenario::rot, and heppy_batch::val.

Referenced by Fit::getInstance().

{ maximum_deviation_derivative = val; }

void ALIUtils::SetOutputAngleDimensionFactors ( )

static

Definition at line 153 of file ALIUtils.cc.

References _OutputAngleSigmaDimensionFactor, _OutputAngleValueDimensionFactor, CalculateAngleDimensionFactorFromInt(), gather_cfg::cout, debug, GlobalOptionMgr::getInstance(), and GlobalOptionMgr::GlobalOptions().

Referenced by Model::readSystemDescription(), and set_time_now().

                                              {
  //--------------------- if it doesn exist, GlobalOptions is 0
  //---------- Calculate factors to convert to radians
  GlobalOptionMgr* gomgr = GlobalOptionMgr::getInstance();
  ALIint ad = ALIint(gomgr->GlobalOptions()[ALIstring("output_angle_value_dimension")]);
  if (ad == 0)
    ad = ALIint(gomgr->GlobalOptions()[ALIstring("angle_value_dimension")]);
  _OutputAngleValueDimensionFactor = CalculateAngleDimensionFactorFromInt(ad);

  ad = ALIint(gomgr->GlobalOptions()[ALIstring("output_angle_error_dimension")]);
  if (ad == 0)
    ad = ALIint(gomgr->GlobalOptions()[ALIstring("angle_error_dimension")]);
  _OutputAngleSigmaDimensionFactor = CalculateAngleDimensionFactorFromInt(ad);

  //---------- Change factor to convert to error dimensions
  //  _AngleValueDimensionFactor /= _AngleSigmaDimensionFactor;
  //_AngleSigmaDimensionFactor = 1;

  if (ALIUtils::debug >= 9)
    std::cout << _OutputAngleValueDimensionFactor << "Output Angle Dimension Factors"
              << _OutputAngleSigmaDimensionFactor << std::endl;
}

void ALIUtils::SetOutputLengthDimensionFactors ( )

static

Definition at line 126 of file ALIUtils.cc.

References _OutputLengthSigmaDimensionFactor, _OutputLengthValueDimensionFactor, CalculateLengthDimensionFactorFromInt(), gather_cfg::cout, debug, GlobalOptionMgr::getInstance(), and GlobalOptionMgr::GlobalOptions().

Referenced by Model::readSystemDescription(), and set_time_now().

                                               {
  //---------------------------------------- if it doesn exist, GlobalOptions is 0
  //---------- Calculate factors to convert to meters
  GlobalOptionMgr* gomgr = GlobalOptionMgr::getInstance();
  ALIint ad = ALIint(gomgr->GlobalOptions()[ALIstring("output_length_value_dimension")]);
  if (ad == 0)
    ad = ALIint(gomgr->GlobalOptions()[ALIstring("length_value_dimension")]);
  _OutputLengthValueDimensionFactor = CalculateLengthDimensionFactorFromInt(ad);

  ad = ALIint(gomgr->GlobalOptions()[ALIstring("output_length_error_dimension")]);
  if (ad == 0)
    ad = ALIint(gomgr->GlobalOptions()[ALIstring("length_error_dimension")]);
  _OutputLengthSigmaDimensionFactor = CalculateLengthDimensionFactorFromInt(ad);

  //---------- Change factor to convert to error dimensions
  //  _LengthValueDimensionFactor /= _LengthSigmaDimensionFactor;
  //_LengthSigmaDimensionFactor = 1;

  if (ALIUtils::debug >= 6)
    std::cout << _OutputLengthValueDimensionFactor << "Output Length Dimension Factors"
              << _OutputLengthSigmaDimensionFactor << std::endl;
}

static void ALIUtils::setReportVerbosity ( ALIint  val )

inlinestatic

Definition at line 37 of file ALIUtils.h.

References heppy_batch::val.

Referenced by GlobalOptionMgr::setGlobalOption().

{ report = val; }

ALIstring ALIUtils::subQuotes ( const ALIstring &  str )

static

Definition at line 439 of file ALIUtils.cc.

References beam_dqm_sourceclient-live_cfg::cerr, and beamvalidation::exit().

Referenced by val0().

                                                  {
  //---------- Take out leading and trailing '"'
  if (str.find('"') != 0 || str.rfind('"') != str.length() - 1) {
    std::cerr << "!!!EXITING trying to substract quotes from a word that has no quotes " << str << std::endl;
    exit(1);
  }

  //  str = str.strip(ALIstring::both, '\"');
  //---------- Take out leading and trallling '"'
  ALIstring strt = str.substr(1, str.size() - 2);

  //-  std::cout << " subquotes " << str << std::endl;
  //---------- Look for leading spaces
  while (strt[0] == ' ') {
    strt = strt.substr(1, strt.size() - 1);
  }

  //---------- Look for trailing spaces
  while (strt[strt.size() - 1] == ' ') {
    strt = strt.substr(0, strt.size() - 1);
  }

  return strt;
}

static time_t ALIUtils::time_now ( )

inlinestatic

Definition at line 39 of file ALIUtils.h.

References _time_now.

Referenced by Fit::getFitQuality(), Fit::multiplyMatrices(), and Fit::PropagateErrors().

{ return _time_now; }

static ALIdouble ALIUtils::val0 ( ALIdouble  val )

inlinestatic

Definition at line 71 of file ALIUtils.h.

References changeName(), getDimensionValue(), AlCaHLTBitMon_QueryRunRegistry::string, subQuotes(), and heppy_batch::val.

                                       {
    //-std::cout << val << " val " << ( (val <= 1.E-9) ? 0. : val) << std::endl;
    //    return (abs(val) <= 1.E-9) ? 0. : val; }
    if (std::fabs(val) <= 1.E-9) {
      return 0.;
    } else {
      return val;
    };
  }

Member Data Documentation

ALIdouble ALIUtils::_AngleSigmaDimensionFactor = 1.

staticprivate

Definition at line 106 of file ALIUtils.h.

Referenced by AngleSigmaDimensionFactor(), and SetAngleDimensionFactors().

ALIdouble ALIUtils::_AngleValueDimensionFactor = 1.

staticprivate

Definition at line 105 of file ALIUtils.h.

Referenced by AngleValueDimensionFactor(), and SetAngleDimensionFactors().

ALIdouble ALIUtils::_LengthSigmaDimensionFactor = 1.

staticprivate

Definition at line 104 of file ALIUtils.h.

Referenced by LengthSigmaDimensionFactor(), and SetLengthDimensionFactors().

ALIdouble ALIUtils::_LengthValueDimensionFactor = 1.

staticprivate

Definition at line 103 of file ALIUtils.h.

Referenced by LengthValueDimensionFactor(), and SetLengthDimensionFactors().

ALIdouble ALIUtils::_OutputAngleSigmaDimensionFactor = 1.

staticprivate

Definition at line 110 of file ALIUtils.h.

Referenced by dumpDimensions(), OutputAngleSigmaDimensionFactor(), and SetOutputAngleDimensionFactors().

ALIdouble ALIUtils::_OutputAngleValueDimensionFactor = 1.

staticprivate

Definition at line 109 of file ALIUtils.h.

Referenced by dumpDimensions(), OutputAngleValueDimensionFactor(), and SetOutputAngleDimensionFactors().

ALIdouble ALIUtils::_OutputLengthSigmaDimensionFactor = 1.

staticprivate

Definition at line 108 of file ALIUtils.h.

Referenced by dumpDimensions(), OutputLengthSigmaDimensionFactor(), and SetOutputLengthDimensionFactors().

ALIdouble ALIUtils::_OutputLengthValueDimensionFactor = 1.

staticprivate

Definition at line 107 of file ALIUtils.h.

Referenced by dumpDimensions(), OutputLengthValueDimensionFactor(), and SetOutputLengthDimensionFactors().

time_t ALIUtils::_time_now

staticprivate

Definition at line 111 of file ALIUtils.h.

Referenced by set_time_now(), and time_now().

ALIint ALIUtils::debug = -1

static

Definition at line 34 of file ALIUtils.h.

Referenced by Measurement::addAffectingEntriesFromOptO(), Fit::addDaMatrixToEntries(), Entry::addFittedDisplacementToValue(), ParameterMgr::addParameter(), ParameterMgr::addRandomFlatParameter(), ParameterMgr::addRandomGaussParameter(), ALILine::ALILine(), ALIPlane::ALIPlane(), CocoaDBMgr::BuildAlignments(), Model::BuildMeasurementsFromOA(), CocoaDaqReaderRoot::BuildMeasurementsFromOptAlign(), CocoaDBMgr::BuildOpticalAlignments(), Measurement::buildOptOList(), OpticalObject::buildRmFromEntryValuesOriginalOriginal(), OpticalObject::buildWordList(), OptOOpticalSquare::calculateFaces(), OpticalObject::calculateLocalRotationAxisInGlobal(), Measurement::calculateOriginalSimulatedValue(), Fit::calculateSimulatedMeasurementsWithOriginalValues(), MeasurementSensor2D::calculateSimulatedValue(), MeasurementDistancemeter3dim::calculateSimulatedValue(), MeasurementDistancemeter::calculateSimulatedValue(), MeasurementCOPS::calculateSimulatedValue(), MeasurementDiffEntry::calculateSimulatedValue(), MeasurementTiltmeter::calculateSimulatedValue(), Fit::CheckIfFitPossible(), Fit::CheckIfMeasIsProportionalToAnother(), Model::CMSLinkDeleteOptOs(), CocoaDaqReaderRoot::CocoaDaqReaderRoot(), CocoaDaqReaderText::CocoaDaqReaderText(), Measurement::construct(), OpticalObject::construct(), Measurement::constructFromOA(), OpticalObject::constructFromOptAligInfo(), OptOCOPS::convertPointToLocalCoordinates(), OptOSensor2D::convertPointToLocalCoordinates(), OpticalObject::copyData(), Measurement::copyMeas(), MeasurementSensor2D::correctValueAndSigma(), MeasurementCOPS::correctValueAndSigma(), MeasurementDistancemeter::correctValueAndSigma(), MeasurementDistancemeter3dim::correctValueAndSigma(), MeasurementTiltmeter::correctValueAndSigma(), OpticalObject::createComponentOptOs(), OpticalObject::createComponentOptOsFromOptAlignInfo(), Model::createCopyComponentList(), FittedEntriesManager::createFileName(), Fit::CreateMatrices(), OpticalObject::createNewOptO(), OptOLaser::defaultBehaviour(), OptOPinhole::defaultBehaviour(), OptOXLaser::defaultBehaviour(), OptOCOPS::defaultBehaviour(), OptOSource::defaultBehaviour(), Model::deleteOptO(), Measurement::DerivativeRespectEntry(), OptOPlateSplitter::detailedDeviatesLightRay(), OptOMirror::detailedDeviatesLightRay(), OptOCubeSplitter::detailedDeviatesLightRay(), OptOModifiedRhomboidPrism::detailedDeviatesLightRay(), OptOOpticalSquare::detailedDeviatesLightRay(), OptOPlateSplitter::detailedTraversesLightRay(), OptOCubeSplitter::detailedTraversesLightRay(), OptOModifiedRhomboidPrism::detailedTraversesLightRay(), OptOOpticalSquare::detailedTraversesLightRay(), OptOMirror::detailedTraversesLightRay(), OptOSensor2D::detailedTraversesLightRay(), EntryLengthAffCentre::displace(), Entry::displace(), OpticalObject::displaceCentreGlob(), OpticalObject::displaceCentreGlobOriginal(), OpticalObject::displaceCentreGlobOriginalOriginal(), EntryLengthAffCentre::displaceOriginal(), Entry::displaceOriginal(), EntryAngleAffAngles::displaceOriginalOriginal(), EntryLengthAffCentre::displaceOriginalOriginal(), Entry::displaceOriginalOriginal(), OpticalObject::displaceRmGlobAroundGlobal(), OpticalObject::displaceRmGlobAroundLocal(), OpticalObject::displaceRmGlobOriginal(), OpticalObject::displaceRmGlobOriginalOriginal(), rrapi.RRApi::dprint(), util.rrapi.RRApi::dprint(), Measurement::DumpBadOrderOptOs(), CocoaDBMgr::DumpCocoaResults(), FittedEntriesManager::dumpEntriesSubstraction(), Fit::dumpEntryAfterFit(), Fit::dumpEntryCorrelations(), Fit::dumpFittedValues(), MatrixMeschach::EliminateColumns(), MatrixMeschach::EliminateLines(), Entry::Entry(), OptOCubeSplitter::fastDeviatesLightRay(), OptOModifiedRhomboidPrism::fastDeviatesLightRay(), OptOOpticalSquare::fastDeviatesLightRay(), OptOPlateSplitter::fastDeviatesLightRay(), OptOMirror::fastDeviatesLightRay(), OptOCubeSplitter::fastTraversesLightRay(), OptOOpticalSquare::fastTraversesLightRay(), OptOModifiedRhomboidPrism::fastTraversesLightRay(), OptOSensor2D::fastTraversesLightRay(), OptOCOPS::fastTraversesLightRay(), OptOPlateSplitter::fastTraversesLightRay(), OptOMirror::fastTraversesLightRay(), EntryData::fill(), Entry::fill(), Model::fillCopyComponentList(), Measurement::fillData(), FittedEntriesSet::FillEntriesAveragingSets(), OptOSensor2D::fillExtraEntry(), OpticalObject::fillExtraEntry(), Entry::fillFromInputFileQuality(), Entry::fillFromInputFileSigma(), Entry::fillFromInputFileValue(), Entry::fillFromReportOutFileValue(), Fit::FillMatricesWithCalibratedParameters(), Fit::FillMatricesWithMeasurements(), EntryMgr::findEntryByLongName(), OpticalObject::findExtraEntryValue(), Fit::fitNextEvent(), Fit::fitParameters(), FittedEntry::FittedEntry(), pkg.AbstractPkg::generate(), rrapi.RRApi::get(), util.rrapi.RRApi::get(), pkg.AbstractPkg::get_kwds(), CocoaDBMgr::GetAlignInfoErrorFromOptO(), CocoaDBMgr::GetAlignInfoFromOptO(), Model::getComponentOptOs(), OpticalObject::getCoordinateFromOptAlignParam(), FittedEntriesManager::GetDifferentBetweenLasers(), OpticalObject::getDispVec(), Model::getEntryByName(), Fit::getEntryValue(), Fit::getFitQuality(), GlobalOptionMgr::getGlobalOption(), GlobalOptionMgr::getGlobalOptionValue(), Fit::getInstance(), OptOCOPS::getMeasFromInters(), OptOCubeSplitter::getMiddlePlate(), Model::getOptOByName(), OpticalObject::getPlate(), OptOModifiedRhomboidPrism::getRotatedPlate(), getRotationAnglesFromMatrix(), OpticalObject::getRotationAnglesFromMatrix(), OpticalObject::getRotationAnglesInOptOFrame(), Fit::GetSChi2(), OptOCubeSplitter::getUpperPlate(), ALILine::intersect(), LightRay::intersect(), MatrixMeschach::inverse(), LightRay::LightRay(), runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017(), FittedEntriesManager::MakeHistos(), OptOTiltmeter::makeMeasurement(), OptODistancemeter::makeMeasurement(), OptODistancemeter3dim::makeMeasurement(), OptOSensor2D::makeMeasurement(), OptOCOPS::makeMeasurement(), Fit::multiplyMatrices(), Model::nextOptOToCopy(), operator+(), MatrixMeschach::operator=(), OpticalObject::OpticalObject(), Measurement::postConstruct(), Fit::printCentreInOptOFrame(), Fit::PrintChi2(), Fit::PropagateErrors(), OpticalObject::readCoordinates(), EntryMgr::readEntryFromReportOut(), CocoaDaqReaderRoot::ReadEvent(), DeviationsFromFileSensor2D::readFile(), FittedEntriesReader::readFittedEntriesFromFile(), ErrorCorrelationMgr::readFromReportFile(), Model::readMeasurementsFromFile(), CocoaDaqReaderText::ReadNextEvent(), Model::readSystemDescription(), Fit::redoMatrices(), LightRay::reflect(), LightRay::refract(), OpticalObject::rotateItAroundGlobal(), runTauIdMVA.TauIDEmbedder::runTauID(), SetAngleDimensionFactors(), OpticalObject::SetCentreGlobFromCentreLocal(), OpticalObject::SetCentreLocalFromEntryValues(), MatrixMeschach::SetCorrelation(), Fit::setCorrelationsInWMatrix(), Fit::setFittableEntries(), GlobalOptionMgr::setGlobalOption(), OpticalObject::setGlobalRMOriginalOriginal(), SetLengthDimensionFactors(), FittedEntriesSet::SetOptOEntries(), OpticalObject::setOriginalEntryValues(), SetOutputAngleDimensionFactors(), SetOutputLengthDimensionFactors(), OpticalObject::SetRMGlobFromRMLocal(), OpticalObject::SetRMLocalFromEntryValues(), Model::SetValueDisplacementsFromReportOut(), LightRay::shiftAndDeviateWhileTraversing(), Fit::startFit(), LightRay::startLightRay(), crabFunctions.CrabController::submit(), Measurement::Substitute2p(), FittedEntriesReader::substitutePointBySlash(), Fit::substractLastDisplacementToEntries(), OpticalObject::transformCylindrical2Cartesian(), OptOUserDefined::userDefinedBehaviour(), EntryAngleAffAngles::valueDisplaced(), EntryLengthAffCentre::valueDisplaced(), Entry::valueDisplaced(), pkg.AbstractPkg::write(), and Fit::WriteVisualisationFiles().

ALIdouble ALIUtils::deg = 0.017453293

static

Definition at line 35 of file ALIUtils.h.

Referenced by OptOModifiedRhomboidPrism::getRotatedPlate().

ALIbool ALIUtils::firstTime = false

staticprivate

Definition at line 113 of file ALIUtils.h.

Referenced by getFirstTime(), and setFirstTime().

ALIdouble ALIUtils::maximum_deviation_derivative = 1.E-6

staticprivate

Definition at line 115 of file ALIUtils.h.

Referenced by getMaximumDeviationDerivative(), and setMaximumDeviationDerivative().

ALIint ALIUtils::report = 1

static

Definition at line 33 of file ALIUtils.h.

Referenced by Fit::dumpEntryAfterFit(), Fit::fitNextEvent(), Fit::getFitQuality(), dataset.Dataset::getPrimaryDatasetEntries(), Fit::PrintChi2(), runall.testit::run(), and addOnTests.testit::run().