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=true)
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 = nullptr
static Model *	theInstance = nullptr
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 54 of file Model.h.

Constructor & Destructor Documentation

Model()

Model::Model

(

)

Definition at line 90 of file Model.cc.

             {
  //  theMeasurementsTime = clock();
}

~Model()

Model::~Model ( )

inline

Definition at line 58 of file Model.h.

{};

Member Function Documentation

addEntryToList()

static void Model::addEntryToList ( Entry *  entry )

inlinestatic

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

Definition at line 140 of file Model.h.

References mps_splice::entry, and 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;
  }

addMeasurementToList()

static void Model::addMeasurementToList ( Measurement *  measadd )

inlinestatic

Definition at line 145 of file Model.h.

References theMeasurementVector.

Referenced by Measurement::postConstruct().

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

BuildMeasurementsFromOA()

void Model::BuildMeasurementsFromOA

(

OpticalAlignMeasurements &

measList

)

Definition at line 1569 of file Model.cc.

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

                                                                      {
  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 = nullptr;
    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") ) {
      //   meastemp = new MeasurementDiffCentre( 1, measType, measName );
      // } else if ( measType == ALIstring("DIFFANGLE") ) {
      //   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;
      throw std::exception();
    }
    meastemp->constructFromOA(*mite);
  }
}

buildMeasurementsLinksToOptOs()

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 537 of file Model.cc.

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();
  }
}

BuildSystemDescriptionFromOA()

void Model::BuildSystemDescriptionFromOA

(

OpticalAlignments &

optAlig

)

Definition at line 1528 of file Model.cc.

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

                                                                   {
  theOpticalAlignments = optAlig.opticalAlignments();

  OpticalAlignInfo oai_system = FindOptAlignInfoByType("system");

  OpticalObject* OptOsystem = new OpticalObject(nullptr, "system", oai_system.name_, false);

  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);
}

cleanParamFittedSigmaVector()

static void Model::cleanParamFittedSigmaVector ( )

inlinestaticprivate

Definition at line 257 of file Model.h.

References cuy::ii, and theParamFittedSigmaVector.

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

cleanParamFittedValueDisplacementMap()

static void Model::cleanParamFittedValueDisplacementMap ( )

inlinestaticprivate

Definition at line 264 of file Model.h.

References theParamFittedValueDisplacementMap.

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

CMSLinkCleanModel()

void Model::CMSLinkCleanModel ( )

private

Definition at line 856 of file Model.cc.

References cuy::ii.

                              {
  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();
  }
}

CMSLinkDeleteOptOs()

void Model::CMSLinkDeleteOptOs ( )

staticprivate

Definition at line 868 of file Model.cc.

References cms::cuda::assert(), 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 {
  }
}

CMSLinkFit()

void Model::CMSLinkFit

(

ALIint

cmslink

)

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

Definition at line 800 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();

  }
  */
}

CMSLinkRecoverParamFittedSigma()

void Model::CMSLinkRecoverParamFittedSigma ( ALIint  cmslink )

staticprivate

Definition at line 1017 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

  }
  */
}

CMSLinkRecoverParamFittedValueDisplacement()

void Model::CMSLinkRecoverParamFittedValueDisplacement ( ALIint  cmslink )

staticprivate

Definition at line 1076 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 {

  }
  */
}

CMSLinkSaveParamFittedSigma()

void Model::CMSLinkSaveParamFittedSigma ( ALIint  cmslink )

staticprivate

Definition at line 929 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

  }
  */
}

CMSLinkSaveParamFittedValueDisplacement()

void Model::CMSLinkSaveParamFittedValueDisplacement ( ALIint  cmslink )

staticprivate

Definition at line 992 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 {

  }
  */
}

copyMeasurements()

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

staticprivate

Definition at line 1411 of file Model.cc.

References ALIUtils::changeName(), Measurement::copyMeas(), gather_cfg::cout, Exception, spr::find(), Measurement::name(), mps_merge::newName, AlCaHLTBitMon_QueryRunRegistry::string, Measurement::type(), and LaserClient_cfi::wl.

                                                           {
  //----- 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 = nullptr;
  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);
    } else {
      throw cms::Exception("LogicError") << "@SUB=Model::copyMeasurements\n"
                                         << "unknown measurement type: " << meas->type();
    }

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

    break;
  }
}

createCopyComponentList()

ALIbool Model::createCopyComponentList ( const ALIstring &  optoname )

static

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

Definition at line 742 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 true;
}

deleteOptO() [1/2]

void Model::deleteOptO ( const ALIstring &  opto_name )

staticprivate

Definition at line 1102 of file Model.cc.

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

deleteOptO() [2/2]

void Model::deleteOptO ( OpticalObject *  opto )

staticprivate

Definition at line 1110 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;
}

EntryList()

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

inlinestatic

Definition at line 86 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; }

fillCopyComponentList()

ALIbool Model::fillCopyComponentList ( const OpticalObject *  opto )

staticprivate

Definition at line 762 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;
}

FindOptAlignInfoByType()

OpticalAlignInfo Model::FindOptAlignInfoByType

(

const ALIstring &

type

)

Definition at line 1543 of file Model.cc.

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

                                                                    {
  OpticalAlignInfo oai;

  ALIbool bFound = false;
  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 = true;
      } 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;
}

getCocoaStatus()

static cocoaStatus Model::getCocoaStatus ( )

inlinestatic

Definition at line 63 of file Model.h.

References theCocoaStatus.

Referenced by Fit::fitParameters().

{ return theCocoaStatus; }

getComponentOptOs()

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 697 of file Model.cc.

References 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 = false;
  for (vocite = OptOList().begin(); vocite != OptOList().end(); ++vocite) {
    if ((*vocite)->parent() != nullptr) {
      //        std::cout << "looping OptOlist" << (*vocite)->name() << " parent " <<(*vocite)->parent()->name() << std::endl;
      if ((*vocite)->parent()->name() == opto_name) {
        opto_found = true;
        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;
}

getComponentOptOTypes()

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

static

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

Definition at line 666 of file Model.cc.

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 true;
  } else {
    return false;
  }
}

getEntryByName() [1/2]

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 615 of file Model.cc.

References DMR_cfg::cerr, gather_cfg::cout, ALIUtils::debug, and beamvalidation::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);
}

getEntryByName() [2/2]

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 120 of file Model.h.

References mps_splice::entry, and 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;
  }

getFittedEntriesReader()

static FittedEntriesReader* Model::getFittedEntriesReader ( )

inlinestatic

Definition at line 254 of file Model.h.

References theFittedEntriesReader.

Referenced by Fit::fitNextEvent().

{ return theFittedEntriesReader; }

getInstance()

Model & Model::getInstance ( )

static

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

Definition at line 80 of file Model.cc.

Referenced by CocoaAnalyzer::runCocoa().

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

getMeasurementByName()

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

static

Definition at line 633 of file Model.cc.

References DMR_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 nullptr;
    }
  }
}

getOpticalAlignments()

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

inlinestatic

Definition at line 137 of file Model.h.

References theOpticalAlignments.

Referenced by OpticalObject::createComponentOptOsFromOptAlignInfo().

{ return theOpticalAlignments; }

getOptOByName()

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

static

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

Definition at line 562 of file Model.cc.

References DMR_cfg::cerr, gather_cfg::cout, ALIUtils::debug, and beamvalidation::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);
  }
}

getOptOByType()

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 593 of file Model.cc.

References DMR_cfg::cerr, and beamvalidation::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);
  }
}

getParameterValue()

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 552 of file Model.cc.

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

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

  return iret;
}

getParamFittedSigmaVectorItem()

ALIdouble Model::getParamFittedSigmaVectorItem ( const ALIuint  position )

static

Definition at line 1267 of file Model.cc.

References DMR_cfg::cerr, beamvalidation::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);
}

MatricesFName()

static ALIstring& Model::MatricesFName ( )

inlinestatic

the name of the File for storing the matrices

Definition at line 102 of file Model.h.

References theMatricesFName.

Referenced by Fit::dumpMatrices().

{ return theMatricesFName; }

MeasFName()

static ALIstring& Model::MeasFName ( )

inlinestatic

the name of the Measurements File

Definition at line 96 of file Model.h.

References theMeasFName.

{ return theMeasFName; }

MeasurementList()

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

inlinestatic

Definition at line 88 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; }

MeasurementsTime()

static struct tm& Model::MeasurementsTime ( )

inlinestatic

Definition at line 135 of file Model.h.

References theMeasurementsTime.

{ return theMeasurementsTime; }

nextOptOToCopy()

OpticalObject * Model::nextOptOToCopy ( )

static

Definition at line 788 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);
}

OptODictionary()

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

inlinestatic

ACCESS STATIC DATA MEMBERS.

Definition at line 82 of file Model.h.

References theOptODictionary.

Referenced by readSystemDescription().

{ return theOptODictionary; }

OptOList()

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

inlinestatic

Definition at line 84 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; }

printCocoaStatus()

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

static

Definition at line 1503 of file Model.cc.

References COCOA_FitCannotImprove, COCOA_FitChi2Worsened, COCOA_FitImproving, COCOA_FitMatrixNonInversable, COCOA_FitOK, COCOA_Init, COCOA_InitFit, COCOA_ReadingModel, callgraph::cs, str, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                      {
  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;
}

readMeasurementsFromFile()

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

static

Definition at line 1281 of file Model.cc.

References DMR_cfg::cerr, Measurement::correctValueAndSigma(), gather_cfg::cout, ALIUtils::debug, Measurement::dim(), ALIUtils::dumpVS(), relativeConstraints::empty, mps_fire::end, ALIFileIn::eof(), ALIFileIn::ErrorInLine(), beamvalidation::exit(), Measurement::fillData(), GlobalOptionMgr::getInstance(), ALIFileIn::getInstance(), ALIFileIn::getWordsInLine(), GlobalOptionMgr::GlobalOptions(), cuy::ii, MeasurementList(), Measurement::measurementsFileName(), heppy_report::oname, 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().empty())
    return true;

  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 false;
  }

  //  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 true;
}

readSystemDescription()

void Model::readSystemDescription ( )

static

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

Definition at line 97 of file Model.cc.

References ParameterMgr::addParameter(), ParameterMgr::addRandomFlatParameter(), ParameterMgr::addRandomGaussParameter(), cms::cuda::assert(), buildMeasurementsLinksToOptOs(), DMR_cfg::cerr, ALIFileIn::close(), COCOA_Init, Measurement::construct(), OpticalObject::construct(), filterCSVwithJSON::copy, gather_cfg::cout, ALIUtils::debug, ALIUtils::dumpVS(), relativeConstraints::empty, mps_fire::end, ALIFileIn::eof(), ALIFileIn::ErrorInLine(), beamvalidation::exit(), spr::find(), ALIUtils::getFloat(), EntryMgr::getInstance(), ParameterMgr::getInstance(), ErrorCorrelationMgr::getInstance(), GlobalOptionMgr::getInstance(), ALIFileIn::getInstance(), ALIUtils::getInt(), ALIFileIn::getWordsInLine(), cuy::ii, 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(), and ParameterMgr::setRandomSeed().

                                  {
  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.empty());

    //----- 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().empty()) {
          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.empty()) {
          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 = false;
          for (vvscite2 = OptODictionary2.begin(); vvscite2 != OptODictionary2.end(); ++vvscite2) {
            if (*((*vvscite).begin()) == *((*vvscite2).begin())) {
              ofound = true;
              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) {
              const 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(nullptr, "system", wordlist[1], false);
        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 = nullptr;
      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("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 = true;
          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;
}

recoverParamFittedSigma()

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

staticprivate

Definition at line 1248 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;
  */
}

reorderOptODictionary()

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 504 of file Model.cc.

                                                                                                             {
  //  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);
    } */
}

ReportFName()

static ALIstring& Model::ReportFName ( )

inlinestatic

the name of the report File

Definition at line 99 of file Model.h.

References theReportFName.

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

{ return theReportFName; }

saveParamFittedCorrelation()

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

staticprivate

Definition at line 1211 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");
  */
}

saveParamFittedSigma()

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

staticprivate

Definition at line 1192 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");
  */
}

SDFName()

static ALIstring& Model::SDFName ( )

inlinestatic

the name of the System Description File

Definition at line 93 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; }

setCocoaStatus()

static void Model::setCocoaStatus ( const cocoaStatus  cs )

inlinestatic

Definition at line 64 of file Model.h.

References callgraph::cs, and theCocoaStatus.

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

{ theCocoaStatus = cs; }

setMatricesFName()

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

inlinestatic

Definition at line 155 of file Model.h.

References Skims_PA_cff::name, and theMatricesFName.

{ theMatricesFName = name; }

setMeasurementsTime()

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

inlinestatic

Definition at line 157 of file Model.h.

References theMeasurementsTime.

{ theMeasurementsTime = tim; }

setReportFName()

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

inlinestatic

Definition at line 153 of file Model.h.

References Skims_PA_cff::name, and theReportFName.

{ theReportFName = name; }

setSDFName()

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

inlinestatic

Definition at line 151 of file Model.h.

References Skims_PA_cff::name, and theSDFName.

{ theSDFName = name; }

SetValueDisplacementsFromReportOut()

void Model::SetValueDisplacementsFromReportOut ( )

staticprivate

Definition at line 1481 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

CMSLinkIteration

ALIint Model::CMSLinkIteration = 0

staticprivate

Definition at line 237 of file Model.h.

Referenced by CMSLinkDeleteOptOs().

CMSLinkRangeDetValue

std::vector< ALIdouble > Model::CMSLinkRangeDetValue

static

Definition at line 223 of file Model.h.

Ncmslinkrange

ALIint Model::Ncmslinkrange = 0

static

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

Definition at line 222 of file Model.h.

theCocoaStatus

cocoaStatus Model::theCocoaStatus = COCOA_Init

staticprivate

Definition at line 272 of file Model.h.

Referenced by getCocoaStatus(), and setCocoaStatus().

theEntryVector

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

staticprivate

std::vector of all Entries

Definition at line 191 of file Model.h.

Referenced by addEntryToList(), and EntryList().

theFittedEntriesReader

FittedEntriesReader * Model::theFittedEntriesReader = nullptr

staticprivate

Definition at line 280 of file Model.h.

Referenced by getFittedEntriesReader().

theInstance

Model * Model::theInstance = nullptr

staticprivate

Definition at line 176 of file Model.h.

theMatricesFName

ALIstring Model::theMatricesFName = "matrices.out"

staticprivate

the name of the File for storing the matrices

Definition at line 203 of file Model.h.

Referenced by MatricesFName(), and setMatricesFName().

theMeasFName

ALIstring Model::theMeasFName = "Measurements.txt"

staticprivate

the name of the Measurements File

Definition at line 199 of file Model.h.

Referenced by MeasFName().

theMeasurementsTime

struct tm Model::theMeasurementsTime

staticprivate

Definition at line 278 of file Model.h.

Referenced by MeasurementsTime(), and setMeasurementsTime().

theMeasurementVector

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

staticprivate

std::vector of all Measurements

Definition at line 194 of file Model.h.

Referenced by addMeasurementToList(), and MeasurementList().

theOpticalAlignments

std::vector< OpticalAlignInfo > Model::theOpticalAlignments

staticprivate

Definition at line 282 of file Model.h.

Referenced by getOpticalAlignments().

theOptODictionary

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

staticprivate

parameters

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

Definition at line 182 of file Model.h.

Referenced by OptODictionary().

theOptOList

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 188 of file Model.h.

Referenced by OptOList().

theOptOsToCopyList

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

staticprivate

Definition at line 216 of file Model.h.

theOptOsToCopyListIterator

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

staticprivate

Definition at line 218 of file Model.h.

theParamFittedSigmaVector

std::vector< ALIdouble > Model::theParamFittedSigmaVector

staticprivate

Definition at line 274 of file Model.h.

Referenced by cleanParamFittedSigmaVector().

theParamFittedValueDisplacementMap

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

staticprivate

Definition at line 276 of file Model.h.

Referenced by cleanParamFittedValueDisplacementMap().

theReportFName

ALIstring Model::theReportFName = "report.out"

staticprivate

the name of the report File

Definition at line 201 of file Model.h.

Referenced by ReportFName(), and setReportFName().

theSDFName

ALIstring Model::theSDFName = "SystemDescription.txt"

staticprivate

the name of the System Description File

Definition at line 197 of file Model.h.

Referenced by SDFName(), and setSDFName().