/data/refman/pasoursint/CMSSW_5_3_10/src/Alignment/CocoaModel/src/FittedEntriesReader.cc

Go to the documentation of this file.

#include "Alignment/CocoaModel/interface/FittedEntriesReader.h"
#include "Alignment/CocoaModel/interface/Model.h"
#include "Alignment/CocoaModel/interface/Entry.h"
#include "Alignment/CocoaModel/interface/OpticalObject.h"
#include "Alignment/CocoaModel/interface/ALIRmDataFromFile.h"
#include "Alignment/CocoaUtilities/interface/ALIUtils.h"
  

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
FittedEntriesReader::FittedEntriesReader( const ALIstring& filename ) 
{
  theFileName = filename;
  theFile = ALIFileIn::getInstance( filename );
  std::vector<ALIstring> wl;
  theFile.getWordsInLine( wl );
  if (wl[0] == ALIstring( "DIMENSIONS:" ) ){
    theLengthDim = ALIUtils::CalculateLengthDimensionFactorFromString( wl[3] );
    theLengthErrorDim = ALIUtils::CalculateLengthDimensionFactorFromString( wl[5] );
    theAngleDim = ALIUtils::CalculateAngleDimensionFactorFromString( wl[8] );
    theAngleErrorDim = ALIUtils::CalculateAngleDimensionFactorFromString( wl[10] );
  } else {
    ALIUtils::dumpVS( wl, "!!! FATAL ERROR FittedEntriesReader: first line is not dimensions " );
    std::exception();
  }

}


//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
ALIbool FittedEntriesReader::readFittedEntriesFromFile()
{
  if( ALIUtils::debug >= 5) std::cout << " readFittedEntriesFromFile " << theFileName << std::endl;
  std::map<OpticalObject*,ALIRmDataFromFile> affAngles;

  std::vector<ALIstring> wl;
  theFile.getWordsInLine( wl );
  unsigned int siz = wl.size();
  for( size_t ii = 1; ii < siz; ii+=3 ){
    ALIstring optOentryName = substitutePointBySlash( wl[ii] );
    Entry* entry = Model::getEntryByName( optOentryName );
    if( ALIUtils::debug >= 5) std::cout << entry->name() << " readFittedEntriesFromFile " << entry->value() << " " << ALIUtils::getFloat( wl[ii+1] ) << std::endl;
    if( entry->name().substr(0,6) != "angles" ) {
      entry->displaceOriginalOriginal(entry->value() - ALIUtils::getFloat( wl[ii+1] )*theLengthDim);
    } else {
      OpticalObject* opto = entry->OptOCurrent();
      if( affAngles.find(opto) == affAngles.end() ){
        affAngles[opto] = ALIRmDataFromFile();
      }
      std::map<OpticalObject*,ALIRmDataFromFile>::iterator ite = affAngles.find(opto);
      (*ite).second.setAngle( optOentryName.substr(optOentryName.size()-1,1), ALIUtils::getFloat( wl[ii+1] )*theAngleDim );
      if( ALIUtils::debug >= 5) std::cout << " setting angle from file " << ALIUtils::getFloat( wl[ii+1] )*theAngleDim << " " << wl[ii+1] << " " << theAngleDim << std::endl;
    }
    entry->setSigma( ALIUtils::getFloat( wl[ii+2] )*theAngleDim );
    //  ar.lass6.laser.angles_X -159.7524 7.2208261 

  }

  ALIstring coordi("XYZ");
  std::map<OpticalObject*,ALIRmDataFromFile>::const_iterator ite;
  for( ite = affAngles.begin(); ite != affAngles.end(); ite++ ){
    ALIRmDataFromFile dff = (*ite).second;
    OpticalObject* opto = (*ite).first;
    for( size_t ii = 0; ii < 3; ii++ ) {
      int ifound =  dff.dataFilled().find( coordi[ii] );
      if( ALIUtils::debug >= 5) std::cout << ii << " dataFilled " << ifound << std::endl; 
      if( ifound == -1 ) { //angles not read from file are taken as the original value
        ALIdouble entval = opto->getEntryRMangle(coordi.substr(ii,1));
        dff.setAngle( coordi.substr(ii,1), entval );
      } 
    }
    dff.constructRm();
    opto->setGlobalRMOriginalOriginal( dff.rm() );
  }
    
  return true; // to avoid warning
  
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
ALIstring FittedEntriesReader::substitutePointBySlash( const ALIstring& nameWithPoints ) const
{
  ALIstring nameWithSlash = nameWithPoints;

  size_t siz = nameWithPoints.length();

  for( size_t ii = 0; ii < siz; ii++ ){
    if( nameWithSlash[ii] == '.' ) nameWithSlash[ii] = '/'; 
  }
  nameWithSlash = "s/" + nameWithSlash;
  if( ALIUtils::debug >= 5) std::cout << " substitutePointBySlash " << nameWithSlash << " " << std::endl;


  return nameWithSlash;
}