Measurement.cc

Go to the documentation of this file.

// COCOA class implementation file
// Id:  Measurement.cc
// CAT: Model
// ---------------------------------------------------------------------------
// History: v1.0
// Authors:
//   Pedro Arce

#include "Alignment/CocoaModel/interface/Model.h"

#include <algorithm>
#include <iomanip>
#include <iostream>
#include <iterator>
//#include <algo.h>
#include <cstdlib>
#include <cmath>  // include floating-point std::abs functions

#include "Alignment/CocoaModel/interface/Entry.h"
#include "Alignment/CocoaUtilities/interface/ALIFileIn.h"
#include "Alignment/CocoaUtilities/interface/GlobalOptionMgr.h"
#include "Alignment/CocoaModel/interface/Measurement.h"
#include "Alignment/CocoaModel/interface/OpticalObject.h"
#include "Alignment/CocoaModel/interface/ParameterMgr.h"
#ifdef COCOA_VIS
#include "Alignment/CocoaVisMgr/interface/ALIVRMLMgr.h"
#include "Alignment/IgCocoaFileWriter/interface/IgCocoaFileMgr.h"
#endif
#include "CondFormats/OptAlignObjects/interface/OpticalAlignMeasurementInfo.h"
#include "CondFormats/OptAlignObjects/interface/OpticalAlignInfo.h"

ALIdouble Measurement::cameraScaleFactor = 1.;
ALIstring Measurement::theMeasurementsFileName = "";
ALIstring Measurement::theCurrentDate = "99/99/99";
ALIstring Measurement::theCurrentTime = "99:99";

ALIbool Measurement::only1 = false;
ALIstring Measurement::only1Date = "";
ALIstring Measurement::only1Time = "";

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@ constructor:
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
Measurement::Measurement(const ALIint measdim, ALIstring& type, ALIstring& name)
    : theDim(measdim), theType(type), theName(name) {
  //  _OptOnames = new ALIstring[theDim];
  theValue = new ALIdouble[theDim];
  theSigma = new ALIdouble[theDim];
  theValueType = new ALIstring[theDim];

  theValueSimulated = new ALIdouble[theDim];
  theValueSimulated_orig = new ALIdouble[theDim];
  theValueIsSimulated = new ALIbool[theDim];
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@ construct (read from file)
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void Measurement::construct() {
  ALIFileIn& filein = ALIFileIn::getInstance(Model::SDFName());

  //---------- Read OptOs that take part in this Measurement
  std::vector<ALIstring> wordlist;
  filein.getWordsInLine(wordlist);

  //--------- Fill the list of names of OptOs that take part in this measurement ( names only )
  buildOptONamesList(wordlist);

  if (ALIUtils::debug >= 3) {
    std::cout << "@@@@ Reading Measurement " << name() << " TYPE= " << type() << std::endl << " MEASURED OPTO NAMES: ";
    std::ostream_iterator<ALIstring> outs(std::cout, " ");
    copy(wordlist.begin(), wordlist.end(), outs);
    std::cout << std::endl;
  }

  //---------- Read the data
  for (unsigned int ii = 0; ii < dim(); ii++) {
    filein.getWordsInLine(wordlist);
    fillData(ii, wordlist);
  }

  if (!valueIsSimulated(0))
    correctValueAndSigma();

  postConstruct();
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void Measurement::constructFromOA(OpticalAlignMeasurementInfo& measInfo) {
  //---- Build wordlist to build object name list
  std::vector<std::string> objNames = measInfo.measObjectNames_;
  std::vector<std::string>::const_iterator site;
  std::vector<ALIstring> wordlist;
  //--- Fill the list of names of OptOs that take part in this measurement ( names only )
  for (site = objNames.begin(); site != objNames.end(); ++site) {
    if (site != objNames.begin())
      wordlist.push_back("&");
    wordlist.push_back(*site);
  }
  buildOptONamesList(wordlist);

  if (ALIUtils::debug >= 3) {
    std::cout << "@@@@ Reading Measurement " << name() << " TYPE= " << type() << " " << measInfo << std::endl
              << " MEASURED OPTO NAMES: ";
    for (size_t ii = 0; ii < _OptONameList.size(); ii++) {
      std::cout << _OptONameList[ii] << " ";
    }
    std::cout << std::endl;
  }

  //---------- No data, set to simulated_value
  for (unsigned int ii = 0; ii < dim(); ii++) {
    wordlist.clear();
    wordlist.push_back((measInfo.values_)[ii].name_);
    char ctmp[20];
    if (measInfo.isSimulatedValue_[ii]) {
      if (ALIUtils::debug >= 5) {
        std::cout << "Measurement::constructFromOA:  meas value " << ii << " " << dim() << " = simulated_value"
                  << std::endl;
      }
      wordlist.push_back("simulated_value");
    } else {
      if (ALIUtils::debug >= 5) {
        std::cout << "Measurement::constructFromOA:  meas value " << ii << " " << dim() << " = "
                  << measInfo.values_.size() << std::endl;
      }
      ALIdouble val =
          (measInfo.values_)[ii].value_ /
          valueDimensionFactor();  //in XML  values are without dimensions, so neutralize multiplying by valueDimensionFactor() in fillData
      gcvt(val, 10, ctmp);
      wordlist.push_back(ctmp);
    }
    ALIdouble err =
        (measInfo.values_)[ii].error_ /
        sigmaDimensionFactor();  //in XML  values are without dimensions, so neutralize multiplying by valueDimensionFactor() in fillData
    gcvt(err, 10, ctmp);
    wordlist.push_back(ctmp);
    std::cout << " sigma " << err << " = " << ctmp << " " << (measInfo.values_)[ii].error_ << std::endl;
    //-    wordlist.push_back( "simulated_value" );
    //-   wordlist.push_back( "1." );
    if (ALIUtils::debug >= 5)
      ALIUtils::dumpVS(wordlist, " Measurement: calling fillData ");
    //-    std::cout << " MEAS INFO " << measInfo << std::endl;
    //-   std::cout << ii << " MEAS INFO PARAM " <<  (measInfo.values_)[ii] << std::endl;
    //- std::cout << ii << " MEAS INFO PARAM VALUE " <<  (measInfo.values_)[ii].value_ << std::endl;
    fillData(ii, wordlist);
  }

  postConstruct();
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void Measurement::postConstruct() {
  //---------- Set name as name of last OptO
  setName();

  //---------- Transform for each Measurement the Measured OptO names to Measured OptO pointers
  buildOptOList();

  //---------- Build list of Entries that affect a Measurement
  buildAffectingEntryList();

  //---------- add this measurement to the global list of measurements
  Model::addMeasurementToList(this);

  if (ALIUtils::debug >= 10) {
    std::cout << Model::MeasurementList().size() << std::endl;
  }
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@ Fills a list of names of OpticalObjects that take part in this measurement
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void Measurement::buildOptONamesList(const std::vector<ALIstring>& wl) {
  int NPairs = (wl.size() + 1) / 2;  // Number of OptO names ( pair of name and '&' )

  //--------- Fill list with names
  for (int ii = 0; ii < NPairs; ii++) {
    _OptONameList.push_back(wl[ii * 2]);
    // Check for separating '&'
    if (ii != NPairs - 1 && wl[2 * ii + 1] != ALIstring("&")) {
      //      ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
      std::cerr << "!!! Measured Optical Objects should be separated by '&', not by" << wl[2 * ii + 1] << std::endl;
      exit(2);
    }
  }
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@ Fill the data of measurement coordinate 'coor' with values in 'wordlist'
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void Measurement::fillData(ALIuint coor, const std::vector<ALIstring>& wordlist) {
  if (ALIUtils::debug >= 3) {
    std::cout << "@@ Filling coordinate " << coor << std::endl;
    //-   ostream_iterator<ALIstring> outs(std::cout," ");
    //-  copy(wordlist.begin(), wordlist.end(), outs);
  }

  ParameterMgr* parmgr = ParameterMgr::getInstance();

  //---------- Check that there are 3 attributes: name, value, error
  if (wordlist.size() != 3) {
    //    ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
    std::cerr << " Incorrect format for Measurement value:" << std::endl;
    std::ostream_iterator<ALIstring> outs(std::cout, " ");
    copy(wordlist.begin(), wordlist.end(), outs);
    std::cout << std::endl << "There should be three words: name value sigma " << std::endl;
    exit(2);
  }

  //---------- check coor value
  if (coor >= theDim) {
    // ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
    std::cerr << "Trying to fill Measurement coordinate No " << coor << " but the dimension is " << theDim << std::endl;
    exit(2);
  }

  //---------- set data members
  //----- Set valueType
  theValueType[coor] = wordlist[0];

  //----- Set value (translate it if a PARAMETER is used)
  ALIdouble val = 0.;
  theValueIsSimulated[coor] = false;
  if (!ALIUtils::IsNumber(wordlist[1])) {
    if (parmgr->getParameterValue(wordlist[1], val) == 0) {
      if (wordlist[1] == ALIstring("simulated_value")) {
        theValueIsSimulated[coor] = true;
      } else {
        //      ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
        std::cerr << "!!! parameter for value not found: " << wordlist[1].c_str() << std::endl;
        exit(2);
      }
    }
    //d val *= valueDimensionFactor();
  } else {
    //d val = DimensionMgr()::getInstance()->extractValue( wordlist[1], ValueDimensionFactor() );
    val = atof(wordlist[1].c_str());
  }
  val *= valueDimensionFactor();
  if (ALIUtils::debug >= 3)
    std::cout << "Meas VALUE= " << val << " (ValueDimensionFactor= " << valueDimensionFactor() << std::endl;

  //----- Set sigma (translate it if a PARAMETER is used)
  ALIdouble sig = 0.;
  if (!ALIUtils::IsNumber(wordlist[2])) {
    if (parmgr->getParameterValue(wordlist[2], sig) == 0) {
      // ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
      std::cerr << "!!! parameter for sigma not found: " << wordlist[2].c_str() << std::endl;
      exit(2);
    }
    //d sig *= sigmaDimensionFactor();
  } else {
    //    sig = DimensionMgr()::getInstance()->extractValue( wordlist[2], ValueDimensionFactor() );
    sig = atof(wordlist[2].c_str());
  }
  sig *= sigmaDimensionFactor();
  if (ALIUtils::debug >= 3)
    std::cout << "SIGMA= " << sig << " (SigmaDimensionFactor= " << sigmaDimensionFactor() << std::endl;

  //----- set theValue & theSigma
  theValue[coor] = val;
  theSigma[coor] = sig;
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void Measurement::fillData(ALIuint coor, OpticalAlignParam* oaParam) {
  if (ALIUtils::debug >= 3) {
    std::cout << "@@ Filling coordinate " << coor << std::endl;
  }

  //  ParameterMgr* parmgr = ParameterMgr::getInstance();

  //---------- check coor value
  if (coor >= theDim) {
    std::cerr << "Trying to fill Measurement coordinate No " << coor << " but the dimension is " << theDim << std::endl;
    exit(2);
  }

  //---------- set data members
  //----- Set value (translate it if a PARAMETER is used)
  ALIdouble val = 0.;
  theValueIsSimulated[coor] = false;
  val = oaParam->value();
  val *= valueDimensionFactor();
  theValue[coor] = val;
  if (ALIUtils::debug >= 3)
    std::cout << "Meas VALUE= " << val << " (ValueDimensionFactor= " << valueDimensionFactor() << std::endl;

  ALIbool sigmaFF = GlobalOptionMgr::getInstance()->GlobalOptions()["measurementErrorFromFile"];
  if (sigmaFF) {
    //----- Set sigma (translate it if a PARAMETER is used)
    ALIdouble sig = 0.;
    sig = oaParam->sigma();  // it is in mm always
    sig *= sigmaDimensionFactor();
    theSigma[coor] = sig;
    if (ALIUtils::debug >= 3)
      std::cout << "SIGMA= " << sig << " (SigmaDimensionFactor= " << sigmaDimensionFactor() << std::endl;
  }
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@ Once the complete list of OptOs is read, convert OptO names to pointers
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void Measurement::buildOptOList() {
  //-  if ( ALIUtils::debug >= 3 ) std::cout << std::endl << " MEASUREMENT: " << " " << this->name() << std::endl;
  ALIstring twopoints("..");  // .. goes one level up in the tree fo OptOs

  //---------------------------------------- Loop OptONameList
  std::vector<ALIstring>::iterator vsite;
  for (vsite = _OptONameList.begin(); vsite != _OptONameList.end(); ++vsite) {
    //----------------------------------- Count how many '..' there are in the name
    ALIuint ii = 0;
    //    ALIuint slen = (*vsite).length();
    ALIuint Ntwopoints = 0;  //---- No '..' in ALIstring
    for (;;) {
      int i2p = (*vsite).find_first_of(twopoints, 3 * ii);  // if it is ., it also finds it!!!
      if (i2p < 0)
        break;
      if (i2p != ALIint(3 * ii)) {
        std::cerr << i2p << "!!! Bad position of '..' in reference ALIstring: " << (*vsite).c_str() << std::endl;
        exit(2);
      } else {
        Ntwopoints++;
        if (ALIUtils::debug >= 9)
          std::cout << "N2p" << Ntwopoints;
      }
      ii++;
    }
    //------ Substitute '..' by reference (the last OptO in list)
    if (Ntwopoints != 0) {
      Substitute2p((*vsite), *(_OptONameList.end() - 1), Ntwopoints);
    }
    //----- Get OpticalObject* that correspond to ALIstring and fill list
    ALIstring referenceOptO = (*vsite);
    //--- a ':' is used in OptOs that have several possible behavious
    ALIint colon = referenceOptO.find(':');
    if (colon != -1) {
      if (ALIUtils::debug >= 99) {
        std::cout << "colon in reference OptO name " << colon << referenceOptO.c_str() << std::endl;
      }
      referenceOptO = referenceOptO.substr(0, colon);
    }
    OpticalObject* OptOitem = Model::getOptOByName(referenceOptO);
    if (ALIUtils::debug >= 3)
      std::cout << "Measurement::buildOptOList: OptO in Measurement: " << OptOitem->name() << std::endl;
    if (OptOitem != (OpticalObject*)nullptr) {
      _OptOList.push_back(OptOitem);
    } else {
      std::cerr << "!!! Error in Measurement: can't find Optical Object " << (*vsite).c_str() << std::endl;
      exit(2);
    }
  }
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@ Build the list of all entries of every OptO that take part in this
//@@ Measurement and also the list of all entries of their OptO ancestors
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void Measurement::buildAffectingEntryList() {
  //---------- Loop OptO MeasuredList
  std::vector<OpticalObject*>::const_iterator vocite;
  for (vocite = _OptOList.begin(); vocite != _OptOList.end(); ++vocite) {
    addAffectingEntriesFromOptO(*vocite);
  }
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@ Get the list of all entries of this OptO that take part in this Measurement
//@@ and of their OptO ancestors
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void Measurement::addAffectingEntriesFromOptO(const OpticalObject* optoP) {
  if (ALIUtils::debug >= 3)
    std::cout << "Measurement::addAffectingEntriesFromOptO: OptO taking part in Measurement: " << optoP->name()
              << std::endl;
  //---------- Loop entries in this OptO
  std::vector<Entry*>::const_iterator vecite;
  std::vector<Entry*>::const_iterator fvecite;
  for (vecite = optoP->CoordinateEntryList().begin(); vecite != optoP->CoordinateEntryList().end(); ++vecite) {
    //T     if( find( theAffectingEntryList.begin(), theAffectingEntryList.end(), (*vecite) ) == theAffectingEntryList.end() ){
    //t      theAffectingEntryList.push_back(*vecite);
    //T    }
    fvecite = find(theAffectingEntryList.begin(), theAffectingEntryList.end(), (*vecite));
    if (fvecite == theAffectingEntryList.end()) {
      theAffectingEntryList.push_back(*vecite);
      if (ALIUtils::debug >= 4)
        std::cout << "Entry that may affect Measurement: " << (*vecite)->name() << std::endl;
    }
  }
  for (vecite = optoP->ExtraEntryList().begin(); vecite != optoP->ExtraEntryList().end(); ++vecite) {
    fvecite = find(theAffectingEntryList.begin(), theAffectingEntryList.end(), (*vecite));
    if (fvecite == theAffectingEntryList.end()) {
      theAffectingEntryList.push_back(*vecite);
      if (ALIUtils::debug >= 4)
        std::cout << "Entry that may affect Measurement: " << (*vecite)->name() << std::endl;
    }
  }
  if (optoP->parent() != nullptr) {
    addAffectingEntriesFromOptO(optoP->parent());
  }
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@ Substitute '..' in 'ref' by name of parent OptO ('firstref')
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void Measurement::Substitute2p(ALIstring& ref, const ALIstring& firstref, int Ntwopoints) {
  // '/' sets hierarchy of OptOs
  ALIstring slash("/");

  int pos1st = firstref.length();
  // Go back an '/' in firstref for each '..' in ref
  for (int ii = 0; ii < Ntwopoints; ii++) {
    pos1st = firstref.find_last_of(slash, pos1st - 1);
    if (ALIUtils::debug >= 9)
      std::cout << "pos1st=" << pos1st;
  }

  if (ALIUtils::debug >= 9)
    std::cout << "before change ref: " << ref << " 1ref " << firstref << std::endl;
  // Substitute name
  ref.replace(0, (Ntwopoints * 3) - 1, firstref, 0, pos1st);
  if (ALIUtils::debug >= 9)
    std::cout << "after change ref: " << ref << " 1ref " << firstref << std::endl;
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void Measurement::printStartCalculateSimulatedValue(const Measurement* meas) {
  std::cout << std::endl
            << "@@@ Start calculation of simulated value of " << meas->type() << " Measurement " << meas->name()
            << std::endl;
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@ Calculate the simulated value of the Measurement propagating the LightRay when all the entries have their original values
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void Measurement::calculateOriginalSimulatedValue() {
  //----------  Calculate the simulated value of the Measurement
  calculateSimulatedValue(true);

#ifdef COCOA_VIS
  if (ALIUtils::getFirstTime()) {
    GlobalOptionMgr* gomgr = GlobalOptionMgr::getInstance();
    if (gomgr->GlobalOptions()["VisWriteVRML"] > 1) {
      ALIVRMLMgr::getInstance().newLightRay();
    }
    /*-    if(Model::GlobalOptions()["VisWriteIguana"] > 1) {
           IgCocoaFileMgr::getInstance().newLightPath( theName );
           } */
  }
#endif

  //---------- Set original simulated values to it
  //-  if(ALIUtils::debug >= 5) std::cout << "MEAS DIMENSION" << dim() << std::endl;
  for (ALIuint ii = 0; ii < dim(); ii++) {
    setValueSimulated_orig(ii, valueSimulated(ii));
    if (ALIUtils::debug >= 4)
      std::cout << "SETsimuvalOriginal" << valueSimulated(ii) << std::endl;
    //----- If Measurement has as value 'simulated_value', set the value to the simulated one
    if (valueIsSimulated(ii) == 1) {
      setValue(ii, valueSimulated(ii));
      //- std::cout << ii << " setting value as simulated " <<  valueSimulated(ii) << " " << value(ii) << this << std::endl;
    }
  }
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void Measurement::DumpBadOrderOptOs() {
  std::cerr << " Detector can not make measurement with these optical objects " << std::endl;
  if (ALIUtils::debug >= 1) {
    //  std::vector<OpticalObject*>::iterator voite;
    //      for ( voite = _OptOList.begin();
    //     voite != _OptOList.end(); voite++) {
    //    std::cout << (*voite)->type() << " : " << (*voite)->name() << std::endl;
    // }
    std::vector<ALIstring>::const_iterator vsite;
    for (vsite = OptONameList().begin(); vsite != OptONameList().end(); ++vsite) {
      std::cerr << (*vsite) << " : ";
    }
    std::cerr << std::endl;
  }
  exit(2);
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
std::vector<ALIdouble> Measurement::DerivativeRespectEntry(Entry* entry) {
  //---------- std::vector of derivatives to return
  std::vector<ALIdouble> deriv;
  ALIdouble sumderiv;

  //---------- displacement to start with
  ALIdouble displacement = entry->startingDisplacement();
  //----- all angles are in radians, so, if displace is not, rescale it before making the displacement
  //-  displacement *= entry->SigmaDimensionFactor();
  if (ALIUtils::debug >= 3)
    std::cout << std::endl
              << "%%% Derivative w.r.t. entry " << entry->name() << ": displacement = " << displacement << std::endl;

  ALIint count_itera = 0;

  //---------- Loop decreasing the displacement a factor 2, until the precision set is reached
  do {
    count_itera++;
    entry->displace(displacement);

    if (ALIUtils::debug >= 5)
      std::cout << "Get simulated value for displacement " << displacement << std::endl;
    calculateSimulatedValue(false);

    //---------- Get sum of derivatives
    sumderiv = 0;
    for (ALIuint ii = 0; ii < theDim; ii++) {
      sumderiv += std::abs(theValueSimulated[ii] - theValueSimulated_orig[ii]);
      if (ALIUtils::debug >= 4) {
        std::cout << "iteration " << count_itera << " COOR " << ii << " difference ="
                  << (theValueSimulated[ii] - theValueSimulated_orig[ii])
                  //-             << "  " << theValueSimulated[ii] << "  " << theValueSimulated_orig[ii]
                  << " derivative = " << (theValueSimulated[ii] - theValueSimulated_orig[ii]) / displacement << " disp "
                  << displacement << " sum derivatives = " << sumderiv << std::endl;
      }
      if (ALIUtils::debug >= 5) {
        std::cout << " new simu value= " << theValueSimulated[ii] << " orig simu value " << theValueSimulated_orig[ii]
                  << std::endl;
      }
    }
    if (count_itera >= 100) {
      std::cerr << "EXITING: too many iterations in derivative, displacement is " << displacement
                << " sum of derivatives is " << sumderiv << std::endl;
      exit(3);
    }
    displacement /= 2.;
    //-    std::cout << "sumderiv " << sumderiv << " maximu " <<  Fit::maximum_deviation_derivative << std::endl;
  } while (sumderiv > ALIUtils::getMaximumDeviationDerivative());
  displacement *= 2;

  //---------- Enough precision reached: pass result
  for (ALIuint ii = 0; ii < theDim; ii++) {
    deriv.push_back((theValueSimulated[ii] - theValueSimulated_orig[ii]) / displacement);
    //----- change it to entry sigma dimensions
    //     deriv[ii] /= entry->SigmaDimensionFactor();
    if (ALIUtils::debug >= 1)
      std::cout << name() << ": " << entry->OptOCurrent()->name() << " " << entry->name() << " " << ii
                << "### DERIVATIVE: " << deriv[ii] << std::endl;
  }
  //-  if(ALIUtils::debug >= 5) std::cout << "END derivative: " << deriv << "disp" << displacement << std::endl;

  //--------------------- Reset _centreGlob and _rmGlob of OptO entry belongs to (and component OptOs)
  entry->OptOCurrent()->resetGlobalCoordinates();

  return deriv;
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@ destructor
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
Measurement::~Measurement() {
  //  delete[] _name;
  delete[] theValue;
  delete[] theSigma;
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@ get the ':X' that determines how the behaviour of the OptO w.r.t. this Measurement
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
ALIstring Measurement::getMeasuringBehaviour(const std::vector<OpticalObject*>::const_iterator vocite) {
  std::vector<ALIstring>::const_iterator vscite =
      _OptONameList.begin() + (vocite - _OptOList.begin());  // point to corresponding name of this OptO
  ALIint colon = (*vscite).find(':');
  ALIstring behav;
  if (colon != -1) {
    behav = (*vscite).substr(colon + 1, (*vscite).size());
  } else {
    behav = " ";
  }
  return behav;
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@ Get the previous OptOs in the list of OptO that take part in this measurement
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
const OpticalObject* Measurement::getPreviousOptO(const OpticalObject* Popto) const {
  //--------- Loop OptOs that take part in this measurement
  std::vector<OpticalObject*>::const_iterator vocite;
  for (vocite = _OptOList.begin(); vocite != _OptOList.end(); ++vocite) {
    if (*vocite == Popto) {
      if (vocite == _OptOList.begin()) {
        std::cerr << " ERROR in  getPreviousOptO of measurement " << name() << std::endl;
        std::cerr << " OptO " << Popto->name() << " is the first one " << std::endl;
        exit(1);
      } else {
        return *(vocite - 1);
      }
    }
  }

  std::cerr << " ERROR in  getPreviousOptO of measurement " << name() << std::endl;
  std::cerr << " OptO " << Popto->name() << " not found " << std::endl;
  exit(1);
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void Measurement::setCurrentDate(const std::vector<ALIstring>& wl) {
  if (wl.size() != 3) {
    std::cerr << "!!!EXITING: reading DATE of measurements set: it must have three words, it is though " << std::endl;
    ALIUtils::dumpVS(wl, " ");
    exit(1);
  } else if (wl[0] != "DATE:") {
    std::cerr << "!!!EXITING: reading DATE of measurements set: first word must be 'DATE:', it is though " << std::endl;
    ALIUtils::dumpVS(wl, " ");
    exit(1);
  } else {
    theCurrentDate = wl[1];
    theCurrentTime = wl[2];
  }
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void Measurement::copyMeas(Measurement* meas, const std::string& subsstr1, const std::string& subsstr2) {
  theDim = meas->dim();
  theType = meas->type();
  theName = ALIUtils::changeName(meas->name(), subsstr1, subsstr2);

  //  _OptOnames = new ALIstring[theDim];
  theValueSimulated = new ALIdouble[theDim];
  theValueSimulated_orig = new ALIdouble[theDim];
  theValueIsSimulated = new ALIbool[theDim];
  theValue = const_cast<ALIdouble*>(meas->value());
  theSigma = const_cast<ALIdouble*>(meas->sigma());

  unsigned int ii;
  for (ii = 0; ii < theDim; ii++) {
    theValueSimulated[ii] = meas->valueSimulated(ii);
    theValueSimulated_orig[ii] = meas->valueSimulated_orig(ii);
    theValueIsSimulated[ii] = meas->valueIsSimulated(ii);
  }

  //--------- Fill the list of names of OptOs that take part in this measurement ( names only )

  std::vector<std::string> wordlist;
  auto& optolist = meas->OptOList();
  ALIuint nOptos = optolist.size();
  for (ALIuint ii = 0; ii < nOptos; ii++) {
    wordlist.push_back(ALIUtils::changeName(optolist[ii]->longName(), subsstr1, subsstr2));
    std::cout << " copymeas " << ALIUtils::changeName(optolist[ii]->longName(), subsstr1, subsstr2) << std::endl;
    if (ii != nOptos - 1)
      wordlist.push_back("&");
  }

  buildOptONamesList(wordlist);

  if (ALIUtils::debug >= 3) {
    std::cout << "@@@@ Reading Measurement " << name() << " TYPE= " << type() << std::endl << " MEASURED OPTO NAMES: ";
    std::ostream_iterator<ALIstring> outs(std::cout, " ");
    copy(wordlist.begin(), wordlist.end(), outs);
    std::cout << std::endl;
  }

  postConstruct();
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void Measurement::setName() {
  // name already set by passing one argument with sensor type
  if (!theName.empty())
    return;
  if (_OptONameList.empty()) {
    std::cerr
        << " !!! Error in your code, you cannot ask for the name of the Measurement before the OptONameList is build "
        << std::endl;
    exit(9);
  }
  std::vector<ALIstring>::iterator vsite = (_OptONameList.end() - 1);
  theName = type() + ":" + (*vsite);
}