MeasurementDiffEntry.cc

Go to the documentation of this file.

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

#include "Alignment/CocoaModel/interface/MeasurementDiffEntry.h"
#include "Alignment/CocoaModel/interface/Model.h"
#include "Alignment/CocoaModel/interface/OpticalObject.h"
#include "Alignment/CocoaModel/interface/Entry.h"
#include "Alignment/CocoaUtilities/interface/ALIFileIn.h"
#include <iostream>
#include <iomanip>
#include <cstdlib>
#include <cmath>        // include floating-point std::abs functions
#ifdef COCOA_VIS
#include "Alignment/CocoaVisMgr/interface/ALIVRMLMgr.h"
#include "Alignment/IgCocoaFileWriter/interface/IgCocoaFileMgr.h"
#include "Alignment/IgCocoaFileWriter/interface/ALIVisLightPath.h"
#endif


//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void MeasurementDiffEntry::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++ ) {
    //--- take out Entry names from object names
    int isl = wl[ii*2].rfind("/");
    AddOptONameListItem( wl[ii*2].substr(0,isl) );
    // 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 entry names
    if( ii == 0 ){
      theEntryNameFirst = wl[ii*2].substr(isl+1,999);
    }else if( ii == 1 ){
      theEntryNameSecond = wl[ii*2].substr(isl+1,999);
    }

  }

}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@ calculate the simulated value propagating the light ray through the OptO that take part in the Measurement
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void MeasurementDiffEntry::calculateSimulatedValue( ALIbool firstTime )
{

  if( ALIUtils::debug >= 2) printStartCalculateSimulatedValue( this ); // important for Examples/FakeMeas

  //---------- Loop list of OptO that take part in measurement
  std::vector<OpticalObject*>::const_iterator vocite =  OptOList().begin();
  if( ALIUtils::debug >= 5) std::cout  << "OptOList size" <<OptOList().size() << std::endl;

  //----- Check that there are only two objects
  if( OptOList().size() == !2 ) {
    std::cerr << "!!! ERROR in MeasurementDiffEntry: " << name() << " There should only be two objects " << std::endl;
    std::cerr << " 1st " << (*vocite)->name() << " 2nd " << (*vocite+1)->name()  << std::endl;
    DumpBadOrderOptOs();
    std::exception();
  }

#ifdef COCOA_VIS
  ALIVisLightPath* vispath = 0;
  if( ALIUtils::getFirstTime() ) {
    GlobalOptionMgr* gomgr = GlobalOptionMgr::getInstance();
    if(gomgr->GlobalOptions()["VisWriteIguana"] > 1) {
      vispath = IgCocoaFileMgr::getInstance().newLightPath( name() );
    }
  }
#endif

  //--- This is a special 'measurement': it represents the fact that you have measured two entries one relative to the other (e.g. relative angle of two objects)
  Entry* entry1 = Model::getEntryByName( (*(OptOList().begin()))->longName(), theEntryNameFirst );
  Entry* entry2 = Model::getEntryByName( (*(OptOList().begin()+1))->longName(), theEntryNameSecond );
  if( ALIUtils::debug >= 5) std::cout << "  entry1 " << (*(OptOList().begin()))->longName() << "/" << entry1->name() << " ->valueDisplaced() " << entry1->valueDisplaced() << " entry2 " <<  (*(OptOList().begin()+1))->longName() << "/" << entry2->name() << " ->valueDisplaced() " << entry2->valueDisplaced() << std::endl;
  setValueSimulated( 0, entry1->valueDisplaced() - entry2->valueDisplaced() );

  if (ALIUtils::debug >= 2) {
    ALIdouble detD = 1000*valueSimulated(0); if(std::abs(detD) <= 1.e-9 ) detD = 0.;
    std::cout << "REAL value: " <<"D: " << 1000.*value()[0] << " (mm)  " << (this)->name()
          << "   DIFF= " << detD-1000*value()[0] << std::endl;
    std::cout << "SIMU value: " << "D: " << detD << " (mm)  " << (this)->name() << std::endl;
  }

  if(ALIUtils::debug >= 5) std::cout << "end calculateSimulatedValue" <<std::endl;

}