Model Class Reference

#include <Model.h>

Public Member Functions
void	BuildMeasurementsFromOA (OpticalAlignMeasurements &measList)
void	BuildSystemDescriptionFromOA (OpticalAlignments &optAlig)
void	CMSLinkFit (ALIint cmslink)
	*************** FOR CMS LINK SYSTEM (to fit it part by part) More...
OpticalAlignInfo	FindOptAlignInfoByType (const ALIstring &type)
	Model ()
	~Model ()

Static Public Member Functions
static void	addEntryToList (Entry *entry)
	***************** SET DATA MEMBERS More...
static void	addMeasurementToList (Measurement *measadd)
static ALIbool	createCopyComponentList (const ALIstring &optoname)
	**************** FOR COPYING AN OPTO More...
static std::vector< Entry * > &	EntryList ()
static cocoaStatus	getCocoaStatus ()
static ALIbool	getComponentOptOs (const ALIstring &opto_name, std::vector< OpticalObject * > &vcomponents)
	--— Get from theOptOList the list of pointers to component OptOs More...
static ALIbool	getComponentOptOTypes (const ALIstring &opto_type, std::vector< ALIstring > &vcomponents)
	--— Get from theOptODictionary the list of component OptO types More...
static Entry *	getEntryByName (const ALIstring &opto_name, const ALIstring &entry_name)
	--— Search an Entry name in the Entry* list and return a pointer to it More...
static Entry *	getEntryByName (const ALIstring &opto_entry_name)
static FittedEntriesReader *	getFittedEntriesReader ()
static Model &	getInstance ()
	-------— Gets the only instance of this class More...
static Measurement *	getMeasurementByName (const ALIstring &name, ALIbool exists=1)
static std::vector < OpticalAlignInfo >	getOpticalAlignments ()
static OpticalObject *	getOptOByName (const ALIstring &opto_name)
	--— Find an OptO name in theOptOList and return a pointer to it More...
static OpticalObject *	getOptOByType (const ALIstring &type)
	--— Find the first OptO of type 'opto_type' in theOptOList and return a pointer to it More...
static int	getParameterValue (const ALIstring &sstr, ALIdouble &val)
	************ ACCESS INFO FROM STATIC DATA More...
static ALIdouble	getParamFittedSigmaVectorItem (const ALIuint position)
static ALIstring &	MatricesFName ()
	the name of the File for storing the matrices More...
static ALIstring &	MeasFName ()
	the name of the Measurements File More...
static std::vector < Measurement * > &	MeasurementList ()
static struct tm &	MeasurementsTime ()
static OpticalObject *	nextOptOToCopy ()
static std::vector < std::vector< ALIstring > > &	OptODictionary ()
	ACCESS STATIC DATA MEMBERS. More...
static std::vector < OpticalObject * > &	OptOList ()
static std::string	printCocoaStatus (const cocoaStatus cs)
static ALIbool	readMeasurementsFromFile (ALIstring only1Date=ALIstring(""), ALIstring only1Time=ALIstring(""))
static void	readSystemDescription ()
	-------— Read the different sections of the SDF and act accordingly More...
static ALIstring &	ReportFName ()
	the name of the report File More...
static ALIstring &	SDFName ()
	the name of the System Description File More...
static void	setCocoaStatus (const cocoaStatus cs)
static void	setMatricesFName (const ALIstring &name)
static void	setMeasurementsTime (struct tm &tim)
static void	setReportFName (const ALIstring &name)
static void	setSDFName (const ALIstring &name)

Static Public Attributes
static std::vector< ALIdouble >	CMSLinkRangeDetValue
static ALIint	Ncmslinkrange = 0
	*************** FOR RANGE STUDIES More...

Private Member Functions
void	CMSLinkCleanModel ()

Static Private Member Functions
static void	buildMeasurementsLinksToOptOs ()
	Read Measurements (to be implemented for reading from an external file the DATA of the measurements) More...
static void	cleanParamFittedSigmaVector ()
static void	cleanParamFittedValueDisplacementMap ()
static void	CMSLinkDeleteOptOs ()
static void	CMSLinkRecoverParamFittedSigma (ALIint cmslink)
static void	CMSLinkRecoverParamFittedValueDisplacement (ALIint cmslink)
static void	CMSLinkSaveParamFittedSigma (ALIint cmslink)
static void	CMSLinkSaveParamFittedValueDisplacement (ALIint cmslink)
static void	copyMeasurements (const std::vector< ALIstring > &wl)
static void	deleteOptO (const ALIstring &opto_name)
static void	deleteOptO (OpticalObject *opto)
static ALIbool	fillCopyComponentList (const OpticalObject *opto)
static void	recoverParamFittedSigma (const ALIstring &opto_name, const ALIstring &entry_name, const ALIuint position)
static void	reorderOptODictionary (const ALIstring &ssearch, std::vector< std::vector< ALIstring > > &OptODictionary2)
	********** private METHODS More...
static void	saveParamFittedCorrelation (const ALIstring &opto_name1, const ALIstring &entry_name1, const ALIstring &opto_name2, const ALIstring &entry_name2)
static void	saveParamFittedSigma (const ALIstring &opto_name, const ALIstring &entry_name)
static void	SetValueDisplacementsFromReportOut ()

Static Private Attributes
static ALIint	CMSLinkIteration = 0
static cocoaStatus	theCocoaStatus = COCOA_Init
static std::vector< Entry * >	theEntryVector
	std::vector of all Entries More...
static FittedEntriesReader *	theFittedEntriesReader = 0
static Model *	theInstance = 0
	Only instance of Model. More...
static ALIstring	theMatricesFName = "matrices.out"
	the name of the File for storing the matrices More...
static ALIstring	theMeasFName = "Measurements.txt"
	the name of the Measurements File More...
static struct tm	theMeasurementsTime
static std::vector< Measurement * >	theMeasurementVector
	std::vector of all Measurements More...
static std::vector < OpticalAlignInfo >	theOpticalAlignments
static std::vector < std::vector< ALIstring > >	theOptODictionary
	parameters More...
static std::vector < OpticalObject * >	theOptOList
	map of OptO*/type of parent OptO, for navigation down the tree structure More...
static std::vector < OpticalObject * >	theOptOsToCopyList
static std::vector < OpticalObject * > ::const_iterator	theOptOsToCopyListIterator
static std::vector< ALIdouble >	theParamFittedSigmaVector
static std::map< ALIstring, ALIdouble, std::less < ALIstring > >	theParamFittedValueDisplacementMap
static ALIstring	theReportFName = "report.out"
	the name of the report File More...
static ALIstring	theSDFName = "SystemDescription.txt"
	the name of the System Description File More...

Detailed Description

Definition at line 38 of file Model.h.

Constructor & Destructor Documentation

Model::Model

(

)

Definition at line 93 of file Model.cc.

{
  //  theMeasurementsTime = clock();
}

Model::~Model ( )

inline

Definition at line 43 of file Model.h.

{};

Member Function Documentation

static void Model::addEntryToList ( Entry *  entry )

inlinestatic

***************** SET DATA MEMBERS

Definition at line 145 of file Model.h.

References theEntryVector.

Referenced by OpticalObject::fillCoordinateEntry(), OpticalObject::fillExtraEntry(), and OpticalObject::setAnglesNull().

                                             {
     theEntryVector.push_back( entry );
     //-     std::cout << entry << entry->OptOCurrent()->name() << "ADDENTRY " << entry->name() << " " << EntryList().size() << std::endl;
  }

static void Model::addMeasurementToList ( Measurement *  measadd )

inlinestatic

Definition at line 150 of file Model.h.

References theMeasurementVector.

Referenced by Measurement::postConstruct().

                                                           {
      theMeasurementVector.push_back( measadd);  
      //   std::cout << "ADD MEASUREMENT" << theMeasurementVector.size() << std::endl ;  
  }

void Model::BuildMeasurementsFromOA

(

OpticalAlignMeasurements &

measList

)

Definition at line 1626 of file Model.cc.

References dtNoiseDBValidation_cfg::cerr, Measurement::constructFromOA(), gather_cfg::cout, ALIUtils::debug, cppFunctionSkipper::exception, and OpticalAlignMeasurements::oaMeasurements_.

Referenced by CocoaAnalyzer::RunCocoa().

{
  std::vector<OpticalAlignMeasurementInfo>::iterator mite;

  if( ALIUtils::debug >= 5 ) std::cout << " BuildMeasurementsFromOA " << std::endl;
  std::vector<OpticalAlignMeasurementInfo> measInfos = measList.oaMeasurements_;
  for( mite = measInfos.begin(); mite != measInfos.end(); mite++ ) {
    std::string measType = (*mite).type_;
    std::string measName = (*mite).name_;
  if( ALIUtils::debug >= 4 ) std::cout << " BuildMeasurementsFromOA measType " << measType << " measName " << measName << std::endl;
    //---------- Create Measurement with appropiate dimension
    Measurement* meastemp = 0;
    if ( measType == ALIstring("SENSOR2D") ) {
      meastemp = new MeasurementSensor2D( 2, measType, measName );  
    } else if ( measType == ALIstring("DISTANCEMETER3DIM") ) {
      meastemp = new MeasurementDistancemeter3dim( 1, measType, measName );
    } else if ( measType == ALIstring("DISTANCEMETER") ||
        measType == ALIstring("DISTANCEMETER1DIM") ) {
      meastemp = new MeasurementDistancemeter( 1, measType, measName );
    } else if ( measType == ALIstring("TILTMETER") ) {
      meastemp = new MeasurementTiltmeter( 1, measType, measName );
    } else if ( measType == ALIstring("COPS") ) {
      meastemp = new MeasurementCOPS( 4, measType, measName );
    } else if ( measType == ALIstring("DIFFCENTRE") ) {
      //t       meastemp = new MeasurementDiffCentre( 1, measType, measName );
    } else if ( measType == ALIstring("DIFFANGLE") ) {
      //t        meastemp = new MeasurementDiffAngle( 2, measType, measName );
    } else if ( measType == ALIstring("DIFFENTRY") ) {
      meastemp = new MeasurementDiffEntry( 1, measType, measName );
    } else {
      std::cerr << " !!! Model::BuildMeasurementsFromOA : measType not found " << measType << std::endl;
    std::exception();
    } 
    meastemp->constructFromOA( *mite );

  }
}

void Model::buildMeasurementsLinksToOptOs ( )

staticprivate

Read Measurements (to be implemented for reading from an external file the DATA of the measurements)

Build for each measuremnt its link to the OptO that take part in it

Definition at line 546 of file Model.cc.

References begin.

Referenced by readSystemDescription().

                                         {

//---------- Loop Measurements
  std::vector< Measurement* >::const_iterator vmcite;
  for ( vmcite = MeasurementList().begin(); vmcite != MeasurementList().end(); vmcite++) {
//---------- Transform for each Measurement the Measured OptO names to Measured OptO pointers
    //     (*vmcite)->buildOptOList();

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

  }
 
}

void Model::BuildSystemDescriptionFromOA

(

OpticalAlignments &

optAlig

)

Definition at line 1580 of file Model.cc.

References OpticalObject::constructFromOptAligInfo(), OpticalAlignInfo::name_, and OpticalAlignments::opticalAlignments().

Referenced by CocoaAnalyzer::RunCocoa().

{
  theOpticalAlignments = optAlig.opticalAlignments();

  OpticalAlignInfo oai_system = FindOptAlignInfoByType( "system" );

  OpticalObject* OptOsystem = new OpticalObject( 0, "system", oai_system.name_, 0 );
 
  OptOsystem->constructFromOptAligInfo( oai_system );

  //-              Model::_OptOtree.insert( std::multimap< ALIstring, OpticalObject*, std::less<ALIstring> >::value_type(OptOsystem->type(), OptOsystem) );
  //              theOptOlist[OptOsystem->name()] = OptOsystem; 
  theOptOList.push_back( OptOsystem ); 
}

static void Model::cleanParamFittedSigmaVector ( )

inlinestaticprivate

Definition at line 271 of file Model.h.

References theParamFittedSigmaVector.

                                            {
    ALIuint pfsv_size = theParamFittedSigmaVector.size();    
    for( ALIuint ii = 0; ii < pfsv_size; ii++) {
      theParamFittedSigmaVector.pop_back();
    }
  }

static void Model::cleanParamFittedValueDisplacementMap ( )

inlinestaticprivate

Definition at line 278 of file Model.h.

References theParamFittedValueDisplacementMap.

                                                     {
    theParamFittedValueDisplacementMap.erase( theParamFittedValueDisplacementMap.begin(), theParamFittedValueDisplacementMap.end() );   
  }

void Model::CMSLinkCleanModel ( )

private

Definition at line 887 of file Model.cc.

{
  deleteOptO("s");
  //---------- Clean OptOdicitionary (in case this is not first reading)
  ALIuint odsize = theOptODictionary.size(); 
  for( ALIuint ii = 0; ii < odsize; ii++) {
      theOptODictionary.pop_back(); 
  }

}

void Model::CMSLinkDeleteOptOs ( )

staticprivate

Definition at line 902 of file Model.cc.

References CMSLinkIteration, gather_cfg::cout, ALIUtils::debug, GlobalOptionMgr::getGlobalOptionValue(), and GlobalOptionMgr::getInstance().

{
  ALIint cmslink_iter = Model::CMSLinkIteration; 
  ALIdouble cmslink_method;

  GlobalOptionMgr* gomgr = GlobalOptionMgr::getInstance();
  assert(gomgr->getGlobalOptionValue("cms_link_method", cmslink_method));
  ALIdouble cmslink_halfplanes;
  assert(gomgr->getGlobalOptionValue("cms_link_halfplanes", cmslink_halfplanes));
  if(ALIUtils::debug >= 2) std::cout << "CMSLinkDeleteOptOs: cms_link_halfplanes " << cmslink_halfplanes <<  cmslink_iter  << std::endl;

  if( cmslink_iter == 1){
    //---------- First fit: delete everything but laser1 and det_tkU
    //    deleteOptO("s/laserboxL/laser2");
    //-    std::cout << "delete mabs" << std::endl;
    deleteOptO("s/mabsL");
    //------- Method 1: detectors at tracker down
    if (cmslink_method == 1){
      deleteOptO("s/tracker/det_trkDL");
      deleteOptO("s/tracker/det_trkDR");
    }

    if( cmslink_halfplanes == 2) { 
      //      deleteOptO("s/laserboxR/laser2");
      deleteOptO("s/mabsR");
    }

  } else if( cmslink_iter == 2) {
    //---------- Second fit (method 1): delete everything but laser1 and det3
    //    deleteOptO("s/laserboxL/laser2");
    deleteOptO("s/mabsL");
    deleteOptO("s/tracker/CST/wheel_trkL/peri/mirror"); //??
    deleteOptO("s/tracker/CST/wheel_trkL/det_trkU");
    //------- Method 1: detectors on CST, Method 2: detectors on tracker
      //not necessary    deleteOptO("s/tracker/CST/det6");

    if( cmslink_halfplanes <= 1) {  
      deleteOptO("s/tracker/CST/wheel_trkR");
    } else if( cmslink_halfplanes == 2) { 
      //      deleteOptO("s/laserboxR/laser2");
      deleteOptO("s/mabsR");
      deleteOptO("s/tracker/CST/wheel_trkR/peri/mirror"); //??
      deleteOptO("s/tracker/CST/wheel_trkR/det_trkU");
    }

  } else if( cmslink_iter == 3) {
    //---------- Third fit: delete everything but laser2 and mabs
    //    deleteOptO("s/laserboxL/laser1");
    deleteOptO("s/tracker");

    if( cmslink_halfplanes == 2) {  
      //      deleteOptO("s/laserboxR/laser1");
    }
  //---------- Do nothing
  } else {

  }

}

void Model::CMSLinkFit

(

ALIint

cmslink

)

*************** FOR CMS LINK SYSTEM (to fit it part by part)

Definition at line 829 of file Model.cc.

{
  /*
  //---------- Get number of fits 
  ALIint cmslinkNoFits = 0;
  if( cmslink == 1) {
    cmslinkNoFits = 3;
  } else if( cmslink == 2) {
    cmslinkNoFits = 2;
  }
  if(ALIUtils::debug >= 4) std::cout << " Start CMS link fit with method " << cmslink << " (number of fits = " << cmslinkNoFits << ")" << std::endl;

  //----------- Loop number of cms_link fits 
  for(ALIint ilink = ALIint(Model::GlobalOptions()["cms_link"]); ilink <= cmslinkNoFits; ilink++) { 

  //----- Iteration 2 of method 2 resembles iteration 3 of method 1
     ALIdouble cmslink_method;
     assert(Model::getGlobalOptionValue("cms_link_method", cmslink_method)); 
     if( cmslink_method == 2 && ilink == 2) ilink = 3;

//---- Set variable CMSLinkIteration, Checked in Fit.C and other     //- std::cout << "ilink" << ilink << std::endl;
    Model::setGlobalOption("cms_link", ilink); 
    Model::CMSLinkIteration = ilink;

    if(ilink > 1)Model::readSystemDescription(); //already read once to fill value Model::GlobalOptions()["cms_link"]

    //---------- Delete the OptO not fitted in this iteration
    //    Model::CMSLinkDeleteOptOs(); 
    // cannot be here because you may recover a parameter saved in previous iteration that now it is deleted (!!MODIFY THIS)  

    //---------- Recover parameters fitted in previous iteration
    Model::CMSLinkRecoverParamFittedSigma( ilink );

    Model::CMSLinkRecoverParamFittedValueDisplacement( ilink );

    //---------- Delete the OptO not fitted in this iteration
    Model::CMSLinkDeleteOptOs();
  
    //---------- Start fit
    Fit::startFit();

    //---------- Save parameters fitted in this iteration (to be used in next one)
    Model::CMSLinkSaveParamFittedSigma( ilink );

    Model::CMSLinkSaveParamFittedValueDisplacement( ilink );

    //---------- Delete whole system to start anew in next iteration
    Model::CMSLinkCleanModel();

  }    
  */
}

void Model::CMSLinkRecoverParamFittedSigma ( ALIint  cmslink )

staticprivate

Definition at line 1058 of file Model.cc.

{
  /*
  ALIdouble cms_link_halfplanes = (GlobalOptions()["cms_link_halfplanes"]);

  if( cms_link == 2) {
  //---------- Second fit: recover laserbox angles
    recoverParamFittedSigma("s/laserboxL","Angles X",0);
    recoverParamFittedSigma("s/laserboxL","Angles Y",1);

    if( cms_link_halfplanes == 2) {  
      recoverParamFittedSigma("s/laserboxR","Angles X",0);
      recoverParamFittedSigma("s/laserboxR","Angles Y",1);
    }

  } else if( cms_link == 3) {
  //---------- Third fit: recover laserbox angles and position and rotate angles to mabs
    recoverParamFittedSigma("s/laserboxL","Centre X",0);
    recoverParamFittedSigma("s/laserboxL","Centre Y",1);
    recoverParamFittedSigma("s/laserboxL","Centre Z",2);
    recoverParamFittedSigma("s/laserboxL","Angles X",3);

    //----- Angle around Y is converted to angle around Z when turning 90 deg
    Entry* slaZ = getEntryByName("s/laserboxL","Angles Z");
    //--- prec_level_laser
    Entry* smaZ = getEntryByName("s/mabsL","Angles Z");
    slaZ->setQuality(0);
    slaZ->setValue( smaZ->value() );
    //    smaZ->setQuality(0); //!!???!!?

    Entry* slaY = getEntryByName("s/laserboxL","Angles Y");
    slaY->setQuality(0); 

    if( cms_link_halfplanes == 2) {   
      recoverParamFittedSigma("s/laserboxR","Centre X",0);
      recoverParamFittedSigma("s/laserboxR","Centre Y",1);
      recoverParamFittedSigma("s/laserboxR","Centre Z",2);
      recoverParamFittedSigma("s/laserboxR","Angles X",3);

      //----- Angle around Y is converted to angle around Z when turning 90 deg
      Entry* slaZ = getEntryByName("s/laserboxR","Angles Z");
      //--- prec_level_laser
      Entry* smaZ = getEntryByName("s/mabsR","Angles Z");
      slaZ->setQuality(0);
      slaZ->setValue( smaZ->value() );
      //    smaZ->setQuality(0); //!!???!!?

      Entry* slaY = getEntryByName("s/laserboxR","Angles Y");
      slaY->setQuality(0); 
    }
  } else {
  //---------- Do nothing

  }
  */
}

void Model::CMSLinkRecoverParamFittedValueDisplacement ( ALIint  cmslink )

staticprivate

Definition at line 1119 of file Model.cc.

{
  /*
  //---------- Second fit: recover laserbox angles
  if( cms_link == 2 || cms_link == 3 ) {

    std::map<ALIstring, ALIdouble, std::less<ALIstring> >::const_iterator vsdmite;
    if ( ALIUtils::debug >= 99) std::cout << "theParamFittedValueDisplacementMap.size " << theParamFittedValueDisplacementMap.size() << std::endl;
    for( vsdmite =  theParamFittedValueDisplacementMap.begin(); vsdmite !=  theParamFittedValueDisplacementMap.end(); vsdmite++) {
      std::cout << "reoverValueDisp" <<  (*vsdmite).first << "  " << (*vsdmite).second << std::endl;
      Entry* this_entry = getEntryByName( (*vsdmite).first);
      this_entry->displaceOriginal( (*vsdmite).second ); 
      this_entry->OptOCurrent()->resetGlobalCoordinates();
      this_entry->setValueDisplacementByFitting(  (*vsdmite).second ); 

    }

  //---------- Do nothing
  } else {

  }
  */
}

void Model::CMSLinkSaveParamFittedSigma ( ALIint  cmslink )

staticprivate

Definition at line 966 of file Model.cc.

{
  /*
  ALIdouble cms_link_halfplanes = (GlobalOptions()["cms_link_halfplanes"]);

  cleanParamFittedSigmaVector(); 

  //---------- First fit: Save laserbox angles fitted at tracker 
  if( cms_link == 1) {
    //?    if (GlobalOptions()["cms_link_method"] < 10){
    saveParamFittedSigma("s/laserboxL","Angles X");   
    saveParamFittedSigma("s/laserboxL","Angles Y");

    saveParamFittedCorrelation("s/laserboxL","Angles X",
     "s/tracker/CST","Centre Y");
    saveParamFittedCorrelation("s/laserboxL","Angles Y",
     "s/tracker/CST","Centre X");
    if( cms_link_halfplanes == 2) {  
      saveParamFittedSigma("s/laserboxR","Angles X");   
      saveParamFittedSigma("s/laserboxR","Angles Y");

      saveParamFittedCorrelation("s/laserboxR","Angles X",
       "s/tracker/CST","Centre Y");
      saveParamFittedCorrelation("s/laserboxR","Angles Y",
       "s/tracker/CST","Centre X");
    }

  } else if( cms_link == 2) {
  //---------- Second fit: Save laserbox angles and position
    saveParamFittedSigma("s/laserboxL","Centre X");
    saveParamFittedSigma("s/laserboxL","Centre Y");
    // Make quality unk to cal
    saveParamFittedSigma("s/laserboxL","Centre Z");
    saveParamFittedSigma("s/laserboxL","Angles X");
    saveParamFittedSigma("s/laserboxL","Angles Y");
   
    saveParamFittedCorrelation("s/laserboxL","Centre X",
        "s/laserboxL","Angles Y");
    saveParamFittedCorrelation("s/laserboxL","Centre Y",
        "s/laserboxL","Angles X");

   if( cms_link_halfplanes == 2) {  
    saveParamFittedSigma("s/laserboxR","Centre X");
    saveParamFittedSigma("s/laserboxR","Centre Y");
    // Make quality unk to cal
    saveParamFittedSigma("s/laserboxR","Centre Z");
    saveParamFittedSigma("s/laserboxR","Angles X");
    saveParamFittedSigma("s/laserboxR","Angles Y");
   
    saveParamFittedCorrelation("s/laserboxR","Centre X",
        "s/laserboxR","Angles Y");
    saveParamFittedCorrelation("s/laserboxR","Centre Y",
        "s/laserboxR","Angles X");
   }  
  } else {
 //---------- Do nothing

  }
  */
}

void Model::CMSLinkSaveParamFittedValueDisplacement ( ALIint  cmslink )

staticprivate

Definition at line 1031 of file Model.cc.

{
  /*
  cleanParamFittedValueDisplacementMap(); 

  //---------- 
  if( cms_link == 1 || cms_link == 2 ) {
    std::vector<Entry*>::const_iterator vecite;
    for( vecite = EntryList().begin(); vecite != EntryList().end(); vecite++) {
      if( (*vecite)->valueDisplacementByFitting() != 0 ) {
        ALIstring names = (*vecite)->OptOCurrent()->name() + "/" + (*vecite)->name();
        std::cout << "saeParamFittedValueDisplacementMap" << names << (*vecite)->valueDisplacementByFitting() << std::endl;
        theParamFittedValueDisplacementMap[ names ] = (*vecite)->valueDisplacementByFitting();
      }
    }

  //---------- Do nothing
  } else {

  }
  */
}

void Model::copyMeasurements ( const std::vector< ALIstring > &  wl )

staticprivate

Definition at line 1458 of file Model.cc.

References ALIUtils::changeName(), Measurement::copyMeas(), gather_cfg::cout, spr::find(), Measurement::name(), and Measurement::type().

{
 
  //----- Check format, e.g.  @copy_measurements _1/_2/*_1 


  //----- get list of Measurement's that satisfy the query in their name
  //t  std::string querystr = wl[1].substr( wl[1].find("/")+1, wl[1].length() );
  std::string subsstr1 = wl[1].substr( 0, wl[1].find("/") );
  std::string subsstr2 = wl[1].substr( wl[1].find("/")+1, wl[1].rfind("/") - wl[1].find("/")-1 );
  std::string querystr = wl[1].substr( wl[1].rfind("/")+1, wl[1].length() );

  std::cout << " Model::copyMeasurements " << " subsstr1 " << subsstr1 << " subsstr2 " << subsstr2 << " querystr " << querystr << std::endl;

  std::vector< Measurement* > measToCopy;
  std::vector< Measurement* >::iterator mite;
  for( mite = theMeasurementVector.begin(); mite != theMeasurementVector.end(); mite++) {
    Measurement* meas = (*mite);
    //improve this
    if( meas->name().find( querystr ) != std::string::npos ) {
      measToCopy.push_back( meas );
    }
  }

  //---- Build new measurements
  Measurement* meastemp = 0;
  for( mite = measToCopy.begin(); mite != measToCopy.end(); mite++) {
    Measurement* meas = (*mite);
    std::vector<ALIstring> wlt;
    wlt.push_back( meas->type() );

    //---- create new name 
    std::string newName = ALIUtils::changeName( meas->name(), subsstr1, subsstr2 );
    std::cout << " newName " << newName << std::endl;
    wlt.push_back( newName );

    ALIstring measType = wlt[0];
    ALIstring measName;
    if( wlt.size() == 2 ) {
      measName = wlt[1];
    } else {
      measName = "";
    }
    if ( meas->type() == ALIstring("SENSOR2D") ) {
      meastemp = new MeasurementSensor2D( 2, measType, measName );      
      //      } else if ( meas->type() == ALIstring("DISTANCEMETER3DIM") ) {
      //        meastemp = new MeasurementDistancemeter3dim( 1, measType, measName );
    } else if ( meas->type() == ALIstring("DISTANCEMETER") ||
        meas->type() == ALIstring("DISTANCEMETER1DIM") ) {
      meastemp = new MeasurementDistancemeter( 1, measType, measName );
    } else if ( meas->type() == ALIstring("TILTMETER") ) {
      meastemp = new MeasurementTiltmeter( 1, measType, measName );
    } else if ( meas->type() == ALIstring("DIFFCENTRE") ) {
 //     meastemp = new MeasurementDiffCentre( 1, measType, measName );
    } else if ( meas->type() == ALIstring("DIFFANGLE") ) {
//      meastemp = new MeasurementDiffAngle( 1, measType, measName );
    } else if ( meas->type() == ALIstring("DIFFENTRY") ) {
      meastemp = new MeasurementDiffEntry( 1, measType, measName );
    } else if ( meas->type() == ALIstring("COPS") ) {
      meastemp = new MeasurementCOPS( 4, measType, measName );
    }

    //later        meastemp->copyConversionFactor( wordlist );
    meastemp->copyMeas(meas, subsstr1, subsstr2);
    
    break;
  }

}

ALIbool Model::createCopyComponentList ( const ALIstring &  optoname )

static

**************** FOR COPYING AN OPTO

Definition at line 769 of file Model.cc.

References gather_cfg::cout, and ALIUtils::debug.

Referenced by OpticalObject::createComponentOptOs().

{
  //---------- Find an OptO with the same type (YOU HAVE TO BE SURE THAT ALL EXISTING OPTOs OF THIS TYPE HAVE SAME COMPONENTS, IF NOT COPYING MAY GIVE YOU UNPREDICTABLE RESULTS)
  if(ALIUtils::debug >= 3) std::cout << "createCopyComponentList " << typ << std::endl; 
  OpticalObject* start_opto = getOptOByType( typ );
  
  //---------- clean list of OptOs to copy
  theOptOsToCopyList.erase( theOptOsToCopyList.begin(), theOptOsToCopyList.end() );
  
  //---------- Fill list of OptOs to copy
  fillCopyComponentList( start_opto );
  //- if(ALIUtils::debug >= 9) std::cout << "createCopyComponentList " << typ << theOptOsToCopyList.size() << std::endl; 
  
  theOptOsToCopyListIterator = theOptOsToCopyList.begin();
  return 1;
}

void Model::deleteOptO ( const ALIstring &  opto_name )

staticprivate

Definition at line 1147 of file Model.cc.

{
  OpticalObject* opto = getOptOByName( opto_name );
  deleteOptO( opto );
}

void Model::deleteOptO ( OpticalObject *  opto )

staticprivate

Definition at line 1157 of file Model.cc.

References OpticalObject::CoordinateEntryList(), gather_cfg::cout, ALIUtils::debug, OpticalObject::ExtraEntryList(), spr::find(), and OpticalObject::name().

{
  if(ALIUtils::debug >= 5 ) std::cout << "DELETING OptO" << opto->name() << std::endl;
  
  //---------- Delete Coordinate Entries of this OptO
  std::vector<Entry*>::const_iterator vecite;
  std::vector<Entry*>::iterator veite2;
  if ( ALIUtils::debug >= 9) std::cout << "SIZE" << theEntryVector.size() <<std::endl;
  for( vecite = opto->CoordinateEntryList().begin(); vecite != opto->CoordinateEntryList().end(); vecite++) {
    //    ALIuint pos = FindItemInVector( (*veite), opto->CoordinateEntryList() );
    veite2 = find( theEntryVector.begin(), theEntryVector.end(), (*vecite) );
    //-  if ( ALIUtils::debug >= 9) std::cout << (*veite2) << "DELETE ENTRY " << (*vecite) <<(*veite2)->OptOCurrent()->name() << (*veite2)->name() << std::endl;
    delete ( (*veite2) );
    theEntryVector.erase( veite2 );
  }

  for( vecite = opto->ExtraEntryList().begin(); vecite != opto->ExtraEntryList().end(); vecite++) {
    //    ALIuint pos = FindItemInVector( (*veite), opto->CoordinateEntryList() );
    veite2 = find( theEntryVector.begin(), theEntryVector.end(), (*vecite) );
    //-    if(ALIUtils::debug >= 9) std::cout << (*veite2) << "DELETE ENTRY " << (*veite2)->OptOCurrent()->name() << (*veite2)->name() << std::endl;
    delete ( (*veite2) );
    theEntryVector.erase( veite2 );
  }

 
  for( vecite = theEntryVector.begin(); vecite != theEntryVector.end(); vecite++ ) {
    //     std::cout << (*vecite) << "ENTReY " << (*vecite)->OptOCurrent()->name() << (*vecite)->name() << std::endl;
  }

  //---------- Delete all Measurement in which opto takes part
  std::vector<Measurement*> MeasToBeDeleted;
  std::vector<Measurement*>::const_iterator vmite;
  std::vector<OpticalObject*>::const_iterator vocite;
  for( vmite = MeasurementList().begin(); vmite != MeasurementList().end(); vmite++ ) {
      if(ALIUtils::debug >= 5) std::cout << "Deleting Measurement" << (*vmite)->name() << std::endl;
    //----- If any of the OptO Measured is opto, delete this Measurement
    for( vocite = (*vmite)->OptOList().begin(); vocite != (*vmite)->OptOList().end(); vocite++) {
      if( (*vocite) == opto ) {
    //-      std::cout << "MEASTBD" << (*vmite) << std::endl;
        MeasToBeDeleted.push_back( *vmite );
    //?       delete (*vmite);
        break;
      }
    } 
  }

 //---------- Delete Measurements from list
  std::vector<Measurement*>::const_iterator vmcite;
  std::vector<Measurement*>::iterator vmite2;
  if ( ALIUtils::debug >= 9) std::cout << "SIZEMEAS" << MeasToBeDeleted.size() <<std::endl;
  for( vmcite = MeasToBeDeleted.begin(); vmcite != MeasToBeDeleted.end(); vmcite++) {
    vmite2 = find( theMeasurementVector.begin(), theMeasurementVector.end(), (*vmcite) );
    //    std::cout << (*vmite2) << "DELETE MSEASU " << (*vmcite) << (*vmite2)->name()[0] << std::endl;
    delete ( (*vmite2) );
    theMeasurementVector.erase( vmite2 );
  }


  //---------- Delete components
  //  std::vector<OpticalObject*>::iterator voite;
  std::vector<OpticalObject*> vopto;
  //  ALIbool opto_found = getComponentOptOs( opto->name(), vopto );
  for( vocite = vopto.begin(); vocite != vopto.end(); vocite++ ) {
    deleteOptO( *vocite );
  }
     
  //---------- Delete this OptO
  //---------- Delete OptO (only from list, to delete it really first delete components)
  /*  map< ALIstring, OpticalObject*, std::less<ALIstring> >::iterator msoite =  
find( theOptOList.begin(), theOptOList.end(), 
map< ALIstring, OpticalObject*, std::less<ALIstring> >::value_type( opto->name(), opto) );*/
  std::vector< OpticalObject* >::iterator dvoite =  
   find( theOptOList.begin(), theOptOList.end(), 
         std::vector< OpticalObject* >::value_type( opto) );
  //-  std::cout << (*dvoite) << "DELETE OPTO " << opto <<"WW" << (*dvoite)->name() << std::endl;
  theOptOList.erase( dvoite );
  delete opto;

}

static std::vector< Entry* >& Model::EntryList ( )

inlinestatic

Definition at line 75 of file Model.h.

References theEntryVector.

Referenced by Fit::addDaMatrixToEntries(), Fit::CheckIfFitPossible(), Fit::CreateMatrices(), Fit::dumpEntryCorrelations(), NtupleManager::FillChi2(), FittedEntriesSet::FillCorrelations(), FittedEntriesSet::FillEntries(), NtupleManager::FillFitParameters(), Fit::FillMatricesWithCalibratedParameters(), Fit::fitNextEvent(), Fit::multiplyMatrices(), Fit::PrintChi2(), Fit::setFittableEntries(), SetValueDisplacementsFromReportOut(), and Fit::substractLastDisplacementToEntries().

                                          {
    return theEntryVector;
  }

ALIbool Model::fillCopyComponentList ( const OpticalObject *  opto )

staticprivate

Definition at line 790 of file Model.cc.

References gather_cfg::cout, ALIUtils::debug, and OpticalObject::name().

{
  if(ALIUtils::debug >= 3) std::cout << "entering fillCopyComponentList(): OptO" << opto->name() << std::endl;
  //---------- Get list of components of 'opto'
  std::vector<OpticalObject*> vopto; 
  ALIbool opto_found = getComponentOptOs( opto->name(), vopto);
  if( !opto_found ) {
    if(ALIUtils::debug >= 5) std::cout << "fillCopyComponentList: NO COMPONENTS TO COPY IN THIS OptO" << opto->name() << std::endl ; 
  }
  
  //---------- Loop list of components of 'opto'
  std::vector<OpticalObject*>::const_iterator vocite;
  for( vocite = vopto.begin(); vocite != vopto.end(); vocite++ ){
    theOptOsToCopyList.push_back( *vocite );
    if(ALIUtils::debug >= 5) std::cout << "fillCopyOptOList " << (*vocite)->type() << " " << (*vocite)->name() << std::endl; 
    //---------- Add components of this component
    fillCopyComponentList( *vocite );
  }
  return opto_found;
}

OpticalAlignInfo Model::FindOptAlignInfoByType

(

const ALIstring &

type

)

Definition at line 1597 of file Model.cc.

References dtNoiseDBValidation_cfg::cerr, cppFunctionSkipper::exception, and OpticalAlignInfo::name_.

{
  OpticalAlignInfo oai;

  ALIbool bFound = 0;
  std::vector<OpticalAlignInfo>::iterator ite;
  for( ite = theOpticalAlignments.begin(); ite != theOpticalAlignments.end(); ite++ ){
    //    std::cout << " Model::FindOptAlignInfoByType " <<  (*ite).type_ << " =? " << type << std::endl;
    if( (*ite).type_ == type ) {
      if( !bFound ){ 
    oai = *ite;
    bFound = 1;
      } else {
    std::cerr << "!! WARNING: Model::FindOptAlignInfoByType more than one objects of type " << type << std::endl;
    std::cerr << " returning object " << oai.name_ << std::endl
          << " skipping object " << (*ite).name_ << std::endl;
      }
    }
  }
  if( !bFound ) { 
    std::cerr << "!! ERROR: Model::FindOptAlignInfoByType object not found, of type " << type << std::endl;
    std::exception();
  }

  return oai;
}

static cocoaStatus Model::getCocoaStatus ( )

inlinestatic

Definition at line 48 of file Model.h.

References theCocoaStatus.

Referenced by Fit::fitParameters().

{ return theCocoaStatus;}

ALIbool Model::getComponentOptOs ( const ALIstring &  opto_name, std::vector< OpticalObject * > &  vcomponents  )

static

--— Get from theOptOList the list of pointers to component OptOs

Definition at line 723 of file Model.cc.

References begin, gather_cfg::cout, ALIUtils::debug, and OpticalObject::name().

Referenced by OpticalObject::displaceCentreGlob(), OpticalObject::displaceCentreGlobOriginal(), OpticalObject::displaceCentreGlobOriginalOriginal(), OpticalObject::displaceRmGlobAroundGlobal(), OpticalObject::displaceRmGlobAroundLocal(), OpticalObject::displaceRmGlobOriginal(), OpticalObject::displaceRmGlobOriginalOriginal(), OpticalObject::propagateGlobalRMOriginalOriginalChangeToChildren(), OpticalObject::resetGlobalCoordinates(), OpticalObject::resetOriginalOriginalCoordinates(), OpticalObject::setGlobalCoordinatesOfComponents(), and OpticalObject::setGlobalRMOriginalOriginal().

{
  //---------- clean std::vector in which you are going to store opto pointers
  std::vector<OpticalObject*>::iterator voite;
  for (voite = vcomponents.begin(); voite != vcomponents.end(); voite++) {
    vcomponents.pop_back();
  }

  //---------- Get OptO corresponding to name 'opto_name'
  OpticalObject* opto = getOptOByName( opto_name );
 
  if ( ALIUtils::debug >= 99) std::cout << opto_name << "getComponentOptOs: opto " << opto << opto->name() << std::endl;
  std::vector<OpticalObject*>::const_iterator vocite;

  if ( ALIUtils::debug >= 99) std::cout << "optolist size " << OptOList().size() << std::endl;
  ALIbool opto_found = 0; 
  for (vocite = OptOList().begin(); vocite != OptOList().end(); vocite++) {
    if( (*vocite)->parent() != 0 ) {
      //        std::cout << "looping OptOlist" << (*vocite)->name() << " parent " <<(*vocite)->parent()->name() << std::endl; 
      if( (*vocite)->parent()->name() == opto_name ) {
        opto_found = 1;
        vcomponents.push_back( (*vocite) );    
      }
    }
  }
    
  /*  std::pair<multimap< ALIstring, OpticalObject*, std::less<ALIstring> >::iterator,
       std::multimap< ALIstring, OpticalObject*, std::less<ALIstring> >::iterator> 
  pmmao =  _OptOtree.equal_range(opto_name);

  if( pmmao.first != _OptOtree.end()) {
    std::multimap< ALIstring, OpticalObject*, std::less<ALIstring> >::const_iterator socite;
    for (socite = pmmao.first; socite != (pmmao.second); socite++) {
         vcomponents.push_back( (*socite).second );
    }
  }
  */
  return opto_found;

}

ALIbool Model::getComponentOptOTypes ( const ALIstring &  opto_type, std::vector< ALIstring > &  vcomponents  )

static

--— Get from theOptODictionary the list of component OptO types

Definition at line 690 of file Model.cc.

References begin.

Referenced by OpticalObject::createComponentOptOs().

{
  //---------- clean std::vector in which you are going to store opto types
  std::vector<ALIstring>::iterator vsite;
  for (vsite = vcomponents.begin(); vsite != vcomponents.end(); vsite++) {
    vcomponents.pop_back();
  }
 
  //---------- Looks the theOptODictionary item that has 'opto_type' as the first ALIstring,
  ALIint ALIstring_found = 0;
  std::vector< std::vector<ALIstring> >::iterator vvsite; 
  for( vvsite = OptODictionary().begin(); vvsite != OptODictionary().end(); ++vvsite) {
    if( *( (*vvsite).begin() ) == opto_type ) {
      ALIstring_found = 1;
      //tt  copies this item to 'vcomponents', substracting the first ALIstring, that is the opto_type itself,
      vcomponents = *vvsite;
      vcomponents.erase( vcomponents.begin() );
      break;
    }
  }
  
  if ( ALIstring_found ) {
    return 1;
  } else {
    return 0;
  } 
} 

Entry * Model::getEntryByName ( const ALIstring &  opto_name, const ALIstring &  entry_name  )

static

--— Search an Entry name in the Entry* list and return a pointer to it

Definition at line 633 of file Model.cc.

References begin, dtNoiseDBValidation_cfg::cerr, gather_cfg::cout, ALIUtils::debug, and cmsRelvalreport::exit.

Referenced by MeasurementDiffEntry::calculateSimulatedValue(), FittedEntry::FittedEntry(), getEntryByName(), FittedEntriesReader::readFittedEntriesFromFile(), Fit::setCorrelationFromParamFitted(), and FittedEntriesSet::SetOptOEntries().

{
  //---------- Look for Entry name type in EntryList
  std::vector<Entry*>::const_iterator vecite;
  for( vecite = EntryList().begin(); vecite != EntryList().end(); vecite++) {
    if( ALIUtils::debug >= 4 ) std::cout <<  "getEntryByName: " <<(*vecite)->OptOCurrent()->name() 
     << " E " << (*vecite)->name() << " Searching: " << opto_name << " E " << entry_name << std::endl; 
    //-    std::cout << " optoName " << (*vecite)->OptOCurrent()->name()<< " " << (*vecite)->name() << std::endl;
    if( (*vecite)->OptOCurrent()->name() == opto_name &&
    (*vecite)->name() == entry_name ) {
      return *vecite;
    }
  }
  //---------- Entry not found!
  std::cerr << "!!!EXITING at getEntryByName: Entry name not found:" << opto_name << "  " << entry_name << std::endl; 
  exit(1);
}

static Entry* Model::getEntryByName ( const ALIstring &  opto_entry_name )

inlinestatic

--— Search an Entry from the full entry path (first substract the name of the OptO and then look in the Entry* list)

Definition at line 122 of file Model.h.

References getEntryByName().

                                                                  {
    ALIint slash_pos = opto_entry_name.rfind('/');
    ALIint length = opto_entry_name.length();
    ALIstring opto_name = opto_entry_name.substr(0, slash_pos);
    ALIstring entry_name = opto_entry_name.substr(slash_pos+1, length);
    Entry* entry = getEntryByName( opto_name, entry_name);
    return entry;
  }

static FittedEntriesReader* Model::getFittedEntriesReader ( )

inlinestatic

Definition at line 267 of file Model.h.

References theFittedEntriesReader.

Referenced by Fit::fitNextEvent().

                                                      {
    return theFittedEntriesReader; }

Model & Model::getInstance ( )

static

-------— Gets the only instance of this class

Definition at line 81 of file Model.cc.

Referenced by CocoaAnalyzer::RunCocoa().

{
  if(!theInstance) {
    theInstance = new Model;
  }
  return *theInstance;
}

Measurement * Model::getMeasurementByName ( const ALIstring &  name, ALIbool  exists = 1  )

static

Definition at line 653 of file Model.cc.

References dtNoiseDBValidation_cfg::cerr.

{
  //---------- Look for Optical Object name in OptOList
  std::vector< Measurement* >::const_iterator vmcite;
  for( vmcite = theMeasurementVector.begin(); vmcite != theMeasurementVector.end(); vmcite++ ) {
    if( (*vmcite)->name() == meas_name ) break;
  }
  
  if ( vmcite != theMeasurementVector.end() ) {
    //---------- If opto_name found, return pointer to it
    return (*vmcite);
  }else {
    if( exists ) {
      //---------- If opto_name not found, exit
      std::cerr << " LIST OF Measurements " << std::endl;
      for( vmcite = theMeasurementVector.begin(); vmcite != theMeasurementVector.end(); vmcite++ ) {
    std::cerr << (*vmcite)->name() << std::endl;
      }
      std::cerr << "!!EXITING at getMeasurementByName: Measurement " << meas_name << " doesn't exist!!" << std::endl;
      abort();
      //       return (OpticalObject*)0;
    } else {
      return 0;
    }

  } 

  
}

static std::vector<OpticalAlignInfo> Model::getOpticalAlignments ( )

inlinestatic

Definition at line 141 of file Model.h.

References theOpticalAlignments.

Referenced by OpticalObject::createComponentOptOsFromOptAlignInfo().

{ return theOpticalAlignments; }

OpticalObject * Model::getOptOByName ( const ALIstring &  opto_name )

static

--— Find an OptO name in theOptOList and return a pointer to it

Definition at line 578 of file Model.cc.

References begin, dtNoiseDBValidation_cfg::cerr, gather_cfg::cout, ALIUtils::debug, and cmsRelvalreport::exit.

Referenced by Measurement::buildOptOList(), Fit::findEntryFitPosition(), and FittedEntriesSet::SetOptOEntries().

{
  //---------- Look for Optical Object name in OptOList
  std::vector< OpticalObject* >::const_iterator vocite;
  for( vocite = OptOList().begin(); vocite != OptOList().end(); vocite++ ) {
    if( (*vocite)->name() == opto_name ) break;
  }

  if ( vocite == OptOList().end() ) {
    //---------- If opto_name not found, exit
    std::cerr << " LIST OF OpticalObjects " << std::endl;
    for( vocite = OptOList().begin(); vocite != OptOList().end(); vocite++ ) {
      std::cerr <<  (*vocite)->name() << std::endl;
    }
    std::cerr << "!!EXITING at getOptOByName: Optical Object " << opto_name << " doesn't exist!!" << std::endl;
    exit(4);
       //       return (OpticalObject*)0;
  } else {
    //---------- If opto_name found, return pointer to it
    if( ALIUtils::debug>999) {
      std::cout << opto_name.c_str() << "SSOptOitem" << (*vocite) << (*vocite)->name() << "len" << OptOList().size() << std::endl;
    }
    return (*vocite);
  } 

}

OpticalObject * Model::getOptOByType ( const ALIstring &  type )

static

--— Find the first OptO of type 'opto_type' in theOptOList and return a pointer to it

Definition at line 610 of file Model.cc.

References begin, dtNoiseDBValidation_cfg::cerr, and cmsRelvalreport::exit.

{
  //---------- Look for Optical Object type in OptOList
  std::vector< OpticalObject* >::const_iterator vocite;
  for( vocite = OptOList().begin(); vocite != OptOList().end(); vocite++ ) {
   //   std::cout << "OPTOList" << (*msocite).first << std::endl;
    if( (*vocite)->type() == opto_type ) break;
  }

  if ( vocite == OptOList().end() ) {
    //---------- If opto_type not found, exit
    std::cerr << "!!EXITING at getOptOByType: Optical Object " << opto_type << " doesn't exist!!" << std::endl;
    exit(4);
  } else {
    //---------- If opto_type found, return pointer to it
    return (*vocite);
 } 

}

ALIint Model::getParameterValue ( const ALIstring &  sstr, ALIdouble &  val  )

static

************ ACCESS INFO FROM STATIC DATA

--— Search a string in theParameters and return 1 if found

Definition at line 565 of file Model.cc.

References ParameterMgr::getInstance(), and ParameterMgr::getParameterValue().

{
  ParameterMgr* parmgr = ParameterMgr::getInstance();
  ALIint iret = parmgr->getParameterValue( sstr, val );

  return iret;

}

ALIdouble Model::getParamFittedSigmaVectorItem ( const ALIuint  position )

static

Definition at line 1319 of file Model.cc.

References dtNoiseDBValidation_cfg::cerr, cmsRelvalreport::exit, and position.

{
  if( position >= theParamFittedSigmaVector.size() ) {
    std::cerr << "!!EXITING at getParamFittedSigma: position" << position << 
 " bigger than dimension of theParamFittedSigmaVector " << theParamFittedSigmaVector.size() << std::endl;
    exit(3);
  }
  std::vector<ALIdouble>::const_iterator vdcite = theParamFittedSigmaVector.begin() + position;
  return (*vdcite);

}

static ALIstring& Model::MatricesFName ( )

inlinestatic

the name of the File for storing the matrices

Definition at line 101 of file Model.h.

References theMatricesFName.

Referenced by Fit::dumpMatrices().

                                   {
    return theMatricesFName;
  }

static ALIstring& Model::MeasFName ( )

inlinestatic

the name of the Measurements File

Definition at line 91 of file Model.h.

References theMeasFName.

                               {
    return theMeasFName;
  }

static std::vector< Measurement* >& Model::MeasurementList ( )

inlinestatic

Definition at line 79 of file Model.h.

References theMeasurementVector.

Referenced by CocoaDaqReaderRoot::BuildMeasurementsFromOptAlign(), Fit::calculateSimulatedMeasurementsWithOriginalValues(), Fit::CheckIfFitPossible(), Fit::CreateMatrices(), OptOSensor2D::fastTraversesLightRay(), NtupleManager::FillChi2(), Fit::FillMatricesWithCalibratedParameters(), Fit::FillMatricesWithMeasurements(), NtupleManager::FillMeasurements(), Fit::GetMeasurementName(), CocoaToDDLMgr::measurementsAsSpecPars(), Measurement::postConstruct(), Fit::PrintChi2(), readMeasurementsFromFile(), and CocoaDaqReaderText::ReadNextEvent().

                                                      {
     return theMeasurementVector;
  }

static struct tm& Model::MeasurementsTime ( )

inlinestatic

Definition at line 137 of file Model.h.

References theMeasurementsTime.

                                       {
    return theMeasurementsTime;
  }

OpticalObject * Model::nextOptOToCopy ( )

static

Definition at line 816 of file Model.cc.

References gather_cfg::cout, and ALIUtils::debug.

Referenced by OpticalObject::copyData().

{
  if(ALIUtils::debug >= 5) std::cout << "entering nextOptOToCopy() " << std::endl;
  theOptOsToCopyListIterator++;
  //  if(ALIUtils::debug >= 5) std::cout <<" nextOptOToCopy " << (*(theOptOsToCopyListIterator-1))->name() << std::endl;
  return *(theOptOsToCopyListIterator-1);
}

static std::vector< std::vector<ALIstring> >& Model::OptODictionary ( )

inlinestatic

ACCESS STATIC DATA MEMBERS.

Definition at line 67 of file Model.h.

References theOptODictionary.

Referenced by readSystemDescription().

                                                             {
      return theOptODictionary;
  }

static std::vector< OpticalObject* >& Model::OptOList ( )

inlinestatic

Definition at line 71 of file Model.h.

References theOptOList.

Referenced by CocoaDBMgr::BuildAlignments(), CocoaDBMgr::BuildOpticalAlignments(), OpticalObject::createComponentOptOs(), OpticalObject::createComponentOptOsFromOptAlignInfo(), Fit::dumpFittedValues(), Fit::dumpFittedValuesInAllAncestorFrames(), OpticalObjectMgr::dumpOptOs(), NtupleManager::FillOptObjects(), Fit::fitNextEvent(), FittedEntriesManager::GetDifferentBetweenLasers(), CocoaToDDLMgr::writeLogicalVolumes(), CocoaToDDLMgr::writeMaterials(), CocoaToDDLMgr::writePhysicalVolumes(), CocoaToDDLMgr::writeSolids(), CocoaToDDLMgr::writeSpecPars(), and CocoaToDDLMgr::writeSpecParsCocoa().

                                                {
      return theOptOList;
  }

std::string Model::printCocoaStatus ( const cocoaStatus  cs )

static

Definition at line 1553 of file Model.cc.

References COCOA_FitCannotImprove, COCOA_FitChi2Worsened, COCOA_FitImproving, COCOA_FitMatrixNonInversable, COCOA_FitOK, COCOA_Init, COCOA_InitFit, and COCOA_ReadingModel.

{
  std::string str ="";

  if( cs == COCOA_Init ) {
    str = "COCOA_Init ";
  }else if( cs == COCOA_ReadingModel) {
    str = "COCOA_ReadingModel";
  }else if( cs == COCOA_InitFit) {
    str = "COCOA_InitFit";
  }else if( cs == COCOA_FitOK) {
    str = "COCOA_FitOK";
  }else if( cs == COCOA_FitImproving) {
    str = "COCOA_FitImproving";
  }else if( cs == COCOA_FitCannotImprove) {
    str = "COCOA_FitCannotImprove";
  }else if( cs == COCOA_FitChi2Worsened) {
    str = "COCOA_FitChi2Worsened";
  }else if( cs == COCOA_FitMatrixNonInversable) {
    str = "COCOA_FitMatrixNonInversable";
  }

  return str;
}

ALIbool Model::readMeasurementsFromFile ( ALIstring  only1Date = ALIstring(""), ALIstring  only1Time = ALIstring("")  )

static

Definition at line 1336 of file Model.cc.

References begin, dtNoiseDBValidation_cfg::cerr, Measurement::correctValueAndSigma(), gather_cfg::cout, ALIUtils::debug, Measurement::dim(), ALIUtils::dumpVS(), end, ALIFileIn::eof(), ALIFileIn::ErrorInLine(), cmsRelvalreport::exit, Measurement::fillData(), GlobalOptionMgr::getInstance(), ALIFileIn::getInstance(), ALIFileIn::getWordsInLine(), GlobalOptionMgr::GlobalOptions(), MeasurementList(), Measurement::measurementsFileName(), Measurement::only1, Measurement::only1Date, Measurement::only1Time, Measurement::setCurrentDate(), Measurement::setSigma(), Measurement::sigma(), and Measurement::valueType().

{
 if(ALIUtils::debug >= 5) std::cout << " readMeasurementsFromFile " << Measurement::measurementsFileName() << std::endl;
  if( Measurement::measurementsFileName() == "") return 1;

  ALIFileIn& filein = ALIFileIn::getInstance( Measurement::measurementsFileName() );
  std::vector<ALIstring> wordlist;
  //---------- read date
  //  ALIint retfil = filein.getWordsInLine(wordlist);
  // std::cout << "@@@@@@@@@@@@@@@ RETFIL " << retfil << std::endl;
  //if( retfil == 0 ) {
  if( filein.getWordsInLine(wordlist) == 0 ) {
    if(ALIUtils::debug>=4 ) std::cout << "@@@@ No more measurements left" << std::endl;
    return 0; 
  }

  //  struct tm tim;
  //t Model::setMeasurementsTime( tim );

  //if you are looking for only one measurement with a given date and time, loop all measurements until you find it
  if( Measurement::only1 ) {
    for(;;) {
      if( wordlist[0] == "DATE:" && wordlist[1] == Measurement::only1Date && wordlist[2] == Measurement::only1Time ) break;
      filein.getWordsInLine(wordlist); 
      if(filein.eof() ) {
    std::cerr << "!! EXITING date not found in measurements file" << Measurement::only1Date << " " <<  Measurement::only1Time << std::endl;
    exit(1);
      }
    }
  }

  //set date and time of current measurement
  if( wordlist[0] == "DATE:" ) {
    Measurement::setCurrentDate( wordlist ); 
  } 

  //---------- loop measurements
  ALIint nMeas = Model::MeasurementList().size();
  if(ALIUtils::debug >= 4) {
    std::cout << " Reading " << nMeas << " measurements from file " << Measurement::measurementsFileName() 
         << " DATE: " << wordlist[1] << " " << wordlist[1] << std::endl;
  }
  ALIint ii;
  for(ii = 0; ii < nMeas; ii++) {
    filein.getWordsInLine(wordlist);  
    if( wordlist[0] == ALIstring("SENSOR2D") || wordlist[0] == ALIstring("TILTMETER") || wordlist[0] == ALIstring("DISTANCEMETER")  || wordlist[0] == ALIstring("DISTANCEMETER1DIM")  || wordlist[0] == ALIstring("COPS") ) {
      if( wordlist.size() != 2 ) {
    std::cerr << "!!!EXITING Model::readMeasurementsFromFile. number of words should be 2 instead of " << wordlist.size() << std::endl;
    ALIUtils::dumpVS( wordlist, " " );
    exit(1);
      }
      std::vector< Measurement* >::const_iterator vmcite;
      for( vmcite = MeasurementList().begin();  vmcite != MeasurementList().end(); vmcite++ ) {
    //-------- Measurement found, fill data
    /*  ALIint last_slash =  (*vmcite)->name().rfind('/');
    ALIstring oname = (*vmcite)->name();
    if( last_slash != -1 ) {
      oname = oname.substr(last_slash+1, (*vmcite)->name().size()-1);
      } 
    */
    ALIint fcolon = (*vmcite)->name().find(':');
    ALIstring oname = (*vmcite)->name();
    oname = oname.substr(fcolon+1,oname.length());
    //-    std::cout << " measurement name " << (*vmcite)->name() << " short " << oname << std::endl;
    if( oname == wordlist[1] ) {
      //-   std::cout << " measurement name found " << oname << std::endl;
      if( (*vmcite)->type() != wordlist[0] ) {
        std::cerr << "!!! Reading measurement from file: type in file is " 
         << wordlist[0] << " and should be " << (*vmcite)->type() << std::endl;
        exit(1);
      }
      Measurement* meastemp = *vmcite;
      
      GlobalOptionMgr* gomgr = GlobalOptionMgr::getInstance();
      ALIbool sigmaFF = gomgr->GlobalOptions()["measurementErrorFromFile"];
      //---------- Read the data 
      for ( ALIuint ii=0; ii < meastemp->dim(); ii++){
        filein.getWordsInLine( wordlist );
            ALIdouble sigma = 0.;
            if( !sigmaFF ) { 
          // keep the sigma, do not read it from file 
          const ALIdouble* sigmav = meastemp->sigma();
          sigma = sigmav[ii];
        }
        //---- Check measurement value type is OK
        if( meastemp->valueType(ii) != wordlist[0] ) {
          filein.ErrorInLine();
          std::cerr << "!!!FATAL ERROR: Measurement value type is " << wordlist[0] << " while in setup definition was " << meastemp->valueType(ii) << std::endl;
          exit(1);
        }
        meastemp->fillData( ii, wordlist );
            if( !sigmaFF ) { 
          meastemp->setSigma( ii, sigma );
        }
      }
      meastemp->correctValueAndSigma();
      break;
    }
      }
      if( vmcite == MeasurementList().end() ) {
    for( vmcite = MeasurementList().begin(); vmcite != MeasurementList().end(); vmcite++ ) {
      std::cerr << "MEAS: " << (*vmcite)->name() << " " << (*vmcite)->type() << std::endl;
    }
    std::cerr << "!!! Reading measurement from file: measurement not found in list: type in file is "  << wordlist[1]  << std::endl;
    exit(1);
      }
    } else {
      std::cerr << " wrong type of measurement: " << wordlist[0] << std::endl
       << " Available types are SENSOR2D, TILTMETER, DISTANCEMETER, DISTANCEMETER1DIM, COPS" << std::endl;
      exit(1);
    }
  }
  //-  std::cout << " returning readmeasff" << std::endl;

  return 1;
}

void Model::readSystemDescription ( )

static

-------— Read the different sections of the SDF and act accordingly

Definition at line 101 of file Model.cc.

References ParameterMgr::addParameter(), ParameterMgr::addRandomFlatParameter(), ParameterMgr::addRandomGaussParameter(), begin, buildMeasurementsLinksToOptOs(), dtNoiseDBValidation_cfg::cerr, ALIFileIn::close(), COCOA_Init, Measurement::construct(), OpticalObject::construct(), filterCSVwithJSON::copy, gather_cfg::cout, ALIUtils::debug, ALIUtils::dumpVS(), end, ALIFileIn::eof(), ALIFileIn::ErrorInLine(), cmsRelvalreport::exit, spr::find(), ALIUtils::getFloat(), ParameterMgr::getInstance(), EntryMgr::getInstance(), ErrorCorrelationMgr::getInstance(), GlobalOptionMgr::getInstance(), ALIFileIn::getInstance(), ALIUtils::getInt(), ALIFileIn::getWordsInLine(), ALIUtils::IsNumber(), ALIFileIn::nline(), Measurement::only1, Measurement::only1Date, Measurement::only1Time, OptODictionary(), EntryMgr::readEntryFromReportOut(), ErrorCorrelationMgr::readFromReportFile(), SDFName(), sectGlobalOptions, sectMeasurements, sectParameters, sectReportOut, sectSystemTreeData, sectSystemTreeDescription, ALIUtils::SetAngleDimensionFactors(), setCocoaStatus(), Measurement::setConversionFactor(), GlobalOptionMgr::setGlobalOption(), ALIUtils::SetLengthDimensionFactors(), ALIUtils::SetOutputAngleDimensionFactors(), ALIUtils::SetOutputLengthDimensionFactors(), ParameterMgr::setRandomSeed(), and findQualityFiles::size.

{
  Model::setCocoaStatus( COCOA_Init );
  
  ALIint data_reading = 0;  // variable to identify the first line of section SYSTEM_TREE_DATA 
  
  //---------- Open the system description file
  ALIFileIn& filein = ALIFileIn::getInstance( Model::SDFName());
 
  //----------- Set section titles
  std::vector<ALIstring> SectionTitle;
  SectionTitle.push_back(ALIstring("GLOBAL_OPTIONS"));
  SectionTitle.push_back(ALIstring("PARAMETERS"));
  SectionTitle.push_back(ALIstring("SYSTEM_TREE_DESCRIPTION"));
  SectionTitle.push_back(ALIstring("SYSTEM_TREE_DATA"));
  SectionTitle.push_back(ALIstring("MEASUREMENTS"));
  SectionTitle.push_back(ALIstring("REPORT.OUT"));
  std::vector<ALIstring>::iterator SectionTitleIterator;

  //---------------------------------------- Loops lines in SDF file  
  std::vector<ALIstring> wordlist;
  ALIint InSectionNo = -1; 
  ALIint currentSectionNo = -1; 
  while (!filein.eof()) {

    if( !filein.getWordsInLine(wordlist) ) break;  //----- Read line
    assert( wordlist.size() != 0 );

    //----- checking
    if( ALIUtils::debug > 99) {
      ALIUtils::dumpVS( wordlist, " ", std::cout );
    }

    //---------- Get in which section the current line is and act accordingly
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@@ ---------- Check if word is start of section 
    SectionTitleIterator = find(SectionTitle.begin(), SectionTitle.end(), *wordlist.begin());
    if( SectionTitleIterator != SectionTitle.end() ) { 
      //---------- Check that previous sections are correct
      currentSectionNo = SectionTitleIterator - SectionTitle.begin();
      if ( currentSectionNo != InSectionNo + 1 ) {
    if( currentSectionNo != sectReportOut ) {
      ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
      std::cerr << "BAD ORDER OF SECTIONS, reading section " 
           << *SectionTitleIterator << std::endl 
           << " currentSectionNo = " << currentSectionNo << " InSectionNo = " << InSectionNo << std::endl
           << " ---------  Please see documentation  ---------- " << std::endl;
      exit(1);
    }
      } else {
    if( currentSectionNo != sectReportOut ) {
      InSectionNo++;
    }
      }
      if( currentSectionNo == sectMeasurements ) {
    SetValueDisplacementsFromReportOut();
      }

      if( ALIUtils::debug >= 4) std::cout << std::endl << "START OF SECTION: " << currentSectionNo << " " << *SectionTitleIterator << std::endl;

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@@ ---------- Reading GLOBAL_OPTIONS section
    } else if( currentSectionNo == sectGlobalOptions ) {
      //-       std::cout << "utilsdebug " << ALIUtils::debug << std::endl;
      //-------- Check format of line read
      //----- Two words
      if( wordlist.size() == 2 ) {
    //----- Second word is number
    int isnumber = ALIUtils::IsNumber( wordlist[1] );
    if (!isnumber && wordlist[0] != ALIstring("external_meas")) {
      ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
      std::cerr << ": EXPECTING A NUMBER, FOUND: " << wordlist[1] << std::endl;
      exit(2);
    }
    
    GlobalOptionMgr* gomgr = GlobalOptionMgr::getInstance();
    gomgr->setGlobalOption( wordlist[0], ALIUtils::getFloat(wordlist[1]), ALIFileIn::getInstance(Model::SDFName()) );
 
      //-or    GlobalOptions.insert(std::pair<const ALIstring, ALIdouble>(wordlist[0] , atof(wordlist[1].c_str()) ) );
    
    if( ALIUtils::debug >= 1) {
      ALIUtils::dumpVS( wordlist, "GLOBAL_OPTION:  ", std::cout );
    }

      } else {
    std::cout << "error < 1" << std::endl;
    ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
        std::cerr << ": IN GLOBAL_OPTIONS section TWO-WORD LINES ARE EXPECTED " << std::endl;
        exit(2);
      }
    
      //------- Set dimension factors for lengths and angles
      ALIUtils::SetLengthDimensionFactors();
      ALIUtils::SetAngleDimensionFactors();
      ALIUtils::SetOutputLengthDimensionFactors();
      ALIUtils::SetOutputAngleDimensionFactors();

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@@ ---------- Reading PARAMETERS section
    } else if( currentSectionNo == sectParameters ||
           currentSectionNo == -1 ) {  // Allow parameters in first lines (for easy study of different parameter values)

      //---------- Check format of line read
      //----- Two words
      if( wordlist.size() == 2 ) {
    /* now checked in ParameterMgr
    //----- Second is number
    int isnumber = ALIUtils::IsNumber( wordlist[1] );
    if( !isnumber ) {
      ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
      std::cerr << ": EXPECTING A NUMBER, FOUND: " << wordlist[1] << std::endl;
      exit(2);
      } */
    
    //old---------- Create parameter with value read (or overwrite existing value)
    //---------- Create parameter with value read if parameter does not exist yet
    ParameterMgr* parmgr = ParameterMgr::getInstance();
    parmgr->addParameter( wordlist[0], wordlist[1] );
    
      } else if( wordlist.size() == 3 ) {
    if( wordlist[1] != "seed" ) { 
      ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
      std::cerr << ": For a three-word parameter line, second has to be 'seed', it is  " << wordlist[1] << std::endl;
      exit(1);
    } 

    if( wordlist[0] == "gauss" || wordlist[0] == "flat" ){ 
      ParameterMgr::getInstance()->setRandomSeed( ALIUtils::getInt( wordlist[2] ));
    } else {
      ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
      std::cerr << ": For a three-word parameter line, first has to be 'gauss' or 'flat', it is  " << wordlist[0] << std::endl;
      exit(1);
    } 
      } else if( wordlist.size() == 4 ) {
    if( wordlist[0] == "gauss" ) {  
      ParameterMgr::getInstance()->addRandomGaussParameter( wordlist[1], wordlist[2], wordlist[3] ); 
    } else if( wordlist[0] == "flat" ) {
     ParameterMgr::getInstance()->addRandomFlatParameter( wordlist[1], wordlist[2], wordlist[3] ); 
    } else {
      ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
      std::cerr << ": For a four-word parameter line, first has to be 'gauss' or 'flat', it is  " << wordlist[0] << std::endl;
      exit(1);
    } 
      } else {
    ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
        std::cerr << ": IN PARAMETERS section TWO-WORD-LINES ARE EXPECTED " << std::endl;
    ALIUtils::dumpVS( wordlist , " " );
        exit(2);
      } 

      //print it out
      if( ALIUtils::debug >= 1) { 
    ALIUtils::dumpVS( wordlist, "PARAMETERS:  ", std::cout );
      }
    
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@@ ---------- Reading SYSTEM_TREE_DESCRIPTION section
    } else if( currentSectionNo == sectSystemTreeDescription ) {

      //---------- Convert the numbers found in list of components (e.g. 2 laser -> laser laser )
      //----- Backup wordlist and erase it
      std::vector<ALIstring> wordlist2 = wordlist;
      std::vector<ALIstring>::iterator vsite;
      ALIint wsiz = wordlist.size();
      for ( ALIint ii = 0; ii < wsiz; ii++) {
        wordlist.pop_back();
      }
      //----- Loop objects looking for numbers to convert
      for ( vsite=wordlist2.begin(); vsite!=wordlist2.end(); vsite++ ) {
        if( ALIUtils::IsNumber( *vsite ) ) {
          int nOpticalObjects = atoi( (*vsite).c_str() );
      //----- If number is 1 it is not necessary
          if( nOpticalObjects == 1 ) {
            if( ALIUtils::debug >= 0) std::cerr << "WARNING: in line " << ALIFileIn::getInstance( Model::SDFName() ).nline() 
                 << " number of repeated OpticalObjects = 1. Please avoid the numbering " << std::endl; 
        //-    wordlist.erase( wordlist.begin() + (vsite-wordlist2.begin()) );
      } else {
        //----- The number cannot be the last in the list and you cannot put two numbers together
        if( vsite+1 == wordlist.end() || ALIUtils::IsNumber(*(vsite+1)) ) {
          ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
              std::cerr << "NUMBER SHOULD BE FOLLOWED BY an OpticalObject type" << std::endl;
              exit(2);
            }
        //----- Right format: convert
            //--- Substitute the number by the object type in wordlist
        //-   *( wordlist.begin() + (vsite-wordlist2.begin()) ) = *(vsite+1);
        //--- Add n-1 object types to wordlist (the nth object will be added as the object taht goes after the number)
            for ( ALIint ii=0; ii<nOpticalObjects-1; ii++ ) {
          //-std::cout << ii << "inserting in wordlist " << *(vsite+1) << std::endl;
              wordlist.push_back( *(vsite+1) ); 
            }
          }
        } else {
      //----- Not number, add it to wordlist
      wordlist.push_back( *vsite );
    }
      }

      //---------- Dump system structure
      if ( ALIUtils::debug >= 1 ) {
    ALIUtils::dumpVS( wordlist, "SYSTEM TREE DESCRIPTION: before ordering, OBJECT: ", std::cout );
      }

      //---------- Fill the list of Optical Object with components (theOptODictionary)
      //----- First word is 'object': new OptO
      if ( wordlist[0] == ALIstring("object") ) {
      //----- Check out repeated objects
        std::vector< std::vector<ALIstring> >::iterator vvsite; 
        for ( vvsite = theOptODictionary.begin(); vvsite != theOptODictionary.end(); vvsite++) {
      //-     std::cout << " system" << vvsite << std::endl;

      if( *( (*vvsite).begin() ) == wordlist[1] ) {
        ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
        std::cerr << "SYSTEM_TREE_DESCRIPTION: REPEATED object " << *((*vvsite).begin()) << " ( NOT ALLOWED NEITHER WITH EQUAL NOR WITH DIFFERENT COMPONENTS)" << std::endl;
        exit(1);
      }
        }
    //------- Add an item to theOptODictionary
        std::vector<ALIstring> vstemp;
        copy( wordlist.begin()+1, wordlist.end(),back_inserter( vstemp ) );
        Model::OptODictionary().push_back(vstemp);
      } else {
        //----- First word is not 'object': add to previous OptO
        if(OptODictionary().size() == 0) {
      ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
      std::cerr << "SYSTEM_TREE_DESCRIPTION section: FIRST LINE SHOULD START WITH 'object'" << std::endl;
          exit(2);
        }
        copy(wordlist.begin(), wordlist.end(),
             back_inserter( *(OptODictionary().end()-1) ) );
      }
 
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//---------------------------------- Reading SYSTEM_TREE_DATA section
    } else if( currentSectionNo == sectSystemTreeData ) {
      //---------- First line of data:
      if ( !data_reading ) {
    //        ALIint dictsize = theOptODictionary.size();
        //------- Reorder theOptODictionary
        std::vector< std::vector<ALIstring> > OptODictionary2; 
        reorderOptODictionary( "system", OptODictionary2);
        if(OptODictionary2.size() == 0) {
      std::cerr << "SYSTEM_TREE_DESCRIPTION section: no object 'system' found " << std::endl;
      exit(9);
    }
    //------- We start reordering with 'system', therefore if 'system' is not the top most object, the objects not hanging from it would not be considered
    //---- Check if all the objects are here
        std::vector< std::vector<ALIstring> >::const_iterator vvscite, vvscite2; 
    //  ALIint dictsizen = 0;
        for( vvscite = theOptODictionary.begin(); vvscite != theOptODictionary.end(); vvscite++) {
          ALIbool ofound = 0;
      for( vvscite2 = OptODictionary2.begin(); vvscite2 != OptODictionary2.end(); vvscite2++) {
        if( *( (*vvscite).begin() ) == *( (*vvscite2).begin() ) ) {
              ofound = 1;
          break;
        }
      }
      if( !ofound  ) {
        std::cerr << "!!!SYSTEM_TREE_DESCRIPTION section: object " << *( (*vvscite).begin()) << " is not hanging from object 'system' " << std::endl;
        for (vvscite=OptODictionary().begin();vvscite!=OptODictionary().end();vvscite++){
          std::vector<ALIstring> ptemp = *vvscite;
          ALIUtils::dumpVS( ptemp, "OBJECT ", std::cerr );
        }       
        exit(9);
      }
    }
        theOptODictionary = OptODictionary2;

        data_reading = 1;

        //------- Dump ordered OptOs
        if ( ALIUtils::debug >= 3) {
          std::vector< std::vector<ALIstring> >::iterator itevs;
          for (itevs=OptODictionary().begin();itevs!=OptODictionary().end();itevs++){
            std::vector<ALIstring> ptemp = *itevs;
        ALIUtils::dumpVS( ptemp, " SYSTEM TREE DESCRIPTION: after ordering: OBJECT ", std::cout );
          }
        }

        //---------- Create OpticalObject 'system' (first OpticalObject object): 
    //---------- it will create its components and recursively all the System Tree of Optical Objects
        if(wordlist[0] != "system" ) {
          std::cerr << "SYSTEM_TREE_DATA section: object 'system' is not the first one " << std::endl;
      exit(9);
    }

        OpticalObject* OptOsystem = new OpticalObject( 0, "system", wordlist[1], 0 );
        OptOsystem->construct();
        //-              Model::_OptOtree.insert( std::multimap< ALIstring, OpticalObject*, std::less<ALIstring> >::value_type(OptOsystem->type(), OptOsystem) );
        //              theOptOlist[OptOsystem->name()] = OptOsystem; 
        theOptOList.push_back( OptOsystem ); 
       
      } else {      
        //----------- All system is read by the Optical Objects, it should not reach here
    ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
    std::cerr << " STILL SOME LINES AFTER ALL SYSTEM TREE IS READ!!!"<< std::endl;
        exit(9);
      }
 
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//----------------------------------- Reading MEASUREMENTS section
    } else if( currentSectionNo == sectMeasurements ) {
      //---------- Create Measurement with appropiate dimension
      Measurement* meastemp = 0;
      ALIstring measType = wordlist[0];
      ALIstring measName;
      if( wordlist.size() == 2 ) {
    measName = wordlist[1];
        wordlist.pop_back();
      } else {
    measName = "";
      }
      if ( measType == ALIstring("SENSOR2D") ) {
        meastemp = new MeasurementSensor2D( 2, measType, measName );      
        meastemp->setConversionFactor( wordlist );
    meastemp->construct();
      } else if ( measType == ALIstring("DISTANCEMETER3DIM") ) {
    meastemp = new MeasurementDistancemeter3dim( 1, measType, measName );
    meastemp->setConversionFactor( wordlist );
    meastemp->construct();
      } else if ( measType == ALIstring("DISTANCEMETER") ||
          measType == ALIstring("DISTANCEMETER1DIM") ) {
        meastemp = new MeasurementDistancemeter( 1, measType, measName );
        meastemp->setConversionFactor( wordlist );
    meastemp->construct();
      } else if ( measType == ALIstring("TILTMETER") ) {
        meastemp = new MeasurementTiltmeter( 1, measType, measName );
        meastemp->setConversionFactor( wordlist );
    meastemp->construct();
      } else if ( measType == ALIstring("COPS") ) {
        meastemp = new MeasurementCOPS( 4, measType, measName );
        meastemp->setConversionFactor( wordlist );
    meastemp->construct();
      } else if ( measType == ALIstring("DIFFCENTRE") ) {
 //t       meastemp = new MeasurementDiffCentre( 1, measType, measName );
    meastemp->construct();
      } else if ( measType == ALIstring("DIFFANGLE") ) {
//t        meastemp = new MeasurementDiffAngle( 2, measType, measName );
    meastemp->construct();
      } else if ( measType == ALIstring("DIFFENTRY") ) {
        meastemp = new MeasurementDiffEntry( 1, measType, measName );
    meastemp->construct();
      } else if ( measType == ALIstring("measurements_from_file") || measType == ALIstring("@measurements_from_file") ) {
    new CocoaDaqReaderText( wordlist[1] );
    //m Measurement::setMeasurementsFileName( wordlist[1] );
    //m if ( ALIUtils::debug >= 2) std::cout << " setting measurements_from_file " << measType << " == " << Measurement::measurementsFileName() << std::endl;
        if( wordlist.size() == 4) {
      Measurement::only1 = 1; 
      Measurement::only1Date = wordlist[2]; 
      Measurement::only1Time = wordlist[3]; 
      //-      std::cout << " setting Measurement::only1" <<  Measurement::only1 << std::endl;
    }
      } else if ( measType == ALIstring("measurements_from_file_ROOT") || measType == ALIstring("@measurements_from_file") ) {
    new CocoaDaqReaderRoot( wordlist[1] );
      } else if ( wordlist[0] == ALIstring("correlations_from_file") || wordlist[0] == ALIstring("@correlations_from_file") ) {
    ErrorCorrelationMgr::getInstance()->readFromReportFile( wordlist[1] );
      } else if ( wordlist[0] == ALIstring("copy_measurements") || wordlist[0] == ALIstring("@copy_measurements") ) {
    copyMeasurements( wordlist );
    //      } else if( wordlist[0] == "scan" || wordlist[0] == "@scan" ) {
    //  ScanMgr::getInstance()->addOptOEntry( wordlist );
      } else if ( wordlist[0] == ALIstring("fittedEntries_from_file") ) {
    theFittedEntriesReader = new FittedEntriesReader( wordlist[1] );
     if ( ALIUtils::debug >= 2) std::cout << " setting fittedEntries_from_file " << wordlist[0] << " == " << wordlist[1] << std::endl;
      } else {
        std::cerr << "Measurement:" << std::endl;
    ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
    std::cerr << "!!! type of measurement not allowed: " << wordlist[0] 
         << std::endl;
    std::cerr << " Allowed types: SENSOR2D, DISTANCEMETER, DISTANCEMETER1DIM, TILTMETER, COPS, DIFFENTRY " << std::endl; 
        exit(2);
      }
//-      meastemp->setGlobalName( wordlist[0] );
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@@ ---------- Reading REPORT OUT  section
    } else if( currentSectionNo == sectReportOut ) {
      //----- It must be after global options section
      if ( InSectionNo+1 != sectParameters ) {
    ALIFileIn::getInstance( Model::SDFName() ).ErrorInLine();
    std::cerr << "BAD ORDER OF SECTIONS, reading section " 
         << *SectionTitleIterator << std::endl 
         << " currentSectionNo = " << currentSectionNo << " InSectionNo = " << InSectionNo << std::endl
         << " ---------  Please see documentation  ---------- " << std::endl;
    exit(1);
      }

      EntryMgr::getInstance()->readEntryFromReportOut( wordlist );
    }
  }


//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@@ All system read
  //---------- make Measurements links to their OptOs 
  if ( ALIUtils::debug >= 2) std::cout << std::endl << "@@@@ Building Measurements links to OptOs" << std::endl;
  Model::buildMeasurementsLinksToOptOs();

  if ( ALIUtils::debug >= 1) {
      std::cout << "----------  SYSTEM SUCCESFULLY READ ----------" << std::endl << std::endl;
  }  
  filein.close();

  return;
  
}

void Model::recoverParamFittedSigma ( const ALIstring &  opto_name, const ALIstring &  entry_name, const ALIuint  position  )

staticprivate

Definition at line 1298 of file Model.cc.

{
  /*
  if( position >= theParamFittedSigmaVector.size() ) {
    std::cerr << "!!EXITING at recoverParamFittedSigma: position" << position << 
 " bigger than dimension of theParamFittedSigmaVector " << theParamFittedSigmaVector.size() << std::endl;
    exit(3);
  }
  ALIdouble sigma = getParamFittedSigmaVectorItem( position );

  Entry* entry = getEntryByName(opto_name, entry_name);
  entry->setSigma( sigma );
  entry->setQuality( 1 );
  std::cout << "recover " << opto_name << entry_name << entry->sigma() <<std::endl;
  */
}

void Model::reorderOptODictionary ( const ALIstring &  ssearch, std::vector< std::vector< ALIstring > > &  OptODictionary2  )

staticprivate

********** private METHODS

Reorder the list of OptOs in a hierarchical structure (tree-like)

Definition at line 511 of file Model.cc.

References begin.

                                                                                                              {

  //  int ALIstring_found = 0;
  std::vector< std::vector<ALIstring> >::iterator vvsite; 
  std::vector< ALIstring >::iterator vsite; 

//---------------------------------------- Look for 'ssearch' as the first ALIstring of an OptODictionary item
  for( vvsite = OptODictionary().begin(); vvsite != OptODictionary().end(); ++vvsite) {
   if( *( (*vvsite).begin() ) == ssearch ) {
     //     ALIstring_found = 1;
     OptODictionary2.push_back(*vvsite);
 
     //-    std::cout << "VVSITE";
     //-   ostream_iterator<ALIstring> outs(std::cout,"&");
     //-   copy( (*vvsite).begin(), (*vvsite).end(), outs);
 
//---------------------------------- look for components of this _OptODictionary item
     for(  vsite = (*vvsite).begin()+1; vsite != (*vvsite).end(); ++vsite ) {
       reorderOptODictionary( *vsite, OptODictionary2 );
     }
     break;
   }
  }
 
  /*  //------- object 'system' should exist
  if(!ALIstring_found && ssearch == "system") {
    std::cerr << "SYSTEM_TREE_DATA section: no 'object system' found " << std::endl;
    exit(9);
    } */

} 

static ALIstring& Model::ReportFName ( )

inlinestatic

the name of the report File

Definition at line 96 of file Model.h.

References theReportFName.

Referenced by Fit::fitNextEvent(), Fit::getFitQuality(), Fit::PrintChi2(), and CocoaAnalyzer::RunCocoa().

                                 {
    return theReportFName;
  }

void Model::saveParamFittedCorrelation ( const ALIstring &  opto_name1, const ALIstring &  entry_name1, const ALIstring &  opto_name2, const ALIstring &  entry_name2  )

staticprivate

Definition at line 1262 of file Model.cc.

{
  /*
//---------- Get sigma of param Optical
  Entry* entry1 = getEntryByName(opto_name1, entry_name1);
  Entry* entry2 = getEntryByName(opto_name2, entry_name2);
  ALIuint fit_pos1 = entry1->fitPos();
  ALIuint fit_pos2 = entry2->fitPos();

  //---------- Get correlation if entry has been fitted
  ALIuint pmsize = Fit::propagationMatrix().NoLines();
  if( fit_pos1 >= 0 && fit_pos1 < pmsize && fit_pos2 >= 0 && fit_pos2 < pmsize ) {
    ALIdouble error1 = sqrt( Fit::propagationMatrix()( fit_pos1, fit_pos1) );
    ALIdouble error2 = sqrt( Fit::propagationMatrix()( fit_pos2, fit_pos2) );
    ALIdouble correl = Fit::propagationMatrix()( fit_pos1, fit_pos2) / error1 / error2;
    theParamFittedSigmaVector.push_back( correl );
    if(ALIUtils::debug>=9) {
      std::cout  << "saveParamFittedCorre" << opto_name1 << entry_name1 << fit_pos1 << 
        opto_name2 << entry_name2 << fit_pos2 << "MM " << correl << std::endl;
    }
  } else {
    if(ALIUtils::debug>=9) {
      std::cout  << "NOsaveParamFittedCorre" << opto_name1 << entry_name1 << fit_pos1 << 
        opto_name2 << entry_name2 << fit_pos2 << "MM " << std::endl;
    theParamFittedSigmaVector.push_back( 0. );
    }

  }
  //-    Fit::propagationMatrix().Dump(" the 2PropagationMatrix");
  */
}

void Model::saveParamFittedSigma ( const ALIstring &  opto_name, const ALIstring &  entry_name  )

staticprivate

Definition at line 1241 of file Model.cc.

{
  /*
//---------- Get sigma of param Optical
  Entry* entry = getEntryByName(opto_name, entry_name);
  ALIuint fit_pos = entry->fitPos();
  if( fit_pos < 0 || fit_pos >= Fit::propagationMatrix().NoLines()) {
    std::cerr << "!!EXITING at saveParamFittedSigma: fit position incorrect " <<
      fit_pos << "propagationMatrix size =" << Fit::propagationMatrix().NoLines() << opto_name << std::endl;
    exit(3);
  }
  std::cout << entry_name << "saveParamFittedSigma" << fit_pos << sqrt(Fit::propagationMatrix()( fit_pos, fit_pos)) << std::endl;
  theParamFittedSigmaVector.push_back( sqrt(Fit::propagationMatrix()( fit_pos, fit_pos)) );
  //-    Fit::propagationMatrix().Dump(" the 2PropagationMatrix");
  */
}

static ALIstring& Model::SDFName ( )

inlinestatic

the name of the System Description File

Definition at line 86 of file Model.h.

References theSDFName.

Referenced by MeasurementDiffEntry::buildOptONamesList(), Measurement::construct(), OpticalObject::construct(), Entry::fill(), Entry::fillFromInputFileQuality(), Entry::fillFromInputFileSigma(), Entry::fillFromInputFileValue(), and readSystemDescription().

                             {
    return theSDFName;
  }

static void Model::setCocoaStatus ( const cocoaStatus  cs )

inlinestatic

Definition at line 49 of file Model.h.

References fwrapper::cs, and theCocoaStatus.

Referenced by Fit::fitNextEvent(), and readSystemDescription().

{ theCocoaStatus = cs; }

static void Model::setMatricesFName ( const ALIstring &  name )

inlinestatic

Definition at line 164 of file Model.h.

References mergeVDriftHistosByStation::name, and theMatricesFName.

                                                        {
    theMatricesFName = name;
  }

static void Model::setMeasurementsTime ( struct tm &  tim )

inlinestatic

Definition at line 168 of file Model.h.

References theMeasurementsTime.

                                                    {
    theMeasurementsTime = tim;
  }

static void Model::setReportFName ( const ALIstring &  name )

inlinestatic

Definition at line 160 of file Model.h.

References mergeVDriftHistosByStation::name, and theReportFName.

                                                      {
    theReportFName = name;
  }

static void Model::setSDFName ( const ALIstring &  name )

inlinestatic

Definition at line 156 of file Model.h.

References mergeVDriftHistosByStation::name, and theSDFName.

                                                  {
    theSDFName = name;
  }

void Model::SetValueDisplacementsFromReportOut ( )

staticprivate

Definition at line 1530 of file Model.cc.

References gather_cfg::cout, ALIUtils::debug, EntryList(), EntryMgr::findEntryByLongName(), EntryMgr::getInstance(), EntryMgr::numberOfEntries(), and EntryData::valueDisplacement().

{
  if( ALIUtils::debug >= 3 ) std::cout << " Model::SetValueDisplacementsFromReportOut() " << std::endl;
  
  EntryMgr* entryMgr = EntryMgr::getInstance();
  
  if( entryMgr->numberOfEntries() != 0 ){
    EntryData* entryData;
    
    std::vector<Entry*>::const_iterator vecite; 
    for ( vecite = Model::EntryList().begin();
      vecite != Model::EntryList().end(); vecite++ ) {
      //----- Find the EntryData corresponding to this entry
      entryData = entryMgr->findEntryByLongName( (*vecite)->OptOCurrent()->longName(), (*vecite)->name() );
      if( ALIUtils::debug >= 3 ) std::cout << "SetValueDisplacementsFromReportOut "<< (*vecite)->OptOCurrent()->longName() << " " <<  (*vecite)->name() << " " << entryData->valueDisplacement() << std::endl;
      (*vecite)->addFittedDisplacementToValue( entryData->valueDisplacement() );
    }   
  }
  
}

Member Data Documentation

ALIint Model::CMSLinkIteration = 0

staticprivate

Definition at line 253 of file Model.h.

Referenced by CMSLinkDeleteOptOs().

std::vector< ALIdouble > Model::CMSLinkRangeDetValue

static

Definition at line 239 of file Model.h.

ALIint Model::Ncmslinkrange = 0

static

*************** FOR RANGE STUDIES

Definition at line 238 of file Model.h.

cocoaStatus Model::theCocoaStatus = COCOA_Init

staticprivate

Definition at line 286 of file Model.h.

Referenced by getCocoaStatus(), and setCocoaStatus().

std::vector< Entry * > Model::theEntryVector

staticprivate

std::vector of all Entries

Definition at line 206 of file Model.h.

Referenced by addEntryToList(), and EntryList().

FittedEntriesReader * Model::theFittedEntriesReader = 0

staticprivate

Definition at line 294 of file Model.h.

Referenced by getFittedEntriesReader().

Model * Model::theInstance = 0

staticprivate

Only instance of Model.

Definition at line 191 of file Model.h.

ALIstring Model::theMatricesFName = "matrices.out"

staticprivate

the name of the File for storing the matrices

Definition at line 218 of file Model.h.

Referenced by MatricesFName(), and setMatricesFName().

ALIstring Model::theMeasFName = "Measurements.txt"

staticprivate

the name of the Measurements File

Definition at line 214 of file Model.h.

Referenced by MeasFName().

struct tm Model::theMeasurementsTime

staticprivate

Definition at line 292 of file Model.h.

Referenced by MeasurementsTime(), and setMeasurementsTime().

std::vector< Measurement * > Model::theMeasurementVector

staticprivate

std::vector of all Measurements

Definition at line 209 of file Model.h.

Referenced by addMeasurementToList(), and MeasurementList().

std::vector< OpticalAlignInfo > Model::theOpticalAlignments

staticprivate

Definition at line 296 of file Model.h.

Referenced by getOpticalAlignments().

std::vector< std::vector< ALIstring > > Model::theOptODictionary

staticprivate

parameters

std::vector of OptOs with components (in tree structure)

Definition at line 197 of file Model.h.

Referenced by OptODictionary().

std::vector< OpticalObject * > Model::theOptOList

staticprivate

map of OptO*/type of parent OptO, for navigation down the tree structure

map of OptO*/name of OptO for quick search based on name

Definition at line 203 of file Model.h.

Referenced by OptOList().

std::vector< OpticalObject * > Model::theOptOsToCopyList

staticprivate

Definition at line 230 of file Model.h.

std::vector< OpticalObject * >::const_iterator Model::theOptOsToCopyListIterator

staticprivate

Definition at line 232 of file Model.h.

std::vector< ALIdouble > Model::theParamFittedSigmaVector

staticprivate

Definition at line 288 of file Model.h.

Referenced by cleanParamFittedSigmaVector().

std::map< ALIstring, ALIdouble, std::less< ALIstring > > Model::theParamFittedValueDisplacementMap

staticprivate

Definition at line 290 of file Model.h.

Referenced by cleanParamFittedValueDisplacementMap().

ALIstring Model::theReportFName = "report.out"

staticprivate

the name of the report File

Definition at line 216 of file Model.h.

Referenced by ReportFName(), and setReportFName().

ALIstring Model::theSDFName = "SystemDescription.txt"

staticprivate

the name of the System Description File

Definition at line 212 of file Model.h.

Referenced by SDFName(), and setSDFName().